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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
5 P6 ?5 F7 @3 ^' R' g这个是我当时毕业论文的核心内容~~
( X9 B1 v. C3 T$ Z! P1 x$ d$ G$ d需要做信号非线性的学弟学妹 可以参考下~~
) k* v% [$ ~- Q9 ?* l+ h& e0 O

$ J# m* k; n+ h; q7 x3 H0 {




  a- O6 }" g: ^" {! Q0 D
function [shang]=jss(xdate)
m=2;
: H. Y/ N' ~) T1 W+ yn=length(xdate);
r=0.2*std(xdate);
) O4 k/ B4 B5 z% e, t% M9 p6 d2 Fcr=[];
# O6 [% ^$ t) w+ bgn=1;
gnmax=m;
; B0 E3 a6 f. @while gn<=gnmax
8 E2 K" v; Q3 Yx2m=zeros(n-m+1,m);%存放变换后的向量
* a' {: Z6 o$ Y0 [d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
k=1;
, U* q1 p% ~/ J6 cfor i=1:n-m+1

for j=1:m

, B  H5 V4 @6 D7 e% j1 lx2m(i,j)=xdate(i+j-1);

end
* Q% g6 z3 H3 H0 x& C) uend
x2m;
2 a9 e: ~/ H7 ^! g

3 t# M1 Z# t/ i& n" Y4 K  Jfor i=1:n-m+1
, D6 E& D! a: e9 q
/ T2 L* n( E( F# w7 F/ L) Xfor j=1:n-m+1
7 r7 M4 j2 I/ x
! Y- i& v8 \7 qif i~=j
1 p; ?. V5 A0 l% p* S, Z+ q- t9 G2 ]
; B; w8 _9 }# E) P9 r' a" H; g7 T' v* A% @d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

# A9 Q- E& [; Mk=k+1;
; H- E7 \8 V$ Q! A( `& d
6 M/ ^  A( @6 z" h: d; qend
, I# L* m  h8 ?) |end
6 o& a# Z7 }7 W# Q! M
) ~( \2 X* ?2 u: Lk=1;
end
d;
/ z) M, V& i" H7 m, c
for i=1:n-m+1
! i" \. `2 p' \1 O4 z7 w
% k: m0 l, l$ B# `$ e& K1 r[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
# s( F. m1 i, e3 o1 W9 s! d
( z! N& {% C+ L/ Z, a  S& `cr1(1,i)=l;
end
5 u" l) O3 H2 C! ]! ?7 i+ s" s/ n4 dcr1;
# F( k+ n' `% A- R/ G* g! Y
3 `' N' A' d  Z0 tcr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1

if cr1(i)~=0
' }( m- z- N/ y: m9 M6 d
sum1=sum1+log(cr1(i));

end
end
cr1=1/(n-m+1)*sum1;
: U2 e6 y4 b8 Z# l2 ]cr(1,gn)=cr1;
3 @. s5 E3 w# L. N/ Cgn=gn+1;
m=m+1;
; R8 X4 Z6 I' [9 Iend
, @& N* I8 ~( E$ j/ Z' O+ a# Fcr;

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

* J4 B; }3 M( m" a
' l. U2 a& }5 hfunction [shang]=ybs(xdate)
. E1 N1 B2 Q/ u. j; m0 xm=2;
n=length(xdate);
r=0.2*std(xdate);
cr=[];
gn=1;
gnmax=m;
while gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
  e$ w! U0 \9 _( k) {+ Ad=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
k=1;
for i=1:n-m+1

3 z! L& l0 n3 L- P( T) wfor j=1:m

2 [8 _5 c: z7 r; Sx2m(i,j)=xdate(i+j-1);

end
8 A  F8 P& w8 rend
: ~. C4 s. ]1 E$ x8 s/ x# S, V* ux2m;

/ b3 Z4 i' r" V* j& j# m
for i=1:n-m+1

for j=1:n-m+1

if i~=j

d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
9 @6 I2 Q! x$ P3 m7 i) c
" D& e  @4 D3 _3 Z8 y/ C3 l- Dk=k+1;
. o( c( H4 z0 y( t7 Z) `
4 n: L! J8 X$ w& c, tend
end

k=1;
2 y6 Z/ ^% W4 b5 H; G$ Tend
; O$ M  N+ V) h' \9 ?d;

0 d- I" X2 u8 x' S3 wfor i=1:n-m+1

7 c6 `1 |' u* C, n0 L[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
' V4 g# g& z# k% f
cr1(1,i)=l;
end
cr1;
+ Z1 Z% ?0 W/ o" {# ]! x1 M; g
cr1=(1/(n-m))*cr1;
( j4 @: i; y* y- r1 ?sum1=0;
2 Y0 ~& W7 x' L) Cfor i=1:n-m+1
6 q0 R. g9 s  Y% P: W: ]
1 r. p  T$ a/ j9 \: [8 bsum1=sum1+cr1(i);
0 s. E- y4 t; B2 Z
+ @  C, ]: n( n6 u) G' [2 b  Xend
end
; L4 m5 n4 c8 Zcr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
! s* C) N/ I, _8 ngn=gn+1;
m=m+1;
5 u( N6 ^' _) `) ^  J; [1 U, Iend
cr;
shang=-log(cr(1,1)/cr(1,2));

/ v( k0 Y/ D) e5 F6 o[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

其中圆脸为
8 R6 p8 C8 M( i: U0 N: ) 中间没有空格