论坛 › 信号处理方法 › 数据进行频谱分析的实现程序

大食人鱼 发表于 2013-11-5 20:30

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

程序1 对地震记录数据进行频谱分析的实现程序load i.txt               %调入原始数据文件（加载地震波数据）Xt=i;                              %得到原始信号序列Fs=50;                           %采样频率为50Hzdt=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 onsubplot(3,1,2);                        %频域坐标方框图plot(f,abs(Xf));             %画出频域中的FFT波形图xlabel('频率/Hz');                      %X轴标示ylabel('频谱密度');                     % Y轴标示title('频域波形图');                     %加注标题 xlim();               %频域之画出采样频率的一半grid onsubplot(3,1,3);                         %频域坐标方框图semilogy(f,abs(Xf));    %画出频域中的FFT波形图，Y轴为对数xlabel('频率/Hz');                      %X轴标示ylabel('频谱密度取对数');               % Y轴标示title('频域波形图');                     %加注标题 xlim();                  %频域之画出采样频率的一半grid on

ChaChing 发表于 2013-11-5 21:22

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

查看完整版本: 数据进行频谱分析的实现程序