Articulo de referencia

Filtro elíptico

Un filtro elíptico (también conocido como filtro de Cauer , en honor a Wilhelm Cauer , o filtro de Zolotarev , en honor a Yegor Zolotarev ) es un filtro de procesamiento de seña...

Un filtro elíptico (también conocido como filtro de Cauer , en honor a Wilhelm Cauer , o filtro de Zolotarev , en honor a Yegor Zolotarev ) es un filtro de procesamiento de señales con comportamiento de rizado ecualizado (equiripple) tanto en la banda de paso como en la banda de rechazo . La cantidad de rizado en cada banda se puede ajustar de forma independiente, y ningún otro filtro de igual orden puede tener una transición de ganancia más rápida entre la banda de paso y la banda de rechazo , para los valores de rizado dados (ya sea que el rizado esté ecualizado o no). Alternativamente, se puede renunciar a la capacidad de ajustar de forma independiente el rizado de la banda de paso y la banda de rechazo, y en su lugar diseñar un filtro que sea lo más insensible posible a las variaciones de los componentes.

Cuando la ondulación en la banda de rechazo se aproxima a cero, el filtro se convierte en un filtro Chebyshev de tipo I. Cuando la ondulación en la banda de paso se aproxima a cero, el filtro se convierte en un filtro Chebyshev de tipo II y, finalmente, cuando ambos valores de ondulación se aproximan a cero, el filtro se convierte en un filtro Butterworth .

La ganancia de un filtro elíptico de paso bajo en función de la frecuencia angular ω viene dada por:

GRAMOnorte(ω)=11+ϵ2Rnorte2(ξ,ω/ω0){\displaystyle G_{n}(\omega )={1 \over {\sqrt {1+\epsilon ^{2}R_{n}^{2}(\xi ,\omega /\omega _{0})}}}}

donde R n es la función racional elíptica de orden n (a veces conocida como función racional de Chebyshev) y

ω0{\displaystyle \omega _{0}} es la frecuencia de corte
ϵ{\displaystyle \epsilon }es el factor de propagación
ξ{\displaystyle \xi }es el factor de selectividad

El valor del factor de rizado especifica el rizado de la banda de paso, mientras que la combinación del factor de rizado y el factor de selectividad especifica el rizado de la banda de rechazo.

Propiedades

Respuesta en frecuencia de un filtro paso bajo elíptico de cuarto orden con ε  =  0,5 y ξ  =  1,05. También se muestran la ganancia mínima en la banda de paso y la ganancia máxima en la banda de rechazo, así como la región de transición entre la frecuencia normalizada 1 y ξ.
Primer plano de la región de transición del gráfico anterior.
  • En la banda de paso, la función racional elíptica varía entre cero y uno. Por lo tanto, la ganancia de la banda de paso variará entre 1 y1/1+ϵ2{\displaystyle 1/{\sqrt {1+\epsilon ^{2}}}}.
  • En la banda de rechazo, la función racional elíptica varía entre el infinito y el factor de discriminación.Lnorte{\displaystyle L_{n}}que se define como:
Lnorte=Rnorte(ξ,ξ){\displaystyle L_{n}=R_{n}(\xi ,\xi )\,}
Por lo tanto, la ganancia de la banda de rechazo variará entre 0 y1/1+ϵ2Lnorte2{\displaystyle 1/{\sqrt {1+\epsilon ^{2}L_{n}^{2}}}}.
  • En el límite deξ{\displaystyle \xi \rightarrow \infty }La función racional elíptica se convierte en un polinomio de Chebyshev y, por lo tanto, el filtro se convierte en un filtro de Chebyshev de tipo I , con factor de rizado ε.
  • Dado que el filtro de Butterworth es una forma límite del filtro de Chebyshev, se deduce que en el límite deξ{\displaystyle \xi \rightarrow \infty },ω00{\displaystyle \omega _ {0}\rightarrow 0}yϵ0{\displaystyle \epsilon \rightarrow 0}de tal manera queϵRnorte(ξ,1/ω0)=1{\displaystyle \epsilon \,R_{n}(\xi ,1/\omega _{0})=1}El filtro se convierte en un filtro Butterworth.
  • En el límite deξ{\displaystyle \xi \rightarrow \infty },ϵ0{\displaystyle \epsilon \rightarrow 0}y ω00{\displaystyle \omega _ {0}\rightarrow 0}de tal manera queξω0=1{\displaystyle \xi \omega _ {0}=1}yϵLnorte=α{\displaystyle \epsilon L_{n}=\alpha }, el filtro se convierte en un filtro Chebyshev tipo II con ganancia
GRAMO(ω)=11+1α2Tnorte2(1/ω){\displaystyle G(\omega )={\frac {1}{\sqrt {1+{\frac {1}{\alpha ^{2}T_{n}^{2}(1/\omega )}}}}}}

Polos y ceros

Logaritmo del valor absoluto de la ganancia de un filtro elíptico de octavo orden en el espacio de frecuencia compleja (s = σ + jω) con ε = 0,5, ξ = 1,05 y ω₀ = 1. Los puntos blancos representan polos y los puntos negros, ceros. Hay un total de 16 polos y 8 ceros dobles. Lo que parece ser un solo polo y un cero cerca de la región de transición son en realidad cuatro polos y dos ceros dobles, como se muestra en la vista ampliada a continuación. En esta imagen, el negro corresponde a una ganancia de 0,0001 o menos, y el blanco a una ganancia de 10 o más.
Una vista ampliada de la región de transición de la imagen anterior, que permite distinguir los cuatro polos y los dos ceros dobles.

Los ceros de la ganancia de un filtro elíptico coincidirán con los polos de la función racional elíptica, que se derivan en el artículo sobre funciones racionales elípticas .

