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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
7 T) U3 p& x& U& r1 b) C  C+ r这个是我当时毕业论文的核心内容~~
需要做信号非线性的学弟学妹 可以参考下~~
9 k" v$ k# b, C: e4 E$ w



- ^6 F7 X* C% m! c0 m5 A, M



function [shang]=jss(xdate)
m=2;
( I) e/ |/ E4 [2 y* [n=length(xdate);
7 D6 x! w  {  n; F* J, sr=0.2*std(xdate);
cr=[];
7 B5 s# |0 O4 w8 Ugn=1;
2 R, F2 j* x1 u1 Dgnmax=m;
& N; `  n7 J' y' E" j) zwhile gn<=gnmax
3 ~: z  ?! M- d- Y6 f9 Jx2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
4 q" [. a6 V  ycr1=zeros(1,n-m+1);%存放
k=1;
( T! o5 u, q  a$ d! zfor i=1:n-m+1
6 x5 X0 M$ R3 L
for j=1:m

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

end
end
4 T3 p. n8 S2 }1 w5 @, Ix2m;


9 c; {  D0 q: C) S7 K# G' B2 cfor i=1:n-m+1
1 X- N) X# O& C4 a
' V& y" s1 ^2 K, q8 q: w/ t5 K" M) ?for j=1:n-m+1

if i~=j
; u8 K" ~; Z* I! Q
d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

k=k+1;
* U, |' _! e6 a
end
% b4 t* Z0 B" g, @end
: t5 ?% V  k: \" n: s
' ^& J' ]4 b! h' e. c* J" Kk=1;
9 m$ p1 T* |+ Q; O+ X5 Lend
8 d2 U4 \5 a: R; @d;
# a) {8 H: E! v9 U; N2 b' @
; b) w/ i7 C# K: ^+ N7 Cfor i=1:n-m+1

# \, m; z# R/ l6 K  a# u1 K[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

0 G5 I# r! t4 x! ~. ycr1(1,i)=l;
! _+ i) S5 l7 A! r+ jend
cr1;
& P, b9 k- a( h
cr1=(1/(n-m))*cr1;
sum1=0;
% v! ~! G6 c) T% O. T. _for i=1:n-m+1
# K8 s9 o1 v9 z% n7 ?8 q
if cr1(i)~=0

* v5 p; H  k2 t; T# x' q- s$ fsum1=sum1+log(cr1(i));
/ f% Q1 w) N& U. C
end
end
cr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
gn=gn+1;
  ]1 p. w! h7 hm=m+1;
8 I( K( o# Q6 V  }$ Nend
. W; M  ?( F& r+ y- H; x8 Pcr;

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

7 r' \9 H) P2 r+ t5 @4 b
& K% L5 X9 q. o0 S! H' {function [shang]=ybs(xdate)
' a0 i3 l+ s6 A, B; p6 v4 `/ t8 ^6 l2 Am=2;
' j, ~5 |# P" a4 d  u' {5 rn=length(xdate);
r=0.2*std(xdate);
3 ~* ^* e# c4 c8 O# ~cr=[];
# A0 n- Z' N7 _3 c: ^1 bgn=1;
gnmax=m;
; @$ m" H2 }. r1 d: j/ j( ?while gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
5 |) @" w3 N9 a( M; |% U! xcr1=zeros(1,n-m+1);%存放
, ^4 L3 B" ^% d  U1 b6 o  @1 C1 Ok=1;
for i=1:n-m+1
8 G$ {9 h& T* b3 r5 Z" d
for j=1:m
! z2 J) T# j$ O
x2m(i,j)=xdate(i+j-1);

end
5 C9 _: M% J2 v9 n+ eend
/ j' k+ a% G3 s6 z- ~x2m;

. i3 M1 y: D% o; C( ]  j$ g5 T
for i=1:n-m+1

for j=1:n-m+1
% O. J2 H1 K9 o. w) D1 @9 _2 V
- |4 [: @/ T3 b  U. U% ~if i~=j
8 Y1 O0 T$ t, v6 D, B
6 _% {  I  c! f, o) M$ Sd(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
3 |3 i8 i2 H$ v/ t3 s8 W/ V( B, f
, R2 _" S& V! z# nk=k+1;
( z6 o+ c6 i4 E4 q
end
end

k=1;
$ w' A8 S0 d; T4 r9 q4 {5 I. G6 }5 send
d;
$ P. [& Z# Z1 U6 a+ V  F
for i=1:n-m+1

3 _5 L  L$ c  a7 ][k,l]=size(find(d(i,<r));%将比R小的个数传送给L
8 v; n6 @  c: r0 N6 Y% W4 F
cr1(1,i)=l;
& ^2 J3 f- ^. J5 O9 A4 k8 r& Uend
" E2 m0 v( L' j. w. g8 |9 _cr1;

cr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1
3 T- P. t" f7 x  a: c
sum1=sum1+cr1(i);
- @) t# l0 Z$ j8 C+ S
end
" d6 }& R5 j9 bend
cr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
gn=gn+1;
" y0 I4 m$ Z3 `m=m+1;
end
cr;
: p" h+ M' J- x7 j  M/ {$ k3 Eshang=-log(cr(1,1)/cr(1,2));
+ Y0 c* U, n5 y, `  X
[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

其中圆脸为
% k. \( x# z  S4 N9 v7 k: ) 中间没有空格