Articulo de referencia

nodos de Chebyshev

Los ceros de Chebyshev (puntos sólidos, líneas rojas) y los extremos (cuadrados huecos, líneas azules) son la proyección de dos conjuntos de puntos equidistantes en el círculo u...

Los ceros de Chebyshev (puntos sólidos, líneas rojas) y los extremos (cuadrados huecos, líneas azules) son la proyección de dos conjuntos de puntos equidistantes en el círculo unitario sobre el eje x . 2n puntos equidistantes en el círculo se proyectan sobre n ceros de Chebyshev o n+1 extremos de Chebyshev. (Aquí n = 5 ).
Los ceros de Chebyshev (puntos sólidos) son raíces de un polinomio de Chebyshev de primera especie (rojo). Los extremos de Chebyshev (cuadrados huecos) son raíces de un polinomio de Chebyshev de segunda especie (azul), y también los extremos (cruces) de un polinomio de Chebyshev de primera especie.

En análisis numérico , los nodos de Chebyshev (también llamados puntos de Chebyshev o cuadrícula de Chebyshev ) son un conjunto de números algebraicos específicos que se utilizan como nodos para la interpolación polinómica y la integración numérica . Son la proyección de un conjunto de puntos equidistantes en el círculo unitario sobre el intervalo real.[1,1]{\displaystyle [-1,1]}, el diámetro del círculo .

Hay dos tipos de nodos de Chebyshev. Elnorte{\displaystyle n}Los nodos de Chebyshev de primera especie , también llamados nodos de Chebyshev-Gauss [ 1 ] o ceros de Chebyshev , son los ceros de un polinomio de Chebyshev de primera especie ,Tnorte{\displaystyle T_{n}} . El correspondientenorte+1{\displaystyle n+1} Los nodos de Chebyshev de segundo tipo , también llamados nodos de Chebyshev-Lobatto [ 2 ] o extremos de Chebyshev , son los extremos deTnorte{\displaystyle T_{n}} , que también son los ceros de un polinomio de Chebyshev de segundo tipo,Unorte1{\displaystyle U_{n-1}} , junto con los dos extremos del intervalo. Ambos tipos de números se conocen comúnmente como nodos de Chebyshev o puntos de Chebyshev en la literatura. [ 3 ] Reciben su nombre del matemático ruso del siglo XIX Pafnuty Chebyshev , quien introdujo por primera vez los polinomios de Chebyshev.

A diferencia de otros nodos de interpolación, los nodos de Chebyshev se "anidan": los nodos existentes se conservan al duplicar el número de nodos, reduciendo a la mitad el cálculo para cada refinamiento de la malla. Los interpolantes polinomiales construidos a partir de nodos de Chebyshev minimizan el efecto del fenómeno de Runge . [ 4 ] Se pueden convertir fácilmente a una representación como una suma ponderada de polinomios de Chebyshev utilizando la transformada rápida de Fourier .

Definición

Nodos de Chebyshev de ambos tipos desdenorte=2{\displaystyle n=2}anorte=50{\displaystyle n=50}.

Para un entero positivo dadonorte{\displaystyle n}, elnorte{\displaystyle n}Los nodos de Chebyshev de primera especie están dados por

incógnitak=porque(k+12)πnorte,k=0,,norte1.{\displaystyle x_{k}=\cos {\frac {{\bigl (}k+{\tfrac {1}{2}}{\bigr )}\pi }{n}},\quad k=0,\ldots ,n-1.}

Esta es la proyección de2norte{\displaystyle 2n}Puntos equidistantes en el círculo unitario sobre el intervalo[1,1]{\displaystyle [-1,1]} , el diámetro del círculo. Estos puntos son también las raíces deTnorte{\displaystyle T_{n}} , el polinomio de Chebyshev de primera especie con gradonorte{\displaystyle n} .

Elnorte+1{\displaystyle n+1}Los nodos de Chebyshev de segundo tipo vienen dados por

incógnitak=porquekπnorte,k=0,,norte.{\displaystyle x_{k}=\cos {\frac {k\pi }{n}},\quad k=0,\ldots ,n.}

