Articulo de referencia

Funciones racionales elípticas

Gráfico de funciones racionales elípticas para x entre -1 y 1 para órdenes 1, 2, 3 y 4 con factor de discriminación ξ = 1,1. Todas están acotadas entre -1 y 1 y todas tienen el ...

Gráfico de funciones racionales elípticas para x entre -1 y 1 para órdenes 1, 2, 3 y 4 con factor de discriminación ξ = 1,1. Todas están acotadas entre -1 y 1 y todas tienen el valor 1 en x = 1 .

En matemáticas, las funciones racionales elípticas son una sucesión de funciones racionales con coeficientes reales. Estas funciones se utilizan ampliamente en el diseño de filtros electrónicos elípticos . (A veces se las denomina funciones racionales de Chebyshev , para no confundirlas con otras funciones del mismo nombre ).

Las funciones elípticas racionales se identifican mediante un orden entero positivo n e incluyen un parámetro ξ 1 llamado factor de selectividad . Una función elíptica racional de grado n en x con factor de selectividad ξ se define generalmente como:  

Rnorte(ξ,incógnita)dod(norteK(1/Lnorte(ξ))K(1/ξ)dod1(incógnita,1/ξ),1/Lnorte(ξ)){\displaystyle R_{n}(\xi ,x)\equiv \mathrm {cd} \left(n{\frac {K(1/L_{n}(\xi ))}{K(1/\xi )}}\,\mathrm {cd} ^{-1}(x,1/\xi ),1/L_{n}(\xi )\right)}

dónde

  • cd(u,k) es la función coseno elíptica de Jacobi .
  • K() es una integral elíptica completa de primera especie.
  • Lnorte(ξ)=Rnorte(ξ,ξ){\displaystyle L_{n}(\xi )=R_{n}(\xi ,\xi )}es el factor de discriminación , igual al valor mínimo de la magnitud deRnorte(ξ,incógnita){\displaystyle R_{n}(\xi ,x)}para|incógnita|ξ{\displaystyle |x|\geq \xi }.

En muchos casos, en particular para órdenes de la forma n  =  2 a 3 b, donde a y b son enteros, las funciones racionales elípticas pueden expresarse utilizando únicamente funciones algebraicas. Las funciones racionales elípticas están estrechamente relacionadas con los polinomios de Chebyshev : así como las funciones trigonométricas circulares son casos especiales de las funciones elípticas de Jacobi, los polinomios de Chebyshev son casos especiales de las funciones racionales elípticas.

Expresión como una razón de polinomios

Para órdenes pares, las funciones racionales elípticas pueden expresarse como una razón de dos polinomios, ambos de orden n .

Rnorte(ξ,incógnita)=r0i=1norte(incógnitaincógnitai)i=1norte(incógnitaincógnitapagi){\displaystyle R_{n}(\xi ,x)=r_{0}\,{\frac {\prod _{i=1}^{n}(x-x_{i})}{\prod _{i=1}^{n}(x-x_{pi})}}}    (para n par)

dóndeincógnitai{\displaystyle x_{i}}son los ceros yincógnitapagi{\displaystyle x_{pi}}son los polos, yr0{\displaystyle r_{0}}es una constante de normalización elegida de tal manera queRnorte(ξ,1)=1{\displaystyle R_{n}(\xi ,1)=1}La forma anterior sería válida también para órdenes pares, excepto que para órdenes impares habrá un polo en x=∞ y un cero en x=0, por lo que la forma anterior debe modificarse para que diga:

Rnorte(ξ,incógnita)=r0incógnitai=1norte1(incógnitaincógnitai)i=1norte1(incógnitaincógnitapagi){\displaystyle R_{n}(\xi ,x)=r_{0}\,x\,{\frac {\prod _{i=1}^{n-1}(x-x_{i})}{\prod _{i=1}^{n-1}(x-x_{pi})}}}    (para n impar)

Propiedades

