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

Pseudo-spectral forward Modeling Algorithm Written in VC++ Realization of Solving Wave Equation

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

The pseudo-spectrum method uses Fourier transform and finite difference to solve wave equation.It is one of the important methods used to simulate the seismic wave field.Using VC++ language to realize this system,we must solve the question that disperse the second order spatial differential coefficient and temporal differential coefficient.We use the character of Fourier transform and Inverse Fourier transform (conjugate character) to solve the difficult problem and temporal differential coefficient, and we use the method of second order center finite difference to solve the question.At the end ,the result of simulating is satisfactory.