数值分析 李庆杨 Cht9常微分方程初值问题数值解法.pdfVIP

第9章 常微分方程初值问题数值解法 §1 引 言 考虑一阶常微分方程初值问题 dy  f (x ,y ) (1.1) dx  y (x ) y (1.2)  0 0 在x0  x  X 上的数值解法。 假定f (x ,y )在R {x0  x  X ,  y  }中连续, 并且 关于y 满足Lipschitz条件： | f (x ,y ) f (x ,y ) | L | y y |, (x ,y ), (x ,y ) R ， 1 2 1 2 1 2 其中L  0为Lipschitz常数,则该初值问题的解在x0  x  X 上存在且唯一，而且连续依赖于初始条件及右端函数。 解析方法与数值解法. 数值解法就是求y (x)在离散节点x1   xn  上的 近似值y ,,y ,. 1 n 常用步进式的方法. 单步法，k步法. 研究计算公式的局部截断误差和阶，数值解y 与准确解 n y (x )的误差估计及收敛性，计算公式的稳定性等. n §2 简单的数值方法与基本概念 一、欧拉法和后退欧拉法 积分曲线. 方向场. 欧拉(Euler)公式： y n1 y n  hf (xn ,y n ), (n 0,1,2,). (2.1) 例1 取h 0.1,利用Euler 公式求解   2x y y  ,（0  x  1)   y (2.2)  y (0) 1,  2x 0.2x 解：y n1 y n  h(y n  n ) 1.1y n  n . 计算及结果如下 y n y n y = 1 x = 0 clear; y=1, x=0, %初始化 y = 1.1000 x = 0.1000 for n=1:10 y = 1.1918 x = 0.2000 y=1.1*y-0.2*x/y, x=x+0.1, y = 1.2774 x = 0.3000 end y = 1.3582 x = 0.4000 y = 1.4351 x = 0.5000 几何意义