fir2iir.m

来自「Matlab时频分析工具箱,希望能对大家有所帮助啊」· M 代码 · 共 47 行

function g=fir2iir(g,Liir);
%FIR2IIR   Extend FIR window to IIR
%   Usage:  g=fir2iir(g,Liir);
%
%   FIR2IIR(g,Liir) will extend the FIR window g to a length Liir window by
%   inserting zeros. Note that this is a slightly different behaviour
%   than MIDDLEPAD.
%
%   FIR2IIR can also be used to extend a FIR window to a longer FIR
%   window, for instance in order to satisfy that the window length is
%   divisable by the number of channels.
%
%   SEE ALSO:  MIDDLEPAD

% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
% 
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
% 
% You should have received a copy of the GNU General Public License
% along with this program.  If not, see <http://www.gnu.org/licenses/>.

error(nargchk(2,2,nargin));

Lfir=length(g);


if Lfir>Liir
  error('Liir must be larger than length of window.');
end;

if rem(Lfir,2)==0
  % HPE middlepad works the same way as the FIR extension (e.g. just
  % inserting zeros) for even-length signals.
  g=middlepad(g,Liir,[],.5);
else
  % WPE middlepad works the same way as the FIR extension (e.g. just
  % inserting zeros) for odd-length signals.
  g=middlepad(g,Liir);
end;