fig7-36.m

数字信号处理（机械工业出版社）的源码

% fig7-36
N = 8;                      	% Order of the filter (number of poles and zeros).
[Bh,Ah] = yulewalk(N,f,H);  	% Working, please wait.....
n = 256;
hh = freqz(Bh,Ah,n);        	% compute complex frequency response
hy  = abs(hh);             	% compute magnitude
ff  = fs/(2*n) * (0:n-1);
plot(fhz,H,ff,hy)
title('Actual vs. Desired Frequency Response')
xlabel('Frequency (Hz)')
ylabel('Magnitude')