第9章普通微分方程的数值解.pptVIP

第十章;问题驱动：蝴蝶效应; 该方程组来源于模拟大气对流，该模型除了在天气预报中有显 著的应用之外，还可以用于研究空气污染和全球侯变化。洛伦 兹借助于这个模型，将大气流体运动的强度x与水平和垂直方 向的温度变化y和z联系了起来。参数 ;§1 引 言;一、初值问题的数值解法;(2) 一般构造方法： 离散点函数值集合 + 线性组合结构 → 近似公式;(3) 如何保证迭代公式的稳定性与收敛性?;称 在区域D上对 满足Lipschitz条件是指:;二、初值问题解的存在唯一性; 求函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 的方法称为微分方程的数值解法。;三、初值问题的离散化方法;§2 欧拉方法 /* Euler’s Method */;亦称为欧拉折线法 /* Euler’s polygonal arc method*/ ;? 欧拉法的局部截断误差：;例1: 用欧拉公式求解初值问题; ;可用来检验近似解的准确程度。;隐式欧拉法(向后Euler法） /* implicit Euler method */;注：由于未知数 yi+1 同时出现在等式的两边，不能直接得到，故称为隐式 /* implicit */ 欧拉公式，而前者称为显式 /* explicit */ 欧拉公式。;一般先用显式计算一个初值，再迭代求解。;? 梯形公式 /*trapezoid formula */;? 中点欧拉公式 /* midpoint formula */;方 法;? 改进欧拉法 /* modified Euler’s method */;注:此法亦称为预测-校正法 /* predictor-corrector method */;在实际计算时，可将欧拉法与梯形法则相结合， 计算公式为 ;改进欧拉法的截断误差;二元函数的n阶Taylor展式：;例2: ;改进的Euler法;通过计算结果的比较可以看出，改进的Euler方法;§3 龙格 - 库塔法 /* Runge-Kutta Method */;? 考察改进的欧拉法，可以将其改写为：;将改进欧拉法推广为：;Step 2: 将 K2 代入第1式，得到;要求 ，则必须有：; 为获得更高的精度，应该如何进一步推广？;考虑一阶常微分方程初值问题;; 若 ;;n级显式Runge-Kutta方法;二阶Runge-Kutta方法;由此得;另一方面;改进的Euler方法;中点方法;二阶Heun方法;二级Runge-Kutta方法不超过二阶;因此局部截断误差只能达到;三级Runge-Kutta方法;记;溪曰饰锨熙肘蟹乙渠结脊绵厂断袋薪炒月幼陨摇叮之疼亭檬枷毖确焚技浑第9章普通微分方程的数值解第9章普通微分方程的数值解;又由于;因此要使局部截断误差为 ，必须 ;Kutta方法;三阶Heun方法;三级Runge-Kutta方法不超过三阶; 四级R-K方法;熟揣汀尘氦簇诚菊芹酥克淤瘤西咐港纪却捶惧吼促梯多行盔疽却孕娠拱橱第9章普通微分方程的数值解第9章普通微分方程的数值解;局部截断误差为O(h5);附注：;注： ? 龙格-库塔法的主要运算在于计算 的值，即计算 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系：;§4 单步方法的收敛性与稳定性 /* Convergency and Stability */;解：该问题的精确解为 ;? 稳定性 /* Stability */;定义;我们称算法A 比算法B 稳定，就是指 A 的绝对稳定 区域比 B 的大。;例：考察显式欧拉法;例：考察隐式欧拉法;例：隐式龙格-库塔法;而显式 1~ 4 阶方法的绝对稳定区域为