examp1.m

来自「基于MATLAB_Simulink的系统仿真技术与应用大家可以学习下」· M 代码 · 共 8 行

function examp1
[t,y]=ode45(@lorenzeq,[0,100],[0;0;1e-10]);
plot3(y(:,1), y(:,2), y(:,3))  % 绘制三维相轨迹
function xdot = lorenzeq(t,x) % 下面的函数直接对应于给出的微分方程
xdot=[-8/3*x(1)+x(2)*x(3);
      -10*x(2)+10*x(3);
      -x(1)*x(2)+28*x(2)-x(3)];