[sissigong]

:unsure:
小妹正在用九点格式的有限差分做频率域的数值模拟。但程序调了许久都没有得到理想的图形。感觉LU分解后解方程出来的波场值在震源附近比较理想，距离变大后就迅速衰减，几乎为零。

如果哪位达人曾经做过类似的模拟实验，可否与我讨论一下。多谢啦！ :Djavascript:emoticon(':D')

msn: gloria_gng@hotmail.com

[blueocean]

对这个问题的具体编程不是很熟, 你是用的pseudo-spectral方法吗? 记得83年的Geophysics有一篇文章介绍这个方法. 我只能问问, 你的远波场的距离多大?

[sissigong]

我用的是有限差分法，理论依据是R.G.Pratt 1990年在Geophysics上的文章。

由于只是实验，我并没有给很大的波场范围，水平方向2000m,垂直方向800m，空间间隔都是10m。