Polynômes orthogonaux - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

- Introduction - Introduction - Propriétés - Exemple : les polynômes de Legendre - Tableau des polynômes orthogonaux classiques - Équations différentielles conduisant à des polynômes orthogonaux

Introduction

En mathématiques, une suite de polynômes orthogonaux est une suite infinie de polynômes p₀(x), p₁(x), p₂(x) ... à coefficients réels, dans laquelle chaque p_n(x) est de degré n, et telle que les polynômes de la suite sont orthogonaux deux à deux pour un produit scalaire de fonctions donné.

Introduction

Le produit scalaire de fonctions le plus simple est l'intégrale du produit de ces fonctions, sur un intervalle borné :

Plus généralement, on peut introduire une "fonction poids" W(x) dans l'intégrale (sur l'intervalle d'intégration $] a, b [$ , W doit être à valeurs finies et strictement positives, et l'intégrale du produit de la fonction poids par un polynôme doit être finie ; les bornes $a, b$ peuvent être infinies) :

Avec cette définition du produit scalaire, deux fonctions sont orthogonales entre elles si leur produit scalaire est égal à zéro (de la même manière que deux vecteurs sont orthogonaux (perpendiculaires) si leur produit scalaire égale zéro). On introduit alors la norme associée : $||f||=\sqrt{\langle f,f \rangle}$ ; le produit scalaire fait de l'ensemble de toutes les fonctions de norme finie un espace de Hilbert.

L'intervalle d'intégration est appelé intervalle d'orthogonalité.

Le domaine des polynômes orthogonaux a été développé durant le XIX^ème siècle par Stieltjes, comme outil de la théorie analytique des fractions continues. De multiples applications en ont découlé, en mathématiques et en physique.

Propriétés

Toute suite de polynômes $p_0, p_1 \dots$ , où chaque $\ p_k$ est de degré k, est une base de l'espace vectoriel $\R[x]$ (de dimension infinie) de tous les polynômes, "adaptée au drapeau $(\mathbb{R}_n[x])_{n\in\N}$ ". Une suite de polynômes orthogonaux est simplement une telle base qui est de plus orthogonale pour un certain produit scalaire. Ce produit scalaire étant fixé, une telle suite est presque unique (unique à produit près de ses vecteurs par des scalaires non nuls), et peut s'obtenir à partir de la base canonique $(1,x,x^2,\ldots)$ (non orthogonale en général), par le procédé de Gram-Schmidt.

Quand on construit une base orthogonale, on peut être tenté de la rendre orthonormale, c'est-à-dire telle que $\langle p_n, p_n \rangle\ =\ 1$ pour tout n, en divisant chaque $p n$ par sa norme. Dans le cas des polynômes, on préfère ne pas imposer cette condition supplémentaire car il en résulterait souvent des coefficients contenant des racines carrées. On préfère souvent choisir un multiplicateur tel que les coefficients restent rationnels, et donnent des formules aussi simples que possible. On appelle cela standardisation. Les polynômes "classiques" énumérés ci-dessous ont été ainsi standardisés ; typiquement, le coefficient de leur terme de plus haut degré ou leur valeur en un point ont été mis à une quantité donnée (pour les polynômes de Legendre, $P n (1) = 1$ ). Cette standardisation n'a pas de signification mathématique, c'est juste une convention, qui pourrait aussi parfois être obtenue par une mise à l'échelle de la fonction poids correspondante. Notons $h_n=\langle p_n,\ p_n \rangle$ (la norme de $\ p_n$ est la racine carrée de $\ h_n$ ). Les valeurs de $\ h_n$ pour les polynômes standardisés sont énumérées dans le tableau ci-dessous. Nous avons

;

où $δ m n h n$ est le delta de Kronecker.

Toute suite $(p k)$ de polynômes orthogonaux possède un grand nombre de propriétés remarquables. Pour commencer :

Lemme 1 : $(p_0,\ldots,p_n)$ est une base de $\R_n[x]$

Lemme 2 : $p n$ est orthogonal à $\R_{n-1}[x]$ .

Le lemme 1 est dû au fait que $p k$ est de degré k. Le lemme 2 vient de ce que, de plus, les $p k$ sont orthogonaux deux à deux.

Relation de récurrence

Pour toute suite de polynômes orthogonaux, il existe une relation de récurrence relativement à trois polynômes consécutifs.

Les coefficients $a n, b n, c n$ sont donnés par

où $k j$ et $k j'$ désignent les deux premiers coefficients de $p j$ :

