以下是一个使用有限差分法模拟三维地震波的Python代码:
import numpy as np
定义空间和时间步长
dx = 0.1 dt = 0.001
定义空间范围
xmin = 0. xmax = 1.
定义时间范围
tmin = 0. tmax = 1.
定义空间和时间网格
x = np.arange(xmin, xmax+dx, dx) t = np.arange(tmin, tmax+dt, dt)
定义地震波速度
v = 1.
定义空间和时间步长
dx2 = dxdx dt2 = dtdt
定义空间和时间网格
u = np.ze