%Appl7_2
%用快速Fourier變換（FFT）對(duì)數(shù)字信號(hào)進(jìn)行頻譜分析的程序
load grbx3.txt;                           %讀取數(shù)據(jù)序列
Xt=grbx3;                                 %把數(shù)據(jù)賦值給變量
Fs=50;                                   %設(shè)定采樣率 單位(Hz)
dt=1/Fs;                                 %求采樣間隔 單位(s)
N=length(Xt);                            %得到序列的長(zhǎng)度
Xf=fft(Xt);                              %對(duì)信號(hào)進(jìn)行快速Fourier變換(FFT)

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

subplot(2,1,2),plot([0:N-1]/(N*dt),abs(Xf)*2/N);%繪制信號(hào)的振幅譜
xlabel('頻率/Hz'),title('頻率域');
ylabel('振幅');
xlim([0 Fs/2]);                          %頻率軸只畫(huà)出Nyquist頻率之前的部分
grid on;