Convert comlpex pole/zero into frequency/Q pole/zero representation. [f0, q]= pri2fq(c) M Hewitson 2003-03-12 $Id: pri2fq.html,v 1.14 2008/03/31 10:27:36 hewitson Exp $
0001 function [f0, q]= pri2fq(c) 0002 0003 % Convert comlpex pole/zero into frequency/Q pole/zero representation. 0004 % 0005 % [f0, q]= pri2fq(c) 0006 % 0007 % M Hewitson 2003-03-12 0008 % 0009 % 0010 % $Id: pri2fq.html,v 1.14 2008/03/31 10:27:36 hewitson Exp $ 0011 % 0012 0013 0014 if(nargin==0) 0015 disp('usage: [f0, q]= ri2fq(c)'); 0016 return 0017 end 0018 0019 a = real(c(1)); 0020 b = imag(c(1)); 0021 0022 % f0 = -a*sqrt(b^2/a^2 + 1) / 2 / pi; 0023 % q = sqrt(b^2/a^2 +1)/2; 0024 0025 k = -a*sqrt(b^2/a^2 + 1) * 2 * pi; 0026 f0 = 1/k; 0027 q = sqrt(b^2/a^2 +1)/2; 0028 0029 % END