ATLAS Ouvrir dans l'atlas →

ExerciceMaths Sup – MPSI, chapitre 10 — Analyse asymptotique et développements limités · Exercices de synthèse

Exercice 51

DémonstrationCalcul

Une approximation rationnelle du sinus. Pour tout couple (a,b)(a, b) de réels, on considère la fonction fa,bf_{a,b}, définie au voisinage de 00 par fa,b(x)=sin⁡x−x+ax31+bx2.f_{a,b}(x) = \sin x - \frac{x + ax^{3}}{1 + bx^{2}}.

1. Déterminer le développement limité de fa,bf_{a,b} à l'ordre 66 en 00.

2. Montrer qu'il existe un unique couple (a0,b0)(a_{0}, b_{0}) tel que fa0,b0(x)=o(x6)f_{a_{0},b_{0}}(x) = o(x^{6}) au voisinage de 00, et donner un équivalent de fa0,b0f_{a_{0},b_{0}} en 00.

3. Parmi les fonctions fa,bf_{a,b}, y en a-t-il une qui soit négligeable en 00 devant toutes les autres ?

4. On pose R(x)=x+a0x31+b0x2R(x) = \dfrac{x + a_{0}x^{3}}{1 + b_{0}x^{2}} et T(x)=x−x36+x5120T(x) = x - \dfrac{x^{3}}{6} + \dfrac{x^{5}}{120} : chacune de ces deux approximations du sinus dépend de trois coefficients. Montrer que R(x)<sin⁡x<T(x)R(x) < \sin x < T(x) pour tout x>0x > 0 assez petit, puis comparer RR et TT au voisinage de +∞+\infty.

Corrigé

1. Comme bx2bx^{2} tend vers 00, le développement 11+y=1−y+y2−y3+o(y3)\frac{1}{1 + y} = 1 - y + y^{2} - y^{3} + o(y^{3}), composé avec y=bx2y = bx^{2}, donne 11+bx2=1−bx2+b2x4+o(x5)\frac{1}{1 + bx^{2}} = 1 - bx^{2} + b^{2}x^{4} + o(x^{5}), le terme en x6x^{6} étant absorbé par le reste. Ainsi x+ax31+bx2=(x+ax3)(1−bx2+b2x4+o(x5))=x+(a−b)x3+(b2−ab)x5+o(x6),\frac{x + ax^{3}}{1 + bx^{2}} = (x + ax^{3})\big(1 - bx^{2} + b^{2}x^{4} + o(x^{5})\big) = x + (a - b)x^{3} + (b^{2} - ab)x^{5} + o(x^{6}), le terme ab2x7ab^{2}x^{7} étant absorbé par le reste. Avec sin⁡x=x−x36+x5120+o(x6)\sin x = x - \frac{x^{3}}{6} + \frac{x^{5}}{120} + o(x^{6}) (théorème 5), fa,b(x)=(b−a−16)x3+(1120+ab−b2)x5+o(x6)\boxed{f_{a,b}(x) = \Big(b - a - \frac{1}{6}\Big)x^{3} + \Big(\frac{1}{120} + ab - b^{2}\Big)x^{5} + o(x^{6})} L'absence de termes de degré pair était prévisible : fa,bf_{a,b} est impaire (proposition 9).

2. Par unicité du développement limité (proposition 8), fa,b(x)=o(x6)f_{a,b}(x) = o(x^{6}) si et seulement si les deux coefficients sont nuls, c'est-à-dire si b−a=16b - a = \frac{1}{6} et b(b−a)=1120b(b - a) = \frac{1}{120}. La première égalité, reportée dans la seconde, donne b6=1120\frac{b}{6} = \frac{1}{120}, d'où un unique couple : a0=−760,b0=120\boxed{a_{0} = -\frac{7}{60}, \qquad b_{0} = \frac{1}{20}} Pour l'équivalent, on pousse le calcul d'un ordre. Le développement 11+b0x2=1−b0x2+b02x4−b03x6+o(x7)\frac{1}{1 + b_{0}x^{2}} = 1 - b_{0}x^{2} + b_{0}^{2}x^{4} - b_{0}^{3}x^{6} + o(x^{7}) ajoute au quotient le terme (a0b02−b03)x7=b02(a0−b0)x7=−x72400(a_{0}b_{0}^{2} - b_{0}^{3})x^{7} = b_{0}^{2}(a_{0} - b_{0})x^{7} = -\frac{x^{7}}{2400}, et le sinus contribue −x75040-\frac{x^{7}}{5040}. Comme −15040+12400=−10+2150 400-\frac{1}{5040} + \frac{1}{2400} = \frac{-10 + 21}{50\,400}, la proposition 13 donne fa0,b0(x)∼x→01150 400 x7\boxed{f_{a_{0},b_{0}}(x) \underset{x \to 0}{\sim} \frac{11}{50\,400}\,x^{7}}

