一类四阶微积分方程的Legendre-Galerkin谱逼近.pdf

2013年 5月 计 算 数 学 第 35卷第 2期 May,2013 MATHEMATICA NUMERICA SINICA Vo1．35，No．2 一 类四阶微积分方程的 Legendre—Galerkin谱逼近木1) 任全伟 庄清渠 (华侨大学数学科学学院，福建泉州362021) 摘 要 针对研究吊桥模型而建立的四阶微积分方程，提出Legendre谱逼近法进行求解．构造迭代 算法来求解得到的线性系统，证明了迭代格式的收敛性，对问题进行了误差分析．数值算例验证了 迭代的收敛性和方法的高精度． 关键词：四阶微积分方程；Legendre谱逼近；迭代算法；误差分析 MR (2000)主题分类：65M70，45305 1．引 言 科学和工程中的很多问题，最终可归结为各种各样的微分方程，为这些微分方程寻求快速 有效的解法，已经成为当今计算数学工作者的主要工作．谱方法的指数阶收敛性和应用的灵 活性越来越受到计算工作者的青睐．在过去的几十年中，随着计算机的发展和快速 Fourier变 换的出现，谱方法的应用越来越广泛．谱方法在单区域上的研究已趋于成熟，如 1【,2】等对谱 方法的几种基本形式做了详细的介绍．谱方法在多区域上的发展也趋于完善，多区域谱方法 的研究可参考 3【,4]等．谱方法是解决很多微分方程的重要工具，近年来也被广泛运用于解决 微积分方程，如 Volterra和 Fredholm类型积分方程 的求解 _5，6J等． 在研究长度为L的吊桥模型中，当负载量为p(x)时可得如下的方程 )_ yr!+ L )如 0 L， (1_1) I(0)= ()=Y”(0)=Y”(L)=0． 这里 和 都是正常数．假设p(x) 0(对于p(x) 0我们将得到类似的结论)，且为区间 [0，L]上的连续函数．方程的推导过程可以参考 [7]． 该四阶微积分方程已经得到部分工作者的关注，如文献 8『1给出了求解该方程的有限元方 法及相应的数值结果，9【】介绍了混合有限元方法的误差估计．文献 1【0]则介绍了该四阶微积 分方程的有限元方法和迭代算法．关于更多四阶微积分方程的研究和应用，可参考文献 [11，12] 等．本文将利用 Legendre—Galerkin谱方法对该微积分方程进行求解． 本文余下内容安排如下：在第 2节，我们给出问题的Legendre—Galerkin谱逼近形式及迭 代格式；迭代格式的收敛性证明将在第 3节给出；第 4节对 Legendre．Galerkin谱逼近形式进 行误差分析；最后在第 5节通过数值例子来验证理论分析的正确性． 2012年 4月 24日收到． )基金项目：国家 自然科学基金项目(No、福建省自然科学基金项目(No．2o11Jo5oo5)、中央高校 基本科研业务费专项资金和华侨大学侨办科研基金资助项 目(10QZR21)． 126 计 算 数 学 2013位 2．Legendre-Galerkin谱逼近 我们将在这一节中给出问题 (1．1)的Legendre—Galerkin谱逼近形式及相应的离散线性系 统．首先记A=(0，)．L (A)表示平方可积函数空间，该空间的内积和范数定义如下： u(，”)A=／u()()d，I钆『IA『：(u，)． 对任意的非负整数r，H (A)是通常意义下的Sobolev空间，其上的半范数和全范数定义为： --rA，=”《--A，--”A=(塞一；A，)． 其