et $h j$ le produit scalaire de $p j$ par lui-même :

(Par convention, $c_0,\ p_{-1},\ k'_0$ sont nuls.)

Avec les valeurs données pour $a n$ et $b n$ , le polynôme $(a n x + b n) p n - p n + 1$ est de degré <n (les termes de degrés n+1 et n s'éliminent). On peut donc l'exprimer sous forme d'une combinaison linéaire des éléments de la base $(p_j)_{j=0}^{n-1}$ de $\R_{n-1}[x]$ :

avec

(car pour j<n, p_j est orthogonal à p_n et p_n₊₁).

De plus, de par la forme intégrale du produit scalaire,

Pour j<n-1, ce produit scalaire est nul car $x p j$ est de degré <n.

Pour j=n-1, il est égal à $\frac{h_n}{a_{n-1}}$ car (par le même raisonnement qu'au début) $a n - 1 x p n - 1 - p n$ est de degré <n.

On peut conclure :

avec

Noyau de Christoffel-Darboux

Dans l'espace L2 associé à W, notons $S n$ la projection orthogonale sur $\R_n[x]$ : pour toute fonction $f\$ telle que $\int_a^bf^2(x)W(x)dx<\infty$ ,

où $K n$ est le noyau de Christoffel-Darboux, défini par :

La relation de récurrence précédente permet alors de montrer :

Démontrons la première de ces deux formules (la seconde s'en déduit en faisant tendre y vers x), par récurrence sur n. Pour n=-1 elle est vraie (par convention, K_-1=0). Supposons-la vraie au rang n-1 et prouvons-la au rang n. En remplaçant $p n + 1$ on obtient

avec, par hypothèse de récurrence,

d'où

Existence de racines réelles

Tout polynôme d'une suite de polynômes orthogonaux dont le degré n est supérieur ou égal à 1 admet n racines distinctes, toutes réelles, et situées strictement à l'intérieur de l'intervalle d'intégration (c'est une propriété remarquable : il est rare, pour un polynôme de degré élevé dont les coefficients ont été choisis au hasard, d'avoir toutes ses racines réelles)

Soit m le nombre des points où p_n change de signe à l'intérieur de l'intervalle d'orthogonalité ; notons $x_1 \dots x_m$ ces points. Ce sont les racines de p_n d'ordre impair appartenant à l'intervalle d'orthogonalité. D'après le théorème fondamental de l'algèbre, m ≤ n. On va montrer m = n. Soit $S(x) = \prod_{j=1}^m (x - x_j)$ ; c'est un polynôme de degré m qui change de signe en chaque point x_j ; S(x)p_n(x) est donc strictement positif, ou strictement négatif, partout sur l'intervalle d'intégration sauf aux points x_j , et il en est donc de même de S(x)p_n(x)W(x). Ainsi, $\langle S, p_n \rangle$ , l'intégrale de ce produit, est non nul. Mais, d'après le Lemme 2, p_n est orthogonal à tous les polynômes de degré inférieur, donc le degré de S doit être n.

Position des racines

Les racines des polynômes se trouvent strictement entre les racines du polynôme de degré supérieur dans la suite.

On met d'abord tous les polynômes sous une forme standardisée telle que le coefficient dominant soit positif (ce qui ne change pas les racines), puis on effectue une récurrence sur n. Pour n=0 il n'y a rien à démontrer. Supposons la propriété acquise jusqu'au rang n. Notons $x_1<\ldots<x_n$ les racines de $p n$ et $y_0<\ldots<y_n$ celles de $p n + 1$ . La relation de récurrence donne $p n + 1 (x j) = - c n p n - 1 (x j)$ avec (d'après le choix de standardisation) $c n > 0$ . Or par hypothèse de récurrence, $( - 1) n - j p n - 1 (x j) > 0$ . On en déduit $( - 1) n + 1 - j p n + 1 (x j) > 0$ . En outre, $\forall x>y_n,\ p_{n+1}(x)>0$ et $\forall x<y_0,\ (-1)^{n+1}p_{n+1}(x)>0$ . Ceci permet de conclure : $y_0<x_1<y_1<\ldots<x_n<y_n$ .

Une autre méthode de démonstration est de prouver (par récurrence, ou plus simplement en utilisant le noyau de Christoffel-Darboux) que pour tout n et tout x, $p_{n+1}'(x)\ p_n(x) > p_{n+1}(x)\ p_n'(x)\,$ , pour en déduire que $p_{n+1}'(y_j)\,$ et $p_n(y_j)\,$ ont même signe, si bien que $( - 1) n - j p n (y j) > 0$ , ce qui permet de conclure que p_n s'annule entre les $y j$ .

Exemple : les polynômes de Legendre

Des chercheurs ont intercepté et décodé "facilement" des messages sous-marins depuis les airs ✈️

Voici pourquoi les femmes ressentent plus de douleur chronique que les hommes 🧐

Un immense iceberg se détache et révèle ce surprenant écosystème caché depuis des siècles 🦑

Un colibri imite une chenille pour survivre: découvrez cette stratégie étonnante 🐦

Cette forme de vie géante ne ressemble à rien d'autre de connu 🌎

Recherche: cette nouvelle molécule se montre terriblement efficace contre la calvitie 🔬

Des tardigrades génétiquement modifiés pour survivre... sur le Soleil ☀️

Or, cuivre... pourquoi ces métaux précieux se trouvent dans les zones de subduction ? 🌋

Une avancée majeure dans la supraconductivité à haute température 🌡️

Les chimpanzés fabriquent leurs outils comme des ingénieurs 🛠️

Mach 16: des surprises physiques découvertes dans les vols hypersoniques 🚀

Des scientifiques ont réussi à filmer la séparation des brins de l'ADN 🧬

Comment 13 degrés peuvent transformer une foule en chaos ? 🚶‍♂️

Découverte: les phoques savent évaluer la quantité d'oxygène dans leur sang 🤿

Les trous noirs pourraient-ils protéger la vie dans l'Univers ? 🌱

Implants mammaires en silicone: une toxicité remise en question

Découverte d'un "broyeur d'étoiles" dans notre galaxie 🌀

Neandertal et Sapiens: une collaboration inattendue il y a 100 000 ans 🤝

Que sont ces sphères noires découvertes sur Mars ? 🔴

Découverte: les champignons contre... la grippe 🍄

Découverte inédite: voici comment notre cerveau attache les souvenirs entre eux 🖇️

Découverte: les requins émettent des sons ! 🦈

Des scientifiques inventent un ciment qui aspire le CO2 🌎

Comment le cerveau décode le temps: du récit aux réseaux neuronaux 🕰️

Horloge atomique dans son smartphone et GPS 1000 fois plus précis ⏱️

L'ocytocine: l'hormone de l'amour ou simple mythe ? 💕

Cette nouvelle méthode permettrait d'identifier une vie extraterrestre inattendue 👽

Une araignée géante filmée dans les abysses de l'Antarctique 🕷️

Trois familles d'objets planétaires identifiées dans le Système Solaire externe 🔭

Découverte d'une griffe géante à 2 doigts appartenant à une nouvelle espèce de dinosaure 🦖

Un nouvelle technologie de vaisseaux nucléaires pour explorer le Système solaire 🚀

Votre cerveau se consomme lui-même pour rechercher de l'énergie, pendant un exercice intense 🏃

Des vestiges organiques découverts sur Mars 🧐

Pourquoi les armes nucléaires sont-elles si difficiles et dangereuses à produire ? 💥

MRAM: la mémoire informatique du futur enfin là ? ⚡

Vidange d'un lac glaciaire: la catastrophe du Sikkim en 2023 expliquée 🌊

Qu'est-ce qui a bien pu creuser ces tunnels il y a plusieurs millions d'années ? ⛏️

Les secrets biologiques de Maria Branyas Morera, la supercentenaire 🕰️

Exploiter la rotation terrestre pour produire de l'électricité: l'énergie verte de demain ? ⚡

Que nous apprennent ces émissions radioactives de météorites lunaires et martiennes ? ☢️

La Chine dévoile un processeur quantique fiable à 99.9% 🚀

Voici les emblématiques espèces éteintes que la science pourrait bientôt ressusciter 🦣

Un poisson de 15 millions d'années incroyablement préservé 🐟

Existence démontrée d'un lien entre muscles, cerveau et fertilité 💪

Prévision météo: une petite IA surclasse les plus gros supercalculateurs 🌤️

Cette fosse commune révèle une histoire de violence extrême et de respect ⚔️

Cette nouvelle cellule solaire dure 10 fois plus longtemps 🌞

Des poulets aux plumes de dinosaures: une expérience génétique étonnante 🐔

Vidéo: une pieuvre chevauche un requin, une rencontre marine insolite 🐙

Une explication aux magnétars à faible champ 🧲

Page générée en 0.138 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales
Version anglaise | Version allemande | Version espagnole | Version portugaise