Finite difference modeling of acoustic waves in Matlab

Carrie F. Youzwishen, Gary F. Margrave

A Matlab toolkit,called the AFD package,has been written to model waves using acoustic finite differences.It uses central finite difference schemes to approximate derivatives to the scalar wave equation. Both a second order or 5 point approximation,and a fourth order or 9 point approximation,to the Laplacian are included.The fourth order approximation is slower,but is more accurate,and results in a broader temporal bandwidth.The AFD package is also equipped with absorbing boundary conditions to suppress reflections from the edges of the grid.The toolkit is able to create velocity models,shot records,exploding reflector models,as well as snapshots and movies of the wavefield propagating in depth.