原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
非奇异， 。线性方程组 可
以转化为 。可以利用迭代法求
解线性方程组。
选定初始向量
;使用迭代法收敛条件：
（1）
（2） ，因为
; 如何转化为 ？
----- A=M+N
----- M x = N x + b, M分裂矩阵，选不同的M即
可构造出不同的
;M选对角矩阵D：
即称为雅可比迭代法。
;当M取对角矩阵D和下三角矩阵L之差：
称为高斯-塞德尔迭代法;Jacobi 和Gauss-Seidel的计算：
Jacobi：
Gauss-Seidel：
区别在于计算 时充分利用最新的迭代值，即
代替
;Jacobi 和Gauss-Seidel收敛性：
Ax=b：
（1）A为严格对角占优矩阵。都收敛
（2）A弱对角占优矩阵，且不可约矩阵。都收敛
（3）A对称正定时：
雅克比收敛： 充要条件A，2D-A正定；
Guass收敛： 成分条件，A正定。
;;;第七章 非线性方程与方程组的数值解法 ;;;;迭代法一：牛顿法
迭代函数（无重根）：
二次收敛。
有m重根：
知道重根数：
不知道重根数：
;;第8章 矩阵特征值问题计算 ;;;;;;;2.
3.将红色部分进行约化，直至A变换为上海森伯格矩阵的形式
;;第9章 常微分方程初值问题数值解法; 微分方程数值计算的一般思路：
算出精确解y(x)在区间?a,b?上的一系列离散节点
处的函数值
的近似值
相邻两个节点的间距 称为步长
微分方程数值解关心的问题：
（1）局部的截断误差和阶数；
（2）数值解Yn的误差估计和收敛性；
（3）递推公式的稳定性；;;;;;四：龙格库塔法
在前面的指导思想下，在一个步长内，如果选取多个点的斜率的一个组
合就可以得到更高的精度。龙格库塔法就是基于此提出的，如何组合这
些斜率使得局部截断误差最小，阶数最高。
以三阶为例：取三点处Xi ， ， 的斜
率做加权平均，计算格式为
利用二元泰勒展开式，得到消去 ， 的条件即可得到p，q， 应满足
的条件：给出一种形式;算法稳定性问题
数值法在计算时候伴随着各种误差（计算机舍入等）。这些
误差在以后的计算过程中是否会放大，还是可以忽略就是稳
定性问题：
比如：在某次计算中出现δ的扰动，在以后的节点
计算中偏差都不大于δ时称为稳定：
讨论稳定性以下面的形式为标准，都可以转化：
;以欧拉法为例：
模型方程：
欧拉法递推公式 ：
反复递推后：
为了使误差不扩大，应该有 ，所以用欧拉法稳定
的充要条件就是：
其他方法的稳定性判断类似。 ;算法收敛性问题
递推公式中用yi来近似替代y(xi)，这种近似替代是否合理，
还要进行判断即： ， 是否成立，也就
是整体的误差要趋于0的。
同样以欧拉法为例：
欧拉法局部截断误差：
则总体截断误差；
;根据李普希茨条件有
即：
反复递推后有：
初值误差为0所以：
当h-0,误差-0，欧拉法整体阶段误差为 ，算法精度为一阶
文档评论（0）