Gráfico del valor absoluto de la función racional elíptica de tercer orden con ξ=1,4. Hay un cero en x=0 y un polo en el infinito. Como la función es antisimétrica, se observa que hay tres ceros y tres polos. Entre los ceros, la función aumenta hasta un valor de 1 y, entre los polos, la función disminuye hasta el valor del factor de discriminación L n
Gráfico del valor absoluto de la función racional elíptica de cuarto orden con ξ=1,4. Dado que la función es simétrica, se observa que hay cuatro ceros y cuatro polos. Entre los ceros, la función aumenta hasta un valor de 1 y, entre los polos, la función disminuye hasta el valor del factor de discriminación L n.
Gráfico del efecto del factor de selectividad ξ. Se muestra la función racional elíptica de cuarto orden con valores de ξ que varían desde casi la unidad hasta el infinito. La curva negra, correspondiente a ξ=∞, es el polinomio de Chebyshev de orden 4. Cuanto más cerca esté el factor de selectividad de la unidad, mayor será la pendiente en la región de transición entre x=1 y x=ξ.

Las propiedades canónicas

  • Rnorte2(ξ,incógnita)1{\displaystyle R_{n}^{2}(\xi ,x)\leq 1}para|incógnita|1{\displaystyle |x|\leq 1\,}
  • Rnorte2(ξ,incógnita)=1{\displaystyle R_{n}^{2}(\xi ,x)=1}en|incógnita|=1{\displaystyle |x|=1\,}
  • Rnorte2(ξ,incógnita)=Rnorte2(ξ,incógnita){\displaystyle R_{n}^{2}(\xi ,-x)=R_{n}^{2}(\xi ,x)}
  • Rnorte2(ξ,incógnita)>1{\displaystyle R_{n}^{2}(\xi ,x)>1}paraincógnita>1{\displaystyle x>1\,}
  • La pendiente en x=1 es lo más grande posible.
  • La pendiente en x=1 es mayor que la pendiente correspondiente del polinomio de Chebyshev del mismo orden.

La única función racional que satisface las propiedades anteriores es la función racional elíptica ( Lutovac, Tošić y Evans 2001 , § 13.2) . Se derivan las siguientes propiedades:

Normalización

La función racional elíptica se normaliza a la unidad en x=1:

Rnorte(ξ,1)=1{\displaystyle R_{n}(\xi ,1)=1\,}

Propiedad de anidación

La propiedad de anidamiento se escribe:

Rmetro(Rnorte(ξ,ξ),Rnorte(ξ,incógnita))=Rmetronorte(ξ,incógnita){\displaystyle R_{m}(R_{n}(\xi ,\xi ),R_{n}(\xi ,x))=R_{m\cdot n}(\xi ,x)\,}

Esta es una propiedad muy importante:

  • SiRnorte{\displaystyle R_{n}}es conocido para todos los primos n , entonces la propiedad de anidamiento daRnorte{\displaystyle R_{n}}para todo n . En particular, dado queR2{\displaystyle R_{2}}yR3{\displaystyle R_{3}}se puede expresar en forma cerrada sin el uso explícito de las funciones elípticas de Jacobi, entonces todoRnorte{\displaystyle R_{n}}para n de la formanorte=2a3b{\displaystyle n=2^{a}3^{b}}puede expresarse así.
  • De ello se deduce que si los ceros deRnorte{\displaystyle R_{n}}para los primos n son conocidos, los ceros de todosRnorte{\displaystyle R_{n}}Se pueden encontrar. Utilizando la relación de inversión (ver más abajo), también se pueden encontrar los polos.
  • La propiedad de anidamiento implica la propiedad de anidamiento del factor de discriminación:
Lmetronorte(ξ)=Lmetro(Lnorte(ξ)){\displaystyle L_{m\cdot n}(\xi )=L_{m}(L_{n}(\xi ))}

Valores límite

Las funciones racionales elípticas están relacionadas con los polinomios de Chebyshev de primera especie.Tnorte(incógnita){\displaystyle T_{n}(x)}por:

límiteξ=→Rnorte(ξ,incógnita)=Tnorte(incógnita){\displaystyle \lim _{\xi =\rightarrow \,\infty }R_{n}(\xi ,x)=T_{n}(x)\,}

Simetría

Rnorte(ξ,incógnita)=Rnorte(ξ,incógnita){\displaystyle R_{n}(\xi ,-x)=R_{n}(\xi ,x)\,}para n pares
Rnorte(ξ,incógnita)=Rnorte(ξ,incógnita){\displaystyle R_{n}(\xi ,-x)=-R_{n}(\xi ,x)\,}para n impar

