第九章-常微方程数值解法 1.pptVIP

第章 第9章 常微分方程数值解法 第9章目录 第8章 序 第9章 序 微分方程 数值解 §1 欧拉（Euler）法 1.1 Euler法及其简单改进 Euler公式的推导(续1） Euler公式的推导（续2） Euler公式的推导（续4） Euler公式的推导（续5） 2）以右矩形面积代替曲边梯形(后退的梯形公式）：如图8-3 Euler公式注释 Euler公式注释（续） 关于Euler法的整体截断误差注释 结 论 结 论（续） 1.2 改进的Euler法 改进的Euler法(续） 预测──校正型公式 改进的Euler公式的局部截断误差分析 改进的Euler公式举例 表格8-1 表格8-2 表8-2计算结果说明（续） §5 收敛性与稳定性 5.1 收敛性 5.2 稳定性 稳定性（续1） 第9章 常微分方程数值解法 第9章 常微分方程数值解法 W Y §1 欧拉（Euler）方法 1.1 Euler法及其简单改进 1.2 改进的Euler法 §2 龙格库塔（Runge-kutta）方法 2.1 龙格-库塔方法的基本思想 2.2 二阶龙格-库塔公式 2.3 高阶R-K公式 2.4 变步长R-K法 §3 线性多步法 §4 一阶方程组与高阶方程初值问题 §5 收敛性与稳定性 许多科学技术问题，例如天文学中的星体运动，空间 技术中的物体飞行，自动控制中的系统分析，力学中的振 动，工程问题中的电路分析等，都可归结为常微分方程的 初值问题。 所谓初值问题，是函数及其必要的导数在积分的起始 点为已知的一类问题，一般形式为： 我们先介绍 简单的一阶问题： 由常微分方程的理论可知：上述问题的解唯一存在。 常微分方程求解求什么？应求一满足初值问题（8—1） 的解函数y = y(x) ，如对下列微分方程： 《高等数学》中，微分方程求解，如对一阶微分方程： y? =f(x,y)是求解解函数y = y(x) ，使满足上述方程。但能够 求出准确的解析函数y(x)的微分方程是很少的，《高数》 中研究微分方程的求解，是分门别类讨论，对不同类型的 微分方程，求解方法不一样，因此，要求解微分方程，首 先必须认清类型。 而常微分方程 初值问题的数值解法，是要寻求解函数y(x) 在一系列点y(xi ) （离散点）: 上 y(xi )的近似值yi（ i=1,2,…,n），并且还可由这些（xi,yi) （i=1,2,…,n）构造插值函数作为近似函数。上述离散点相 邻两点间的距离hi=xi-1-xi 称为步长，若hi 都相等为一定数 h, 则称为定步长，否则为变步长。 由于在实际问题和科学研究中遇到的微分方程往往很 复杂，绝大多数很难，甚至不可能求出解析函数y(x)，因 此只能考虑求其数值解。 本章重点讨论如下 一阶微分方程： 在此基础上介绍一阶微分方程组与 高阶微分方程的数值解法。 以Euler法及其改进方法为例,说明 常微分方程初值问题数值解法的一般概 念，Euler法很简单，准确度也不高， 介绍此方法的目的，是由于对它的分析 讨论能够比较清楚地显示出方法的一些 特点，而这些特点及基本方法反映了其 它方法的特点。 Euler法用于求 解一阶微分方 程初值问题： Euler公 式为： 由x0出发?x1,x2,…,xN，而利用此式可算出对应的 y1,y2,…,yN，式（8-2）称为差分方程（序列{yn}满足的方 程）。 下面是Euler公式的推导 ： 一、从几何意义出发：y?=f (x,y)的解函数y=y (x) 在xoy平 面上是一族解曲线， 而初值问题则是其中一条积分曲线假定y = y(x)的曲线如图8-1从给定的点P0(x0,y0) 出发，以P0为切点，作切线，切线斜率为曲线y(x)的切线斜率 y? =f (x0,y0)，因此可 得切线：（点斜式） P1 P2 y(x) P0 x2 x1 x0 几何意义：用折线近似解曲线y = y(x)，折线不会偏离太远 ，因为每项以f (x, y)（斜率）修正。 切线与x = x1交于P1(x1,y1)，在[x0,x1]上以切线 近似曲线， 二、利用Taylor级数：将y(x)在xn处展开： 四、利用数值积分公式：在[xn,xn+1]上对y?(x)=f (x,y(x)) 积分 对右端积分项采用不同的数值积分公式，便可得到各种 不同的求解dE初值