版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
单步法的收敛性和稳定性第一页,共十一页,2022年,8月28日8.3.1单步法的收敛性
数值解法的基本思想就是要通过某种离散化方法,将微分方程转化为某种差分方程(例如,(8.1.8)式)来求解。这种转化是否合理,还要看差分方程的解,是否收敛到微分方程的准确解。
定义8.3
对于任意固定的,若对于初值问题(8.1.1)的显式单步法(8.1.8)产生的近似解,均有,则称该方法是收敛的。在定义中,是固定的点,当时有,n不是固定的。显然,若方法是收敛的,则在固定点处的整体截断误差趋于零。下面给出方法收敛的条件。
定理8.1
设初值问题(8.1.11)的单步法(8.1.8)是p阶的(),且函数满足对y的Lipschitz条件即存在常数,使第二页,共十一页,2022年,8月28日对一切成立,则方法(8.1.8)收敛,且。因为(8.1.8)是p阶的,所以存在,当时有。再用的Lipschitz条件有为了方便,记,即有。由此可推得证仍记,根据局部截断误差的定义将此式与(8.1.8)相减得第三页,共十一页,2022年,8月28日利用关系式可以得到现在取,有,于是有。定理得证。
容易证明,如果(8.1.1)的满足Lipschitz条件是,且初值是正确的,则显示Euler法、改进的Euler法和R-K方法是收敛的。由定理8.1说明,f关于y满足Lipschitz条件是使单步收敛的充分条件,而且,还说明一个方法的整体截断误差比局部截断误差低一阶。所以,常常通过求出局部截断误差去了解整体截断误差的大小。
单步法的显式形式(8.1.8)可写成
(8.3.1)第四页,共十一页,2022年,8月28日称为增量函数。对于收敛的方法,固定,有从而。对于(8.3.1),我们自然要考虑
是否成立。这就是相容性问题。则称方法(8.1.8)与初值问题(8.1.1)是相容的。
相容性说明数值计算的差分方程(8.3.1)趋于(8.1.1)中微分方程。我们本章讨论的数值方法都是与原初值问题相容的。定义8.4若方法(8.1.8)的增量函数满足第五页,共十一页,2022年,8月28日8.3.2单步法的稳定性
对于一种收敛的相容的差分方程,由于计算过程中舍入误差总会存在,我们需要讨论其数值稳定性。一个不稳定的差分方程会使计算解失真或计算失败。
为了讨论方便起见。将(8.1.1)中的在解域内某一点作Taylor展开并局部线性化,即令利用线性化的关系,可得。因此,我们通过如下的试验方程第六页,共十一页,2022年,8月28日(8.3.2)讨论数值方法的稳定性。当某一步有舍入误差时,若以后的计算中不会逐步扩大,则称这种稳定性为绝对稳定性。
现在讨论显式Euler法的稳定性。将显式Euler法用于试验方程(8.3.2),有
。当有舍入误差时,其近似值为,从而有
。令,得到误差传播方程。令,只要,则显式Euler方法的解和误差都不会恶性发展,即时,显式Euler方法是稳定的,即是条件稳定的。
对于梯形方法,应用于试验方程后,有第七页,共十一页,2022年,8月28日同理,有误差方程,其中。因此当时,梯形方法是稳定的。
一般地,在试验方程(8.3.2)中,我们只考虑的情形,而对的情形,我们认为微分方程是不稳定的。比如,将显式Euler方法用于(8.1.1)中的方程,有当时,有。对于每一种单步法应用于试验方程(8.3.2),可得
(8.3.3)然而,对于不同的单步法,有不同的表达式。第八页,共十一页,2022年,8月28日定义8.5
若(8.3.3)式中的,则称对应的单步法是绝对稳定的。在复平面上,满足的区域,称为方法的绝对稳定区域,它与实轴的交称为绝对稳定区间。
一些单步法的表达式和它们的绝对稳定区间列于表8-4。从表中可见,隐式方法比显式方法的绝对稳定性好。表8-4
方法绝对稳定区间Euler法改进的Euler法三阶R-K法四阶R-K法隐式Euler法梯形式第九页,共十一页,2022年,8月28日例8.4分别取h=1,2,4,用经典R-K方法计算其准确解为。
解本题分别为-1,-2,-4。有表8-4可知,当时,该方法才稳定,计算结果列于表8-5h=1的解h=2的解h=4的解准确解表8-553.63943.67305.47153.638997.63237.636716.82917.63221311.632111.632657.617
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 《2024年 信息技术与课程整合的新趋向-基于iPad+Apps的课程设计研究》范文
- 《 具有串并混联结构腿的四足机器人设计》范文
- DB21-T 3996.2-2024 农田主要杂草防控技术规程 第2部分:玉米田
- 出版行业数字化出版流程优化与管理
- 红星照耀中国优美词汇摘抄
- 农业行业智能农业大数据方案
- 苏教版六年级数学下册教学设计-【六 正比例和反比例】
- 【人教版】二年级数学下册导学案设计《第1课时 简单的推理(1)》
- 生产作业场所危害因素检测制度
- 电工基础常识概要
- GB/T 14100-1993燃气轮机验收试验
- 银行个人信贷业务档案管理办法
- 摄影测量与遥感-课件
- 第23课《范进中举》课件(41张PPT) 部编版语文九年级上册
- 重大隐患限期整改计划和防范措施(3篇)
- 配电线路工(题库)
- 项目一-发售车票课件
- 2023年四川省自然资源投资集团有限责任公司招聘笔试题库及答案解析
- 机械制造工艺学-课后答案
- 结业证书word模板
- 脊柱侧弯科普(脊柱侧弯那点事)课件
评论
0/150
提交评论