ccd2.m

来自「用matlab实现利用统计混沌方法解决非线性系统时间序列预测的问题」· M 代码 · 共 28 行

function dat = ccd2(nbr,len,wid,stm,std);

if nargin<5,
        std = 0.25;
end;
if nargin<4,
        stm = 1;
end;
if nargin<3,
        wid = 1;
end;
if nargin<2,
        len = 3;
end;
if nargin<1,
        nbr = 10000;
end;

pm1 = rand(nbr,stm);
pm2 = rand(nbr,stm);

dat = zeros(nbr,len*wid);

for l = 1:len,
        for w = 1:wid,
                dat(:,(l-1)*wid+(w-1)+1) = sum(exp(-0.5/std.^2*( ((l-0.5)/len-pm1).^2 + ((w-0.5)/wid-pm2).^2 )),2);
        end;
end;