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

刹那公子

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

& O! m8 F2 H. L8 J) A- r
* L2 r( l  T  m

: V3 R+ P$ G9 }! T1 ~6 ~9 S- L# @' F
- Y7 k* `0 ]; R6 p6 ?# C% s6 e) ]
4 s( i% Z" E: f

function [shang]=jss(xdate)
3 q2 h0 [3 [% n6 @! dm=2;
/ a0 v( A/ F' kn=length(xdate);
7 X2 x5 m, X3 \8 E& Y. ^r=0.2*std(xdate);
cr=[];
1 T6 [$ X- A, Y# c$ _/ zgn=1;
, J, K2 T& F% _$ b2 sgnmax=m;
while gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
( k& u) s% s9 T2 X3 a, ~cr1=zeros(1,n-m+1);%存放
k=1;
for i=1:n-m+1
: D1 @. L6 g# `4 p. d
for j=1:m

x2m(i,j)=xdate(i+j-1);
; O  N# ^4 H6 P0 J9 v% I( S1 ~6 X4 q
end
# ]; d% }8 g& x, {/ s6 Oend
  c( N& D4 x) P' ?x2m;


for i=1:n-m+1
! S1 e+ u+ k9 v9 R+ B
for j=1:n-m+1

) U* I/ X: T2 n0 g9 `0 P5 R5 W9 zif i~=j
' F( x  i) E# t; P: ^  b" m
d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

% V$ t% M1 [5 G- g' Z$ _k=k+1;
0 z% p- l5 f- h& S# `
end
1 B* r1 ~! x! \/ ?3 r; \end

k=1;
$ J$ h0 R5 f* T& G9 W; T6 a+ S+ I+ ]end
d;
+ k* X$ Z" |$ u* C. m
for i=1:n-m+1

4 Q/ L3 j# T7 t+ f; i[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
- F- z5 f7 q! t. y1 j4 R7 H
cr1(1,i)=l;
end
+ T# ~8 V& j* z; f- u) _! acr1;
2 o7 H! N/ |8 {- b
cr1=(1/(n-m))*cr1;
$ z8 L% R. A7 Z! l: }; ~sum1=0;
) s9 i- k7 M" L7 j! Ofor i=1:n-m+1
9 G9 d4 G: J3 @$ E, d: p# j: z
, |! N5 \3 i& U6 Y9 Hif cr1(i)~=0

7 z" x% B2 _( L; D4 G3 ~! xsum1=sum1+log(cr1(i));
; h3 q) E( @$ M, _- f, ?
end
end
# d- `, m4 C  z* V4 c, acr1=1/(n-m+1)*sum1;
! ^1 c, M0 i; K7 j- e4 j: T# L. hcr(1,gn)=cr1;
8 J5 u9 s5 Q2 q4 Q1 Ugn=gn+1;
m=m+1;
end
cr;

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

% ?5 \. M( @% V0 c


) [, l' i* |, O0 T' o4 k
function [shang]=ybs(xdate)
% E" Z& Q( {& g2 H4 C( Im=2;
: w% S* r" D: e  hn=length(xdate);
# o0 _$ |; T& W- y4 K! Kr=0.2*std(xdate);
cr=[];
9 {6 P9 _9 @7 hgn=1;
gnmax=m;
while gn<=gnmax
8 _4 S6 I5 i( _8 ^" px2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
% a9 d# l' w% c: x8 ucr1=zeros(1,n-m+1);%存放
k=1;
1 y: O5 r/ c1 V/ s# tfor i=1:n-m+1

for j=1:m

# O7 B; G) _+ @2 }# dx2m(i,j)=xdate(i+j-1);
6 X. ?- F/ h1 v" {  {' j0 `# U
5 O1 ^7 D+ T$ Iend
$ c0 X9 n( S2 Send
) ~) j. n3 T3 W( ~) yx2m;

4 o8 c$ Y# `8 n$ O" d' q1 v
& X9 I) X# {1 `" K  pfor i=1:n-m+1

& w+ M0 D' H, X0 l# gfor j=1:n-m+1
- Z& ~/ l" g9 T; U, k; |
$ |+ J" M2 h' J6 q% R# g" u+ x1 l. bif i~=j
5 M  J$ ^/ @& P7 o  p) ]
& O" H: A% D& g9 E. @6 jd(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
, y4 n" K7 p0 ^+ E9 P
k=k+1;
2 ^# l9 }8 R: y3 n
end
end
6 m& H; t2 G4 {# l8 I
k=1;
end
3 i" f0 r7 \9 c2 Y5 F' K! ]2 R, td;
5 A$ p# |8 R" \( b; }' O! ?
for i=1:n-m+1
& n' U# I' a4 Z  s' e
[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
9 z8 w0 C8 Y- d8 F- d% ^! ?- o. |
cr1(1,i)=l;
$ L/ U0 o% K% w  L) aend
! j! T- j9 p/ Rcr1;
+ R2 T  D6 c. {
cr1=(1/(n-m))*cr1;
  v7 m3 L6 y! Lsum1=0;
for i=1:n-m+1
, ~  b: C9 W/ k- F- C  M7 n" m
7 K; g. _7 r+ P3 Wsum1=sum1+cr1(i);

end
end
0 v$ ?+ w: C2 Mcr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
gn=gn+1;
, `) H  k1 J! ]2 W; C% ]: Jm=m+1;
end
1 i7 Z8 ]1 b' j! |8 k  r# \0 Ycr;
shang=-log(cr(1,1)/cr(1,2));

[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

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