Introduction - If you have any usage issues, please Google them yourself
MATLABScript for obtaning the analytical solution for 2D wave propagation in a viscoelastic medium, based on the Appendix B of paper of Carcione et al. (1988)(1), corrected in (2).
The source used is a Ricker wavelet in the y velocity component.
A figure of the source temporal function and its inverse Fourier transform are provided to verify the good working of the method. Also works for elastic variables by substituting the values of vp and vs by real constant ones (vp_0 & vs_0).
Attenuation is implemented by GMB-EK rheology presented in (4) and thoroughfully explained in (3).
In order to make a proper comparison with numerical data it is recommended to either use the same rheology in both the numerical and the analytical. It is also possible to get an almost constant Q value by increasing a lot the number of mechanisms used in both the numerical and the analytical solutions, thus minimizing the effects of the Q fitting in the comparison.