Esta es también la proyección de2norte{\displaystyle 2n} puntos equidistantes en el círculo unitario sobre[1,1]{\displaystyle [-1,1]} , esta vez incluyendo los extremos del intervalo, cada uno de los cuales es solo la proyección de un punto en el círculo en lugar de dos. Estos puntos son también los extremos deTnorte{\displaystyle T_{n}}en[1,1]{\displaystyle [-1,1]} , los lugares donde toma el valor±1{\displaystyle \pm 1} . [ 5 ] Los puntos interiores entre los nodos, sin incluir los extremos, son también los ceros deUnorte1{\displaystyle U_{n-1}} , un polinomio de Chebyshev de segundo tipo, un reescalamiento de la derivada deTnorte{\displaystyle T_{n}} .

Para nodos en un intervalo arbitrario[a,b]{\displaystyle [a,b]}una transformación afín de[1,1]{\displaystyle [-1,1]}se puede utilizar: incógnita~k=12(a+b)+12(ba)incógnitak.{\displaystyle {\tilde {x}}_{k}={\tfrac {1}{2}}(a+b)+{\tfrac {1}{2}}(ba)x_{k}.}

Propiedades

Ambos tipos de nodos son siempre simétricos con respecto a cero, el punto medio del intervalo.

Ejemplos

Los conjuntos de nodos para los primeros números enterosnorte{\displaystyle n}son: raíces(T0)={},raíces(U0)={},extremos(T1)={1,+1},raíces(T1)={0},raíces(U1)={0},extremos(T2)={1,0,+1},raíces(T2)={1/2,+1/2},raíces(U2)={1/2,+1/2},extremos(T3)={1,1/2,+1/2,+1}{\displaystyle {\begin{aligned}{\text{raíces}}(T_{0})&=\{\},&{\text{raíces}}(U_{0})&=\{\},&{\text{extremo}}(T_{1})&=\{-1,+1\},\\{\text{raíces}}(T_{1})&=\{0\},&{\text{raíces}}(U_{1})&=\{0\},&{\text{extremo}}(T_{2})&=\{-1,0,+1\},\\{\text{raíces}}(T_{2})&=\{-1/{\sqrt {2}},+1/{\sqrt {2}}\},&{\text{raíces}}(U_{2})&=\{-1/2,+1/2\},&{\text{extremos}}(T_{3})&=\{-1,-1/2,+1/2,+1\}\\\end{aligned}}}

Si bien estos conjuntos están ordenados por valores ascendentes, las fórmulas de definición dadas anteriormente generan los nodos de Chebyshev en orden inverso, del mayor al menor.

Aproximación

Los nodos de Chebyshev son importantes en la teoría de la aproximación porque forman un conjunto de nodos particularmente bueno para la interpolación polinómica . Dada una función f en el intervalo[1,+1]{\displaystyle [-1,+1]}ynorte{\displaystyle n}agujasincógnita1,incógnita2,,incógnitanorte,{\displaystyle x_{1},x_{2},\ldots ,x_{n},}En ese intervalo, el polinomio de interpolación es ese polinomio único.PAGnorte1{\displaystyle P_{n-1}}de grado como máximonorte1{\displaystyle n-1}que tiene valorF(incógnitai){\displaystyle f(x_{i})}en cada puntoincógnitai{\displaystyle x_{i}}. El error de interpolación enincógnita{\displaystyle x}es F(incógnita)PAGnorte1(incógnita)=F(norte)(ξ)norte¡i=1norte(incógnitaincógnitai){\displaystyle f(x)-P_{n-1}(x)={\frac {f^{(n)}(\xi )}{n!}}\prod _{i=1}^{n}(x-x_{i})} para algunosξ{\displaystyle \xi }(dependiendo de x ) en [ −1, 1 ] . [ 6 ] Por lo tanto, es lógico intentar minimizar máximoincógnita[1,1]|i=1norte(incógnitaincógnitai)|.{\displaystyle \max _{x\in [-1,1]}{\biggl |}\prod _{i=1}^{n}(x-x_{i}){\biggr |}.}

