Approximant de Padé

Le concept de l'article doit son nom au mathématicien français Henri Padé (1863-1953).

En mathématiques, et plus précisément en analyse complexe, l'approximant de Padé est une méthode d'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 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 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 prolonger des fonctions analytiques et d'étudier certains aspects de la question des 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'irrationalité du [[e (nombre)|nombre Modèle:Math]]. Une stratégie plus élaborée permet à Jean-Henri Lambert de montrer celle de [[pi|Modèle: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

Charles Hermite utilise ce qui est maintenant appelé un approximant de Padé pour démontrer la transcendance de Modèle:Math en 1873.

Il est utile d'approcher une fonction donnée par une suite de fonctions aisément calculables. Cette démarche est à l'origine de la théorie des séries entières consistant à approcher de plus en plus précisément une fonction analytique à l'aide de la suite de ses 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 fraction rationnelle. « Le mieux possible » signifie ici que pour un 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 fractions continues (généralisées). Dans le cas des polynômes, les développements limités sont indexés par un 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|Modèle:Math]] ou [[E (nombre)|Modèle:Math]] est plus simplement mise en valeur. Les preuves de l'irrationalité de Modèle:Math et de Modèle:Math, comme celles de leur transcendance, se fondent naturellement sur cet outil. 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ô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

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 tModèle:Ind est un point de domaine de définition de f et t un scalaire de module strictement inférieur à r, le rayon de convergence de f au point tModèle:Ind, il est possible de développer f de la manière suivante : f(tModèle:Ind + t) = αModèle:Ind + αModèle:IndtModèle:Exp + …. Ici, αModèle:Ind ainsi que les coefficients suivants, sont choisis non nuls. Pour cette raison, nModèle:Ind n'est pas nécessairement égal à 1. Il existe une fonction analytique fModèle:Ind définie en 0 et telle que

$f(t_0+t)=a_0+\frac{t^{n_1}}{f_1(t)}\text{avec}{a}_0={\alpha }_0.$

Un traitement analogue sur la fonction fModèle:Ind, puis fModèle:Ind, montre qu'il est encore possible d'écrire la fonction f de la manière suivante :

$f(t_0+t)=a_0+\frac{t^{n_1}}{a_1+\frac{t^{n_2}}{f_2(t)}}=a_0+\frac{t^{n_1}}{a_1+\frac{t^{n_2}}{a_2+\frac{t^{n_3}}{f_3(t)}}}=\cdots$

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 fModèle:Ind 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 nombres rationnels pour la fraction continue.

On utilise dans cet article la notation de Pringsheim :

$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$

Si le processus est arrêté à la pModèle:E étape, on obtient une fraction rationnelle qui « approche » la fonction initiale f et qui est un exemple d'approximant de Padé.

Pour une raison de simplicité, dans la suite de l'article, la valeur tModèle:Ind est choisie nulle.

Fonction tangente

Pour illustrer cette démarche, considérons la fonction fonction tangente, l'un des premiers exemples historiques de fonction analytique dont le développement en fraction continue est calculé, par Lambert :

$\tan (t)=\frac{t}{1-\frac{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$

Appliquons l'algorithme naïf du paragraphe précédent au point t₀ = 0. La valeur a₀ est égale à 0 et n₁ à 1. On obtient une première fraction rationnelle h₀(t)/k₀(t) approchant la fonction tangente :

$\tan (t)=0+\frac{t}{f_1(t)},f_1(t)=\frac{t}{\tan (t)}\text{et}{h}_0(t)=0,k_0(t)=1.$

Prolongeons le calcul :

$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)}.$

On obtient une nouvelle expression plus précise de la fonction tangente :

$\tan (t)=0+\frac{t\mid }{\mid 1}-\frac{t^2\mid }{\mid f_2(t)},f_2(t)=\frac{t^2\tan t}{\tan t-t}\text{et}{h}_1(t)=t,k_1(t)=1.$

Par une autre méthode, Lambert découvre une formule de récurrence permettant de calculer plus rapidement h_n(t)/k_n(t) :

$\mathrm{\forall }n\ge 2h_n(t)=(2n-1)h_{n-1}(t)-t^2{h}_{n-2}(t)\text{et}{k}_n(t)=(2n-1)k_{n-1}-t^2{k}_{n-2}(t).$

Modèle:Démonstration/début Lambert transpose de la manière suivante (aux notations près) l'algorithme d'Euclide usuellement appliqué au calcul des fractions continues de réels : il « divise » d'abord RModèle:Ind := cos(t) par Modèle:Nobr ce qui lui donne un « quotient » Modèle:Nobr et un « reste » –RModèle:Ind :

$\cos t=\frac{\sin t}{t}-R_1\text{ donc }\tan t=\frac{t}{1-\frac{t{R}_1}{R_0}}.$

Il prévoit qu'indéfiniment, RModèle:Ind divisé par RModèle:Ind donne un quotient QModèle:Ind = (2n – 1)/t (et un reste –RModèle:Ind) :

$R_{n-2}=\frac{2n-1}{t}{R}_{n-1}-R_n\text{ donc }\frac{t{R}_{n-1}}{R_{n-2}}=\frac{t^2}{2n-1-\frac{t{R}_n}{R_{n-1}}}.$

Cela fournit le développement annoncé, sous réserve de prouver que la valuation des RModèle:Ind, définis par

$R_{-1}=\cos ,R_0=\sin \mathrm{e}\mathrm{t}\mathrm{\forall }n>0R_n=\frac{2n-1}{t}{R}_{n-1}-R_{n-2},$

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 :

$\mathrm{\forall }n\in \mathbb{N}{R}_n=t^{n+1}{2}^n\sum _{k=0}^{\mathrm{\infty }}\frac{(-1{)}^k{t}^{2k}(n+k)!}{k!(2n+2k+1)!}$

puis de vérifier cette égalité par une récurrence d'ordre 2 : si elle est vraie aux ordres n – 2 et n – 1, on trouve bien, en remplaçant RModèle:Ind et RModèle:Ind par leurs expressions et en regroupant les termes :

$R_n=\frac{2n-1}{t}{R}_{n-1}-R_{n-2}=t^{n-1}{2}^n\sum _{k=0}^{\mathrm{\infty }}\frac{k(-1{)}^{k-1}{t}^{2k}(n+k-1)!}{k!(2n+2k-1)!}=t^{n+1}{2}^n\sum _{k=1}^{\mathrm{\infty }}\frac{(-1{)}^{k-1}{t}^{2(k-1)}(n+k-1)!}{(k-1)!(2n+2(k-1)+1)!}.$

La relation de récurrence sur les numérateurs hModèle:Ind et les dénominateurs kModèle:Ind des fractions partielles est alors une conséquence immédiate de l'algorithme d'Euclide étendu. Modèle:Démonstration/fin

Il démontre de plus la [[#Illustration par l'exemple de la fonction tangente|convergence des quotients hModèle:Ind/kModèle:Ind vers la fonction tangente]]. Il remarque même que Modèle:Citation, ce qui est la clé de sa démonstration de l'irrationalité de Modèle:Math.

Modèle:Démonstration/début Par parité des kModèle:Ind et imparité de tan et des hModèle:Ind, il suffit de prouver la convergence sur l'intervalle [0, 1].

La méthode utilisée ici consiste à étudier la suite des fonctions ψ_n(t) définies sur [0, 1] par

${\psi }_n(t)=k_n(t)\sin t-h_n(t)\cos t,$

autrement dit, par

${\psi }_0(t)=\sin t,{\psi }_1(t)=\sin t-t\cos t\text{et}\mathrm{\forall }n\ge 2{\psi }_n(t)=(2n-1){\psi }_{n-1}(t)-t^2{\psi }_{n-2}(t).$

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 récurrence d'ordre 2) que pour tout entier n > 0, la dérivée de ψ_n(t) est égale à tψModèle:Ind(t) puis, que chaque fonction ψModèle:Ind est nulle en 0 et strictement croissante (donc positive). On en déduit alors que la suite des ψ_n tend uniformément vers zéro et même (par récurrence, en utilisant que ψ_n est la primitive de tψModèle:Ind(t) qui s'annule en 0) que

$\mathrm{\forall }t\in [0,1]{\psi }_n(t)\le \frac{t^{2n+1}}{2^{n}n!}.$

Pour en déduire que ψ_n/(kModèle:Indcos) converge aussi vers 0, on a besoin de minorer kModèle:Ind. Une récurrence montre que pour tout entier n > 0,

$\mathrm{\forall }t\in [0,1]2^{n-1}(n-1)!\le k_n(t)\le \frac{(2n)!}{2^{n}n!}.$

On dispose donc de l'encadrement suivant :

$\mathrm{\forall }t\in [0,1]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)},$

ce qui montre non seulement que la suite (hModèle:Ind/kModèle:Ind) converge uniformément sur [0, 1] vers la fonction tangente mais aussi que, pour n fixé, hModèle:Ind/kModèle:Ind et tangente ont même développement limité en 0 à l'ordre 2n. Modèle:Démonstration/fin

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 Modèle:Math, où k est un entier, hModèle:Ind(t)/kModèle:Ind(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 à Modèle:Sqrt au lieu de Modèle:Math. 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 Modèle:Math. Elle possède deux pôles sur les réels positifs. La huitième approximation, en vert, « colle à la fonction tangente » sur l'intervalle Modèle:Math 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 converge uniformément sur la réunion des intervalles Modèle:Math et Modèle:Math où j varie de 1 à n.

Modèle:Refnec 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 Modèle:Math. Modèle:Référence insuffisante l'étude de la fonction zêta de Riemann.

Fonction exponentielle

Modèle:Article détaillé 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 :

$\exp (x)=1+x+\frac{\frac{1}{2}{x}^2\mid }{\mid 1-\frac{1}{3}x}+\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$

L'algorithme utilisé est ici un peu différent. Les numérateurs ne sont plus des constantes mais des fonctions affines. En revanche une propriété reste commune : le développement limité des fractions h_p/k_p à 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é pour introduire sa théorie.

Fragments d'histoire

Origine

Lagrange énonce une propriété sur des réduites de fractions continues correspondant à notre définition moderne d'approximant de Padé.

Padé place l'origine de ses approximants aux travaux d'Euler (1707-1783) avec sa démonstration de l'irrationalité de Modèle:Math. Il est vrai qu'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 Modèle:Math]]). Lambert (1707-1777) franchit un pas fondamental en développant des fonctions, ce qui lui permet de démontrer que si un nombre réel non nul t est rationnel, alors ni tan(t), ni exp(t) ne l'est, démontrant par là l'irrationalité de Modèle:Math. Il développe en fraction continue non seulement tan(t) mais tanh(t), arctan(t), log(1 + t)Modèle: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 (1736-1813) généralise l'emploi des fractions continues à l'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'é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 Modèle:Citation, ce qui correspond à une définition moderne des approximants de Padé. Carl Friedrich Gauss (1777-1855) utilise les résultats de Lagrange et généralise l'usage des fractions continues au plan complexe ; les fractions continues de Gauss permettent de retrouver des développements en fraction continue de nombreux quotients de 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

