nlp2bs.mws

syntfil[NLP2BS]  - frequency transformation of transfer function poles and zeros of NLP to a particular biquadratic transfer function parameters of band-rejection filter

Calling sequence:

      NLP2BS(fp, f_p, Gc, poles);

      NLP2BS(fp, f_p, Gc, poles, zeros);

Parameters:

      fp       - lower cutoff frequency [Hz]

      f_p      - upper cutoff frequency [Hz]

      Gc       - leading coefficient of denominator polynomial of transfer function H [-]

      poles    - 1D arrays of transfer function poles

      zeros    - 1D arrays of transfer function zeros (only for Cauer and Inverse Chebyshev transformations)

Parameters   fp , f_p   and   Gc   must be positive and fp  < f_p .

Description:

The function returns a table  with description of filter type, value of gain constant and parameters omega[0] , Q , ( omega[n] ) of particular passband biquadratic transfer functions - sections . Pole-zero pairing and section arranging are made so that dynamic balancing and noise properties of final filter were optimal. Section number determines sequence of the section in cascade. The filter type is cascade_BS . This specification is necessary for transfer function compilation and for following filter synthesis. Biquadratic transfer function with gain constant h[0] = 1  is denoted by:

Info level:

Setting of variable infolevel[syntfil]  can be used to get more detailed results.

infolevel[syntfil] =

      2  - print of gain constant and parameters omega[0] , Q , ( omega[n] ) of particular biquadratic transfer functions in the order of Q-factor magnitudes

      3  - print of filter type + print of level 2

Examples:

>    with(syntfil):

`Syntfil version 1.53 loaded`

`You can set infolevel[syntfil] variable to 2..5 to get more detailed results!`

>    x:=BS2NLP(1000,1200,3000,4000,3,30):

>    N:=ChebyshevNLPOrder(x);

>    Gc,poles:=ChebyshevPoles(N);

N := 5, 1.406250000, 3

Gc, poles := 15.96205352, vector([-.5485987094e-1+.9659274757*I, -.5485987094e-1-.9659274757*I, -.1436250067+.5969760105*I, -.1436250067-.5969760105*I, -.1775302716])
Gc, poles := 15.96205352, vector([-.5485987094e-1+.9659274757*I, -.5485987094e-1-.9659274757*I, -.1436250067+.5969760105*I, -.1436250067-.5969760105*I, -.1775302716])

>    NLP2BS(1000,4000,Gc,poles);

