计算机方法-常微分方程初值

6.1基本离散方法,6.2Runge-Kutta方法,6.3线性多步法,6.4收敛性与稳定性,6常微分方程数值解法,考虑一阶常微分方程的初值问题:,例如：,其解析解为：,6.1基本离散方法,但是,只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解。因此，只能依赖于数值方法去获得微分方程的数值解。,例如：,其解析解为：,很难得到其解析解,例如：,其解析解为,只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解。因此，只能依赖于数值方法去获得微分方程的数值解。,要计算出解函数y(x)在一系列节点a=x0x1xn=b处的近似值,通常取节点间距为步长，通常采用等距节点，即取hi=h(常数)。它适合计算机求解,应用广泛,具有应用价值。,Taylor级数复习,函数在点作Taylor级数展开：,这里x,x0都可以是任意一点。,Taylor级数复习,则：,6.1.2Euler公式,利用等距分割，数值微分来代替导数项，建立差分格式。,1、向前差商公式,所以，可以构造差分方程,称为局部截断误差。显然，这个误差在逐步计算过程中会传播，积累。因此还要估计这种积累,定义:在假设yi=y(xi)，即第i步计算是精确的前提下，考虑的截断误差Ri=y(xi+1)yi+1称为局部截断误差。,定义:若某算法的局部截断误差为O(hp+1)，则称该算法有p阶精度。,收敛性：,考察局部误差的传播和积累,2、向后差商公式,是隐格式，要迭代求解,可以由向前差商公式求出,3、中心差商公式,是多步，2阶格式，该格式不稳定,对微分方程,积分有:,类似，可以算出其误差估计式：,2阶的方法,所以，有,是个隐式的方法，要用迭代法求解,局部截断误差,4、梯形公式,5、欧拉公式的改进：,隐式欧拉法,由于未知数yi+1同时出现在等式的两边，不能直接得到，故称为隐式欧拉公式，而前者称为显式欧拉公式。,中点欧拉公式（欧拉二步法）,假设，则可以导出即中点公式具有2阶精度。,需要2个初值y0和y1来启动递推过程，这样的算法称为双步法/*double-stepmethod*/，而前面的三种算法都是单步法/*single-stepmethod*/。,简单,精度低,稳定性最好,精度低,计算量大,精度提高,计算量大,精度提高,显式,多一个初值,可能影响精度,改进欧拉法,Step1:先用显式欧拉公式作预测，算出,Step2:再将代入隐式梯形公式的右边作校正，得到,Euler方法、隐式Euler方法、梯形方法与单步法计算公式的,显式单步法,对应关系,隐式单步法,显式Euler方法,隐式Euler方法,梯形方法(隐式),6.1.3总结,算例：分别用Euler公式和改进的Euler公式求解：,取步长，计算y(0.5)的近似值,解：欧拉公式：,改进的Euler公式：,算例分别用显式Euler方法，梯形方法和预估校正Euler方法初值问题,解：,取h=0.1，,(1)Euler方法为:,续,算例分别用显式Euler方法，梯形方法和预估校正Euler方法解初值问题,解：,取h=0.1，,梯形方法为：,续,算例分别用显式Euler方法，梯形方法和预估校正Euler方法解初值问题,解：,取h=0.1，,梯形方法为：,预估校正Euler方法：,续,数值例子表明，梯形方法和预估校正Euler方法比显式Euler方法有更好的精度。,续,基本思想:根据微分中值定理有：,6.2龙格-库塔方法,建立高精度的单步递推格式。,6.2RungeKutta法,由Taylor展开,记为,所以，可以构造格式,这种格式使用到了各阶偏导数，使用不便。,从另一个角度看，,取(x,y)及其附近的点做线性组合，表示F，问题就好办了。当然，要求此时的展开精度相同。这种方法称为RungeKutta法,在(x,y)处展开有,而,以2阶为例，设,比较对应系数,有：,1、改进的Euler公式,2、Heun公式,一般的RungeKutta法构造,常见的为3阶，4阶公式,下列公式是三阶公式中的一个典型例子：,下列公式是经典的四阶R-K方法(古典的R-K方法),算例：分别用Euler公式，改进的Euler公式，经典4阶R-K公式计算一阶常微分方程初值问题。,并与准确解比较。,解：Euler公式,改进的Euler公式取步长h=0.1,经典4阶R-K公式取步长h=0.2。,4阶R-K公式：,计算结果见下表：,6.3线性多步法,思想:用若干节点处的y及y值的线性组合来近似y(xn+1)。,其通式可写为：,当10时，为隐式公式;1=0则为显式公式。,Adams方法是线性多步法的一个代表，它是利用插值多项式进,行积分得出来的，这样构造线性多步法的方法称为数值求积法，,它是构造线性多步法的一种途径，另外还有Taylor法。,构造线性多步法的途径：,1数值积分法,显式Adams方法,隐式Adams方法,其它方法,Simpson方法（二步法）,方程(1.1),由Simpson求积公式，得,用近似，得Simpson方法,局部截断误差为,结论：Simpson方法是四阶隐式方法。,例,Milne法,用过点的插值多项式近似,得Milne方法,局部截断误差为,结论：Milne方法是四阶显式线性多步法。,例：设,确定式中待定系数0,1,2,0,1,2,3,使得公式具有4阶精度。,解：,/*y(xi)=yi*/,2Taylor展开法,解：,个未知数个方程,7,5,令1=2=0,以yi+1取代yi1，并取1=2=0,取1=1,2=0得到辛甫生公式,辛甫生公式,例：设,确定式中待定系数0,1,2,0,1,2,3,使得公式具有4阶精度。,例：确定下列公式：,中的待定系数，使公式具有3阶精度。,由于：,只需：,具有3阶精度，但可以进一步验证公式具有4阶精度,例：确定下列公式：,中的待定系数，使公式具有3阶精度。,6.4收敛性与稳定性,收敛性,例：就初值问题考察欧拉显式格式的收敛性。,解：该问题的精确解为,欧拉公式为,对任意固定的x=xi=ih，有,稳定性,例：考察初值问题在区间0,0.5上的解。分别用欧拉显、隐式格式和改进的欧拉格式计算数值解。,1.00002.50001016.25001021.56251023.90631039.7656104,1.00002.50006.25001.56261013.90631019.7656101,1.00004.97871022.47881031.23411046.14421063.0590107,1.00002.00004.00008.00001.60001013.2000101,一般分