| 查看: 1780 | 回复: 5 | ||
[求助]
求信息熵的程序
|
哪位高手知道求混沌信息k熵的程序，也就是给出一段混沌时间序列的数据，求出它的k熵。如果会的话给我一下，小弟将不胜感激。 [ Last edited by 李三石 on 2012-11-5 at 22:11 ] |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
bohrzqc
新虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1201.9
- 散金: 1026
- 红花: 4
- 帖子: 731
- 在线: 878.8小时
- 虫号: 705017
- 注册: 2009-02-20
- 性别: GG
- 专业: 激光
2楼2012-11-06 08:42:52
3楼2012-11-06 11:53:53
bohrzqc
新虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1201.9
- 散金: 1026
- 红花: 4
- 帖子: 731
- 在线: 878.8小时
- 虫号: 705017
- 注册: 2009-02-20
- 性别: GG
- 专业: 激光
4楼2012-11-06 17:23:56
5楼2012-11-06 21:04:50
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
李三石: 金币+50 2013-10-23 18:49:58
nono2009: 请认真阅读此帖，http://emuch.net/bbs/viewthread.php?tid=3711903。 并与两日内解释您和ID 李三石 之间的关系。 2013-10-24 12:50:04
nono2009: 金币-50, 马甲存档, 马甲转移金币行为：http://emuch.net/bbs/viewthread.php?tid=1996653 2013-10-28 14:43:23
nono2009: 金币-40, 扣除转移金币 2013-10-28 14:43:37
李三石: 金币+50 2013-10-23 18:49:58
nono2009: 请认真阅读此帖，http://emuch.net/bbs/viewthread.php?tid=3711903。 并与两日内解释您和ID 李三石 之间的关系。 2013-10-24 12:50:04
nono2009: 金币-50, 马甲存档, 马甲转移金币行为：http://emuch.net/bbs/viewthread.php?tid=1996653 2013-10-28 14:43:23
nono2009: 金币-40, 扣除转移金币 2013-10-28 14:43:37
|
%———————求一维离散序列信息熵matlab代码 function Hx=yyshang(y,duan) %不以原信号为参考的时间域的信号熵 %输入：maxf:原信号的能量谱中能量最大的点 %y:待求信息熵的序列 %duan:待求信息熵的序列要被分块的块数 %Hx:y的信息熵 %duan=10;%将序列按duan数等分，如果duan=10,就将序列分为10等份 x_min=min(y); x_max=max(y); maxf(1)=abs(x_max-x_min); maxf(2)=x_min; duan_t=1.0/duan; jiange=maxf(1)*duan_t; % for i=1:10 % pnum(i)=length(find((y_p>=(i-1)*jiange)&(y_p<i*jiange))); % end pnum(1)=length(find(y<maxf(2)+jiange)); for i=2:duan-1 pnum(i)=length(find((y>=maxf(2)+(i-1)*jiange)&(y<maxf(2)+i*jiange))); end pnum(duan)=length(find(y>=maxf(2)+(duan-1)*jiange)); %sum(pnum) ppnum=pnum/sum(pnum);%每段出现的概率 %sum(ppnum) Hx=0; for i=1:duan if ppnum(i)==0 Hi=0; else Hi=-ppnum(i)*log2(ppnum(i)); end Hx=Hx+Hi; end end %---------------- |
6楼2013-10-23 18:48:30
20