常微分方程初值问题9.pptVIP

北京科技大学数理学院

卫宏儒

;

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

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

;

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

反之，假设y(x)满足积分方程〔1.4〕，可验证它满足〔1.1〕和〔1.2〕，所以〔1.4〕式与初值问题〔1.1〕，〔1.2〕等价,这说明可用积分方程构造初值问题的数值解法。;;为了研究数值方法的绝对稳定性，下面给出常系数线性差分方程的有关概念。;;讨论线性多步法的绝对稳定性条件;其对应的特征方程为：;;2、数值解法的构造途径;它称为〔向前〕欧拉(Euler)公式。〔类似地可以用向后差商、中心差商代替导数产生相应的欧拉(Euler)公式〕;〔2〕数值积分法

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

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

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

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

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

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

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

改进欧拉法的预-校公式 ;Euler公式的几何意义;例题：用Euler公式和改进的Euler公式分别求以下初值问题的数值解(取步长h=0.1计算到y3)：

y′=-2xy2

y(0)=1

解：由欧拉公式

yn+1=yn+hf(xn,yn)=yn-2hxnyn2

计算如下

y1=y0-2hx0y02=1-2·0.1·0·12=1

y2=y1-2hx1y12=1-2·0.1·0.1·12=0.98

y3=y2-2hx2y22=0.98-2·0.1·0.2·0.982=0.9416;用改进欧拉法的预-校公式计算如下：;3、Runge-Kutta方法;对照标准形式yn+1=yn+hф(xn,yn,h)。假设取

ф(x,y,h)=y(x)+(h/2!)y(x)+......+(hp-1/p!)y(p)(x)

并以yn代替y(xn)，那么得到一个p阶近似公式

yn+1=yn+hф(xn,yn,h)

(n=0,1,2,......)(*);R-K方法不是直接使用Taylor级数,而是利用它的思想,即计算f(x,y)在不同结点的函数值,然后作这些函数值的线性组合,构造近似公式,式中有一些可供选择的参数。将近似公式与Taylor展开式相比较,使前面的假设干项密合,从而使近似公式到达一定的精度。

下面以二级二阶R-K方法为例说明这一方法的根本思想。;

在[xn,xn+1]上,取f(x,y)在两个点的函数值作线性组合,即得到二级R-K方法:

yn+1=yn+h(c1K1+c2K2)

K1=f(xn,yn) 〔**〕 K2=f(xn+a2h,yn+b21hK1)

其中c1,c2,a2,b21为待定参数。对照式(*)有：

ф(x,y,h)=c1f(x,y)+c2f(x+a2h,y+b21hf(x,y)); ф(xn,y(xn);h)=(c1+c2)y(xn)+c2(a2hfx+b21hfyf)+O(h2) 因为y(xn+1)在xn处的Taylor展开为

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

由显式单步法在xn+1的局部截断误差定义有：

Tn+1=y(xn+1)-y(xn)-hф(xn,y(xn),h)

=h(1-c1-c2)y(xn)+h2[(1/2-a2c2