1、首先我们要把微分方程组用MATLAB表示出来function dydt = lorenz_diff(t,y)dydt = [ -10*(y(1像粜杵泳)-y(2)); 30*y(1)-y(1)*y(3)-y(2) -8/3*y(3)+y(1)*y(2)]; 这就是lorzen系统的微分方程组。
2、用龙格库塔法解上述微分方程组，画出相图clearclcy0 = [0,2,9]; [t,y] = ode45('lorenz_diff',缪梨痤刻[0,200],y0);plot(y(:,1),y(:,3),'.');grid on;
3、结果展示（一）xz相图
4、结果展示（一）xy相图
5、结果展示（一）yz相图
