第8章 常微分方法的数值解法 ----8.1 Euler 方法.PDFVIP

第八章常微分方程数值解法 第8章常微分方法的数值解法 8.1 Euler 方法 8.1.1 Euler 方法及其有关的方法 8.1.2 局部误差和方法的阶 第八章常微分方程数值解法 第8章常微分方法的数值解法 教学目的 1. 掌握解常微分方程的单步法：Euler方法、 Taylor方法和Runge-Kutta方法； 2. 掌握解常微分方程的多步法：Adams步法、 Simpson方法和Milne方法等； 3. 了解单步法的收敛性、相容性与稳定性；多步 法的稳定性。 教学重点及难点 重点是解常微分方程的单步法：Euler方法、 Taylor方法和Runge-Kutta方法和解常微分方程的 多步法：Adams步法、Simpson方法和Milne方法 等； 难点是理解单步法的收敛性、相容性与稳定性及多 步法的稳定性。 第八章常微分方程数值解法 第8章常微分方法的数值解法 科学技术与工程问题常常需要建立微分方程形式的数学模型， 下面是这类问题的例子。 设N （t ）为某物种的数量， 为该物种的的出生率与死亡率之 α β 差， 为生物的食物供给及它们所占空间的限制，描述该物种增长 率的数学模型是 dN 2 α β N (t ) − N (t ) dt N (t ) N 。 0 0 设Q是电容器上的带电量，C为电容，R为电阻，E为电源的 电动势，描述该电容器充电过程的数学模型是 dQ Q (t ) E − ， dt RC Q (t 0 ) Q0 。 第八章常微分方程数值解法 以上两个例子是常微分方程初值问题，下面是一个两点边值问 题的例子。 设一跟长为L的矩形截面的梁，两端固定。E是弹性模量，S是端 点作用力，I （x ）是惯性矩，q是均匀荷载强度，梁的桡度y （x ）满 足如下方程 d 2 y S qx 2 y (x ) + (x −l )， dx EI (x ) 2EI (x ) ( ) y （0） y L 0。 针对实际问题建立的数学模型，要找出模型解的解析表达式往往 是困难的，甚至是不可能的。因此，需要研究和掌握微分方程的数值 解法，即计算解域内离散点上的近似值的方法。本章讨论常微分方程 数值解的基本方法和理论。 第八章常微分方程数值解法 8.1 Euler 方法 8.1.1 Euler 方法及其有关的方法 考虑一阶常微分方程初值的问题： ⎧y f ( x , y )