ccd1.m

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

function dat = ccd1(nbr,len,stm,std);

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

prm = rand(nbr,stm);

dat = zeros(nbr,len);

for l = 1:len,
        dat(:,l) = prod(exp(-0.5/std.^2*((l-0.5)/len-prm).^2),2);
end;