%Appl7_5
clear all;
wp=2/25;N=400;    %通帶邊界頻率（歸一化頻率）和濾波器階數
dt=0.02;   %中國數字地震臺網的采樣間隔為0.02s,采樣頻率為50Hz
b=fir1(N,wp,hamming(N+1));           %設計FIR低通濾波器
figure(1)
[H,f]=freqz(b,1,512,1/dt);        %求出頻率特性
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;
load cz02.txt   %加載數據
x=cz02'; %輸入信號
y=filtfilt(b,1,x); %產生零相位輸出，采用filtfilt 對輸入信號濾波
t=[0:length(x)-1]*dt;
figure(2)
subplot (2,1,1), plot (t,x),title ('輸入信號')%繪出輸入信號波形
subplot (2,1,2),plot (t,y)                  %繪出輸出信號波形
title('輸出信號'),xlabel('時間/s')