%% ode45解决微分方程 %无附加参数 %Lorenz 模型状态方程 F=@(t,x)[-8/3*x(1) x(2)*x(3); %dx1 -10*x(2) 10*x(3); -x(1)*x(2) 28*x(2)-x(3) tn=100；%终止时间 x0=[0;0;1e-10]; -10=10^(-10); [t,x]=ode45(F,[0,tn],x0); figure() plot(t,x,'Linewidth',1.3) figure() plot3(x(:,1),x(:,2),x(三)Linewidth',1.3) figure() comet3(x(:,1),x(:,2),x(3)%三维彗星图