Voir la source de Approximant de Padé

[[Image:Henri Padé.jpeg|vignette|Le concept de l'article doit son nom au [[mathématicien]] [[France|français]] [[Henri Padé]] ([[1863 en science|1863]]-[[1953 en science|1953]]).]]
En [[mathématiques]], et plus précisément en [[analyse complexe]], l''''approximant de Padé''' est une méthode d'[[approximation de fonction|approximation]] d'une [[fonction analytique]] par une [[fonction rationnelle]]. En ce sens, elle est un peu analogue à un [[développement limité]] qui approche la fonction selon les mêmes critères à l'aide d'un [[polynôme]].

De même que les développements limités forment une suite appelée [[série entière]], convergeant vers la fonction initiale, les approximants de Padé apparaissent comme les réduites de diverses [[Fraction continue généralisée|fractions continues (généralisées)]] dont la limite est aussi la fonction initiale. En ce sens, ces approximants font partie de la vaste théorie des [[Fraction continue|fractions continues]].

En [[analyse complexe]], les approximants offrent un développement dont le domaine de convergence est parfois plus large que celui d'une série entière. Ils permettent ainsi de [[Prolongement analytique|prolonger des fonctions analytiques]] et d'étudier certains aspects de la question des [[Série divergente|séries divergentes]]. En [[théorie analytique des nombres]], l'approximant permet de mettre en évidence la nature d'un nombre ou d'une fonction [[arithmétique]] comme celle de la [[fonction zêta de Riemann]]. Dans le domaine du [[calcul numérique]], l'approximant joue un rôle, par exemple, pour évaluer le comportement d'une solution d'un [[système dynamique]] à l'aide de la [[théorie des perturbations]].

Le développement d'une fonction en fraction continue est utilisé pour la première fois par [[Leonhard Euler]], pour démontrer l'[[Nombre irrationnel|irrationalité]] du [[e (nombre)|nombre {{math|e}}]]. Une stratégie plus élaborée permet à [[Jean-Henri Lambert]] de montrer celle de [[pi|{{math|π}}]]. Cette notion est développée plus systématiquement par [[Henri Padé]] et érigée en théorie à part entière.

== Préambule ==
=== Présentation du concept ===
[[Image:Charles Hermite circa 1887.jpg|redresse|vignette|[[Charles Hermite]] utilise ce qui est maintenant appelé un approximant de Padé pour démontrer la transcendance de {{math|e}} en 1873.]]
Il est utile d'approcher une fonction donnée par une [[Suite et série de fonctions|suite de fonctions]] aisément calculables. Cette démarche est à l'origine de la théorie des [[série entière|séries entières]] consistant à approcher de plus en plus précisément une [[fonction analytique]] à l'aide de la suite de ses [[développement limité|développements limités]]. Les développements limités, puis les séries entières offrent de nombreuses possibilités comme le calcul d'une limite ou la résolution d'une [[équation différentielle]].

La démarche de Padé est, à beaucoup d'égards, analogue à celle des séries entières. L'objectif est d'approcher « le mieux possible » une fonction analytique à l'aide d'une [[fonction rationnelle|fraction rationnelle]]. « Le mieux possible » signifie ici que pour un [[Couple (mathématiques)|couple]] (''p'', ''q'') d'entiers positifs donné, la fraction rationnelle ''h''/''k'' est telle que les polynômes ''h'' et ''k'' sont de degrés majorés respectivement par ''p'' et ''q'' et son développement limité « colle » le mieux possible à celui de la fonction ''f''. Si les deux développements limités coïncident jusqu'à l'ordre ''p + q'', une telle fraction est un approximant de Padé de ''f''.

De même qu'il est intéressant de considérer la suite des développements limités, il est judicieux d'étudier les suites d'approximants, de plus en plus précises et convergeant vers la fonction ''f''. Ces suites peuvent toujours s'exprimer sous forme de [[Fraction continue généralisée|fractions continues (généralisées)]]. Dans le cas des polynômes, les développements limités sont indexés par un [[Entier naturel|entier positif]], ce qui engendre une suite naturelle : la série entière. Les approximants de Padé sont indexés par un couple d'entiers positifs, engendrant une infinité plus complexe de suites. D'autres causes rendent l'approche plus complexe. La somme, le produit, la dérivée ou la primitive d'approximant de Padé n'ont aucune raison d'être encore un approximant de Padé.

Ces suites offrent néanmoins quelques avantages que n'ont pas les séries entières. La nature [[arithmétique]] d'un nombre comme [[Pi|{{math|π}}]] ou [[E (nombre)|{{math|e}}]] est plus simplement mise en valeur. Les preuves de l'irrationalité de {{math|π}} et de {{math|e}}, comme celles de leur [[Nombre transcendant|transcendance]], se fondent naturellement sur cet outil<ref name=FCAD/>. Les suites d'approximants de Padé possèdent parfois un autre avantage. Certaines fonctions, comme la [[fonction zêta de Riemann]], sont initialement définies par des séries dont le domaine de convergence est limité par la présence des [[Pôle (mathématiques)|pôles]], c'est-à-dire des espèces de montagnes qui grimpent jusqu'à l'infini. La série entière ne permet pas de voir plus loin qu'un pôle. Les approximants de Padé possèdent le mérite de permettre de voir l'autre versant des pôles. Le caractère parfois plus vaste du domaine de convergence d'une suite d'approximants de Padé par rapport à une expression sous forme de série entière est donc un atout non négligeable. La physique ou l'astronomie offrent des exemples d'application. L'étude d'un [[système dynamique]] un peu complexe peut nécessiter le recours à la [[théorie des perturbations]], les solutions sont souvent des fonctions analytiques exprimées sous forme de séries entières. Or la présence de pôles impose parfois un [[rayon de convergence]] trop petit pour permettre un calcul effectif intéressant. Ce phénomène se traduit par le besoin de calculer la limite d'une [[série divergente]]. Trouver une bonne suite d'approximants de Padé permet d'élargir le domaine de convergence, le comportement de la fonction analytique est ainsi connu sur une plus vaste étendue et le calcul effectif de la limite est possible.

=== Approche par la fraction continue ===
{{Article détaillé|Fraction continue}}

Une [[fraction continue]] permet d'approcher un [[nombre réel]] avec une précision aussi bonne que souhaitée. Une démarche analogue peut être appliquée aux fonctions analytiques. Une fonction analytique ''f'' se développe en série entière ; si ''t''{{ind|0}} est un point de [[Ensemble de définition|domaine de définition]] de ''f'' et ''t'' un scalaire de [[module d'un nombre complexe|module]] strictement inférieur à ''r'', le rayon de convergence de ''f'' au point ''t''{{ind|0}}, il est possible de développer ''f'' de la manière suivante : ''f''(''t''{{ind|0}} + ''t'') = α{{ind|0}} + α{{ind|''n''{{ind|1}}}}''t''{{exp|''n''{{ind|1}}}} + …. Ici, α{{ind|''n''{{ind|1}}}} ainsi que les coefficients suivants, sont choisis non nuls. Pour cette raison, ''n''{{ind|1}} n'est pas nécessairement égal à 1. Il existe une fonction analytique ''f''{{ind|1}} définie en 0 et telle que
<center><math>f(t_0 + t) = a_0 + \frac {t^{n_1}}{f_1(t)}\quad\text{avec}\quad a_0 = \alpha_0.</math></center>

Un traitement analogue sur la fonction ''f''{{ind|1}}, puis ''f''{{ind|2}}, montre qu'il est encore possible d'écrire la fonction ''f'' de la manière suivante :
<center><math>f(t_0 + t) = a_0 + \cfrac {t^{n_1}}{a_1 + \cfrac {t^{n_2}}{f_2(t)}} = a_0 + \cfrac {t^{n_1}}{a_1 + \cfrac {t^{n_2}}{a_2 + \cfrac {t^{n_3}}{f_3(t)}}} = \cdots</math></center>

Cet algorithme, bien que naïf et impliquant des calculs laborieux, permet d'approcher la fonction analytique ''f'' par une fraction rationnelle. Si ''f'' est elle-même une fraction rationnelle, il existe une valeur ''p'' tel que ''f{{ind|p}}'' soit une [[fonction constante]] et le processus s'arrête ; dans le cas contraire, le processus continue indéfiniment. Cette configuration est analogue à celle des [[Nombre rationnel|nombres rationnels]] pour la fraction continue.

On utilise dans cet article la notation de [[Alfred Pringsheim|Pringsheim]] :

<center><math>f(t_0 + t) = a_0 + \frac{t^{n_1}\mid}{\mid a_1} + \frac{t^{n_2}\mid}{\mid a_2} + \frac{t^{n_3}\mid}{\mid a_3} +\dots </math></center>

Si le processus est arrêté à la ''p''{{e}} étape, on obtient une fraction rationnelle qui « approche » la fonction initiale ''f'' et qui est un exemple d'[[#Définitions|approximant de Padé]]<ref>{{Ouvrage|langue=en|nom1=[[Lloyd N. Trefethen]]|titre=Approximation Theory and Approximation Practice|éditeur=[[Society for Industrial and Applied Mathematics|SIAM]]|année=2013|pages totales=295|page=248-249|isbn=978-1-61197-239-9|présentation en ligne=https://books.google.com/books?id=En41UGQ6YXsC&pg=PA248}}, {{Lang|en|Exercise}} 27.7.</ref>.

Pour une raison de simplicité, dans la suite de l'article, la valeur ''t''{{ind|0}} est choisie nulle.

=== Fonction tangente ===
Pour illustrer cette démarche, considérons la fonction [[fonction trigonométrique|fonction tangente]], l'un des premiers exemples historiques de fonction analytique dont le développement en fraction continue est calculé, par Lambert<ref name=FCAD>{{Cf.|maj}} [[Fraction continue et approximation diophantienne]].</ref> :

<center><math>\tan (t) = \cfrac t{1 - \cfrac {t^2}{3 - \frac {t^2} {5 - \frac {t^2}{7 - \cdots}}}} = \frac{t \mid}{\mid 1} - \frac{t^2 \mid}{\mid 3} - \frac{t^2 \mid}{\mid 5} - \frac{t^2 \mid}{\mid 7} + \cdots </math></center>

