《常微分方程及其数值解法》实验指导书.docVIP

《常微分方程及其数值解法》 实验指导书 姓名 ________________________ 班级 ________________________ 学号 _______________________ 东北电力大学 理学院 2010年 《常微分方程及其数值解法》课程实验教学大纲 一、课程类别及学时 专业必修课，总学时6学时 二、教学对象 信息与计算科学、数学与应用数学专业二年级学生 三、教学目的 《常微分方程及其数值解法》实验是通过，使学生能够，培养实际动手能力和提高学生分析问题、解决问题的能力为以后从事现代数学科研工作以及工程和科学计算实践打下良好的基础。（1）； Euler方法及其改进方法.................................................4 实验二 Runge-Kutta方法..............................................................7 实验三 数值计算方法稳定域、收敛域检验………………….....9 实验一 Euler方法及其改进方法 一、实验目的与要求 1．通过用Matlab编程运用Euler方法及其改进方法求解常微分方程初值问题，更进一步掌握常微分方程及其数值解法课程的理论内容，加深对数值解法的理解。 2．熟悉Matlab编程环境。 二、实验内容和指导（步骤） 【实验题目】运用Euler方法及其改进方法求解常微分方程初值问题，具体问题在教材课后习题中学生任意选取。 【实验原理】 Euler方法： Euler改进方法： 【设计思想】 选择教材中能够获得精确解的题目，分别应用Euler方法、Euler改进方法求解数值结果，将所得结果列表或者画图，参照精确结果，对Euler方法、Euler改进方法的计算误差进行分析。 源代码： %fun为目标函数字符串 %x0为自变量初始值。 %y0为fun(x0); %bou=[a,b]自变量区间 %h为步长 fun=1/x^2-y^2; bou=[1,6]; a=bou(1); b=bou(2); x0=1; y0=0; h=0.1; n=ceil((b-a)/h); xx=linspace(a,b,n+1); yy=zeros(1,n+1); lengthx=length(xx); xx(1)=x0;yy(1)=y0; for i=2:n+1 x=xx(i-1);y=yy(i-1); k=eval(fun); yy(i)=yy(i-1)+h*k; end Ys=dsolve(Dy=1/x^2-y^2,y(1)=0,x); for i=1:lengthx x=xx(i); exacty(i)=eval(Ys); end YY=exacty; yend=[xx,yy,YY] p=plot(xx,yend(:,2),k-o,LineWidth,1,... MarkerEdgeColor,k,... MarkerFaceColor,g,... MarkerSize,4); hold on; pp=plot(xx,yend(:,3),r-.+,LineWidth,0.8,... MarkerEdgeColor,r,... MarkerFaceColor,m,... MarkerSize,6); legend([p,pp],Eula,jiequejie); %fun为目标函数字符串 %x0为自变量初始值。 %y0为fun(x0); %bou=[a,b]自变量区间 %h为步长 fun=1/x^2-y^2; bou=[1,6]; a=bou(1); b=bou(2); x0=1; y0=0; h=0.1; n=ceil((b-a)/h); xx=linspace(a,b,n+1); yy=zeros(1,n+1); lengthx=length(xx); xx(1)=x0;yy(1)=y0; for i =2:n+1 x=xx(i-1);y=yy(i-1); k1=eval(fun); x=xx(i);y=yy(i-1)+h*k1; k2=eval(fun);