Los polos de la ganancia de un filtro elíptico se pueden derivar de una manera muy similar a la derivación de los polos de la ganancia de un filtro Chebyshev de tipo I. Para simplificar, supongamos que la frecuencia de corte es igual a la unidad. Los polos(ωpagmetro){\displaystyle (\omega _ {pm})}de la ganancia del filtro elíptico serán los ceros del denominador de la ganancia. Usando la frecuencia complejas=σ+jω{\displaystyle s=\sigma +j\omega }Esto significa que:

1+ϵ2Rnorte2(js,ξ)=0{\displaystyle 1+\epsilon ^{2}R_{n}^{2}(-js,\xi )=0\,}

Definiciónjs=dod(w,1/ξ){\displaystyle -js=\mathrm {cd} (w,1/\xi)}donde cd() es la función coseno elíptica de Jacobi y, utilizando la definición de las funciones racionales elípticas, se obtiene:

1+ϵ2dod2(nortewKnorteK,1Lnorte)=0{\displaystyle 1+\epsilon ^{2}\mathrm {cd} ^{2}\left({\frac {nwK_{n}}{K}},{\frac {1}{L_{n}}}\right)=0\,}

dóndeK=K(1/ξ){\displaystyle K=K(1/\xi )}yKnorte=K(1/Lnorte){\displaystyle K_{n}=K(1/L_{n})}. Resolviendo para w

w=KnorteKnortedod1(±jϵ,1Lnorte)+metroKnorte{\displaystyle w={\frac {K}{nK_{n}}}\mathrm {cd} ^{-1}\left({\frac {\pm j}{\epsilon }},{\frac {1}{L_{n}}}\right)+{\frac {mK}{n}}}

donde los múltiples valores de la función inversa cd() se hacen explícitos usando el índice entero m .

Los polos de la función de ganancia elíptica son entonces:

spagmetro=idod(w,1/ξ){\displaystyle s_{pm}=i\,\mathrm {cd} (w,1/\xi )\,}

Como ocurre con los polinomios de Chebyshev, esto puede expresarse en forma explícitamente compleja ( Lutovac y otros, 2001 , § 12.8).

spagmetro=a+jbdo{\displaystyle s_{pm}={\frac {a+jb}{c}}}
a=ζnorte1ζnorte21incógnitametro21incógnitametro2/ξ2{\displaystyle a=-\zeta _{n}{\sqrt {1-\zeta _{n}^{2}}}{\sqrt {1-x_{m}^{2}}}{\sqrt {1-x_{m}^{2}/\xi ^{2}}}}
b=incógnitametro1ζnorte2(11/ξ2){\displaystyle b=x_{m}{\sqrt {1-\zeta _{n}^{2}(1-1/\xi ^{2})}}}
do=1ζnorte2+incógnitai2ζnorte2/ξ2{\displaystyle c=1-\zeta _ {n}^{2}+x_ {i}^{2}\zeta _ {n}^{2}/\xi ^{2}}

dóndeζnorte{\displaystyle \zeta _{n}}es una función denorte,ϵ{\displaystyle n,\,\epsilon }yξ{\displaystyle \xi }yincógnitametro{\displaystyle x_{m}}son los ceros de la función racional elíptica.ζnorte{\displaystyle \zeta _{n}}es expresable para todo n en términos de funciones elípticas de Jacobi, o algebraicamente para algunos órdenes, especialmente los órdenes 1, 2 y 3. Para los órdenes 1 y 2 tenemos

ζ1=11+ϵ2{\displaystyle \zeta _{1}={\frac {1}{\sqrt {1+\epsilon ^{2}}}}}
ζ2=2(1+t)1+ϵ2+(1t)2+ϵ2(1+t)2{\displaystyle \zeta _{2}={\frac {2}{(1+t){\sqrt {1+\epsilon ^{2}}}+{\sqrt {(1-t)^{2}+\epsilon ^{2}(1+t)^{2}}}}}}

dónde

t=11/ξ2{\displaystyle t={\sqrt {1-1/\xi ^{2}}}}

La expresión algebraica paraζ3{\displaystyle \zeta _{3}}es bastante complejo (véase Lutovac y otros (2001 , § 12.8.1) ).

La propiedad de anidamiento de las funciones racionales elípticas se puede utilizar para construir expresiones de orden superior paraζnorte{\displaystyle \zeta _{n}}:

ζmetronorte(ξ,ϵ)=ζmetro(ξ,1ζnorte2(Lmetro,ϵ)1){\displaystyle \zeta _{m\cdot n}(\xi ,\epsilon )=\zeta _{m}\left(\xi ,{\sqrt {{\frac {1}{\zeta _{n}^{2}(L_{m},\epsilon )}}-1}}\right)}

dóndeLmetro=Rmetro(ξ,ξ){\displaystyle L_{m}=R_{m}(\xi ,\xi )}.

Consideraciones de diseño

Factores Q normalizados de los polos de un filtro elíptico de octavo orden con ξ  =  1,1 en función del factor de rizado ε . Cada curva representa cuatro polos, ya que los pares de polos complejos conjugados y los pares de polos positivo-negativo tienen el mismo factor Q. (Las curvas azul y cian casi coinciden). El factor Q de todos los polos se minimiza simultáneamente en ε Qmin  =  1  / L n = 0,02323...   

Véase Lutovac y et al. (2001 , § 12.11, 13.14) .

