%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;