%Appl7_4
figure(1)
load fjdos0351.txt;                   %讀取數據序列
Xt= fjdos0351;                   %把數據賦值給變量
Fs=50; dt=0.02;              %設定采樣率(單位Hz)和采樣間隔(單位s)
N=length(Xt); t=[0:N-1]*dt;      %得到序列的長度和時間序列
subplot(2,1,1),plot(t,Xt);        %繪制原始值序列
title('濾波前');
ylabel('振幅');
Fcp=1;  Fcs=2;   %設置的通帶邊界頻率，
Wp=Fcp/(Fs/2);  Ws=Fcs/(Fs/2);              %轉換為標準頻率
Wn=(Wp+Ws)/2;                    %給出設計時用的邊界頻率
wdelta=Ws-Wp;               %過渡帶寬
N=ceil(8*pi/wdelta);           %按Hamming窗求最小階數
pa=(N-1)/2/Fs;                 %計算相位延遲
b=fir1(N,Wn);           % 設計FIR帶阻濾波器
Yt=filtfilt(b,1,Xt);                  %對信號進行濾波
t=t-pa;                            %計算相位延遲
subplot(2,1,2),plot(t,Yt);           %繪制濾波后波形序列
xlim([0 300])
xlabel('時間/s');ylabel('振幅');title('濾波后信號');
figure(2)
[H,f]=freqz(b,1,128,50);
subplot(2,1,1),plot(f,20*log10(abs(H)))
xlabel('頻率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('頻率/Hz');ylabel('相位/^o');grid on;