Matlab生成雷克子波

前言

雷克子波为地震数据处理中一个基本波

代码如下

close all;
fpeak=25;
n=100;
dt=0.001;
t=(0:n-1).*dt;
ricker=(1-2*(pi*fpeak*t).^2).*exp(-(pi.*fpeak.*t).^2);
num = 1000;
DataT = zeros(1,num);
for k=1:n
    DataT(k) = ricker(k);
end
for k=2:n
    DataT(num+2-k) = ricker(k);
end

subplot(2,1,1);
plot(DataT);

DataF = fft(DataT);
subplot(2,1,2);
plot(real(DataF));

生成效果

地震图像中众多雷克子波叠加后的案例