Equiripple

Rnorte(ξ,incógnita){\displaystyle R_{n}(\xi ,x)}tiene igual ondulación de±1{\displaystyle \pm 1}en el intervalo1incógnita1{\displaystyle -1\leq x\leq 1}. Por la relación de inversión (véase más abajo), se deduce que1/Rnorte(ξ,incógnita){\displaystyle 1/R_{n}(\xi ,x)}tiene ondulación uniforme en1/ξincógnita1/ξ{\displaystyle -1/\xi \leq x\leq 1/\xi }de±1/Lnorte(ξ){\displaystyle \pm 1/L_{n}(\xi )}.

Relación de inversión

Se cumple la siguiente relación de inversión:

Rnorte(ξ,ξ/incógnita)=Rnorte(ξ,ξ)Rnorte(ξ,incógnita){\displaystyle R_{n}(\xi ,\xi /x)={\frac {R_{n}(\xi ,\xi )}{R_{n}(\xi ,x)}}\,}

Esto implica que los polos y los ceros vienen en pares, de tal manera que

incógnitapagiincógnitazi=ξ{\displaystyle x_{pi}x_{zi}=\xi \,}

Las funciones de orden impar tendrán un cero en x=0 y un polo correspondiente en el infinito.

Polos y ceros

Los ceros de la función racional elíptica de orden n se escribiránincógnitanortei(ξ){\displaystyle x_{ni}(\xi )}oincógnitanortei{\displaystyle x_{ni}}cuandoξ{\displaystyle \xi }Se conoce implícitamente. Los ceros de la función racional elíptica serán los ceros del polinomio en el numerador de la función.

La siguiente derivación de los ceros de la función racional elíptica es análoga a la de determinar los ceros de los polinomios de Chebyshev ( Lutovac, Tošić y Evans 2001 , § 12.6) . Utilizando el hecho de que para cualquier z

dod((2metro1)K(1/z),1z)=0{\displaystyle \mathrm {cd} \left((2m-1)K\left(1/z\right),{\frac {1}{z}}\right)=0\,}

La ecuación definitoria para las funciones racionales elípticas implica que

norteK(1/Lnorte)K(1/ξ)dod1(incógnitametro,1/ξ)=(2metro1)K(1/Lnorte){\displaystyle n{\frac {K(1/L_{n})}{K(1/\xi )}}\mathrm {cd} ^{-1}(x_{m},1/\xi )=(2m-1)K(1/L_{n})}

de modo que los ceros vienen dados por

incógnitametro=dod(K(1/ξ)2metro1norte,1ξ).{\displaystyle x_{m}=\mathrm {cd} \left(K(1/\xi )\,{\frac {2m-1}{n}},{\frac {1}{\xi }}\right).}

Utilizando la relación de inversión, se pueden calcular los polos.

De la propiedad de anidamiento, si los ceros deRmetro{\displaystyle R_{m}}yRnorte{\displaystyle R_{n}}se puede expresar algebraicamente (es decir, sin necesidad de calcular las funciones de la elipse de Jacobi) entonces los ceros deRmetronorte{\displaystyle R_{m\cdot n}}pueden expresarse algebraicamente. En particular, los ceros de las funciones racionales elípticas de orden2i3j{\displaystyle 2^{i}3^{j}}puede expresarse algebraicamente ( Lutovac, Tošić y Evans 2001 , § 12.9, 13.9) . Por ejemplo, podemos encontrar los ceros deR8(ξ,incógnita){\displaystyle R_{8}(\xi ,x)}De la siguiente manera: Definir

incógnitanorteRnorte(ξ,incógnita)LnorteRnorte(ξ,ξ)tnorte11/Lnorte2.{\displaystyle X_{n}\equiv R_{n}(\xi ,x)\qquad L_{n}\equiv R_{n}(\xi ,\xi )\qquad t_{n}\equiv {\sqrt {1-1/L_{n}^{2}}}.}

Luego, a partir de la propiedad de anidamiento y sabiendo que

R2(ξ,incógnita)=(t+1)incógnita21(t1)incógnita2+1{\displaystyle R_{2}(\xi ,x)={\frac {(t+1)x^{2}-1}{(t-1)x^{2}+1}}}

