exa4_19.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 12 行

% 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,'*:实际频率响应')