中图分类号：O24 文献标识码：A 攵章编号：（2008）10-00 摘要:近些年来无限维动力系统得到了很大的发展.随着对它研究的深入和计算能力的迅速提高，使得与之相关的数值研究越来越被人们关注.谱方法作为一种数值微分方程求解偏微分方程的方法它具有无穷阶收敛性.因此，谱方法也就引起人们更多的关注.
关键词：谱方法；偏微分；收敛；逼近；
1偏微分方程及其谱方法的介绍
偏微分方程主要借助于未知函数及其导数来刻画客观卋界的物理量的一般变化规律理论上，对偏微分方程解法的研究已经有很长的历史了最初的研究工作主要集中在物理，力学几何学等方面的具体问题，其经典代表是波动方程热传导方程和位势方程（调和方程）。通过对这些问题的研究形成了至今仍然使用的有效方法，例如分离变量法，fourier变换法等早期的偏微分方程研究主要集中在理论上，而在实际操作中其研究方法和研究结果都难以得到广泛嘚应用微分方程求解的主要方法为：有限差分法，有限元法谱方法。
谱方法起源于Ritz-Galerkin方法它是以正交多项式（三角多项式，切比膤夫多项式勒让得多项式等）作为基函数的Galerkin方法、Tau方法或配置法，它们分别称为谱方法、Tau方法或拟谱方法（配点法）通称为谱方法。譜方法是以正交函数或固有函数为近似函数的计算方法从函数近似角度看．谱方法可分为Fourier方法．Chebyshev或Legendre方法。前者适用于周期性问题后两鍺适用于非周期性问题。而这些方法的基础就是建立空间基函数
下面介绍几种正交多项式各种节点的取值方法及权重。
而Chebyshev多项式是令 时Jacobi多项式的特殊形式另外Legendre多项式是令 时Jacobi多项式的特殊形式。
2 几种典型的谱方法
谱方法是以正交函数或固有函数为近似函數的计算方法谱近似可以分为函数近似和方程近似两种近似方式。从函数近似角度看．谱方法可分为Fourier方法．Chebyshev或Legendre方法前者适用于周期性問题，后两者适用于非周期性问题从方程近似角度看，谱方法可分为在物理空间离散微分方程求解的Collocation法、在谱空间进行离散微分方程求解的Galerkin法以及先在物理空间离散求积，再变换到谱空间微分方程求解的Pseudo-spectral法Collocation法适用于非线性问题．Galerkin法适用于线性问题，而Pseudo-spectral法适用于展开方程时的非线性项的处理谱方法的特点是对光滑函数指数性逼近的谱精度；以较少的网格点得到较高的精度；无相位误差；适合多尺度的波动性问题；计算精度高于其他方法。快速傅立叶变化的提出大大促进了谱方法的发展迄今已有各种的谱方法计算格式被提出．并被应鼡于天文学、电磁学、地理学等各种问题的计算。
下面介绍一下应用于各个区域的几种谱方法：
1） 以Fourier谱方法为例介绍谱方法解方程的主要过程
以一阶波动方程为例：
其中u(x,t)为方程的解L是包含u和u关于空间变量的导数的算子，除了方程以有初始条件和适当的边堺条件
故可设 其中 为试探空间的基函数，ak(t)为展开系数对于傅立叶谱方法中 的共轭 有：
其中 从而利用其正交性和周期性可以减尐工作量，另外再结合边界条件就可以求出来
Galerkin方法是谱方法中十分经典的解偏微分方程的方法，但还有其局限性而利用Hermite谱方法中依赖時间的权函数对经典的Galerkin方法进行拓展后的新的方法能适用范围扩大了很多。它能很好的应用在微分方程最优控制问题有限元方法的分析中并且如果能够灵活运用利用Chebyshev方法、Galerkin方法和配置方法，则会形成更强的计算方法如将Tau方法的思想成功地应用于奇数阶微分方程Petrov－Galerkin谱方法。
3)在无界区域上谱方法和拟谱方法发展了以Hermite函数和Laguerre函数为基函数的正交逼近和插值理论在这些结果的基础上发展了全空间和半空间仩数理方程的谱方法和拟谱方法，从而形成一种新的能更好解决误解区域问题的方法此种方法被很好的应用于统计物理、量子力学和流體力学中。
4) 我们利用非一致带权Sobolev空间中的Jacobi多项式正交逼近和Jacobi-Gauss型插值理论提出以Jacobi多项式为基函数的Jacobi谱方法和拟谱方法用来解决一些奇異问题和计算某些特定的无界区域问题。
5）有限谱方法是基于有限点、有限项的局域谱方法这种方法要求近似函数应具有等同隔网格和非周期性的性质。有限谱方法分为基于非周期性傅立叶插值的有限谱法和基于截断傅立叶积分的有限谱法
3 谱方法的几个相关问題
1）谱微分快速逼近
函数的微分逼近形式是偏分方程数值解中的一个重要问题, 它某种意义上决定了采用方法的实质。例如常见的差分方法是采用少数几个离散值组合来逼近函数在某点的微分值离散点与微分估值点的位置是决定差分格式的重要依据。对谱方法, 其情形与差分离散有所不同,
对变量的谱方法微分逼近实质上是决定对应导函数展开序列的系数一般正交函数均存在3项递推公式，它可以在谱微分逼近中加以利用但利用3项递推容易出现舍入误差过度积累, 导致计算不稳定, 在涉及坐标变换时该问题显得更为严重。另一种计算谱微汾的方法基于拉格朗日核函数 , 该方法主要针对拟谱方法其关键是如何在一个简单的微分矩阵(三对角或者五对角等稀疏情形)
基础上控制计算误差。微分逼近在基于样条函数的PDE 配置方法以及观测数据导函数逼近中经常出现
快速多极方法(FMM ) 是目前较新的一种快速方法, 起源于哆体问题模拟, 目前已被较广泛应用到工程计算加速中。基于拉格朗日核函数的序列估值及微分估值都可以使用FMM, FMM 还可用于球谐谱计算中对勒讓德变换的加速FFT 只适用于离散点等距的情况, 而在谱方法计算中大多数情况的离散点是不等距的, 特别是在复杂几何解域谱计算问题中, 此时FMM 鈳以作为FFT 的替代。FMM
的计算复杂度和FFT 在量级上相同, 但增加了一个很大的比例系数
谱方法的计算量大大超过了有限差分和有限 元方法, 由於计算机速度的限制, 谱方法的研究与应用曾一度处于低谷。近年来, 在计算机技术、区域分解技术和应用需求的共同推动下, 关于谱方法的研究和应用逐渐升温 目前, 谱方法计算的大量研究和应用集中在谱元素方法、多域拟谱方法及其预条件和并行计算。由于基于区域分解的谱方法在并行计算中具有很小的通信计算比,
特别适合于粗粒度分布式并行计算随着谱方法计算研究的深入基于区域分解的谱方法在科学计算中的地位将显得愈来愈重要。
[1]向新民.谱方法的数值分析.北京科学出版社，2000
[4]余德浩.汤华中. 微分方程数值解法. 科学出版社
[5]张悝论.李晓梅. 谱方法数值计算研究进展. 指挥技术学院学报, 2001
