Seismology and MCMC

Keywords:

MATLAB, Markov Chain Monte Carlo, Data Visualization

In this project, I aimed to fit the seismology related data into a MCMC model, by using three parameters from the dataset - SeismoEx.mat. The three variables I used to construct the model are: x (location), s (standard deviation, a vector full of ones) and y (the velocity data). I used the function inside SeisModel.m which implements the following model: $$m(v, D, x_0) = \frac{v}{\pi}arctan\left(\frac{x - x_0}{D}\right)$$
The project was designed as individual work. During the first half of the project, I solved the problem set by myself in MATLAB and later gathered together with my group to give out a 30 minutes presentation and Q&A session to the entire class collaboratively.

The following file is our project presentation: