%Appl3_3
load hns.dat ;                           %讀取數據序列
Xt=hns;                                 %把數據賦值給變量
Fs=50;                                   %設定采樣率 單位(Hz)
dt=1/Fs;                                 %求采樣間隔 單位(s)
N=length(Xt);                            %得到序列的長度
Xf=fft(Xt);                              %對信號進行快速Fourier變換(FFT)

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

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


figure(2)
load hns1.dat ;                           %讀取數據序列
Xt=hns1;                                 %把數據賦值給變量
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('時間域');
grid on;

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

figure(3)
load hns2.dat ;                           %讀取數據序列
Xt=hns2;                                 %把數據賦值給變量
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('時間域');
grid on;

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

figure(4)
load hns3.dat ;                           %讀取數據序列
Xt=hns3;                                 %把數據賦值給變量
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('時間域');
grid on;

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