第十常微分方程数值解.ppt
《第十常微分方程数值解.ppt》由会员分享,可在线阅读,更多相关《第十常微分方程数值解.ppt(34页珍藏版)》请在三一文库上搜索。
1、第十章 常微分方程数值解,第一节 求解初值问题数值方法的基本原理,第二节 高精度的单步法,第三节 线性多步法,第四节 一阶微分方程组的解法,第五节 边值问题的打靶法和差分法,考虑一阶常微分方程的初值问题 /* Initial-Value Problem */:,只要 f (x, y) 在a, b R1 上连续,且关于 y 满足 Lipschitz 条件,即存在与 x, y 无关的常数 L 使 对任意定义在 a, b 上的 y1(x) 和 y2(x) 都成立,则上述IVP存在唯一解。,要计算出解函数 y(x) 在一系列节点 a = x0 x1 xn= b 处的近似值,节点间距 为步长,通常采用等
2、距节点,即取 hi = h (常数)。,第一节 求解初值问题数值方法的基本原理,数值解,(10-1),一、初值问题的数值解,求解(10-1)最基本的方法是单步法,单步法:从初值 开始,依次求出 ,后一步的值 只依靠前一步的,典型的单步法是Euler(欧拉)方法,其计算格式是:,例:求解常微分方程初值问题,由此可见,Euler公式的近似值接近方程的精确值.,二、构造初值问题数值方法的基本途径,以Euler法为例说明构造IVP问题数值方法的三种基本途径,1. 数值微分法,用差商代替微商,亦称为欧拉折线法,2. Taylor展开法,忽略高阶项,取近似值可得到Euler公式,3. 数值积分法区间,将
3、区间 积分,隐式欧拉法 /* implicit Euler method */,由于未知数 yi+1 同时出现在等式的两边,不能直接得到,故称为隐式 /* implicit */ 欧拉公式,而前者称为显式 /* explicit */ 欧拉公式。,一般先用显式计算一个初值,再迭代求解。,三、Euler法的改进及梯形公式,梯形公式 /* trapezoid formula */, 显、隐式两种算法的平均,中点欧拉公式 /* midpoint formula */,改进欧拉法 /* modified Eulers method */,注:此法亦称为预测-校正法 /* predictor-corre
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第十 微分方程 数值
链接地址:https://www.31doc.com/p-2565358.html