[工学]7常微分方程数值解法.pptVIP

第七章 常微分方程数值解法 第一节 问题的提出 含有一个或多个导数及其函数的方程式称为微分方程，在工程中常遇到求解微分方程的问题。 很多微分方程的解不能用初等函数来表示，有时即使能够用解析式表示其解，但计算量太大而不实用(表达式过于复杂)。 需要用数值方法来求解，一般只要求得到若干个点上的近似值或者解的简单的近似表达式(精度要求满足即可)。 ?常微分方程（一个自变量） 微分方程? ?偏微分方程（一个以上自变量） 一. 定解问题 实际中求解常微分方程的所谓定解问题有两类： 初值问题和边值问题 定解指已知因变量和/或其导数在某些点上是已知的(约束条件)。 1. 边值问题 约束条件为已知，在自变量的任一非初值上，已知函数值和/或其导数值，如 ?y ” = f (x,y,y’) ? 求解y ?y(a)=?，y(b)=? 常常可以将边值问题转化为初值问题求解。 2. 初值问题 约束条件为在自变量的初值上已知函数值，如： ?y’= f (x,y) ?dy/dx = f (x,y) x?x0 ? ?? ?y(x0) = y0 ? y(x0) = y0 求解y(x)，以满足上述两式，即在 a?x0?x1?……? xn?b上的y(xi)的 近似值 yi (i=0,1,2,…,n)。 通常取等距节点，即h= xi+1-xi ，有 xi = x0+ih (i=0,1,2,…,n) 初值问题的数值解法特点：按节点顺序依次推进，由已知的y0 , y1 , …, yi ，求出yi+1，这可以通过递推公式得到。 ?单步法：利用前一个单步的信息 ? (一个点)，在y=f(x) 上找下 ? 一点yi， 有欧拉法， ? 龙格－库格法。 ? 初值问题的常见解法? ? ?预测校正法：多步法，利用一个以 ? 上的前点信息求f(x)的 ? 下一个yi，常用迭代法 ? 如改进欧拉法、 ? 阿当姆斯法。 第二节 欧拉法 ?y’= f (x,y) ? 求解计算y(xi) (i=1,2,…,n)的近似值yi ?y(x0) = y0 一、欧拉折线法 1. 解一阶方程初值问题的几何意义 y’= f (x,y) (x,y)?R，R={a?x?b，-??y?+?} 也即解y = y(x)所表示的积分曲线 y =∫ f(x,y) dx+c 上，每一点(x，y)的切线斜率等于函数f(x,y)在该点的值，即： y ’( xk) = f (xk，yk) y(x0)=y0，表示积分曲线从P0（x0，y0）点出发且在 P0（x0，y0）的切线斜率为f（x0，y0），故设想积分曲线y = y (x)在x = x0 附近可用切线近似代替。 2. 欧拉折线法 在 （ x0，y0 ）点上的切线方程为 y=y0+ f (x0 , y0)(x-x0) 设方程与x = x1交点P1（x1，y1）纵坐标y1取为y(x1)的近似值，则有 y(x1) ? y1 = y0+ f (x0,y0)(x1-x0) = y0+ h f (x0,y0) 同理 ：在(x1，y1 )上的切线方程 y = y1+ f (x1,y1)(x-x1)与x=x2交点P2(x2,y2)的纵坐标y2取为y(x2)的近似值有 y(x2) ? y2 = y1+ f (x1,y1)(x2-x1) = y1+ h f (x1,y1) 上述的h = x2－x1= x1－x0 过Pi(xi,yi)作斜率为f(xi,yi)的切线方程与x = xi+1的交点， 有欧拉公式 ? y(xi+1) ? yi+1 = yi+ h f (xi，yi) (i=0,1,…,n-1) ? ?xi = x0 +i