Sine-Gordon方程的三次配点法.pdfVIP

第31卷 第2期 工 程 数 学 学 报 Vol. 31 No. 2 2014年04月 CHINESE JOURNAL OF ENGINEERING MATHEMATICS Apr. 2014 doi: 10.3969/j.issn.1005-3085.2014.02.011 文章编号: 1005-3085(2014)02-0254-13 方程的三次配点法∗ 李 宏, 黄春霞, 何斯日古楞, 赵 猛, 刘 洋 (内蒙古大学数学科学学院，呼和浩特 010021) 摘 要: Sine-Gordon 方程在非线性光学和生物物理等多个物理问题中有着广泛的应用．本 文研究一维Sine-Gordon 方程的三次配点法，利用复合高斯求积公式近似内积的一 种离散化的H –Galerkin 方法建立半离散和全离散格式．采用先验估计方法推理 了L , H 和H 模最优估计结果．通过Matlab 软件编程计算，获得了数值解和真实 解的对比结果及误差估计数据． 关键词: Sine–Gordon 方程；三次配点法；高斯求积公式；误差估计 分类号: (2000) 65M05; 65M10 中图分类号: O241 文献标识码: A 引引引言言言 考虑具有边界条件u → 0, x → ±∞ 的一维Sine-Gordon 方程 u + βu = u − sin(u), (1) 其中β ≥ 0 是阻尼系数．Sine-Gordon 方程是应用科学中非常重要的非线性方程之一，在 凝聚态物理，非线性光学和生物物理等领域中有广泛运用，因此得到学术界的广泛关注． 近年来出现了许多求解Sine-Gordon 方程的数值方法，如有限元方法 ，有限差分法 ， 无网格法 ，混合元方法 等． 在偏微分方程的数值求解问题中，寻找一种较高效便捷的数值方法至关重要，其 中三次配点法就是一种数值计算方便有效的高阶数值方法．该方法可以看成离散的H – Galerkin 方法，也可以被认为是以三次样条函数为试验函数空间、以分段线性函数为 试探函数空间，用复合两点高斯求积公式近似内积的Petrov-Galerkin 方法．该方法与 正交三次样条配点法相比较有如下一些优点：格式的收敛精度高；该方法所需未知数 仅为正交的三次样条配点法需要的未知数一半，从而降低了计算成本．Sloan 于1988 年 在文献[5] 中提出了三次配点法，并在文献[6] 中研究二阶两点边值问题的三次配点法； 后来，Jones 和Pani 等人将该方法推广到二阶半线性两点边值问题 ，一维抛物型方 程的初边值问题 ，Stefan 问题 和Burgers 方程 等．文献[11] 利用此方法数值求解 收稿日期: 2012-06-26. 作者简介: 李宏(1973年12月生)，女，博士，教授. 研究方向：微分方程数值解. 基金项目: 国家自然科学基金；内蒙古自治区高等学校科学研究项目(NJ10006; NJZZ12011) ；内蒙古自然科学基金(2012MS0106; 2012MS0108)；内蒙古大学高层次人才项目(30105-125132; 125119). 第2期