基于matlab的Lorenz系統仿真可視化

原文連結：http://tecdat.cn/?p=12307

我使用MATLAB解決以下Lorenz初始值問題：

我編寫了一個函數LorenzRK4IVP（），該函數将三個微分方程組作為輸入，并使用 帶有步長的Runge-Kutta方法求解該系統。我使用MATLAB生成了解決方案的GIF。

%   Inputs:%      f1,f2,f3 = y'(t,y) as a string%            y0 = initial condition%         inter = interval%             h = step size% call function to solve Lorenz equationsL=LorenzRK('-10*y1+10*y2','-y1*y3+28*y1-y2','y1*y2-(8/3)*y3',[0,50],[5,5,5],.01)% plot Lorenz solutionsplot3(L(:,2),L(:,3),L(:,4))% create GIF source: L=LorenzRK('-10*y1+10*y2','-y1*y3+28*y1-y2','y1*y2-(8/3)*y3',[0,50],[5,5,5],.01)           

複制