c语言合成地震记录

matlab代码

clc;
clear all;
close all;
%1. generate the plane waves recorded at y = 300;
dt=0.002;
f=30;

[wavelet,tt1]= ricker(f,dt);

nt=450;
t=(0:nt-1)*dt;

nx = 100;dx=5;
d1 = zeros(nt, nx);

n1=40+50;
n2 = 100+50;
n3 = 250+50;
v1 = 1000;
v2 = 4000;
v3 = -5000;

x=[0:nx - 1]*dx;
data = zeros(nt,nx);
for ix = 1:nx
    g = zeros(nt,1);
    m1 = n1 + round( (ix-1)*dx/v1/dt);
    g(m1) = g(m1) + 1;
    
    m2 = n2 + round( (ix-1)*dx/v2/dt);
    g(m2) = g(m2) - .8;

    m3 = n3 + round( (ix-1)*dx/v3/dt);
    g(m3) = g(m3) + .5;
    
    s1 = conv(g, wavelet);
    data(:,ix) = s1(1:nt);
end
figure;imagesc(x,t, data)
set(gca,'Fontsize',16,'Fontweight','bold');
xlabel('X (m)');;ylabel('Time (s)');;
colorbar;
caxis([-1.5,1.5]);

怎么把这个改成c语言的