下沙论坛

 找回密码
 注册论坛(EC通行证)

QQ登录

QQ登录

下沙大学生网QQ群8(千人群)
群号:6490324 ,验证:下沙大学生网。
用手机发布本地信息严禁群发,各种宣传贴请发表在下沙信息版块有问必答,欢迎提问 提升会员等级,助你宣传
新会员必读 大学生的论坛下沙新生必读下沙币获得方法及使用
查看: 6154|回复: 3
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1
发表于 2008-12-27 22:27:54 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
又是大家开始忙碌着毕业论文的时候了~~
% `, h2 ~2 o' L0 f9 _- s这个是我当时毕业论文的核心内容~~. z, a9 {6 F( c" z" C! m
需要做信号非线性的学弟学妹 可以参考下~~
- u7 x3 e4 q% t: B" N8 y9 ]2 m( B5 V, Z3 q

  }$ A$ g. l* k( v7 ^5 a- ^7 @2 E, A( L; a7 S+ r

8 a- ~1 P0 n) V' A
4 {# k' d+ g: t- Z) n0 Q# M  Y! o; i9 {& @# X5 W- |& v* K, A
* e; m% p4 v5 X: c

" I* Y- q; M  L5 Q8 O6 l. v+ W( yfunction [shang]=jss(xdate)( n' }& [- J, E5 L# L
m=2;
. K. o1 y2 p$ }& _0 bn=length(xdate);
) P2 ?* R: K' l& V& Mr=0.2*std(xdate);
6 U  H: G/ l: W4 H' O' zcr=[];0 q% \' l! a4 n/ |1 F5 A
gn=1;
! v; q/ g& E4 v1 G  M' r! Hgnmax=m;
8 q, `( b* z! `2 c, |while gn<=gnmax
+ f! h+ N2 E" N( P7 Dx2m=zeros(n-m+1,m);%存放变换后的向量9 ^9 @3 l. m( c+ T
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵/ ]7 t% x( }2 L; b: [; ^# b. d
cr1=zeros(1,n-m+1);%存放
8 T8 K2 J, Z( O- nk=1;
  t4 D; p: z9 L8 x+ F* o. }, sfor i=1:n-m+1
0 U5 {( w  }8 k9 t5 j; e# o; B
3 F0 n! c, h  Y. e  |4 }5 N( sfor j=1:m

% J7 x& `9 t  l3 d& E" g6 E- F, I: R% U
x2m(i,j)=xdate(i+j-1);
) R% F* y' r' [. h- C, m

& m. D7 f4 i# e- L9 i3 x- x# _end
& W! h: H, W  O% |
end
3 t% b3 d7 M5 [x2m;7 H+ ?8 c/ t" y$ r( f7 E
2 U1 z% f9 f& p: @0 R. H

4 e, U% S9 U8 o5 P$ Mfor i=1:n-m+15 v" E2 K4 b0 Y) j7 G

6 S  R, J* U$ u6 J/ Z" Jfor j=1:n-m+1
! ]2 Y7 J2 N, |! {. V# A

  x& s% f6 q. W: n& _3 W9 iif i~=j
' h" a) f% }8 ~( P8 V+ J
, a( Q! F* |7 s8 ~2 M
d(i,k)=max(abs(x2m(i,-x2m(j,));%
计算各个元素和响应元素的距离
1 A0 [' C8 Y8 ^, b
6 o/ D' x: c' L. _- Ik=k+1;

: T7 v. b6 c' \/ ^- \, J. v' a" `
3 W6 [: i; X) f; B8 N$ ]end
, v7 }5 E1 W/ t9 z' a- P2 p9 _# G
end* Z8 X! \* J6 O1 q0 d, i3 P- d+ X- M& x

1 G$ i6 L' ]; C3 E9 Vk=1;

8 h- O! h3 q/ K5 h$ |( ~3 P9 dend
/ z8 l  C! A& s& G$ {" md;+ a# r" K& x  y1 M

5 L. v# J. R! k- vfor i=1:n-m+1
' u; X, L4 q$ @, t! h- u. m
  |# e- G$ i; I: U$ m5 c[k,l]=size(find(d(i,<r));%
将比R小的个数传送给L
9 j3 H% N4 Q1 h
' U+ S, K2 d' ccr1(1,i)=l;

4 `( A1 ?' b; {end
; m7 G+ R  T$ J$ p$ D- ccr1;: f# N) Q$ v9 L" g1 t

) y. g- E: E* n5 t# t2 m2 Kcr1=(1/(n-m))*cr1;! d7 b9 ]3 |* b* v7 J
sum1=0;9 K" R4 ]' q" f( a+ e  q
for i=1:n-m+15 P4 C+ h1 C- c8 Q

* b5 v2 {8 _5 w) ^$ g8 `if cr1(i)~=0

$ I( H9 g. t) W* R9 @4 x4 k" [0 j! {# B
sum1=sum1+log(cr1(i));

$ G+ b3 X, C' w+ ~! V, A: Z
% Q, a3 e9 F7 C0 s5 d! Vend

( ~6 f1 g# b% j, ?+ ~end
! L( C$ e3 H! I& d6 J" h  n& scr1=1/(n-m+1)*sum1;! W, P* R) l* W5 ~6 p( `/ x
cr(1,gn)=cr1;% q' P( s8 l- f3 E% h
gn=gn+1;
6 r) J: E' P8 l2 G0 \m=m+1;
  \: g: }* K' ~# G' C2 iend% \' v/ X$ P9 @
cr;
" M% L5 R: P* ~. p
shang=cr(1,1)-cr(1,2);
; v. A6 S  q' m0 Q  d9 m% m
/ E5 s+ @: m3 Q7 N: T

5 ~% a! P+ _5 Y& \  E0 S! @! m" r% i8 u8 T5 e: k
$ Y" X2 A* b$ \( u& |2 k
function [shang]=ybs(xdate)8 `! J+ a( p* U7 T
m=2;1 q+ M0 M" j4 g! g* D8 R$ u
n=length(xdate);
% k# i1 ?1 ?0 A4 ?r=0.2*std(xdate);
1 V# q' ]: X; M; F' I0 y7 P7 Ucr=[];: {4 D9 i0 y: S  Z5 P& Z9 `/ e. J
gn=1;6 j8 N# q, p% Y# V: N
gnmax=m;
$ o" x! X9 G  cwhile gn<=gnmax$ p9 V% K$ p+ D7 z2 K, ]
x2m=zeros(n-m+1,m);%存放变换后的向量
, L$ p) r" t: A2 Vd=zeros(n-m+1,n-m);% 存放距离结果的矩阵
$ V: @/ q; p6 o& ?cr1=zeros(1,n-m+1);%存放0 h* q% w* Q$ p" c5 P* g" q6 `9 y
k=1;
' `& {1 {8 `& X: dfor i=1:n-m+1; i% c& ~) ?* [7 \: d

. \4 G8 H8 v& N: H) f- Pfor j=1:m
) d; @# o8 j8 ]7 O

) b7 F) Z6 C( n3 ^; K3 Tx2m(i,j)=xdate(i+j-1);
0 z- A0 c( @) R' D$ a8 Z- B8 |
- P% U8 T9 u% o1 _$ h/ |7 J% K
end

4 {% b4 V2 h  S6 Y# p! Tend
1 i5 Y- d; x5 A2 i' T0 h; ]x2m;' E! X6 |& w$ ^4 B' D' A
. l/ q5 }# K& I

2 {  u* R2 E8 e; [for i=1:n-m+10 r4 f- A, ]  h. h; S8 }; V  N/ {. E

0 p/ z6 E' Z5 h$ Kfor j=1:n-m+1
2 M% X) L2 `6 t" u; m6 ]( v/ M

: l) f' e( |+ ?8 y0 _4 l- Bif i~=j

1 e; E1 B: Q- h3 M* X
6 n+ S  l4 ~, }6 l/ Jd(i,k)=max(abs(x2m(i,-x2m(j,));%
计算各个元素和响应元素的距离
% e1 ~# c  `. R# B, u  A1 J+ @( K( J# F  N
k=k+1;
$ ?0 A; ^8 p$ Q2 [
" x  @( I* T7 ^$ P/ }
end

/ h8 j5 p7 Y6 a8 ]* Bend
+ u  [* s  m7 K" d8 ]' \
0 a' \) l7 r5 k. m9 qk=1;
  {0 a: s* |6 f1 C1 N# G( ]4 H
end3 x( H3 r: p" X0 J  ]
d;
; {5 d+ b- n# b1 u) [
" ]. M0 {! k/ e0 \6 n3 j/ h
for i=1:n-m+1
" S0 v+ }- d% V$ t( l5 J+ T/ i6 W
[k,l]=size(find(d(i,<r));%
将比R小的个数传送给L
+ V5 b4 ]' w$ ^" p
& G, ^/ ?' g% o* O) |cr1(1,i)=l;
' S. t6 W5 a" _
end1 I* o" _3 D( g0 t5 l1 y. w! S
cr1;
0 m4 B* |, A/ c4 a' ^4 m3 {; |9 c
' o3 a. ]0 f; \7 Jcr1=(1/(n-m))*cr1;
( O( w3 }; V" @0 W# psum1=0;
2 `6 ~4 p% c1 `) @) kfor i=1:n-m+1
7 Z2 G: z7 |+ p0 D) l/ k. }3 |9 A* b/ E5 N: J
sum1=sum1+cr1(i);
! U# m7 b- z5 \2 y6 F

9 V0 R3 h" R6 b, f6 l/ eend

8 N9 t' g( }3 B  Bend
) f/ s, D1 ^& Y1 B% Scr1=1/(n-m+1)*sum1;
4 F8 x2 ^3 Z6 g1 p$ O6 hcr(1,gn)=cr1;0 K/ P, `: c2 A) w5 o1 s
gn=gn+1;
) Q+ s  v5 g% Um=m+1;, Y! p+ E# x  t& W
end
8 ]' [$ v1 x- s) [cr;4 S* e2 V& X8 E% X7 y6 }
shang=-log(cr(1,1)/cr(1,2));
6 e! m5 y5 y) }) Q* K$ L
1 y# M) D0 J# I: n, _# Q[ 本帖最后由 刹那公子 于 2008-12-27 22:30 编辑 ]
分享到:  QQ好友和群QQ好友和群 QQ空间QQ空间 腾讯微博腾讯微博 腾讯朋友腾讯朋友
收藏收藏 分享分享 顶 踩

该用户从未签到

2
 楼主| 发表于 2008-12-27 22:31:56 | 只看该作者
其中圆脸为% Q. k6 T' M" t. p
: ) 中间没有空格
回复 支持 反对

使用道具 举报

本版积分规则

关闭

下沙大学生网推荐上一条 /1 下一条

快速回复 返回顶部 返回列表