Riemann met en évidence un des intérêts des approximants de Padé : ils permettent parfois de prolonger des 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 (1826-1866). 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 Modèle:Citation. Riemann étudie la convergence des fractions continues de Gauss. Modèle:Refnec démonstration du théorème d'Apéry, prouvé en 1977 et énonçant que l'image de 3 par la fonction zêta de Riemann est un nombre irrationnel.

Charles Hermite, le directeur de thèse de Padé, é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 Modèle:Math, le premier nombre démontré transcendant qui n'a pas été spécialement conçu pour une démonstration (comme ceux de Liouville). Sur la base de ces idées, Ferdinand von Lindemann démontre la transcendance de Modèle:Math en 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 : Modèle:Citation é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é.

Il est donc naturel que la question des fractions continues, construites à l'aide de fonctions rationnelles, soient à l'ordre du jour à la fin du Modèle:XIXe siècle. Cependant, en 1890, Poincaré décrit la théorie des fractions continues algébriques comme une Modèle:Citation dont Modèle:Citation. 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, 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, 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 convergence des différentes fractions continues que l'on peut construire. En 1894, Thomas-Joannes Stieltjes est d'un apport considérable sur cette question. Il établit exactement, pour une vaste famille de fractions continues, le domaine de convergence ; 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. Son jury est composé par des grands noms de l'époque comme Poincaré ou Picard. Il est remporté par Padé qui dit présenter Modèle:Citation, tout en Modèle:Citation

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ⁿ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 :

