声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1192|回复: 1

[FFT] 数据进行频谱分析的实现程序

[复制链接]
发表于 2013-11-5 20:30 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
程序1 对地震记录数据进行频谱分析的实现程序
load i.txt               %调入原始数据文件(加载地震波数据)
Xt=i;                              %得到原始信号序列
Fs=50;                             %采样频率为50Hz
dt=1/Fs;                            %采样间隔(单位为s)
N=length(Xt);                       %原始信号序列长度
Xf=fft(Xt);         %对原始波形数据进行快速Fourier变换
n=0:N-1;
t=n*dt;                              %得到时间序列
f=n/(N*dt);                           %得到频率序列
subplot(3,1,1);                         %时域坐标方框图
plot(t, Xt);                       %画出时域中的原始波形图
xlabel('时间/s');                       %X轴标示
ylabel('振幅/counts');                   %Y轴标示
title('时域波形图');                     %加注标题
grid on
subplot(3,1,2);                        %频域坐标方框图
plot(f,abs(Xf));             %画出频域中的FFT波形图
xlabel('频率/Hz');                      %X轴标示
ylabel('频谱密度');                     % Y轴标示
title('频域波形图');                     %加注标题
xlim([0  Fs/2]);                 %频域之画出采样频率的一半
grid on
subplot(3,1,3);                         %频域坐标方框图
semilogy(f,abs(Xf));    %画出频域中的FFT波形图,Y轴为对数
xlabel('频率/Hz');                      %X轴标示
ylabel('频谱密度取对数');               % Y轴标示
title('频域波形图');                     %加注标题
xlim([0 Fs/2]);                  %频域之画出采样频率的一半
grid on

回复
分享到:

使用道具 举报

发表于 2013-11-5 21:22 | 显示全部楼层
1.分享或讨论问题?
2.abs(Xf))应该还不能算频谱密度吧!?
3.好奇问下,"研究方向 震颤信号", 什麼是震颤信号?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-6 09:48 , Processed in 0.067723 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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