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:
donde R n es la función racional elíptica de orden n (a veces conocida como función racional de Chebyshev) y
- es la frecuencia de corte
- es el factor de propagación
- 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


- 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 y.
- En la banda de rechazo, la función racional elíptica varía entre el infinito y el factor de discriminación.que se define como:
- Por lo tanto, la ganancia de la banda de rechazo variará entre 0 y.
- En el límite deLa 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,yde tal manera queEl filtro se convierte en un filtro Butterworth.
- En el límite de,y de tal manera quey, el filtro se convierte en un filtro Chebyshev tipo II con ganancia
Polos y ceros


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 polosde la ganancia del filtro elíptico serán los ceros del denominador de la ganancia. Usando la frecuencia complejaEsto significa que:
Definicióndonde cd() es la función coseno elíptica de Jacobi y, utilizando la definición de las funciones racionales elípticas, se obtiene:
dóndey. Resolviendo para w
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:
Como ocurre con los polinomios de Chebyshev, esto puede expresarse en forma explícitamente compleja ( Lutovac y otros, 2001 , § 12.8).
dóndees una función deyyson los ceros de la función racional elíptica.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
dónde
La expresión algebraica paraes 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:
dónde.
Consideraciones de diseño

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 auxiliarPara 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 ]Esa ecuación surge de una aproximación a la función integral elíptica completa K; una formulación exacta es [ 2 ].dóndeOtra 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:
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:
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,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ípticoque 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.. [ 3 ] El factor de escala utilizado,, 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,, se pueden derivar de las ecuaciones que establecen el orden mínimo, n, arriba.
elrelación,puede derivarse trabajando el problema de orden mínimo, n , anterior hacia atrás desde n para encontrar. [ 1 ]
Los polinomios característicos,calculado a partir dey los requisitos de atenuación, pueden luego traducirse a los polinomios de la función de transferencia,con la traducción clásica,dóndey 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.
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 paraera necesario diseñar una banda de rechazo con una atenuación de exactamente 40 dB,se encuentra que es 1,2186824.
La función de inversión escalada polinómica se puede realizar trasladando cada raíz, s , a, lo cual se puede lograr fácilmente invirtiendo el polinomio y escalándolo por, como se muestra.
Los pasos de diseño elíptico son los siguientes: [ 3 ]
- Diseñar filtros Chebyshev con una ondulación de banda de paso de 1 dB.
- Invierte todos los ceros de reflexión sobrepara crear ceros de transmisión
- 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.
- 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 queSin embargo, este paso puede posponerse hasta la última iteración, si se desea.
Para encontrar elfunción de transferencia, haga lo siguiente. [ 3 ]
Para obtenerDesde 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.con las raíces restantes. [ 3 ] [ 4 ] Generalmente, normalizara 1 en.
Para confirmar que el ejemploes correcto, el gráfico dea lo largo deSe 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.

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 eny los valores finitos de S12 que existen en. 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 ay la transmisión de orden par más alta de cero amanteniendo 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 ay 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 elnumerador, y el cero de transmisión de frecuencia más alta se puede encontrar factorizando eldenominador.
Para trasladar los ceros de reflexión, se aplica la siguiente ecuación a todos los polos y ceros de. [ 5 ] Si bien en teoría, las operaciones de traducción pueden realizarse en cualquiera de los doso, los ceros de reflexión deben extraerse de, por lo que generalmente es más eficiente realizar las operaciones de traducción en.
Dónde:
es el cero o polo de la función elíptica original
es el cero o polo mapeado para la función de transferencia de orden par modificada.
es el cero de reflexión de frecuencia más baja en la banda de paso.
El signo del componente imaginario deestá determinado por el signo del original.
Para trasladar los ceros de transmisión, se aplica la siguiente ecuación a todos los polos y ceros de. [ 5 ] Si bien en teoría, las operaciones de traducción pueden realizarse en cualquiera de los doso, si los ceros de reflexión deben extraerse de, puede ser más eficiente realizar las operaciones de traducción en.
Dónde:
es el cero o polo de la función elíptica original
es el cero o polo mapeado para la función de transferencia de orden par modificada.
es el cero de transmisión de frecuencia más alta en la banda de paso.
El signo del componente imaginario deestá determinado por el signo del original. Si se opera enel signo del componente real dedebe 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 ]
CuandoUna vez completado, se crea una función de transferencia de ondulación uniforme con valores de matriz de dispersión para S12 de 1.y 0 en, 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 amanteniendo una respuesta en frecuencia de banda de paso y banda de rechazo con ondulación uniforme .

Ely 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 ycálculos. Por lo tanto, es importante aplicar modificaciones de orden par después de todosLas 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 unarequisito, el realserá ligeramente más grande que el diseño. 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.

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. [ 7 ] Junto con la conocida magnitud dB a traducción aritmética,, la manipulación algebraica produce el siguiente requisito de atenuación de banda de paso calculado.
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
Sies la función de transferencia de reloj de arena para encontrar la frecuencia de 3,01 dB, yes la frecuencia de 3 dB a encontrar, los pasos a continuación se pueden utilizar para encontrarla
- Sino está ya disponible, multiplicarporpara obtener.
- negar todos los términos decuandoes divisible porEso sería,,y así sucesivamente. La función modificada se llamará y esta modificación permitirá el uso de números reales en lugar de números complejos al evaluar el polinomio y su derivada.ahora se puede utilizar en lugar del complejo
- Convierta la atenuación deseada en dB,, a un valor de ganancia aritmética al cuadrado,, mediante el usoPor ejemplo, 3,010 dB se convierte en 0,5, 1 dB se convierte en 0,79432823 y así sucesivamente.
- Calcular el modificadoen el método de Newton usando el valor real,Siempre tome el valor absoluto.
- Calcula la derivada modificadacon respecto al valor real,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:
utilizando un valor real parasin necesidad de aritmética compleja. El movimiento dedebe limitarse para evitar que se vuelva negativo al principio de las iteraciones para aumentar la fiabilidad. Cuando la convergencia esté completa,puede utilizarse para elque se puede utilizar para escalar el originalDenominador de la función de transferencia. La atenuación de la modificadaEl 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
Desdeno 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 porpara eliminar todos los poderes extraños de, lo que a su vez obligaser real en todas las frecuencias, y luego encontrar la frecuencia que resulta en el cuadrado de la atención deseada.
- Sino está ya disponible, multiplicarporpara obtener.
- Convierta la atenuación deseada en dB,, a un valor de ganancia aritmética al cuadrado,, mediante el usoPor ejemplo, 3,010 dB se convierte en 0,5, 1 dB se convierte en 0,79432823 y así sucesivamente.
- Encontrar
- Encuentra las raíces de P(S) usando un algoritmo de búsqueda de raíces.
- 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..
Escalado de la función de transferencia
CuandoUna vez determinado, el polinomio de la función de transferencia de Hourglass se puede escalar de la siguiente manera:
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 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.
- ↑ 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.
- 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.
- 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.
- ^ 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.
- ↑ 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.
- ↑ 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.
- Filtros lineales
- Filtros de síntesis de red
- Diseño electrónico