• 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 + t². 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 + tModèle: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_[p,q].

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 premiers entre eux, puis rendre égal à 1 le terme constant de k (qui est non nul puisque 0 n'est pas un pôle) : Modèle:Théorème

• 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)) :
  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 noyau est donc non réduit au vecteur nul. Soit t^ωk(t) un polynôme non nul de cet espace, avec ω le plus grand possible et k(t) de terme constant égal à 1.
  Soit h₁(t) le développement à l'ordre p de t^ωk(t)f(t). Le développement à l'ordre p + q de t^ωk(t)f(t) – h₁(t) est nul par construction, ce qui montre que h₁(t) est un multiple de t^ω. Soit h(t) le polynôme h₁(t) divisé par t^ω. 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^ω+ω'(k(t)/d(t)) appartient au noyau ci-dessus donc ω' = 0.
• Pour un tel (h, k), toute réduite est égale à h/k :
  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.

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é

Modèle:Voir 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 trivialement pas le cas pour une fonction f 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(tModèle:2). Voici par exemple la table de la fonction arctan(t)/t :

Table de Padé pour $\arctan (t)/t$

$\phantom{p}\begin{array}{c}\\ q\end{array}p$ $0$ $1$ $2$ $3$ $4$ $5$ $0$ $1$ $1$ $1-\frac{1}{3}{t}^2$ $1-\frac{1}{3}{t}^2$ $1-\frac{1}{3}{t}^2+\frac{1}{5}{t}^4$ $1-\frac{1}{3}{t}^2+\frac{1}{5}{t}^4$ $1$ $1$ $1$ $1-\frac{1}{3}{t}^2$ $1-\frac{1}{3}{t}^2$ $1-\frac{1}{3}{t}^2+\frac{1}{5}{t}^4$ $1-\frac{1}{3}{t}^2+\frac{1}{5}{t}^4$ $2$ $\frac{1}{1+\frac{1}{3}{t}^2}$ $\frac{1}{1+\frac{1}{3}{t}^2}$ $\frac{1+\frac{4}{15}{t}^2}{1+\frac{9}{15}{t}^2}$ $\frac{1+\frac{4}{15}{t}^2}{1+\frac{9}{15}{t}^2}$ $\frac{1+\frac{8}{21}{t}^2+\frac{4}{105}{t}^4}{1+\frac{5}{7}{t}^2}$ $\frac{1+\frac{8}{21}{t}^2+\frac{4}{105}{t}^4}{1+\frac{5}{7}{t}^2}$ $3$ $\frac{1}{1+\frac{1}{3}{t}^2}$ $\frac{1}{1+\frac{1}{3}{t}^2}$ $\frac{1+\frac{4}{15}{t}^2}{1+\frac{9}{15}{t}^2}$ $\frac{1+\frac{4}{15}{t}^2}{1+\frac{9}{15}{t}^2}$ $\frac{1+\frac{8}{21}{t}^2+\frac{4}{105}{t}^4}{1+\frac{5}{7}{t}^2}$ $\frac{1+\frac{8}{21}{t}^2+\frac{4}{105}{t}^4}{1+\frac{5}{7}{t}^2}$ $4$ $\frac{1}{1+\frac{1}{3}{t}^2-\frac{4}{45}{t}^4}$ $\frac{1}{1+\frac{1}{3}{t}^2-\frac{4}{45}{t}^4}$ $\frac{1+\frac{11}{21}{t}^2}{1+\frac{6}{7}{t}^2+\frac{3}{35}{t}^4}$ $\frac{1+\frac{11}{21}{t}^2}{1+\frac{6}{7}{t}^2+\frac{3}{35}{t}^4}$ $\frac{1+\frac{7}{9}{t}^2+\frac{64}{945}{t}^4}{1+\frac{10}{9}{t}^2+\frac{5}{21}{t}^4}$ $\frac{1+\frac{7}{9}{t}^2+\frac{64}{945}{t}^4}{1+\frac{10}{9}{t}^2+\frac{5}{21}{t}^4}$ $5$ $\frac{1}{1+\frac{1}{3}{t}^2-\frac{4}{45}{t}^4}$ $\frac{1}{1+\frac{1}{3}{t}^2-\frac{4}{45}{t}^4}$ $\frac{1+\frac{11}{21}{t}^2}{1+\frac{6}{7}{t}^2+\frac{3}{35}{t}^4}$ $\frac{1+\frac{11}{21}{t}^2}{1+\frac{6}{7}{t}^2+\frac{3}{35}{t}^4}$ $\frac{1+\frac{7}{9}{t}^2+\frac{64}{945}{t}^4}{1+\frac{10}{9}{t}^2+\frac{5}{21}{t}^4}$ $\frac{1+\frac{7}{9}{t}^2+\frac{64}{945}{t}^4}{1+\frac{10}{9}{t}^2+\frac{5}{21}{t}^4}$

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 :

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 :

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. 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 :

