数学实验微分方程组的解.doc

新乡学院数学与信息科学系 实验报告 实验项目名称 微分方程（组）的解 所属课程名称 数学实验 实 验 类 型 微分方程实验 实 验 日 期 2013-5-10 班 级 学 号 姓 名 成 绩 一、实验概述： 【实验目的】 熟悉matlab函数文件编写 灵活运用循环语编写程序 用作图来观察ode45、ode23、的差别 掌握matlab文件操作和画图操作． 【实验原理】 1. 求解析解：y=dsolve(eq1,eq2, ... ,cond1,cond2, ... ,v) ；其中 y 为输出， eq1、eq2、...为微分方程，cond1、cond2、...为初值条件，v 为自变量。 2．求数值解： ode113、ode23t、ode15s、 ode23s、ode23tb 3. 初值问题的Euler折线法 具体步骤：分割求解区间，差商代替微商，解代数方程 【实验环境】 Matlab 7.0.1 Microsoft Windows 7 二、实验内容： 【实验方案】 1．利用matlab求微分方程数值解。并通过ode23、ode45这两种不同的求解器来画出图，通过图来形象的观察来观察出两种求解器的不同。 2．运用Euler折线法求解微分方程初值。 【实验过程】（实验步骤、记录、数据、分析） 1 求微分方程组 其中，分别用 ode23、ode45 求题中的微分方程初值问题的数值解(近似解)，求解区间为t∈[0,0.5]．利用画图来比较两种求解器之间的差异． function dy=myfun1(t,y); dy=[-y(1)-y(2);y(2)-y(1)]; [t,y]=ode23(myfun1,[0,0.5],[1,0]) y1=y(:,1); plot(t,y1,r*) hold on [t,x]=ode45(myfun1,[0,0.5],[1,0]) y1=y(:,1); plot(t,y1,b*); 2 用 Euler 折线法求解微分方程初值问题 的数值解(步长h取0.001)，求解范围为区间[0,2]． f=sym(y-12*x^2/y^3); a=0; b=2; h=0.001; n=(b-a)/h+1; x=0; y=1; szj=[x,y]; for i=1:n-1 y=y+h*subs(f,{x,y},{x,y}); x=x+h; szj=[szj;x,y]; end szj plot(szj(:,1),szj(:,2)) 【实验结论】（结果） 1 t = 0 0.0001 0.0005 0.0025 0.0125 0.0625 0.1125 0.1625 0.2125 0.2625 0.3125 0.3625 0.4125 0.4625 0.5000 x = 1.0000 0 0.9999 -0.0001 0.9995 -0.0005 0.9975 -0.0025 0.9877 -0.0125 0.9413 -0.0626 0.8997 -0.1130 0.8626 -0.1639 0.8298 -0.2157 0.8011 -0.2685 0.7765 -0.3227 0.7557 -0.3786 0.7387 -0.4363 0.7254 -0.4962 0.7179 -0.5427 通过作图可以明显的看出ode23、ode45这两种求解器之间的差别，明显的ode23比较适用于精确度比较低的数值求解，相对误差比较大，而de45则适用于求精确度比较高的数值求解，相对于ode23来说，误差比较小。 2 【实验小结】（收获体会） 很好的掌握用Matlab编写M文件的方法，能够熟练的掌握用Matlab来求出微分方程的数值解以及通