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

刹那公子

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

; n! i0 |4 B5 P9 |

9 ]  {$ M; p0 I5 Y1 `9 @2 N
) w; }& h: v5 m



function [shang]=jss(xdate)
" m& \0 k, w- }5 C4 p; Fm=2;
n=length(xdate);
r=0.2*std(xdate);
" E, C" Y* P+ |, {; gcr=[];
gn=1;
/ C& d0 l9 v9 A; }gnmax=m;
9 P8 U8 H8 k3 U3 Nwhile gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
) r* {, b$ ]/ Yk=1;
for i=1:n-m+1

' l; S# ?9 b! o4 O. ^1 Mfor j=1:m

7 i  @; M; e+ C6 @- ?6 zx2m(i,j)=xdate(i+j-1);
6 a1 ~+ H$ U. W
, `5 d% d+ ^6 ], r. Wend
end
  o* T1 |" a0 Mx2m;

- i) f+ q$ V7 d7 G( P4 r) G
' ~( s# f6 Y( ]/ C* ofor i=1:n-m+1

for j=1:n-m+1

7 |( e( X- p9 l/ S  d% ^+ W' @if i~=j

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

' Y& n$ n8 J% }- ]: T. a/ C1 Vk=k+1;
) Y. Q* i; t* [7 o7 K6 ]2 S
4 o3 k# C9 Z7 d$ H3 K& H( |+ F9 Y4 Bend
end

/ r; m; s3 a! S0 [8 K! Qk=1;
end
' |$ O6 m% A3 e8 }0 id;

  h; z1 U( r; K8 Hfor i=1:n-m+1
& m6 o. y" t8 @4 r3 m% X
( y3 Q; U4 k& ~; j. p- G6 ?[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
! L; b" Z8 H# D9 o8 ^6 D  [
  V+ J9 g% a  |, x* e4 Bcr1(1,i)=l;
end
5 H# ?6 x: J: \2 P$ y* F  m9 ncr1;
8 v" z) l  Y; z8 I) Y3 A) q
cr1=(1/(n-m))*cr1;
# D4 S- w2 n4 a/ K- {+ S; Z5 V3 a6 ssum1=0;
for i=1:n-m+1

; }' a4 v. U: B4 Q6 f1 g8 Hif cr1(i)~=0

2 l, O9 o0 w; f; F/ [* Y) \sum1=sum1+log(cr1(i));

" d: B8 X4 ~# Z. G& S" z+ `4 a0 [end
end
5 o- q2 |. W2 y5 mcr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
! K: j# h8 u# p0 Pgn=gn+1;
m=m+1;
' b5 M* K' d1 Xend
/ U/ d$ Y( P1 P  H4 U: icr;
; M( [% x8 |+ }; Q( ~

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

. d' n5 N% T; a8 d" xfunction [shang]=ybs(xdate)
m=2;
n=length(xdate);
r=0.2*std(xdate);
) w# b, S( n3 }) m9 ~cr=[];
gn=1;
gnmax=m;
! w4 b" Z, M7 O  Q" `! {) n. U- Y+ Iwhile gn<=gnmax
3 n5 A1 e/ e1 b6 f% t) N4 c5 _x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
- M' N7 d' C3 E9 Z% q# V+ V3 l2 B7 ycr1=zeros(1,n-m+1);%存放
! ?1 z4 q9 e; f* P* j. B% Sk=1;
for i=1:n-m+1

for j=1:m
, d, v/ N- u3 D. {0 V* E
3 E: j9 r% B! t9 h% v) I( b9 Gx2m(i,j)=xdate(i+j-1);

end
end
) K! h0 G1 h1 L9 f9 o5 Q4 m" Wx2m;
* f0 H2 O5 g7 y
/ e4 o* e, b1 J6 W: R
4 ]* H. ]# W' b- O7 O+ v5 sfor i=1:n-m+1
( Q0 s. F$ {3 o( u* c0 J
for j=1:n-m+1
" |& [9 \+ h) O8 Y" x0 [3 h) \
if i~=j

! J* K% ~& R# u0 O" kd(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

% u/ ~- V# o: q: Wk=k+1;

end
/ Z  l- q6 D: n& O/ Xend
: x7 Z. A* j" s2 D
k=1;
end
d;

for i=1:n-m+1

[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
1 o) Q; q( k7 U" N' M" B. E. T
cr1(1,i)=l;
5 d5 l7 c: t# V% h  h" Xend
! L& s  i( z2 s6 j3 R" Z* U' dcr1;

2 M; g5 m  E7 [! w7 ycr1=(1/(n-m))*cr1;
sum1=0;
2 R. {" s0 N2 |: Z; Y+ `3 nfor i=1:n-m+1

# N% E1 y8 W8 x1 k6 m5 {2 K, Isum1=sum1+cr1(i);

2 [, J. r. ~# e8 Z& uend
end
cr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
gn=gn+1;
1 I+ W2 U- ^# dm=m+1;
end
, @5 @1 a& ]0 Y. {5 L$ j$ [$ h4 j! Mcr;
shang=-log(cr(1,1)/cr(1,2));

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

刹那公子

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