论文部分内容阅读
波动方程数值模拟的有限差分格式有隐式差分格式和显式差分格式两种，各有优点和缺点。针对高维波动方程提出了一种新的隐式分裂有限差分格式。其基本原理是：首先将高维波动方程按传播方向分解为一系列的一维波动问题，然后分别沿各方向隐式求解。该格式包含了X，Y，Z三个方向相互独立的一维隐式差分格式，每个方向的一维格式在数值离散后归结为一个三对角矩阵问题，可以用追赶法快速地求解。将该格式从时间一空间域变换至时间一波数域，证明此格式可以通过适当地选取参数来提高计算精度，保证计算过程的稳定性和与八阶显式差分格式同样的频散特性