Filtro "L" óptimo

Comparación de la magnitud del filtro entre los filtros tipo 1 de Butterworth, Legendre y Chebyshev

El filtro Optimum "L" (también conocido como filtro Legendre–Papoulis ) fue propuesto por Athanasios Papoulis en 1958. Tiene la tasa de caída máxima para un orden de filtro dado mientras mantiene una respuesta de frecuencia monótona . Proporciona un compromiso entre el filtro Butterworth , que es monótono pero tiene una caída más lenta, y el filtro Chebyshev , que tiene una caída más rápida pero tiene ondulación en la banda de paso o en la banda de rechazo . El diseño del filtro se basa en polinomios de Legendre , que es la razón de su nombre alternativo y de la "L" en Optimum "L".

Sintetizando los polinomios característicos

La solución a la síntesis polinomial característica del filtro L óptimo de orden N surge de la resolución del polinomio característico, , dadas las siguientes restricciones y definiciones. [1] L N ( ω 2 ) {\displaystyle L_{N}(\omega ^{2})}

L N ( 0 ) = 0 L ( 1 ) = 1 d L N ( ω 2 ) d ω  0 for 0  ω 1 d L N ( ω 2 ) d ω | ω = 1  is maximum {\displaystyle {\begin{aligned}&L_{N}(0)=0\\&L(1)=1\\&{dL_{N}(\omega ^{2}) \over d\omega }\geq {\text{ 0 for 0 }}\leq \omega \leq 1\\&{dL_{N}(\omega ^{2}) \over d\omega }{\biggr |}_{\omega =1}{\text{ is maximum}}\\\end{aligned}}}

El caso de orden impar [2] y el caso de orden par [1] pueden resolverse utilizando polinomios de Legendre de la siguiente manera.

