二维粘弹性介质中波的传播的解析形式。应用Carcione等在1988年发表的《WAVE PROPAGATION SIMULATION IN A LINEAR VISCOELASTIC MEDIUM》一文中的方法，更加详细的介绍见英文描述。该代码对于做波场模拟的同学有很好的参考价值。- MATLABscr ipt 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.