Los filtros elípticos generalmente se especifican requiriendo un valor particular para la ondulación de la banda de paso, α p , la ondulación de la banda de rechazo, α s, y la nitidez del corte. Estos especifican un orden mínimo n para la función elíptica racional que describe el comportamiento del filtro. Específicamente, consideremos un filtro que deja pasar frecuencias por debajo de ω p y rechaza frecuencias por encima de ω s , y definamos la función auxiliar=1211(ωpagωs)241+1(ωpagωs)24{\displaystyle u={\frac {1}{2}}\cdot {\frac {1-{\sqrt[{4}]{1-\left({\frac {\omega _{p}}{\omega _{s}}}\right)^{2}}}}{1+{\sqrt[{4}]{1-\left({\frac {\omega _{p}}{\omega _{s}}}\right)^{2}}}}}}Para restringir la ganancia de la banda de paso a [1- G p , 1] y la ganancia de la banda de rechazo a [0, 1- G s ] , el orden n debe satisfacer [ 1 ]norteregistro(16GRAMOsGRAMOpag)registro(+25+159+15013)Dónde:GRAMOs=10αs/101GRAMOpag=10αpag/101{\displaystyle {\begin{aligned}&n\geq {\frac {\log {\left(16{\frac {G_{s}}{G_{p}}}\right)}}{-\log(u+2u^{5}+15u^{9}+150u^{13})}}\\&{\text{Where:}}\\&Gs=10^{\alpha _{s/10}}-1\\&Gp=10^{\alpha _{p/10}}-1\\\end{aligned}}}Esa ecuación surge de una aproximación a la función integral elíptica completa K; una formulación exacta es [ 2 ].norteK(1τ12)K(τ2)K(τ1)K(1τ22){\displaystyle n\geq {\frac {K\left({\sqrt {1-\tau _{1}^{2}}}\right)K(\tau _{2})}{K(\tau _{1})K\left({\sqrt {1-\tau _{2}^{2}}}\right)}}}dóndeτ1=GRAMOpagGRAMOsτ2=ωpagωs{\displaystyle {\begin{aligned}\tau _{1}&={\sqrt {\frac {G_{p}}{G_{s}}}}\\\tau _{2}&={\frac {\omega _{p}}{\omega _{s}}}\end{aligned}}}Otra consideración de diseño es la sensibilidad de la función de ganancia a los valores de los componentes electrónicos utilizados para construir el filtro. Esta sensibilidad es inversamente proporcional al factor de calidad ( factor Q ) de los polos de la función de transferencia del filtro . El factor Q de un polo se define como:

Q=|spagmetro|2Rmi(spagmetro)=12porque(arg(spagmetro)){\displaystyle Q=-{\frac {|s_{pm}|}{2\mathrm {Re} (s_{pm})}}=-{\frac {1}{2\cos(\arg(s_{pm}))}}}

y es una medida de la influencia del polo en la función de ganancia. Para un filtro elíptico, sucede que, para un orden dado, existe una relación entre el factor de rizado y el factor de selectividad que minimiza simultáneamente el factor Q de todos los polos en la función de transferencia:

ϵQmetroinorte=1Lnorte(ξ){\displaystyle \epsilon _{Qmin}={\frac {1}{\sqrt {L_{n}(\xi )}}}}

Esto da como resultado un filtro que es mínimamente insensible a las variaciones de los componentes, pero se pierde la capacidad de especificar de forma independiente las ondulaciones de la banda de paso y la banda de rechazo. Para dichos filtros, a medida que aumenta el orden, la ondulación en ambas bandas disminuirá y la tasa de corte aumentará. Si se decide utilizar un filtro elíptico de Q mínimo para lograr una ondulación mínima particular en las bandas del filtro junto con una tasa de corte particular, el orden necesario generalmente será mayor que el orden que se necesitaría de otro modo sin la restricción de Q mínimo. Una imagen del valor absoluto de la ganancia se verá muy similar a la imagen de la sección anterior excepto que los polos están dispuestos en un círculo en lugar de una elipse. No estarán espaciados uniformemente y habrá ceros en el eje ω, a diferencia del filtro Butterworth , cuyos polos están dispuestos en un círculo espaciado uniformemente sin ceros.

Comparación con otros filtros lineales

Aquí se muestra una imagen del filtro elíptico junto a otros tipos de filtros comunes obtenidos con el mismo número de coeficientes:

Como se puede apreciar en la imagen, los filtros elípticos son más nítidos que todos los demás, pero muestran ondulaciones en todo el ancho de banda.

Construcción a partir de ceros de transmisión de Chebyshev

Las bandas de rechazo de filtros elípticos son esencialmente filtros de Chebyshev con ceros de transmisión donde los ceros de transmisión están dispuestos de manera que producen una banda de rechazo de rizado uniforme. Dado esto, es posible convertir una ecuación característica de filtro de Chebyshev,K(s){\displaystyle K(s)}que contiene los ceros de reflexión de Chebyshev en el numerador y ningún cero de transmisión en el denominador, a un filtro elípticoK(s){\displaystyle K(s)}que contiene los ceros de reflexión elíptica en el numerador y los ceros de transmisión elíptica en el denominador, creando iterativamente ceros de transmisión a partir del inverso escalado de los ceros de reflexión de Chebyshev, y luego restableciendo una banda de paso de Chebyshev de ondulación uniforme a partir de los ceros de transmisión, y repitiendo hasta que las iteraciones no produzcan más cambios significativos.K(s){\displaystyle K(s)}. [ 3 ] El factor de escala utilizado,Ωdo{\displaystyle \Omega _{c}}, es la relación de frecuencia de corte de la banda de rechazo a la banda de paso y también se conoce como el inverso del "factor de selectividad". [ 1 ] Dado que los diseños elípticos generalmente se especifican a partir de los requisitos de atenuación de la banda de rechazo,Ωdo{\displaystyle \Omega _{c}}, se pueden derivar de las ecuaciones que establecen el orden mínimo, n, arriba.

elωs/ωpag{\displaystyle \omega _{s}/\omega _{p}}relación,Ωdo{\displaystyle \Omega _{c}}puede derivarse trabajando el problema de orden mínimo, n , anterior hacia atrás desde n para encontrarΩdo{\displaystyle \Omega _{c}}. [ 1 ]