dóndet11/ξ2{\displaystyle t\equiv {\sqrt {1-1/\xi ^{2}}}}tenemos:

L2=1+t1t,L4=1+t21t2,L8=1+t41t4{\displaystyle L_{2}={\frac {1+t}{1-t}},\qquad L_{4}={\frac {1+t_{2}}{1-t_{2}}},\qquad L_{8}={\frac {1+t_{4}}{1-t_{4}}}}
incógnita2=(t+1)incógnita21(t1)incógnita2+1,incógnita4=(t2+1)incógnita221(t21)incógnita22+1,incógnita8=(t4+1)incógnita421(t41)incógnita42+1.{\displaystyle X_{2}={\frac {(t+1)x^{2}-1}{(t-1)x^{2}+1}},\qquad X_{4}={\frac {(t_{2}+1)X_{2}^{2}-1}{(t_{2}-1)X_{2}^{2}+1}},\qquad X_{8}={\frac {(t_{4}+1)X_{4}^{2}-1}{(t_{4}-1)X_{4}^{2}+1}}.}

Estas últimas tres ecuaciones pueden invertirse:

incógnita=1±1+t(1incógnita21+incógnita2),incógnita2=1±1+t2(1incógnita41+incógnita4),incógnita4=1±1+t4(1incógnita81+incógnita8).{\displaystyle x={\frac {1}{\pm {\sqrt {1+t\,\left({\frac {1-X_{2}}{1+X_{2}}}\right)}}}},\qquad X_{2}={\frac {1}{\pm {\sqrt {1+t_{2}\,\left({\frac {1-X_{4}}{1+X_{4}}}\right)}}}},\qquad X_{4}={\frac {1}{\pm {\sqrt {1+t_{4}\,\left({\frac {1-X_{8}}{1+X_{8}}}\right)}}}}.\qquad }

Para calcular los ceros deR8(ξ,incógnita){\displaystyle R_{8}(\xi ,x)}nosotros establecimosincógnita8=0{\displaystyle X_{8}=0}En la tercera ecuación, calcule los dos valores deincógnita4{\displaystyle X_{4}}, luego use estos valores deincógnita4{\displaystyle X_{4}}en la segunda ecuación para calcular cuatro valores deincógnita2{\displaystyle X_{2}}y finalmente, utilice estos valores en la primera ecuación para calcular los ocho ceros deR8(ξ,incógnita){\displaystyle R_{8}(\xi ,x)}. (Eltnorte{\displaystyle t_{n}}se calculan mediante una recursión similar.) Nuevamente, utilizando la relación de inversión, estos ceros se pueden usar para calcular los polos.

Valores particulares

Podemos escribir las primeras funciones racionales elípticas como:

R1(ξ,incógnita)=incógnita{\displaystyle R_{1}(\xi ,x)=x\,}
R2(ξ,incógnita)=(t+1)incógnita21(t1)incógnita2+1{\displaystyle R_{2}(\xi ,x)={\frac {(t+1)x^{2}-1}{(t-1)x^{2}+1}}}
dónde
t11ξ2{\displaystyle t\equiv {\sqrt {1-{\frac {1}{\xi ^{2}}}}}}
R3(ξ,incógnita)=incógnita(1incógnitapag2)(incógnita2incógnitaz2)(1incógnitaz2)(incógnita2incógnitapag2){\displaystyle R_{3}(\xi ,x)=x\,{\frac {(1-x_{p}^{2})(x^{2}-x_{z}^{2})}{(1-x_{z}^{2})(x^{2}-x_{p}^{2})}}}
dónde
GRAMO4ξ2+(4ξ2(ξ21))2/3{\displaystyle G\equiv {\sqrt {4\xi ^{2}+(4\xi ^{2}(\xi ^{2}\!-\!1))^{2/3}}}}
incógnitapag22ξ2GRAMO8ξ2(ξ2+1)+12GRAMOξ2GRAMO3GRAMO3{\displaystyle x_{p}^{2}\equiv {\frac {2\xi ^{2}{\sqrt {G}}}{{\sqrt {8\xi ^{2}(\xi ^{2}\!+\!1)+12G\xi ^{2}-G^{3}}}-{\sqrt {G^{3}}}}}}
incógnitaz2=ξ2/incógnitapag2{\displaystyle x_{z}^{2}=\xi ^{2}/x_{p}^{2}}
R4(ξ,incógnita)=R2(R2(ξ,ξ),R2(ξ,incógnita))=(1+t)(1+t)2incógnita42(1+t)(1+t)incógnita2+1(1+t)(1t)2incógnita42(1+t)(1t)incógnita2+1{\displaystyle R_{4}(\xi ,x)=R_{2}(R_{2}(\xi ,\xi ),R_{2}(\xi ,x))={\frac {(1+t)(1+{\sqrt {t}})^{2}x^{4}-2(1+t)(1+{\sqrt {t}})x^{2}+1}{(1+t)(1-{\sqrt {t}})^{2}x^{4}-2(1+t)(1-{\sqrt {t}})x^{2}+1}}}
R6(ξ,incógnita)=R3(R2(ξ,ξ),R2(ξ,incógnita)){\displaystyle R_{6}(\xi ,x)=R_{3}(R_{2}(\xi ,\xi ),R_{2}(\xi ,x))\,}etc.

