计算方法第六章精要.pptVIP

* } 第6章 常微分方程数值解法 6.1 欧拉法和改进的欧拉法 6.2 龙格-库塔法 6.3 线性多步法 引言 } 一阶常微分方程的初值问题 当f(x,y)连续,且满足Lipschitz(李普希兹)条件 ? |f(x,y1)–f(x,y2)|≤L|y1–y2|, x∈[a,b], y1,y2∈R ?则上述初值问题存在连续可微的解函数 y = y(x). 例如初值问题 可求出方程y′=1+ex的通解为 y=x+ex+c, 将初值条件 x=0, y=2 代入得 2=1+c, 故 c=1, 所以初值问题的解为 y=x+ex+1 引言 } 又如初值问题 可求出它的解为 但要进一步计算指定点的函数值，还需要用数值积分方法。有些微分方程的解是隐函数，例如 要求函数值还需要解超越方程。应用中所处理的微分方程往往很复杂且大都得不出一般解，所以一般用数值解法。 数值解法： 给定节点a=x0x1…xN=b, 将初值问题离散化为差分方程,求出解函数 y(x) 在这些点的近似值y1 ,y2 ,…,yn 。 所求得的近似值称为数值解。 } 6.1 欧拉法和改进的欧拉法 6.1.1 欧拉法及其截断误差 6.1.2 改进的欧拉法及预测-校正公式 } 6.1.1 欧拉法及其截断误差 1. 公式的构造思想：用差商代替导数 设 等距，步长为 令x=xi , x+h=xi+1 , y(xi )≈yi ，y(xi+1 ) ≈yi+1 ,初值问题离散化为 初值问题 (欧拉公式) } 例 6.1 取步长 h=0.1，用欧拉法求解初值问题 (欧拉公式) 解 y1=y0+h f(x0 ,y0 )=1+0.1(0 + 1 )=1.1 y2=y1+h f(x1 ,y1 )=1.1 + 0.1(0.1 + 1.1 )=1.22 y3=y2+h f(x2 ,y2 )=1.22+0.1(0.2+ 1.22)=1.22+0.142=1.362 ……… y10=y9+h f(x9 ,y9 )=y9+0.1(x9 + y9 )=3.18748 } 2. 几何意义——用折线代替曲线计算解函数的近似值。 准确解为 } 3. 数值公式的误差来源。 (1) 局部截断误差（简称截断误差）：假设 yi=y(xi )是准确的 ， 计算yi+1所产生的误差 y(xi+1 ) - yi+1 若局部截断误差可以表示为O(hk+1), k为正整数， 则称公式是k阶公式。 （2）由于实际上yi不是准确值，因此它的误差会传播下去。 （3）实际计算时，每一步都可能产生舍入误差。 } 4. 欧拉公式的截断误差是O(h2)，公式是1 阶的。 因为 （泰勒公式） 两式相减，由设 yi=y(xi ) ，有 } 6.1.2 改进的欧拉法及预测-校正公式 对微分方程y′=f(x,y) 两边求xi 到xi+1 的定积分，有 利用梯形公式计算积分，有 将y(xi ) 、y(xi+1 )分别用yi 、yi+1 代替，构造相应的数值公式： （改进的欧拉公式） 1. 公式的构造 } 2. 截断误差 假设 yi=y(xi), 解函数在x=xi处的泰勒公式为 在改进的欧拉公式中， 设 则有 求出在h=0处的泰勒公式，整理后得 上式h 和h2 项的乘数应为零，于是 } 因而改进的欧拉法是二阶的。 } 3 改进的欧拉法的具体使用格式。 改进的欧拉法是隐式公式 ,计算时常用迭代法。一般每一步先由欧拉公式计算出yi+1 的初始值yi+1(0),再迭代计算yi+1。 当满足 时，取 可证明当f(x,y)满足一定条件时，迭代是收敛的。 } 改进的欧拉法的预测－校正公式 可证明预测－校正公式的截断误差也为 O(h3)。 } 例 取步长h=0.2,用改进的欧拉法的预测－校正公式 求解初值问题 的数值解y1 , y2 . 解 预测-校正公式具体是 } } 6.2 龙格-库塔法 6.2.1 二阶龙格-库塔（Runge-Kutta)公式 6.2.2 四阶龙格-库塔公式 引言 与其他算法的实例比较 } 思想：从泰勒公式出发，寻找更高阶的数值公式。 例如，泰勒公式计算到二阶可得 令 则 ，略去余项，得出一个二阶的 理论上按此方式可以得到更高阶的公式。但需要计算复合函数的高阶导数，使算法复杂而不实用。 因 数值公式为 } 龙格和库塔的思想（间接地运用泰勒公式）： 利用y(x)在若干个点上的函数值和导数值，作出一个适当的线性组合，使这个线性