Lorenz系统MATLAB仿真与可视化分析

我使用MATLAB解决以下Lorenz初始值问题：

我编写了一个函数，该函数将三个微分方程组作为输入，并使用 带有步长 的Runge-Kutta方法求解该系统。我使用MATLAB生成了解决方法的GIF。

%   输入：
%      f1,f2,f3 = y'(t,y) as a string
%            y0 = initial condition
%         inter = interval
%             h = step size

 
% 调用函数求解Lorenz方程
L=LorenzRK('-10*y1+10*y2','-y1*y3+28*y1-y2','y1*y2-(8/3)*y3',[0,50],[5,5,5],.01)
 
% 绘制Lorenz解
plot3(L(:,2),L(:,3),L(:,4))
     
% 创建GIF：

 
L=LorenzRK('-10*y1+10*y2','-y1*y3+28*y1-y2','y1*y2-(8/3)*y3',[0,50],[5,5,5],.01)

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删