重庆大学《数值分析》课件-第9章.pdf

第九章常微分方程初值问题的数值解法

§9.1引言

科学研究和工程技术中的问题往往归结为求某个常微分方

程的定解问题.

常微分方程的理论指出，很多方程的定解虽然存在，但可

能十分复杂难于计算，也可能不能用简单的初等函数表示,因

此常求其能满足精度要求的近似解.

常微分方程的数值解法常用来求近似解，由于它提供的算

法能通过计算机便捷地实现，因此近年来得到迅速的发展和

广泛的应用常微分方程数值解法的特点是：对求解区间进行

剖分，然后把微分方程离散成在节点上的近似公式或近似方

程，最后结合定解条件求出近似解.因此数值解法得到的近似

解是一个离散的函数表.

1结束

9.1.1基本知识复习

一阶常微方程的初值问题：

yf(x,y)axb(1)



y(a)y0(2)

其中f(x,y)是已知函数，(2)是定解条件，常微分方程的理论指

出：当f(x,y)定义在区域G(a≤x≤b,｜y｜＜∞)，若存在正的常

数L，使：

｜f(x,y)–f(x,y)｜≤L｜y-y|

1212

对任意的x∈[a,b]和任意的y,y成立，则称f(x,y)对y满足李普

12

希兹(Lipschitz)条件,L称李普希兹常数.若f(x,y)在区域G连续，

关于y满足李普希兹条件，则上述一阶常微分方程的初值问题

的解存在且唯一.我们以下的讨论，都在满足上述条件下进行.

2结束

9.1.2一阶常微分方程组和高阶常微分方程

一阶常微分方程组常表述为：

y1f1(x,y1,,ym)





yf(x,y,,y)

mm1m

(axb)

y1(a)1





ym(a)m



写成向量形式为

yf(x,y)axb

T

y(a)y,y(,,)

001m

3结束

还有高阶常微分方程定解问题，如二阶定解问题：

yf(x,y,y)



y(a)



y(a)

可化作一阶常微方程组的定解问题：

yu(x,y)



uf(x,y,u)



y(a),u(a)

本章着重讨论一阶常微方程初值问题的各种数值解法.这

些解法都可以写成向量形式而用于一阶常微分方程组的初值

问题.这样，也就解决了高阶方程的定解问题.

4结束

§9.2欧拉(Euler)方法

对问题

yf(x,y)