norte= número de polos (orden del filtro)αpag y αs=atenuación de la banda de paso y de la banda de rechazo, respectivamente.ωpag y ωs=frecuencias de banda de paso y banda de rechazo, respectivamenteD=GRAMOsGRAMOpag=10αs/10110αpag/101q=(16D)(1/norte)0=q++25+159+15013= raíz real de la ecuación anteriork= el factor de selectividad =1(121+2)4Ωdo=ωsωpag=1k=11(121+2)4{\displaystyle {\begin{aligned}n&={\text{ number of poles (order of the filter)}}\\\alpha _{p}{\text{ and }}\alpha _{s}&={\text{pass band and stop band attenuation, respectively}}\\\omega _{p}{\text{ and }}\omega _{s}&={\text{pass band and stop band frequencies, respectively}}\\D&={\frac {G_{s}}{G_{p}}}={\frac {10^{\alpha _{s}/10}-1}{10^{\alpha _{p}/10}-1}}\\q&=(16D)^{(-1/n)}\\0&=-q+u+2u^{5}+15u^{9}+150u^{13}\\u&={\text{ real root of above equation}}\\k&={\text{ the selectivity factor }}={\sqrt {1-{\bigg (}{\frac {1-2u}{1+2u}}{\bigg )}^{4}}}\\\Omega _{c}&={\frac {\omega _{s}}{\omega _{p}}}={\frac {1}{k}}={\frac {1}{\sqrt {1-{\bigg (}{\frac {1-2u}{1+2u}}{\bigg )}^{4}}}}\\\end{aligned}}}

Los polinomios característicos,K(s){\displaystyle K(s)}calculado a partir deΩdo{\displaystyle \Omega _{c}}y los requisitos de atenuación, pueden luego traducirse a los polinomios de la función de transferencia,GRAMO(s){\displaystyle G(s)}con la traducción clásica,GRAMO(s)=1/(1+ε2K(s)K(s))|postes LHP{\displaystyle G(s)={\sqrt {1/(1+\varepsilon ^{2}K(s)K(-s))}}{\bigg |}_{\text{LHP poles}}}dóndeε2=10Apag/10.1.{\displaystyle \varepsilon ^{2}=10^{Ap/10.}-1.}yApag{\displaystyle A_{p}} es la ondulación de la banda de paso. [ 3 ] [ 4 ]

Ejemplo sencillo

Diseñar un filtro elíptico con una ondulación de banda de paso de 1  dB de 0 a 1 rad/seg y una ondulación de banda de rechazo de 40  dB de al menos 1,25 rad/seg a{\displaystyle \infty }.

Aplicando los cálculos anteriores para el valor de n antes de aplicar la función ceil() , se encuentra que n es 4.83721900 redondeado al siguiente entero, 5, al aplicar la función ceil() , lo que significa que se requiere un filtro elíptico de 5 polos para cumplir con los requisitos de diseño especificados. Aplicando los cálculos anteriores paraΩdo{\displaystyle \Omega _{c}}era necesario diseñar una banda de rechazo con  una atenuación de exactamente 40 dB,Ωdo{\displaystyle \Omega _{c}}se encuentra que es 1,2186824.

La función de inversión escalada polinómica se puede realizar trasladando cada raíz, s , aΩdo/s{\displaystyle \Omega _{c}/s}, lo cual se puede lograr fácilmente invirtiendo el polinomio y escalándolo porΩdo{\displaystyle \Omega _{c}}, como se muestra.

asnorte+bsnorte1dos2+ds1+mis0(miΩdonorte)snorte+(dΩdonorte1)snorte1(doΩdo2)s2+(bΩdo1)s1+(aΩdo0)s0{\displaystyle {\begin{aligned}&as^{n}+bs^{n-1}\dots cs^{2}+ds^{1}+es^{0}\Longrightarrow ({\frac {e}{\Omega _{c}^{n}}})s^{n}+({\frac {d}{\Omega _{c}^{n-1}}})s^{n-1}\dots ({\frac {c}{\Omega _{c}^{2}}})s^{2}+({\frac {b}{\Omega _{c}^{1}}})s^{1}+({\frac {a}{\Omega _{c}^{0}}})s^{0}\\\end{aligned}}}

Los pasos de diseño elíptico son los siguientes: [ 3 ]

  1. Diseñar filtros Chebyshev con una  ondulación de banda de paso de 1 dB.
  2. Invierte todos los ceros de reflexión sobreΩdo{\displaystyle \Omega _{c}}para crear ceros de transmisión
  3. Cree una banda de paso de rizado uniforme a partir de los ceros de transmisión utilizando el proceso descrito en los ceros de transmisión de Chebyshev.
  4. Repita los pasos 2 y 3 hasta que tanto la banda de paso como la banda de rechazo ya no varíen de forma apreciable. Normalmente, entre 15 y 25 iteraciones producen diferencias en los coeficientes del orden de 1e-15.

Para ilustrar los pasos, las siguientes ecuaciones K(s) comienzan con una K(s) de Chebyshev estándar y luego iteran a través del proceso. Se observan diferencias visibles en las tres primeras iteraciones. Al llegar a las 18 iteraciones, las diferencias en K(s) se vuelven insignificantes. Las iteraciones pueden interrumpirse cuando el cambio en los coeficientes K(s) sea suficientemente pequeño como para cumplir con los requisitos de precisión de diseño. Todas las siguientes iteraciones de K(s) se han normalizado de tal manera que|K(j)|=1{\displaystyle |K(j)|=1}Sin embargo, este paso puede posponerse hasta la última iteración, si se desea.

