| 查看: 1776 | 回复: 5 | |||
[求助]
求信息熵的程序
|
哪位高手知道求混沌信息k熵的程序，也就是给出一段混沌时间序列的数据，求出它的k熵。如果会的话给我一下，小弟将不胜感激。 [ Last edited by 李三石 on 2012-11-5 at 22:11 ] |
回复此楼» 猜你喜欢
面上意见出来了
已经有21人回复
-
能否申诉？
已经有8人回复
-
梦想
已经有5人回复
-
中青基了要发朋友圈吗？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
基金不中，共勉
已经有12人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
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
10