% exa4-19_yulewalk.m , for example 4-19
% Design Yule-Walk DF

clear all;
f = [0 0.6 0.6 1];
m = [1 1 0 0];
[b,a] = yulewalk(10,f,m);
[h,w] = freqz(b,a,256);
plot(f,m,'-',w/pi,abs(h),'*')
text(0.7,1.2,'-:理想頻率響應')
text(0.7,1.1,'*:實際頻率響應')