appl7_2.m

MATLAB7.x数字信号处理 光盘内容

%Appl7_2
%用快速Fourier变换（FFT）对数字信号进行频谱分析的程序
load grbx3.txt;                           %读取数据序列
Xt=grbx3;                                 %把数据赋值给变量
Fs=50;                                   %设定采样率 单位(Hz)
dt=1/Fs;                                 %求采样间隔 单位(s)
N=length(Xt);                            %得到序列的长度
Xf=fft(Xt);                              %对信号进行快速Fourier变换(FFT)

subplot(2,1,1),plot([0:N-1]/Fs,Xt);        %绘制原始值序列
xlabel('时间/s'),title('时间域');
ylabel('振幅');
grid on;

subplot(2,1,2),plot([0:N-1]/(N*dt),abs(Xf)*2/N);%绘制信号的振幅谱
xlabel('频率/Hz'),title('频率域');
ylabel('振幅');
xlim([0 Fs/2]);                          %频率轴只画出Nyquist频率之前的部分
grid on;