Nodes de Txebixev

En anàlisi numèrica, els nodes de Txebixev són una distribució de nodes que permeten fer una interpolació més estable numèricament. Entre altres avantatges, aquests nodes permeten interpolar amb error proper a la màquina funcions que pateixen el fenomen de Runge.

Introducció

Els nodes equiespaiats tenen desavantatges, entre d'altres que en els extrems de l'interval d'interpolació l'error és molt acusat. Per això cal utilitzar distribucions de nodes més condensades en els extrems, que permeten solucionar aquest problema. La manera de trobar aquests nodes està estretament relacionada amb els polinomis de Txebixev de primer i segon tipus.

Polinomis de Txebixev

Definim el polinomi de Txebixev de grau n de primer tipus com:

$\, T_n(x) = \cos(n \arccos(x))$

Per altra banda definim el polinomi de Txebixev de grau n de segon tipus com:

$U_n(x) = \frac{\sin[(n+1) \arccos(x)]}{\sin(\arccos x)}$

Definició

Hi ha dues classes de nodes de Txebixev diferents, però que es defineixen de forma similar.

• El primer cas és el següent:

Definim el conjunt de n+1 nodes d'interpolació de Txebixev com les abscisses de les arrels del polinomi n+1 de Txebixev de primer tipus. És a dir: $x_j = -\cos \frac{(2j+1)\pi}{2(n+1)} \quad j = 0, 1, \ldots, n$ Més concretament, aquests són els anomenats nodes de Gauss-Txebixev.

Demostració

Tenim el polinomi de Txebixev de primer tipus $\, T_{n+1}(x) = \cos[(n+1)\arccos x]$ i l'igualem a zero per trobar les arrels: $\cos[(n+1)\arccos x] = 0 \quad \Rightarrow \quad (n+1)\arccos x = (2j+1)\pi, \quad j = 0, 1, \ldots, n$ Acabant d'aïllar aquesta expressió: $x = \cos \frac{(2j+1)\pi}{2(n+1)}, \quad j = 0, 1, \ldots, n$ Com que volem que el primer node sigui el més proper a -1 i l'últim el més proper a 1, afegim un signe menys i obtenim: $x_j = -\cos \frac{(2j+1)\pi}{2(n+1)}, \quad j = 0, 1, \ldots, n$

• El segon cas és:

Definim el conjunt de n+1 nodes d'interpolació de Txebixev com les abscisses dels zeros del polinomi de Txebixev de segon tipus de grau n-1. És a dir: $x_j = - \cos\left(\frac{\pi j}{n}\right) \quad j = 0,1,\ldots, n$ Més concretament, aquests són els anomenats nodes de Gauss-Lobatto-Txebixev.

Demostració

Tenim el n-1 polinomi de Txebixev de segon tipus $\, U_{n-1}(x) = \frac{\sin(n\arccos x)}{\sin(\arccos x)}$ i l'igualem a zero per trobar les arrels: $\frac{\sin(n\arccos x)}{\sin(\arccos x)} = 0 \quad \Rightarrow \quad n \arccos x = \pi j, \quad j = 0, 1, \ldots, n$ Acabant d'aïllar aquesta expressió: $x = \cos\left(\frac{\pi j}{n}\right) \quad j = 0,1,\ldots, n$ Com que volem que el primer node sigui igual a -1 i l'últim a 1, afegim un signe menys i obtenim: $x_j = -\cos\left(\frac{\pi j}{n}\right) \quad j = 0,1,\ldots, n$

S'ha comentat al principi que ambdós conjunts de nodes estan relacionats. Fixem-nos, doncs, què passa si definim els nodes de Txebixeb de primer tipus com:

$x_j = -\cos\left(\frac{\pi j}{n}\right), \quad j = \frac{1}{2}, \, \ldots \, , \, \frac{1}{2}+i \, , \, \ldots \, , n - \frac{1}{2}.$

La qual cosa queda com:

$x_i = -\cos\left(\frac{\pi(i + \frac{1}{2})}{n}\right) = -\cos\left(\frac{\pi(2i + 1)}{2n}\right), \quad i = 0, \ldots, n-1.$

Que de fet, és el conjunt de n (en lloc de n+1) nodes de primera classe tal i com els hem definit abans.

Vegeu també

• Polinomis de Txebixev
• Constant de Lebesgue
• Interpolació polinòmica
• Fenomen de Runge