数值分析课件-第八章 常微分方程数值解.pdfVIP

第八章 常微分方程数值解 8.1 引言(基本求解公式) § 在工程和科学技术的实际问题中，常需要求解微分方程 只有简单的和典型的微分方程可以求出解析解 而在实际问题中的微分方程往往无法求出解析解 在高等数学中我们见过以下常微分方程： y f (x,y) axb y(a)y0 (1) y f (x,y,y ) axb    (2) y(a)y ,y (a)  0 y f (x,y,y ) axb   (3) y(a)y ,y(b)y0 n (1),(2)式称为初值问题,(3)式称为边值问题 另外,在实际应用中还经常需要求解常微分方程组: y f (x,y ,y ) y (x )y  1 1 1 2 1 0 10 (4) y f (x,y ,y ) y (x )y  2 2 1 2 2 0 20 本课程主要研究问题(1)的数值解法,对(2)~(4)只作简单介绍 我们首先介绍初值问题(1)的解存在的条件 定理1. 如果连续函数f (x,y)满足Lipschitz条件，即 正数L, 使得x[ab, ]，均 | f (x,y ) f (x,y ) |L |y y | 1 2 1 2 则初值问题(1)的解存在且唯一. 对于问题(1),要求它的数值解 就是求未知函数y(x)在区间[a,b]上的一系列离散点(节点) ax x x x b 0 1 2 n 上函数值y(x )的近似值y (k1,2,,n)k k 而y (k1,2,,n)就是问题(1)的数值解 k 从(1)的表达式 y f (x,y) axb  (1) y(a)y0 可以看出,求它的数值解的关键在于 y (x)数值计算问题 x 或者它的等价的积分方程 yx( ) y  f (t, yt( ))dt 中 0 a x 积分 a f (t,yt( ))dt 的数值计算问题 数值微分 求解微分方程的数值方法 数值积分 而数值积分问题我们已经学习过, 下考虑数值微分方法 微积分中，关于导数的定义如下： f (x h) f (x) f (x) lim