-
大小: 746B文件类型: .m金币: 1下载: 0 次发布日期: 2021-06-26
- 语言: Matlab
- 标签: Mackey-Glass Runge-Kutta
资源简介
用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
评论
共有 条评论