Single-file vectorized implementation of acoustic wave propagation in isotropic medium in MATLAB. The program is solving second-order wave equation in displacement formulation. Modelling area is surrounded by simple absorbing sponge boundaries with exponential decay (Cerjan, 1985).
- Finite-Differences in Time Domain (FDTD)
- Explicit time stepping
- O(2,2)
- Conventional stencils derived from Taylor series:
- in space [1: -2 :1]/dx2
- in time [1: -2 :1]/dt2
- Isotropic acoustic (vp)
- Sponge or reflecting boundaries
Run acoustic_3D_FDTD_O22.m in MATLAB