偏微分方程数值解法抛物型方程差分法2.pptVIP

下载本文档

13
0
约3.67千字
约 17页
2020-05-06 发布于天津
举报
版权申诉

偏微分方程数值解法抛物型方程差分法2.ppt

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/17 《偏微分方程数值解法》 7 —— 抛物型方程差分法 2 差分格式稳定性概念显、隐格式稳定性分析稳定性分析的矩阵方法 ? ? ? 2/17 x t 2 / h r ? ? 抛物型方程 k j k j x k j k j f u h a u u ? ? ? ? 2 2 2 1 ? ? k j k j k j k j k j f u u ra u ra u ? ? ? ? ? ? ? ? ? ) ( ) 2 1 ( 1 1 2 2 1 ) , ( 2 2 2 t x f x u a t u ? ? ? ? ? ? 简单显式差分格式 ? 在实际应用时 , 取逐层计算形式 . 当初始层数据有误差时 , 误差会逐层传播 , 影响以后各层的解 . 记的误差为 , 设无误差 , 则有 k j u k j ? k j f ) ( ) 2 1 ( 1 1 2 2 1 k j k j k j k j ra ra ? ? ? ? ? ? ? ? ? ? ? 3/17 ) ( 2 1 1 1 1 k j k j k j ? ? ? ? ? ? ? ? 2 / 1 2 ? ra 取 ) ( ) 2 1 ( 1 1 2 2 1 k j k j k j k j ra ra ? ? ? ? ? ? ? ? ? ? ? 设初始层上 , 仅有 , 其它点处无误差 ? ? ? 0 0 j ? 在各计算层上 , 误差传播得到控制 4/17 ) ( ) 2 1 ( 1 1 2 2 1 k j k j k j k j ra ra ? ? ? ? ? ? ? ? ? ? ? 1 2 ? ra 取 ? ) ( 1 1 1 k j k j k j k j ? ? ? ? ? ? ? ? ? ? ? 设初始层上 , 仅有 , 其它点处无误差 ? ? ? 0 0 j 在各计算层上 , 误差传播没有得到控制 5/17 无穷大范数定义双层差分格式记矩阵双层格式的矩阵形式 k k k k k f u B u A ? ? ? ? ) ( 1 ) ( 双层差分格式初值稳定概念 : k k k k B A ? ? ) ( 1 ) ( ? ? 任意解都满足 || || || || 0 k k M ? ? ? 其中 M 与无关 . k k 0 ? | | max || || 1 k j n j k u u ? ? ? ? k j n m k m k jm n m k m k jm f u u ? ? ? ? ? ? ? ? ? ? 1 ) ( 1 1 ) ( n n k jm k A ? ? ) ( ) ( ) ( ? n n k jm k B ? ? ) ( ) ( ) ( ? 6/17 简单显式差分格式 ) ( ) 2 1 ( 1 1 2 2 1 k j k j k j k j u u ra u ra u ? ? ? ? ? ? ? ) ( 0 j j x u ? ? 0 1 0 ? ? ? k n k u u ) , , 2 , 1 ( n j ? ? 稳定性分析 , 设 2 / 1 2 ? ra |) | | (| | | ) 2 1 ( | | 1 1 2 2 1 k j k j k j k j u u ra u ra u ? ? ? ? ? ? ? ? ? ? ? ? || || 2 || || ) 2 1 ( 2 2 k k u ra u ra ? ? ? || || | | 1 k k j u u ? ? ? ? ? || || || || 1 k C k u u 此时差分格式稳定 7/17 k k k k B A ? ? ) ( 1 ) ( ? ? 设齐次方程系数矩阵可逆 ? k k k k B A ? ? ) ( 1 ) ( 1 ] [ ? ? ? 记称之为过渡矩阵 ? ) ( 1 ) ( ) ( ] [ k k k B A H ? ? k k k H ? ? ) ( 1 ? ? k k H ? ? ? ? 1 常系数差分格式 ? H 的谱半径 : | ) ( | max ) ( 1 H H j n j ? ? ? ? ? ? ? 1 1 ) ( M H ? ? 定理 : 双层差分格式稳定的必要条件是 , 存在与无关的常数 M 1 , 使得 ? 8/17 定理若 H = A -1 B 为正规矩阵 , 即 HH * = H * H , 则条件 ? ? 1 1 ) ( M H ? ? 是双层差分格式按欧氏范数稳定的充分条件注：欧氏范数 ( 或离散 L 2 范数 ) 2 / 1 1 2 0 ] ) ( [ || || ? ? ? n j k j k u h u 9/17 简单显式差分格式矩阵形式 k k k f u C ra