样本熵 近似熵 代码 应用于matlab

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
这个是我当时毕业论文的核心内容~~
需要做信号非线性的学弟学妹 可以参考下~~
. @0 D* M( V9 v4 x* Y" E* ~1 Q

  K) B/ r& m% \6 u


+ ~. A1 e- y: Q

1 k- s. ~+ b, E
. d" Q3 |+ W6 efunction [shang]=jss(xdate)
m=2;
n=length(xdate);
% t1 y: G1 Q. H, E; X1 e1 U/ mr=0.2*std(xdate);
cr=[];
gn=1;
gnmax=m;
while gn<=gnmax
. o/ T1 T* i2 D' ?3 n4 m6 Z& p- Wx2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
. n: l. Z! r2 O! h0 l2 }cr1=zeros(1,n-m+1);%存放
0 E6 N  f1 z: m8 I3 dk=1;
for i=1:n-m+1

for j=1:m

5 e+ W8 C0 X0 g3 B3 K/ Zx2m(i,j)=xdate(i+j-1);
8 q- V  c) i; v: E: T. c
9 i* o5 q( l- q  p0 Zend
3 |; }# l0 D# `. y: @" uend
x2m;

* k4 Q% e" [& P' m" a
for i=1:n-m+1
6 Y2 Q8 A6 [! o9 z4 f
for j=1:n-m+1
% H$ t4 z( \7 p
" p9 ~7 i8 @/ ?8 f5 Hif i~=j

d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
0 I+ m. A$ B0 b
( [  Y1 w0 G% ^k=k+1;
1 |. I& B1 z/ ]' b+ G! a
end
end
0 H5 l: G; t" W' d, Y/ B  B
/ ?" X$ ?3 ]# _* U7 Wk=1;
end
d;
3 `( k( D; ^( T
for i=1:n-m+1

# S/ O. j) G0 B9 Z6 o8 _[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

cr1(1,i)=l;
end
& h/ t0 \- Y! j2 pcr1;
3 s( N& f1 x4 C& j$ P0 {; x
- q  r0 X, f% {9 |  T! kcr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1

. J2 _& H' t% _& U& B7 i3 K/ iif cr1(i)~=0
( \, @) A9 s% m( A2 i
sum1=sum1+log(cr1(i));
- ?6 @+ ?  ?& l
end
& j6 F5 H- ^1 o6 j! n3 r7 j) y" dend
cr1=1/(n-m+1)*sum1;
1 _' ^9 d1 U2 Zcr(1,gn)=cr1;
gn=gn+1;
m=m+1;
end
cr;

shang=cr(1,1)-cr(1,2);

5 t. |) b/ X8 ?' Y
, S7 ^4 b: [0 Q5 C- d% H

1 V1 ~" e; W9 V' \6 s3 r  xfunction [shang]=ybs(xdate)
# X% N) N! T( @7 Pm=2;
- L; ]/ Z1 H7 D8 e1 G' un=length(xdate);
r=0.2*std(xdate);
cr=[];
gn=1;
gnmax=m;
1 g7 G: [) m4 Y1 }2 c( W4 dwhile gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
6 b4 _+ B6 d) T/ G* H  ?, b+ xd=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
# ?$ s5 e3 p, ]% z# Dk=1;
8 E+ ?  T# ^1 V1 w9 nfor i=1:n-m+1

for j=1:m

x2m(i,j)=xdate(i+j-1);
, v2 N  ]& f, U, P1 ~8 x
end
end
+ a6 @5 b; j& \2 ?4 sx2m;
+ Q! j" r" @" u+ \& V

for i=1:n-m+1

% M' G# T0 ~; h, z: Lfor j=1:n-m+1
1 k4 A+ g+ ^3 _' j: A$ A
: Q1 u% t0 ^: k* @& @if i~=j
7 n; Y1 T4 U2 L6 M6 D, |4 s8 r
  M+ `( g5 O6 a" z! \8 bd(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

3 U: ~2 o6 e) |8 yk=k+1;

8 l+ I! k1 H+ m' Tend
end
5 Q6 d, K5 A# F( G  D0 Q
k=1;
9 M7 i* j  [; Y5 M- n% }2 Jend
d;

for i=1:n-m+1

[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

cr1(1,i)=l;
" I# [) O7 D5 |( s* @3 ], K( ^0 oend
cr1;
0 C1 R; f% ~% d/ d3 n% S
cr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1

sum1=sum1+cr1(i);
9 X. h7 e& {' V% {7 O
# y, e" q1 r0 O7 }: W# h1 B9 ~/ aend
end
1 E- n7 `, O. O/ vcr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
gn=gn+1;
5 T1 R3 o: m/ X. u& O$ Q( A  Km=m+1;
3 Y: r2 ^9 M$ K; w& a; k$ Cend
+ l+ a3 [9 ]- R2 F& O% [3 K5 Ccr;
9 P1 @4 g3 R" l8 n( r, U: Y, Jshang=-log(cr(1,1)/cr(1,2));

! p" t; A. ]# D  H[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

其中圆脸为
: ) 中间没有空格