sigrand2.m

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

function psmp = sigrand2(nsig,nsmp);

psmp = zeros(nsig,nsmp);
accu = 5 * (0.2 + randn(nsig,nsmp+1)) / nsmp;
freq = 5 * (2 + rand(nsig,1));

for smp = 1:nsmp,
        accu(:,smp+1) = accu(:,smp+1) + accu(:,smp);
        psmp(:,smp) = sin(2*pi*accu(:,smp+1).*freq).^2;
end;