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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
9 S9 v7 ^% |6 G: M$ w4 n这个是我当时毕业论文的核心内容~~
& V% b' x; A  r% ]需要做信号非线性的学弟学妹 可以参考下~~




! w' C% p* L9 i$ R6 M. x6 @


# T4 u; p' M- o4 f( h6 L
7 n4 e$ v' |; ?* t6 J* hfunction [shang]=jss(xdate)
+ N, u, L7 F- c! {5 [+ _m=2;
" b6 c6 o# n& n; R1 x4 H  _' i" Nn=length(xdate);
1 H: S+ }9 {& H+ ?# a  p5 v/ }r=0.2*std(xdate);
. b- N5 t: z: T) X* Scr=[];
# g5 M. B7 A( _- agn=1;
gnmax=m;
while gn<=gnmax
; l! E5 [; H! v6 A& b0 }x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
k=1;
  ]2 X6 c& z' j8 nfor i=1:n-m+1

, ^) {. {) Z9 T1 l2 g" efor j=1:m

) W1 c) @% u+ |x2m(i,j)=xdate(i+j-1);
9 ~% r3 O3 P, e
( Y' T  Z: i( g5 e# B; }end
end
- z$ S2 j1 p! }4 w9 X2 b9 \& i8 tx2m;
% @: O# j% ^' R8 Q

2 ]% t" G8 _+ g/ \+ d! dfor i=1:n-m+1

4 k* J- s1 L/ M0 v5 ^for j=1:n-m+1
2 X& o$ y% x4 n6 P) B  {8 m
5 B3 n7 M. W, z( H! Pif i~=j

d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

k=k+1;

1 ]0 G/ f. F2 e  o; a; Send
end

k=1;
end
' U: f( S4 v4 }d;

for i=1:n-m+1
: |/ P* `  d, U5 K
2 u' M  V7 W% G6 z+ o0 U! [+ H, h/ T[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
5 Z+ O) |5 z- t8 v  M2 j* T# k6 C1 J
, G) U: {" j+ h4 P& U$ m7 Qcr1(1,i)=l;
end
' ]' Q% o& s7 y4 Ocr1;

cr1=(1/(n-m))*cr1;
sum1=0;
* u8 `! f& [. s3 J7 W6 `; Gfor i=1:n-m+1
" E; T% ?( _: v( w
2 w2 l' k- d2 `7 eif cr1(i)~=0

sum1=sum1+log(cr1(i));
7 z+ l" ?% y0 h1 p- z/ F0 j* A
end
end
$ P9 a4 X, W( Q! H- ^, p: ~. ?cr1=1/(n-m+1)*sum1;
' I' F& u/ p. x' E' f6 f6 zcr(1,gn)=cr1;
gn=gn+1;
m=m+1;
end
- ^- D. j- S0 y2 b) v: {2 Ncr;

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

5 U& a6 l# s& r+ `6 l/ q/ ?! q
2 E9 R# ~, O! t- o% N( @& Q
1 V# y; ?" L2 t: T2 y
function [shang]=ybs(xdate)
. w$ [. P" z0 G) x! lm=2;
n=length(xdate);
r=0.2*std(xdate);
/ M* _* E" s3 }* Dcr=[];
4 i/ d4 @$ a( ^1 \$ c1 |+ jgn=1;
+ `" x8 ?2 d( P/ e* L( `+ dgnmax=m;
" b6 [" x* T5 _7 cwhile gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
9 T- ]+ Z- L5 s% S2 i5 Yd=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
k=1;
for i=1:n-m+1
1 D0 ]* Y+ f/ M% I) ~  W
for j=1:m
# ]1 v0 b$ i1 H+ C& K, N
x2m(i,j)=xdate(i+j-1);
# i  b( `2 B. \8 X7 A) l
end
: L" O1 X1 x. ~9 G8 ^) xend
x2m;
) ^! H% |! K8 y

) B: N2 i0 _- xfor i=1:n-m+1
  L& ^. ~0 |" R; W* F( i, i' v5 G
for j=1:n-m+1
1 \6 H  u8 y( C' a& Y
9 A+ {3 T4 [  I. T6 hif i~=j

d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

k=k+1;
+ U" e! }8 T# M
5 B8 D. @& _& Mend
end
, N  i& G# ^8 M5 f) ]4 \+ u9 o
' C4 j. g( g4 |& d3 yk=1;
end
8 @+ s# p4 U5 yd;

for i=1:n-m+1

/ u1 h3 Q8 `% d8 K[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
6 \" O9 G+ ~1 I3 u4 A2 @% B
cr1(1,i)=l;
end
) j9 F4 X9 j4 d2 c3 Fcr1;

cr1=(1/(n-m))*cr1;
sum1=0;
$ {  H9 z# Q8 M1 k/ hfor i=1:n-m+1

sum1=sum1+cr1(i);

/ s6 T% `& @, K- o# y/ |8 D$ wend
5 \+ ~" P3 M! v5 H  P2 Fend
cr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
$ R% ?7 Q$ J* |* vgn=gn+1;
( [) [; ]+ D5 f* h3 v7 {. bm=m+1;
/ M% V* }# ^1 ?- p- E" c, fend
' F: K! y) |* ]cr;
8 ?& _5 c; H& f& N2 ?. Vshang=-log(cr(1,1)/cr(1,2));
7 z" x5 m6 T  ?% f3 m
[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

其中圆脸为
/ i7 M7 l# c; G8 Q5 e: ) 中间没有空格