LTPDA_RI2FQ Convert complex pole/zero into frequency/Q pole/zero representation. [f0, q]= ltpda_ri2fq(c) M Hewitson 26-01-07 $Id: ltpda_ri2fq.m,v 1.4 2007/01/26 14:03:24 hewitson Exp $
0001 function [f0, q]= ltpda_ri2fq(c) 0002 0003 % LTPDA_RI2FQ Convert complex pole/zero into frequency/Q pole/zero representation. 0004 % 0005 % [f0, q]= ltpda_ri2fq(c) 0006 % 0007 % M Hewitson 26-01-07 0008 % 0009 % $Id: ltpda_ri2fq.m,v 1.4 2007/01/26 14:03:24 hewitson Exp $ 0010 % 0011 0012 0013 if(nargin==0) 0014 disp('usage: [f0, q]= ltpda_ri2fq(c)'); 0015 return 0016 end 0017 0018 a = real(c); 0019 b = imag(c); 0020 0021 f0 = -a*sqrt(b^2/a^2 + 1) / 2 / pi; 0022 q = sqrt(b^2/a^2 +1)/2; 0023 0024 % END