Il suffit en effet de définir (α_n) et (β_n) par la relation de récurrence :

${\alpha }_{n+2}=-\frac{h_{n+1}{k}_{n+2}-h_{n+2}{k}_{n+1}}{h_n{k}_{n+1}-h_{n+1}{k}_n},{\beta }_{n+2}=\frac{h_n{k}_{n+2}-h_{n+2}{k}_n}{h_n{k}_{n+1}-h_{n+1}{k}_n}$

et l'initialisation (en supposant, dans le cas Modèle:Math, que kModèle:Ind = 1) :

$(I){\beta }_0=h_0,{\alpha }_1=h_1-h_0{k}_1,{\beta }_1=k_1,$ $(II){\alpha }_0=h_0,{\beta }_0=k_0,{\alpha }_1=(h_0{k}_1-h_1{k}_0)/h_0,{\beta }_1=h_1/h_0.$

Les (α_n) et (β_n) ainsi construites sont en général seulement des suites de fractions rationnelles. On peut bien sûr, dans ce cas, les convertir (sans modifier les réduites) en des suites de polynômes, mais l'inconvénient est que les nouveaux numérateurs AModèle:Ind et dénominateurs BModèle:Ind des réduites ne coïncideront alors plus avec les anciens, hModèle:Ind et kModèle:Ind : on aura seulement égalité des fractions rationnelles, AModèle:Ind/BModèle:Ind = hModèle:Ind/kModèle:Ind, 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 αModèle:Ind et βModèle:Ind, définies comme ci-dessus à partir des hModèle:Ind et kModèle:Ind, soient déjà des polynômes. Padé limite l'étude des suites (hModèle:Ind/kModèle:Ind) pour lesquelles cette condition est vérifiée, en ajoutant des contraintes :

