Ch13A：数值计算方法之常微分方程数值解法.pptVIP

第13章：常微分方程的数值解 在生产实践和科学研究过程中，我们经常要求解常微分方程或常微分方程组的定解问题。在本章中，我们主要研究一阶常微分方程（组）的初值问题，其数学形式为： 为了避免讨论解的存在性和唯一性之类的问题，我们约定f(x,y)是适当光滑的函数，而且关于y满足李普希兹 (Lipshitz)条件，从而保证了解的存在性和唯一性。 作为实际问题，我们还可以假定要求的是某个区间[a,b]上的解，此时我们总是假定x0=a。 13.1 问题求解的基本方法 假如y=y(x)为问题的理论解，如果无法求得，或者所得结果不是初等函数，我们可以直接求它的数值解，也就是求出y=y(x)在一些离散点处的近似值。 求数值解的基本步骤是： 第一步：将问题离散化； 第二步：建立递推格式； 第三步：按步进方式前进。 其中最主要的工作还是建立递推格式，基本目的是减小计算量，提高计算精度。 1.将问题离散化 求微分方程的数值解，实际就是求理论结果y=y(x)在[a,b]上的一些离散点x0,x1,…,xn处的函数值y(x0),y(x1),…,y(xn)的近似值y0,y1,…,yn,接下来就可以利用插值的方法得到相应的插值函数或插值多项式作为问题的近似解。 为了将问题适当简化，我们总是把区间划分为n等份，并记 x0=a,xn=b,h=(b-a)/n,xk=x0+kh,k=0,1,…,n. 作为变通方法，我们也可以先确定步长h和步数n,然后得到b=x0+nh. 2.建立递推关系 最简单的递推关系就是yk+1与yk之间的函数关系，由此所形成的方法称为单步法；较复杂一些的就是yk+1与yk-r,yk-r+1,…,yk之间的函数关系，由此所形成的方法称为多步法。 通常我们用下面一些方法来建立递推关系: 利用差商替代微商 利用求积公式替代积分 利用泰勒展式 下面分别介绍这些方法。 (1) 利用差商替代微商 利用 在这个公式中，我们可以直接利用xk,yk计算yk+1,所以称为显式格式。 （2）利用求积公式替代积分 注意到问题在区间[xk-1,xk]上的理论解可以表为 我们这里得到的yk+1是一个非线性方程的解，或者说是一个隐函数，我们称为是一个隐式格式。直观感觉是隐式格式更精确一些，但要解一个非线性方程，好在并不难。 （3）利用泰勒展式 利用xk+1=xk+h,我们可以把y(xk+1)在xk处展成泰勒级数，从而有 利用泰勒展式可以得到更为准确的解，不过公式推导却非常麻烦，这个问题留到后面再详细讨论。 3.算法框架 为了计算诸y(xk)的近似值yk,我们可以利用y(x0)=y0以及所建立的yk+1与yk之间的递推关系，依次令k=0,1,2,…从而可以依次计算出y1,y2,y3,…这就是按步进的方式进行计算的含义，也可以看成是一个算法框架。 求微分方程的数值解的算法框架都是相同的，不同方法的差别只在于所引用的递推形式不同，所以我们很容易对不同的方法进行比较。 术语：人们通常把yk+1,与yk, yk-1,…之间的递推关系，或者说函数关系，称为差分方程。 4.记号约定 为了求问题（1）在某个区间[a,b]上的数值解，我们通常根据精度要求把区间[a,b]划分为某个n等分,记h=(b-a)/n，称之为步长；我们总是约定x0=a,xk=x0+kh,k=0,1,…,n,它们都是精确值； 记y=y(x)为问题的解，尽管有时候我们写不出它的表达式来，但总是有定义的；记y(xk)表示当x=xk时y(x)的精确值；记yk为我们采用某个算法求得的y(xk)的近似值，记fk=f(xk,yk)。 提示：对于上面的记号，我们一般地有yk≠y(xx)，从而还有fk≠f(xx,y(xk))。 5.局部截断误差与整体截断误差 假如我们利用递推格式yk+1=yk+hf(xk,yk)或其他形式的递推关系来计算yk+1的值，那么它与真值y(xk+1)之间的误差既与步长h有关，也与yk的值有关。为了误差分析的方便，我们把仅与步长h有关的误差称为局部截断误差，而与yk和h都有关的误差称为整体截断误差。简言之，记 为局部截断误差和整体截断误差。 6.收敛性定义 假如我们求微分方程在某个区间[a,b]上的数值解，记x0,x1,…,xn为计算基点，其中x0=a,xk=a+kh,,k=0,1,…,n, 记ε1(h),ε2(h),…,εn(h),和e1(h),e2(h),…,en(h)分别为y1,y2,…,yn的整体截断误差和局部截断误差，如果 则称数值解法是收敛的。 注释：在一般情况下，我们所讨论的数值方法对于我们所假定的问题都是收敛的，所以在我们的课程中并不专门讨论数值解法的收敛性问题