二阶常微分返慕程的数值求解.pptVIP

二阶常微分方程的数值求解; 考虑如下的二阶微分方程初值问题;利用Euler方法求解上述方程组可得如下数值格式;利用四阶R-K方法求解上述方程组可得如下数值格式;例1：用 Euler 法求解如下初值问题;clc;clear; h=0.1; a=0;b=2; x=a:h:b; y(1)=1; z(1)=-1; for i=1:length(x)-1 y(i+1)=y(i)+h*z(i); z(i+1)=z(i)+h*y(i); end plot(x,y,r+,x,exp(-x),k-); xlabel(Variable x); ylabel(Variable y); ;数值解与真解如下图;例2：利用4阶R-K方法求解例1，并与Euler方法进行比较。;function w=rightf_sys1(x,y,z) w=y;;L2=rightf_sys1(x(i)+0.5*h,RK_y(i)+0.5*h*K1,RK_z(i)+0.5*h*L1); % K2 and L2 K3=RK_z(i)+0.5*h*L2; L3=rightf_sys1(x(i)+0.5*h,RK_y(i)+0.5*h*K2,RK_z(i)+0.5*h*L2); % K3 and L3 K4=RK_z(i)+h*L3; L4=rightf_sys1(x(i)+h,RK_y(i)+h*K3,RK_z(i)+h*L3); % K4 and L4 RK_y(i+1)=RK_y(i)+1/6*h*(K1+2*K2+2*K3+K4); RK_z(i+1)=RK_z(i)+1/6*h*(L1+2*L2+2*L3+L4); end plot(x,Euler_y,r+,x,exp(-x),k-,x,RK_y,b*); xlabel(Variable x); ylabel(Variable y); ;例3：分别用 Euler 法和R-K4求解如下初值问题;当 h=0.1，即 n=20 时，Matlab 源程序见 RK_sys2.m, 数值结果如下图;function w=rightf_sys2(x,y,z) w=-y+2*exp(-x)*(x-1);;K2=RK_z(i)+0.5*h*L1; L2=rightf_sys2(x(i)+0.5*h,RK_y(i)+0.5*h*K1,RK_z(i)+0.5*h*L1); % K2 and L2 K3=RK_z(i)+0.5*h*L2; L3=rightf_sys2(x(i)+0.5*h,RK_y(i)+0.5*h*K2,RK_z(i)+0.5*h*L2); % K3 and L3 K4=RK_z(i)+h*L3; L4=rightf_sys2(x(i)+h,RK_y(i)+h*K3,RK_z(i)+h*L3); % K4 and L4 RK_y(i+1)=RK_y(i)+1/6*h*(K1+2*K2+2*K3+K4); RK_z(i+1)=RK_z(i)+1/6*h*(L1+2*L2+2*L3+L4); end plot(x,Euler_y,r+,x,cos(x)+x.*exp(-x),k-,x,RK_y,b*); xlabel(Variable x); ylabel(Variable y); ; dsolve 的调用格式; 几点说明;例 5：求微分方程 在初值条件 下的特解，并画出解函数的图形。;例 6;作业