资源描述：
《实验资料报告材料七 常微分方程初值问题地数值解法.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、大学城市学院实验报告课程名称数值计算方法实验项目名称常微分方程初值问题的数值解法实验成绩指导老师（签名）日期2015/12/16一.实验目的和要求1．用Matlab软件掌握求微分方程数值解的欧拉方法和龙格－库塔方法；2．通过实例学习用微分方程模型解决简化的实际问题。二.实验容和原理编程题2-1要求写出Matlab源程序(m文件)，并有适当的注释语句；分析应用题2-2，2-3，2-4，2-5要求将问题的分析过程、Matlab源程序和运行结果和结果的解释、算法的分析写在实验报告上。2-1编程编写用向前欧拉公式和改进欧拉公式求
2、微分方程数值解的Matlab程序，问题如下：在区间个等距点处，逼近下列初值问题的解，并对程序的每一句添上注释语句。Euler法y=euler(a,b,n,y0,f,f1,b1)改进Euler法y=eulerpro(a,b,n,y0,f,f1,b1)2-2分析应用题假设等分区间数，用欧拉法和改进欧拉法在区间求解初值问题并作出解的曲线图形，同时将方程的解析解也画在同一图上，并作比较，分析这两种方法的精度。2-3分析应用题用以下三种不同的方法求下述微分方程的数值解，取画出解的图形，与精确值比较并进行分析。1）欧拉法；2）改进欧
3、拉法；3）龙格－库塔方法；2-1分析应用题考虑一个涉及到社会上与众不同的人的繁衍问题模型。假设在时刻(单位为年)，社会上有人口人，又假设所有与众不同的人与别的与众不同的人结婚后所生后代也是与众不同的人。而固定比例为的所有其他的后代也是与众不同的人。如果对所有人来说出生率假定为常数，又如果普通的人和与众不同的人的婚配是任意的，则此问题可以用微分方程表示为：其中变量表示在时刻社会上与众不同的人的比例，表示在时刻人口中与众不同的人的数量。1）假定和，当步长为年时，求从到解的近似值，并作出近似解的曲线图形。2）精确求出微分方程的
4、解，并将你当时在分题(b)中得到的结果与此时的精确值进行比较。【MATLAB相关函数】n求微分方程的解析解及其数值的代入dsolve(‘egn1’，‘egn2’，‘’)subs(expr,{x,y,…},{x1,y1,…})其中‘egn’表示第个方程，‘’表示微分方程中的自变量，默认时自变量为。subs命令中的expr、x、y为符合型表达式，x、y分别用数值x1、x2代入。>>symsxyz>>subs('x+y+z',{x,y,z},{1,2,3})ans=6>>symsx>>subs('x^2',x,2)ans=4s
5、=dsolve(‘’,‘’,‘’)ans=>>symsx>>subs(s,x,2)ans=-0.3721n右端函数的自动生成f=inline(‘expr’,’var1’,‘var2’,……)其中’expr’表示函数的表达式，’var1’,‘var2’表示函数表达式中的变量，运行该函数，生成一个新的函数表达式为f(var1,var2,……)。>>f=inline('x+3*y','x','y')f=Inlinefunction:f(x,y)=x+3*y>>f(2,3)ans=11n4，5阶龙格－库塔方法求解微分方程数值解[
6、t,x]=ode45(f,ts,x0,options)其中f是由待解方程写成的m文件名；x0为函数的初值；t,x分别为输出的自变量和函数值(列向量)，t的步长是程序根据误差限自动选定的。若ts=[t0,t1,t2,…,tf]，则输出在自变量指定值，等步长时用ts=t0:k:tf，输出在等分点；options用于设定误差限(可以缺省，缺省时设定为相对误差，绝对误差)，程序为：options=odeset(‘reltol’,rt,’abstol’,at)，这里rt,at分别为设定的相对误差和绝对误差。常用选项见下表。选项名功
7、能可选值省缺值AbsTol设定绝对误差正数RelTol设定相对误差正数InitialStep设定初始步长正数自动MaxStep设定步长上界正数MaxOrder设定ode15s的最高阶数1,2,3,4,55Stats显示计算成本统计on,offoffBDF设定ode15s是否用反向差分on,offoff例：解微分方程在命令窗口执行=(‘’,‘’,‘’)；；ans=01.00000.0502 1.04900.10051.09590.15071.14083.85072.95033.90052.96723.95022.98394
8、.00003.0006plot(,,‘o-’,)%解函数图形表示%不用输出变量，则直接输出图形；ans=01.00001.0000 1.73212.00002.23613.00002.64584.00003.0006三.操作方法与实验步骤（包括实验数据记录和处理）2-1编程编写用向前欧拉公式和改进欧拉公式求微分方程数