iteración 0: K(s)=16s5+20s3+5s1iteración 1: K(s)=9.2965947s5+12.999133s3+4.0025668s0,14167325s4+0,84164496s2+1iteración 2: K(s)=8.6496472s5+12.270597s3+3.8746611s0,19518773s4+0,94147634s2+1iteración 17: K(s)=8.550086786383502s5+12.157269873073034s3+3.854163602012615s0,2043607336740334s4+0,9573802183509494s2+1iteración 18: K(s)=8.550086786383422s5+12.157269873072942s3+3.854163602012599s0,2043607336740334s4+0,9573802183509494s2+1{\displaystyle {\begin{aligned}{\text{iteration 0: }}K(s)&={\frac {16s^{5}+20s^{3}+5s}{1}}\\{\text{iteration 1: }}K(s)&={\frac {9.2965947s^{5}+12.999133s^{3}+4.0025668s}{0.14167325s^{4}+0.84164496s^{2}+1}}\\{\text{iteration 2: }}K(s)&={\frac {8.6496472s^{5}+12.270597s^{3}+3.8746611s}{0.19518773s^{4}+0.94147634s^{2}+1}}\\&\vdots \\{\text{iteration 17: }}K(s)&={\frac {8.550086786383502s^{5}+12.157269873073034s^{3}+3.854163602012615s}{0.2043607336740334s^{4}+0.9573802183509494s^{2}+1}}\\{\text{iteration 18: }}K(s)&={\frac {8.550086786383422s^{5}+12.157269873072942s^{3}+3.854163602012599s}{0.2043607336740334s^{4}+0.9573802183509494s^{2}+1}}\\\end{aligned}}}

Para encontrar elGRAMO(s){\displaystyle G(s)}función de transferencia, haga lo siguiente. [ 3 ]

ε2=101dB/10.1.=.25892541GRAMO(s)=GRAMO(s)GRAMO(s)|postes LHP=11+ε2K(s)K(s)|postes LHP=K(s)dminorteK(s)dminorteK(s)dminorteK(s)dminorte+ε2K(s)nortemetroK(s)nortemetro|postes LHP=0,20436073s4+0,95738022s2+118.928479s1053.78661s854.942632s622.939175s41.931467+1|postes LHP{\displaystyle {\begin{aligned}\varepsilon ^{2}&=10^{1dB/10.}-1.=.25892541\\G(s)&={\sqrt {G(s)G(-s)}}{\bigg |}_{\text{LHP poles}}={\sqrt {\frac {1}{1+\varepsilon ^{2}K(s)K(-s)}}}{\bigg |}_{\text{LHP poles}}={\sqrt {\frac {K(s)_{den}K(-s)_{den}}{K(s)_{den}K(-s)_{den}+\varepsilon ^{2}K(s)_{num}K(-s)_{num}}}}{\bigg |}_{\text{LHP poles}}\\&={\frac {0.20436073s^{4}+0.95738022s^{2}+1}{{\sqrt {-18.928479s^{10}-53.78661s^{8}-54.942632s^{6}-22.939175s^{4}-1.931467+1}}{\bigg |}_{\text{LHP poles}}}}\\\end{aligned}}}

Para obtenerGRAMO(s){\displaystyle G(s)}Desde el semiplano izquierdo, factoriza el numerador y el denominador para obtener las raíces usando un algoritmo de búsqueda de raíces . Descarta todas las raíces del semiplano derecho del denominador, reduce a la mitad las raíces repetidas en el numerador y reconstruye.GRAMO(s){\displaystyle G(s)}con las raíces restantes. [ 3 ] [ 4 ] Generalmente, normalizar|GRAMO(s)|{\displaystyle |G(s)|}a 1 ens=0{\displaystyle s=0}.

GRAMO(s)=0,20436073s4+0,95738022s2+14.3506872s5+4.0174213s4+8.0362343s3+4.9129149s2+3.4288915s+1{\displaystyle {\begin{aligned}&G(s)={\frac {0.20436073s^{4}+0.95738022s^{2}+1}{4.3506872s^{5}+4.0174213s^{4}+8.0362343s^{3}+4.9129149s^{2}+3.4288915s+1}}\\\end{aligned}}}

Para confirmar que el ejemploGRAMO(s){\displaystyle G(s)}es correcto, el gráfico deGRAMO(s){\displaystyle G(s)}a lo largo dejω{\displaystyle j\omega }Se muestra a continuación con una ondulación de banda de paso de 1  dB, una frecuencia de corte de 1 rad/seg y una atenuación de banda de parada de 40  dB que comienza en 1,21868 rad/seg.

Simulación elíptica de cinco polos

Incluso modificaciones de pedidos

Los filtros elípticos de orden par implementados con elementos pasivos, típicamente inductores, capacitores y líneas de transmisión, con terminaciones de igual valor en cada lado, no pueden implementarse con la función de transferencia elíptica tradicional sin el uso de bobinas acopladas, lo cual puede no ser deseable o factible. Esto se debe a la incapacidad física de acomodar los ceros de reflexión y transmisión de Chebyshev de orden par que dan como resultado valores de la matriz de dispersión S12 que exceden el valor de S12 enω=0{\displaystyle \omega =0}y los valores finitos de S12 que existen enω={\displaystyle \omega =\infty }. Si no es factible diseñar el filtro con una de las terminaciones aumentada o disminuida para acomodar la banda de paso S12, entonces la función de transferencia elíptica debe modificarse para mover el cero de reflexión de orden par más bajo aω=0{\displaystyle \omega =0}y la transmisión de orden par más alta de cero aω={\displaystyle \omega =\infty }manteniendo la respuesta de rizado uniforme de la banda de paso y la banda de rechazo. [ 5 ]

La modificación necesaria implica mapear cada polo y cero de la función de transferencia elíptica de manera que el cero de reflexión de frecuencia más baja se mapee a cero, el cero de transmisión de frecuencia más alta se mapee a{\displaystyle \infty }y los polos y ceros restantes según sea necesario para mantener la banda de paso y la banda de rechazo de rizado uniforme. El cero de reflexión de frecuencia más baja se puede encontrar factorizando elK(s){\displaystyle K(s)}numerador, y el cero de transmisión de frecuencia más alta se puede encontrar factorizando elK(s){\displaystyle K(s)}denominador.

