oder=10;
f=0:0.1:1;
m=[0 0 1 1 0 0 1 1 1 0 0];
[b,a]=yulewalk(oder,f,m)
[h,w]=freqz(b,a,128)
axes('position',[0.2 0.2 0.4 0.4]);
plot(f,m,'b-',w/pi,abs(h),'m--');
xlabel('frequency(pi)');
ylabel('magnitude');
title('direct IIR design-yulewalk');
legend('理想圖形','實(shí)際圖形',1);
grid