加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

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

(2012-11-04 22:44:55)
标签:

杂谈

分类: 地球物理专题
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
    f(i)=10*(i-1);%计算频率
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。
雷克子波及其振幅谱的matlab输出

0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有