Este producto es un polinomio mónico de grado n . Se puede demostrar que el valor absoluto máximo (norma máxima) de cualquier polinomio de este tipo está acotado inferiormente por 2 1− n . Esta cota se alcanza mediante los polinomios de Chebyshev escalados 2 1− n T n , que también son mónicos. (Recordemos que | T n ( x )| ≤ 1 para x ∈ [−1, 1] . [ 7 ] ) Por lo tanto, cuando los nodos de interpolación x i son las raíces de T n , el error satisface |F(incógnita)PAGnorte1(incógnita)|12norte1norte¡máximoξ[1,1]|F(norte)(ξ)|.{\displaystyle \left|f(x)-P_{n-1}(x)\right|\leq {\frac {1}{2^{n-1}n!}}\max _{\xi \in [-1,1]}\left|f^{(n)}(\xi )\right|.} Para un intervalo arbitrario [ a , b ], un cambio de variable muestra que |F(incógnita)PAGnorte1(incógnita)|12norte1norte¡(ba2)nortemáximoξ[a,b]|F(norte)(ξ)|.{\displaystyle \left|f(x)-P_{n-1}(x)\right|\leq {\frac {1}{2^{n-1}n!}}\left({\frac {b-a}{2}}\right)^{n}\max _{\xi \in [a,b]}\left|f^{(n)}(\xi )\right|.}

Nodos de orden par modificados

Algunas aplicaciones para nodos de interpolación, como el diseño de filtros Chebyshev pasivos con terminación igual , no pueden usar directamente nodos Chebyshev de orden par debido a la falta de una raíz en 0. En cambio, los nodos Chebyshev se pueden mover hacia cero, con una raíz doble en cero directamente, usando una transformación: [ 8 ]

incógnita~k=sgn(incógnitak)incógnitak2incógnitanorte/221incógnitanorte/22{\displaystyle {\tilde {x}}_{k}=\operatorname {sgn} (x_{k}){\sqrt {\frac {x_{k}^{2}-x_{n/2}^{2}}{1-x_{n/2}^{2}}}}}

Por ejemplo, los nodos de Chebyshev de primera especie de orden 4 son0,9239,0,3827,0,3827,0,9239{\displaystyle {0.9239,0.3827,-0.3827,-0.9239}}, conincógnitanorte/2=0,382683{\displaystyle x_{n/2}=0.382683}Al aplicar la transformación se obtienen nuevos nodos.0,910180,0,0,0,910180{\displaystyle {0.910180,0,0,-0.910180}}Los nodos de orden par modificados ahora incluyen el cero dos veces.

Véase también

Notas

  1. El nombre nodos de Chebyshev-Gauss proviene del uso de ceros de Chebyshev en la integración numérica, que puede verse como una variante de la cuadratura gaussiana .
  2. El nombre de nodos de Chebyshev-Lobatto proviene de Rehuel Lobatto , quien creó una variante de la cuadratura gaussiana, conocida como cuadratura de Lobatto , cuyos nodos incluían los extremos del intervalo, una característica compartida por los extremos de Chebyshev.
  3. Trefethen 2013 , págs. 7 
  4. Fink y Mathews 1999 , págs. 236–238 
  5. Trefethen 2013 , pág. 14 
  6. Stewart 1996 , (20.3)
  7. Stewart 1996 , Lección 20, §14
  8. ^ Saal, Rudolf (enero de 1979). Manual de diseño de filtros (en inglés y alemán) (1ª  ed.). Múnich, Alemania: Allgemeine Elektricitäts-Gesellschaft. págs.25  , 26, 56-61 , 116, 117. ISBN 3-87087-070-2.

Referencias

  • Fink, Kurtis D.; Mathews, John H. (1999). Métodos numéricos con MATLAB (3.ª  ed.). Upper Saddle River, NJ: Prentice Hall.
  • Stewart, Gilbert W. (1996). Notas finales sobre análisis numérico . SIAM . ISBN 978-0-89871-362-6.
  • Trefethen, Lloyd N. (2013), Teoría de la aproximación y práctica de la aproximación , SIAM

Lecturas adicionales

  • Burden, Richard L.; Faires, J. Douglas: Análisis numérico , 8.ª ed., páginas 503–512, ISBN 0-534-39200-8.