eg10_13.m

来自「matlab7.0从入门到精通(光盘内容)」· M 代码 · 共 13 行

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('理想图形','实际图形',1);
grid