3. Soit (a,b)≠(a0,b0)(a, b) \neq (a_{0}, b_{0}). Par la question 2, l'un des coefficients de x3x^{3} ou de x5x^{5} dans le développement de fa,bf_{a,b} est non nul, et la proposition 13 donne fa,b(x)∼c xpf_{a,b}(x) \sim c\,x^{p}, avec c≠0c \neq 0 et p∈{3,5}p \in \{3, 5\}. Comme x7=o(xp)x^{7} = o(x^{p}), on obtient fa0,b0=o(fa,b)f_{a_{0},b_{0}} = o(f_{a,b}). Aucune autre fonction ne convient : si fa,bf_{a,b}, avec (a,b)≠(a0,b0)(a, b) \neq (a_{0}, b_{0}), était négligeable devant fa0,b0f_{a_{0},b_{0}}, elle serait négligeable devant elle-même par transitivité, donc nulle au voisinage de 00, ce que contredit fa,b(x)∼c xpf_{a,b}(x) \sim c\,x^{p}. oui : fa0,b0 est la seule, avec x+a0x31+b0x2=x(60−7x2)3(20+x2)\boxed{\text{oui : } f_{a_{0},b_{0}} \text{ est la seule, avec } \frac{x + a_{0}x^{3}}{1 + b_{0}x^{2}} = \frac{x(60 - 7x^{2})}{3(20 + x^{2})}}

4. Le développement du sinus à l'ordre 88 donne sin⁡x−T(x)∼−x75040\sin x - T(x) \sim -\frac{x^{7}}{5040}, et la question 2 donne sin⁡x−R(x)∼11x750 400\sin x - R(x) \sim \frac{11x^{7}}{50\,400}. Par la proposition 13, ces deux différences ont respectivement le signe de −x7-x^{7} et de x7x^{7} au voisinage de 00 : R(x)<sin⁡x<T(x)pour tout x>0 assez petit\boxed{R(x) < \sin x < T(x) \quad\text{pour tout } x > 0 \text{ assez petit}} Au voisinage de +∞+\infty, une fonction rationnelle est équivalente au quotient de ses termes de plus haut degré : R(x)∼−7x3/60x2/20R(x) \sim \frac{-7x^{3}/60}{x^{2}/20} et T(x)∼x5120T(x) \sim \frac{x^{5}}{120}, soit R(x)∼x→+∞−7x3,T(x)∼x→+∞x5120\boxed{R(x) \underset{x \to +\infty}{\sim} -\frac{7x}{3}, \qquad T(x) \underset{x \to +\infty}{\sim} \frac{x^{5}}{120}} Aucune des deux ne reste bornée comme le sinus, mais RR s'en écarte bien plus lentement. Numériquement, R(1)=5363≈0,841 270R(1) = \frac{53}{63} \approx 0{,}841\,270 et T(1)=101120≈0,841 667T(1) = \frac{101}{120} \approx 0{,}841\,667 encadrent sin⁡1≈0,841 471\sin 1 \approx 0{,}841\,471 avec des écarts presque égaux, mais R(π)≈−0,319R(\pi) \approx -0{,}319 est plus proche de sin⁡π=0\sin\pi = 0 que T(π)≈0,524T(\pi) \approx 0{,}524.

Piège. Être négligeable devant toutes les autres est une propriété locale, qui ne dit pas que RR est la meilleure approximation. Près de 00, c'est même TT qui l'emporte de peu, puisque 15040≈1,98⋅10−4\frac{1}{5040} \approx 1{,}98 \cdot 10^{-4} est inférieur à 1150 400≈2,18⋅10−4\frac{11}{50\,400} \approx 2{,}18 \cdot 10^{-4} ; l'avantage de RR n'apparaît que loin de 00. Ces fractions rationnelles, qui coïncident avec une fonction à un ordre aussi élevé que possible, ont été étudiées systématiquement par Henri Padé dans sa thèse de 1892, et portent son nom ; les approximations rationnelles de ce type servent au calcul numérique des fonctions usuelles.

Ouvrir dans l'atlas interactif →