Polynôme de Tchebychev
# Polynôme de Tchebychev

En mathématiques, un polynôme de Tchebychev est un terme de l'une des deux suites de polynômes orthogonaux particulières reliées à la formule de Moivre. Les polynômes de Tchebychev sont nommés ainsi en l'honneur du mathématicien russe Pafnouti Lvovitch Tchebychev.

Il existe deux suites de polynômes de Tchebychev, l'une nommée polynômes de Tchebychev de première espèce et notée $T_n$ et l'autre nommée polynômes de Tchebychev de seconde espèce et notée $U_n$ (dans les deux cas, l'entier naturel n correspond au degré).

Ces deux suites peuvent être définies par la relation de récurrence :

    $\forall n\in \mathbb {N} \quad P_{n+2}=2X~P_{n+1}-P_{n}$

et les deux premiers termes :

    $T_{0}=1,~T_{1}=X{\text{ pour la suite }}T$ et
    $U_{0}=1,~U_{1}=2X{\text{ pour la suite }}U$.

Chacune est une suite de polynômes orthogonaux par rapport à un produit scalaire de fonctions, associé à la fonction poids $w(x)=\left(1-x^{2}\right)^{-1/2}$ sur [–1, 1]. Ces polynômes constituent un cas particulier des polynômes ultrasphériques.

Une définition alternative de ces polynômes peut être donnée par les relations trigonométriques :

    $T_{n}\left(\cos \theta \right)=\cos \left(n\theta \right),\quad {\text{soit :}}\quad T_{n}(x)=\cos \left(n\arccos x\right)\quad {\text{et}}\quad U_{n}\left(\cos \theta \right)={\frac {\sin \left((n+1)\theta \right)}{\sin \theta }}$,

ce qui revient, par exemple, à considérer $T_n(cos$ θ) comme le développement de cos(nθ) sous forme de polynôme en cos θ.

Contrairement à d'autres familles de polynômes orthogonaux, tels ceux de Legendre, d'Hermite ou de Laguerre, les polynômes de Tchebychev n'ont pratiquement pas d'application directe en physique. En revanche, ils sont particulièrement utiles en analyse numérique pour l'interpolation polynomiale de fonctions. En premier lieu, en ce qui concerne le choix des points d'interpolation, comme les zéros de $T_n(x)$ ou abscisses de Tchebychev, en vue de limiter le phénomène de Runge. Également, ils constituent une base alternative de polynômes par rapport à la base canonique $X^n$ de $\mathbb {R} [X]$ des polynômes de Lagrange, ce qui permet d'améliorer sensiblement la convergence. Ils sont notamment utilisés pour le calcul des éphémérides astrononomiques

## Polynômes de Tchebychev de première espèce

Il existe plusieurs possibilités pour définir cette famille de polynômes. La plus simple est par la relation de récurrence, qui permet de générer rapidement l'expression des différents polynômes. Toutefois, une telle définition ne permet guère d'établir les propriétés générales de ces polynômes, en premier lieu leur orthogonalité, aussi une autre définition, à partir des propriétés des fonctions trigonométriques, doit être envisagée.

### Définition par la relation de récurrence

La définition classique des polynômes de Tchebychev de première espèce est le plus souvent donnée par la relation de récurrence suivante :

    $T_{n+1}=2XT_{n}-T_{n-1},\ \ \forall n\geq 1,\ {\textrm {avec}}\ T_{0}=1\ {\textrm {et}}\ T_{1}=X$.

Par récurrence, $T_n$ est un polynôme de degré n.

Les premiers polynômes de Tchebychev de première espèce sont :

    $T_{0}=1$
    $T_{1}=X$
    $T_{2}=2X^{2}-1$
    $T_{3}=4X^{3}-3X$
    $T_{4}=8X^{4}-8X^{2}+1$
    $T_{5}=16X^{5}-20X^{3}+5X$
    $T_{6}=32X^{6}-48X^{4}+18X^{2}-1$
    $T_{7}=64X^{7}-112X^{5}+56X^{3}-7X$
    $T_{8}=128X^{8}-256X^{6}+160X^{4}-32X^{2}+1$
    $T_{9}=256X^{9}-576X^{7}+432X^{5}-120X^{3}+9X$.

### Définition trigonométrique

On démontre que pour tout entier naturel n,

    $\forall \theta \in \mathbb {R} \quad \cos(n\theta )=T_{n}(\cos \theta )$,

ce qui peut servir de définition alternative des polynômes $T_n,$ vus comme fonctions polynomiales définies sur l'intervalle réel [–1, 1].

L'une des démonstrations se fait par récurrence d'ordre 2, à l'aide de l'identité trigonométrique de Simpson suivante :

    $\cos(a+b)+\cos(a-b)=2\cos a\cos b$.

Le caractère orthogonal des polynômes $T_n$ découle alors directement de celui des fonctions cos(nθ). Plus précisément, cette formule de Simpson montre de plus que les polynômes $T_n$ sont orthogonaux par rapport à la fonction poids $w(x)=\left(1-x^{2}\right)^{-{\tfrac {1}{2}}}$. En effet, pour deux entiers naturels n et p et avec le changement de variable x = cos θ, il vient

    $\int _{0}^{\pi }{\cos(n\theta )\cos(p\theta )\,\mathrm {d} \theta }=\int _{-1}^{+1}{T_{n}(x)T_{p}(x)\left(1-x^{2}\right)^{-1/2}\,\mathrm {d} x}$

Puis, à l'aide de la formule de Simpson :

    $\int _{-1}^{+1}{T_{n}(x)T_{p}(x)\left(1-x^{2}\right)^{-1/2}\,\mathrm {d} x}=\int _{0}^{\pi }{{\tfrac {\cos((n+p)\theta )+\cos((n-p)\theta )}{2}}\,\mathrm {d} \theta }={\begin{cases}0&{\text{si }}n\neq p\\\pi &{\text{si }}n=p=0\\\pi /2&{\text{si }}n=p\neq 0\end{cases}}$.

### Équation différentielle

Pour tout n, la fonction $\theta \mapsto \cos(n\theta )=T_{n}(\cos \theta )$ est solution de l'équation différentielle linéaire homogène d'ordre 2 à coefficients constants :

    $y''+n^{2}y=0$.

Par suite, les polynômes de Tchebychev sont solutions de l'équation différentielle formelle :

    $\left(1-X^{2}\right)T''_{n}-XT'_{n}+n^{2}T_{n}=0$.

Celle-ci peut aussi se mettre sous la forme d'une équation différentielle de Sturm-Liouville :

    ${\frac {\mathrm {d} }{\mathrm {d} x}}\left(\left(1-x^{2}\right)^{\tfrac {1}{2}}{\frac {\mathrm {d} T_{n}}{\mathrm {d} x}}\right)+n^{2}\left(1-x^{2}\right)^{-{\tfrac {1}{2}}}T_{n}=0$.

### Autres propriétés

* Pour tout entier naturel n,$T_{n}=\sum _{k=0}^{\left\lfloor {\frac {n}{2}}\right\rfloor }{\binom {n}{2k}}(X^{2}-1)^{k}X^{n-2k}$où ${\binom {n}{2k}}$ est un coefficient binomial et $\left\lfloor \cdot \right\rfloor$ note la partie entière.
* Pour tout entier n strictement positif,

$T_{n}={\frac {n}{2}}\sum _{k=0}^{\left\lfloor {\frac {n}{2}}\right\rfloor }(-1)^{k}{\frac {(n-k-1)!}{k!(n-2k)!}}(2X)^{n-2k}$.

* Pour tout entier naturel n,$T_{n}(1)=1$.

Cette propriété se démontre aisément en considérant la forme trigonométrique de $T_n,$ le cas x = 1 correspondant à θ = 0. On a aussi $T_{n}(-1)=(-1)^{n}$, qui découle de la symétrie $T_{n}(-x)=(-1)^{n}T_{n}(x)$.

* Quels que soient les entiers naturels m et n,$T_{n}\circ T_{m}=T_{mn}$ et $2T_{m}T_{n}=T_{m+n}+T_{|m-n|}$.
* Pour tout entier n strictement positif, le coefficient dominant de $T_n$ est $2^{n–1}$ et ses n racines sont

$a_{k}^{(n)}=\cos \left({\frac {(2k-1)\pi }{2n}}\right),\quad \forall k\in \{1,\ldots ,n\}$.

* Pour tout entier n > 0, les extremums de $T_n$ sur l'intervalle [–1, 1] sont atteints en$e_{k}^{(n)}=\cos {\frac {k\pi }{n}},\quad \forall k\in \{0,\ldots ,n\}$

(ce sont –1, 1 et les racines de $U_{n–1}),$ et $T_{n}(e_{k}^{(n)})=(-1)^{k}$.

* La parité dépend de n : $T_{n}(-X)=(-1)^{n}T_{n}(X)$.
* Représentation intégrale :$T_{n}(x)={\frac {1}{4i\pi }}\int _{C}{\frac {1}{z^{n}}}{\frac {1-z^{2}}{z(1-2xz+z^{2})}}\,\mathrm {d} z$ où C est un contour du plan complexe parcouru dans le sens trigonométrique, contenant zéro et excluant les zéros de 1 – 2xz + $z^2.$
* Séries génératrices
   * ordinaire : $\sum _{n=0}^{\infty }T_{n}t^{n}={\frac {1-tX}{1-2tX+t^{2}}}$,
   * exponentielle : $\sum _{n=0}^{\infty }T_{n}{\frac {t^{n}}{n!}}={\frac {1}{2}}\left({\rm {e}}^{(X-{\sqrt {X^{2}-1}})t}+{\rm {e}}^{(X+{\sqrt {X^{2}-1}})t}\right),$
   * $\sum \limits _{n=1}^{\infty }T_{n}(X){\frac {t^{n}}{n}}=-{\frac {1}{2}}\ln \left(1-2tX+t^{2}\right),$ pertinente en particulier en théorie du potentiel.

## Polynômes de Tchebychev de seconde espèce

### Définition par récurrence

Les polynômes de seconde espèce $U_n$ peuvent se définir par la même relation de récurrence que ceux de première espèce, avec des premiers termes différents :

    $U_{n+1}=2XU_{n}-U_{n-1},\ \ \forall n\geq 1\ {\textrm {avec}}\ U_{0}=1\ {\textrm {et}}\ U_{1}=2X$.

Par récurrence, $U_n$ est un polynôme de degré n.

Les premiers polynômes de Tchebychev de seconde espèce sont :

    $U_{0}=1$
    $U_{1}=2X$
    $U_{2}=4X^{2}-1$
    $U_{3}=8X^{3}-4X$
    $U_{4}=16X^{4}-12X^{2}+1$
    $U_{5}=32X^{5}-32X^{3}+6X$
    Échec de l’analyse (SVG (MathML peut être activé via une extension du navigateur) : réponse non valide(« Math extension cannot connect to Restbase. ») du serveur « http://localhost:6011/fr.wikipedia.org/v1/ » :): {\displaystyle U_6 = 64X^6 - 80X^4 + 24X^2 - 1}
    $U_{7}=128X^{7}-192X^{5}+80X^{3}-8X$
    $U_{8}=256X^{8}-448X^{6}+240X^{4}-40X^{2}+1$
    $U_{9}=512X^{9}-1024X^{7}+672X^{5}-160X^{3}+10X$.

### Définition trigonométrique

De la même façon que pour ceux de première espèce, les polynômes $U_n$ peuvent se définir alternativement par la forme trigonométrique de leur fonction polynomiale associée sur ]–1 ; 1[. On montre en effet que pour tout n :

    $\forall \theta \in \mathbb {R} \setminus \pi \mathbb {Z} \quad {\frac {\sin \left((n+1)\theta \right)}{\sin \theta }}=U_{n}(\cos \theta )$.

Là encore, le caractère orthogonal des polynômes $U_n$ découle directement de celui des fonctions $\sin n\theta$. Plus précisément, comme $2\sin a\sin b=\cos(a-b)-\cos(a+b)$ :

    $\int _{0}^{\pi }{U_{p}(\cos \theta )U_{n}(\cos \theta )\sin ^{2}(\theta )\,\mathrm {d} \theta }=\int _{0}^{\pi }{\sin \left((p+1)\theta \right)\sin \left((n+1)\theta \right)\,\mathrm {d} \theta }=0{\text{ si }}p\neq n$.

En exprimant cette intégrale en fonction de la variable x = cos θ, on en déduit que les polynômes $U_n$ sont orthogonaux par rapport à la fonction poids $w(x)=\left(1-x^{2}\right)^{1/2}$ :

    $\int _{-1}^{+1}{U_{n}(x)U_{p}(x)\left(1-x^{2}\right)^{1/2}\,\mathrm {d} x}=0{\text{ si }}p\neq n$.

### Équation différentielle

Pour tout n, la fonction $\theta \mapsto \sin(n\theta )=U_{n-1}(\cos \theta )\sin \theta$ est solution de l'équation différentielle linéaire homogène d'ordre 2 à coefficients constants :

    $y''+n^{2}y=0$.

Par suite, les polynômes de Tchebychev de seconde espèce sont solutions de l'équation différentielle formelle :

    $(1-X^{2})U_{n}''-3XU_{n}'+n(n+2)U_{n}=0$.

### Autres propriétés

* Pour tout entier naturel n,$U_{n}=\sum _{k=0}^{\left\lfloor {\frac {n}{2}}\right\rfloor }{\binom {n+1}{2k+1}}(X^{2}-1)^{k}X^{n-2k}=X^{n}\sum _{k=0}^{\left\lfloor {\frac {n}{2}}\right\rfloor }{\binom {n+1}{2k+1}}(1-X^{-2})^{k}=\sum _{k=0}^{\left\lfloor {\frac {n}{2}}\right\rfloor }(-1)^{k}{\binom {n-k}{k}}(2X)^{n-2k}$.
* Les $U_n$ sont orthogonaux pour le produit scalaire associé à la pondération ${\sqrt {1-x^{2}}}$ sur l'intervalle [–1 ; 1]. Plus précisément :

    $\int _{-1}^{1}U_{n}(x)U_{m}(x){\sqrt {1-x^{2}}}\,\mathrm {d} x={\begin{cases}0&{\text{si }}n\neq m\\\pi /2&{\text{si }}n=m.\end{cases}}$
* Pour tout entier naturel n,

    $U_{n}(1)=n+1$.
* Si $m\geq n$, $U_{m}U_{n}=\sum _{k=0}^{n}U_{m-n+2k}$.
* Pour tout entier n strictement positif, les n racines de $U_n$ sont

    $e_{k}^{(n+1)}=\cos {\frac {k\pi }{n+1}},\quad \forall k\in \{1,\ldots ,n\}$.
* La parité dépend de n :

    $U_{n}(-X)=(-1)^{n}U_{n}(X)$.
* Représentation intégrale :

    $U_{n}(x)={\frac {1}{2i\pi }}\int _{C}{\frac {1}{z^{n}}}{\frac {1}{z(1-2xz+z^{2})}}\,\mathrm {d} z$ où C est un contour du plan complexe parcouru dans le sens trigonométrique, contenant zéro et excluant les zéros de 1 – 2xz + $z^2.$
* Série génératrice :

    $\sum _{n=0}^{\infty }U_{n}(X)t^{n}={\frac {1}{1-2tX+t^{2}}}$.

## Quelques relations avec d'autres fonctions spéciales

* $T_{n}=U_{n}-XU_{n-1},\quad T_{n+1}=XT_{n}-(1-X^{2})U_{n-1}{\text{ et }}T_{n}'=nU_{n-1}$,
* $T_{n}={\frac {n}{2}}C_{n}^{(0)}{\text{ et }}U_{n}=C_{n}^{(1)}$ où les $C^{(k)}$ $_n$ sont les polynômes de Gegenbauer et
* $T_{n}(x)=F\left(-n,n;{\frac {1}{2}};{\frac {1-x}{2}}\right){\text{ et }}U_{n}(x)=(n+1)F\left(-n,n+2;{\frac {3}{2}};{\frac {1-x}{2}}\right)$ où F est la fonction hypergéométrique.

## Historique

Tchebychev a découvert ces familles en travaillant sur le problème de convergence des interpolations de Lagrange. On peut démontrer qu'en choisissant les racines des polynômes de Tchebychev comme points d'interpolation, on minimise les écarts (cf. phénomène de Runge). Dans ce contexte, les a(n) k indiqués ci-dessus, éventuellement ajustés à un autre intervalle d'interpolation [a, b] (par une transformation affine x ↦ b – a/2x + b + a/2), sont appelés les abscisses de Tchebychev.

En effet, on peut montrer que l'erreur entre la fonction interpolée et le polynôme d'interpolation aux points $x_0,...,x_n$ sur [a, b] s'exprime en

    $\max _{x\in [a,b]}\left|(x-x_{0})\dots (x-x_{n})\right|{\frac {\|f^{(n+1)}\|_{\infty ,[a,b]}}{(n+1)!}}$.

L'idée fut donc de minimiser $L=\max _{x\in [a,b]}\left|(x-x_{0})\dots (x-x_{n})\right|$ pour n points donnés. Tchebychev montra que dans le cas où l'intervalle est [–1, 1] et la répartition des points est symétrique, le polynôme optimal prend les valeurs –L et +L alternativement et n + 1 fois exactement (on dit que le polynôme présente une alternance de Tchebychev). C'est cette propriété qui permet de déduire que les abscisses de Tchebychev sont les meilleurs points d'interpolation pour minimiser les oscillations du polynôme d'interpolation et donc obtenir la meilleure convergence possible.

## Applications

Les polynômes de Tchebychev permettent de démontrer le théorème de Weierstrass selon lequel toute fonction continue sur un segment est limite uniforme d'une suite de fonctions polynomiales.

Ils sont également impliqués dans le calcul de filtres en électronique analogique, les filtres de Tchebychev.

Enfin, ils permettent une explication théorique de l'efficacité supérieure de la transformée en cosinus discrète dans le cadre de l'interpolation d'un signal numérique échantillonné, par rapport à d'autres méthodes comme le « zéro-padding + filtrage passe-bande ».

