常微分方程初值问题的数值解法.doc

《常微分方程初值问题的数值解法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第九章常微分方程初值问题的数值解法第一部分内容提要一、数值解的一般概念常微分方程初值问题的数值解是指通过一定的近似方法得出准确解在一列离散点上的近似值。数值解的特征是步进式，即在点的近似值是由等若干点处的近似值的信息给出的递推公式。若依赖于前面步的值，则称为步法；称为单步法。b5E2RGbCAP利用在的精确解借助某种算法计算出，则称为该方法的局部截断误差。如果一个算法的局部截断误差是，则称该方法是阶的；而利用数值解得到的与微分方程的精确解之差称为整体截断误差，即是该数值方法的误差。p1EanqFDPw对于固定的，取，用某种算法得到，如有=0，则称该方法是收敛的。注意，因是固

2、定的，随着，数值解的步数。在实际计算时由于舍入误差不可避免，实际得到数值解是，稳定性即研究是否随着计算步骤的增加而增加。通常所提的稳定性是通过模型方程来讨论的。若当某一步有舍入误差时，在以后的计算中误差不会逐步扩大，则称这种稳定性为绝对稳定性。DXDiTa9E3d二、简单单步法及其收敛性、稳定性6/6Euler法的局部截断误差为，整体截断误差为，即一阶收敛。对于模型问题，当时，Euler法是数值稳定的。RTCrpUDGiT隐式Euler法的误差与Euler法相同，但是无条件稳定：即对任意步长，隐式Euler法都是稳定的。梯形法的误差比Euler法高一阶，也是无条件稳定的。改

3、进Euler法是一种预测－校正方法：——Euler法预测——梯形法校正它保持了梯形法的误差阶数，但不是无条件稳定的。一、龙格-库塔方法龙格-库塔类算法采用区间内若干点的斜率的加权平均来近似整个区间的平均斜率，一般形式为如经典的4级4阶<局部截断误差为）Runge-Kutta公式为6/6<1）四、线性多步法线性步法具有一般形式<2）为隐式公式；为显式公式。构造多步法公式有基于数值积分和Taylor展开两种途径。多步法<2）的局部截断误差为5PCzVD7HxA利用原微分方程后，成为因此利用Taylor公式，分别对和作Taylor展开，可确定线性多步<2）中的待定参数和，使它达到

4、最高阶精度<或指定精度）。jLBHrnAILg预测－校正格式：不论单步法还是多步法，隐式公式比显式公式的稳定性好，但隐式公式的计算比较困难。预测－校正格式是用显式公式进行预测，再用隐式公式进行校正。xHAQX74J0X五、高阶方程和一阶常微分方程组所有应用于一阶常微分方程初值问题的数值方法都可以直接推广到一阶常微分方程组，只需把公式中的未知函数改为向量形式。高阶方程总可以降为一阶方程组，进而用数值方法求解。第二部分例题精选一、对初值问题，，用梯形公式求解，求数值解的表达式<写成步长的函数），并证明数值解的收敛性。分析：微分方程用数值方法离散后即变成差分方程6/6，单步法导出

5、的差分方程通常是关于数值解序列的一个递推公式，因此问题变为已知数列首项和递推公式求数列通项公式，收敛性即是数列极限。LDAYtRyKfE解：梯形公式对本问题于是梯形公式的解即由于，易得显然初值问题的准确解为对于固定的点，准确解在该点的值为而数值解为令固定，上式对取极限得：这证明了数值解的收敛性。二、用隐式单步法<该方法也属于隐式Runge-Kutta方法）求解微分方程初值问题，时，试推出其绝对稳定区间。分析：将格式应用于所给方程可导出误差传播方程，从而求出绝对稳定区间解：记舍入误差为。该隐式单步法应用于方程时6/6从而误差传播方程为解不等式得故，该方法应用到所给方程的绝对稳

6、定区间为。三、用Taylor展开原理构造形如的两步法，试确定系数使方法具有二阶精度，并推导其局部截断误差主项。分析：本题考察构造多步法的方法；二阶精度即局部截断误差为。解：由多步法局部截断误差的定义将，和分别在点展开得：将此三式代入局部截断误差(3>式，我们有要使方法具有二阶精度，必需6/6解得。此时，所以局部截断误差主项为。三、取，试用Euler法求解初值问题。分析：先将二阶方程写为一阶方程组，再用Euler法求解。解：引进，并记向量函数，则原二阶方程变为一阶方程组，它满足初值条件。将Euler法应用到该方程组得：即，用Euler法求得原二阶方程的数值解为。申明：所有资料

7、为本人收集整理，仅限个人学习使用，勿做商业用途。6/6