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

刹那公子

又是大家开始忙碌着毕业论文的时候了~~
+ k, ~/ T1 k5 w1 x; Q这个是我当时毕业论文的核心内容~~
, J3 g4 d" y. X! i# K1 u2 |% N, m需要做信号非线性的学弟学妹 可以参考下~~
# w4 \5 U' h- N4 F5 N# K7 n/ S' e* Y
1 h) C4 n) Y9 O7 q% w
/ X# L& F8 Q5 l, B) Q3 l0 W
% S; z8 L/ z7 o5 Z4 J' `) {; d
9 G& w& J* d3 J
! q" J! N4 h/ A' d  e! m/ ]

- J" @3 }) a/ H( `7 G, |2 q' G
function [shang]=jss(xdate)
m=2;
4 G' z5 V$ A" r5 {3 Cn=length(xdate);
r=0.2*std(xdate);
cr=[];
gn=1;
gnmax=m;
# K, g# x7 P, v7 g" ?" ^while gn<=gnmax
  X9 G/ y9 g9 `x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
( w/ z: A; N2 r7 A# Xcr1=zeros(1,n-m+1);%存放
k=1;
0 n# C# |+ F  `& I' D, X& U6 M& Xfor i=1:n-m+1
9 y- S& {; q1 p- \- P4 \
: Q% d  P1 B. N, A" l/ wfor j=1:m

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

end
end
x2m;


$ k$ t, M2 A3 B1 g& Qfor i=1:n-m+1

for j=1:n-m+1

$ q) x! K. ^/ sif i~=j
) l$ ^: a  C/ K$ i% L9 w0 Y' z
! G* h5 T* X8 a( J: D; Id(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

k=k+1;

end
end

k=1;
, B/ M+ `) m) q" v) C  ^/ F1 Uend
( |0 B; @! H4 S% Q1 {( dd;
. w+ @. D) C( }0 f5 ?+ H- ?, c
for i=1:n-m+1
4 |5 [1 W7 a& B5 _1 N
0 x0 h- T$ j  I5 s7 I8 f3 K[k,l]=size(find(d(i,<r));%将比R小的个数传送给L

  d( m! V; h5 x, @$ _6 p2 X# tcr1(1,i)=l;
end
cr1;

. w$ H8 g+ z/ n5 @" E2 _/ R* Scr1=(1/(n-m))*cr1;
- `3 U1 v! K8 ^" K, \, vsum1=0;
; M/ r# q% [( G' [6 afor i=1:n-m+1

if cr1(i)~=0

7 o# v3 K2 k% b$ R' Osum1=sum1+log(cr1(i));
, I/ u0 w+ C. Z2 @1 V7 v. k
+ Y: W8 V$ Z" @$ K' [end
+ ]4 z3 f( W5 N3 ?5 C" D* Uend
cr1=1/(n-m+1)*sum1;
" P) U; _8 W* ~cr(1,gn)=cr1;
  n/ o+ F3 {; u$ qgn=gn+1;
! K8 s$ n6 d% H  H7 Nm=m+1;
end
cr;

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

3 N+ O! h- o# a( q- S+ k

- w' ^& f4 \  |$ n4 ~

" H  d7 q* }* C+ U. R1 v7 @function [shang]=ybs(xdate)
m=2;
9 _4 b( I5 x+ i' b9 tn=length(xdate);
r=0.2*std(xdate);
cr=[];
; x1 E! D$ u5 h7 j3 H' p6 C5 Xgn=1;
% x+ y& f1 \. K- y) q; x" Y, Ngnmax=m;
while gn<=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
k=1;
for i=1:n-m+1
8 f. E) h8 g& y- Y, i9 h3 y
for j=1:m

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

end
! C+ ~9 v! \( E8 Iend
& K4 d3 l' o' |( g+ q/ T5 Yx2m;


0 n& _6 @9 z0 j* s1 Jfor i=1:n-m+1

for j=1:n-m+1

3 {# E3 z" H" {- y9 I1 pif i~=j

+ Y; l% K5 r* H9 a6 @d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离

- \/ q, O5 C0 ~7 E8 [1 E3 Xk=k+1;

/ q" p3 _1 M3 Q9 c9 Jend
end
/ A7 I9 R- k) O+ c: O, |
k=1;
' |; J* B# Q. u) `5 M, Rend
d;

for i=1:n-m+1
  {; j5 Q- }: i- }6 w6 W8 ~7 _4 \; j
5 @' Q: I! x$ [8 L1 w* l7 U[k,l]=size(find(d(i,<r));%将比R小的个数传送给L
0 X, c0 |+ }$ M2 k: K- G& F
cr1(1,i)=l;
  j( Y0 {* l8 Jend
: ~: T3 G* n' a) H3 o2 a7 ?8 A; |cr1;

cr1=(1/(n-m))*cr1;
sum1=0;
8 h; K' e  {9 q* ?$ X6 Q" Cfor i=1:n-m+1
* a4 U* M( i0 M+ I% Y5 [
5 g2 p; A5 D$ G, wsum1=sum1+cr1(i);
$ o% ^1 q/ B! S- H' ]
end
7 b% N7 x3 e$ D! X$ G6 C4 Nend
6 Y; L- a4 \8 O! {) X! D0 Fcr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
9 Q$ }2 D( K; [gn=gn+1;
7 F; m, o( Q! p# @. _% am=m+1;
% N3 l# A( S: h8 T( C6 eend
cr;
5 J* T" }& _' V3 A% M! g- `3 |* dshang=-log(cr(1,1)/cr(1,2));

- a# l- i6 L7 W4 ]  `[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]

刹那公子

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