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

刹那公子

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





: {$ @0 v' r$ p) S; V$ P
7 f9 }& x6 M. c  z

9 O6 d& J2 X+ C3 C: @function [shang]=jss(xdate)
$ V# x6 F' j, }8 tm=2;
* f+ U1 }3 \6 A3 rn=length(xdate);
3 X' c4 ^& M2 {r=0.2*std(xdate);
cr=[];
gn=1;
& c) f$ Z. P" v" Ugnmax=m;
. A, o6 M: _. P" n! {while gn<=gnmax
8 `- n5 a/ B  ^" d& Zx2m=zeros(n-m+1,m);%存放变换后的向量
0 G* G) s  b4 k7 o# Gd=zeros(n-m+1,n-m);% 存放距离结果的矩阵
+ z; N& D) i" m  g) Z# K4 |cr1=zeros(1,n-m+1);%存放
; c$ {2 D7 m0 B9 h# qk=1;
6 K7 ]- z5 h. |+ Y+ afor i=1:n-m+1
5 a( F3 a# {3 ^3 l) X0 s% Z4 m
9 e- ?+ C6 e7 z/ }for j=1:m

6 b. }( K) R8 I" o3 L+ T# ]% cx2m(i,j)=xdate(i+j-1);
: U: U/ P& s8 t$ V# K* b5 ]! |
end
end
" |$ w! E4 @( K6 ax2m;
8 l1 i6 q% C  \5 f9 H

% j9 ~  e" G* P. T- ?for i=1:n-m+1
. L; v  X) j  E; ~$ q3 m* w
' _  `1 u8 K5 B8 }4 gfor j=1:n-m+1
. F4 ~$ P2 p# K; ?/ X1 y" V9 ^8 A
, B3 T% @" y; q4 K( F6 tif i~=j
. |# [7 D& u$ q  ?
1 ]4 v3 M) j% u& T8 s0 ~% ad(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
- q8 o, Q. R" l/ Z; }8 J5 w
k=k+1;
+ c4 \6 N* @! }; g6 F% O
* q7 U6 R  d" p: Fend
end
8 a7 k3 o9 D$ [6 z0 k5 E0 t
k=1;
5 _) ^# R1 M8 @& x- z6 gend
; H4 T) V& n: W; a% e9 {d;

5 L& a& Q! B6 Hfor i=1:n-m+1

1 e8 X. K6 S' z* T% ~7 c7 ?- _[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
5 `  C8 U- [9 f6 w: Y7 i% Z) E
0 n+ d$ T& ]7 {cr1(1,i)=l;
end
cr1;

2 c7 @4 h$ L: O- s$ @cr1=(1/(n-m))*cr1;
* q+ Q+ m0 B  c9 X  Usum1=0;
for i=1:n-m+1
  c; T- A3 e! A7 {, H- Q
4 g/ r$ |  Q4 Q6 L  k9 x. Yif cr1(i)~=0
1 u8 Z( J% v, w% h4 v
' O: P- T6 M6 K# p/ Usum1=sum1+log(cr1(i));
8 T- A( o. O- U3 L% w0 j7 z
) o9 K; W0 c3 n/ i5 N$ L1 Dend
; L' ?+ O# N+ kend
' r' A+ V; E# ocr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
0 Z( N5 [) k3 N+ f3 h$ u, U& @gn=gn+1;
m=m+1;
end
& I4 t1 s0 m3 }: O+ O- |cr;
/ _; H+ r" k5 h. w4 j1 X  O: B

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

0 t& `5 N* u9 A: y, S! Z' z
/ R3 U+ C; w2 O4 u

function [shang]=ybs(xdate)
7 K, B. C5 s6 n% tm=2;
n=length(xdate);
r=0.2*std(xdate);
; [% _" B4 R6 g0 B6 q" \cr=[];
( y& U4 I4 o- j, [gn=1;
' s' m% p, V3 G$ `2 Ognmax=m;
$ a5 s; H$ k0 u/ B% Iwhile gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
- d0 Y8 f( I8 @" Y8 Acr1=zeros(1,n-m+1);%存放
k=1;
for i=1:n-m+1

, m1 |! l) e1 K) p* W" f5 Rfor j=1:m

x2m(i,j)=xdate(i+j-1);

4 V& V% l4 g& |! o- i4 Pend
9 E+ \# T/ s2 h. }5 I6 G- \! Y/ Jend
* [4 Y% Z$ {3 o; l8 [x2m;


for i=1:n-m+1

for j=1:n-m+1

if i~=j
4 Z; ]: i% d! c
) g3 g9 N9 \0 ~! w* T1 v. H: @d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

k=k+1;
( @8 x* B# y  ?
end
8 A) Z, C" Y( Y3 Nend

k=1;
end
d;

for i=1:n-m+1

[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

. A" r. G0 l6 `9 {cr1(1,i)=l;
end
cr1;
( D" Q% O( F) S9 }7 x9 k- h! N% N8 {, o
8 V+ u8 B. }8 J* \, fcr1=(1/(n-m))*cr1;
2 V# W0 P" J% c5 p! Qsum1=0;
" @7 }' k7 ], g" qfor i=1:n-m+1

sum1=sum1+cr1(i);
$ z+ }2 O  e% o. ?( Q  h
end
7 J" L7 M  U7 D1 y# tend
cr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
gn=gn+1;
m=m+1;
end
: x4 b, B& Z3 }. |3 G; J* Scr;
1 g# ]+ U; ~0 D* D( @6 y7 _% Xshang=-log(cr(1,1)/cr(1,2));

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

刹那公子

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