%Samp6_12
Order=10;   %濾波器的階數
f=0:0.1:1;   %歸一化頻率點
m=[0 0 1 1 0 0  1 1 1 0 0];  %幅度點
[b,a]=yulewalk(Order,f,m);  %設計濾波器
[h,w]=freqz(b,a,128);   %計算128個點的頻率特性
figure(1)
plot(f,m,'b-',w/pi,abs(h),'k:');  %繪制理想濾波器和設計濾波器的幅頻特性
xlabel('歸一化頻率');ylabel('振幅');
title('運用yulewalk方法設計IIR濾波器');
legend('理想特性','實際設計',1);  %給定圖例
figure(2)
Fs=50;   %信號采樣頻率
f1=6;f2=17;  %信號的頻率成分
N=100; %數據點數
dt=1/Fs;n=0:N-1;t=n*dt; %時間序列
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t);  %輸入信號
subplot(2,1,1),plot(t,x),title('輸入信號') %繪制輸入信號
y=filtfilt(b,a,x);  %對信號進行濾波
subplot(2,1,2),plot(t,y)   %繪制輸出信號
title('輸出信號'),xlabel('時間/s')