SineGordon方程的次配点法.pdfVIP

第31卷 第2期 工 程 数 学 学 报 Vo1．31N0．2 2014~z04f] CHINESE JOURNALOFENGINEERING MATHEMATICS Apr．2014 doi：10．3969／j．issn．1005—3085．2014．02．011 文章编号：1005—3085(2014)02—0254—13 Sine—Gordon方程的三次配点法冰 李 宏，黄春霞，何斯 日古楞，赵 猛，刘 洋 (内蒙古大学数学科学学院，呼和浩特 010021) 摘 要：Sine—Gordon方程在非线性光学和生物物理等多个物理 问题 中有着广泛 的应用 ．本 文研究一维Sine—Gordon方程的三次配点法 ，利用复合高斯求积公式近似 内积 的一 种离散化的HiGalerkin方法建立半离散和全离散格式 ．采用先验估计方法推理 了 ，日 和 日0模最优估计结果．通过Matlab软件编程计算，获得 了数值解和真实 解的对比结果及误差估计数据． 关键词：Sine-Gordon方程；三次配点法；高斯求积公式；误差估计 分类号：AMS(2000)65M05；65M10 中图分类号：O241 文献标识码：A 1 引言 考虑具有边界条件u 0， Jzcx~的一维Sine-Gordon方程 乱tt+ 钍t=，“一sin(u)， (1) 其中 0是阻尼系数．Sine—Gordon方程是应用科学中非常重要的非线性方程之一，在 凝聚态物理，非线性光学和生物物理等领域 中有广泛运用，因此得到学术界的广泛关注 ． 近年来出现了许多求解Sine—Gordon方程的数值方法，如有限元方法[1】’有限差分法2[]， 无网格法[3】，混合元方法 【J等． 在偏微分方程 的数值求解 问题 中，寻找一种较高效便捷 的数值方法至关重要，其 中三次配点法就是一种数值计算方便有效的高阶数值方法 ．该方法可 以看成离散的H Galerkin方法，也可 以被认为是 以三次样条函数为试验函数空间、以分段线性函数为 试探 函数空间，用复合两点高斯求积公式近似 内积 的Petrov—Galerkin方法．该方法与 正交三次样条配点法相 比较有如下一些优点：格式的收敛精度高；该方法所需未知数 仅为正交的三次样条配点法需要的未知数一半，从而降低 了计算成本 ．Sloan于 1988年 在文献[5】中提出了三次配点法，并在文献f61中研究二阶两点边值 问题的三次配点法； 后来，Jones和Pani等人将该方法推广到二阶半线性两点边值 问题 7【】一维抛物型方 程的初边值问题8【】Stefan问题l。1和Burgers方程1【o]等．文献 1『l1利用此方法数值求解 收稿 日期：2012—06—26． 作者简介：李宏 (1973年12月生)，女，博士，教授．研究方向：微分方程数值解． 基金项目：国家自然科学基金 11361035；内蒙古自治区高等学校科学研究项目(NJ10006 NJZZ12011)：内蒙古自然科学基金 (2012MS0106；2012MS0108)：内蒙古大学高层次人才项目(30105—125132； 125119)． 第2期 李宏，等：Sine-Gordon方程的三次配点法 了RLW—Burgers方程 ，并分析了半离散和全离散格式的最优误差估计 ．在此基础上， 本文利用该方法研究了具有物理边界条件 +0，X ~oc的一类Sine—Gordon方程， 证明半离散格式的 ()， (日)，L (。)， o。(日)和 o。(日。)模以及全离散格式 的 ，日 和 日。模最优误差估计． 2 半离散格式及最优误差估计 在数值计算中，我们不妨设X∈I三(a，b)．考虑具有边界条件 u(a，t)=0，u(b，t)=0，t0， 和初始条件 札(，0)=札0(z)，ut(，0)：9(z)，