2
f (xi , yi )
0
言
Numerical Analysis
M. H. Xu
二 理论基础 定理 7. 1 若f (x, y )在区 域 D = {(x, y )|a ≤ x ≤ b, |y | < ∞} 内连 续, 关于 y 满足 Lipschitz条件 , 即 存在 常数 L > 0, 对 任意 y1 , y2 , 不 等 式 |f (x, y1 ) − f (x, y2 )| ≤ L|y1 − y2 | 对所 有的 x ∈ [a, b]成立 , 则 初值 问题 dy = f (x, y ), a ≤ x ≤ b dx y (a) = y
Numerical Analysis
M. H. Xu
§7. 2 Euler方法与梯 形方 法 一 方法导出 由微分方程知 y (xi ) = f (xi , y (xi )), 用差商
y0 y yi y=y(x) yi+1 yn
y (xi+1 ) − y (xi ) h 近似 导数 y (xi )可 得 y (xi+1 ) ≈ y (xi ) + hf (xi , y (xi ))
0
在区 间[a, b]上有 唯 一解 y (x), 并且 y (x)为 连续 可微 的, 解函 数y (x) 连续 地依 赖于 初值 及f (x, y ).
Numerical Analysis
M. H. Xu
三 数值解法的基本步骤 第一 步: 把区 间[a, b]进行 划分 , 通 常进 行n等 分, 节点 xi = a + ih, i = 0, 1, 2 · · · , n, 其中 h = (b − a)/n; 第二 步: 求y (x)在节 点xi 处 函数 值y (xi )的近 似值 yi , 得 一列 表 函数 ; 第三 步: 根据 需要 可由 插值 方法求得 函数 y (x)在 x处的 近似 值, 或 由 列表 函数 求得 y = y (x)的近 似函 数. 说明 : 数值 解法 的关 键在 于如 何由 y0 得到y (x1 )的近 似值 y1 , 一 般地 , 如何 由y (xi )的近 似值 yi 得到y (xi+1 )的近 似值 yi+1 .