找回密码
 注册
Simdroid-非首页
查看: 272|回复: 6

请教MATLAB中如何画地震波的傅里叶幅值谱?

[复制链接]
发表于 2008-11-15 20:17:39 | 显示全部楼层 |阅读模式 来自 上海
请教如题:MATLAB中如何画地震波的傅里叶幅值谱?
发表于 2008-11-17 09:02:41 | 显示全部楼层 来自 北京
Simdroid开发平台
原帖由 zzxzd1984 于 2008-11-15 20:17 发表
请教如题:MATLAB中如何画地震波的傅里叶幅值谱?


H=FFT(x,n)
然后
plot(abs(H))

具体含义自己在帮助中看一下

命令窗中输入
help fft

即可看到如何使用

评分

1

查看全部评分

回复 不支持

使用道具 举报

 楼主| 发表于 2008-12-30 22:46:06 | 显示全部楼层 来自 上海

回复 2# ljelly 的帖子

谢谢!
回复 不支持

使用道具 举报

发表于 2013-9-18 10:30:33 | 显示全部楼层 来自 江西南昌
来晚了!!

评分

1

查看全部评分

回复 不支持

使用道具 举报

发表于 2013-9-23 08:11:55 | 显示全部楼层 来自 台湾
ljelly 发表于 2008-11-17 09:02
H=FFT(x,n)
然后
plot(abs(H))

一般傅里叶幅值谱还需 *2/n
但针对地震波是否有特殊定义, 个人不确定
不过以为应该不会, 不然幅值将随n而不同

评分

1

查看全部评分

回复 不支持

使用道具 举报

发表于 2013-12-5 00:24:02 | 显示全部楼层 来自 江苏南京
ChaChing 发表于 2013-9-23 08:11
一般傅里叶幅值谱还需 *2/n
但针对地震波是否有特殊定义, 个人不确定
不过以为应该不会, 不然幅值将随n而 ...

说得没错,对地震波也一样。
求地震波x的傅里叶幅值谱程序:
x=        ;   %地震波数据
dt=0.02       ; %采样时间间隔, s,
N=length(x);
X=2/N*abs(fft(x));
k=1:fix(N/2)+1;
fk=k/(N*dt);
plot(fk,X(1:length(k)))
xlabel('频率'),ylabel('傅里叶振幅')

评分

1

查看全部评分

回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 22:30 , Processed in 0.049365 second(s), 20 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表