原文链接:tecdat.cn/?p=12307
原文出处:拓端数据部落公众号
我使用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)
最受欢迎的见解
3.Python数据可视化-seaborn Iris鸢尾花数据
7.R语言动态可视化:制作历史全球平均温度的累积动态折线图动画gif视频图