求解抛物型方程的一种高精度紧致差分格式.PDF

第３８卷第２期 湖北大学学报（自然科学版） Ｖｏｌ．３８　 Ｎｏ．２　 ２０１６年３月 Ｊｏｕｒｎａｌ ｏｆ Ｈｕｂｅｉ Ｕｎｉｖｅｒｓｉｔｙ（Ｎａｔｕｒａｌ Ｓｃｉｅｎｃｅ） Ｍａｒ．，２０１６　 文章编号：１０００ ２３７５（２０１６）０１ ０１６０ ０８ 求解抛物型方程的一种高精度紧致差分格式 杨晓佳，魏剑英 （宁夏大学数学与计算机学院，宁夏 银川 ７５００２１） 摘要：利用四阶Ｐａｄé逼近公式和扩展的１／ ３⁃Ｓｉｍｐｓｏｎ公式，构造一种求解一维抛物型方程的高精度紧致隐式差分格 ４ ４ ＋ 式，其截断误差为Ｏ（τ ｈ ）． 然后通过理论分析证明此格式是无条件稳定的，并通过数值实验验证本文中格式的精确性 和可靠性． 关键词：抛物型方程；高精度紧致格式；无条件稳定；扩展的１／ ３⁃Ｓｉｍｐｓｏｎ公式；Ｐａｄé逼近 中图分类号：Ｏ２４１．８２　 　 文献标志码：Ａ　 　 ＤＯＩ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００⁃２３７５．２０１６．０２．０１３ Ａ ｈｉｇｈ⁃ｏｒｄｅｒ ｃｏｍｐａｃｔ ｄｉｆｆｅｒｅｎｃｅ ｓｃｈｅｍｅ ｆｏｒ ｔｈｅ ｐａｒａｂｏｌｉｃ ｅｑｕａｔｉｏｎ ＹＡＮＧ Ｘｉａｏｊｉａ，ＷＥＩ Ｊｉａｎｙｉｎｇ （Ｃｏｌｌｅｇｅ ｏｆ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｃｏｍｐｕｔｅｒ Ｓｃｉｅｎｃｅ，Ｎｉｎｇｘｉａ Ｕｎｉｖｅｒｓｉｔｙ，Ｙｉｎｃｈｕａｎ７５００２１，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ａ ｈｉｇｈ ａｃｃｕｒａｔｅ ｉｍｐｌｉｃｉｔ ｃｏｍｐａｃｔ ｄｉｆｆｅｒｅｎｃｅ ｓｃｈｅｍｅ ｉｓ ｃｏｎｓｔｒｕｃｔｅｄ ｆｏｒ ｓｏｌｖｉｎｇ ｔｈｅ ｏｎｅ ｄｉｍｅｎｔｉｏｎａｌ ｐａｒａｂｏｌｉｃ ｅｑｕａｔｉｏｎ ｂｙ ｔｈｅ ｆｏｕｒｔｈ⁃ｏｒｄｅｒ Ｐａｄé ｆｏｒｍｕｌａ ｃｏｍｂｉｎｅｄ ｗｉｔｈ ｔｉｍｅ ｅｘｔｅｎｓｉｏｎ ｏｆ ｔｈｅ １／ ３⁃ ４ ４ ＋ Ｓｉｍｐｓｏｎ ｆｏｒｍｕｌａｓ． Ｔｈｅ ｔｒｕｎｃａｔｉｏｎ ｅｒｒｏｒ ｏｆ ｔｈｅ ｓｃｈｅｍｅ ｉｓ Ｏ（τ ｈ ）．Ｔｈｅ ｓｔａｂｉｌｉｔｙ ｏｆ ｔｈｅ ｓｃｈｅｍｅ ｉｓ ａｎａｌｙｚｅｄ ｂｙ ｔｈｅｏｒｉｔｉｃａｌ ｍｅｔｈｏｄ． Ｉｔｉｓｕｎｃｏｎｄｉｔｉｏｎａｌｌｙ ｓｔａｂｌｅ．Ｎｕｍｅｒｉｃａｌｅｘｐｅｒｉｍｅｎｔｓａｒｅｃｏｎｄｕｃｔｅｄｔｏｖｅｒｉｆｙｔｈｅａｃｃｕｒａｃｙ ａｎｄ ｒｅｌｉａｂｉｌｉｔｙ ｏｆ ｔｈｅ ｐｒｅｓｅｎｔ ｓｃｈｅｍｅ． Ｋｅｙ ｗｏｒｄｓ：ｐａｒａｂｏｌｉｃ ｅｑｕａｔｉｏｎ；ｈｉｇｈ⁃ｏｒｄｅｒ ｃｏｍｐａｃｔ ｄｉｆｆｅｒｅｎｃｅ ｓｃｈｅｍｅ；ｕｎｃｏｎｄｉｔｉｏｎａｌｌｙ ｓｔａｂｌｅ；ｅｘｔｅｎｓｉｏｎ ｏｆ ｔｈｅ １／ ３⁃Ｓｉｍｐｓｏｎ ｆｏｒｍｕｌａｓ；Ｐａｄé ａｐｐｒｏｘｉｍａｔｉｏｎ ０　 引言 抛物型方程的数值解法一直是科研工作者们关注的热点问题之一，该类方程可描述众多物理现象， 如大气、海洋、江河及地下水污染物的扩散过程，因此对此类方程模型进行数值研究具有重要意义． 数值 求解该类方程的关键问题是如何在计算工作量尽可能小的前提下，使计算精度尽可能得到提高． 目前， ［１⁃９］ 已经发展了很多求解该方程的数值方法 ，其中Ｃｒａｎｋ⁃Ｎｉｃｏｌｓｏｎ方法是具有二阶精度的绝对稳定的经