伪谱法求解波动方程算法的VC++实现

伪谱法利用傅主叶变换和有限差分求解波动方程，是三维波场正演模拟的重要方法之一。用VC 语言实现这一系统必须解决二阶空间偏导数、二阶时间偏导数的离散化问题。我们利用实信号的傅立叶变换争反傅立叶变换的基本性质即共轭性解决了空间偏导数的离散化问题，对于时间偏导数，我们则采用二阶中心差分的方法解决。最后模拟结果证明是有效的。