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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
这个是我当时毕业论文的核心内容~~
需要做信号非线性的学弟学妹 可以参考下~~

7 W7 m3 ?  \4 I) \' r" W


3 J2 F' h% a" }



function [shang]=jss(xdate)
m=2;
n=length(xdate);
r=0.2*std(xdate);
cr=[];
gn=1;
4 N  j# G7 c) ]! o( ]$ p! Y* \/ J+ Hgnmax=m;
; E0 u9 C7 y- c. e0 {* |% }' mwhile gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
1 k  c+ q+ {" s% P. ]. Qk=1;
for i=1:n-m+1
# ~! C  G4 F' z3 q9 r1 g$ s4 p
for j=1:m

+ u) T2 j! y7 F# H$ d$ ox2m(i,j)=xdate(i+j-1);
; L5 T* y- m6 {1 V; _
; |& y3 J9 j& e% N* v+ `1 l9 bend
, D. A' W2 _# u) E8 v# \" Hend
7 S; @; Q6 E' }) E# O. F; c0 v6 M& lx2m;
5 ?* v! w- f0 @0 r6 h, l+ |' j
: b# c: T: a0 o
0 n! D; K- |% q0 K# Zfor i=1:n-m+1
9 f0 E+ W( M0 J% o" x) O
% o( @( V* O$ X1 d0 j" P- j) D  V) Rfor j=1:n-m+1

if i~=j

6 U, ~9 b* G$ u: cd(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

k=k+1;

/ X" X2 h& ^) V# N/ _+ Fend
end
/ ?4 Q' K: n1 A+ p, i9 a# ^" K, v4 x
k=1;
end
d;
% |' r. h3 ~; a+ o
for i=1:n-m+1

1 z8 V8 N2 r  i, e/ M3 C+ W[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
# C2 |. ^9 E( m1 E
) d" W2 d) q/ p' }) p+ Ccr1(1,i)=l;
: J" _1 x( E% i( ^, L9 Z; F: ^2 jend
6 J* H: t% D" N/ Y. @7 Q+ Y4 X" `cr1;
! I, o9 p$ t* V, `. \' r
0 [* a4 T( k( @% N( Qcr1=(1/(n-m))*cr1;
sum1=0;
3 M* E6 n2 |* W$ {for i=1:n-m+1

if cr1(i)~=0
) Z& x* _" H, d% S, y- V
( `/ B0 R7 V8 s+ ?5 S. R9 Asum1=sum1+log(cr1(i));
/ \: `7 p$ Z, J4 N9 ]% q' B8 _3 L
end
end
0 ^6 w/ o2 y. Acr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
9 C: R8 b" _7 e' @' F8 s. Mgn=gn+1;
m=m+1;
end
cr;
7 X! |) `) a' [

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

* l2 m3 x; y( B7 m
3 k5 O" k: P- j6 w  w

6 o8 r: A, A2 }* u# Z4 \
' u% g6 a% w/ P$ e. \function [shang]=ybs(xdate)
m=2;
9 ?5 B$ K% T% x- i( sn=length(xdate);
8 a" E  B6 p+ ^; ^: Q: Mr=0.2*std(xdate);
. L( a. ?( f- V4 Mcr=[];
gn=1;
" [* N$ A/ C  i: u6 T% f% A# U, O3 e. Tgnmax=m;
while gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
3 ^' _/ M) S! n8 y" n& \7 @! L+ p& ed=zeros(n-m+1,n-m);% 存放距离结果的矩阵
! ?5 ?7 f& X: R9 x  Z( }- _cr1=zeros(1,n-m+1);%存放
k=1;
6 p+ d' d) q* f; ~- ^# _3 rfor i=1:n-m+1

for j=1:m

x2m(i,j)=xdate(i+j-1);
: g1 i1 Z9 }# N' w
end
end
x2m;

# p5 V' d# I4 s+ f
. Z% z" z, R( ^7 o8 [& P% U- a5 ?for i=1:n-m+1

for j=1:n-m+1

  H9 A( u; H. ?* Z: ]9 M7 c2 [- Gif i~=j
3 t; ^" g- u4 o' v" n- ?
d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
! Q2 F' k- d8 ]* X1 V
k=k+1;

  R* N- j* u" v+ Y& Qend
end

8 d+ G0 {9 n0 t6 d$ ik=1;
end
d;
3 W1 I9 B+ [* ~7 Z5 ~
1 }/ ^# x' p+ X+ k8 Z! S+ a3 Efor i=1:n-m+1

$ _: |: p% J* p& H. |1 u[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
: [, O7 p+ n8 _) S
cr1(1,i)=l;
0 O# J7 V3 ]( ?9 hend
cr1;
( t( b5 P$ s8 j
cr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1
. g' p* N0 i8 V$ x
sum1=sum1+cr1(i);
( F$ N1 o4 l! h2 f+ [, z
end
4 D* V( J. S% r" jend
9 H+ V1 H& K: {% S- u& \3 u, a* }, Lcr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
8 t, p9 S8 e1 v& f) Ygn=gn+1;
  }9 y5 X( L9 r2 v2 Bm=m+1;
end
/ _7 n# h# d+ j1 K7 zcr;
shang=-log(cr(1,1)/cr(1,2));
* `6 ^+ Z# \2 _$ x
+ J. ]  k8 w' Y( j[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

其中圆脸为
1 t  j* J3 N& \  E! \1 }* [# r: ) 中间没有空格