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语言的
