3-常微分方程及其求解实例.pptVIP

第4章 常微分方程及其求解 例： 自变函数 t_final=100; x0=[0;0;1e-10]; % t_final为设定的仿真终止时间 [t,x]=ode45(lorenzeq,[0,t_final],x0); plot(t,x), figure; % 打开新图形窗口 plot3(x(:,1),x(:,2),x(:,3)); axis([10 42 -20 20 -20 25]); % 根据实际数值手动设置坐标系 可采用comet3( )函数绘制动画式的轨迹。 comet3(x(:,1),x(:,2),x(:,3)) 描述微分方程是常微分方程初值问题数值求解的关键。 f1=inline([[-8/3*x(1)+x(2)*x(3); -10*x(2)+10*x(3);... -x(1)*x(2)+28*x(2)-x(3)]],t,x); t_final=100; x0=[0;0;1e-10]; [t,x]=ode45(f1,[0,t_final],x0); plot(t,x), figure; plot3(x(:,1),x(:,2),x(:,3)); axis([10 42 -20 20 -20 25]); 得出完全一致的结果。 （3） 延迟微分方程求解 sol:结构体数据，sol.x:时间向量t, sol.y:状态向量。