N Odd: L N ( ω 2 ) = 2 N + 1 1 2 ω 2 1 ( i = 0 i = k a i P i ( x ) ) 2 d x Where P i ( x )  is the Legendre polynomial of the first kind of order i k = N 1 2 a i = 2 i + 1 2 ( k + 1 ) N Even: L N ( ω 2 ) = 1 2 ω 2 1 ( x + 1 ) ( i = 0 i = k a i P i ( x ) ) 2 d x Where k = N 2 2 a i = { 2 i + 1 ( k + 2 ) ( k + 1 ) , if  i  is odd and  k  is odd OR  i  is even and  k  is even 0 , otherwise {\displaystyle {\begin{aligned}&{\text{N Odd:}}\\&L_{N}(\omega ^{2})={\frac {2}{N+1}}\int _{-1}^{2\omega ^{2}-1}{\bigg (}\sum _{i=0}^{i=k}a_{i}P_{i}(x){\bigg )}^{2}dx\\&{\text{Where}}\\&P_{i}(x){\text{ is the Legendre polynomial of the first kind of order i}}\\&k={\frac {N-1}{2}}\\&a_{i}={\frac {2i+1}{\sqrt {2(k+1)}}}\\&\\&\\&{\text{N Even:}}\\&L_{N}(\omega ^{2})=\int _{-1}^{2\omega ^{2}-1}(x+1){\bigg (}\sum _{i=0}^{i=k}a_{i}P_{i}(x){\bigg )}^{2}dx\\&{\text{Where}}\\&k={\frac {N-2}{2}}\\&a_{i}={\begin{cases}{\frac {2i+1}{\sqrt {(k+2)(k+1)}}},&{\text{if }}i{\text{ is odd and }}k{\text{ is odd OR }}i{\text{ is even and }}k{\text{ is even}}\\0,&{\text{otherwise}}\\\end{cases}}\\\end{aligned}}}

Respuesta de frecuencia y función de transferencia

La magnitud de la frecuencia de magnitud se crea utilizando la siguiente fórmula. Dado que la función característica "L" óptima ya está en forma cuadrada, no se la debe elevar al cuadrado nuevamente como se hace para otros tipos de filtros, como los filtros Chebyshev y los filtros Butterworth .

T ( ω ) = 1 1 + ϵ 2 L N ( ω 2 ) ϵ 2 = 10 | δ | / 10 1. δ = magnitude attenuation of the passband in dB, usually 3.0103 {\displaystyle {\begin{aligned}&T(\omega )={\sqrt {\frac {1}{1+\epsilon ^{2}L_{N}(\omega ^{2})}}}\\&\epsilon ^{2}=10^{|\delta |/10}-1.\\&\delta ={\text{magnitude attenuation of the passband in dB, usually 3.0103}}\\\end{aligned}}}

Para obtener la función de transferencia, , haga que todos los coeficientes sean positivos para tener en cuenta el eje de frecuencia y luego use los polos del semiplano izquierdo para construir . Tenga en cuenta que es +1 para N par y -1 para N impar (consulte la tabla a continuación). El signo de debe factorizarse en las ecuaciones para las siguientes. [3] [4] T ( j ω ) {\displaystyle T(j\omega )} L N ( ω 2 ) {\displaystyle L_{N}(\omega ^{2})} j ω {\displaystyle j\omega } T ( j ω ) {\displaystyle T(j\omega )} L N ( ( j ω ) 2 ) {\displaystyle L_{N}((j\omega )^{2})} L N ( ω 2 ) {\displaystyle L_{N}(\omega ^{2})} L N ( ( j ω ) 2 ) {\displaystyle L_{N}((j\omega )^{2})} T ( j ω ) {\displaystyle T(j\omega )}

T ( j ω ) = 1 a + ϵ 2 L N ( ( j ω ) 2 ) | Left half plane Where: a = { 1 , if  N  is even 1 , if  N  is odd ϵ 2 = 10 | δ | / 10 1. δ = magnitude attenuation of the passband in dB, usually 3.010 {\displaystyle {\begin{aligned}&T(j\omega )={\sqrt {\frac {1}{a+\epsilon ^{2}L_{N}((j\omega )^{2})}}}{\bigg |}_{\text{Left half plane}}\\&{\text{Where:}}\\&a={\begin{cases}1,&{\text{if }}N{\text{ is even}}\\-1,&{\text{if }}N{\text{ is odd}}\end{cases}}\\&\epsilon ^{2}=10^{|\delta |/10}-1.\\&\delta ={\text{magnitude attenuation of the passband in dB, usually 3.010}}\\\end{aligned}}}

La restricción del "Semiplano izquierdo" se refiere a encontrar las raíces en todos los polinomios contenidos entre corchetes, seleccionar solo las raíces en el semiplano izquierdo y recrear los polinomios a partir de esas raíces.

Ejemplo: Función de transferencia de cuarto orden

N = 4 (cuarto orden), atenuación de banda de paso = -3,010 a 1 r/s.

Un filtro de cuarto orden tiene un valor para k de 1, que es impar, por lo que la suma utiliza solo valores impares de i para y , que incluye solo el término i = 1 en la suma. a i {\displaystyle a_{i}} P i ( x ) {\displaystyle P_{i}(x)}

La función de transferencia, , puede derivarse de la siguiente manera: T 4 ( j ω ) {\displaystyle T_{4}(j\omega )}

k = N 2 2 = 1  ( k  is odd) a 1 = 2 ( 1 ) + 1 ( ( 1 ) + 2 ) ( ( 1 ) + 1 ) = 1.2247449 P 1 ( x ) = x ( x + 1 ) ( i = 0 i = k a i P i ( x ) ) 2 = ( x + 1 ) ( 1.2247449 ( x ) ) 2 = 1.5 x 3 + 1.5 x 2 L 4 ( x 2 ) = 1 2 x 2 1 1.5 x 3 + 1.5 x 2   d x = 6 x 8 8 x 6 + 3 x 4 L 4 ( x 2 ) = 6 x 8 8 x 6 + 3 x 4 L 4 ( j ω 2 ) = 6 ( j ω ) 8 + 8 ( j ω ) 6 + 3 ( j ω ) 4 e c h o = ϵ 3.0103 / 10 1 = 1 T 4 ( j ω ) = [ 1 1 + 1 2 ( 6 ( j ω ) 8 + 8 ( j ω ) 6 + 3 ( j ω ) 4 ) ] Left Half Plane T 4 ( j ω ) = 1 2.4494897 ( j ω ) 4 + 3.8282201 ( j ω ) 3 + 4.6244874 ( j ω ) 2 + 3.0412127 ( j ω ) + 1 {\displaystyle {\begin{aligned}&k={\frac {N-2}{2}}=1{\text{ (}}k{\text{ is odd)}}\\&a_{1}={\frac {2(1)+1}{\sqrt {((1)+2)((1)+1)}}}=1.2247449\\&P_{1}(x)=x\\&(x+1){\bigg (}\sum _{i=0}^{i=k}a_{i}P_{i}(x){\bigg )}^{2}=(x+1){\bigr (}1.2247449(x){\bigr )}^{2}=1.5x^{3}+1.5x^{2}\\&L_{4}(x^{2})=\int _{-1}^{2x^{2}-1}1.5x^{3}+1.5x^{2}{\text{ }}dx=6x^{8}-8x^{6}+3x^{4}\\&L_{4}(x^{2})=6x^{8}-8x^{6}+3x^{4}\\&L_{4}(j\omega ^{2})=6(j\omega )^{8}+8(j\omega )^{6}+3(j\omega )^{4}\\&echo={\sqrt {\epsilon ^{3.0103/10}-1}}=1\\&T_{4}(j\omega )={\bigg [}{\frac {1}{1+1^{2}(6(j\omega )^{8}+8(j\omega )^{6}+3(j\omega )^{4})}}{\bigg ]}_{\text{Left Half Plane}}\\&\\&T_{4}(j\omega )={\frac {1}{2.4494897(j\omega )^{4}+3.8282201(j\omega )^{3}+4.6244874(j\omega )^{2}+3.0412127(j\omega )+1}}\end{aligned}}}

Una rápida comprobación de cordura calcula un valor de -3,0103 dB, que es lo esperado. T 4 ( j ) {\displaystyle T_{4}(j)}

Tabla de los primeros 10 polinomios característicos

norte L N ( ω 2 ) {\displaystyle L_{N}(\omega ^{2})}
1 ω 2 {\textstyle \omega ^{2}}
2 ω 4 {\textstyle \omega ^{4}}
3 3 ω 6 3 ω 4 + ω 2 {\textstyle 3\omega ^{6}-3\omega ^{4}+\omega ^{2}}
4 6 ω 8 8 ω 6 + 3 ω 4 {\textstyle 6\omega ^{8}-8\omega ^{6}+3\omega ^{4}}
5 20 ω 10 40 ω 8 + 28 ω 6 8 ω 4 + ω 2 {\textstyle 20\omega ^{10}-40\omega ^{8}+28\omega ^{6}-8\omega ^{4}+\omega ^{2}}
6 50 ω 12 120 ω 10 + 105 ω 8 40 ω 6 + 6 ω 4 {\textstyle 50\omega ^{12}-120\omega ^{10}+105\omega ^{8}-40\omega ^{6}+6\omega ^{4}}
7 175 ω 14 525 ω 12 + 615 ω 10 355 ω 8 + 105 ω 6 15 ω 4 + ω 2 {\textstyle 175\omega ^{14}-525\omega ^{12}+615\omega ^{10}-355\omega ^{8}+105\omega ^{6}-15\omega ^{4}+\omega ^{2}}
8 490 ω 16 1668 ω 14 + 2310 ω 12 1624 ω 10 + 615 ω 8 120 ω 6 + 10 ω 4 {\textstyle 490\omega ^{16}-1668\omega ^{14}+2310\omega ^{12}-1624\omega ^{10}+615\omega ^{8}-120\omega ^{6}+10\omega ^{4}}
9 1764 ω 18 7056 ω 16 + 11704 ω 14 10416 ω 12 + 5376 ω 10 1624 ω 8 + 276 ω 6 24 ω 4 + ω 2 {\textstyle 1764\omega ^{18}-7056\omega ^{16}+11704\omega ^{14}-10416\omega ^{12}+5376\omega ^{10}-1624\omega ^{8}+276\omega ^{6}-24\omega ^{4}+\omega ^{2}}
10 5292 ω 20 23520 ω 18 + 44100 ω 16 45360 ω 14 + 27860 ω 12 10416 ω 10 + 2310 ω 8 280 ω 6 + 15 ω 4 {\textstyle 5292\omega ^{20}-23520\omega ^{18}+44100\omega ^{16}-45360\omega ^{14}+27860\omega ^{12}-10416\omega ^{10}+2310\omega ^{8}-280\omega ^{6}+15\omega ^{4}}

La tabla se calcula a partir de las ecuaciones anteriores para L N ( ω 2 ) {\displaystyle L_{N}(\omega ^{2})}

Véase también

Referencias

  1. ^ ab Fukada, Minoru (septiembre de 1959). "Filtros óptimos de órdenes pares con respuesta monótona". IRE Transactions on Circuit Theory . 6 (3): 277–281. doi :10.1109/TCT.1959.1086558 – vía IEEE Xplore.
  2. ^ Papoulis, Athanasios (marzo de 1958). "Filtros óptimos con respuesta monótona". Actas del IRE . 46 (3): 606–609. doi :10.1109/JRPROC.1958.286876 – vía IEEE Xplore.
  3. ^ Notas de la conferencia sobre diseño de filtros del Dr. Byron Bennett, 1985, Universidad Estatal de Montana, Departamento de Ingeniería Eléctrica, Bozeman , Montana, EE. UU.
  4. ^ Sedra, Adel S.; Brackett, Peter O. (1978). Teoría y diseño de filtros: activos y pasivos. Beaverton, Oegon, EE. UU.: Matrix Publishers, Inc., págs. 45-73. ISBN 978-0916460143.{{cite book}}: CS1 maint: date and year (link)
  • Papoulis, Athanasios (1958). "Filtros óptimos con respuesta monótona". Proc. IRE . 46 (marzo): 606–609. doi :10.1109/JRPROC.1958.286876.
  • Kuo, Franklin F. (1966). Análisis y síntesis de redes . Wiley. ISBN 0-471-51118-8.Segunda edición.
  • Filtros “L” óptimos: polinomios, polos y elementos de circuitos, de C. Bond, 2004
  • Notas sobre filtros “L” (óptimos) por C. Bond, 2011


Retrieved from "https://en.wikipedia.org/w/index.php?title=Optimum_%22L%22_filter&oldid=1223005113"