资源简介

用Runge-Kutta方法求解Mackey-Glass混沌时间序列的matlab源程序

资源截图

代码片段和文件信息

N=10000;
t=zeros(N1);
x=zeros(N1);
x(1)=1.2;
t(1)=0;
h=0.1;
ta=17;
for k=1:N-1
    t(k+1)=t(k)+h;
    if t(k)        k1=-0.1*x(k);
        k2=-0.1*(x(k)+h*k1/3);
        k3=-0.1*(x(k)+(-k1+3*k2)*h/3);
        k4=-0.1*(x(k)+(k1-k2+k3)*h);
        k5=-0.1*(x(k)+(k1+3*k2+3*k3+k4)*h/8);
        x(k+1)=x(k)+(k1+3*k2+3*k3+k4)*h/8;
    else
        n=floo

评论

共有 条评论