Chebyshev、Fourier伪谱方法二阶、四阶有限差分法求解波动方程精度比较.pdfVIP

http :// Chebyshev、Fourier伪谱方法与二阶、四阶有限差分 法求解波动方程的精度比较1 李本文，耿佃桥，赫冀成 东北大学材料电磁过程教育部重点实验室(110004) E-mail: heatli@ 摘 要：本文针对具有精确解的波动方程u tt a 2 u xx ，分别采用Chebyshev、Fourier伪谱方 法，和二阶、四阶有限差分法求解，并对各种解法得到的结果进行了分析。结果表明：在相 同的较少节点数时，伪谱方法的结果精度明显高于差分法的结果；当节点数增加到一定数时， 4种方法的结果精度基本一致。总体而言，Chebyshev、Fourier伪谱方法的计算效率高于差 分法。 关键词：波动方程；伪谱方法；有限差分法；数值计算 1．引言 近二十年来，有关谱方法求解微分或偏微分方程的研究已经引起国内外学者的广泛关 ， ， 注[1-6] ，主要原因是谱方法能够提供很高的计算精度[1 3 5] 。大多数数值方法的计算精度随节 点数的增加线形提高，而谱方法的精度随节点数的增加呈指数趋势提高。本文运用 Chebyshev 、Fourier 伪谱方法，和二阶、四阶有限差分法求解具有精确解的波动方程 utt a 2 uxx ，并分析各种解法的结果以比较谱方法与有限差分法的优劣。结果表明：在相同 的较少节点数时，伪谱方法精度明显高于差分法；当节点数增加到一定数时，4 种方法的结 果基本一致。在计算同样的问题时，使用谱方法可以在计算量较少的情况下得到精度较高的 结果。 2. 波动方程及其精确解 [7] 采用文献 给出的波动方程及相应的边界和初始条件 ⎧ 2 u a u x l t , 0 , 0 ≤ ≤ ≥ ⎪tt xx ( ) ( ) ( ) (1) u x x l u x ⎨ ,0 sinπ / , t ,0 0 ⎪( ) ( ) u t u l t ⎩ 0, , 0 为了方便计算，这里取a 2 1 。 将初始和边界条件代入达伦贝尔公式[8]，分析可得如下精确解 ( ) ( ) ( ) (2) u x,t sin Nπx / l ⋅cos Nπt / l 1 本课题得到教育部高等学校博士学科点专项科研基金项目（2001145023）以及国家自然基金和上海宝钢 集团公司联合资助项目）资助。 http :// 3. 波动方程的伪谱方法求解 谱方法中被求解的变量用特定函数的有限项级数展开表示。经常采用的特定函数有 Fourier多项式、Chebyshev多项式、Bessel多项式等[1] 。只要级数中的每一项的系数被确定， 该被求解变量的近似形式也就确定了。级数的项越多，其精度也越高。与其它数值方法不同 的是[1-6]，在谱方法中要建立的代数方程可以是关于这些系数的代数方程，也可以是关于节 点上被求变量的代数方程。本文将采用在谱空间求解系数方程的方法，以充分利用快速傅立 叶变换和快速余弦变换高效的优点。 3.1 Chebyshev 伪谱方法求解[1,3,5] ˆ (1) 对任意求解变量φ及其系数φ ，引入 Chebyshev 多项式 ⎧ N π 1