实验4 常微分方程的数值解几道解方程的题——重要.PDF

实验 4   常微分方程的数值解  几道解方程的题——重要  2 2 2 1.   x y’’+xy’+(x ‐n )y=0    2 2     y ( ) =2,y’ ( ) =‐      (贝赛尔方程，令 n=0.5)，精确解y=sin x   2 2  x 解：首先将高阶方程装降阶化简为一阶常微分方程组  令 y1=y’ ，令y2=y1’=y’’  将 n=0.5 代入，则原方程转化为：  y ’=y   1 2 y 2 1 (1 )y y 2    2 1  x 4x   2 y1 ( ) =2, y2 ( ) =‐   2 2  function dy=fangcheng(x,y) dy=[y(2);-y(2)/x-(1-1/(4*x^2))*y(1)];   xs=pi/2:0.1:100; y0=[2,-2/pi]; [x,y]=ode45(@fangcheng,xs,y0); [x,y] plot(x,y(:,1)),grid   2001    x  2． 设y (x) y (x) sin x ye 0, y (0) 1, y (0) 0,   用数值解法算出    y(1)=               （精确到4 位小数）, 你用的方法是                    ，调用的  Matlab 命令是                          ，算法精度为          。  function dy=fangcheng2(x,y) dy=[y(2);y(2)*sin(x)-y(1)*exp(x)]; xs=0:0.01:1; y0=[1,0]; [x,y]=ode45(@fangcheng2,xs,y0); [x,y]   2000    1.  设y (x ) y (x ) sin x 0, y (0) 1, y (0) 0, 用数值解法算出 y(1)=              ,你 用的方法是                    ，调用的  Matlab 命令是                          ，算 法精度为          。  function dy=fangcheng1(x,y) dy=[y(2);y(1)*sin(x)]; xs=[0,1]; y0=[1,0]; [x,y]=ode45(@fangcheng1,xs,y0); [x,y]   题目 1  2 2 2 1.   x y’’+xy’+(x ‐n )y=0    2 2       y ( ) =2,y’ ( ) =‐      (贝赛尔方程，令 n=0.5)，精确解y=sin x 2 2  x 解：首先将高阶方程装降阶化简为一阶常微分方程组  令 y1=y’ ，令y2=y1’=y