三阶线性常微分方程Sinc方程组的结构预处理方法.pdfVIP

2013年8月 计算数学 第35卷第3期 Aug．，2013 MATHEMATICANUMERICASINICA V01．35，No．3 三阶线性常微分方程Sine方程组 的结构预处理方法术 任志茹 (中国科学院数学与系统科学研究院，计算数学与科学工程计算研究所 科学与工程计算国家重点实验室，北京100190) 捅 要 三阶线性常微分方程在天文学和流体力学等学科的研究中有着广泛的应用．本文介绍求解三 阶线性常微分方程由Sine方法离散所得到的线性方程组的结构预处理方法，首先，我们利用Sine 方法对三阶线性常微分方程进行离散，证明了离散解以指数阶收敛到原问题的精确解．针对离散 后线性方程组的系数矩阵的特殊结构，提出了结构化的带状预处理子，并证明了预处理矩阵的特 征值位于复平面上的一个矩形区域之内．然后，我们引入新的变量将三阶线性常微分方程等价地 转化为由两个二阶线性常微分方程构成的常微分方程组，并利用Sine方法对降阶后的常微分方 程组进行离散．离散后线性方程组的系数矩阵是分块2×2的，且每一块都是Toeplitz矩阵与对 角矩阵的组合．为了利用Krylov子空间方法有效地求解离散后的线性方程组，我们给出了块对角 预处理子，并分析了预处理矩阵的性质．最后，我们对降阶后二阶线性常微分方程组进行了一些比 较研究．数值结果证实了Sine方法能够有效地求解三阶线性常微分方程． 关键词：三阶线性常微分方程；Sine方法；收敛性分析；预处理子；特征值估计 MR(2000)主题分类：65F08，65F10，65F15，65L10，65T10 1．引 言 考虑三阶线性常微分方程边值问题 (1．1) I y(a)=0，y(b)=0，Y(a)=0， 文学、流体力学等学科的研究中有着广泛的应用．比如，流体力学中用来模拟水面上孤波传播 问题[13-151等通过线性化以后都需要求解三阶线性常微分方程(1．1)．这类常微分方程的一个 显著特点是其最高项是三次的，从而使得离散后的线性方程组的系数矩阵是非对称部分占优 的，且是高度病态的，因此在求解此线性方程组时会有相当大的困难． 目前，关于采用Sine方法离散求解三阶线性常微分方程的数值方法的研究并不多见[18,27】 由于Sine方法具有较高的精度和能够处理奇异问题，它在求解常微分方程时是非常有效的． 中，sine函数的定义如下： sinc@)：—sinF(rt)，一oot。。． ’2013年1月31日收到 万方数据 计算数学 2013缸 Sinc方法是基于sinc函数的一类数值离散方法．我们用 蛳m)：=喘缘铲，一。。t。。，触， (1．2) 的第J个Sinc基函数． 题(1．1)的离散解以指数阶收敛到其精确解[5J-由于离散后线性方程组的系数矩阵A是一 是，对这种系数矩阵为Toeplitz矩阵和对角矩阵组合的线性方程组并没有快速直接法．究其原 矩阵的快速直接法并不适用于这类线性方程组． logn)的运算量来实 然而，对于任意一个佗维向量q，矩阵与向量的乘积Aq可以用O(n 空间方法的收敛速度，需要为矩阵A构造一个有效的预处理矩阵．本文中，我们针对离散后 所得到线性方程组的系数矩阵的特殊结构，构造出了结构化的预处理子【7_9，19121】，并证明了 预处理之后的矩阵的特征值是一致有界的，且位于复平面上一个与矩阵维数无关的矩形区域 内【5】． 由于采用Sinc方法