




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、五、线性多步法,第九章 微分方程初值问题的数值解法,1. 线性多步法的一般理论,线性多步法的基本思想: 如果充分利用前面多步的信息来预测 yn+k , 则可以期望获得较高的精度.,记 y(xn)的近似值为 yn , xn = x0+ nh , 并记 fn = f (xn , yn ), 则 k 步线性多步法的一般形式为:,其中 均为常数, 不全为零. 由于上式两端可同乘一常数, 故一般可设 , 则上式可写成:,或写成,注: k=1时, 即为一种单步法. 如k=1, 0= -1, 0=1 ,1=0 就是Euler方法; k=1, 0= -1, 0= 1=1/2 就是梯形方法.,其中 任意给出,
2、s = 0,1,2,迭代到满足给定的精度要求为止.,对于隐式情形的公式( ), 由于f (x , y)一般是非线性函数, 很难求得 yn+k 的显式表达式, 故常用迭代法求解:,可以证明, 当f (x , y)关于 y 满足Lipschitz条件(或 )时,若 , 则为隐式方法; 若 , 则为显式方法. 上式称为k 步线性法.,只要 , 迭代关系式就是收敛的.,Def : 对于线性多步法,定义 xn+k 处的局部截断误差为:,或在 yn+j = y ( xn+j ) ( j = 0,1,2, k -1)的假定下, 定义 xn+k 处的局部截断误差为Rn+k = y ( xn+k ) - yn+
3、k .,Rn+k 按 h 展开的首项称为主局部截断误差.,Def : 若 , 则称线性多步法为 p 阶方法.,若线性多步法为 p 阶方法, 则,设 yn+j = y ( xn+j ), j = 0,1,2, k -1, 则有,其中Cp+1称为主局部截断误差系数.,又,上两式相减, 得:,由微分中值定理,即主局部截断误差为:,其中 介于 y ( xn+k )与 yn+k 之间, 因此,故在 y ( xn+j ) = yn+j ( j = 0,1,2, k -1) 假设下, 若 (显式公式), 则:,若 (隐式公式), 且 , 即为 p 阶方法, 则当 y(x)充分可微时,即Rn+k 的首项与 y
4、 ( xn+k ) - yn+k 的首项相同, 因此两种方法局部截断误差的定义是相同的.,2. 线性多步法的构造,注2: 可以证明, 显式线性多步法的整体截断误差比局部截断误差低一阶., 基于数值积分的构造方法,注1: Euler方法的主局部截断误差为: 梯形法的主局部截断误差为:,将 两端从 yn -j 到 yn+k 积分, 得:,构造 p 次Lagrange插值多项式:,其中,对 f ( t , y (t ), 取等距节点 , 对应的函数值为,代入 式有,用 yn - i 代替 y( xn -i ), fn-i 表示 f ( xn-i , yn-i ), 则得线性多步法显式公式:,所以,若
5、令 t = xn + s h , 则,(1) 取 k=1, j=0可得Adams显式公式:,对k , j 和 p 选择不同的值, 则可得到不同类型的具体公式.,其中,其中,再用n+1代替 n , 得到:,最常用的是 p=3 的情形:,(2) 取 k=0, j=1可得Adams 隐式公式:,最常用的也是 p=3 的情形:,(3) 取 k=1, j=1可得Nystrm 显式公式:,于是线性多步法,(4) 取 k=0, j=2可得Milne 隐式公式(用n+1代替n ):,由于,的局部截断误差为:,对于Adams显式与隐式公式, 由于 在0,1( j=0, k=1, 显式公式) 上恒为正; 在0,1
6、( j = 0, k =1, 隐式公式) 上恒为负. 由微分中值定理可得:,当p=3时,其中 为某中间点, E表示Explicit, I表示Implicit .,由此可知, 在 f (x , y(x)具有p+1阶连续偏导数的条件下, p+1阶Adams显式方法和隐式方法的局部截断误差与 同阶, 即它们都是p+1阶方法. 特别地, 当p=3时, 它们都是4阶方法.,上式右边近似认为,Taylor展开法更具一般性, 下面构造4步方法, 即k=4 . 用Taylor展开法构造下述公式, 使其为4阶方法, 并求其局部截断误差的主项., 基于Taylor展开的构造方法,将函数 及 在点 xn 处按 h的
7、幂展开, 则得:,(1),比较 (1), (2) 两式中h 的同次幂, 可得:,另一方面, 有,(2),如果要求前6个方程成立, 得到的公式的局部截断误差为O( h6 ). 也可要求前面的几个方程成立.,(1) Simpson公式:,要求前4个方程成立, 令 可得:,公式为:,即:,此时, 第五个方程恰好也成立. (2) - (1) 得:,得隐式的Adams公式:,(2) 隐式的Adams公式:,要求前5个方程成立, 令 , 可得:,将 (2) - (1) 得:,(3) 显式的Adams公式:,要求前5个方程成立, 令 , 可得:,得显式的Adams公式:,将 (2) - (1) 得:,(4)
8、 Milne 方法(四步四阶):,要求前4个方程成立, 令 , 可得:,此时, 第五个方程恰好成立.,得Milne 公式:,将 (2) - (1) 得:,(5) Hamming方法:,要求前5个方程成立, 令 , 可得:,即:,得Hamming 公式:,将 (2) - (1) 得:,Remark 1: 用Taylor展开法比用数值积分法推导线性多步法更加灵活, 推导出的方法也比积分法更多. 用数值积分法可以得到的方法, 用Taylor展开方法都可以得到; 但用Taylor展开法得到的方法, 数值积分法却不一定能够得到.,例: 求微分方程初值问题,在0,1上的解, 取步长 h = 0.1 .,Remark 2: 应用线性多步法求解初值问题时, 开始几点处的函数值要用别的方法先计算出来. 一般选用与多步法同阶或更高阶的单步法. 如Runge-Kutta方法、Taylor展开方法等.,Remark 3: 对线性隐式多步法, 除开始几点处的函数值需要单独计算外,
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 农村电子商务创业扶持政策手册
- 迪庆2025年云南迪庆州事业单位招聘工作人员291人笔试历年参考题库附带答案详解
- 茂名广东茂名市自然资源局电白分局招聘劳务派遣工作人员10人笔试历年参考题库附带答案详解
- 盐城2025年江苏盐城射阳县教育局下属事业单位赴高校招聘教师22人笔试历年参考题库附带答案详解
- 甘肃2025年甘肃交通职业技术学院招聘博士研究生笔试历年参考题库附带答案详解
- 贵州省部分高中2024-2025学年高一上学期1月期末考试化学试题(解析版)
- 2025年贵州黔南布依族苗族自治州高三二模政治试卷(2月)【含答案】
- 2022年中级经济师考试《建筑与房地产经济专业知识与实务》真题及答案
- 2022年高级经济实务《人力资源管理》真题及解析
- 阿奇霉素与妇科千金胶囊联合用药应用于慢性附件炎作用分析
- 腹腔引流管护理技术操作规程及评分标准
- 少年管乐团管理制度
- 中学生睡眠质量研究性学习报告
- 酒店水单账单范本
- 《思想道德与法治》第三章
- 空压机(储气罐)日常安全检查表
- 桥梁预应力结构张拉压浆智能化施工成套技术
- 谐波减速器仿真优化
- 11 我是一只小虫子(第二课时一等奖创新教案)
- 压力容器年度检查表
- 2023-2024学年湖北省洪湖市初中语文九年级下册期末通关试题
评论
0/150
提交评论