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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
这个是我当时毕业论文的核心内容~~
; r9 g1 i3 @6 X' E1 M* S: l需要做信号非线性的学弟学妹 可以参考下~~


( _& K: V/ ^8 Z; Y4 x: ^3 k
' j6 P7 ~" E3 k& Y; p3 Q& o
% e! b- C- O, H! e6 Y9 V' H
; W. m) }$ u2 N3 Z, \
) z' k1 n* `  M/ o1 j$ L% T1 D

function [shang]=jss(xdate)
m=2;
n=length(xdate);
r=0.2*std(xdate);
2 v% B9 z3 A- L: e, qcr=[];
gn=1;
gnmax=m;
while gn<=gnmax
3 t! q/ Q- g* b9 y$ Vx2m=zeros(n-m+1,m);%存放变换后的向量
2 F  g5 _. B" R) s5 k( T5 ?2 ld=zeros(n-m+1,n-m);% 存放距离结果的矩阵
; T2 H. E) l! l6 y) \/ y3 pcr1=zeros(1,n-m+1);%存放
k=1;
: x  g0 ]+ Q4 c. l& ~! e4 Rfor i=1:n-m+1
$ ~% a. m. E# P+ q/ c
& G' f+ @( r9 P$ o5 x" A! sfor j=1:m
% ?0 ?2 @/ _# U0 J# c( x- ~
1 l3 B' o  {2 g8 T8 c9 Vx2m(i,j)=xdate(i+j-1);

2 j' R+ J, l, u+ t/ {, Zend
2 s- c" _2 R5 e5 [* N: {end
: A* q9 G/ W" F& [0 kx2m;
3 K6 H) w  o! a
& ?1 T& F' A1 R
for i=1:n-m+1
0 n8 N- u+ b: w) v1 Y& T( P0 |
4 X; |- ?7 t" b* z, Ifor j=1:n-m+1

0 K' T; F, g, m! Hif i~=j

7 P$ i- U5 p( p- V* y" L( F1 ]d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
9 F% J$ Z. t) Q1 v$ l
7 C, O7 \: N  n) Rk=k+1;
% h; u! D& A% C# r. ?
  b; y  a1 p. z/ D8 R9 y$ Q- l8 pend
. ^% c2 P$ v8 a, p( _# n& Uend
( M' Z( b0 X9 v* O1 m" U5 m- T
7 ^, k+ ?8 p: Q4 X2 j0 b1 J' Qk=1;
8 n/ z( {, a9 h# E) yend
& H  Q" u4 [/ b( C; G3 _3 b  md;

. B% ^$ Z' h: _- Q# e( n) h! Z  Tfor i=1:n-m+1
# k- |( G& R" O  J' y
# M* }. v: _0 M) c+ u" @, n4 u% }[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
, C- e2 q* b: S- q
' V& M$ J6 l) \0 Rcr1(1,i)=l;
end
9 R) I6 d1 @- s: t6 I# T( Jcr1;

cr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1
  R8 b! W0 |8 R" M1 s1 F. ^
' }- h2 y* _1 Lif cr1(i)~=0
% z" o  O8 P8 O" x
. u; E; b3 t) V! q5 c2 u/ \sum1=sum1+log(cr1(i));

end
end
, I7 v$ Q, L: f+ r8 d$ |& C% B* Pcr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
gn=gn+1;
m=m+1;
& H% V  c, z/ k/ E$ W6 tend
cr;
& J+ `" \2 `+ I$ w& r! B9 w

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

0 l1 g9 A7 K4 v

9 r/ v( l0 C* p

function [shang]=ybs(xdate)
m=2;
n=length(xdate);
r=0.2*std(xdate);
5 v2 [) m4 j' Y5 l! m6 _# E* m* ]cr=[];
gn=1;
gnmax=m;
2 o- \% }% ?1 t5 z! Fwhile gn<=gnmax
# X$ ?$ Y4 D+ ]x2m=zeros(n-m+1,m);%存放变换后的向量
5 K# ~6 v5 @2 U3 r) [7 ed=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
/ E% s, w% O8 R4 u% p. mk=1;
for i=1:n-m+1
' S% C& J' i% K' T0 g7 `( r
for j=1:m

; L4 X/ V) A# }0 P. c$ j: Wx2m(i,j)=xdate(i+j-1);
$ I9 M& G: ^$ t# c- @
9 |8 t- \- d9 l( t9 J- Bend
end
x2m;
" }, Q0 S. ^, y6 G8 W  [

% |/ T" N  ]# S! k$ L; e4 K  Dfor i=1:n-m+1

for j=1:n-m+1

if i~=j
& q1 a( a  r' o! f; D: k# C
d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
5 o$ ^* b, Z6 e
* i  G3 M& K/ P# `k=k+1;
( D' Y. U( p$ Z. i" {1 O# Q: x
end
end
( h# W# z- s8 w( t5 p
k=1;
0 }0 ]7 Z7 b5 n$ k+ Yend
d;
# \- Q) v# I# {: }9 G, s) H' |
for i=1:n-m+1

( l+ `9 c3 l5 W& i[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

cr1(1,i)=l;
2 O, r5 ^( n% ]+ W4 jend
, O  ~. a6 O/ B6 qcr1;

( k# @2 _/ Q! @+ F- ?( X, l' ?7 jcr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1

sum1=sum1+cr1(i);

, p- h: O3 k3 R# [5 ^end
3 J4 O; Z3 v+ Zend
cr1=1/(n-m+1)*sum1;
0 b3 x0 s4 X0 {2 }  @& g# [cr(1,gn)=cr1;
9 _2 z* x8 E, o2 `1 W0 H9 H' Zgn=gn+1;
m=m+1;
5 Q, X' W+ g3 [) g! Dend
cr;
  B, Y. {; t$ d' J* vshang=-log(cr(1,1)/cr(1,2));
2 i; f0 I! m/ h- e6 q  M
[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

其中圆脸为
" [. x8 j: h! t4 N& g7 Q: ) 中间没有空格