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

刹那公子

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


) b8 m  l3 f; H' i% l) g

$ E1 E7 D7 F1 L7 P/ G' N" M( U2 i


* F- I* L$ J# P' Q: q
function [shang]=jss(xdate)
( [% m+ \- s- x1 K! B2 B7 H  {% Rm=2;
  |8 q& v& C1 y9 M0 y& a. Tn=length(xdate);
r=0.2*std(xdate);
- @9 ]- ?5 x) X2 x2 m- v6 ycr=[];
gn=1;
5 S& F3 W. K& a; A" q8 tgnmax=m;
5 C6 d8 H9 x0 v+ ]- h$ iwhile gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
8 J' Z. Y) B$ ~8 q( Z/ D) _4 Icr1=zeros(1,n-m+1);%存放
* k8 V) q2 \' H5 \  W! lk=1;
for i=1:n-m+1

for j=1:m
. A) _/ ]9 ]" V. H% P; n3 Z! F* H7 v
x2m(i,j)=xdate(i+j-1);
/ O* `: l4 O+ j  O# L0 e
end
5 a# K# L8 w) \( f, f2 jend
& o' v! A/ v. kx2m;


! Z! k! I! ^5 r9 rfor i=1:n-m+1

for j=1:n-m+1
# T  D) _: k) e$ _, r3 l
/ c  U% o$ n6 O, _4 \if i~=j
) h- E  o, O5 P' Y$ x  L
$ |+ s6 b# Q' M0 e8 id(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
/ H9 J0 b* T1 ~7 H6 z1 ]! {5 J, }4 w
  E2 Y+ W' P3 ^/ c! hk=k+1;

6 b7 ]% m& Q) Aend
end
" @# c/ Q4 W: P
k=1;
# D0 W0 b" k- F1 k  s$ ]end
6 k) y( f1 o; r: Y+ V9 x3 ?1 id;
3 N+ r' a  a9 \+ O; t! Y3 ]5 D
% Q3 d/ d1 X# P8 w6 ufor i=1:n-m+1

8 l; F  N7 `7 C' [& I0 C[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

cr1(1,i)=l;
end
/ U5 M7 z6 _. W3 l3 m, Rcr1;

8 w; `8 M2 m* e; Ucr1=(1/(n-m))*cr1;
2 I: z9 {+ Z, }7 I$ A' zsum1=0;
for i=1:n-m+1
2 y7 l" F3 {8 W& @; d% N; F) Z
& x" `, m  I5 _if cr1(i)~=0

( m( n: E$ C6 r/ s% \$ A! l9 gsum1=sum1+log(cr1(i));

end
end
+ E, J: w8 x9 j" s: ~cr1=1/(n-m+1)*sum1;
" z* H/ {' i: J, o9 Y* ncr(1,gn)=cr1;
gn=gn+1;
m=m+1;
3 w% D% i0 }3 O; ^/ A# Vend
cr;
! L$ D  ?7 q( N( p8 L4 j/ q) Y

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

function [shang]=ybs(xdate)
0 i6 K7 b9 d6 Sm=2;
n=length(xdate);
# m* t6 h- ]+ N) T6 ~) yr=0.2*std(xdate);
: {; X6 ~& x6 d/ z* ycr=[];
gn=1;
gnmax=m;
" Y* P7 i2 P3 z/ Qwhile gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
: c; `2 l4 M  g  q8 vd=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
k=1;
+ H5 u8 R0 f# P7 y6 f2 gfor i=1:n-m+1

for j=1:m

x2m(i,j)=xdate(i+j-1);
: y2 ^7 ^6 v' }
" d0 E! D8 v2 j; a1 T6 _end
/ r; ]# ^: @1 Z( Mend
x2m;


1 `0 M2 h* ^  Y* d8 z! Cfor i=1:n-m+1
6 }; f4 v6 b) F- D, J4 }
/ S; p1 e; h$ H& h$ Tfor j=1:n-m+1
5 r) e  W1 J3 K
if i~=j

7 C0 c6 y& K% W5 [8 Z# `d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

k=k+1;
( x+ i8 y" ^+ U
( z8 k; Y& ^) o; y/ Y8 W9 Dend
end

0 ^7 o, z1 Y( M! s2 ~+ g1 a; Nk=1;
# A7 E* V2 C1 W3 _- }8 Aend
d;
% _$ v& u7 ^, ^. B9 q2 w! p
for i=1:n-m+1

[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
, Q2 z# L$ C" ^! ^7 }* A
9 o, A7 c# k( U2 e' ?  Fcr1(1,i)=l;
end
& U* K: C* z& s/ O# V' |cr1;
" R' D7 v6 ~: I$ F$ ]
# _: {% Q/ h, wcr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1

sum1=sum1+cr1(i);

end
end
cr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
6 O/ v1 M) D* O7 M, j: Sgn=gn+1;
m=m+1;
  f% z7 W4 q! @* [' t1 vend
cr;
shang=-log(cr(1,1)/cr(1,2));
' p5 S; `3 [, q
[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

其中圆脸为
  [0 x- l' B5 x; K) S4 w: ) 中间没有空格