数学实验8.docVIP

数学实验(八) 实验内容：微分方程的求解方法及其应用 实验目的：掌握微分方程的求解方法及其应用 实验重点：微分方程的求解方法及其应用 实验难点：微分方程的求解应用 实验关键：微分方程的求解的适应条件 实验过程： 一、微分方程基本求解格式 1、Dsolve[equ,y,x],求微分方程的解析解（精确解）DSolve[方程，y[x],x]给出独立变量为x的微分方程的一般解y[x].DSolve[方程，y,x]给出微分方程的一般解y，这里把解表示为列表内的纯粹函（就是不引用特定变量进行定义函数，这时定义的变量依次用#1，#2等均表示。 函数的格式为： 2、NDSolve[equ,y,{x,,}],给出微分方程的数值解（近似解）求解常微分方程数值解的函数NDSolve的一般形式是： NDSolve[eqns,y,{x,xmin,xmax}]对常微方程或方程组eqns，求函数y关于x在[xmin,xmax]范围内的数值解NDSolve[{eqn1,eqn2,……},{y1,y2,……},{x,xmin,xmax}] 对常微方程组，求函数y1,y2,……[xmin,xmax]范围内的数值解 NDSolve[eqns,y,{x,xmin,xmax}，{t,tmin,tmax}] 二、微分方程的求解 例1:求解微分方程+2xy=x并加以验证 说明： ①方程中的函数应写成完整形式y[x]，以表明y是x的函数； ②方程必须写成…==…的形式，如果是多个方程，还必须写成表的形式，这些方程是表中的元素； ③第一行：是经常要做的事，以清除变量的以前定义，减少疑义，保证程序正常运行； ④程序的运行结果使一个表，表中的元素以纯函数的形式给出微分方程的解； 所谓的纯函数相当于数学中的算子，仅定义了函数名本身，而没有具体给出谁使自变量。 ⑤第三行：一个代入规则，同样没有具体给出谁是自变量，结果为 ⑥第四行：具体给出x是自变量，结果为 ⑦第五行：使用所得的解，这里将解代入原微分方程，结果应该为0，但这里给出： ⑧第六行：用Simplify[……]函数化简，结果为0，表明y=y[x]的确是微分方程的解。 例2:求微分方程x+y-=0在初始条件y(1)=2e下的特解，并画出解函数的图形。 微分方程的解为： 其特解为：y= 例4:求微分方程组在初始条件x,y=0下的特解，并画出解函数的图形。 解： 求解本问题的Mathematica程序为: 其特解为： 解函数的图形如下： 例5：解初值问题 （0）并计算y(0)和y(0.7)的值。 三、Euler 折线法的数值求解 例1：用Euler 折线法求解微分方程初值问题 的数值解（步长h取0.1），求解范围为区间[0,1] 解 本问题的差分方程为 h=0.1 k=0,1,2……，n-1 (其中：f[x,y]=y+) 相应的Mathematica程序为： 结果表示为矩阵： 5