我使用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)