原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
那么，
2 2 2
E 2 ≤ E F ≤εn( T1 F U1 F + T2 F U2 F ) ≤εn y 2 ( 1+ x 2 + x 2 ) ≤εn y 2 (1+2 x 2 )
F ≤ε n T −1 。
2 2
考虑所有的误差边界，经过计算，可得
~
−1 −1 ~ −1
T −T 2 ≤n(2ε +nε) y 2 (1+2 x 2 ) +ε n T 2 。
由Ty e1 和 Tx f ，可得 y 2 ≤ T −1 2 和 x 2 ≤ T −1 2 f 2 ，因此，相对误差为
~
−1 −1
T −T
2 ~ −1
−1 ≤n(2ε +nε)(1+2 T 2 f 2 ) +ε n 。
T
2
因为T 是良态的，所以 T −1 是有界的。而 f 也是有界的。
2 2
因此，公式T −1 T U +T U 是向前稳定的。
1 1 2 2
2.4.3 算法分析
现在我们来考虑一下算法，对于公式T −1 T U +T U ，关键就是要解线性方
1 1 2 2
程组Tx f 和Ty e 。如果用 Zohar 算法求解 Toeplitz 线性方程组Tx f 和Ty e ，
1 1
则要8n 2 −4n −4 次乘除法运算，8n 2 −14n +6 次加减法运算；如果用快速傅里叶
（FFT ）方法，则运算复杂性为O(n log n) 。这些方法将在第四章里做详细的阐述。
因此，如果用快速傅里叶方法，得到分解式T −1 T U +T U 只需要O(n log n) 的运
1 1 2 2
算量。
15
第三章 特殊 Toeplitz 矩阵的逆矩阵的研究
特殊的 Toeplitz 矩阵有广泛的应用前景，特别在生物学、物理学、数学、社会
科学中的许多问题都和 Toeplitz 矩阵的理论有着密切的关系。在数理统计、石油、
地震物探及其它应用科学中，尤其是图象、数字信号处理中，常会遇到循环矩阵
这类特殊矩阵。带状 Toeplitz 方程组广泛应用于科学计算和工程计算中，尤其是在
双曲偏微分方程的数值求解中，随着科学技术的不断发展，问题的规模越来越大，
适应于大规模计算的并行算法也被陆续提出来。因此，对于这类特殊的 T
文档评论（0）