• Une fraction continue (de type Modèle:Math ou Modèle:Math ci-dessus) est dite simple lorsque ses numérateurs partiels α_n sont des monômes non constants, sauf αModèle:Ind (pour le type Modèle:Math) qui doit être une constante non nulle et lorsque de plus, chacun de ses dénominateurs partiels β_n est un polynôme de terme constant non nul.

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 premiers entre eux, d'après la caractérisation suivante :

Cette propriété générale indique comment choisir une suite d'approximants de Padé f_{[pn, qn]} 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 pModèle:Ind ou qModèle:Ind 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 :

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 Modèle:Nobr forment les sommets du carré contenant toutes les occurrences de l'approximant de Padé f_{[p, q]}, comme indiqué à la fin du § « Table de Padé ». La propriété suivante permet de trouver les candidats pour la construction d'une fraction continue simple :

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.

• Pour une fraction continue simple :
  1. hModèle:Ind(0) et les kModèle:Ind(0) sont non nuls : immédiat pour hModèle:Ind, kModèle:Ind et kModèle:Ind (pour les deux types) et s'en déduit pour les autres kModèle:Ind, d'après leur définition par récurrence (kModèle:Ind = βModèle:IndkModèle:Ind + αModèle:IndkModèle:Ind) et le fait qu'en 0, les α_n sont nuls (pour n > 0) mais pas les β_n (on peut montrer de même que tous les hModèle:Ind(0) sont non nuls) ;
  2. les hModèle:IndkModèle:Ind – hModèle:IndkModèle:Ind sont des monômes non nuls : d'après leur expression comme produits de monômes –α_i ;
  3. la suite de leurs degrés est strictement croissante : puisque le quotient par ce monôme du monôme suivant est –αModèle:Ind ;
  4. kModèle:Ind = 1 (type Modèle:Math) ou hModèle:Ind est constant (type Modèle:Math) : immédiat ;
  5. hModèle:Ind et kModèle:Ind sont premiers entre eux : conséquence directe des propriétés 1 et 2 ;