Appliquons l'algorithme naïf du paragraphe précédent au point ''t''<sub>0</sub> = 0. La valeur ''a''<sub>0</sub> est égale à 0 et ''n''<sub>1</sub> à 1. On obtient une première fraction rationnelle ''h''<sub>0</sub>(''t'')/''k''<sub>0</sub>(''t'') approchant la fonction tangente :
<center><math>\tan(t) = 0 + \frac t{f_1(t)},\quad f_1(t) = \frac t{\tan(t)}\quad\text{et}\quad h_0(t) = 0,\; k_0(t) = 1.</math></center>

Prolongeons le calcul :

<center><math>f_1(t)=\frac t{\tan(t)}=\frac t{t+\frac{t^3}3+O(t^5)}=1-\frac{t^2}{3+O(t^2)}=1-\frac {t^2}{f_2(t)}.</math></center>

On obtient une nouvelle expression plus précise de la fonction tangente :
<center><math>\tan(t) = 0 + \frac{t \mid}{\mid 1} - \frac{t^2 \mid}{\mid f_2(t)},\quad f_2(t) =\frac{t^2\tan t}{\tan t-t}\quad\text{et}\quad h_1(t) = t,\; k_1(t) = 1.</math></center>

Par une autre méthode, Lambert découvre<ref>{{Article|prénom=J.-H.|nom=Lambert|titre=Mémoire sur quelques propriétés remarquables des quantités transcendentes [sic] circulaires et logarithmiques|revue=Mémoires de l'[[Académie royale des sciences de Prusse|Académie royale des sciences et belles-lettres]]|lieu=Berlin|vol=17|année=1761|p.=265-322|url=http://www.bibnum.education.fr/mathematiques/theorie-des-nombres/lambert-et-l-irrationalite-de-p-1761}}, § 21.</ref> une formule de récurrence permettant de calculer plus rapidement ''h<sub>n</sub>''(''t'')/''k<sub>n</sub>''(''t'') :
<center><math>\forall n \ge 2\quad h_n(t) = (2n-1)h_{n-1}(t) -t^2h_{n-2}(t) \quad \text{et}\quad k_n(t) = (2n-1)k_{n-1} -t^2k_{n-2}(t).</math></center>

{{Démonstration/début|titre=Démonstration historique de Lambert}}
Lambert transpose de la manière suivante (aux notations près) l'[[algorithme d'Euclide]] usuellement appliqué au calcul des [[Fraction continue|fractions continues de réels]] : il « [[Division suivant les puissances croissantes|divise]] » d'abord ''R''{{ind|–1}} := cos(''t'') par {{nobr|1=''R''{{ind|0}} := sin(''t''),}} ce qui lui donne un « quotient » {{nobr|1=''Q''{{ind|1}} = 1/''t''}} et un « reste » –''R''{{ind|1}} :
<center><math>\cos t=\frac{\sin t}t-R_1\text{ donc }\tan t=\frac t{1-\frac{tR_1}{R_0}}.</math></center>
Il prévoit qu'indéfiniment, ''R''{{ind|''n''–2}} divisé par ''R''{{ind|''n''–1}} donne un quotient ''Q{{ind|n}}'' = (2''n'' – 1)/''t'' (et un reste –''R{{ind|n}}'') :
<center><math>R_{n-2}=\frac{2n-1}tR_{n-1}-R_n\text{ donc }\frac{tR_{n-1}}{R_{n-2}}=\frac{t^2}{2n-1-\frac{tR_n}{R_{n-1}}}.</math></center>
Cela fournit le développement annoncé, sous réserve de prouver que la [[Valuation#Ordre d'annulation d'un polynôme|valuation]] des ''R{{ind|n}}'', définis par
<center><math>R_{-1}=\cos,\quad R_0=\sin\quad{\rm et}\quad\forall n>0\quad R_n=\frac{2n-1}tR_{n-1}-R_{n-2},</math></center>
croît avec ''n''. Il lui suffit pour cela d'expliciter la valeur des restes qu'il observe pour les premières valeurs de ''n'' et qu'il prévoit pour tout ''n'' :
<center><math>\forall n\in\N\quad R_n=t^{n+1}2^n\sum_{k=0}^{\infty}\frac{(-1)^kt^{2k}(n+k)!}{k!(2n+2k+1)!}</math></center>
puis de vérifier cette égalité par une [[Raisonnement par récurrence#Récurrence d’ordre 2|récurrence d'ordre 2]] : si elle est vraie aux ordres ''n'' – 2 et ''n'' – 1, on trouve bien, en remplaçant ''R''{{ind|''n''–2}} et ''R''{{ind|''n''–1}} par leurs expressions et en regroupant les termes :
<center><math>R_n=\frac{2n-1}tR_{n-1}-R_{n-2}=t^{n-1}2^n\sum_{k=0}^{\infty}\frac{k(-1)^{k-1}t^{2k}(n+k-1)!}{k!(2n+2k-1)!}=t^{n+1}2^n\sum_{k=1}^{\infty}\frac{(-1)^{k-1}t^{2(k-1)}(n+k-1)!}{(k-1)!(2n+2(k-1)+1)!}.</math></center>

La relation de récurrence sur les numérateurs ''h{{ind|n}}'' et les dénominateurs ''k{{ind|n}}'' des fractions partielles est alors une conséquence immédiate de l'[[algorithme d'Euclide étendu]].
{{Démonstration/fin}}

Il démontre de plus la [[#Illustration par l'exemple de la fonction tangente|convergence des quotients ''h{{ind|n}}''/''k{{ind|n}}'' vers la fonction tangente]]<ref>{{Harvsp|Lambert|1761|loc=§ 32}} : {{Citation|Il s'agit de déterminer la loi, suivant laquelle les fractions […] approchent de la valeur de la tangente.}}</ref>. Il remarque même que {{Citation|toutes ces suites sont plus convergentes, que ne l'est aucune progression géométrique décroissante<ref><!--cf. analyse de Juhel sur Bibnum-->{{Harvsp|Lambert|1761}}, § 34.</ref>}}, ce qui est la clé de sa démonstration de l'irrationalité de {{math|π}}<ref name=FCAD/>.

{{Démonstration/début|titre=Une démonstration de convergence dans un cas simple : sur [–1, 1]<ref>La convergence est étudiée sur un domaine plus vaste en fin d'article.</ref>}}
Par parité des ''k{{ind|n}}'' et imparité de tan et des ''h{{ind|n}}'', il suffit de prouver la convergence sur l'intervalle [0, 1].

La méthode utilisée ici consiste à étudier la suite des fonctions ψ<sub>''n''</sub>(''t'') définies sur [0, 1] par
<center><math>\psi_n(t) = k_n(t)\sin t - h_n(t) \cos t,</math></center>
autrement dit, par
<center><math>\psi_0(t)=\sin t,\quad\psi_1(t)=\sin t-t\cos t\quad\text{et}\quad\forall n\ge2\quad\psi_n(t)=(2n-1)\psi_{n-1}(t)-t^2\psi_{n-2}(t).</math></center>
Pour montrer que la suite de ces fonctions tend uniformément vers 0, on déduit d'abord facilement de la définition (par une [[Raisonnement par récurrence#Récurrence d’ordre 2|récurrence d'ordre 2]]) que pour tout entier ''n'' > 0, la dérivée de ψ<sub>''n''</sub>(''t'') est égale à ''t''ψ{{ind|''n''–1}}(''t'') puis, que chaque fonction ψ{{ind|''n''}} est nulle en 0 et strictement croissante (donc positive). On en déduit alors que la suite des ψ<sub>''n''</sub> tend uniformément vers zéro et même (par récurrence, en utilisant que ψ<sub>''n''</sub> est la primitive de ''t''ψ{{ind|''n''–1}}(''t'') qui s'annule en 0) que
<center><math>\forall t\in[0,1]\quad\psi_n(t)\le\frac{t^{2n+1}}{2^nn!}.</math></center>
Pour en déduire que ψ<sub>''n''</sub>/(''k{{ind|n}}''cos) converge aussi vers 0, on a besoin de minorer ''k{{ind|n}}''. Une récurrence montre que pour tout entier ''n'' > 0,
<center><math> \forall t \in [0,1]\quad 2^{n-1}(n-1)!\le k_n(t)\le\frac{(2n)!}{2^nn!}.</math></center>
On dispose donc de l'encadrement suivant :
<center><math> \forall t \in [0,1]\quad 0\le\tan (t) - \frac {h_n(t)}{k_n(t)}  = \frac {\psi_n(t)}{\cos(t)k_n(t)}\le\frac{t^{2n+1}}{2^{2n-1}n!(n-1)!\cos(1)},</math></center>
ce qui montre non seulement que la suite (''h{{ind|n}}''/''k{{ind|n}}'') converge uniformément sur [0, 1] vers la fonction tangente mais aussi que, pour ''n'' fixé, ''h{{ind|n}}''/''k{{ind|n}}'' et tangente ont même développement limité en 0 à l'ordre 2''n''.
{{Démonstration/fin}}

[[Image:Padé Tangente.jpg|vignette|upright=2|Dans cet exemple, la suite des approximants de Padé « avale » les pôles et fournit une approximation sur tous les nombres réels, à l'exception des pôles.]]

Une analyse plus précise montre que si ''t'' n'est pas de la forme {{math|''k''π + π/2}}, où ''k'' est un entier, ''h{{ind|n}}''(''t'')/''k{{ind|n}}''(''t'') tend vers tan(''t'').

Les approximations successives permettent d'« avaler » les pôles de la fonction tangente. Sur les réels positifs, la deuxième approximation, en violet sur la figure, simule le premier pôle avec une [[asymptote]] à {{sqrt|3}} au lieu de {{math|π/2}}. La quatrième approximation, en bleu, ne diffère pas de manière visible de la fonction tangente, en rouge sur le graphique, sur l'intervalle {{math|[0,&nbsp;π]}}. Elle possède deux pôles sur les réels positifs. La huitième approximation, en vert, « colle à la fonction tangente » sur l'intervalle {{math|[0,&nbsp;2π]}} et possède 4 pôles positifs dont 3 sont visibles sur la figure. De manière plus générale, si ''n'' est un entier strictement positif et ε un réel strictement positif, la suite des approximants [[convergence uniforme|converge uniformément]] sur la réunion des intervalles {{math|[0,&nbsp;π/2 – ε]}} et {{math|[(2''j''–1)π/2 – ε, (2''j''+1)π/2 + ε]}} où ''j'' varie de 1 à ''n''.

{{refnec|date=avril 2014|Cette capacité à « avaler les pôles » est l'un des attraits des approximants de Padé,}} d'autant plus fort que ce résultat reste valable si les fonctions sont considérées comme complexes de la variable complexe. À la différence des séries entières, les approximants de Padé fournissent des informations sur la fonction tangente en dehors du disque de rayon {{math|π/2}}. {{référence insuffisante|date=mars 2014|Cette propriété est utilisée pour}} l'étude de la [[fonction zêta de Riemann]]<ref>{{Article|titre=Nombres d'Euler, approximants de Padé et constante de Catalan|revue=Ramanujan J.|vol=11|année=2006|p.=199-214|url=http://www-fourier.ujf-grenoble.fr/~rivoal/articles/gnew6.pdf|prénom=Tanguy|nom=Rivoal}}.</ref>.

=== Fonction exponentielle ===
{{Article détaillé|Approximant de Padé de la fonction exponentielle}}
Il est possible d'appliquer un [[algorithme]] un peu de même manière avec la [[fonction exponentielle]]. On peut par exemple, établir une expression en fraction continue :
<center><math>\exp(x) = 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>

L'algorithme utilisé est ici un peu différent. Les numérateurs ne sont plus des constantes mais des [[Fonction affine|fonctions affines]]. En revanche une propriété reste commune : le [[développement limité]] des fractions ''h''<sub>p</sub>/''k''<sub>p</sub> à l'ordre égal à la somme des degrés du numérateur et du dénominateur est identique à celui de la fonction exponentielle.

Il existe de multiples expressions différentes de la fonction exponentielle sous forme de fraction continue, ce qui amène une série de questions d'ordre général sur les approximations par des fractions rationnelles d'une fonction analytique. Quatre sont particulièrement importantes aux yeux de Padé : pour un couple d'entiers strictement positifs (''p'', ''q''), existe-t-il une fraction rationnelle ''h''/''k'' telle que ''h'' soit un polynôme de degré ''p'', ''k'' un polynôme de degré ''q'' et telle que la fraction ait le même [[développement limité]] à l'ordre ''p'' + ''q'' que la fonction exponentielle ? Existe-t-il des relations de récurrence permettant de passer d'un approximant à un autre à l'aide de formules analogues à celles présentées dans cet exemple ? Ces formules de récurrence permettent-elles d'établir des fractions continues ? Enfin, ces fractions continues convergent-elles vers la fonction cible ?

Ces réponses, toutes positives pour la fonction exponentielle, sont l'objet de l'article détaillé. Cet exemple est choisi par Padé<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|vol=16|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> pour introduire sa théorie.

== Fragments d'histoire ==
=== Origine ===
[[Image:Joseph Louis Lagrange.jpg|vignette|[[Joseph-Louis Lagrange|Lagrange]] énonce une propriété sur des réduites de fractions continues correspondant à notre définition moderne d'approximant de Padé.]]
Padé place<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]), {{p.|38}}.</ref> l'origine de ses approximants aux travaux d'Euler <small>(1707-1783)</small> avec sa démonstration de l'irrationalité de {{math|e}}. Il est vrai qu'[[Formule de fraction continue d'Euler#Fonction exponentielle|une formule d'Euler permet d'exprimer la fonction exponentielle]] sous forme de fraction continue, mais d'une forme très particulière : chaque réduite est non seulement une fraction rationnelle mais un polynôme (la somme partielle de la [[Exponentielle complexe|série entière {{math|exp}}]]). Lambert <small>(1707-1777)</small> franchit un pas fondamental en développant des fonctions<!--cf. commentaire de Juhel sur Bibnum-->, ce qui lui permet de démontrer que si un nombre réel non nul ''t'' est rationnel, alors ni [[Fonction trigonométrique|tan]](''t''), ni exp(''t'') ne l'est<ref name=FCAD/>, démontrant par là l'irrationalité de {{math|π}}. Il développe en fraction continue non seulement tan(''t'') mais [[Tangente hyperbolique|tanh]](''t'')<!--1761-->, [[arctangente|arctan]](''t'')<!--http://librarum.org/book/10065/144-->, [[Logarithme naturel|log]](1 + ''t'')<!--http://librarum.org/book/10065/144-->{{etc.}} L'application de son approche reste néanmoins limitée à la recherche d'expressions sous forme de fraction continue généralisée de nombres spécifiques ; Lambert ne cherche pas ''a priori'' à mieux connaître le comportement d'une fonction.

[[Joseph-Louis Lagrange]] <small>([[1736]]-[[1813 en science|1813]])</small> généralise l'emploi des fractions continues à l'[[Analyse (mathématiques)|analyse]]. Son optique est plus proche de celle de l'article. La fraction continue d'une fonction est utilisée comme un outil de résolution d'[[équation différentielle|équations différentielles]] ; il est nécessaire de réduire la fraction continue qui s'exprime alors comme une fonction rationnelle et Lagrange remarque que dans l'un de ses exemples, les développements limités de la solution et de la fraction rationnelle coïncident {{Citation|jusqu'à la puissance de ''x'' inclusivement qui est la somme des deux plus hautes puissances de ''x'' dans le numérateur et dans le dénominateur<ref>J. L. Lagrange, ''[https://gallica.bnf.fr/ark:/12148/bpt6k229223s/f334.image Sur l'usage des fractions continues dans le calcul intégral]'', Nouveaux mémoires de l'Académie royale des sciences et Belles lettres de Berlin, 1776.</ref>}}, ce qui correspond à une définition moderne des approximants de Padé. [[Carl Friedrich Gauss]] <small>([[1777 en science|1777]]-[[1855 en science|1855]])</small> utilise les résultats de Lagrange et généralise l'usage des fractions continues au plan complexe ; les [[Fraction continue de Gauss|fractions continues de Gauss]] permettent de retrouver des développements en fraction continue de nombreux quotients de [[série hypergéométrique|fonctions hypergéométriques]]. L'approche reste néanmoins empirique. Rien ne permet de déterminer de manière exhaustive toutes les fractions continues associées à une fonction analytique donnée ou même encore d'ébaucher une classification des fractions continues.

=== Formalisation ===
[[Image:Georg Friedrich Bernhard Riemann.jpeg|vignette|gauche|[[Bernhard Riemann|Riemann]] met en évidence un des intérêts des approximants de Padé : ils permettent parfois de prolonger des [[analyse complexe|fonctions analytiques complexes]].]]
Une triple motivation modifie à la fois l'approche et les questions relatives à ces fractions continues un peu spéciales. La première est la conséquence des travaux de [[Bernhard Riemann]] <small>([[1826 en science|1826]]-[[1866 en science|1866]])</small>. Elles mettent en évidence l'importance de l'[[analyse complexe]], dans cet univers, les fonctions sont souvent définies par des séries entières. Or ces séries ne sont pas convergentes à l'extérieur d'un disque donné, la question du prolongement est parfois cruciale. Les approximants de Padé n'ont pas toujours cette limitation, comme le montre l'exemple de la tangente. Leur étude est naturellement {{citation|devenue à l'ordre du jour, de la légitimité de l'emploi, dans le calcul, des séries divergentes<ref>{{Harvsp|Padé|1899|p=395}}.</ref>}}. Riemann étudie la convergence des [[Fraction continue de Gauss|fractions continues de Gauss]]. {{Refnec|date=mars 2014|Cette convergence est à la base d'une}} démonstration du [[théorème d'Apéry]], prouvé en [[1977 en science|1977]] et énonçant que l'image de 3 par la [[fonction zêta de Riemann]] est un nombre irrationnel<ref>{{Article|prénom=T.|nom=Rivoal|titre=Séries hypergéométriques et irrationalité des valeurs de la fonction zêta de Riemann|revue=Journal de théorie des nombres de Bordeaux|vol=15|numéro=1|année=2003|p.=351-365|url=http://www-fourier.ujf-grenoble.fr/~rivoal/articles/zetaja.pdf}}.</ref>.

[[Charles Hermite]], le [[directeur de thèse]] de Padé<ref name=MacTutor>Voir par exemple {{MacTutor|id=Pade|title=Henri Eugène Padé}}.</ref>, étudie la question de l'interpolation d'une fonction par une fraction rationnelle ainsi que les propriétés de fractions continues associées à l'exponentielle, ce qui lui permet de montrer la transcendance de {{math|e}}<ref>{{Article|prénom=C.|nom=Hermite|titre=Sur la fonction exponentielle|lien périodique=CRAS|revue=Comptes rendus hebdomadaires des séances de l'Académie des sciences|p.=18-24|année=1873|vol=77|url=https://gallica.bnf.fr/ark:/12148/bpt6k3034n/f18}}, [http://www.bibnum.education.fr/mathematiques/theorie-des-nombres/la-demonstration-de-la-transcendance-de-e présenté et analysé] par [[Michel Waldschmidt]] sur le site [[Bibnum]].</ref>, le premier nombre démontré transcendant qui n'a pas été spécialement conçu pour une démonstration (comme ceux de [[nombre de Liouville|Liouville]]). Sur la base de ces idées, [[Ferdinand von Lindemann]] démontre la transcendance de {{math|π}} en [[1882 en science|1882]] clôturant ainsi la question millénaire de la [[quadrature du cercle]]. Ce problème, que Charles Hermite croyait hors de portée et dont il disait : {{Citation|Je ne me hasarderai point à la recherche d’une démonstration de la transcendance du nombre {{math|π}}. Que d’autres tentent l’entreprise […] mais croyez-m’en, mon cher ami, il ne laissera pas que de leur en coûter quelques efforts<ref>{{Article|revue=[[Journal für die reine und angewandte Mathematik|J. reine angew. Math.]]|titre=Extrait d'une lettre de Mr. ''Ch. Hermite'' à Mr. ''[[Karl Wilhelm Borchardt|Borchardt]]''|vol=76|année=1873|p.=342-344|url=http://www.digizeitschriften.de/dms/img/?PPN=PPN243919689_0076&DMDID=dmdlog27}}. Voir aussi C. Brezinski, ''Histoires de sciences : Inventions, découvertes et savants'', [[Éditions L'Harmattan|L'Harmattan]], 2006 {{ISBN|978-2-29600350-7}}, [https://books.google.fr/books?id=TZAIalOhq3YC&pg=PA14 {{p.|14}}].</ref>}} était résoluble avec les techniques qu'il avait développées. Enfin, l'étude de systèmes dynamiques complexes comme ceux étudiés par [[Henri Poincaré]] pour démontrer la statibilité du [[système solaire]] impose l'usage de fractions continues de type Padé<ref>Voir par exemple H. Poincaré, ''Méthodes nouvelles de la mécanique céleste'', 3 vol., Gauthier-Villars, Paris, 1892-1899.</ref>.

Il est donc naturel que la question des fractions continues, construites à l'aide de fonctions rationnelles, soient à l'ordre du jour à la fin du {{XIXe siècle}}. Cependant, en [[1890]], Poincaré décrit la théorie des fractions continues algébriques comme une {{citation|sorte de ''terra incognita''}} dont {{Citation|la carte est encore presque blanche<ref>H. Poincaré, « Notice sur Halphen », ''Journal de l'École Polytechnique'', {{60e|cahier}}, 1890, {{p.|137-161}}.</ref>}}. Henri Padé étudie la question sous un jour nouveau. Il cherche systématiquement la fraction rationnelle dont le numérateur et le dénominateur ont des degrés respectifs ne dépassant pas les valeurs d'un couple (''p'', ''q'') et qui approxime localement le mieux une fonction analytique donnée. Son travail de thèse<ref>{{Harvsp|Padé|1892}}.</ref>, soutenu à l'université de la [[Sorbonne]] en 1892, consiste à étudier la table de ces approximants et à en établir une théorie très générale. Il détermine les lois qui permettent de passer de cette table aux fractions continues, dont il parvient à classifier les représentants importants. En [[1898]], il applique sa théorie à la fonction exponentielle<ref name=Pade1899/>, indique comment construire toutes les fractions continues régulières de cette fonction. Avant son mémoire, seuls cinq exemples étaient connus ; son travail montre qu'ils ne correspondent qu'à des cas particuliers d'une triple infinité de fractions.

Une question reste encore largement ouverte : la [[Problème de convergence|convergence des différentes fractions continues]] que l'on peut construire. En [[1894 en science|1894]], [[Thomas-Joannes Stieltjes]] est d'un apport considérable sur cette question<ref>{{Article|prénom=T.-J.|nom=Stieltjes|titre=Recherches sur les fractions continues|lien périodique=Annales de la Faculté des sciences de Toulouse|revue=Annales de la Faculté des sciences de Toulouse|année=1894|url=http://www.numdam.org/numdam-bin/item?id=AFST_1894_1_8_4_J1_0|série=6|vol=8|numéro=4|p.=J1-J122}}.</ref>. Il établit exactement, pour une vaste famille de fractions continues, le domaine de convergence<ref>{{Article|prénom=T.-J.|nom=Stieltjes|titre=Recherches sur les fractions continues [suite et fin]|revue=Annales de la Faculté des sciences de Toulouse|année=1895|url=http://archive.numdam.org/ARCHIVE/AFST/AFST_1995_6_4_4/AFST_1995_6_4_4_A5_0/AFST_1995_6_4_4_A5_0.pdf|vol=9|p.=A5-A47}}.</ref> ; ce sont celles qui s'avèrent nécessaires pour démontrer la stabilité du système solaire étudiée par Poincaré. Cette question de la convergence est jugée suffisamment importante pour faire l'objet du ''Grand prix de l'Académie des sciences de Paris'' de [[1906 en science|1906]]<ref name=MacTutor/>. Son jury est composé par des grands noms de l'époque comme Poincaré ou [[Émile Picard|Picard]]. Il est remporté par Padé qui dit présenter<ref>{{Article|prénom=H.|nom=Padé|titre=Recherches sur la convergence des développements en fractions continues d'une certaine catégorie de fonction|revue=ASENS|série=3|vol=24|année=1907|url=http://archive.numdam.org/ARCHIVE/ASENS/ASENS_1907_3_24_/ASENS_1907_3_24__341_0/ASENS_1907_3_24__341_0.pdf|p.=341-400}}.</ref> {{Citation|une sorte de synthèse et très large généralisation des résultats […] de [[Ludwig Wilhelm Thomé|Thomé]], [[Edmond Laguerre|Laguerre]] et [[Robert de Montessus de Ballore|M. de Montessus]]}}, tout en {{Citation|ne parvenant pas à achever […] la transformation et la généralisation des résultats […] de Stieltjes, [… ayant interrompu ces recherches à la suite de] la publication des beaux résultats obtenus dans cette voie par M. [[Edward Burr Van Vleck|Van Vleck]]<ref>{{Article|lang=en|prénom=Edward Burr|nom=Van Vleck|titre=On an extension of the 1894 Memoir of Stieltjes|année=1903|revue=[[Transactions of the American Mathematical Society|Trans. Amer. Math. Soc.]]|vol=4|p.=297-332}}.</ref>.}}

== Définitions et premières propriétés ==
=== Définitions ===
Dans toute la suite de l'article, ''f''(''t'') désigne une [[fonction analytique]] en 0 et qui ne s'annule pas en 0. Le fait que cette fonction ne s'annule pas en 0 ne limite en rien la généralité des propos. Si ''g''(''t'') est une fonction analytique en 0, il existe nécessairement une valeur ''n'' telle que ''g''(''t'') soit égal à ''t<sup>n</sup>f''(''t'') où ''f'' est une fonction qui satisfait les hypothèses précédentes. Les lettres ''p'' et ''q'' désignent deux entiers positifs.

Pour Padé, la première question associée à sa théorie est : Existe-t-il une fraction rationnelle ''h''/''k'' telle que ''h'' soit un polynôme de degré ''p'', ''k'' un polynôme de degré ''q'' et telle que la fraction ait le même [[développement limité]] à l'ordre ''p'' + ''q'' que ''f'' ? Elle amène la définition suivante<ref name=GHL/> :

* Un '''approximant de Padé d'indice (''p'', ''q'') de la fonction ''f''''' désigne une fraction rationnelle ''h''/''k'' telle que les degrés des polynômes ''h'' et ''k'' soient inférieurs ou égaux respectivement à ''p'' et ''q'' et que ''h''/''k'' approxime ''f'' en 0 au moins à l'ordre ''p + q''.

Si cette définition est souvent reprise, elle n'est pas totalement satisfaisante. Pour s'en rendre compte le plus simple est de considérer le cas où ''f''(''t'') est égal à 1&nbsp;+&nbsp;''t''<sup>2</sup>. Recherchons l'approximation d'indice (1,1). Si ''a'' + ''bt'' désigne le dénominateur, son produit avec la fonction ''f''(''t'') ne doit pas comporter de terme du deuxième degré, ce qui impose à ''a'' d'être nul. Le numérateur doit alors être égal, lui aussi, à ''bt''. Cependant, la fonction constante ''bt''/''bt'' = 1 n'approxime 1 + t{{2}} qu'à l'ordre 1. Cet exemple montre qu'il n'existe pas d'approximant de Padé d'indice (1,1). Pour cette raison, une deuxième définition s'avère nécessaire. On n'espère plus une fraction qui approche la fonction ''f'' à l'ordre ''p'' + ''q'' et l'on se contente de choisir, parmi les fractions dont les degrés des numérateur et dénominateur sont majorés respectivement par ''p'' et ''q'', une fraction qui approxime ''f'' le mieux possible :

* Les termes de '''fraction réduite''' ou '''réduite d'indice (''p'',''q'') de la fonction ''f''''' désignent une fraction rationnelle ''h''/''k'' telle que les degrés des polynômes ''h'' et ''k'' soient inférieurs ou égaux respectivement à ''p'' et ''q'' et telle que pour tous polynômes ''u'' et ''v'' de degrés inférieurs ou égaux respectivement à ''p'' et ''q'', le développement limité de ''u''/''v'' coïncide avec celui de ''f'' à un ordre moins grand que ne le fait celui de ''h''/''k''. La réduite d'indice (''p'',''q'') est parfois notée ''f''<sub>[p,q]</sub>.

=== Premières propriétés ===
La définition de réduite, moins restrictive que celle d'approximant de Padé pallie l'absence de réponse positive à la première question : par définition même, son existence est garantie. Un résultat essentiel à la théorie affirme non seulement l'existence d'une réduite ''h''/''k'' mais aussi son unicité, sous une certaine forme : on peut toujours rendre ''h'' et ''k'' [[Primalité dans un anneau|premiers entre eux]], puis rendre égal à 1 le terme constant de ''k'' (qui est non nul puisque 0 n'est pas un [[Fraction rationnelle#Construction algébrique|pôle]])<ref>Cette normalisation de la réduite correspond au choix de {{Harvsp|Padé|1892}}.</ref> :
{{Théorème|1=Pour tout couple (''p'', ''q''), la réduite d'indice (''p'', ''q'') de ''f'' est l'unique fraction rationnelle [[fraction irréductible|irréductible]] ''h''/''k'' telle que [[Degré d'un polynôme|deg]](''h'') ≤ ''p'', deg(''k'') ≤ ''q'' et ''h''/''k'' approxime ''f'' à l'ordre (au moins) max(''p'' + deg(''k''), ''q'' + deg(''h'')).}}

{{Démonstration/début|titre=Démonstrations}}
* '''Il existe des polynômes ''h'' et ''k'' premiers entre eux tels que deg(''h'') ≤ ''p'', deg(''k'') ≤ ''q'' et ''h''/''k'' approxime ''f'' à l'ordre max(''p'' + deg(''k''), ''q'' + deg(''h'')) :'''<br />Considérons l'application qui à tout polynôme ''g''(''t'') de degré inférieur ou égal à ''q'' associe la partie de degrés ''p'' + 1 à ''p + q'' du développement en série entière de ''g''(''t'')''f''(''t''). C'est une [[application linéaire]] d'un [[espace vectoriel]] de dimension ''q'' + 1 dans un espace de dimension ''q''. Son [[Application linéaire#Noyau et image|noyau]] est donc non réduit au vecteur nul. Soit ''t''<sup>ω</sup>''k''(''t'') un polynôme non nul de cet espace, avec ω le plus grand possible et ''k''(''t'') de terme constant égal à 1.<br />Soit ''h''<sub>1</sub>(''t'') le développement à l'ordre ''p'' de ''t''<sup>ω</sup>''k''(''t'')''f''(''t''). Le développement à l'ordre ''p'' + ''q'' de ''t''<sup>ω</sup>''k''(''t'')''f''(''t'') – ''h''<sub>1</sub>(''t'') est nul par construction, ce qui montre que ''h''<sub>1</sub>(''t'') est un multiple de ''t''<sup>ω</sup>. Soit ''h''(''t'') le polynôme ''h''<sub>1</sub>(''t'') divisé par ''t''<sup>ω</sup>. La fraction rationnelle ''h''(''t'')/''k''(''t'') possède un développement limité à l'ordre ''p'' + ''q'' – ω égal à celui de ''f''(''t''). De plus, ''h''(''t'') et ''k''(''t'') sont premiers entre eux car pour tout facteur commun ''d''(''t''), de degré ω', le polynôme ''t''<sup>ω+ω'</sup>(''k''(''t'')/''d''(''t'')) appartient au noyau ci-dessus donc ω' = 0.
* '''Pour un tel (''h'', ''k''), toute réduite est égale à ''h''/''k'' :'''<br />Soit ''u''/''v'' une réduite d'indice (''p'', ''q''). Elle approxime ''f'' au moins aussi bien que ne le fait ''h''/''k''. Par conséquent, le numérateur de la fraction rationnelle suivante admet 0 comme racine multiple d'ordre au moins ''p'' + ''q'' – ω : ''h''/''k'' – ''u''/''v'' = (''hv – ku'')/''kv''. Ce numérateur étant de degré inférieur ou égal à ''p'' + ''q'' – ω, il est donc nul, ce qui prouve que ''u''/''v'' = ''h''/''k''.
{{Démonstration/fin}}

Les analogies avec la [[fraction continue]] (comme le premier des trois corollaires de la section suivante) sont multiples, ce qui justifie un vocabulaire commun.

=== Table de Padé ===
{{Voir|contenu=Article détaillé : {{Lien|langue=en|trad=Padé table|fr=Table de Padé}}.}}
Une méthode de présentation des approximants est la '''table de Padé'''. Elle consiste en un tableau à double entrée dont la case de ''colonne p'' et de ''ligne q'' contient la réduite d'indice (''p'', ''q''). L'article [[Approximant de Padé de la fonction exponentielle]] propose celle de l'exponentielle, où chaque case contient une fraction rationnelle distincte. Tel n'est [[trivial]]ement pas le cas pour une fonction ''f'' [[Parité d'une fonction|paire]], puisque sa table n'est qu'un dédoublement des lignes et colonnes de celle de la fonction ''g'' telle que ''f''(''t'') = ''g''(''t''{{2}}). Voici par exemple<ref>Cette table figure, sans commentaire, dans {{Harvsp|Padé|1892|p=15}}.</ref> la table de la fonction [[arctan]](''t'')/''t'' :
<center>
{| class="wikitable"
|+Table de Padé pour <math>\arctan(t)/t</math>
|
{| border="1"
| <center><math>\phantom{p}\quad\begin{matrix} \\ q \end{matrix}\quad \!p\,</math></center> ||<center><math> 0 </math></center>||<center><math> 1 </math></center>||<center><math> 2 </math></center>||<center><math> 3 </math></center> ||<center><math> 4 </math></center>||<center><math> 5 </math></center>
|-
| <center><math> 0 </math></center>
||<center><math>\quad{\color{Blue} 1}\quad</math></center>
||<center><math>\quad{\color{Blue} 1}\quad</math></center>
|| <center><math>{\color{Red} 1-\frac13t^2}</math></center>
|| <center><math>{\color{Red} 1-\frac13t^2}</math></center>
|| <center><math>{\color{PineGreen} 1-\frac13t^2+\frac15t^4}</math></center>
|| <center><math>{\color{PineGreen} 1-\frac13t^2+\frac15t^4}</math></center>
|-
|<center><math> 1 </math></center>
|| <center><math>\quad{\color{Blue} 1}\quad</math></center>
|| <center><math>\quad{\color{Blue} 1}\quad</math></center>
|| <center><math>{\color{Red} 1-\frac13t^2}</math></center>
|| <center><math>{\color{Red} 1-\frac13t^2}</math></center>
|| <center><math>{\color{PineGreen} 1-\frac13t^2+\frac15t^4}</math></center>
|| <center><math>{\color{PineGreen} 1-\frac13t^2+\frac15t^4}</math></center>
|-
|<center><math> 2 </math></center>
|| <center><math>{\color{Maroon} \frac1{1+\frac13t^2}}</math></center>
|| <center><math>{\color{Maroon} \frac1{1+\frac13t^2}}</math></center>
|| <center><math>{\color{Plum} \frac{1+\frac4{15}t^2}{1+\frac9{15}t^2}}</math></center>
|| <center><math>{\color{Plum} \frac{1+\frac4{15}t^2}{1+\frac9{15}t^2}}</math></center>
|| <center><math>{\color{BrickRed} \frac{1+\frac8{21}t^2+\frac4{105}t^4}{1+\frac57t^2}}</math></center>
|| <center><math>{\color{BrickRed} \frac{1+\frac8{21}t^2+\frac4{105}t^4}{1+\frac57t^2}}</math></center>
|-
|<center><math> 3 </math></center>
|| <center><math>{\color{Maroon} \frac1{1+\frac13t^2}}</math></center>
|| <center><math>{\color{Maroon} \frac1{1+\frac13t^2}}</math></center>
|| <center><math>{\color{Plum} \frac{1+\frac4{15}t^2}{1+\frac9{15}t^2}}</math></center>
|| <center><math>{\color{Plum} \frac{1+\frac4{15}t^2}{1+\frac9{15}t^2}}</math></center>
|| <center><math>{\color{BrickRed} \frac{1+\frac8{21}t^2+\frac4{105}t^4}{1+\frac57t^2}}</math></center>
|| <center><math>{\color{BrickRed} \frac{1+\frac8{21}t^2+\frac4{105}t^4}{1+\frac57t^2}}</math></center>
|-
|<center><math> 4 </math></center>
|| <center><math>{\color{ForestGreen} \frac1{1+\frac13t^2-\frac4{45}t^4}}</math></center>
|| <center><math>{\color{ForestGreen} \frac1{1+\frac13t^2-\frac4{45}t^4}}</math></center>
|| <center><math>{\color{Brown} \frac{1+\frac{11}{21}t^2}{1+\frac67t^2+\frac3{35}t^4}}</math></center>
|| <center><math>{\color{Brown} \frac{1+\frac{11}{21}t^2}{1+\frac67t^2+\frac3{35}t^4}}</math></center>
|| <center><math>{\color{Violet} \frac{1+\frac79t^2+\frac{64}{945}t^4}{1+\frac{10}9t^2+\frac5{21}t^4}}</math></center>
|| <center><math>{\color{Violet} \frac{1+\frac79t^2+\frac{64}{945}t^4}{1+\frac{10}9t^2+\frac5{21}t^4}}</math></center>
|-
|<center><math> 5 </math></center>
|| <center><math>{\color{ForestGreen} \frac1{1+\frac13t^2-\frac4{45}t^4}}</math></center>
|| <center><math>{\color{ForestGreen} \frac1{1+\frac13t^2-\frac4{45}t^4}}</math></center>
|| <center><math>{\color{Brown} \frac{1+\frac{11}{21}t^2}{1+\frac67t^2+\frac3{35}t^4}}</math></center>
|| <center><math>{\color{Brown} \frac{1+\frac{11}{21}t^2}{1+\frac67t^2+\frac3{35}t^4}}</math></center>
|| <center><math>{\color{Violet} \frac{1+\frac79t^2+\frac{64}{945}t^4}{1+\frac{10}9t^2+\frac5{21}t^4}}</math></center>
|| <center><math>{\color{Violet} \frac{1+\frac79t^2+\frac{64}{945}t^4}{1+\frac{10}9t^2+\frac5{21}t^4}}</math></center>
|}
|}
</center>

Ici, les mêmes fractions sont affichées à l'aide d'une couleur commune. On remarque qu'elles couvrent des surfaces carrées. On peut préciser ces carrés :

{{Théorème|Soient ''h''/''k'' une fraction rationnelle irréductible, ''p'' le degré de ''h'', ''q'' le degré de ''k'', et ''p + q'' + ω l'ordre exact auquel ''h''/''k'' approxime ''f''. Les indices pour lesquels ''h''/''k'' est une réduite de ''f'' sont les couples (''p + i'', ''q + j'') pour 0 ≤ ''i'' ≤ ω et 0 ≤ ''j'' ≤ ω<ref>{{Harvsp|Padé|1892|p=30}}.</ref>.}}

En effet, d'après le théorème de la section précédente, les indices (''m'', ''n'') pour lesquels ''h''/''k'' est une réduite sont ceux pour lesquels ''p'' ≤ ''m'', ''q'' ≤ ''n'' et max(''m + p'', ''n + q'') ≤ ''p + q'' + ω.

On en déduit le cas particulier ω infini, ainsi que — lorsque l'ensemble de couples ci-dessus est non vide, c'est-à-dire ω ≥ 0 — le lien entre les notions d'approximant de Padé et de réduite :

{{Théorème|Corollaires|2=Soient ''h'', ''k'', ''p'', ''q'' et ω comme ci-dessus.
*La fonction ''f'' coïncide avec ''h''/''k'' sur un voisinage de 0 si et seulement si, pour tous ''m ≥ p'' et ''n ≥ q'', sa  réduite d'indice (''m'', ''n'') est égale à ''h''/''k''.
*Si ''h''/''k'' est un approximant de Padé d'indice (''m'', ''n'') alors il est égal à la réduite d'indice (''m'', ''n'').
*Si ''h''/''k'' fait partie des réduites de ''f'' alors, les indices pour lesquels c'est même un approximant de Padé sont les couples (''p + i'', ''q + j'') pour ''i'', ''j'' ≥ 0 et ''i + j'' ≤ ω.
}}

== Fraction continue ==
=== Généralités ===
L'approximant de Padé est particulièrement utile au sein d'une suite de réduites de plus en plus avancées. Une réduite est dite '''plus avancée''' qu'une autre lorsque la somme des deux coordonnées du centre du carré (constitué de cases de réduites identiques) où se situe la première est strictement plus élevée que celle de l'autre, autrement dit — d'après le théorème précédent — lorsqu'elle constitue (au voisinage de 0) une meilleure approximation<ref>{{Harvsp|Padé|1892|p=31}}.</ref>. Ces suites représentent un peu l'équivalent d'une série entière. Une série entière peut être vue comme une [[suite de polynômes]] qui approxime localement de mieux en mieux la fonction cible, à l'image des suites d'approximants décrites dans ce paragraphe. Une suite de cette nature est toujours la suite des réduites d'une [[fraction continue généralisée]] :

{{Énoncé|Soit (''h{{ind|n}}''/''k{{ind|n}}'') une suite de [[Fraction rationnelle|fractions rationnelles]] dont la première est non nulle et chacune est distincte de la suivante. Il existe deux suites de fractions rationnelles (α<sub>''n''</sub>) et (β<sub>''n''</sub>) telles que les ''h{{ind|n}}'' et ''k{{ind|n}}'' soient exactement les [[Fraction continue généralisée|numérateurs et dénominateurs des réduites]] de<center><math>(I)\quad\beta_0+\frac{\alpha_1\mid}{\mid\beta_1}+\frac{\alpha_2\mid}{\mid\beta_2}+\cdots\quad\text{ou de}\quad(II)\quad\frac{\alpha_0\mid}{\mid\beta_0}+\frac{\alpha_1\mid}{\mid\beta_1}+\cdots.</math></center>
}}

Il suffit en effet de définir (α<sub>''n''</sub>) et (β<sub>''n''</sub>) par<ref>{{Harvsp|Padé|1892|p=47}}.</ref> la [[Suite définie par récurrence|relation de récurrence]] :
<center><math>
\alpha_{n+2}=-\frac{h_{n+1}k_{n+2}-h_{n+2}k_{n+1}}{h_nk_{n+1}-h_{n+1}k_n},\quad
\beta_{n+2}=\frac{h_nk_{n+2}-h_{n+2}k_n}{h_nk_{n+1}-h_{n+1}k_n}</math></center>
et l'initialisation (en supposant, dans le cas {{math|(''I'')}}, que ''k''{{ind|0}} = 1) :
<center><math>(I)\quad\beta_0=h_0,\quad\alpha_1=h_1-h_0k_1,\quad\beta_1=k_1,</math></center>
<center><math>(II)\quad\alpha_0=h_0,\quad\beta_0=k_0,\quad\alpha_1=(h_0k_1-h_1k_0)/h_0,\quad\beta_1=h_1/h_0.</math></center>

Les (α<sub>''n''</sub>) et (β<sub>''n''</sub>) ainsi construites sont en général seulement des suites de fractions rationnelles. On peut bien sûr, dans ce cas, les [[Fraction continue généralisée#Conversions|convertir (sans modifier les réduites)]] en des suites de polynômes, mais l'inconvénient est que les nouveaux numérateurs ''A{{ind|n}}'' et dénominateurs ''B{{ind|n}}'' des réduites ne coïncideront alors plus avec les anciens, ''h{{ind|n}}'' et ''k{{ind|n}}'' : on aura seulement égalité des fractions rationnelles, ''A{{ind|n}}''/''B{{ind|n}} = h{{ind|n}}''/''k{{ind|n}}'', les nouveaux numérateur et dénominateur pouvant être multipliés par un polynôme parasite. Il complexifie inutilement l'expression de la réduite et, pour un calcul effectif, ajoute une instabilité par l'adjonction de singularités factices.

=== Fraction continue simple ===
Pour éviter l'apparition ''par conversion'' des polynômes parasites décrits au paragraphe précédent, il faut et il suffit que les fractions rationnelles α{{ind|''n''}} et β{{ind|''n''}}, définies comme ci-dessus à partir des ''h{{ind|n}}'' et ''k{{ind|n}}'', soient ''déjà'' des polynômes. Padé limite l'étude des suites (''h{{ind|n}}''/''k{{ind|n}}'') pour lesquelles cette condition est vérifiée, en ajoutant des contraintes :

* Une fraction continue (de type {{math|(''I'')}} ou {{math|(''II'')}} ci-dessus) est dite '''simple''' lorsque ses numérateurs partiels α<sub>''n''</sub> sont des monômes non constants, sauf α{{ind|0}} (pour le type {{math|(''II'')}}) qui doit être une constante non nulle et lorsque de plus, chacun de ses dénominateurs partiels β<sub>''n''</sub> est un polynôme de terme constant non nul<ref>{{Harvsp|Padé|1892|p=48}}.</ref>.

Cette condition est un peu équivalente à celle qui sépare les fractions continues de nombres réels, dont les numérateurs partiels sont toujours égaux à 1, d'avec les fractions continues généralisées qui ont des numérateurs partiels quelconques. À l'image de la situation des fractions continues de nombres réels, le numérateur et le dénominateur de chaque réduite sont alors [[Primalité dans un anneau|premiers entre eux]], d'après la caractérisation suivante :

{{Énoncé|1=Une suite de couples de polynômes (''h{{ind|n}}'', ''k{{ind|n}}''), est la suite des numérateurs et dénominateurs des réduites d'une fraction continue simple de type {{math|(''I'')}} (resp. {{math|(''II'')}}) si et seulement si :
#''h''{{ind|0}} et tous les ''k{{ind|n}}'' ont un terme constant non nul ;
#les ''h{{ind|n}}k''{{ind|''n''+1}} – ''h''{{ind|''n''+1}}''k{{ind|n}}'' sont des monômes non nuls ;
#la suite de leurs degrés est strictement croissante ;
#''k''{{ind|0}} = 1 (resp. ''h''{{ind|0}} est constant).
}}

Cette propriété générale indique comment choisir une suite d'approximants de Padé ''f''<sub>[''p<sub>n</sub>'', ''q<sub>n</sub>'']</sub> pour que la fraction continue associée soit simple. En effet, la condition 1 est toujours vérifiée, la 4 l'est si le premier approximant est au bord du tableau — c'est-à-dire ''p''{{ind|0}} ou ''q''{{ind|0}} nul, selon le type — et, sous réserve que la condition 2 soit assurée par d'autres moyens, la propriété suivante garantit la condition 3 :

{{Énoncé|Soit une suite d'approximants de Padé ''f''<sub>[''p<sub>n</sub>'', ''q<sub>n</sub>'']</sub> tels que pour tout entier naturel ''n'', le polynôme suivant soit un monôme non nul :<center><math>h_{p_n,q_n}k_{p_{n+1},q_{n+1}}-h_{p_{n+1},q_{n+1}}k_{p_n,q_n}.</math></center>Pour que la suite des degrés de ces monômes soit strictement croissante, il suffit que les deux suites (''p{{ind|n}}'') ''et'' (''q{{ind|n}}'') soient croissantes.
}}

[[Image:Padé fraction continue simple.jpg|vignette]]

Il reste encore à trouver un moyen pratique pour assurer la condition 2. La connaissance de la structure de la table de Padé permet de répondre à ce besoin. Elle est composée de carrés ne contenant que des réduites égales. Soient ''p'', ''q'' et ω trois entiers positifs tels que les couples (''p'', ''q''), (''p''+ω, ''q''), (''p'', ''q''+ω) et {{nobr|(''p''+ω, ''q''+ω)}} forment les sommets du carré contenant toutes les occurrences de l'approximant de Padé ''f''<sub>[p, q]</sub>, comme indiqué à la fin du [[#Table de Padé|§ « Table de Padé »]]. La propriété suivante permet de trouver les candidats pour la construction d'une fraction continue simple :

{{Énoncé|
*Les deux polynômes suivants sont des monômes non nuls :<center><math>h_{p,q}k_{p+\omega+1,q}-h_{p+\omega+1,q}k_{p,q},\quad h_{p,q}k_{p,q+\omega+1}-h_{p,q+\omega+1}k_{p,q}.</math></center>
*Si ω est égal à 0, le polynôme suivant est un monôme non nul :<center><math>h_{p,q}k_{p+1,q+1}-h_{p+1,q+1}k_{p,q}.</math></center>
}}

La figure de droite illustre cette situation. Toutes les cases du carré bleu contiennent la même réduite. Si l'une de ses cases est un élément de la suite des approximants de Padé et si la suivante est l'une des deux cases rouges, le numérateur partiel de la fraction continue est un monôme non nul et le dénominateur un polynôme de terme constant non nul. Dans le cas où le carré se réduit à une unique case, il est aussi possible de choisir la case diagonale en bas à droite.

{{Démonstration/début|titre=Démonstrations}}
* '''Pour une fraction continue simple :'''
*# '''''h''{{ind|0}}(0) et les ''k{{ind|n}}''(0) sont non nuls :''' immédiat pour ''h''{{ind|0}}, ''k''{{ind|0}} et ''k''{{ind|1}} (pour les deux types) et s'en déduit pour les autres ''k{{ind|n}}'', d'après [[Fraction continue généralisée#Transformations de Möbius|leur définition par récurrence]] (''k{{ind|n}}'' = β''{{ind|n}}k''{{ind|''n''–1}} + α''{{ind|n}}k''{{ind|''n'' – 2}}) et le fait qu'en 0, les α<sub>''n''</sub> sont nuls (pour ''n'' > 0) mais pas les β<sub>''n''</sub> (on peut montrer de même que tous les ''h{{ind|n}}''(0) sont non nuls) ;
*# '''les ''h{{ind|n}}k''{{ind|''n''+1}} – ''h''{{ind|''n''+1}}''k{{ind|n}}'' sont des monômes non nuls :''' d'après [[Fraction continue généralisée#Réduites|leur expression comme produits de monômes –α<sub>''i''</sub>]] ;
*# '''la suite de leurs degrés est strictement croissante''' : puisque le quotient par ce monôme du monôme suivant est –α{{ind|''n''+2}} ;
*# '''''k''{{ind|0}} = 1 (type {{math|(''I'')}}) ou ''h''{{ind|0}} est constant (type {{math|(''II'')}}) :''' immédiat ;
*# '''''h{{ind|n}}'' et ''k{{ind|n}}'' sont premiers entre eux''' : conséquence directe des propriétés 1 et 2 ;
* '''Si des polynômes ''h{{ind|n}}'', ''k{{ind|n}}'', vérifient 1, 2, 3 et 4, alors la fraction continue associée est simple :'''
** Les α{{ind|''n''}} pour ''n'' > 0 sont bien des monômes non constants, comme quotients d'un monôme par un autre, de degré strictement plus petit. Dans le cas du type {{math|(''II'')}}, α{{ind|0}} est la constante non nulle ''h''{{ind|0}}.
** Par construction, β{{ind|0}} et β{{ind|1}} sont bien des polynômes non nuls en 0, et chaque β{{ind|''n''+2}} est le quotient d'un polynôme par un monôme. Or β{{ind|''n''+2}} est aussi égal à (''k''{{ind|''n''+2}} – α{{ind|''n''+2}}''k{{ind|n}}'')/''k''{{ind|''n''+1}}, où les ''k{{ind|i}}''(0) sont non nuls et α{{ind|''n''+2}}(0) est nul. C'est donc un polynôme, non nul en 0.
* '''Si ''f''<sub>[''p<sub>n</sub>'', ''q<sub>n</sub>'']</sub> est une suite d'approximants de Padé telle que (''p{{ind|n}}'') et (''q{{ind|n}}'') sont croissantes et si les'''<center><math>h_{p_n,q_n}k_{p_{n+1},q_{n+1}}-h_{p_{n+1},q_{n+1}}k_{p_n,q_n}</math></center>'''sont des monômes non nuls, alors la suite de leurs degrés est strictement croissante :'''<br />Le (''n'' + 1){{e}} monôme peut se réécrire<center><math>k_{p_{n+1},q_{n+1}}\Big(k_{p_{n+2},q_{n+2}}f-h_{p_{n+2},q_{n+2}}\Big)-k_{p_{n+2},q_{n+2}}\Big(k_{p_{n+1},q_{n+1}}f-h_{p_{n+1},q_{n+1}}\Big).</math></center>Dans cette différence, la série entière de droite n'a que des monômes de degrés au moins égaux à ''p''<sub>''n''+1</sub>&nbsp;+&nbsp;''q''<sub>''n''+1</sub>&nbsp;+&nbsp;1, et celle de gauche aussi d'après l'hypothèse de croissance. Le degré du (''n'' + 1){{e}} monôme est par conséquent au moins égal à ''p''<sub>''n''+1</sub>&nbsp;+&nbsp;''q''<sub>''n''+1</sub>&nbsp;+&nbsp;1. Il est donc bien strictement supérieur à celui du ''n''{{e}}, puisque ce dernier ne peut dépasser le maximum des deux valeurs ''p<sub>n</sub>&nbsp;+&nbsp;q''<sub>''n''+1</sub> et ''p''<sub>''n''+1</sub>&nbsp;+&nbsp;''q<sub>n</sub>''.
* '''Les deux polynômes suivants sont des monômes non nuls :'''<center><math>h_{p,q}k_{p+\omega+1,q}-h_{p+\omega+1,q}k_{p,q},\quad h_{p,q}k_{p,q+\omega+1}-h_{p,q+\omega+1}k_{p,q}.</math></center>Le raisonnement est analogue au précédent. Démontrons la proposition pour le premier polynôme. Il s'exprime aussi comme la différence des deux séries entières suivantes :<center><math>k_{p,q}\Big(k_{p+\omega+1,q}f-h_{p+\omega+1,q}\Big)-k_{p+\omega+1,q}\Big(k_{p,q}f-h_{p,q}\Big),</math></center>ce qui montre qu'il ne contient que des monômes de degrés au moins ''p'' + ''q'' + ω + 1. Or son degré est exactement ''p'' + ''q'' + ω + 1. C'est donc un monôme non nul.
{{Démonstration/fin}}

=== Fraction continue régulière ===
Même en imposant à la fraction continue d'être simple, il reste encore un très vaste choix de suites d'approximants disponibles. Il suffit de partir du bord de la table, c'est-à-dire d'un indice (''p'', 0) ou (0, ''q''), et de sauter d'un carré au carré adjacent situé soit en bas soit à droite. Deux choix au moins se présentent à chaque étape et trois si la largeur du carré est égal à 1. Une nouvelle contrainte peut encore simplifier massivement les calculs :

* Une fraction simple est dite '''régulière''' lorsque<ref>{{Harvsp|Padé|1892|p=54}}, aux notations près.</ref> tous les α<sub>''n''</sub> sont de même degré sauf éventuellement α{{ind|0}} et α{{ind|1}} et tous les β<sub>''n''</sub> sont de même degré sauf éventuellement β{{ind|0}}.

Une telle recherche est plus facile à formuler si les carrés de même réduites de la table de Padé présentent une géométrie régulière. Étudions le cas où ω, la longueur du côté d'un carré quelconque, est toujours égale à 1. Trois déplacements sont possibles : vers le bas, la droite ou en diagonale, en bas à droite, pour assurer la simplicité. Ajouter la régularité impose un déplacement précis. Il en existe de trois types différents.
[[Image:Padé type 1.jpg|vignette|<center>Fraction continue régulière de type I,</center>privée de la première case de sa trajectoire<!--donc image à rectifier-->]]
[[Image:Padé type 2.jpg|vignette|<center>Fraction continue régulière de type II</center>]]

Pour le premier type, on démarre soit par une descente le long du bord vertical, soit par un déplacement à droite le long du bord horizontal<ref>{{Harvsp|Padé|1899|p=411}}.</ref>, puis la suite suit une trajectoire en escalier, illustrée sur la figure de gauche et composée d'une alternance régulière de déplacements une fois vers la droite puis une fois vers le bas. Les numérateurs partiels de la fraction continue sont des monômes du premier degré et les dénominateurs partiels sont égaux à 1 (sauf quelques termes initiaux, {{cf.}} définition ci-dessus). Par exemple, pour la fonction exponentielle, Lagrange démarre par la réduite d'indice (0, 0) et obtient, en suivant la trajectoire illustrée en rouge, le résultat suivant :
<center><math>\exp(t) = 1 + \frac{t\mid}{\mid 1} - \frac{\frac 12 t \mid}{\mid 1} + \frac{\frac 16t \mid}{\mid 1} - \frac{\frac 16t \mid}{\mid 1} + \frac{\frac 1{10}t \mid}{\mid 1} - \frac{\frac 1{10}t \mid}{\mid 1} + \frac{\frac 1{14}t \mid}{\mid 1} - \frac{\frac 1{14}t \mid}{\mid 1}+\cdots </math></center>
L'approche consistant à construire la fraction continue à l'aide d'une relation de récurrence est déjà présente dans celle des « [[Algorithme d'Euclide|divisions successives]] » que faisait Lambert<ref>Méthode que {{Harvsp|Padé|1892|p=39}}, pour sa part, qualifie de {{Citation|fort pénible}}.</ref> et [[#Fonction tangente|son développement de la fonction tangente]], divisée par ''t'' pour être non nulle en 0, fait aussi partie de ce type : le développement de la fonction paire tan(''t'')/''t'' (vue comme une fonction de ''t''{{2}})
<center><math>\frac{\tan t}t=\frac{1\mid}{\mid1}-\frac{t^2\mid}{\mid3}-\frac{t^2\mid}{\mid5}-\frac{t^2\mid}{\mid7}+\cdots</math></center>
correspond à la trajectoire (0, 0), (0, 1), (1, 1), (1, 2)…

Le type II correspond à la figure de droite. Il suppose un déplacement constant soit vers le bas, soit vers la droite. Le numérateur est toujours un monôme de degré 1 mais cette fois, le dénominateur est un polynôme de degré 1. Lorsque ce déplacement se fait le long de la première ligne du tableau (en vert sur la figure), la fraction continue obtenue est celle déduite d'[[Formule de fraction continue d'Euler|une formule d'Euler]] et ses réduites sont les sommes partielles de la série entière. [[Formule de fraction continue d'Euler#Fonction exponentielle|Un premier exemple]] est celle associée à la série exponentielle :

<center><math>\exp(t) = \frac{1\mid}{\mid 1} - \frac{t \mid}{\mid 1+t} - \frac{\frac 12t \mid}{\mid 1+ \frac 12t} - \frac{\frac 13t \mid}{\mid 1+\frac 13 t} - \frac{\frac 14t \mid}{\mid 1+\frac 14t} - \cdots </math></center>
[[Formule de fraction continue d'Euler#Arctangente hyperbolique|Un autre exemple]] est donné par la fonction paire suivante (vue comme fonction de ''t''{{2}}) issue du [[logarithme]] :

<center><math>\frac1t{\rm artanh}(t)=\frac1{2t}{\rm Log}\left(\frac{1+t}{1-t}\right)=\frac{1\mid}{\mid 1} - \frac{\frac 13 t^2 \mid}{\mid 1+\frac 13 t^2} - \frac{\frac 35t^2 \mid}{\mid 1+ \frac 35t^2} - \frac{\frac 57 t^2 \mid}{\mid 1+\frac 57 t^2} - \frac{\frac 79t^2 \mid}{\mid 1+\frac 79t^2} - \cdots
</math></center>
On a longtemps cru que ce type de fraction continue n'était qu'une exception<ref>{{Harvsp|Padé|1892|p=41}}, fait référence à {{de}} [[Eduard Heine]], ''Handbuch der Kugelfunctionnen'', vol. 1, 1878, [https://books.google.fr/books?id=D79hEMl2GM0C&pg=PA266 {{p.|266}} et s.].</ref>. Il existe pourtant une infinité de fractions continues de cette nature — une par ligne et une par colonne — pour chaque fonction analytique, sous réserve que chaque carré soit de côté 1.

[[Image:Padé type 3.jpg|vignette|<center>Fraction continue régulière de type III</center>]]

Le troisième type est celui de l'exemple du [[#Fonction exponentielle|§ « Fonction exponentielle »]] ci-dessus. Le numérateur est un monôme de degré 2 et le dénominateur un polynôme de degré 1, avec une exception possible : il arrive que le monôme dominant du dénominateur se simplifie systématiquement. Un exemple se trouve dans l'article [[Approximant de Padé de la fonction exponentielle]]. Cet article traite exclusivement d'une fonction particulière mais la logique décrite est générale aux tables de Padé dont les carrés sont tous de côtés 1 et le raisonnement se généralise à des tables plus complexes.

== Convergence ==
=== Structure du domaine ===
{{refnec|date=avril 2014|L'un des attraits de l'approximant de Padé est sa capacité à avaler les pôles d'une fonction analytique.}} La [[Problème de convergence|question de la convergence]] est, en conséquence, particulièrement cruciale. Elle est néanmoins difficile, contrairement à ce que laisseraient espérer l'exemple de la fonction exponentielle et le développement particulier de la fonction tangente étudié ici. Sans même aborder la délicate question du lien, en un point où toutes deux sont définies, entre la limite de la fraction continue et la série entière ''f'' de départ, ce paragraphe se borne, comme {{Harv|Padé|1892}}, à quelques remarques élémentaires sur la convergence des fractions continues simples.

Soit (''h{{ind|n}}''/''k{{ind|n}}'') la suite des réduites d'une fraction continue simple de type {{math|(''I'')}}. On a donc, pour tout entier ''n'' > 0, ''h{{ind|n}}''(''t'')''k''{{ind|''n''–1}}(''t'') – ''h''{{ind|''n''–1}}(''t'')''k{{ind|n}}''(''t'') = ''a''{{ind|ν{{ind|''n''}}}}''t''{{exp|ν{{ind|''n''}}}}, où la suite des exposants ν{{ind|''n''}} est strictement croissante, si bien qu'en notant ''d{{ind|n}}'' le polynôme ''k{{ind|n}}k''{{ind|''n''–1}} :
<center><math>\frac{h_n(t)}{k_n(t)}-\frac{h_{n-1}(t)}{k_{n-1}(t)}=\frac{a_{\nu_n}}{d_n(t)}t^{\nu_n}</math></center>et
<center><math>\frac{h_n(t)}{k_n(t)}=h_0(t)+\sum_{i=1}^n\frac {a_{\nu_n}}{d_n(t)}t^{\nu_n}.</math></center>

Notons ''D'' le disque de convergence de la série entière de terme général ''a<sub>k</sub>t<sup>k</sup>'' (en posant ''a{{ind|k}}'' = 0 si ''k'' n'est pas égal à l'un des ν<sub>''n''</sub>).
* Pour ''N'' et ''c'' > 0 fixés, soit ''H'' l'ensemble des ''t'' pour lesquels, à partir de l'indice ''N'', le [[module d'un nombre complexe|module]] de ''d{{ind|n}}''(''t'') reste supérieur à ''c''. Sur l'intersection de ''H'' et du disque ouvert ''D'', la fraction continue converge (elle converge même [[Convergence uniforme|uniformément]] sur tout [[Compacité (mathématiques)|compact]] de ''H''⋂''D'').
* En tout point ''t'' extérieur au disque fermé {{surligner|''D''}} et pour lequel la suite (''d{{ind|n}}''(''t'')) reste [[Suite bornée|bornée]], la fraction continue diverge.

=== Illustration par l'exemple de la fonction tangente ===
Modifions un peu les conventions du § « Fonction tangente » [[#Fonction tangente|ci-dessus]] pour rendre égal à 1 le terme constant de chaque dénominateur ''k<sub>n</sub>'' :
<center><math>\forall n \ge 2\quad h_n(t) = h_{n-1}(t) -\frac {t^2}{(2n-1)(2n-3)}h_{n-2}(t) \quad \text{et}\quad k_n(t) = k_{n-1}(t) -\frac {t^2} {(2n-1)(2n-3)} k_{n-2}(t).</math></center>

Les fractions sont les mêmes à un facteur multiplicatif près du numérateur et du dénominateur, qui ne modifie donc pas leurs valeurs. Elles sont maintenant normalisées. Avec les notations du paragraphe précédent, on a :
<center><math>\forall n\in\N\quad a_{2n}=\prod_{i=1}^n\frac1{2i-1}\quad\text{et}\quad a_{2n+1}=0.</math></center>

Il suffit de minorer 2''i'' – 1 par ''i'' pour conclure que la série entière de coefficients ''a{{ind|n}}'' est majorée terme à terme par exp(''t''<sup>2</sup>). Son [[rayon de convergence]] est donc infini.

[[Image:Padé Tangente dénominateur.jpg|vignette|upright=2|centré|<center>La suite des dénominateurs de la fraction continue converge manifestement vers la fonction cosinus.</center>]]
Il suffit alors de vérifier que la suite des dénominateurs n'a pas tendance à se rapprocher de zéro en dehors des pôles de la fonction tangente. Le graphique ci-dessus laisse prévoir que cette suite converge vers la fonction cosinus. Un calcul ({{cf.}} boîte déroulante ci-dessous) confirme que cette convergence a lieu sur tout le [[plan complexe]] et est uniforme sur tout compact. La suite (''d''<sub>n</sub>(''t'')) converge donc vers cos{{2}}(''t'').

Soit ''K'' un [[Fermé (topologie)|fermé]] du plan complexe, [[Partie bornée|borné]] par un réel ''M'', et ne contenant aucun des zéros de la fonction cosinus (qui sont les (2''j''+1)π/2 pour tout entier ''j''). D'après le paragraphe précédent, la suite des approximants de Padé est uniformément convergente sur ''K''. (En effet, il existe une constante strictement positive 2''c'' telle que le module de cos{{2}}(''t'') soit supérieur à 2''c'' pour tout ''t'' de ''K'' donc il existe un rang ''N'' tel que toute fonction de la suite (''d<sub>n</sub>''(''t'')) d'indice supérieur à ''N'' est, en module, supérieure à ''c'' sur ''K''.)

Cela ne suffit pas à prouver que la limite est bien la fonction tangente. Lambert, lui, montre directement<ref>{{Harvsp|Lambert|1761|loc= § 24-29}}.</ref> non seulement cette convergence des dénominateurs vers cosinus mais aussi des numérateurs vers sinus.

{{Démonstration|titre=Preuve de la convergence des dénominateurs|
La technique consiste à déterminer la suite des coefficients des polynômes du dénominateur. (Une autre méthode utilise le calcul différentiel<ref name=GHL>On en trouve une version sur : M.Gouy, G. Huvent et A. Ladureau, ''[http://gery.huvent.pagesperso-orange.fr/irem/pade.pdf Approximants de Padé]'', [[IREM]] de Lille, 2002.</ref>. Elle n'est pas choisie ici car elle n'illustre pas le caractère structurel du domaine de convergence à l'aide de la logique de Padé et ne montre que la convergence de la fraction continue vers la fonction tangente et non celle spécifique de son dénominateur.)

* '''Calcul de coefficients des ''k<sub>n</sub>'' :'''
Lambert calcule les polynômes ''k{{ind|n}}''(''t'') pour les premières valeurs de ''n'' puis devine la formule générale :
<center><math>k_n(t)=\sum (-1)^ib_{i,n}t^{2i}\quad\text{avec}
\quad b_{i,n}=\begin{cases}\frac1{(2i)!}\frac{(2n-2i)(2n-2i-2)\ldots(2n-4i+2)}{(2n-1)(2n-3)\ldots(2n-2i+1)}&\text{si }0\le 2i\le n,\\0&\text{sinon,}\end{cases}</math></center>
qu'il démonte alors facilement, grâce à la relation de récurrence
<center><math>\forall n \ge 2\quad b_{i,n}=b_{i,n-1}+\frac{b_{i-1,n-2}}{(2n-1)(2n-3)}.</math></center>

* '''Calcul de leurs limites :'''
Dans l'expression de ''b{{ind|i,n}}'', le produit des ''i'' facteurs du numérateur est inférieur, terme à terme, à celui des ''i'' facteurs correspondants du dénominateur, si bien que ''b{{ind|i,n}}'' est majoré par 1/(2i)!.

Lambert affirme que ce majorant est en fait la limite de ''b{{ind|i,n}}'' quand ''n'' tend vers l'infini. Il ne prend pas la peine de le démontrer mais c'est un exercice facile, par exemple en utilisant la [[formule de Stirling]]. En chaque degré, le coefficient de ''k{{ind|n}}''(''t'') a donc pour limite celui de la série entière cos(''t'').

* '''Convergence de la suite des polynômes ''k<sub>n</sub>'' :
La convergence uniforme sur un disque ''E'' de rayon ''R'' est maintenant simple à établir. Le raisonnement est analogue à [[Approximant de Padé de la fonction exponentielle#Convergence des réduites|celui sur la fonction exponentielle]]. Soit ε un réel strictement positif. Fixons un entier naturel ''m'' tel que<center><math>\sum_{i=m}^{\infty}\frac{R^{2i}}{(2i)!}\le\frac{\varepsilon}2.</math></center>Alors, pour tout entier naturel ''n'',
<center><math>\forall t\in E\quad|\cos(t)-k_n(t)|\le\sum_{i=0}^{\infty}\left(\frac1{(2i)!}-b_{i,n}\right)R^{2i}\le R_n+\frac{\varepsilon}2</math></center>
avec
<center><math>R_n=\sum_{0\le i<m}\left(\frac1{(2i)!}-b_{i,n}\right)R^{2i}.</math></center>

La proposition précédente montre qu'il existe ''N'' tel que pour tout ''n ≥ N'', ''R{{ind|n}}'' ≤ ε/2. On en déduit :
<center><math>\forall n\ge N\quad\forall t\in E\quad|\cos(t)-k_n(t)|\le\varepsilon.</math></center>
}}

== Notes et références ==
{{références|colonnes=2}}

== Voir aussi ==
=== Article connexe ===
[[Fonction de Bessel]]
=== Bibliographie ===
* Jean-Étienne Rombaldi, ''Interpolation & approximation, analyse pour l'agrégation'', Vuibert, 2005
*{{Ouvrage|langue=en|auteur1=C. Brezinski|titre=History of Continued Fractions and Padé Approximants|lieu=Berlin|éditeur=Springer-Verlag|année=1991|doi=10.1007/978-3-642-58169-4|lire en ligne=https://books.google.com/books?id=rxzsCAAAQBAJ}}
* {{en}} C. Brezinski et M. Redivo-Zaglia, ''Extrapolation Methods: Theory and Practice'', North-Holland, 1991 {{ISBN|978-0-44488814-3}}
*{{en}} G. A. Baker et P. 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}}
*H. Padé, ''Œuvres'' rassemblées et présentées par Claude Brezinski, Librairie Albert Blanchard, Paris, 1984.

=== Liens externes ===
*D. Vekemans, ''[http://vekemans.free.fr/E1/pade.pdf Approximants de Padé]'', Laboratoire de [[mathématiques pures]] et appliquées Joseph Liouville
*{{MathWorld|nom_url=PadeApproximant|titre=Padé Approximant}}

{{Portail|analyse}}

[[Catégorie:Analyse complexe]]
[[Catégorie:Fraction continue]]
[[Catégorie:Transformation de suite]]