Décomposition d'Adomian
# Décomposition d'Adomian

La décomposition d'Adomian est une méthode semi-analytique de résolution d'équations différentielles développée par le mathématicien américain George Adomian (en) durant la seconde partie du XXe siècle. On rencontre fréquemment l'utilisation d'ADM pour Adomian Decomposition Method.

## Généralités

On considère le problème de Cauchy suivant :

${\dfrac {dy}{dt}}=f(t,y)\quad {\text{et}}\quad y(0)=y_{0}$

Cette équation vérifiée par y est générale dans la mesure où y peut être à valeurs vectorielles et que nous n'avons pas de condition sur f. Il faut bien noter que dans cette méthode, il est plus commode de considérer y comme un vecteur et f comme une fonction pour éviter des confusions :

${\begin{matrix}f:&\mathbb {R} \times \mathbb {R} ^{d}&\to &\mathbb {R} ^{d}\\&(t,y)&\mapsto &f(t,y)\end{matrix}}$

Considérant que f est analytique proche de $y=y_0$ et t=0, résoudre le problème initial revient à résoudre :

$y(t)=y_{0}+\int _{0}^{t}f\left(s,y(s)\right)\,\mathrm {d} s$

## Méthode d'Adomian

La méthode d'Adomian consiste à décomposer y comme une série :

$y=y_{0}+\sum _{n=1}^{\infty }y_{n}$

et à décomposer de la même manière la fonction f :

$f(t,y)=\sum _{n=0}^{\infty }A_{n}(t,y_{0},y_{1},...,y_{n}),$

où les fonctions $A_n$ sont les polynômes d'Adomian, calculés formellement comme-ci :

$A_{n}={\frac {1}{n!}}{\frac {\mathrm {d} ^{n}}{\mathrm {d} \varepsilon ^{n}}}f\left(t,\sum _{i=1}^{\infty }y_{i}\varepsilon ^{i}\right){\Bigg |}_{\varepsilon =0}.$

En injectant les deux premières décompositions dans l'équation intégrale, on peut en déduire une méthode itérative de calcul des $y_n$ :

$y_{n+1}=\int _{0}^{t}A_{n}\left(s,y_{0}(s)...,y_{n}(s)\right)\,\mathrm {d} s$

On part de la formulation intégrale

    $y(t)=y_{0}+\int _{0}^{t}f\left(s,y(s)\right)\,\mathrm {d} s\quad \Rightarrow \quad y(t)-y_{0}=\int _{0}^{t}f\left(s,y(s)\right)\,\mathrm {d} s$

Or la décomposition de y donne :

    $y(t)-y_{0}=\sum _{n=1}^{\infty }y_{n}(t)$

Soit par égalité :

    $\sum _{n=1}^{\infty }y_{n}(t)=\int _{0}^{t}f\left(s,y(s)\right)\,\mathrm {d} s$

Il suffit d'injecter la décomposition de la fonction f :

    $\sum _{n=1}^{\infty }y_{n}=\int _{0}^{t}\sum _{n=0}^{\infty }A_{n}(t,y_{0},y_{1},...,y_{n})\,\mathrm {d} s$

Dans un premier temps on ré-indice la série de gauche. Ensuite sous réserve que la convergence de la série des $A_n$ vers f est uniforme sur le segment [0,t] (ce qui est valide avec l'hypothèse d'analyticité de f), on intervertit somme et intégrale :

    $\sum _{n=0}^{\infty }y_{n+1}=\sum _{n=0}^{\infty }\int _{0}^{t}A_{n}(t,y_{0},y_{1},...,y_{n})\,\mathrm {d} s$

Puis on identifie les termes. Évidemment ce n'est pas l'unique solution mais elle nous satisfait entièrement. En effet, le calcul de $y_n$ en lui-même n'est pas très important, c'est la série qui découlera de ce calcul itératif qui nous intéressera (celle-ci construisant la solution y) :

$y_{n+1}=\int _{0}^{t}A_{n}(t,y_{0},y_{1},...,y_{n})\,\mathrm {d} s$

Pouvant calculer à la suite chaque $y_n,$ on construit au fur et à mesure la solution finale y par sommation.

## Polynômes en dimension 1

Dans cette section on se restreint à d=1. Nous allons expliciter les premiers polynômes d'Adomian (les dérivées correspondent à des dérivées partielles par rapport à y) :

    $A_{0}={\frac {1}{0!}}\ f\left(t,\sum _{i=1}^{\infty }y_{i}\varepsilon ^{i}\right){\Bigg |}_{\varepsilon =0}=f\left(t,y_{0}\right),$
    $A_{1}={\frac {1}{1!}}\ {\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}f\left(t,\sum _{i=1}^{\infty }y_{i}\varepsilon ^{i}\right){\Bigg |}_{\varepsilon =0}=y_{1}{\frac {\partial }{\partial y}}f\left(t,y_{0}\right),$
    $A_{2}={\frac {1}{2!}}\ {\frac {\mathrm {d} ^{2}}{\mathrm {d} \varepsilon ^{2}}}f\left(t,\sum _{i=1}^{\infty }y_{i}\varepsilon ^{i}\right){\Bigg |}_{\varepsilon =0}=y_{2}{\frac {\partial }{\partial y}}f\left(t,y_{0}\right)+{1 \over 2}y_{1}{\frac {\partial ^{2}}{\partial y^{2}}}f\left(t,y_{0}\right).$
