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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
$ ~  T6 ~% ^4 Z4 L7 }! d3 B; y这个是我当时毕业论文的核心内容~~
需要做信号非线性的学弟学妹 可以参考下~~

9 k1 y. u  x' y2 [( c
8 a  Z# c5 H0 n: c3 l1 a
3 c2 y, t5 \# r# [7 d% R6 N6 Q
. I: X& i0 h8 U
; _# f: f1 `* C" F1 S+ e
% Q5 h) f) y# S, m! y' _& z
' B- G4 T4 h  G, v) ?6 X
function [shang]=jss(xdate)
- E* Y5 O/ J( @m=2;
n=length(xdate);
r=0.2*std(xdate);
cr=[];
gn=1;
6 ^* s" {/ G) B, e' ?1 W  o. p/ Ognmax=m;
: W5 k1 k; |. z5 u2 @while gn<=gnmax
6 ^* L# L4 p4 w" n/ A2 S4 |x2m=zeros(n-m+1,m);%存放变换后的向量
* M/ Z9 P; C" Q" s6 O6 d8 Rd=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
k=1;
5 N3 X4 P: h2 V7 p0 N# a/ A. F, Vfor i=1:n-m+1
" Y- P0 Q) N5 s% v4 @
8 ]' z* A2 A& H$ ^. U5 Z% C9 zfor j=1:m
0 a) q6 C7 F2 [/ i9 T# c
x2m(i,j)=xdate(i+j-1);

. Z8 z: {9 }# P9 A. e8 yend
end
9 L3 k% D3 g8 I2 gx2m;
7 I1 j( @5 S% x4 B
) X+ |2 D+ s9 z9 w0 {/ R
# W2 f# n$ h1 e, qfor i=1:n-m+1
8 e: O% j; Z' C7 u
for j=1:n-m+1

if i~=j

0 h7 w- @. T  F) c; z! }( ]d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
' n. ~2 Z& F( K1 J. L& i; ^
: E: e) X/ t/ T1 O% C! E( ck=k+1;
2 R) W3 ?( }9 g
end
/ |/ ~! e6 }% r. ?end

+ N! q  X; F( [) Yk=1;
end
# O! Z5 k; F0 f4 y7 {; h4 Z, Jd;
1 R) }9 t4 ]8 ?% O1 r+ b
for i=1:n-m+1

. [: w; S" f- c6 `% j[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

' \, Z3 |) P5 P& \cr1(1,i)=l;
end
cr1;
% ?2 A0 ]  b/ ^) ~2 J
cr1=(1/(n-m))*cr1;
sum1=0;
/ s# P" g3 I# x" w2 Tfor i=1:n-m+1

if cr1(i)~=0
; S  a. m: y8 Z* a" Y/ u2 o
1 @- U1 c& \! ysum1=sum1+log(cr1(i));

end
end
cr1=1/(n-m+1)*sum1;
$ n& {. o# s6 S: w7 lcr(1,gn)=cr1;
gn=gn+1;
m=m+1;
end
cr;
3 K* d7 r3 l7 H+ ]! E4 o) v6 }

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

. u; D: R) X8 j7 g: u$ P


" V" }5 _/ g# M+ k& T& i1 wfunction [shang]=ybs(xdate)
3 L4 H0 Q. g) f) G. A5 Z" S5 om=2;
! W, Z- U% M2 vn=length(xdate);
, E/ Q7 [5 C+ l* ]/ M7 pr=0.2*std(xdate);
cr=[];
: v3 B4 N! `, L' Rgn=1;
gnmax=m;
" {; g0 N. L- k0 V3 I0 Z# V. {while gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
( q" P  \) Q4 q* T: I+ d+ _- |cr1=zeros(1,n-m+1);%存放
k=1;
for i=1:n-m+1
4 |! d( b  N0 T
for j=1:m

3 S1 U+ }" v* g# hx2m(i,j)=xdate(i+j-1);
0 v& G/ f$ ^, m! w- O& Y
; n* L( J$ }7 J- N1 vend
end
x2m;


4 _6 l7 [: _( B2 ~1 f% K7 A+ c( d- Rfor i=1:n-m+1
( B( v+ I0 \) U* h* m
2 H- H2 |' ?! t& L1 b7 W! h* rfor j=1:n-m+1
5 ?# N, z" t  y& W! W  S
if i~=j

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

% q0 M; T2 r8 q; P0 Ek=k+1;

end
end
! ~$ h; Z, G7 ^, l2 [: Q) \1 ]
8 W  r& ~. N) r# @  wk=1;
end
d;

for i=1:n-m+1
6 o" _6 N$ T  g) g) U. j
  @9 _& I) q- B( {2 n0 M2 }+ U; x[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
9 I6 v: T7 ~- ^& p; F
6 e) d& z, K) M0 J" dcr1(1,i)=l;
; M1 T& p7 p* d. A6 @end
cr1;

% f7 ?4 o% i3 n+ Z) Ecr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1

sum1=sum1+cr1(i);
. s5 S2 `; i) R1 F; k# q
end
end
& P% u' ]: ?+ i) F( \2 {cr1=1/(n-m+1)*sum1;
: w: E4 h. I) h, k. j, u5 Hcr(1,gn)=cr1;
gn=gn+1;
3 L) ~$ ?5 v1 a3 }' E2 zm=m+1;
end
cr;
) }- ^% s, s, Y$ n! s' Ashang=-log(cr(1,1)/cr(1,2));
: a  n( }4 q2 W" i9 H
[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

其中圆脸为
/ q; I+ Y5 @; v! X$ q: O: ) 中间没有空格