examp1_4.m

matlab的大量实例代码

function examp1_4
[t,y]=ode45(@lorenzeq,[0,100],[0;0;1e-10]); % 求解方程
plot3(y(:,1), y(:,2), y(:,3))  % 绘制三维相轨迹
axis([10 42 -20 20 -20 28]);   % 手动设置坐标系位置
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)];