dfdab.m

来自「AFD - Advanced Filter Design using MATLA」· M 代码 · 共 23 行

function [a, b] = dfdab(denD)

% dfdab DFD  Coefficients a and b
%
%   Authors: Dejan V. Tosic, Miroslav D. Lutovac, 1999.03.07
%                 [e-mail]
%                 [e-mail]
%
%   Copyright (c) 1999 by Tosic & Lutovac
%   $Revision: 1.0 $  $Date: 1999/03/07 03:07:42 $
%
%   References:
%        Miroslav D. Lutovac, Dejan V. Tosic, Brian L. Evans
%           Advanced Filter Design for Signal Processing
%                   Using MATLAB and Mathematica

p = sort(roots(denD));
b = abs(p).^2;
r = -2*real(p);
a = r./(b+1);

disp([a b])