微分方程作业.docxVIP

P10习题 1.用Euler法和改进的Euler法求u =?5u (OWtWl), u(0)=l的数值解，步长h=0. 1, 0. 05;并 比较两个算法的精度。 解:function du=Euler funl (t, u) du=-5*u;clear; h=0.1; tend=l;N=l/h;t (l)=0;u(l)=l; t=h. *(0:N); for n=l:N u(n+l) =u(n) +h*Eu 1 er_funl(t(n),u(n)); end plot (t, u, *) ;hold on for n=l:N v(l)=u (n) +h*Euler_funl (t (n), u(n)); for k二1:6 v (k+l)二u (n) +h/2* (Euler_funl (t (n), u(n))+Euler_funl (t (n+1), v (k))); end u(n+l)=v (k+1); end plot (t,叮 o); sol=dsolve(，Du二-5*u, u(0)=T ); u_real=eval(sol); plot (t, unreal, r‘)； 将上述h换为o?05得: 由图像知道： 显然改进的Euler法要比Euler法精确度要高； 3.将u二-u(OWtWl),u(O)=O, u (0)=1化为一阶方程组，并用Euler法和改进的的Euler 法求解，步长h=0.1,0. 05；并比较两个算法的精度。 解： function du=fun31 (y) clu=y; function dy=fun32(u) dy二一u; clear; h=0. 1; tend=l ;N=l/h; t (1)=0; u (1)=0; y(l) = ; t=h. *(0:N); for n=l:N u(n+1) =u(n)+h*y (n); y (n+1) =y (n) +h* (-u (n)); end plot (t,叮 *，);hold on for n=l:N v(l)=u(n)+h*fun31(y(n)); w(1)=y(n)+h*fun32(u(n)); for k二1:6 v(k+l)=u (n)+h/2* (fun31 (y (n)) +fun31 (... w (k))); w (k+1) =y (n) +h/2* (fun32 (u (n)) +fun32 (... v (k))); end u(n+1)二V(k+1); y (n+1)二\v (k+1); end plot (t, u/ o); sol=dsolve(，D2u二-u‘, u(0)=0,, Du(0)=T ; u_real=eval(sol); plot (t, unreal, F ); 将上述h换为o?05得: 由图像可以知道： 显然改进的Euler法要比Euler法精确度要高; 实习题（二） 1.取步长h = QA ,分别用Euler法和改进的Euler法求下列初值问题的解，并与真解 相比较. 9 Y u= a- — , 0 x 1, u u(0) = h 真解 u(x) = Jl + 2x ； 解： function du=funl(x, u) du=u~2*x/u; clear; h 二 0. 1 ;xend二 1 ;N二 l/h;x(l)二0;u(l)二 1; x=h.*(0:N); %——Eluer法——% for n=l:N u (n+l)=u (n) +h*fun 1 (x (n), u(n)); end plot (x, u, ; hold on %——改进的Eluer法——% for n=l:N v(l)=u(n) +h*funl (x(n), u(n)); for k二1:6 v(k+l)=u(n) +h/2*(funl (x(n), u(n))+funl (x(n+1), v(k))); end u(n+1)二v(k+1); end plot (x,叮 o); hold on %——真解——% u real=sqrt(l+2*x); plot (x, unreal, F ); 由图像可以知道： 显然改进的Euler法要比Euler法精 确度要高； ,U X 1 “ ⑵八厂卩心2, M(l) = 2, 真解 u(x) = x(8-31nx)3 ； 解： function du=fun2(x, u) du= (u/x)-x. *2/u. 2; clear; h=0. l;N=l/h;x=l:h:2;x(l)=l;u(l)=2; for n=l:N u(n+l)=u (n) +h*fun2 (x (n), u(n)); end plot (x, u, ; hold on for n=l:N