| 查看: 847 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
[求助]
求混沌或者噪声熵的程序
|
|||
大家好,谁有求混沌或者噪声熵的程序,c或者matlab都可以。跪求 |
» 猜你喜欢
厉娜
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有12人回复
售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急
已经有15人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急
已经有16人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有15人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
【答案】应助回帖
|
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










回复此楼