TABLE([1 = TABLE([omega_n = 12566.37062, Q = 2.744050378, omega_0 = 4527.003370]), 2 = TABLE([omega_n = 12566.37062, Q = 14.38685866, omega_0 = 25625.77433]), 3 = TABLE([omega_n = 12566.37062, Q = 14.3...
TABLE([1 = TABLE([omega_n = 12566.37062, Q = 2.744050378, omega_0 = 4527.003370]), 2 = TABLE([omega_n = 12566.37062, Q = 14.38685866, omega_0 = 25625.77433]), 3 = TABLE([omega_n = 12566.37062, Q = 14.3...
TABLE([1 = TABLE([omega_n = 12566.37062, Q = 2.744050378, omega_0 = 4527.003370]), 2 = TABLE([omega_n = 12566.37062, Q = 14.38685866, omega_0 = 25625.77433]), 3 = TABLE([omega_n = 12566.37062, Q = 14.3...
TABLE([1 = TABLE([omega_n = 12566.37062, Q = 2.744050378, omega_0 = 4527.003370]), 2 = TABLE([omega_n = 12566.37062, Q = 14.38685866, omega_0 = 25625.77433]), 3 = TABLE([omega_n = 12566.37062, Q = 14.3...
TABLE([1 = TABLE([omega_n = 12566.37062, Q = 2.744050378, omega_0 = 4527.003370]), 2 = TABLE([omega_n = 12566.37062, Q = 14.38685866, omega_0 = 25625.77433]), 3 = TABLE([omega_n = 12566.37062, Q = 14.3...
TABLE([1 = TABLE([omega_n = 12566.37062, Q = 2.744050378, omega_0 = 4527.003370]), 2 = TABLE([omega_n = 12566.37062, Q = 14.38685866, omega_0 = 25625.77433]), 3 = TABLE([omega_n = 12566.37062, Q = 14.3...
TABLE([1 = TABLE([omega_n = 12566.37062, Q = 2.744050378, omega_0 = 4527.003370]), 2 = TABLE([omega_n = 12566.37062, Q = 14.38685866, omega_0 = 25625.77433]), 3 = TABLE([omega_n = 12566.37062, Q = 14.3...

 

>    Gc2,poles2,zeros2:=InvChebyshevPolesZeros(N):

>    infolevel[syntfil]:=3:

>    NLP2BS(1000,4000,Gc2,poles2,zeros2);

type = cascade_BS

H0 = 1.947291

omega_0 = 12566.37062, Q = 0.94876, omega_n = 12566.37062

omega_0 =  7925.27924, Q = 1.29955, omega_n =  9229.99680

omega_0 = 19925.31313, Q = 1.29955, omega_n = 17108.74595

omega_0 =  6440.05130, Q = 3.78220, omega_n =  7716.45418

omega_0 = 24520.56091, Q = 3.78220, omega_n = 20464.53808

TABLE([1 = TABLE([omega_n = 9229.996805, Q = 1.299554986, omega_0 = 7925.279235]), 2 = TABLE([omega_n = 20464.53808, Q = 3.782204824, omega_0 = 24520.56091]), 3 = TABLE([omega_n = 7716.454182, Q = 3.78...
TABLE([1 = TABLE([omega_n = 9229.996805, Q = 1.299554986, omega_0 = 7925.279235]), 2 = TABLE([omega_n = 20464.53808, Q = 3.782204824, omega_0 = 24520.56091]), 3 = TABLE([omega_n = 7716.454182, Q = 3.78...
TABLE([1 = TABLE([omega_n = 9229.996805, Q = 1.299554986, omega_0 = 7925.279235]), 2 = TABLE([omega_n = 20464.53808, Q = 3.782204824, omega_0 = 24520.56091]), 3 = TABLE([omega_n = 7716.454182, Q = 3.78...
TABLE([1 = TABLE([omega_n = 9229.996805, Q = 1.299554986, omega_0 = 7925.279235]), 2 = TABLE([omega_n = 20464.53808, Q = 3.782204824, omega_0 = 24520.56091]), 3 = TABLE([omega_n = 7716.454182, Q = 3.78...
TABLE([1 = TABLE([omega_n = 9229.996805, Q = 1.299554986, omega_0 = 7925.279235]), 2 = TABLE([omega_n = 20464.53808, Q = 3.782204824, omega_0 = 24520.56091]), 3 = TABLE([omega_n = 7716.454182, Q = 3.78...
TABLE([1 = TABLE([omega_n = 9229.996805, Q = 1.299554986, omega_0 = 7925.279235]), 2 = TABLE([omega_n = 20464.53808, Q = 3.782204824, omega_0 = 24520.56091]), 3 = TABLE([omega_n = 7716.454182, Q = 3.78...
TABLE([1 = TABLE([omega_n = 9229.996805, Q = 1.299554986, omega_0 = 7925.279235]), 2 = TABLE([omega_n = 20464.53808, Q = 3.782204824, omega_0 = 24520.56091]), 3 = TABLE([omega_n = 7716.454182, Q = 3.78...

>   

See also:

ButterworthPoles, ChebyshevPoles, InvChebyshevPolesZeros, InvChebyshevBPolesZeros, CauerPolesZeros, CauerBPolesZeros, CauerCPolesZeros

NLP2LP, NLP2HP, NLP2BP, NLP2BP2

MakeH, ARCSynt, ARCBlock, ARCBlock1