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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
) f' z$ b. u5 @" c这个是我当时毕业论文的核心内容~~
" I+ x* P9 `. z+ m需要做信号非线性的学弟学妹 可以参考下~~

9 [1 o- ?  G& o) z
7 o0 B0 ]# ]! C% g, G- V* ^
. n& J) K2 }$ {5 E; E

$ m6 A# m. x* j  v0 f+ M
! B4 ]. f5 Y" {2 e6 W3 M

function [shang]=jss(xdate)
m=2;
n=length(xdate);
! G; |& D, A' h6 z+ m$ R1 Ar=0.2*std(xdate);
cr=[];
gn=1;
gnmax=m;
while gn<=gnmax
  h7 r" R: V+ Fx2m=zeros(n-m+1,m);%存放变换后的向量
2 ~& n1 R! y. ^: V  gd=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
, ^, V  f: f* d# n0 m' \; l$ G/ Ak=1;
for i=1:n-m+1

, C) s7 x6 i2 A+ Z) y* @for j=1:m

x2m(i,j)=xdate(i+j-1);
' w) ^+ g- P9 Q, Y. I( x, d
6 w% s# ^( f  jend
0 M) B. b" P$ ~end
, P- K9 w& d+ `" Mx2m;

( Y0 i# k& f7 V! S
5 [7 r% m' O( d2 y$ jfor i=1:n-m+1

! u6 |. v) D) _1 Yfor j=1:n-m+1
) n9 t) a; d. Q* c7 t* T
+ D  u4 _1 i8 j" P6 z) N, Iif i~=j
/ ^( o; Q. K9 W" q1 L
d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
, }7 O! k% V  _) d# v9 a9 }
k=k+1;
4 u0 f+ y# U1 o1 N! q/ e
% r* {+ ]6 x5 K# vend
end

6 _- Y( }! Z+ U" Hk=1;
end
d;
$ O3 l% P5 _  @% [
+ a2 C- S' S1 ?4 I6 _5 o3 @for i=1:n-m+1
& I1 x. t4 p# P0 s
[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

6 Q/ B; H# }! ?cr1(1,i)=l;
end
cr1;
4 A" x, `3 s4 ^& b- r# ~3 R/ _) S5 r
cr1=(1/(n-m))*cr1;
+ L6 l5 W' e$ ysum1=0;
for i=1:n-m+1
1 ~0 y9 C' @/ x) C, s; N/ m
if cr1(i)~=0
8 p- l* P. l2 h* j/ W, ^# q
% v" }# |2 H1 asum1=sum1+log(cr1(i));
- R" h! i; X$ m2 Y! p% M4 R, c' ]9 _7 a
end
4 r- {1 ]( Z  R! j  F8 Nend
cr1=1/(n-m+1)*sum1;
  b( f2 v; |1 ccr(1,gn)=cr1;
; Y1 ]7 ~4 g" o! ~' xgn=gn+1;
m=m+1;
end
cr;

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

$ ?2 C& Z2 |- L0 E% C$ H

+ s" U# j( C3 \& U* O! z0 }6 P6 Y

8 o( }) l; }! {) D, B% bfunction [shang]=ybs(xdate)
m=2;
n=length(xdate);
r=0.2*std(xdate);
8 L& e/ A& r5 ]/ s' t4 ]' ncr=[];
gn=1;
( V% w( }% A, o6 V7 d1 H9 P# `* Ngnmax=m;
9 G7 |& f4 j/ W+ j7 c  c  l3 w5 Jwhile gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
6 \+ e6 j# b: z& u2 `d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
( h2 U9 ^/ n' e& R/ w2 pk=1;
for i=1:n-m+1
# P5 p( n. x. }$ X. m
  Y; Z& r: }' a1 H7 H# ]for j=1:m
! g! ]4 F  e7 G( w+ F
) |( k) J$ M* W% R+ h: Xx2m(i,j)=xdate(i+j-1);
+ }$ W& h8 e2 e
- r( [0 ~5 B2 j: [1 mend
" M3 [2 Q9 ~" t$ hend
; e  `# Z9 U* vx2m;


for i=1:n-m+1

$ N' @2 y7 P; ~1 O2 q  R0 ~7 jfor j=1:n-m+1

: ]+ q( ~/ W; K7 sif i~=j

d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
1 N. e; n% Q) P+ g; J
k=k+1;
% A* F$ {+ S& V. @& ]% x4 I# c
end
. V" f, U% ]) \* i( M2 O+ X8 @end

+ b: c# Q. P, N! l6 `5 d! W( |( Bk=1;
& X  u5 w; B  Y1 U, aend
d;
' Y* F) ]& x. W2 q# B, x
. d; n8 ~0 i1 ~$ v! ?3 z5 @for i=1:n-m+1
% u( `/ {: {) v% s
/ R2 j$ u- v; @& i( k! R. k0 S[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
3 A) i) q2 u1 _
cr1(1,i)=l;
end
/ J& u: Z2 o- Icr1;

cr1=(1/(n-m))*cr1;
sum1=0;
; d! \. q8 a+ L& P4 W" Qfor i=1:n-m+1
9 C; s5 b5 Q/ b3 ]; i, I  Q
% b0 ]" r7 ?+ h5 U9 t7 csum1=sum1+cr1(i);

- V2 Z  ]9 R6 S' zend
" ^* e6 }  V2 B2 C' [" Aend
cr1=1/(n-m+1)*sum1;
8 x  p* j) ~$ S) g: ]0 Ycr(1,gn)=cr1;
gn=gn+1;
m=m+1;
  k/ C! h' o! Z* L7 c* |end
cr;
' y* Y4 U' O* j8 G3 G3 m( }/ Zshang=-log(cr(1,1)/cr(1,2));

5 H2 `- }. R9 O& a3 [+ ?$ r% K. S' r[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

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