Write a matlab code to do the following:

Compare your results to this example.

3-iv) Now, experiment reproducing a noiseless synthetic trace, using a Ricker wavelet and a reflectivity series represented as a series of impulse responses located at the double traveltime value ton, of each reflection event. A Ricker wavelet can be calculated using Where fis the dominant frequency and to is the lag or delay. Write down a code to generate a Ricker wavelet 256 points evenly spaced and delay it so the first break of the wavelet coincides with time zero. Convolve the Ricker wavelet with a reflectivity series, r, (dimension or length equal to 1,024) and whose non-zero values (0.25,-0.12, 0.08, 0.31) are located at the indexes 149, 299, 511, and 724. Use the function r=zeros (N) to create a vector r with N zero values and then add the non- zero values manually. You can also use delta (n-n0) over a sequence of indexes n at the location nO (see command list at the end of this document. Use Matlab function w = conv(u, v) %%Krone ker delta at index 2 for a vector of 11 components >>clear all >> delta-inline ('n 0') >delta (n-2) >stem (n,x) 8 in this case, 2 in the argument of delta correspond to the index at which the Krone ker delta is located at. To place a reflection 88coefficient, rl, at index 2, simply write down; >r1*delta (n-2) 8you can build a whole reflectivity series by adding as many terms as needed at locations nj: >>xrl*delta (n-nl)+ r2*delta (n-n2)+ r3*delta (n-n3)+