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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
这个是我当时毕业论文的核心内容~~
2 n/ E2 [$ V# ]; h! V需要做信号非线性的学弟学妹 可以参考下~~


1 U' B1 L9 Q; j% ]
4 b; p7 v" y$ @) h( o+ y
! c6 b/ J6 |4 t5 Y

9 ]8 p( D2 t# x+ N3 H
2 \0 C  y" S" a
6 K: ~/ Y, v- zfunction [shang]=jss(xdate)
m=2;
2 ?* N3 @! O7 R6 ?$ X( {n=length(xdate);
r=0.2*std(xdate);
cr=[];
' b) S4 R3 x+ ^# A3 @gn=1;
gnmax=m;
while gn<=gnmax
  B) C, ?$ f. u! A: ?7 D2 e, Kx2m=zeros(n-m+1,m);%存放变换后的向量
  w5 h1 O  |" H1 L5 [; z2 id=zeros(n-m+1,n-m);% 存放距离结果的矩阵
6 z- V+ C: [% j9 |$ ncr1=zeros(1,n-m+1);%存放
k=1;
for i=1:n-m+1
: Q$ r! j6 ~8 p7 c. Q6 G9 V" d3 y' ]+ C
for j=1:m
8 M3 l$ r3 f* c) U
x2m(i,j)=xdate(i+j-1);
( q: v6 t6 I# y. R! [3 y
% I+ k# O+ d  P5 k8 ^# Yend
3 y6 H& b+ _( ], o6 dend
! y) Y: X4 R2 Mx2m;
) J, S3 L! f; G: U- b

9 Y. S* a. y: v0 D1 Ffor i=1:n-m+1

for j=1:n-m+1
! C! W8 T% Q1 g% t+ M: a# V
if i~=j
: }6 j3 h7 Q& Y* A8 [
d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

7 a- r1 u0 ^/ Fk=k+1;
  E' [" f; K3 M: s
end
end
+ Q; p2 H8 u4 H! p" r5 \" W
' x1 w) [" [4 W  W2 K6 ~4 t8 |k=1;
end
d;
) d$ j: g/ E5 X
  Q+ E, ?6 i2 r( I" `" c, T* ]5 Cfor i=1:n-m+1

0 t: k/ M) k" u0 v2 R8 R[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

cr1(1,i)=l;
7 x9 G+ Q& u# k' \end
2 d0 C3 ^6 P; S/ L' o% Y  D. Jcr1;
- v$ f4 x6 O" q1 F
cr1=(1/(n-m))*cr1;
sum1=0;
# N$ S$ T0 @) J3 Sfor i=1:n-m+1
' ^; k1 o5 N# K! |) Z
- a' v8 r; @& {0 Qif cr1(i)~=0

sum1=sum1+log(cr1(i));
. e7 i$ {1 q# s' m, W
$ s2 \* p* N) \: ?5 Fend
end
  Z1 \  Y6 k* Y) f$ }cr1=1/(n-m+1)*sum1;
5 T: G1 U+ j- [: Wcr(1,gn)=cr1;
gn=gn+1;
$ g$ A. l2 M# ?  c9 r, v* C7 B; qm=m+1;
end
' t# Y" R% O$ P0 \cr;

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

$ L. u9 t+ O8 S& j$ n- ?0 Z/ j
$ M2 m8 \9 @5 `& g" _

) T; b, R( k0 ]/ Rfunction [shang]=ybs(xdate)
m=2;
# [. S) K! [4 L' d5 Tn=length(xdate);
6 B" r3 D1 e: @, c5 y) T2 Zr=0.2*std(xdate);
cr=[];
5 n8 S3 E' I: N/ F. Zgn=1;
gnmax=m;
1 C# a: e& I: ?% p1 ?  uwhile gn<=gnmax
0 l" r$ y2 V8 \x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
8 J. }; A* E7 Rcr1=zeros(1,n-m+1);%存放
k=1;
+ _" n9 h9 u+ g' X% G) R  mfor i=1:n-m+1

for j=1:m

x2m(i,j)=xdate(i+j-1);

end
0 n- x3 d2 i% Z( X$ u* e7 R" [! Z! Hend
x2m;

+ G% ]5 D0 d  F9 l6 n: u% q, l, Q
/ p5 g, `0 C9 m3 Ffor i=1:n-m+1

- C# u9 d, E& _5 Q4 C6 z3 `for j=1:n-m+1
6 f) u. `6 l4 A: ~" [4 w7 k
- v- ]1 ], O4 r4 Wif i~=j

2 _  e8 x8 J- \/ wd(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
/ m, R& C# q' S8 X# B0 r
/ P6 S8 E7 I. O4 E3 Mk=k+1;

end
end
$ ?4 q- |+ b$ y" Z# q
4 E1 J# ~* t/ S% J7 @! b6 Uk=1;
end
d;
8 p8 \& h/ c- O& u7 \
for i=1:n-m+1
/ ~- L* K# J2 O/ k2 L( K0 F$ t2 I
[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

# r6 E3 f' z' G3 M1 |0 ~: gcr1(1,i)=l;
7 V6 W0 G" _( p! cend
cr1;

: s* V( }6 a) L: I" x; bcr1=(1/(n-m))*cr1;
sum1=0;
$ M$ B& A% L4 F2 C% C: t# {for i=1:n-m+1

' Y) {# e7 C& ?& U1 T; l% p3 q# Vsum1=sum1+cr1(i);

: z: u$ l8 I& @% i" ^) g" Nend
end
7 ?; {, Y3 b9 T# s) Q; Ecr1=1/(n-m+1)*sum1;
# k) w) @0 Y% m% K4 Fcr(1,gn)=cr1;
2 {# w- Y  L1 F: W9 z8 pgn=gn+1;
1 N5 {3 W$ V7 N0 x: o1 Xm=m+1;
end
cr;
0 C0 X/ U$ A  f. x7 a( ~' Gshang=-log(cr(1,1)/cr(1,2));

0 W3 l; u# L. v3 s( ][ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

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