Dear Ehsan Gharebaghi,

Well, I am not an expert on seismology so I do not think I can help you. However, in my opinion,  Matlab is a powerful tool for modeling. Firstly and importantly, you need to show the calculation method, anyway.

Regards, Kai LIU

I advice you explore CREWES MATLAB toolbox in the following URL

https://www.crewes.org/ResearchLinks/FreeSoftware/

There are some m-script to model acoustic and elastic wave ppropagation with finite differences or reflectivity approaches etc... I used it to solve direct problem (trace modelling) also with complex models. You can input a generic source in your model (also in depth) to simulate a source with a definite signature.

I hope this will useful for you