单步法的收敛和稳定.ppt_第1页
单步法的收敛和稳定.ppt_第2页
单步法的收敛和稳定.ppt_第3页
单步法的收敛和稳定.ppt_第4页
单步法的收敛和稳定.ppt_第5页
已阅读5页，还剩6页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

8.3 单步法的收敛性和稳定性,8.3.2 单步法的稳定性,8.3.1 单步法的收敛性,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条件即存在常数 ，使,对一切 成立，则方法（8.1.8）收敛，且 。,因为（8.1.8）是p阶的，所以存在 ，当 时有 。 再用 的Lipschitz条件有,为了方便，记 ，即有 。由此可推得,利用关系式,（8.3.1）,称 为增量函数。对于收敛的方法，固定 ，有 从而 。对于（8.3.1），我们自然要考虑 是否成立。这就是相容性问题。,8.3.2 单步法的稳定性,对于一种收敛的相容的差分方程，由于计算过程中舍入误差总会存在，我们 需要讨论其数值稳定性。一个不稳定的差分方程会使计算解失真或计算失败。 为了讨论方便起见。将（8.1.1）中的 在解域内某一点 作 Taylor展开并局部线性化，即,现在讨论显式Euler法的稳定性。将显式Euler法用于试验方程（8.3.2），有 。当 有舍入误差时，其近似值为 ，从而有 。令 ，得到误差传播方程。,令 ，只要 ，则显式Euler方法的解和误差都不会恶性 发展，即 时，显式Euler方法是稳定的，即是条件稳定的。 对于梯形方法，应用于试验方程后，有,同理，有误差方程 ，其中 。 因此当 时，梯形方法是稳定的。,定义8.5 若（8.3.3）式中的 ，则称对应的单步法是绝对稳定的。在复 平面上， 满足 的区域，称为方法的绝对稳定区域，它与实轴的交 称为绝对稳定区间。,一些单步法的 表达式和它们的绝对稳定区间列于表8-4。从表中可见， 隐式方法比显式方法的绝对稳定性好。,解 本题 分别为-1，-2，-4。有表8-4可知，当 时，该方法 才稳定，计算结果列于表8-5,由表8-5可见，h=1和h=2时，计算结果确实稳定，h=4时，结果发散。此外，h为1的计算精度比h为2的计算精度高。因为h 越小，方法的截断误差越小。但若h过分
人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论