雷克子波及其振幅谱的matlab输出

标签:
杂谈 |
分类: 地球物理专题 |
fm=50;%主频
dt=0.001;%时间域采样间隔
number=100;%采样点数
t=-number/2+1:number/2;
a=(1-2*(pi*fm*t*0.001).^2).*exp(-(pi*fm*t*0.001).^2);
subplot(2,1,1);
plot(t,a);
title('Ricker-雷克子波');
xlabel('时间t(ms)');
ylabel('幅值A');
for i=1:100
end
Y=abs(fft(a));%fourier变换,取振幅谱
subplot(2,1,2)
plot(f,Y);
title('Ricker子波的振幅谱');
xlabel('频率f(hz)');
ylabel('振幅谱');
%注:时间域采样间隔为0.001s,采样点数为100点,总的时间长度为0.1s,则频率域采样间隔为10hz。
前一篇:天天向上(9)