exa4_22.m

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

% exa4-22_stmcb.m , for example 4-22 
% Design Steiglitz-McBride DF

clear all;
[b,a] = butter(6,0.2);
h = filter(b,a,[1 zeros(1,100)]);
freqz(b,a,128)
figure(2)
[bb,aa] = stmcb(h,4,4);
freqz(bb,aa,128)