%Appl7_1
clear all;
wp=[0.032 0.2];N=320;    %通帶邊界頻率（歸一化頻率）和濾波器階數
dt=0.02;   %中國數字地震臺網的采樣間隔為0.02s,采樣頻率為50Hz
b=fir1(N,wp,hanning(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 ChangChun.txt   %加載地震波形記錄
y=filtfilt(b,1,ChangChun); %產生零相位輸出，采用filtfilt 對輸入信號濾波
t=[0:length(ChangChun)-1]*dt;
figure(2)
subplot (2,1,1), plot (t,ChangChun),title ('輸入信號')%繪出輸入信號波形
subplot (2,1,2),plot (t,y)                  %繪出輸出信號波形
title('輸出信號'),xlabel('時間/s')