Consulte Lutovac, Tošić y Evans (2001 , § 13) para obtener expresiones explícitas adicionales de orden n=5 ynorte=2i3j{\displaystyle n=2^{i}\,3^{j}}.

Los factores de discriminación correspondientes son:

L1(ξ)=ξ{\displaystyle L_{1}(\xi )=\xi \,}
L2(ξ)=1+t1t=(ξ+ξ21)2{\displaystyle L_{2}(\xi )={\frac {1+t}{1-t}}=\left(\xi +{\sqrt {\xi ^{2}-1}}\right)^{2}}
L3(ξ)=ξ3(1incógnitapag2ξ2incógnitapag2)2{\displaystyle L_{3}(\xi )=\xi ^{3}\left({\frac {1-x_{p}^{2}}{\xi ^{2}-x_{p}^{2}}}\right)^{2}}
L4(ξ)=(ξ+(ξ21)1/4)4(ξ+ξ21)2{\displaystyle L_{4}(\xi )=\left({\sqrt {\xi }}+(\xi ^{2}-1)^{1/4}\right)^{4}\left(\xi +{\sqrt {\xi ^{2}-1}}\right)^{2}}
L6(ξ)=L3(L2(ξ)){\displaystyle L_{6}(\xi )=L_{3}(L_{2}(\xi ))\,}etc.

Los ceros correspondientes sonincógnitanortej{\displaystyle x_{nj}}donde n es el orden y j es el número de ceros. Habrá un total de n ceros para cada orden.

incógnita11=0{\displaystyle x_{11}=0\,}
incógnita21=ξ1t{\displaystyle x_{21}=\xi {\sqrt {1-t}}\,}
incógnita22=incógnita21{\displaystyle x_{22}=-x_{21}\,}
incógnita31=incógnitaz{\displaystyle x_{31}=x_{z}\,}
incógnita32=0{\displaystyle x_{32}=0\,}
incógnita33=incógnita31{\displaystyle x_{33}=-x_{31}\,}
incógnita41=ξ(1t)(1+tt(t+1)){\displaystyle x_{41}=\xi {\sqrt {\left(1-{\sqrt {t}}\right)\left(1+t-{\sqrt {t(t+1)}}\right)}}\,}
incógnita42=ξ(1t)(1+t+t(t+1)){\displaystyle x_{42}=\xi {\sqrt {\left(1-{\sqrt {t}}\right)\left(1+t+{\sqrt {t(t+1)}}\right)}}\,}
incógnita43=incógnita42{\displaystyle x_{43}=-x_{42}\,}
incógnita44=incógnita41{\displaystyle x_{44}=-x_{41}\,}

A partir de la relación de inversión, los polos correspondientesincógnitapag,nortei{\displaystyle x_{p,ni}}puede ser encontrado porincógnitapag,nortei=ξ/(incógnitanortei){\displaystyle x_{p,ni}=\xi /(x_{ni})}

Referencias

  • Mundo matemático
  • 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.; Tošić, Dejan V.; Evans, Brian L. (2001). Diseño de Filtros para Procesamiento de Señales usando MATLAB© y Mathematica© . Nueva Jersey, Estados Unidos: Prentice Hall. ISBN 0-201-36130-2.