| 查看: 826 | 回复: 2 | |||
[求助]
求混沌或者噪声熵的程序
|
大家好，谁有求混沌或者噪声熵的程序，c或者matlab都可以。跪求 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有34人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有78人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有12人回复
-
FileCode能看出啥？
已经有30人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
-
关于Filecode分析方法
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
李三石: 金币+100, ★★★★★最佳答案, 20 2013-10-23 18:41:30
ben_ladeng: 请认真阅读此帖，http://emuch.net/bbs/viewthread.php?tid=3711903。 并与两日内解释您和ID 李三石 之间的关系。 2013-10-24 09:49:55
nono2009: 金币-50, 马甲存档, 马甲转移金币行为：http://emuch.net/bbs/viewthread.php?tid=1996653 2013-10-28 14:46:00
nono2009: 金币-80, 扣除转移金币 2013-10-28 14:46:14
李三石: 金币+100, ★★★★★最佳答案, 20 2013-10-23 18:41:30
ben_ladeng: 请认真阅读此帖，http://emuch.net/bbs/viewthread.php?tid=3711903。 并与两日内解释您和ID 李三石 之间的关系。 2013-10-24 09:49:55
nono2009: 金币-50, 马甲存档, 马甲转移金币行为：http://emuch.net/bbs/viewthread.php?tid=1996653 2013-10-28 14:46:00
nono2009: 金币-80, 扣除转移金币 2013-10-28 14:46:14
| 好好学学啊！ |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 3chaos.m
2013-10-23 18:36:17, 1.93 K
2楼2013-10-23 18:35:46
【答案】应助回帖
|
function [shang]=jss(xdate) m=2; n=length(xdate); r=0.2*std(xdate); cr=[]; gn=1; gnmax=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 for j=1:m x2m(i,j)=xdate(i+j-1); end end x2m; for i=1:n-m+1 for j=1:n-m+1 if i~=j d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离 k=k+1; end end k=1; end d; for i=1:n-m+1 [k,l]=size(find(d(i,<r));%将比R小的个数传送给L cr1(1,i)=l; end cr1; cr1=(1/(n-m))*cr1; sum1=0; for i=1:n-m+1 if cr1(i)~=0 sum1=sum1+log(cr1(i)); end end cr1=1/(n-m+1)*sum1; cr(1,gn)=cr1; gn=gn+1; m=m+1; end cr; shang=cr(1,1)-cr(1,2); function [shang]=ybs(xdate) m=2; n=length(xdate); r=0.2*std(xdate); cr=[]; gn=1; gnmax=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 for j=1:m x2m(i,j)=xdate(i+j-1); end end x2m; for i=1:n-m+1 for j=1:n-m+1 |
3楼2013-10-23 18:40:08
