【数值计算】基于四阶Runge-Kutta方法实现Mackey_Glass混沌时间序列生成以及相位时差的绘制附matlab代码

1 内容介绍

基于四阶Runge-Kutta方法实现Mackey_Glass混沌时间序列生成以及相位时差的绘制.

2 部分代码

function [x,t]=Mackey_Glass(N,tau)

% 麦克-格拉斯(Mackey-Glass)混沌延迟微分方程 

%N=10000; 

t=zeros(N,1);

x=zeros(N,1);  

x(1)=1.2; t(1)=0; 

a=0.2;b=0.1;h=0.1;

%tau=30; 

for k=1:N-1

   t(k+1)=t(k)+h; 

   if t(k)

       k1=-b*x(k); 

       k2=-b*(x(k)+h*k1/2); 

       k3=-b*(x(k)+k2*h/2); 

       k4=-b*(x(k)+k3*h);

       x(k+1)=x(k)+(k1+2*k2+2*k3+k4)*h/6; 

   else 

       n=floor((t(k)-tau-t(1))/h+1); 

       k1=Df(x(n))-b*x(k); 

       k2=Df(x(n))-b*(x(k)+h*k1/2); 

       k3=Df(x(n))-b*(x(k)+2*k2*h/2); 

       k4=Df(x(n))-b*(x(k)+k3*h); 

       x(k+1)=x(k)+(k1+2*k2+2*k3+k4)*h/6; 

   end 

end

%

%subplot(121);plot(t,x);

%Dt=tau;  % 相图时间差 

%subplot(122);

%plot(x((Dt*10+1):end),x(1:end-10*Dt));

%dlmwrite('x.txt',x,'-append','delimiter', '\t');

3 运行结果

4 参考文献

[1]吴世宝, 要淞洋, 马佳,等. 基于四阶Runge-Kutta方法的柴油密度数值模拟[J]. 现代商贸工业, 2020, 41(6):2.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机、雷达通信、无线传感器等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。