常微分方程讲解.pptVIP

第九章 常微分方程数值解 * * * 第九章 常微分方程数值解法 §1 Euler（折线）法 1.Euler法 2.改进Euler法 3. 预估—校正（ PC）法 §2 Runge—Kutta法 1.二级R—K法 2.三级R—K法 3.四级四阶R—K法 对于常微分方程初值问题 则（9.1）在区间[a,b]上存在唯一解 y=y(x). 如果f(x,y) 在 [a,b]×(-∞,+∞)上连续，且关于 y 满足 Lipschtz 条件： （9.1） |f (x,y1) –f(x,y2)|≤L|y1-y2|（常数L0） （9.2） 对于(9.1)在区间[a,b]上的唯一解 y=y(x)，一般情况下很难求出其解析解，因此只能通过数值解法求其近似解。 也就说，构造适当的数值方法，利用(9.1)求出 y=y(x) 在离散节点 x1,x2, …, xn 处的近似函数值 y1 , y2 , …, yn 。 常用的两种方法：Euler 折线法和 Runge-Kutta 法。 x∈ [a,b] §1 欧拉(折线)法 一、Euler 法 xi=x0+ih , i=0,1,2, …，n 对于初值问题 将区间[a,b] n 等分，步长为 h=(b-a)/n ，得到n+1个分点 已知 y=y(x) 在