• Si des polynômes hModèle:Ind, kModèle:Ind, vérifient 1, 2, 3 et 4, alors la fraction continue associée est simple :
  • Les αModèle:Ind 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 Modèle:Math, αModèle:Ind est la constante non nulle hModèle:Ind.
  • Par construction, βModèle:Ind et βModèle:Ind sont bien des polynômes non nuls en 0, et chaque βModèle:Ind est le quotient d'un polynôme par un monôme. Or βModèle:Ind est aussi égal à (kModèle:Ind – αModèle:IndkModèle:Ind)/kModèle:Ind, où les kModèle:Ind(0) sont non nuls et αModèle:Ind(0) est nul. C'est donc un polynôme, non nul en 0.
• Si f_{[pn, qn]} est une suite d'approximants de Padé telle que (pModèle:Ind) et (qModèle:Ind) sont croissantes et si les$h_{p_n,q_n}{k}_{p_{n+1},q_{n+1}}-h_{p_{n+1},q_{n+1}}{k}_{p_n,q_n}$sont des monômes non nuls, alors la suite de leurs degrés est strictement croissante :
  Le (n + 1)Modèle:E monôme peut se réécrire$k_{p_{n+1},q_{n+1}}(k_{p_{n+2},q_{n+2}}f-h_{p_{n+2},q_{n+2}})-k_{p_{n+2},q_{n+2}}(k_{p_{n+1},q_{n+1}}f-h_{p_{n+1},q_{n+1}}).$Dans cette différence, la série entière de droite n'a que des monômes de degrés au moins égaux à p_n+1 + q_n+1 + 1, et celle de gauche aussi d'après l'hypothèse de croissance. Le degré du (n + 1)Modèle:E monôme est par conséquent au moins égal à p_n+1 + q_n+1 + 1. Il est donc bien strictement supérieur à celui du nModèle:E, puisque ce dernier ne peut dépasser le maximum des deux valeurs p_n + q_n+1 et p_n+1 + q_n.
• Les deux polynômes suivants sont des monômes non nuls :$h_{p,q}{k}_{p+\omega +1,q}-h_{p+\omega +1,q}{k}_{p,q},h_{p,q}{k}_{p,q+\omega +1}-h_{p,q+\omega +1}{k}_{p,q}.$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 :$k_{p,q}(k_{p+\omega +1,q}f-h_{p+\omega +1,q})-k_{p+\omega +1,q}(k_{p,q}f-h_{p,q}),$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.

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 tous les α_n sont de même degré sauf éventuellement αModèle:Ind et αModèle:Ind et tous les β_n sont de même degré sauf éventuellement βModèle:Ind.

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.

Fraction continue régulière de type I,privée de la première case de sa trajectoire

Fraction continue régulière de type II

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, 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, Modèle: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 :

$\exp (t)=1+\frac{t\mid }{\mid 1}-\frac{\frac{1}{2}t\mid }{\mid 1}+\frac{\frac{1}{6}t\mid }{\mid 1}-\frac{\frac{1}{6}t\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$

L'approche consistant à construire la fraction continue à l'aide d'une relation de récurrence est déjà présente dans celle des « divisions successives » que faisait Lambert et 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 tModèle:2)

$\frac{\tan t}{t}=\frac{1\mid }{\mid 1}-\frac{t^2\mid }{\mid 3}-\frac{t^2\mid }{\mid 5}-\frac{t^2\mid }{\mid 7}+\cdots$

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'une formule d'Euler et ses réduites sont les sommes partielles de la série entière. Un premier exemple est celle associée à la série exponentielle :

$\exp (t)=\frac{1\mid }{\mid 1}-\frac{t\mid }{\mid 1+t}-\frac{\frac{1}{2}t\mid }{\mid 1+\frac{1}{2}t}-\frac{\frac{1}{3}t\mid }{\mid 1+\frac{1}{3}t}-\frac{\frac{1}{4}t\mid }{\mid 1+\frac{1}{4}t}-\cdots$

Un autre exemple est donné par la fonction paire suivante (vue comme fonction de tModèle:2) issue du logarithme :

$\frac{1}{t}\mathrm{a}\mathrm{r}\mathrm{t}\mathrm{a}\mathrm{n}\mathrm{h}(t)=\frac{1}{2t}\mathrm{L}\mathrm{o}\mathrm{g}\left(\frac{1+t}{1-t}\right)=\frac{1\mid }{\mid 1}-\frac{\frac{1}{3}{t}^2\mid }{\mid 1+\frac{1}{3}{t}^2}-\frac{\frac{3}{5}{t}^2\mid }{\mid 1+\frac{3}{5}{t}^2}-\frac{\frac{5}{7}{t}^2\mid }{\mid 1+\frac{5}{7}{t}^2}-\frac{\frac{7}{9}{t}^2\mid }{\mid 1+\frac{7}{9}{t}^2}-\cdots$

On a longtemps cru que ce type de fraction continue n'était qu'une exception. 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.

Fraction continue régulière de type III

Le troisième type est celui de l'exemple du § « 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

Modèle:Refnec La 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 Modèle:Harv, à quelques remarques élémentaires sur la convergence des fractions continues simples.

Soit (hModèle:Ind/kModèle:Ind) la suite des réduites d'une fraction continue simple de type Modèle:Math. On a donc, pour tout entier n > 0, hModèle:Ind(t)kModèle:Ind(t) – hModèle:Ind(t)kModèle:Ind(t) = aModèle:IndtModèle:Exp, où la suite des exposants νModèle:Ind est strictement croissante, si bien qu'en notant dModèle:Ind le polynôme kModèle:IndkModèle:Ind :

$\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}$

