[理学]常微分方程的数值解法.docVIP

第六章 常微分方程的数值解法 §6.0 引言 §6.1 算法构造的主要途径 §6.2 Runge-Kutta Method算法 §6.3 线性多步法 §6.4 线性多步法的一般形式 §6.5 算法的稳定性、收敛性 §6.0 引 言 1． 主要考虑如下的一阶常微分方程初值问题的求解: 微分方程的解就是求一个函数y=y(x)，使得该函数满足微分方程并且符合初值条件。 例如微分方程: xy-2y=4x ；初始条件: y(1)=-3。 于是可得一阶常微分方程的初始问题 。 显然函数y(x)=x2-4x满足以上条件,因而是该初始问题 的微分方程的解。 但是,只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解，有的甚至无法用解析表达式来表示。因此，只能依赖于数值方法去获得微分方程的数值解。 4． 微分方程的数值解：设微分方程问题的解y(x)的存在区间是[a,b]，初始点x0=a，将[a,b]进行划分得一系列节点x0 , x1 ,...,xn，其中a= x0 x1… xn =b。 y(x)的解析表达式不容易得到或根本无法得到，我们用数值方法求得y(x)在每个节点xk的近似值y(xk)，即y≈y(xk)，这样y0 , y1 ,...,yn称为微分方程的数值解。如图所示: §6.1 算法构造的主要途径 1 欧拉公式 1.1 构造的思想： 微分方程初值问题: 利用差商代替一阶导数，即 ， 则 。 于是，可求出y(x1)的近似值y1， 同样地，可利用x1处的微分方程可得： 一般地，利用在xn处的微分方程可得： 此式称为欧拉公式。 1.2 几何意义： 对于微分方程y=2(x+1),其通解是y=(x+1)2+c,是一个曲线族，当给定初值条件y(0)=2，其特解为y=(x+1)2+1。 如图所示： 由y(0)=2，过该曲线上一点(0,2)作曲线的切线， 其斜率， ∴切线为： , 因此可计算出y1, 如此，可根据：，… 故欧拉法又称欧拉折线法。 1.3 算例： 例： 解：h=0.2 , xi=1+ih y1= y0+hf(x0，y0)=-1+0.2× y2= y1+hf(x1，y1)=-1+0.2× y3= y2+hf(x2，y2)= -0.9333+0.2× y4= y3+hf(x3，y3)=0.8+0.2× y5= y4+hf(x4，y4)=0.6+0.2× y6= y5+hf(x5，y5)=0.3333+0.2× 精确解为：y=x2-2x xk y(xk) yk ek 1.2 -0.96 -1 0.04 1.4 -0.84 -0.9333 0.0933 1.6 -0.64 -0.8 0.16 1.8 -0.36 -0.6 0.24 2.0 0 -0.3333 0.3333 2.2 0.44 0 0.44 可以看出误差随着计算在积累。 1.4 Euler法的特点和误差 迭代格式： 特点： (1) 单步方法； (2) 显式格式； (3) 局部截断误差为。 局部截断误差： 当时，由按照欧拉方法计算来的的 误差称为局部截断误差。即， 是局部截断误差。 如： 欧拉法得： 因此，局部截断误差是。 2 改进Euler法 2.1方法构造 在微分方程初值问题,对其从到进行定积分得: 将右端的定积分用梯形公式来进行近似计算。 用和来分别代替和得计算格式: 这就是改进欧拉方法。 2.2 显式格式和隐式格式 在欧拉式中 每一步计算已知,直接用格式可以计算出,此类格式称为显式格式。 而在改进欧拉方法中 在每一步计算中是未知,待求的,未知量在中这是一个方程,如是非线形或超越函数,此方程是无法直接解出来(要依靠迭代法才能解出)。这类格式称为隐式格式。 2.3 算例 例: 用改进欧拉方法求解。 解：解得: 注意：由于,是线形函数可以从隐式格式中 解出。 问题的精确解是 欧拉方法 误差 … … … … … 3 预测－校正方法 由于改进Euler法是隐式格式，无法从格式中直接求出必须要解方程。下面用预测－校正方法来求隐式格式中的。 预测值: 校正值: 此式相当于对隐式格式求时采用迭代的方法,用欧拉格式得到的作为初始值迭代公式迭代一次而已,此公式代入后得: 如改写成平均的形式为: §6.2 龙格-库塔法Runge-Kutta Method 1 龙格-库塔法的思想 1.1 考虑微分方程的初值问题: 根据微分中值定理有： ，其中0θ1。 于是 即 我们称为y(x)在区间[xk, xk+1]上的平均斜率，记作K，其中，θ是存在但是未知的