Home > m > helper > ltpda_fq2ri2.m

ltpda_fq2ri2

PURPOSE ^

LTPDA_FQ2RI2 Convert frequency/Q pole/zero representation into real

SYNOPSIS ^

function ri= ltpda_fq2ri2(f0, Q)

DESCRIPTION ^

 LTPDA_FQ2RI2 Convert frequency/Q pole/zero representation into real
 and imaginary conjugate pairs. Returns [ri conj(ri)]
 
 ri= ltpda_fq2ri2(f0, Q)
 
 M Hewitson 26-01-07
 

 $Id: ltpda_fq2ri2.html,v 1.1 2007/06/08 14:15:10 hewitson Exp $

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function ri= ltpda_fq2ri2(f0, Q)
0002 
0003 % LTPDA_FQ2RI2 Convert frequency/Q pole/zero representation into real
0004 % and imaginary conjugate pairs. Returns [ri conj(ri)]
0005 %
0006 % ri= ltpda_fq2ri2(f0, Q)
0007 %
0008 % M Hewitson 26-01-07
0009 %
0010 %
0011 % $Id: ltpda_fq2ri2.html,v 1.1 2007/06/08 14:15:10 hewitson Exp $
0012 %
0013 
0014 if(nargin==0)
0015   disp('usage: ri = ltpda_fq2ri2(f0, Q)')
0016   return
0017 end
0018 
0019 % a = Q - 0.5 ;
0020 % re = 2*pi*f0 / sqrt(1+tan(a)^2);
0021 % im = re * tan(a);
0022 
0023 w0 = 2*pi*f0;
0024 re = -w0/(2*Q);
0025 im = w0*sqrt(4*Q*Q-1)/(2*Q);
0026 G = 1/(w0*w0);
0027 tri = complex(re,im);
0028 ri = [tri conj(tri)]';
0029 
0030 % END

Generated on Fri 08-Jun-2007 16:09:11 by m2html © 2003