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

常微分方程初值问题的数值解法;;关于初值问题解的存在，唯一及对初始条件的连续依赖性，有以下定理：

定理1.1设f(x,y)是在D={(x,y)|a≤x≤b,y∈R}上的连续函数，其中a,b为有限实数，而且f(x,y)满足对y的lipschitz条件，那么对(x0,y0)∈D,初值问题〔1.1〕，〔1.2〕在[a,b]存在唯一的连续函数微分解y(x).

定理1.2设f(x,y)是在D={(x,y)|a≤x≤b,y∈R}上的连续函数，满足lipschitz条件〔1.3〕并设y(x;s)是初值问题

dy/dx=f(x,y)，y(x0)=S

的解，那么当s取s1和s2时，有估计

｜y(x;s1)-y(x;s2)｜≤EXP〔L｜x-x0｜｜s1-s2｜)

定理1.2的意义：假设初值问题〔1.1〕，〔1.2〕中，初始值有一微小扰动，那么解的扰动也是微小的，也即解连续依赖初始条件.;

假设y(x)是式〔1.1〕，〔1.2〕的解，从方程〔1.1〕两边积分，再利用式〔1.2〕可得积分方程

反之，假设y(x)满足积分方程〔1.4〕，可验证它满足〔1.1〕和〔1.2〕，所以〔1.4〕式与初值问题〔1.1〕，〔1.2〕等价,这说明可用积分方程构造初值问题的数值解法。;

求初值问题，是给出它的解在某些节点数值的近似值，这称为数值离散方法，假设要求出在[a,b}上的离散解，引入点列{xK}，其中

xk=xk-1+hk-1，k=1，2，3…..

xk称为节点，hk称为步长。

通常，步长h不变，取为等距步长hk=h=〔b-a〕/N，N为等份区间[a,b]分割数,节点变为等距节点，此时有

xk=xk-1+h=x0+kh,

记式〔1.1〕，〔1.2〕的准确解y(x)在节点xk之值为y(xk)，而记y(xk)的近似值为yk，又记fk=f〔xk，yk〕，

;求数值解一般逐步进行，有单步法与多步法之分：

单步法：在计算yk+1之时只用到yk

多步法：在计算yk+1之时不仅用到yk，还要用yk-1，yk-2，…

一般m步法要用到yk,yk-1,yk-2，…yk-m+1，二者都有显式方法和隐式方法之分.

单步显式计算公式可写成:

yk+1=yk+hΦf〔xk，yk;h〕〔*〕

单步法隐式格式yk+1=yk+hΦf〔xk，yk,yk+1;h〕

它每步求解yk+1需要解一个隐式方程。

上式函数Φ是连续函数，称为增量函数。;多步法也可写出类似形式，不过这里只讨论线性多步法。

m步线性多步法的一般格式为：

定义1：设y(x)是初始问题??1.1〕的精确解，yn+1表示在

yn=y(xn)的假设下，由单步显式公式〔*〕算出的数值解，

Tn+1=y(xn+1)-yn+1

=y(xn+1)-y(xn)-h?f(xn,y(xn);h)

称为显式单步法在xn+1的局部截断误差。;; 绝对稳定性与绝对稳定域

求解初值问题的数值方法，当给定不同步长计算时结果的舍入误差影响差异很大，如果舍入误差不增长算法就是数值稳定的，假设舍入误差增长很快算法就不稳定。

定义4：用一个数值方法求解微分方程

y′=?y?是复数(1.5)

对给定的步长h,在计算yn时引起的误差?n，假设这个误差在计算后面的yn+k中所引起的误差?n+k满足:

|?n+k|≤|?n|(k=1,2,…)

就说这个数值方法对步长h和复数?是绝对稳定的,使得数值方法是绝对稳定的H=?h在复平面上的允许范围称为数值方法的绝对稳定域.;2、解一阶微分方程的初值问题的离散化方法;数值积分方法

把y′=f(x,y)在[xk,xk+1]积分，得

对右端的定积分用数值积分方法做离散化，可得计算公式，如用矩形公式可得欧拉公式，假设用梯形公式可得改进的欧拉公式，它也称为梯形公式：;Taylor展开法

设f(x,y)充分光滑,将y(xn+1)在xn点作Taylor展开:

y(xn+1)=y(xn)+hy(xn)+(h2/2!)y(xn)+O(h3)

取其关于h的线性局部，并用yk代替y(xk),就得到Euler公式。

易知Euler公式的局部截断误差为

T1=(h2/2!)y(xn)+O(h3)=O(h2)