exa4_20.m

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

% exa4-20_prony.m , for example 4-20
% Design prony DF

clear
[b,a]=butter(8,0.2)  % Butterworth滤波器设计
h=filter(b,a,[1,zeros(1,25)]); % 滤波器的冲激响应
[bp,ap]=prony(h,8,8)