Para trasladar los ceros de reflexión, se aplica la siguiente ecuación a todos los polos y ceros deK(s){\displaystyle K(s)}. [ 5 ] Si bien en teoría, las operaciones de traducción pueden realizarse en cualquiera de los dosK(s){\displaystyle K(s)}oGRAMO(s){\displaystyle G(s)}, los ceros de reflexión deben extraerse deK(s){\displaystyle K(s)}, por lo que generalmente es más eficiente realizar las operaciones de traducción enK(s){\displaystyle K(s)}.

Ri=Ri2+ωLO21ωLO2{\displaystyle R_{i}'={\sqrt {\frac {R_{i}^{2}+\omega _{LO}^{2}}{1-\omega _{LO}^{2}}}}}

Dónde:

Ri{\displaystyle R_{i}}es el cero o polo de la función elíptica original

Ri{\displaystyle R_{i}'}es el cero o polo mapeado para la función de transferencia de orden par modificada.

ωLO{\displaystyle \omega _{LO}}es el cero de reflexión de frecuencia más baja en la banda de paso.

El signo del componente imaginario deRi{\displaystyle R_{i}'}está determinado por el signo del originalRi{\displaystyle R_{i}}.

Para trasladar los ceros de transmisión, se aplica la siguiente ecuación a todos los polos y ceros deK(s){\displaystyle K(s)}. [ 5 ] Si bien en teoría, las operaciones de traducción pueden realizarse en cualquiera de los dosK(s){\displaystyle K(s)}oGRAMO(s){\displaystyle G(s)}, si los ceros de reflexión deben extraerse deK(s){\displaystyle K(s)}, puede ser más eficiente realizar las operaciones de traducción enK(s){\displaystyle K(s)}.

Ri=(ωHI21)Ri2ωHI2+Ri2{\displaystyle R_{i}'={\sqrt {\frac {(\omega _{HI}^{2}-1)R_{i}^{2}}{\omega _{HI}^{2}+R_{i}^{2}}}}}

Dónde:

Ri{\displaystyle R_{i}}es el cero o polo de la función elíptica original

Ri{\displaystyle R_{i}'}es el cero o polo mapeado para la función de transferencia de orden par modificada.

ωHI{\displaystyle \omega _{HI}}es el cero de transmisión de frecuencia más alta en la banda de paso.

El signo del componente imaginario deRi{\displaystyle R_{i}'}está determinado por el signo del originalRi{\displaystyle R_{i}}. Si se opera enGRAMO(s){\displaystyle G(s)}el signo del componente real deRi{\displaystyle R_{i}'}debe ser negativo para cumplir con el requisito del semiplano izquierdo.

Es importante señalar que todas las aplicaciones requieren traducciones tanto de paso como de parada. Los diplexores de red pasivos, por ejemplo, solo requieren traducciones de banda de parada de orden par y funcionan de manera más eficiente con bandas de paso de orden par sin traducir. [ 5 ]

CuandoGRAMO(s){\displaystyle G(s)}Una vez completado, se crea una función de transferencia de ondulación uniforme con valores de matriz de dispersión para S12 de 1.ω=0{\displaystyle \omega =0}y 0 enω={\displaystyle \omega =\infty }, que puede implementarse con redes pasivas terminadas de forma igual.

La ilustración a continuación muestra un filtro elíptico de octavo orden modificado para admitir redes pasivas terminadas de orden par mediante la reubicación del cero de reflexión de frecuencia más baja de una frecuencia finita a 0 y el cero de transmisión de frecuencia más alta a{\displaystyle \infty }manteniendo una respuesta en frecuencia de banda de paso y banda de rechazo con ondulación uniforme .

Ilustración elíptica modificada de orden par
Ilustración elíptica modificada de orden par

ElΩdo{\displaystyle \Omega _{c}}y el cálculo del orden en el párrafo de construcción elíptica anterior son solo para filtros elípticos no modificados. Aunque incluso las modificaciones del orden no tienen efecto en la atenuación de la banda de paso o de la banda de rechazo, se esperan pequeños errores en el orden yΩdo{\displaystyle \Omega _{c}}cálculos. Por lo tanto, es importante aplicar modificaciones de orden par después de todosK(s){\displaystyle K(s)}Las iteraciones se completan si se desea preservar las atenuaciones de la banda de paso y de parada. Si la función elíptica modificada de orden par se crea a partir de unaΩdo{\displaystyle \Omega _{c}}requisito, el realΩdo{\displaystyle \Omega _{c}}será ligeramente más grande que el diseñoΩdo{\displaystyle \Omega _{c}}. Asimismo, un cálculo de orden, n , puede dar como resultado un valor menor que el orden realmente requerido.

Implementación del reloj de arena

Un filtro de reloj de arena es un caso especial de filtro donde los ceros de reflexión son el recíproco de los ceros de transmisión alrededor de una frecuencia de atenuación de corte normalizada de 3,01  dB de 1 rad/seg, lo que resulta en que todos los polos del filtro residen en el círculo unitario . [ 6 ] La implementación elíptica de reloj de arena tiene una ventaja sobre un filtro de Chebyshev inverso en que la banda de paso es más plana, y tiene una ventaja sobre los filtros elípticos tradicionales en que el retardo de grupo tiene un pico menos pronunciado en la frecuencia de corte.

Respuesta de frecuencia recíproca S11 y S12 en forma de reloj de arena
Respuesta en frecuencia recíproca S11 y S12 de reloj de arena de 7 polos

Proceso de síntesis

La forma más sencilla de sintetizar un filtro de reloj de arena es diseñar un filtro elíptico con una atenuación de banda de parada de diseño especificada, As , y una atenuación de banda de paso calculada que cumpla con el requisito de red de dos puertos sin pérdidas de que los parámetros de dispersión|S11|2+|S12|2=1{\displaystyle |S_{11}|^{2}+|S_{12}|^{2}=1}. [ 7 ] Junto con la conocida magnitud dB a traducción aritmética,(Sij)dB=20logramo10(|Sij|arith){\displaystyle (S_{ij})_{dB}=20log_{10}(|S_{ij}|_{arith})}, la manipulación algebraica produce el siguiente requisito de atenuación de banda de paso calculado.

Apag=10registro10(1.10(As/10)){\displaystyle A_{p}=-10\log _{10}{(1.-10^{(-A_{s}/10)})}}

El parámetro A p , definido anteriormente, producirá ceros de reflexión y transmisión recíprocos alrededor de una  frecuencia de corte de 3,01 dB aún desconocida. Para diseñar un filtro elíptico con una frecuencia de banda de paso de 1 rad/s,  es necesario determinar la frecuencia de atenuación de 3,01 dB y utilizarla para escalar inversamente los polinomios de diseño elíptico. El resultado serán polinomios con una atenuación de 3,01  dB a una frecuencia normalizada de 1 rad/s. Se puede utilizar el método de Newton o resolver las ecuaciones directamente con un algoritmo de búsqueda de raíces para determinar la  frecuencia de atenuación de 3,01 dB.

Escalado de frecuencias con el método de Newton

SiGRAMO(s){\displaystyle G(s)}es la función de transferencia de reloj de arena para encontrar la  frecuencia de 3,01 dB, yωdo{\displaystyle \omega _{c}}es la  frecuencia de 3 dB a encontrar, los pasos a continuación se pueden utilizar para encontrarlaωdo{\displaystyle \omega _{c}}

  1. SiGRAMO(s)GRAMO(s){\displaystyle G(s)G(-s)}no está ya disponible, multiplicarGRAMO(s){\displaystyle G(s)}porGRAMO(s){\displaystyle G(-s)}para obtenerGRAMO(s)GRAMO(s){\displaystyle G(s)G(-s)}.
  2. negar todos los términos desnorte{\displaystyle s^{n}}cuando(norte+2){\displaystyle (n+2)}es divisible por4{\displaystyle 4}Eso serías2{\displaystyle s^{2}},s6{\displaystyle s^{6}},s10{\displaystyle s^{10}}y así sucesivamente. La función modificada se llamará GRAMO2(s)GRAMO2(s){\displaystyle G_{2}(s)G_{2}(-s)}y esta modificación permitirá el uso de números reales en lugar de números complejos al evaluar el polinomio y su derivada.ωa{\displaystyle \omega _{a}}ahora se puede utilizar en lugar del complejojωa{\displaystyle j\omega _{a}}
  3. Convierta la atenuación deseada en dB,AdB{\displaystyle A_{dB}}, a un valor de ganancia aritmética al cuadrado,Barith2{\displaystyle B_{arith}^{2}}, mediante el usoBarith2=10AdB/10{\displaystyle B_{arith}^{2}=10^{A_{dB}/10}}Por ejemplo, 3,010  dB se convierte en 0,5, 1  dB se convierte en 0,79432823 y así sucesivamente.
  4. Calcular el modificado|GRAMO2(s)GRAMO2(s)|{\displaystyle |G_{2}(s)G_{2}(-s)|}en el método de Newton usando el valor real,ωa{\displaystyle \omega _{a}}Siempre tome el valor absoluto.
  5. Calcula la derivada modificadaGRAMO2(ωa)GRAMO2(ωa){\displaystyle G_{2}(\omega _{a})G_{2}(-\omega _{a})}con respecto al valor real,ωa{\displaystyle \omega _{a}}NO tome el valor absoluto de la derivada.

Cuando se completan los pasos 1) al 4), la expresión que involucra el método de Newton se puede escribir como:

ωa=ωa([GRAMO2(ωa)GRAMO2(ωa)|B2)/(d[GRAMO2(ωa)GRAMO2(ωa)]/dωa){\displaystyle \omega _{a}=\omega _{a}-([G_{2}(\omega _{a})G_{2}(-\omega _{a})|-B^{2})/(d[G_{2}(\omega _{a})G_{2}(-\omega _{a})]/d\omega _{a})}

utilizando un valor real paraωa{\displaystyle \omega _{a}}sin necesidad de aritmética compleja. El movimiento deωa{\displaystyle \omega _{a}}debe limitarse para evitar que se vuelva negativo al principio de las iteraciones para aumentar la fiabilidad. Cuando la convergencia esté completa,ωa{\displaystyle \omega _{a}}puede utilizarse para elωdo{\displaystyle \omega _{c}}que se puede utilizar para escalar el originalGRAMO(s){\displaystyle G(s)}Denominador de la función de transferencia. La atenuación de la modificadaGRAMO(s){\displaystyle G(s)}El valor será prácticamente el deseado a 1 rad/s. Si se realiza correctamente, solo se necesitan unas pocas iteraciones para ajustar la atenuación en un amplio rango de valores deseados, tanto para filtros de orden pequeño como de orden muy grande.

Escalamiento de frecuencia con búsqueda de raíces

Desde|GRAMO(jωa)|{\displaystyle |G(j\omega _{a})|}no contiene ninguna información de fase, factorizar directamente la función de transferencia no producirá resultados utilizables. Sin embargo, la función de transferencia puede modificarse multiplicándola porGRAMO(s){\displaystyle G(-s)}para eliminar todos los poderes extraños deGRAMO(jωa){\displaystyle G(j\omega a)}, lo que a su vez obligaGRAMO(jωa){\displaystyle G(j\omega a)}ser real en todas las frecuencias, y luego encontrar la frecuencia que resulta en el cuadrado de la atención deseada.

  1. SiGRAMO(s)GRAMO(s){\displaystyle G(s)G(-s)}no está ya disponible, multiplicarGRAMO(s){\displaystyle G(s)}porGRAMO(s){\displaystyle G(-s)}para obtenerGRAMO(s)GRAMO(s){\displaystyle G(s)G(-s)}.
  2. Convierta la atenuación deseada en dB,AdB{\displaystyle A_{dB}}, a un valor de ganancia aritmética al cuadrado,Barith2{\displaystyle B_{arith}^{2}}, mediante el usoBarith2=10AdB/10{\displaystyle B_{arith}^{2}=10^{A_{dB}/10}}Por ejemplo, 3,010  dB se convierte en 0,5, 1  dB se convierte en 0,79432823 y así sucesivamente.
  3. EncontrarPAG(S)=GRAMOnortemetro(S)GRAMOnortemetro(S)Barith2GRAMOdminorte(S)GRAMOdminorte(S){\displaystyle P(S)=G_{num}(S)G_{num}(-S)-B_{arith}^{2}G_{den}(S)G_{den}(-S)}
  4. Encuentra las raíces de P(S) usando un algoritmo de búsqueda de raíces.
  5. Del conjunto de raíces anteriores, seleccione la raíz imaginaria positiva para todos los filtros de orden y la raíz real positiva para los filtros de orden par.ωdo{\displaystyle \omega _{c}}.

Escalado de la función de transferencia

Cuandoωdo{\displaystyle \omega _{c}}Una vez determinado, el polinomio de la función de transferencia de Hourglass se puede escalar de la siguiente manera:

GRAMO(s)original=nortenortenortesnortenorte++norte2s2+norte1s+norte0Dnortenortesnorted++D2s2+D1s+D0 (polinomios de función de transferencia originales sin escalar)GRAMO(s) escamoso=nortenortenorte(ωdos)nortenorte++norte2(ωdos)2+norte1ωdos+norte0Dnortenorte(ωdos)norted++D2(ωdos)2+D1ωdos+D0 (3.01 dB a 1 rad/seg (polinomios de función de transferencia escalados)ωdo=|GRAMO(s)| Frecuencia de atenuación de 3,01 dB nortenorte,norted= orden del numerador y del denominador, respectivamentenorte,D= coeficientes del numerador y del denominador, respectivamente{\displaystyle {\begin{aligned}G(s)_{\text{orig}}&={\frac {N_{nn}s^{nn}+\dots +N_{2}s^{2}+N_{1}s+N_{0}}{D_{nn}s^{nd}+\dots +D_{2}s^{2}+D_{1}s+D_{0}}}{\text{ (original unscaled transfer function polynomials)}}\\G(s)_{\text{ scaled}}&={\frac {N_{nn}(\omega _{c}s)^{nn}+\dots +N_{2}(\omega _{c}s)^{2}+N_{1}\omega _{c}s+N_{0}}{D_{nn}(\omega _{c}s)^{nd}+\dots +D_{2}(\omega _{c}s)^{2}+D_{1}\omega _{c}s+D_{0}}}{\text{ (}}3.01{\text{ dB at 1 rad/sec scaled transfer function polynomials)}}\\\omega _{c}&=|G(s)|{\text{ 3.01 dB attenuation frequency }}\\nn,nd&={\text{ order of numerator and denominator, respectively}}\\N,D&={\text{ coefficients of numerator and denominator, respectively}}\\\end{aligned}}}

Incluso modificaciones de pedidos

Los filtros de reloj de arena de orden par presentan las mismas limitaciones respecto a las redes pasivas con terminación igual que otros filtros elípticos. Las mismas modificaciones de orden par que resuelven el problema con los filtros elípticos también lo resuelven con los filtros de reloj de arena.

Referencias

  1. 1 2 3 Rorabaugh, C. Britton (1 de enero de 1993). Manual del diseñador de filtros digitales (  Edición reimpresa). Blue Ridge Summit, PA, EE. UU.: Tab Books, División de McGraw-Hill, Inc. págs.  93 a 95. ISBN 978-0830644315.
  2. Paarmann, Larry D. (2001). Diseño y análisis de filtros analógicos: una perspectiva de procesamiento de señales . Norwell, Massachusetts, EE. UU.: Kluwer Academic Publishers. págs. 182 a 198. ISBN  0-7923-7373-1.
  3. 1 2 3 4 5 Apuntes de clase sobre diseño de filtros del Dr. Byron Bennett, archivados el 23/04/2024 en la Wayback Machine , 1985, Universidad Estatal de Montana . Archivados el 28/03/2023 en la Wayback Machine , Departamento de Ingeniería Eléctrica . Archivados el 15/04/2024 en la Wayback Machine , Bozeman , Montana, EE. UU.
  4. 1 2 Sedra, Adel S.; Brackett, Peter O. (1978). Filter Theory and Design: Active and Passive . Beaverton, Ohio, EE. UU.: Matrix Publishers, Inc. págs. 45–73 . ISBN  978-0916460143.
  5. ^ Saal , Rudolf (enero de 1979 ) . Manual de diseño de filtros (en inglés y alemán) (1ª ed.). Múnich, Alemania: Allgemeine Elektricitais-Gesellschaft. págs.25 , 26, 56–61 , 116, 117. ISBN   3-87087-070-2.
  6. Bennett, Byron J. (diciembre de 1988). "Una nueva técnica de síntesis de filtros: el reloj de arena" . IEEE Transactions on Circuits and Systems . 35 (12): 1469– 1477. doi : 10.1109/31.9910 vía IEEE.
  7. Matthaei, George L.; Young, Leo; Jones, EMT (1984). Filtros de microondas, redes de adaptación de impedancia y estructuras de acoplamiento . 610 Washington Street, Dedham, Massachusetts, EE. UU.: Artech House, Inc. (publicado en 1985). pág. 44. ISBN  0-89006-099-1.{{cite book}}: CS1 mantenimiento: ubicación ( enlace )
  • Daniels, Richard W. (1974). Métodos de aproximación para el diseño de filtros electrónicos . Nueva York: McGraw-Hill. ISBN 0-07-015308-6.
  • Lutovac, Miroslav D.; Tosic, Dejan V.; Evans, Brian L. (2001). Diseño de filtros para el procesamiento de señales con MATLAB y Mathematica . Nueva Jersey, EE. UU.: Prentice Hall. ISBN 0-201-36130-2.