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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
. b. I8 w- F6 B# r" D这个是我当时毕业论文的核心内容~~
需要做信号非线性的学弟学妹 可以参考下~~

" ]& j$ m  {% m6 _. i

4 @9 R7 j: n' Y* I  N' \, p
9 d. ?$ R% x5 d6 ~
6 w1 \6 s/ A/ @
$ m+ R, [# r2 H7 Z5 k! t  [

0 x$ w* Y/ b# I3 n7 |- Gfunction [shang]=jss(xdate)
m=2;
9 ]7 B1 @* e) N! [, mn=length(xdate);
  [  ^2 K. b# w3 Q: sr=0.2*std(xdate);
cr=[];
gn=1;
5 Y: {: g) s" X% S. S4 ]8 z7 S, @gnmax=m;
while gn<=gnmax
' p* ^# H- w: @, W; d( O) X  Qx2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
: S% e& S7 B- x5 w6 _1 ~k=1;
for i=1:n-m+1

for j=1:m

x2m(i,j)=xdate(i+j-1);
* F& r; h/ m' p/ f" o9 t$ ?
end
end
) E3 t2 t7 L8 W. }/ y2 {8 d. Gx2m;
8 G! o  b  }  ~5 Q7 w  j' n, A

6 X8 R9 }9 V0 l/ K. }/ q3 T+ }for i=1:n-m+1

! V0 C' \8 G& ofor j=1:n-m+1

$ I! J  G9 o, M& E, h! ?! dif i~=j
- f3 _" J  f# v7 |: S1 s$ z
% |6 i5 d. d% I- f& o- pd(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
0 ^2 F1 ?1 c4 p' Y% ]/ r0 g
k=k+1;
" G  ~* q" I* A& ^; v1 i  G! n
end
: h/ b9 w8 o$ a8 e5 t, S+ e% Lend

k=1;
end
" L, D5 ^/ D3 N3 G/ H/ j7 Wd;

: }) Z6 v1 g6 a/ V3 Q) Ufor i=1:n-m+1
7 O; V, p+ Z( I7 c$ H
[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
7 j$ s2 Y' u6 b* K- R$ l
cr1(1,i)=l;
- `; n$ ?% F8 Z7 i, h- T' Dend
# w+ P) i8 S7 a- V; Q, Vcr1;

1 Y- i: \% T/ P/ u! [cr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1

/ }, v  O1 W: ]( \. o# ?! p$ eif cr1(i)~=0

3 _) V* V' u5 V* m4 i' P2 m9 q4 _sum1=sum1+log(cr1(i));
# }* g# Y9 K0 y) D, h
end
end
% f" x. z3 r- ^+ S$ R4 m  t5 Rcr1=1/(n-m+1)*sum1;
. N" T* i" w- ^. tcr(1,gn)=cr1;
* S5 a! z: l9 {" i* i0 F* A, Wgn=gn+1;
m=m+1;
end
cr;
* k* P% P, |" f

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

/ m$ t. F) x0 Q

. t: C, V% w4 `4 {& d

' s" ?* p1 h1 j) X% ?; ~function [shang]=ybs(xdate)
# }. W8 e% J8 |3 t* a$ vm=2;
n=length(xdate);
r=0.2*std(xdate);
2 k$ x' B3 B* N7 }. L) B  ]  A- J, ccr=[];
gn=1;
gnmax=m;
while gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
3 ^8 {' b7 v  G( P- x( e4 Ycr1=zeros(1,n-m+1);%存放
' }2 L9 U+ ~" g: Bk=1;
9 G% o) x7 m- Y, f7 c* ]% p* Ofor i=1:n-m+1

9 V8 c3 f5 x7 t6 a% D- I: Vfor j=1:m
3 \, Z' M, t7 S; Z. _* h7 M, x* U, @/ l
6 b' V3 A4 \4 ^; s& nx2m(i,j)=xdate(i+j-1);
; N" D' G6 ~( U& I7 r7 Q6 S
/ \9 Y4 e$ ^+ `3 \) h- Oend
end
6 `) ^" ^; P. W4 d/ t5 Hx2m;
* F9 ]  p# ]; F/ v) Q+ D
& X, T1 Z5 ^0 L/ {
1 _% V3 C- n! I* Z. A/ Sfor i=1:n-m+1

8 V9 `- S5 _3 e7 d; i: o: _for j=1:n-m+1

( a" ]- w. X. n! w* E6 hif i~=j
" h( Y( E2 O1 I$ [0 r* `: I7 L
d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

3 K5 k) K! W/ ?, l2 {k=k+1;

  y8 z1 S5 ?& q4 ^end
end
8 c; r& h' y( V8 X3 o) p
k=1;
, z- s, x- s& i" S. \+ }- rend
d;
% R+ k+ i0 ?$ j# R0 T3 @2 X
for i=1:n-m+1

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

  b' p% |4 R" q* n, f# Mcr1(1,i)=l;
$ o2 D: ~+ X; \0 D" \1 T- G4 G: W6 Send
cr1;

% u; O. ]% e/ V6 ecr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1
" d# o3 h% I6 u  P
' |6 ~( [1 o: B1 msum1=sum1+cr1(i);

& p" ]. j+ K7 pend
end
6 h+ R1 h2 ^( E! i- d" kcr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
+ y2 m% g* E0 r* ]3 s# w9 Sgn=gn+1;
m=m+1;
end
  v: S: p. Q1 b3 P" p0 v$ [- B2 kcr;
shang=-log(cr(1,1)/cr(1,2));
. _# O3 h; }. f; p. `' Y: I
- Y" G: z, z. G/ w+ h- S- m[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

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