Voir la source de Approximant de Padé de la fonction exponentielle

[[Image:Henri Padé.jpeg|thumb|Les travaux présentés dans cet article sont l'œuvre du [[mathématicien]] [[France|français]] [[Henri Padé]] ([[1863 en science|1863]]-[[1953 en science|1953]]).]]
En [[mathématiques]], un '''[[approximant de Padé]] de la [[fonction exponentielle]]''' est une [[fraction rationnelle]] ''h''(''x'')/''k''(''x''), où ''h''(''x'') désigne un [[polynôme]] de degré ''p'' et ''k''(''x'') de degré ''q'', telle que le [[développement limité]] de la fraction à l'ordre ''p'' + ''q'' soit identique à celui de l'exponentielle. L'étude de cette question est l'exemple introductif choisi par [[Henri Padé]] pour la théorie des approximants portant son nom.

L'existence de suites de fractions rationnelles ayant pour limite l'exponentielle est une question déjà abordée avant les travaux de Padé. [[Leonhard Euler]] ouvre le bal avec deux expressions dont l'une fournit le [[Fraction continue et approximation diophantienne#Exemple : le nombre e|développement en fraction continue infinie de {{math|e}}]], prouvant ainsi son [[nombre irrationnel|irrationalité]] ; [[Jean-Henri Lambert]] démontre plus rigoureusement ce développement, [[Joseph-Louis Lagrange]] en trouve deux autres<ref>J.-L. Lagrange, ''Sur l’usage des fractions continues dans le Calcul intégral'', Nouveaux Mémoires de l’Académie Royale des Sciences et Belles-Lettres, t. 4, {{p.|301}} (1779).</ref> et [[Carl Friedrich Gauss]] encore un<ref>{{Article|langue=la|prénom=C. F.|nom=Gauss|titre=Disquisitiones generales circa seriem infinitam…|revue=Commentationes Societatis Regiae Scientiarum Gottingensis recentiores|année=1813|url=http://gdz.sub.uni-goettingen.de/dms/load/img/?PPN=PPN35283028X_0002_2NS&DMDID=DMDLOG_0024&LOGID=LOG_0024&PHYSID=PHYS_0239|sous-titre=Sectio secunda — Fractiones continuae|p.=16}}.</ref>.

Le travail de Padé consiste à généraliser ces travaux précédents en vue d'illustrer par un exemple une théorie générale s'appliquant à toute [[fonction analytique]]. Il traite cette question sous quatre aspects : il montre l'existence d'un approximant de Padé d'indice (''p'', ''q''), établit les [[Suite définie par récurrence|relations de récurrence]] permettant de déterminer un approximant d'ordre supérieur, en déduit les différentes expressions sous forme de [[fraction continue généralisée|fractions continues (généralisées)]] de l'exponentielle et montre la [[convergence uniforme]] de certaines suites d'approximants.

La présentation ici correspond à une reformulation<ref name=Pade1899>{{Article|prénom=H.|nom=Padé|titre=Mémoire sur les développements en fractions continues de la fonction exponentielle|revue=[[Annales scientifiques de l'École normale supérieure|ASENS]]|série=3|année=1899|url=http://archive.numdam.org/ARCHIVE/ASENS/ASENS_1899_3_16_/ASENS_1899_3_16__395_0/ASENS_1899_3_16__395_0.pdf|p.=395-426}}.</ref> de [[1899 en science|1899]] et un enrichissement d'une partie de son travail de thèse<ref>{{Article|prénom=Henri |nom=Padé|titre=Sur la représentation approchée d'une fonction par des fractions rationnelles|revue=ASENS|série=3|vol=9|année=1892|url=http://archive.numdam.org/ARCHIVE/ASENS/ASENS_1892_3_9_/ASENS_1892_3_9__S3_0/ASENS_1892_3_9__S3_0.pdf|p.=3-93}} ([https://archive.org/details/surlareprsentat00padgoog thèse]).</ref>.

== Réduite de la fonction exponentielle ==
Dans la suite de l'article, ''p'' et ''q'' désignent deux [[Entier naturel|entiers positifs]]. Le premier résultat établit l'existence et l'[[Unicité (mathématiques)|unicité]], [[À quelque chose près|à un facteur multiplicatif près]], de deux polynômes ''h<sub>p,q</sub>'' et ''k<sub>p,q</sub>'' de degrés respectifs ''p'' et ''q'', tels que :

* ''Les développements en [[série entière|séries entières]] au point'' 0 ''de la fraction rationnelle h<sub>p,q</sub>/k<sub>p,q</sub> et de la fonction exponentielle coïncident sur les p + q'' + 1 ''premiers termes.''

Une telle fraction rationnelle est appelée « approximant de Padé ». Pour établir ce résultat, le mathématicien se fonde sur l'expression suivante d'une [[primitive]] de e<sup>''tx''</sup>''f''(''x'') où ''t'' est un paramètre, ''x'' la variable et {{math|''f''}} un polynôme dont le degré est noté ''n'' :
<center><math>\int \exp (tx)f(x){\rm d}x= \exp (tx)\sum_{i=0}^n (-1)^i\frac {f^{(i)}(x)}{t^{i+1}}.</math></center>

Il en déduit les expressions suivantes :
<center><math>h_{p,q}(t)=\sum_{i=0}^p\frac{p!/(p-i)!}{(p+q)!/(p+q-i)!}\frac{t^i}{i!}\quad\text{et}\quad k_{p,q}(t)=\sum_{j=0}^q (-1)^j\frac{q!/(q-j)!}{(p+q)!/(p+q-j)!}\frac{t^j}{j!}.</math></center>

La configuration présente de nombreuses analogies avec les [[fraction continue|fractions continues]], ce qui justifie la définition suivante :

: La fraction rationnelle ''h<sub>p,q</sub>''/''k<sub>p,q</sub>'' est appelée '''réduite d'ordre ou d'indice (''p'', ''q'') '''de la fonction exponentielle.

On dispose de propriétés comme :
*''si'' (''r'', ''s'') ''est égal à'' (''p'' + 1, ''q''), (''p'', ''q'' + 1) ''ou'' (''p'' + 1, ''q'' + 1), ''alors h''{{sub|''p'',''q''}}''k''{{sub|''r'', ''s''}} – ''h''{{sub|''r'', ''s''}}''k''{{sub|''p'',''q''}} ''est un [[Monôme (mathématiques)|monôme]] (non nul) de degré p'' + ''q'' + 1 ;
* ''les deux polynômes h<sub>p,q</sub> et k<sub>p,q</sub> sont uniques (à [[Polynôme unitaire|normalisation]] près) et premiers entre eux.''

{{Démonstration/début|titre=Démonstrations}}
Soit ''n'' l'entier ''p + q''. 
* '''Pour tout polynôme {{math|''f''}} de degré ''n'' :'''<center><math>\int \exp (tx)f(x){\rm d}x=\exp(tx)\sum_{k=0}^n (-1)^k\frac {f^{(k)}(x)}{t^{k+1}}.</math></center>La vérification est immédiate, en dérivant par rapport à {{math|''x''}} le membre de droite.
*Soient ''a'' et ''b'' deux réels et {{math|''f''}}, {{math|α}} et {{math|β}} les polynômes définis par <center><math>f(x)=(x-a)^q(x-b)^p,\quad\alpha(t)=\sum_{i=0}^p(-1)^if^{(n-i)}(a)t^i,\quad\beta(t)=\sum_{j=0}^q(-1)^jf^{(n-j)}(b)t^j.</math></center>
**'''L'égalité suivante est vérifiée :'''<center><math>\beta(t)\exp(bt)-\alpha(t)\exp(at)=(-1)^nt^{n+1}\int_a^b \exp(tx)f(x){\rm d}x.</math></center>En effet, puisqu'ici ''f''<sup>(''k'')</sup>(''a'') est nul si ''k'' est strictement inférieur à ''q'' et ''f''<sup>(''k'')</sup>(''b'') est nul si ''k'' est strictement inférieur à ''p'', l'expression précédente, appliquée entre ''a'' et ''b'', montre que<center><math>\begin{align}t^{n+1}\int_a^b \exp (tx)f(x){\rm d}x&=\exp (tb)\sum_{k=p}^n (-1)^kf^{(k)}(b)t^{n-k}-\exp (ta)\sum_{k=q}^n (-1)^kf^{(k)}(a)t^{n-k}\\
&=(-1)^n\Big(\exp(tb)\beta(t)-\exp(ta)\alpha(t)\Big).\end{align}</math></center>
**{{Citation|Les '''expressions explicites des coefficients de α et β''' se déduisent aisément des [[Formule de Taylor|développements de {{math|''f''}} suivant les puissances]], soit de ''x – a'', soit de ''x – b''<ref name=Pade1899/>.}}<br />Par exemple pour α :<center><math>f(x)=(x-a)^q\left((x-a)+(a-b)\right)^p=(x-a)^q\sum_{i=0}^p{p\choose i}(x-a)^{p-i}(a-b)^i,</math></center>si bien que<center><math>f^{(n-i)}(a)=(n-i)!{p\choose i}(a-b)^i\quad{\rm  et}\quad\alpha(t)=\sum_{i=0}^p(n-i)!{p\choose i}(b-a)^it^i.</math></center>
* '''Si ''h<sub>p,q</sub>'' et ''k<sub>p,q</sub>'' sont les deux fonctions polynomiales suivantes alors le premier monôme non nul de la [[fonction entière]] ''k<sub>p,q</sub>''e<sup>''t''</sup> – ''h<sub>p,q</sub>'' est de degré ''p'' + ''q'' + 1 :'''<center><math>h_{p,q}(t)=\sum_{i=0}^p\frac{(p+q-i)!p!}{(p+q)!(p-i)!i!}t^i\quad\text{et}\quad k_{p,q}(t)=\sum_{j=0}^q (-1)^j\frac{(p+q-j)!q!}{(p+q)!(q-j)!j!}t^j.</math></center>En effet, ''h<sub>p,q</sub>'' et ''k<sub>p,q</sub>'' sont les quotients par ''n''! des expressions ci-dessus des polynômes α et β, pour ''a'' = 0 et ''b'' = 1. Pour ces valeurs de ''a'' et ''b'', l'égalité précédente devient donc :<center><math>k_{p,q}(t)\exp(t)-h_{p,q}(t)=\frac{(-1)^n}{n!}t^{n+1}\int_0^1 \exp(tx)f(x){\rm d}x</math></center>or la fonction entière qui à ''t'' associe l'intégrale de 0 à 1 de {{math|''x'' ↦ exp(''tx'')''f''(''x'')}} est non nulle en ''t'' = 0.

Une fois prouvée l'existence de la fraction réduite et déterminée son expression, il reste à démontrer son unicité. On le fait grâce à un lemme dont la démonstration est de même nature que celle montrant que le numérateur et dénominateur d'une fraction réduite sont premiers entre eux.
*'''Si (''r'', ''s'') est égal à (''p'' + 1, ''q'') ou (''p'', ''q'' + 1) ou (''p'' + 1, ''q'' + 1), alors ''h''{{sub|''p'',''q''}}''k''{{sub|''r'', ''s''}} – ''h''{{sub|''r'', ''s''}}''k''{{sub|''p'',''q''}} est un monôme (non nul) de degré ''p'' + ''q'' + 1 :'''<br />En effet, d'une part ce polynôme est de degré inférieur ou égal à ''p'' + ''q'' + 1 et d'autre part, d'après l'égalité<center><math>h_{p,q}k_{r,s}-h_{r,s}k_{p,q} =k_{p,q}\Big(k_{r,s}\exp-h_{r,s}\Big)-k_{r,s}\Big(k_{p,q}\exp-h_{p,q}\Big),</math></center>son premier monôme non nul est de degré ''p'' + ''q'' + 1.

* '''Les polynômes ''h<sub>p,q</sub>'' et ''k<sub>p,q</sub>'' sont premiers entre eux :'''<br />En effet, d'après la propriété précédente, tout polynôme divisant à la fois ''h''{{sub|''p'',''q''}} et ''k''{{sub|''p'',''q''}} divise un monôme ; on en déduit que les seuls facteurs communs aux deux polynômes sont des puissances de ''t''. Comme ces deux polynômes ne sont pas divisibles par ''t'', ils sont bien premiers entre eux.
* '''Les polynômes ''h<sub>p,q</sub>'' et ''k<sub>p,q</sub>'' sont uniques :'''<br />Soient ''h'' et ''k'' deux polynômes, de degrés respectivement majorés par ''p'' et ''q'', tels que le développement en série entière en 0 de la fraction ''h''/''k'' coïncide jusqu'au degré ''n'' avec celui la fonction exponentielle. Par le même raisonnement que ci-dessus, ''h''{{sub|''p'',''q''}}''k'' – ''hk''{{sub|''p'',''q''}} est un polynôme de degré inférieur ou égal à ''n'' sans monômes de degrés inférieurs ou égaux à ''n''. Cette différence est donc nulle. Comme les polynômes ''h<sub>p,q</sub>'' et ''k<sub>p,q</sub>'' sont premiers entre eux et de degrés respectifs ''p'' et ''q'', le couple (''h'', ''k'') est proportionnel au couple (''h<sub>p,q</sub>'', ''k<sub>p,q</sub>'').
{{Démonstration/fin}}

== Table de Padé ==
Ainsi, la fonction exponentielle s'approxime par des fractions rationnelles, de manière un peu analogue à l'approximation par des polynômes avec les séries entières. Si les polynômes forment une suite, les approximants de Padé définissent un tableau à double entrée appelé {{Lien|trad=Padé table|table de Padé}}, dans lequel l'approximant d'indice (''p'', ''q'') figure dans la ''colonne p'' et la ''ligne q''. Les premiers termes sont les suivants<ref>Cette table est extraite de {{Harvsp|Padé|1892|p=16}}.</ref> :

{| class="wikitable" style="text-align: center; width: 100%"
! scope="col" | Table de Padé 
! scope="col" | 0
! scope="col" | 1
! scope="col" | 2
! scope="col" | 3
! scope="col" | …
|-
! scope="row" | 0
| <math>1\;</math>
| <math> 1 + x\; </math>
| <math> 1 + x + \frac {x^2}2\; </math>
| <math> 1 + x + \frac {x^2}2 + \frac {x^3}6\; </math>
| …
|-
! scope="row" | 1
| <math> \frac 1{1 - x}</math>
| <math> \frac {1 + \frac 12x}{1 - \frac 12x}</math>
| <math> \frac {1 + \frac 23 x + \frac 16 x^2}{1 - \frac 13x}</math>
| <math> \frac {1 + \frac 34x + \frac 14x^2 + \frac 1{24}x^3}{1 - \frac 14x}</math>
| …
|-
! scope="row" | 2
| <math> \frac 1{1 - x + \frac {x^2}2}</math>
| <math> \frac {1 + \frac 13x}{1 - \frac 23x+ \frac 16x^2}</math>
| <math> \frac {1 + \frac 12x + \frac 1{12}x^2}{1 - \frac 12x + \frac 1{12}x^2}</math>
| <math> \frac {1 + \frac 35 x + \frac 3{20}x^2 + \frac 1{60}x^3}{1 - \frac 25x + \frac 1{20}x^2}</math>
| …
|-
! scope="row" | 3
| <math> \frac 1{1 - x + \frac 12x^2 - \frac 16 x^3}</math>
| <math> \frac {1 + \frac 14x}{1 - \frac 34x+ \frac 14x^2 - \frac 1{24}x^3}</math>
| <math> \frac {1 + \frac 25x + \frac 1{20}x^2}{1 - \frac 35x + \frac 3{20}x^2 - \frac 1{60}x^3}</math>
| <math> \frac {1 + \frac 12 x + \frac 1{10}x^2 + \frac 1{120}x^3}{1 - \frac 12x + \frac 1{10}x^2 - \frac 1{120}x^3}</math>
| …
|-
! scope="row" | …
| …
| …
| …
| …
| …
|}
</center>

La première ligne horizontale correspond au développement en série entière. Les lignes diagonales définies par l'égalité ''p'' + ''q'' = ''C'', où ''C'' est une constante positive donnée, correspondent à un ensemble de fractions rationnelles dite ''droite d'égale approximation''. On y trouve par exemple, si ''C'' est égal à 2, les couples (2,0), (1,1) et (0,2). Une fraction rationnelle de la table est dite ''plus avancée'' qu'une autre lorsque son coefficient ''C'' est plus élevé.

Le graphique suivant illustre la convergence de différentes suites extraites de la table de Padé. Ici la notation <code>Exp[p,q]</code> désigne l'approximant d'indice (''p'', ''q''). La fonction exponentielle, notée <code>Exp</code>, est illustrée en rouge. La première ligne du tableau correspond à la suite des polynômes de la série entière. Elle est illustrée en bleu et correspond à la suite <code>Exp[1,0]</code>, <code>Exp[2,0]</code>, <code>Exp[3,0]</code>, etc. En vert et pointillé est illustrée la suite correspondant à la deuxième ligne de la table, elle est formée des approximants <code>Exp[0,1]</code>, <code>Exp[1,1]</code>, <code>Exp[2,1]</code>, <code>Exp[3,1]</code>, etc. En violet, on trouve les fractions de la diagonale : <code>Exp[1,1]</code>, <code>Exp[2,2]</code>, <code>Exp[3,3]</code>, etc.

[[Image:Padé Expontielle.jpg|thumb|upright=2|center|Différentes suites extraites de la table de Padé convergent uniformément vers la fonction exponentielle sur tout intervalle borné.]]

== Formule de récurrence ==
[[Image:Table de Padé Expontielle.jpg|thumb|upright=1.5|right|Illustration d'une suite de couples choisis par Lagrange, ''commençant par le couple (0, 0)'', pour construire une fraction continue approchant la fonction exponentielle]]
Padé cherche ensuite à obtenir des suites à partir de sa table. Si (''p<sub>n</sub>'') et (''q<sub>n</sub>'') sont deux suites croissantes de couples d'indices, l'objectif est d'exprimer par récurrence la suite des « réduites » — déjà calculées — d'indice (''p<sub>n</sub>'', ''q<sub>n</sub>''), pour comprendre le cas général. Cette récurrence est aussi une étape pour exprimer la fonction exponentielle sous forme de fractions continues généralisées, justifiant ainsi le nom de « [[Fraction continue généralisée#Réduites|réduites]] » pour ces approximants. L'exemple illustré par la figure de droite est utilisé par Lagrange pour obtenir un résultat de cette nature. Dans cet exemple, si ''f''<sub>n</sub>(''x'') désigne le ''n''-ième terme de la suite, on a :
<center><math>f_0(x) = \exp_{[0,0]}=1,\quad f_1(x) = \exp_{[1,0]} = 1 + x,\quad f_2(x) = \exp_{[1,1]} = \frac {1 + \frac 12x}{1 - \frac 12x}= 1 + \frac x{1 - \frac 12x}</math></center>
puis :
<center><math>f_3(x) = \exp_{[2,1]} = \frac {1 + \frac 23 x + \frac 16 x^2}{1 - \frac 13x} = 1 + \frac x{1 - \cfrac{\frac 12x}{1 + \frac 16x}},\quad\text{etc.} </math></center>

Padé n'étudie l'incrémentation des indices (''p''<sub>''n''+1</sub> – ''p<sub>n</sub>'', ''q''<sub>''n''+1</sub> – ''q<sub>n</sub>'') que dans les trois cas les plus simples : (1,0), (0,1) et (1,1). Le premier cas correspond à un déplacement horizontal d'une case, le deuxième à un déplacement vertical d'une case et le troisième à la conjonction des deux. Cette restriction lui permet d'obtenir des relations de récurrence :

* ''Soit'' (''f<sub>n</sub>'') ''une suite de réduites dont le couple de la différence des indices entre f''<sub>''n''+1</sub> ''et f<sub>n</sub> correspond toujours à l'un des trois cas'' (1,0), (0,1) ''ou'' (1,1), ''alors il existe une relation de récurrence de type :''
<center><math>h_{(p_{n+2},q_{n+2})} =\alpha_{n+2}h_{(p_n,q_n)} + \beta_{n+2}h_{(p_{n+1},q_{n+1})}\quad\text{et}\quad k_{(p_{n+2},q_{n+2})} =\alpha_{n+2}k_{(p_n,q_n)} + \beta_{n+2}k_{(p_{n+1},q_{n+1})}</math></center>
:''où ''α<sub>''n''+2</sub> ''désigne un monôme (non nul) de degré'' 1 ''ou'' 2, ''et ''β<sub>''n''+2</sub> ''un polynôme de degré'' 0 ''ou'' 1 ''et de terme constant égal à'' 1.

Le tableau suivant résume les trois configurations possibles :
[[Image:Padé type 1,2 et 3.jpg|600px|center]]
Les flèches rouges indiquent les incréments utilisés pour passer de ''f<sub>n</sub>''(''x'') à  ''f''<sub>''n''+1</sub>(''x''). Pour les deux premières configurations, ils correspondent à (1,0) ou (0,1), c'est-à-dire que soit le degré du [[numérateur]], soit celui du dénominateur est incrémenté de 1 ; pour la dernière, ils sont tous deux incrémentés de 1. Les flèches vertes indiquent les incréments utilisés pour passer de ''f''<sub>''n''+1</sub> à ''f''<sub>''n''+2</sub>. L'illustration précédente indique le degré des deux polynômes α<sub>''n''+2</sub>(''x'') et β<sub>''n''+2</sub>(''x'') utilisés pour exprimer la formule de récurrence. Pour la configuration 3, le seul cas où β<sub>''n''</sub>(''x'') est une constante est celui où ''f<sub>n</sub>''(''x''), ''f''<sub>''n''+1</sub>(''x'') et ''f''<sub>''n''+2</sub>(''x'') se trouvent sur la [[diagonale principale]], d'indice un couple (''p'', ''p'').

{{Démonstration/début}}
Pour un ''n'' fixé, on notera pour simplifier (''p'', ''q''), (''r'', ''s'') et (''u'', ''v'') les couples correspondant aux indices ''n'', ''n'' + 1 et ''n'' +2.
* '''Il existe deux polynômes α et β tels que les égalités suivantes soient vérifiées. Les degrés de ces polynômes sont donnés par la figure du paragraphe :'''<center><math>\begin{cases}\alpha h_{p,q}+\beta h_{r,s}&=h_{u,v},\\\alpha k_{p,q}+\beta k_{r,s}&=k_{u,v}.\end{cases}</math></center>Ce [[système de Cramer]] a pour déterminant Δ = ''h''{{sub|''p'',''q''}}''k''{{sub|''r'', ''s''}} – ''h''{{sub|''r'', ''s''}}''k''{{sub|''p'',''q''}} et sa solution est :<center><math>\alpha=\frac{h_{u,v}k_{r,s}-h_{r,s}k_{u,v}}\Delta,\quad \beta=\frac {h_{p,q}k_{u,v}-h_{u,v}k_{p,q}}\Delta.</math></center>D'après la section précédente, Δ est un monôme de degré ''p'' + ''q'' + 1. Il en résulte — par le même argument pour son numérateur — que α est un monôme de degré 1 pour les configurations 1 et 2 et de degré 2 pour la configuration 3. De même, β est un monôme de degré 0 dans la configuration 1 et un polynôme de degré inférieur ou égal à 1 dans les configurations 2 et 3 (de degré 1 en général, mais 0 si les trois couples d'indices sont sur la diagonale principale).
* '''Le coefficient constant de β vaut 1 :'''<br />Immédiat, par identification des termes constants dans l'égalité<center><math>\alpha k_{p,q}+\beta k_{r,s}=k_{u,v}.</math></center>
{{Démonstration/fin}}

== Fraction continue ==
{{Voir|Fraction continue généralisée#Réduites{{!}}§ « Réduites » de l'article sur les fractions continues généralisées}}
La relation de récurrence permet d'écrire un approximant de Padé sous forme de fraction continue, en partant d'un bord du tableau. Plus précisément :
<center><math>\begin{matrix}\text{si }q_0=0 :&f_n =\beta_0 + \frac{\alpha_1\mid}{\mid \beta_1} + \frac{\alpha_2\mid}{\mid \beta_2} + \cdots + \frac{\alpha_n\mid}{\mid \beta_n},\\
\text{si }p_0=0 :&f_n =\frac{1\mid}{\mid \beta_0} +\frac{\alpha_1\mid}{\mid \beta_1}+ \cdots + \frac{\alpha_n\mid}{\mid \beta_n},\end{matrix}</math></center>
où les polynômes α<sub>''n''</sub> et β<sub>''n''</sub> sont définis pour ''n'' ≥ 2 par les relations du paragraphe précédent et pour ''n'' égal à 0 ou 1, par les choix suivants :
<center><math>\begin{matrix}\text{si }q_0=0 :&\beta_0= h_{p_0,0},&\alpha_1=h_{p_1,q_1}-h_{p_0,0}k_{p_1,q_1},&\beta_1=k_{p_1,q_1}~;\\
\text{si }p_0=0 :&\beta_0=k_{0,q_0},&\alpha_1 =k_{p_1,q_1}-k_{0,q_0}h_{p_1,q_1},&\beta_1=h_{p_1,q_1}.\end{matrix}</math></center>
(Si ''p''{{ind|0}} et ''q''{{ind|0}} sont tous les deux nuls, on fait le choix menant à l'expression la plus simple : {{cf.}} exemples ci-dessous.)

La fraction continue est dite régulière si les polynômes α<sub>''j''</sub> pour ''j'' > 1 sont tous de même degré, ainsi que les polynômes β<sub>''j''</sub>. Les tableaux du paragraphe précédent montrent qu'il en existe de trois types différents.

=== Fraction continue du premier type ===
[[Image:Padé type 1.jpg|thumb|left|Les fractions continue du premier type s'obtiennent en considérant une suite extraite de la table de Padé selon la première configuration du paragraphe précédent.]]

Les fractions continues du premier type s'obtiennent avec des polynômes β égaux à 1 et des monômes α du premier degré. L'étude des relations de récurrence montrent qu'elles s'obtiennent nécessairement à l'aide d'une suite de réduites correspondant à la première configuration du paragraphe précédent. Quitte à remonter cette suite jusqu'au bord du tableau, on remarque qu'il en existe une pour chaque case à la frontière du tableau (ou deux, pour la case (0, 0)). Lagrange a développé celle de case initiale (0, 0) avec pour première réduite β{{ind|0}} = ''h''{{ind|0,0}} = 1 et suivie des couples rouges sur la figure, obtenant :
<center><math>1 + \frac{x\mid}{\mid 1} - \frac{\frac 12 x \mid}{\mid 1} + \frac{\frac 16x \mid}{\mid 1} - \frac{\frac 16x \mid}{\mid 1} + \frac{\frac 1{10}x \mid}{\mid 1} - \frac{\frac 1{10}x \mid}{\mid 1} + \frac{\frac 1{14}x \mid}{\mid 1} - \frac{\frac 1{14}x \mid}{\mid 1}+\cdots</math></center>

En remarquant que {{math|1=e{{exp|''x''}} = 1/e{{exp|–''x''}}}}, elle équivaut à sa symétrique par rapport à la diagonale principale, démontrée par Gauss comme exemple de [[Fraction continue de Gauss|ses fractions continues pour les fonctions hypergéométriques]] :<center><math>\frac{1\mid}{\mid 1} - \frac{x \mid}{\mid 1} + \frac{\frac 12x \mid}{\mid 1} - \frac{\frac 16x \mid}{\mid 1} + \frac{\frac 1{6}x \mid}{\mid 1} - \frac{\frac 1{10}x \mid}{\mid 1} + \frac{\frac 1{10}x \mid}{\mid 1} - \frac{\frac 1{14}x \mid}{\mid 1}+\cdots.</math></center>

Il est possible d'en construire d'autres de cette nature, par exemple à l'aide de la série bleue illustrée sur la figure, éventuellement précédée (pour une expression plus simple) du couple (0, 2).

=== Fraction continue du deuxième type ===
[[Image:Padé type 2.jpg|thumb|left|Les fractions continue du deuxième type s'obtiennent comme précédemment, mais en suivant un déplacement soit toujours vertical soit toujours horizontal.]]

Les fractions continues de la deuxième catégorie correspondent à la configuration 2. Elles s'obtiennent à l'aide d'un déplacement continu soit vers le bas, soit vers la droite. Les numérateurs α sont des monômes du premier degré et les dénominateurs β, des polynômes du premier degré de terme constant égal à 1. Trois exemples sont illustrés sur la figure de gauche. Le plus simple, qui correspond à suite illustrée en vert, a pour réduite d'indice ''n'' la somme partielle de la série exp(''x''), [[Formule de fraction continue d'Euler#Fonction exponentielle|selon une formule générale d'Euler]]<ref>{{la}} L. Euler, ''[[Introductio in analysin infinitorum]]'', t. 1, [https://gallica.bnf.fr/ark:/12148/bpt6k69587.f392 1748, par. 368-373].</ref> :
<center><math>\frac{1\mid}{\mid 1} - \frac{x \mid}{\mid 1+x} - \frac{\frac 12x \mid}{\mid 1+ \frac 12x} - \frac{\frac 13x \mid}{\mid 1+\frac 13 x} - \frac{\frac 14x \mid}{\mid 1+\frac 14x} - \cdots </math></center>
Padé montre que l'exemple d'Euler ne correspond qu'à un cas particulier de fraction continue de cette nature. Il en existe en fait une infinité, exactement une par case de la frontière de la table de Padé (ou deux, pour la case (0, 0)).

=== Fraction continue du troisième type ===
[[Image:Padé type 3.jpg|thumb|left|Les fractions continue du troisième type se caractérisent par des dénominateurs égaux à des monômes du deuxième degré.]]

Le dernier type de configuration correspond à la troisième décrite dans le paragraphe sur les relations de récurrence. Ici, le passage d'une réduite à une autre s'obtient par un déplacement diagonale. Les monômes α sont toujours du second degré. Les polynômes β sont du premier degré, à l'exception de la fraction continue associée à la diagonale principale, en rouge sur la figure (dans ce cas, β = 1).

Lagrange découvre celle associée à la série bleue sur la figure. Elle fournit l'expression suivante :
<center><math>1 + x + \frac{\frac 12x^2\mid}{\mid 1-\frac 13x} + \frac{\frac 1{36}x^2 \mid}{\mid 1-\frac 1{15}x} + \frac{\frac 1{100}x^2 \mid}{\mid 1- \frac 1{35}x} + \cdots</math></center>

Les formules exactes sont :
<center><math> \forall n\ge2\quad \alpha_n(x) = \frac {x^2}{4(2n-3)^2}\quad\text{et}\quad \beta_n(x) = 1 - \frac x{(2n-1)(2n-3)}.</math></center>

Celle correspondant à la diagonale principale<ref>{{harvsp|Padé|1899|p=416}}, l'attribue à Lagrange et ne mentionne même plus Lambert. Dans les {{p.|38-39}} de {{Article|prénom=H.|nom=Padé|titre=Sur les fractions approchées d'une fonction par des fractions rationnelles|revue=ASENS|série=3|vol=9|année=1892|url=http://archive.numdam.org/ARCHIVE/ASENS/ASENS_1892_3_9_/ASENS_1892_3_9__S3_0/ASENS_1892_3_9__S3_0.pdf|p.=3-93}}, il citait seulement ce précurseur {{citation|comme ayant donné un développement en fraction continue de quelques fonctions}}, retrouvées par Lagrange {{citation|sous la forme la plus convenable pour nous}}.</ref> [[Fraction continue et approximation diophantienne#Exemple : le nombre e|a été découverte par Euler]], à partir de son développement de la fonction [[tangente hyperbolique]] (que [[Fraction continue et approximation diophantienne#Exponentielle d'un rationnel|Lambert démontrera]] en toute rigueur en 1761, puis qui deviendra un autre [[Fraction continue de Gauss#Applications à 0F1|exemple de fraction continue de Gauss]]). Elle n'est pas à proprement parler du troisième type car les dénominateurs sont des constantes à partir de la valeur 2 de l'indice ; c'est l'unique exception.
<center><math>1 + \frac{x\mid}{\mid 1-\frac 12x} + \frac{\frac 1{12}x^2 \mid}{\mid 1} + \frac{\frac 1{60}x^2 \mid}{\mid 1} + \frac{\frac 1{140}x^2 \mid}{\mid 1}\cdots\quad\text{avec}\quad \forall n\ge 2 \quad \alpha_n(x) = \frac {x^2}{4(2n-1)(2n-3)}.</math></center>

{{Démonstration/début}}
* '''Calcul pour la première fraction continue :'''
Le monôme α<sub>''n''</sub> est de la forme ''a''<sub>''n''</sub>''x''<sup>2</sup> et le polynôme β<sub>''n''</sub>, de la forme 1 + ''b<sub>n</sub>x''. La relation de récurrence montre que pour tout entier naturel ''n'', on a pour les numérateurs :
<center><math>h_{n+2,n+1} =a_{n+2}x^2h_{n,n-1}+(1+ b_{n+2}x)h_{n+1,n}</math></center>
et de même pour les dénominateurs.

Par identification des coefficients dominants, au numérateur et dénominateur,
<center><math>\begin{cases}\frac {(n+1)!}{(2n+3)!}&= b_{n+2}\frac {n!}{(2n+1)!} + a_{n+2}\frac {(n-1)!}{(2n-1)!},\\\frac {(n+2)!}{(2n+3)!}&= -b_{n+2}\frac {(n+1)!}{(2n+1)!} + a_{n+2}\frac {n!}{(2n-1)!}.\end{cases}</math></center>
En [[Système d'équations linéaires|résolvant le système]], on trouve

<center><math>a_{n+2} = \frac 1{4(2n+1)^2}\quad\text{et}\quad b_{n+2} = \frac {-1}{(2n+1)(2n+3)}.</math></center>

* '''Calcul pour la deuxième fraction continue :'''
Avec les mêmes notations que pour le calcul précédent mais compte tenu du fait qu'ici, β{{ind|''n''+2}} = 1, on a
<center><math>h_{n+2,n+2}=a_{n+2}x^2h_{n,n}+h_{n+1,n+1}.</math></center>
Par identification des coefficients dominants, on trouve
<center><math>a_{n+2} =\frac 1{4(2n+1)(2n+3)}.</math></center>
{{Démonstration/fin}}

== Convergence des réduites ==
Les fractions continues des trois types précédents sont toutes [[convergence uniforme|uniformément convergentes]] sur toute [[partie bornée]] du [[plan complexe]]. Plus précisément :
*'' Soient ''(''p<sub>n</sub>'') ''et'' (''q<sub>n</sub>'')'' deux suites d'entiers positifs dont l'une au moins tend vers l'infini. Si le rapport p<sub>n</sub>/q<sub>n</sub> tend vers une limite, finie ou infinie, la suite de fractions rationnelles h''<sub>(''p<sub>n</sub>, q<sub>n</sub>'')</sub>/''k''<sub>(''p<sub>n</sub>, q<sub>n</sub>'')</sub>'' converge uniformément sur tout ensemble borné.''

Le numérateur et le dénominateur tendent tous les deux vers des séries entières. Si ''p<sub>n</sub>''/''q<sub>n</sub>'' a pour limite ω, alors :
<center><math>\lim_{n \to \infty} h_{(p_n,q_n)}(x) = \exp \left(\frac {\omega x}{1 + \omega}\right)\quad \text{et}\quad 
\lim_{n \to \infty} k_{(p_n,q_n)}(x) = \exp \left(\frac {-x}{1 + \omega}\right),</math></center>
ces deux formules signifiant, dans le cas où ω est infinie, que le numérateur tend vers la fonction exponentielle et le dénominateur vers la [[fonction constante]] 1.

{{Démonstration/début}}
Si la limite ω de la suite (ω<sub>''n''</sub>) = (''p<sub>n</sub>''/''q<sub>n</sub>'') est infinie, on se ramène au cas où elle est nulle en remplaçant ''x'' par –''x'' et en intervertissant les lettres ''p'' et ''q'' ainsi que ''h'' et ''k''. Supposons désormais que ω est finie. La suite (''q<sub>n</sub>'') tend alors vers l'infini.

* '''Si ''a<sub>k,n</sub>'' et ''b<sub>k,n</sub>'' désignent les coefficients des termes du ''k''-ième degré des polynômes ''h''<sub>(''p<sub>n</sub>'', ''q<sub>n</sub>'')</sub>  et ''k''<sub>(''p<sub>n</sub>, q<sub>n</sub>'')</sub>, alors :'''<center><math>\lim_{n \to \infty} a_{k,n}=\frac1{k!}\left(\frac {\omega}{1+\omega}\right)^k \quad\text{et}\quad \lim_{n \to \infty} b_{k,n} = \frac 1{k!}\left(\frac {-1}{1+\omega}\right)^k.</math></center>En effet, par définition,<center><math>a_{k,n} =\begin{Bmatrix}\frac {p_n!(p_n+q_n -k)!}{(p_n-k)!(p_n+q_n)!k!}&\text{ si }p_n\ge k\\0&\text{ sinon}\end{Bmatrix}=\frac 1{k!} \prod_{i=0}^{k-1} \frac {p_n -i}{p_n + q_n -i}</math></center>et<center><math>b_{k,n}=\begin{Bmatrix}\frac{(-1)^kq_n!(p_n+q_n -k)!}{(q_n-k)!(p_n+q_n)!k!}&\text{ si }q_n\ge k\\0&\text{ sinon}\end{Bmatrix}=\frac{(-1)^k}{k!}\prod_{i=0}^{k-1}\frac {q_n-i}{p_n+q_n -i}.</math></center>On conclut en remarquant que pour chacune des ''k'' valeurs de ''i'' considérées,<center><math>\frac{p_n-i}{p_n+q_n-i}=\frac{\omega_n-i/q_n}{\omega_n+1-i/q_n}\to\frac{\omega}{\omega+1}\quad\text{et}\quad\frac{q_n-i}{p_n+q_n-i}=\frac{1-i/q_n}{\omega_n+1-i/q_n}\to\frac1{\omega+1}.</math></center>

Soit ''E'' un ensemble borné, ''A'' un majorant de la borne de ''E'' et ε un réel strictement positif.

*'''Les deux suites (''h''<sub>(''p<sub>n</sub>, q<sub>n</sub>'')</sub>(''x'')) et (''k''<sub>(''p<sub>n</sub>, q<sub>n</sub>'')</sub>(''x'')) convergent uniformément sur ''E'' ainsi que leur rapport<ref name=Pade1899/> :'''<br />Fixons un entier naturel ''m'' tel que<center><math>\sum_{k=m}^{\infty}\frac{A^k}{k!}\le\frac{\varepsilon}3.</math></center>Alors, pour tout entier naturel ''n'' (en utilisant que ''a<sub>k,n</sub>'' ≤ 1/''k''! et ω/(1+ ω) ≤ 1),<center><math>\forall x\in E\quad\left| h_{(p_n,q_n)}(x) - \exp\left(\frac {\omega x}{1 + \omega}\right)\right|\le\sum_{k=0}^{\infty}A^k\left| a_{k,n}-\frac 1{k!}\left(\frac{\omega}{1 + \omega}\right)^k\right| \le R_n+\frac{2\varepsilon}3</math></center>avec<center><math>R_n=\sum_{0\le k<m}A^k\left|a_{k,n}-\frac 1{k!}\left(\frac{\omega}{1+\omega}\right)^k\right|.</math></center>La proposition précédente montre qu'il existe ''N'' tel que pour tout ''n ≥ N'', ''R<sub>n</sub>'' ≤ ε/3. On en déduit :<center><math>\forall n\ge N\quad\forall x\in E\quad\left| h_{(p_n,q_n)}(x) - \exp\left(\frac {\omega x}{1 + \omega}\right)\right|\le\varepsilon.</math></center>On a donc prouvé la première des trois convergences annoncées, uniformément sur ''E''. La deuxième se montre de manière analogue et la troisième s'en déduit, puisque |exp(–''x''/(1+ ω))| est minorée sur ''E'' par exp(–''A''/(1+ ω)) > 0.
{{Démonstration/fin}}

== Notes et références ==
<references/>

== Voir aussi ==
=== Lien externe ===
S. Khémira, [http://tel.archives-ouvertes.fr/docs/00/04/82/06/PDF/tel-00009653.pdf Approximants de Hermite-Padé, déterminants d’interpolation et approximation diophantienne], thèse de doctorat de l'[[université Paris 6]], 2005

=== Bibliographie ===
*{{en}} Claude Brezinski et Michela Redivo-Zaglia, ''Extrapolation Methods: Theory and Practice'', North-Holland, 1991 {{ISBN|978-0-44488814-3}}
*{{en}} George A. Baker et Peter Graves-Morris, ''Padé Approximants'', coll. « Encyclopedia of Mathematics and its Applications » (n° 59), 1996, {{2e}} éd. {{lire en ligne|url=https://books.google.fr/books?id=Vkk4JNLKbLoC}} {{ISBN|978-0-521-45007-2}}

{{Portail|analyse}}

[[Catégorie:Fraction continue]]
[[Catégorie:Analyse complexe]]
[[Catégorie:Leonhard Euler]]