China-bub.com 第11章积分和微分方程组 157 下载 ■例11.2 (a)求解下面的ODE: x/=-x2 x(0)=1 创建函数xprim1,将此函数保存在M文件xpriml.m中： function xprim xprimi(t,x) xprim -x.2; 然后，调用MATLAB的ODE算法求解方程，最后画出解的图形： [t,x]=ode45('xprim1',[01],1); p1ot(t,x,’-’,t,x,’0); xlabel('time to =0,tt =1'); ylabel('x values x(0)=1'); 得到图11-2。MATLAB计算出的解用圆圈标记。在13.1节中介绍绘图命令p1ot。 05 05 0.7 a85 06 Q10.2030.4·050.6070.809 m80=0量=1 图1l-2由函数xprim1定义的ODE解的图形 (b)解下面的ODE过程是等价的： [x'=x2 x(0)=1 首先创建函数xprim.2,并将其保存在M文件xprim2.m中： function xprim xprim2(t,x) xprim x.2; 然后调用ODE的求解方程并画出其解的图形：■ 例11 . 2 (a) 求解下面的O D E： 创建函数x p r i m1，将此函数保存在M文件x p r i m 1 . m中： 然后，调用M AT L A B的O D E算法求解方程，最后画出解的图形： 得到图11 - 2。M AT L A B计算出的解用圆圈标记。在 1 3 . 1节中介绍绘图命令p l o t。 图11-2 由函数x p r i m 1定义的O D E解的图形 (b) 解下面的O D E过程是等价的： 首先创建函数x p r i m2，并将其保存在M文件x p r i m 2 . m中： 然后调用O D E的求解方程并画出其解的图形： 第11章 积分和微分方程组 1 5 7 下载