et

$\frac{h_n(t)}{k_n(t)}=h_0(t)+\sum _{i=1}^n\frac{a_{{\nu }_n}}{d_n(t)}{t}^{{\nu }_n}.$

Notons D le disque de convergence de la série entière de terme général a_kt^k (en posant aModèle:Ind = 0 si k n'est pas égal à l'un des ν_n).

• Pour N et c > 0 fixés, soit H l'ensemble des t pour lesquels, à partir de l'indice N, le module de dModèle:Ind(t) reste supérieur à c. Sur l'intersection de H et du disque ouvert D, la fraction continue converge (elle converge même uniformément sur tout compact de H⋂D).
• En tout point t extérieur au disque fermé Modèle:Surligner et pour lequel la suite (dModèle:Ind(t)) reste bornée, la fraction continue diverge.

Illustration par l'exemple de la fonction tangente

Modifions un peu les conventions du § « Fonction tangente » ci-dessus pour rendre égal à 1 le terme constant de chaque dénominateur k_n :

$\mathrm{\forall }n\ge 2h_n(t)=h_{n-1}(t)-\frac{t^2}{(2n-1)(2n-3)}{h}_{n-2}(t)\text{et}{k}_n(t)=k_{n-1}(t)-\frac{t^2}{(2n-1)(2n-3)}{k}_{n-2}(t).$

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 :

$\mathrm{\forall }n\in \mathbb{N}{a}_{2n}=\prod _{i=1}^n\frac{1}{2i-1}\text{et}{a}_{2n+1}=0.$

Il suffit de minorer 2i – 1 par i pour conclure que la série entière de coefficients aModèle:Ind est majorée terme à terme par exp(t²). Son rayon de convergence est donc infini.

La suite des dénominateurs de la fraction continue converge manifestement vers la fonction cosinus.

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 (Modèle: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_n(t)) converge donc vers cosModèle:2(t).

Soit K un fermé du plan complexe, borné par un réel M, et ne contenant aucun des zéros de la fonction cosinus (qui sont les (2j+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 2c telle que le module de cosModèle:2(t) soit supérieur à 2c pour tout t de K donc il existe un rang N tel que toute fonction de la suite (d_n(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 non seulement cette convergence des dénominateurs vers cosinus mais aussi des numérateurs vers sinus.

Voir aussi

Article connexe

Fonction de Bessel

Bibliographie

• Jean-Étienne Rombaldi, Interpolation & approximation, analyse pour l'agrégation, Vuibert, 2005
• Modèle:Ouvrage
• Modèle:En C. Brezinski et M. Redivo-Zaglia, Extrapolation Methods: Theory and Practice, North-Holland, 1991 Modèle:ISBN
• Modèle:En G. A. Baker et P. Graves-Morris, Padé Approximants, coll. « Encyclopedia of Mathematics and its Applications » (n° 59), 1996, Modèle:2e éd. Modèle:Lire en ligne Modèle:ISBN
• H. Padé, Œuvres rassemblées et présentées par Claude Brezinski, Librairie Albert Blanchard, Paris, 1984.

Liens externes

• D. Vekemans, Approximants de Padé, Laboratoire de mathématiques pures et appliquées Joseph Liouville
• Modèle:MathWorld