常微分方程数值解法.pptVIP

常微分方程数值解法第一页，共七十二页，2022年，8月28日 §1 引言1.0 基本概念1. 常微分方程的初值问题：称为具有初值(1.2)的常微分方程.① 若f(x,y)在{a?x?b, |y|+?}上连续，且关于y满足Lip条件：?常数L使| f(x, y1) –f(x, y2)| ? L|y1 – y2|则初值问题(1.1)(1.2)存在唯一连续可微解y(x).注：以下总假设f 满足Lip条件.第二页，共七十二页，2022年，8月28日 §1 引言1.0 基本概念1. 常微分方程的初值问题：称为具有初值(1.2)的常微分方程.② (1.1)(1.2)等价于微分方程： (1.3)注：一般无初等解(解析解)，即使有形式也复杂.第三页，共七十二页，2022年，8月28日 §1 引言1.0 基本概念2. 初值问题的数值解 设(1.1)(1.2)的解y(x)在节点xi处的近似解值为 yi ? y(xi), a x1 x2 … xn = b则称yi (i = 1, 2, …, n)为(1.1)(1.2)的数值解，又称y(xi)的计算值.第四页，共七十二页，2022年，8月28日 §1 引言1.0 基本概念3. 数值方法① 两种转化： 由微分出发的数值方法. 由积分 出发的数值方法.② 计算方法 步进法：从初始条件出发，逐步求y1, y2, …, yn. 又有两种：单步法，多步法.注：采用等距节点：第五页，共七十二页，2022年，8月28日 §1 引言1.1 基于数值微分的求解公式. (1.6)第六页，共七十二页，2022年，8月28日 §1 引言1.1 基于数值微分的求解公式.1. 前进欧拉公式 (1.6)的前半部分为：令 yi+1 = yi + hf(xi, yi) (1.7)其中yi = y(xi) , 则yi+1 ? y(xi+1)第七页，共七十二页，2022年，8月28日 §1 引言1.1 基于数值微分的求解公式.1. 前进欧拉公式 令 yi+1 = yi + hf(xi, yi) (1.7)其中yi = y(xi) , 则yi+1 ? y(xi+1)记 (1.8)则称(1.7)为前进欧拉求解公式. 简称为欧拉公式或欧拉法. (1.8)称为欧拉公式的余项：ei+1(h) = y(xi+1) –yi+1 第八页，共七十二页，2022年，8月28日 §1 引言1.1 基于数值微分的求解公式.2. 后退欧拉公式 (1.6)的后半部分令 yi+1 = yi + hf(xi+1, yi+1) (1.9)其中yi = y(xi), 则yi+1 ? y(xi+1) 第九页，共七十二页，2022年，8月28日 §1 引言1.1 基于数值微分的求解公式.2. 后退欧拉公式令 yi+1 = yi + hf(xi+1, yi+1) (1.9)其中yi = y(xi), 则yi+1 ? y(xi+1) 注：①(1.9)中f(xi+1, yi+1) ? f(xi+1, y(xi+1)) ∴余项 (1.10)第十页，共七十二页，2022年，8月28日 §1 引言1.1 基于数值微分的求解公式.2. 后退欧拉公式令 yi+1 = yi + hf(xi+1, yi+1) (1.9)其中yi = y(xi), 则yi+1 ? y(xi+1) 注：② 称(1.9)为后退欧拉公式(后退欧拉法). 称(1.10)为后退欧拉法的误差近似值.③ 欧拉法与后退欧拉公式的区别：(1.7)为直接计算公式称显式公式.(1.9)