example7_4.m

来自「用dsp解压mp3程序的算法」· M 代码 · 共 11 行

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Example7_4.m - Program for example 7-4
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

[b,a] = ellip(4,1,50,.3); % design an ellipe filter
h = impz(b,a,30);         % impulse response

[bb,aa] = prony(h,4,4);   % design using Prony method
h1 = impz(bb,aa,30);      % impulse response