Série de Neumann

En analyse fonctionnelle, une série de Neumann est une série d'opérateurs de la forme

$\sum _{k=0}^{\mathrm{\infty }}{T}^k.$

où T est un opérateur et T^k désigne une itération de T répétée k fois. Elle étend l'idée de série géométrique.

La série est nommée d'après le mathématicien Carl Neumann, qui l'a utilisé en 1877 dans le cadre de la théorie du potentiel. Elle est également centrale dans l'étude du spectre d'opérateurs bornés.

Propriétés

Supposons que T est un opérateur linéaire borné dans un espace vectoriel normé X. Si la série de Neumann converge pour la norme d'opérateur, alors l'opérateur Modèle:Math est inversible et son inverse est l'opérateur somme de la série :

$(\mathrm{I}\mathrm{d}-T{)}^{-1}=\sum _{k=0}^{\mathrm{\infty }}{T}^k$,

avec Modèle:Math est l'identité sur X. Pour s'en convaincre, on peut regarder les sommes partielles

$S_n:=\sum _{k=0}^n{T}^k$.

On a alors, par télescopage,

$\underset{n\to \mathrm{\infty }}{\lim }(\mathrm{I}\mathrm{d}-T)S_n=\underset{n\to \mathrm{\infty }}{\lim }(\sum _{k=0}^n{T}^k-\sum _{k=0}^n{T}^{k+1})=\underset{n\to \mathrm{\infty }}{\lim }(\mathrm{I}\mathrm{d}-T^{n+1})=\mathrm{I}\mathrm{d}.$

On reconnait le résultat analogue pour la série géométrique sur la droite réelle :

$\frac{1}{1-x}=1+x+x^2+\cdots$

Un cas où la convergence est garantie est si X est un espace de Banach et |T| < 1 pour la norme d'opérateur. Il existe cependant des résultats montrant la convergence pour des conditions plus faibles.

Exemples

Opérateur matriciel

Soit $C\in {ℝ}^{n\times n}$ une matrice carrée réelle. On peut voir que la série de Neumann liée à Modèle:Mvar converge pour une norme matricielle si :

$||C|{|}_{\mathrm{\infty }}=\underset{i}{\max }\sum _j|c_{ij}|<1.$

et la série converge vers la matrice inverse de Modèle:Math, qui existe.

Ce résultat peut s'appliquer dans l'étude de convergence et de stabilité d'un schéma numérique : pour un problème de la forme

${𝐲}^{n+1}=A{𝐲}^n+b$

si la série de Neumann de Modèle:Mvar converge, alors le schéma converge vers une solution du problème

$𝐲=A𝐲+b.$

Opérateur intégral (série de Liouville-Neumann)

Pour une équation intégrale de Fredholm du second type

$\varphi (x)-\lambda \underset{a}{\overset{b}{\int }}K(x,t)\varphi (t)\mathrm{d}t=f(x)$

on peut définir la série, dite de Liouville-Neumann, pour des valeurs de λ suffisamment petites pour assurer la convergence de la série :

$\varphi (x)=\sum _{n=0}^{+\mathrm{\infty }}{\lambda }^n{\varphi }_n(x)$

avec

${\varphi }_0(x)=f(x),\text{ et }\mathrm{\forall }n\in \mathbb{N},{\varphi }_{n+1}(x)=\underset{a}{\overset{b}{\int }}K(x,t){\varphi }_n(t)\mathrm{d}t=\underset{a}{\overset{b}{\int }}{K}_{n+1}(x,t)f(t)\mathrm{d}t.$

Le noyau Modèle:Mvar est défini par la convolution itérée du noyau Modèle:Mvar avec lui-même n-1 fois :

$K_n(x,z)=\underset{x_{n-1}=a}{\overset{b}{\int }}\cdots \underset{x_1=a}{\overset{b}{\int }}K(x,x_1)K(x_1,x_2)\dots K(x_{n-1},z)\mathrm{d}{x}_1\mathrm{d}{x}_2\dots \mathrm{d}{x}_{n-1}.$

La série de Liouville-Neumann est l'unique solution continue de l'équation intégrale de Fredholm.

Le concept peut être étendu pour la résolution des équations intégrales de Volterra.

L'ensemble des opérateurs inversibles est ouvert

Un corollaire est que l'ensemble des opérateurs inversibles entre deux espaces de Banach B et BModèle:' est ouvert pour la topologie induite par la norme d'opérateur. En effet, soit S : B → BModèle:' un opérateur inversible et T: B → B' un autre opérateur.

Si |S – T | < |S⁻¹|⁻¹, alors T est également inversible. En effet, on peut remarquer que

$T=S(\mathrm{I}\mathrm{d}-(\mathrm{I}\mathrm{d}-S^{-1}T)).$

On en tire alors que la norme de T⁻¹ est majorée par

$|T^{-1}|\le \frac{1}{1-q}|S^{-1}|\text{avec}q=|S-T||S^{-1}|.$

Références

• Modèle:Ouvrage