%Appl6_4
Fcp=4;Fcs=6;   %設置邊界頻率為4Hz，6Hz
load iir_signaldata.txt; %調入數據文件
Xt=iir_signaldata; %原始波形數據
Fs=50; %采樣頻率50Hz
wp=Fcp*2/Fs;ws=Fcs*2/Fs; %將通帶臨界頻率Fcp和阻帶臨界頻率Fcs轉換為歸一化頻率
Rp=1;Rs=25; %通帶衰減和阻帶衰減
Nn=128; %繪頻譜圖所用點數
[N,Wn]=buttord(wp,ws,Rp,Rs);  %求濾波器的最小階數和臨界頻率（歸一化頻率）
[b,a]=butter(N,Wn);  %求Butterworth數字濾波器傳遞函數的分子b和分母a
figure(1); %圖形(一)
[H,f]=freqz(b,a,Nn,Fs); %濾波器特性圖
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;
n=0:length(Xt)-1; t=n/Fs; %轉換為時間序列
figure(2); %圖形(二)
subplot(2,1,1),plot(t,Xt);title('輸入信號'); %繪制輸入信號
xlabel('時間/s');ylabel('振幅'); %坐標軸標識
Yt=filter(b,a,Xt); %對輸入信號進行濾波
subplot(2,1,2),plot(t,Yt),title('輸出信號');%繪制輸出信號
xlabel('時間/s');ylabel('振幅'); %坐標軸標識