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

大食人鱼

程序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

ChaChing

1.分享或讨论问题?
2.abs(Xf))应该还不能算频谱密度吧!?
3.好奇问下,"研究方向 震颤信号", 什麼是震颤信号?