exa4_23.m

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

% exa4-23_invfreqz.m , for example 4-23
% Design Steiglitz-McBride DF

clear all;
[b,a] = butter(5,0.2)
[h,w] = freqz(b,a,64);
[bv,av] = invfreqz(h,w,5,5)