M12 Mathématiques 2
M12 Mathématiques 2
M12 Mathématiques 2
Mathématiques 2
Licence Sciences et Techniques
Aide-mémoire.
L1 INFO MATH SI
Gloria Faccanoni
i http://faccanoni.univ-tln.fr/enseignements.html
1 Nombres complexes 5
2 Fonctions polynomiales 11
3 Primitives et Intégrales 29
L’objet de ce aide-mémoire est de proposer une explication succincte des concepts vu en cours. De nombreux livres, parfois
très fournis, existent. Ici on a cherché, compte tenu des contraintes de volume horaire, des acquis des étudiants au lycée et
des exigences pour la suite du cursus, à dégager les points clés permettant de structurer le travail personnel de l’étudiant
voire de faciliter la lecture d’autres ouvrages. Ce polycopiée ne dispense pas des séances de cours et de TD ni de prendre des
notes complémentaires. Il est d’ailleurs important de comprendre et apprendre le cours au fur et à mesure. Ce polycopié
est là pour éviter un travail de copie qui empêche parfois de se concentrer sur les explications données oralement mais ce
n’est pas un livre auto-suffisant (il est loin d’être exhaustif ) ! De plus, ne vous étonnez pas si vous découvrez des erreurs.
Malgré de très nombreuses relectures, il restera toujours des fautes, ce polycopié est donc fournit sans garanties ! N’hésitez
pas à me signaler les erreurs que vous remarquez.
On a inclus dans ce texte nombreux exercices corrigés. Ceux-ci, de difficulté variée, répondent à une double nécessitée.
Il est important de jongler avec les différents concepts introduits en cours et même de faire certaines erreurs une fois
pour bien identifier les pièges. Les exercices permettent d’orienter les raisonnements vers d’autres domaines (physique,
économie, etc.), cela afin d’exhiber l’intérêt et l’omniprésence du calcul différentiel. Cependant, veuillez noter que vous
n’obtiendrez pas grande chose si vous vous limitez à choisir un exercice, y réfléchir une minute et aller vite voir le début de
la correction en passant tout le temps à essayer de comprendre la correction qui va paraitre incompréhensible. Pour que la
méthode d’étude soit vraiment efficace, il faut d’abord vraiment essayer de chercher la solution. En particulier, il faut avoir
un papier brouillon à coté de soi et un crayon. La première étape consiste alors à traduire l’énoncé (pas le recopier), en
particulier s’il est constitué de beaucoup de jargon mathématique. Ensuite il faut essayer de rapprocher les hypothèses de
la conclusion souhaitée, et pour cela faire quelques calculs ou transformer les hypothèses pour appliquer un théorème
dont on aura vérifier que les hypothèses sont bien satisfaite. C’est ici que l’intuition joue un grand rôle et il ne faut pas
hésiter à remplir des pages pour s’apercevoir que l’idée qu’on a eu n’est pas la bonne. Elle pourra toujours resservir dans
une autre situation. Quand finalement on pense tenir le bon bout, il faut rédiger soigneusement en s’interrogeant à chaque
pas sur la validité (logique, mathématique) de ce qu’on a écrit. Si l’étape précédente ne donne rien, il faut chercher de l’aide
(voir le début de la correction, en parler à un autre étudiant, interroger les tuteurs).
Gloria FACCANONI
Intervalles
Inégalité(s) Ensemble Représentations graphique
a b a b
a ≤x ≤b [a, b]
a b a b
a <x <b ]a, b[
a b a b
a ≤x <b [a, b[
a b a b
a <x ≤b ]a, b]
a a
x≥a [a, +∞[
a a
x>a ]a, +∞[
b b
x ≤b ] − ∞, b]
b b
x <b ] − ∞, b[
−a a −a a
|x| ≤ a avec a ≥ 0 [−a, a]
−a a −a a
|x| < a avec a ≥ 0 ] − a, a[
−a a −a a
|x| ≥ a avec a ≥ 0 ] − ∞, −a] ∪ [a, +∞[
−a a −a a
|x| > a avec a ≥ 0 ] − ∞, −a[∪]a, +∞[
∀x ∈ R ] − ∞, +∞[
a a
x 6= a ] − ∞, a[∪]a, +∞[= R \ { a }
exercice de base
exercice
exemple
curiosité
3
Notations Mercredi 4 septembre 2019
4 © 2019-2020 G. Faccanoni
Nombres complexes
1
Forme algébrique Tout nombre z ∈ C s’écrit, de manière unique, sous la forme algébrique z = x + i y où x, y ∈ R et i est tel
que i 2 = −1. Le réel x s’appelle “partie réelle de z” et se note ℜ(z). Le réel y s’appelle “partie imaginaire de z” et se
note ℑ(Z ). Si y = 0 alors z ∈ R, si x = 0 alors z est un imaginaire pur.
Égalité Deux nombres complexes z = x + i y et z 0 = x 0 + i y 0 sont égaux ssi (x, y) = (x 0 , y 0 ).
Attention : il n’y a pas d’inégalités en C.
Opérations Soient deux nombres complexes z = x + i y et z 0 = x 0 + i y 0 .
z + z 0 = (x + x 0 ) + i (y + y 0 ), zz 0 = (xx 0 − y y 0 ) + i (x y 0 + x 0 y).
y
Argument On le note arg(z) et il est défini, modulo 2π, par cos ϑ = %x et sin ϑ = % . On a (à 2kπ près, k ∈ Z)
Forme trigonométrique Tout nombre complexe non nul z s’écrit sous forme trigonométrique z = %(cos ϑ + sin ϑ) avec
% > 0 le module de z et ϑ un argument de z.
Forme exponentielle Tout nombre complexe non nul z s’écrit sous forme exponentielle z = %e i ϑ .
? Formule de Moivre : pour tout ϑ ∈ R et pour tout k ∈ Z
e i x + e −i x e i x − e −i x
cos(x) = , sin(x) = ,
2 2i
e i nx + e −i nx e i nx − e −i nx
cos(nx) = , sin(nx) = .
2 2i
5
1 Nombres complexes Mercredi 4 septembre 2019
z = %e i ϑ = x + i y
y
%
ϑ
−x % x ℜ
−ϑ
−z = %e i (ϑ+π) = −x − i y −y z̄ = %e −i ϑ = x − i y
Racines n -ièmes Tout nombre complexe non nul z = %e i ϑ possède n racines n-ièmes
p ϑ+2kπ
wk = n
%e i n avec k ∈ {0, 1, . . . , n − 1}.
ℑ ℑ
z
z
ρ
p
ρ
w1 ϑ
w1 ϑ w2
pρ w0 p
6 ρ
w0
ϑ/3 ϑ/6
ℜ w3 ℜ
w5
w2 w4
6 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 1 Nombres complexes
ℑ ℑ
y z z
ϑ ϑ
pr w0 w̄ 1 w0
β
−α ϑ/2 ϑ/2
α x ℜ ℜ
−ϑ/2
−β
w1 w1 w̄ 0
−ϑ
z̄
1.1 Proposition
L’équation az 2 +bz +c = 0, où a 6= 0, b et c appartiennent à C, admet deux solutions dans C (pas nécessairement distinctes) :
−b + δ −b − δ
z1 = , z2 = ,
2a 2a
En effet
b 2
¶ µ ¶2
b c b c
µ ¶ ·µ ¸
0 = az 2 + bz + c = a z 2 + z + =a z+ − +
a a 2a 2a a
donc ¶2 ¶2
b b c b 2 − 4ac
µ µ
z+ = − = .
2a 2a a 4a 2
Exercice 1.1
Calculer les racines carrées de
a) 1 b) i c) 3 + 4i d) 8 − 6i
1+i
e) 7 + 24i f) 3 − 4i g) 24 − 10i h) p
2
Correction (
α=1
(a) z = 1 + 0i =⇒ x = 1, y = 0 =⇒ =⇒ w 0 = 1 et w 1 = −1
β=0
α = p1
(
(b) z = i =⇒ x = 0, y = 1 =⇒ 2 =⇒ w 0 = p1 + i p1 = e i π/4 et w 1 = − p1 − i p1 = e i 5π/4
β = p1 2 2 2 2
2
(
α=2
(c) z = 3 + 4i =⇒ x = 3, y = 4 =⇒ =⇒ w 0 = 2 + i et w 1 = −2 − i
β=1
(
α=3
(d) z = 8 − 6i =⇒ x = 8, y = −6 =⇒ =⇒ w 0 = 3 − i et w 1 = −3 + i : on aurait pu tout simplement remarquer
β = −1
que z est le conjugué du complexe au point (c)
(
α=4
(e) z = 7 + 24i =⇒ x = 7, y = 24 =⇒ =⇒ w 0 = 4 + 3i et w 1 = −4 − 3i
β=3
(
α=2
(f) z = 3 − 4i =⇒ x = 3, y = −4 =⇒ =⇒ w 0 = 2 − i et w 1 = −2 + i
β = −1
© 2019-2020 G. Faccanoni 7
1 Nombres complexes Mercredi 4 septembre 2019
(
α=5
(g) z = 24 − 10i =⇒ x = 24, y = −10 =⇒ =⇒ w 0 = 5 − i et w 1 = −5 + i
β = −1
(h) z = 1+i
p = 12 e i π/4 =⇒ w 0 = p1 e i π/8 et w 1 = p1 e i (π/8+π)
2 2 2
p p p p
(i) z = 2 = 2e i 0 =⇒ w 0 = 2e i 0 = 2 et w 1 = 2e i π = − 2
p p
(j) z = 3i = 3e i π =⇒ w 0 = 3e i π/4 et w 1 = 3e i 5π/4
q p
α = 2+ 13 q p q p
(k) z = 2 + 3i =⇒ x = 2, y = 3 =⇒ q 2p =⇒ w 0 = 2+ 2 13 + i −2+2 13 et w 1 = −w 0
β = −2+2 13
q p
α = 1+ 1+b 2
(l) z = 1 + bi (b est un paramètre réel) =⇒ x = 1, y = b =⇒ 2 q =⇒
p
2
β = signe(b) −1+ 21+b
q p q p
1+ 1+b 2 −1+ 1+b 2
w = + i si b > 0,
0
q p2 2
p
q
1+ 1+b 2 −1+ 1+b 2
w 0 = 2 − i 2 si b < 0,
w = 1 si b = 0,
0
et w 1 = −w 0
Exercice 1.2
Trouver (sous forme algébrique et trigonométrique) toutes les racines complexes de l’équation
z 2 − (1 + 2i )z + i − 1 = 0.
Correction
−b+δ0
A = 1, b = −1 − 2i , c = i − 1, ∆ = b 2 − 4ac = (−(1 + 2i ))2 − 4(i − 1) = 1 donc δ0 = 1 et δ1 = −δ0 et finalement z 0 = 2 =
p ¡
1 + i = 2 cos π4 + i sin π4 et z 1 = −b−δ = i = cos π2 + i sin π2 .
¡ ¢ ¡ ¢¢ 0
¡ ¢ ¡ ¢
2
Exercice 1.3
Trouver toutes les solutions complexes des équations
a) z 2 − 2i z + 2 − 4i = 0, b) z 2 + 2z + (1 − 2i ) = 0, c) z 2 + 2(1 + i )z − 5(1 + 2i ) = 0,
p
d) z 2 − (5 + 3i )z + 7i + 4 = 0, e) z 2 − 3z − i = 0, f ) z 2 + z + 1 = 0,
j) 4z 2 − 2z + 1 = 0, k) z 2 − (11 − 5i )z + 24 − 27i = 0, l) z 4 + 2z 2 + 4 = 0.
Correction
2i +2+4i
1. z 2 −2i z+2−4i = 0 =⇒ δ2 = (−2i )2 −4(2−4i ) = −12+16i =⇒ δ0 = 2+4i =⇒ z 0 = 2 = 1+3i et z 1 = 2i −2−4i
2 = −1−i
−2+2+2i −2−2−2i
2. z 2 + 2z + (1 − 2i ) = 0 =⇒ δ2 = 22 − 4(1 − 2i ) = 8i =⇒ δ0 = 2 + 2i =⇒ z 0 = 2 = i et z 1 = 2 = −2 − i
−2−2i +6+4i
3. z 2 + 2(1 + i )z − 5(1 + 2i ) = 0 =⇒ δ2 = (2(1 + i ))2 − 4(−5(1 + 2i )) = 20 + 48i =⇒ δ0 = 6 + 4i =⇒ z 0 = 2 = 2+i
et z 1 = −2−2i2−6−4i = −4 − 3i
4. z 2 −(5+3i )z +7i +4 = 0 =⇒ δ2 = (−(5+3i ))2 −4(4+7i ) = 2i =⇒ δ0 = 1+i =⇒ z 0 = 5+3i2+1+i = 3+2i et z 1 = 5+3i −1−i
2 =
2+i
p p p p
5. z 2 − 3z − i = 0 =⇒ δ2 = (− 3)2 − 4(−i ) = 3 + 4i =⇒ δ0 = 2 + i =⇒ z 0 = 3+2+i
2 et z 1 = 3−2−i
2
p p p
2 2 −1+i 3 −1−i 3
6. z + z + 1 = 0 =⇒ δ = −3 =⇒ δ0 = i 3 =⇒ z 0 = 2 et z 1 = 2
1+2i −1 1+2i +1
7. z 2 − (1 + 2i )z + i − 1 = 0 =⇒ δ2 = (−(1 + 2i ))2 − 4(i − 1) = 1 =⇒ δ0 = 1 =⇒ z 0 = 2 = i et z 1 = 2 = 1+i
2 2
8. z − (5 − 14i )z − 2(5i + 12) = 0 =⇒ δ = −75 − 100i =⇒ δ0 = 5 − 10i =⇒ z 0 = 5 − 12i et z 1 = −2i
9. z 2 − (3 + 4i )z − 1 + 5i = 0 =⇒ δ2 = −3 + 4i =⇒ δ0 = 1 + 2i =⇒ z 0 = 2 + 3i et z 1 = 1 + i
p p p
10. 4z 2 − 2z + 1 = 0 =⇒ δ2 = −12 =⇒ δ0 = 2 3i =⇒ z 0 = 1−4 3i et z 1 = 1+4 3i
11. z 2 − (11 − 5i )z + 24 − 27i = 0 =⇒ δ2 = −2i =⇒ δ0 = 1 − i =⇒ z 0 = 6 − 3i et z 1 = 5 − 2i
8 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 1 Nombres complexes
p p p
12. On pose w = z 2 et on résout d’abord w 2 +2w +4 = 0 =⇒ δ2 = −12 =⇒ δ0 = 2 3i =⇒ w 0 = −1+ 3i et w 1 = −1− 3i .
On résout maintenant les deux équations
α = p1
(
p p
2
? z = w 0 =⇒ p2 =⇒ z 0 = p1 + i p3 et z 1 = − p1 − i p3
β = p3 2 2 2 2
2
α = p1
(
p p
2
? z = w 1 =⇒ 2p
3
=⇒ z 3 = p1 − i p3 et z 4 = − p1 + i p3
β = −p 2 2 2 2
2
© 2019-2020 G. Faccanoni 9
Fonctions polynomiales
2
2.1 Définition (Fonctions polynomiales)
Une application P de C dans C est polynomiale s’il existe une famille finie {a n , . . . , a 1 , a 0 } de nombres complexes telle que
pour tout x ∈ C,
P (x) = a 0 + a 1 x + · · · + a n x n .
Pour faire bref, nous ferons un abus de langage en utilisant le terme “polynôme” pour “application polynomiale”.
Remarque
Le polynôme nul n’a pas de degré, mais par convention do 0 < do g pour tout polynôme non nul g .
On a donc n 1 + · · · + n k = do P .
n (x − x 0 )k
P (k) (x 0 )
X
P (x) =
k=0 k!
11
2 Fonctions polynomiales Mercredi 4 septembre 2019
2.8 Théorème
Soient P un polynôme et x 0 une racine de P . Elle est de multiplicité n si et seulement si
Astuce
Pour factoriser les polynômes on pourra utiliser Wolfram Alpha www.wolframalpha.com et taper factor(P (x)).
Correction
1.
7x 4 −x 3 +2x −4 x2 −3x +5
25x 2 −98x −4
−23x −129
−x 5 +2x 4 −x 3 x3 +x 2 +x +1
x4 −x 3 −x +1
−x 4 +2x 3 −x 2
x3 −x 2 −x +1
−x 3 +2x 2 −x
x2 −2x +1
−x 2 +2x −1
12 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 2 Fonctions polynomiales
Correction
1. On remarque que z 3 − i est factorisable par (z + i ) car z 3 − i = z 3 + i 3 = (z + i )(z 2 − 2i z − 1) donc
Exercice 2.3
Soit p un polynôme de degré 4. On sait qu’il assume en x = 2 et x = 3 sa valeur maximale 3 et qu’il s’annule en x = 1. Calculer
p(4).
Correction
? p(2) = p(3) = 3 et p a degré 4 donc p(x) − 3 = (x − 2)(x − 3)q(x) avec q un polynôme de degré 2.
? x = 2 et x = 3 sont des maxima pour p, donc p 0 (2) = p 0 (3) = 0. Comme p 0 (x) = (2x − 5)q(x) + (x − 2)(x − 3)q 0 (x), alors
q(2) = 0 et q(3) = 0, ce qui implique q(x) = a(x − 2)(x − 3). On obtient alors p(x) = 3 + a(x − 2)2 (x − 3)2 .
? Comme p(1) = 0 alors a = −3/4 et on conclut que p(4) = 0.
Exercice 2.4
Soit n ∈ N∗ . Démontrer que P = x n+1 + x n − 2x n−1 + nx − n est divisible par (x − 1).
Correction
Un polynôme P est divisible par (x − 1) si x = 1 est racine de P , donc il suffit de vérifier que P (1) = 0.
Exercice 2.5
Soit n ∈ N∗ . Démontrer que P = (x − 2)2n + (x − 1)n − 1 est divisible par x 2 − 3x + 2.
Correction
Comme x 2 − 3x + 2 = (x − 1)(x − 2), il suffit de montrer que x = 1 et x = 2 sont racines du polynôme P , autrement dit que
P (1) = P (2) = 0.
Exercice 2.6
Soit n ∈ N, n ≥ 2. On considère le polynôme P défini par P = ax n+1 + bx n + 1. Déterminer les réels a et b pour que P soit
divisible par (x − 1)2 .
Correction
Si P est divisible par (x − 1)2 , cela signifie que x = 1 est racine au moins double, donc il faut que P (1) = P 0 (1) = 0. Comme
P (1) = a + b + 1 et P 0 (1) = a(n + 1) + bn, on conclut que a = n et b = −n − 1.
© 2019-2020 G. Faccanoni 13
2 Fonctions polynomiales Mercredi 4 septembre 2019
Correction
? x 0 = 1 a multiplicité 3 car
P (x) = x 4 − x 3 − 3x 2 + 5x − 2 P (1) = 0
0 3 2
P (x) = 4x − 3x − 6x + 5 P 0 (1) = 0
P 00 (x) = 12x 2 − 6x − 6 P 00 (1) = 0
P 000 (x) = 24x − 6 P 000 (1) = 18
P (x) = x 3 − i x 2 + x − i P (1) = 0
0 2
P (x) = 3x − 2i x + 1 P 0 (1) = 0
P 00 (x) = 6x − 2i P 00 (1) = 4i
Exercice 2.8
On considère le polynôme
P (x) := x 5 − 2x 4 + x 3 + x 2 − 2x + 1.
1. Écrire la formule de Taylor de P en 1. Montrer que 1 est une racine de P et en déterminer son ordre de multiplicité.
2. Factoriser P sur R.
3. Factoriser P sur C.
Correction
1. Pour un polynôme P de degré 5, la formule de Taylor en 1 s’écrit
k=5 (x − 1)k
P (k) (1)
X
P (x) = .
k=0 k!
P (x) = x 5 − 2x 4 + x 3 + x 2 − 2x + 1 P (1) = 0
0 4 3 2
P (x) = 5x − 8x + 3x + 2x − 2 P 0 (1) = 0
P 00 (x) = 20x 3 − 24x 2 + 6x + 2 P 00 (1) = 4
P 000 (x) = 60x 2 − 48x + 6 P 000 (1) = 18
P I V (x) = 120x − 48 P I V (1) = 72
P V (x) = 120 P V (1) = 120
14 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 2 Fonctions polynomiales
2. On pose Q(x) = x 3 + 1. Une racine réelle évidente de Q est −1. En effectuant la division euclidienne de Q par (x + 1)
on obtient
Q(x) = (x + 1)(x 2 − x + 1).
Étant donné que x 2 − x + 1 ne se factorise pas sur R, on conclut que la factorisation de P sur R est
3. Pour factoriser P sur C il ne reste que factoriser le polynôme x 2 − x + 1 sur C. Avec les notations usuelles
p
du cours,
p 1 3 π
on a ∆ = −3 et une racine carrée de ∆ est δ = i 3 d’où les deux racines complexes conjuguées x 1 = 2 − i 2 = e i 3 et
p 5π
3
x 2 = 12 + i 2 = ei 3 . On conclut que la factorisation de P sur C est
π 5π
P (x) = (x − 1)2 (x + 1)(x − e i 3 )(x − e i 3 ).
Exercice 2.9
On considère le polynôme
P (x) := x 4 − x 3 − 3x 2 + 5x − 2.
1. Montrer que x = 1 est une racine du polynôme et en déterminer son ordre de multiplicité.
2. Factoriser P à l’aide de la formule de Taylor.
Correction
1. On a
P (x) = x 4 − x 3 − 3x 2 + 5x − 2 P (1) = 0
0 3 2
P (x) = 4x − 3x − 6x + 5 P 0 (1) = 0
P 00 (x) = 12x 2 − 6x − 6 P (1) = 0
000 000
P (x) = 24x − 6 P (1) = 18
IV
P (x) = 24 P I V (1) = 24
k=4 (x − 1)k
P (k) (1)
X
P (x) = .
k=0 k!
(x − 1)3 (x − 1)4 ³ ´
P (x) = 18 + 24 = 3(x − 1)3 + (x − 1)4 = (x − 1)3 x + 2 .
3! 4!
Exercice 2.10
On considère le polynôme
P (x) = x 5 − x 4 − 2x 3 + 5x 2 − 5x + 2.
1. Écrire la formule de Taylor de P en 1 ; montrer que 1 est une racine de P et déterminer son ordre de multiplicité.
2. Factoriser P sur R en sachant qu’il possède une autre racine réelle évidente (i.e. elle est un diviseur du terme de degré
0).
3. Factoriser P sur C.
Correction
1. Pour un polynôme P de degré 5, la formule de Taylor en 1 s’écrit
k=5 (x − 1)k
P (k) (1)
X
P (x) = .
k=0 k!
P (x) = x 5 − x 4 − 2x 3 + 5x 2 − 5x + 2 P (1) = 0
© 2019-2020 G. Faccanoni 15
2 Fonctions polynomiales Mercredi 4 septembre 2019
3. Pour factoriser P sur C il ne reste que factoriser le polynôme x 2 − x + 1 sur C. Avec les notations usuelles
p
du cours,
p π
on a ∆ = −3 et une racine carrée de ∆ est δ = i 3 d’où les deux racines complexes conjuguées x 1 = 12 − i 23 = e i 3 et
p 5π
3
x 2 = 12 + i 2 = ei 3 . On conclut que la factorisation de P sur C est
π 5π
P (x) = (x − 1)2 (x + 2)(x − e i 3 )(x − e i 3 ).
Exercice 2.11
On considère le polynôme
P (x) = x 5 − x 4 − 2x 3 + x 2 − x − 2.
1. Écrire la formule de Taylor de P en −1 ; montrer que −1 est une racine de P et déterminer son ordre de multiplicité.
2. Factoriser P sur R en sachant qu’il possède une autre racine réelle évidente (i.e. elle est un diviseur du terme de degré
0).
3. Factoriser P sur C.
Correction
1. Pour un polynôme P de degré 5, la formule de Taylor en −1 s’écrit
k=5 (x + 1)k
P (k) (−1)
X
P (x) = .
k=0 k!
P (x) = x 5 − x 4 − 2x 3 + x 2 − x − 2 P (−1) = 0
0 4 3 2
P (x) = 5x − 4x − 6x + 2x − 1 P 0 (−1) = 0
P 00 (x) = 20x 3 − 12x 2 − 12x + 2 P 00 (−1) = −18
P 000 (x) = 60x 2 − 24x − 12 P 000 (−1) = 72
P I V (x) = 120x − 24 P I V (−1) = −144
P V (x) = 120 P V (−1) = 120
16 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 2 Fonctions polynomiales
3. Pour factoriser P sur C il ne reste que factoriser le polynôme x 2 − x + 1 sur C. Avec les notations usuelles
p
du cours,
p 1 3 π
on a ∆ = −3 et une racine carrée de ∆ est δ = i 3 d’où les deux racines complexes conjuguées x 1 = 2 − i 2 = e i 3 et
p 5π
3
x 2 = 12 + i 2 = ei 3 . On conclut que la factorisation de P sur C est
π 5π
P (x) = (x + 1)2 (x − 2)(x − e i 3 )(x − e i 3 ).
Exercice 2.12
On considère le polynôme
P (x) = x 5 − 5x 4 + 10x 3 − 11x 2 + 7x − 2.
1. Écrire la formule de Taylor de P en 1 ; montrer que 1 est une racine de P et déterminer son ordre de multiplicité.
2. Factoriser P sur R en sachant qu’il possède une autre racine réelle évidente (i.e. elle est un diviseur du terme de degré
0).
3. Factoriser P sur C.
Correction
1. Pour un polynôme P de degré 5, la formule de Taylor en 1 s’écrit
k=5 (x − 1)k
P (k) (1)
X
P (x) = .
k=0 k!
(x − 1)2 (x − 1)5
P (x) = −2 + 120
2³ 120 ´
= (x − 1)2 x 3 − 3x 2 + 3x − 2
© 2019-2020 G. Faccanoni 17
2 Fonctions polynomiales Mercredi 4 septembre 2019
3. Pour factoriser P sur C il ne reste que factoriser le polynôme x 2 − x + 1 sur C. Avec les notations usuelles
p
du cours,
p π
on a ∆ = −3 et une racine carrée de ∆ est δ = i 3 d’où les deux racines complexes conjuguées x 1 = 12 − i 23 = e i 3 et
p 5π
3
x 2 = 12 + i 2 = ei 3 . On conclut que la factorisation de P sur C est
π 5π
P (x) = (x − 1)2 (x − 2)(x − e i 3 )(x − e i 3 ).
Exercice 2.13
On considère le polynôme
P (x) = x 6 − 5x 5 + 10x 4 − 12x 3 + 11x 2 − 7x + 2.
1. Écrire la formule de Taylor de P en 1 ; montrer que 1 est une racine de P et déterminer son ordre de multiplicité.
2. Factoriser P sur R en sachant qu’il possède une autre racine réelle évidente (i.e. elle est un diviseur du terme de degré
0).
3. Factoriser P sur C.
Correction
1. Pour un polynôme P de degré 6, la formule de Taylor en 1 s’écrit
k=6 (x − 1)k
P (k) (1)
X
P (x) = .
k=0 k!
3. Pour factoriser P sur C il ne reste que factoriser le polynôme x 2 + 1 sur C ; on a les deux racines complexes conjuguées
x 1 = −i et x 2 = +i . On conclut que la factorisation de P sur C est
Exercice 2.14
On considère le polynôme
P (x) = x 8 + 5x 7 + 9x 6 + 7x 5 + 3x 4 + 5x 3 + 9x 2 + 7x + 2.
1. Écrire la formule de Taylor de P en −1 ; montrer que −1 est une racine de P et déterminer son ordre de multiplicité.
18 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 2 Fonctions polynomiales
2. Factoriser P sur R en sachant qu’il possède une autre racine réelle évidente (i.e. elle est un diviseur du terme de degré
0).
3. Factoriser P sur C.
Correction
1. Pour un polynôme P de degré 8, la formule de Taylor en −1 s’écrit
k=8 (x + 1)k
P (k) (−1)
X
P (x) = .
k=0 k!
P (x) = x 8 + 5x 7 + 9x 6 + 7x 5 + 3x 4 + 5x 3 + 9x 2 + 7x + 2 P (−1) = 0
0 7 6 5 4 3 2
P (x) = 8x + 35x + 54x + 35x + 12x + 15x + 18x + 7 P 0 (−1) = 0
P 00 (x) = 56x 6 + 210x 5 + 270x 4 + 140x 3 + 36x 2 + 30x + 18 P 00 (−1) = 0
P 000 (x) = 336x 5 + 1050x 4 + 1080x 3 + 420x 2 + 72x + 30 P 000 (−1) = 12
P I V (x) = 1680x 4 + 4200x 3 + 3240x 2 + 840x + 72 P I V (−1) = −48
P V (x) = 6720x 3 + 12600x 2 + 6480x + 840 P V (−1) = 240
P V I (x) = 20160x 2 + 25200x + 6480 P V I (−1) = 1440
P V I I (x) = 40320x + 25200 P V I I (−1) = −15120
P V I I I (x) = 40320 P V I I I (−1) = 40320
P (x) = 2(x + 1)3 − 2(x + 1)4 + 2(x + 1)5 + 2(x + 1)6 − 3(x + 1)7 + (x + 1)8 =
= (x + 1)3 2 − 2(x + 1) + 2(x + 1)2 + 2(x + 1)3 − 3(x + 1)4 + (x + 1)5 =
£ ¤
= (x + 1)3 (x 5 + 2x 4 + x + 2)
1+i 1−i
3. Pour factoriser P sur C il ne reste que factoriser le polynôme x 4 + 1 sur C ; on a les quatre racines x 1 = p , x2 = p ,
2 2
−1+i −1−i
x3 = p et x 4 = p . On conclut que la factorisation de P sur C est
2 2
1+i 1−i −1 + i −1 − i
µ ¶µ ¶µ ¶µ ¶
P (x) = (x + 1)3 (x + 2) x − p x− p x− p x− p .
2 2 2 2
Exercice 2.15
On considère le polynôme
P (x) = x 6 − x 5 − x 4 + x 2 + x − 1.
1. Écrire la formule de Taylor de P en 1 ; montrer que 1 est une racine de P et déterminer son ordre de multiplicité.
2. Factoriser P sur R en sachant qu’il possède une autre racine réelle évidente (i.e. elle est un diviseur du terme de degré
0).
3. Factoriser P sur C.
© 2019-2020 G. Faccanoni 19
2 Fonctions polynomiales Mercredi 4 septembre 2019
Correction
1. Pour un polynôme P de degré 6, la formule de Taylor en 1 s’écrit
k=6 (x − 1)k
P (k) (1)
X
P (x) = .
k=0 k!
P (x) = x 6 − x 5 − x 4 + x 2 + x − 1 P (1) = 0
P 0 (x) = 6x 5 − 5x 4 − 4x 3 + 2x + 1 P 0 (1) = 0
P 00 (x) = 30x 4 − 20x 3 − 12x 2 + 2 P 00 (1) = 0
P 000 (x) = 120x 3 − 60x 2 − 24x P 000 (1) = 36
P I V (x) = 360x 2 − 120x − 24 P I V (1) = 216
P V (x) = 720x − 120 P V (1) = 600
P V I (x) = 720 P V I (1) = 720
3. Pour factoriser P sur C il ne reste que factoriser le polynôme x 2 + x + 1 sur C. Avec les notations usuelles du
p
cours, on
p 2
a ∆ = −3 et une racine carrée de ∆ est δ = i 3 d’où les deux racines complexes conjuguées x 1 = − 12 + i 23 = e i 3 π et
p 4
x 2 = − 21 − i 2
3
= e i 3 π . On conclut que la factorisation de P sur C est
2 4
P (x) = (x − 1)3 (x + 1)(x − e i 3 π )(x − e i 3 π ).
Exercice 2.16
On considère le polynôme
P (x) = x 6 + 3x 5 + 3x 4 − 3x 2 − 3x − 1.
1. Écrire la formule de Taylor de P en −1 ; montrer que −1 est une racine de P et déterminer son ordre de multiplicité.
2. Factoriser P sur R en sachant qu’il possède une autre racine réelle évidente (i.e. elle est un diviseur du terme de degré
0).
3. Factoriser P sur C.
Correction
1. Pour un polynôme P de degré 6, la formule de Taylor en −1 s’écrit
k=6 (x + 1)k
P (k) (−1)
X
P (x) = .
k=0 k!
P (x) = x 6 + 3x 5 + 3x 4 − 3x 2 − 3x − 1 P (−1) = 0
20 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 2 Fonctions polynomiales
3. Pour factoriser P sur C il ne reste que factoriser le polynôme x 2 + x + 1 sur C. Avec les notations usuelles du
p
cours, on
p 2
a ∆ = −3 et une racine carrée de ∆ est δ = i 3 d’où les deux racines complexes conjuguées x 1 = − 2 + i 2 = e i 3 π et
1 3
p 4
x 2 = − 12 − i 2
3
= e i 3 π . On conclut que la factorisation de P sur C est
2 4
P (x) = (x + 1)3 (x − 1)(x − e i 3 π )(x − e i 3 π ).
Exercice 2.17
On considère le polynôme
P (x) = x 5 − x 3 − x 2 + 1.
1. Écrire la formule de Taylor de P en 1 ; montrer que 1 est une racine de P et déterminer son ordre de multiplicité.
2. Factoriser P sur R en sachant qu’il possède une autre racine réelle évidente (i.e. elle est un diviseur du terme de degré
0).
3. Factoriser P sur C.
Correction
1. Pour un polynôme P de degré 5, la formule de Taylor en 1 s’écrit
k=5 (x − 1)k
P (k) (1)
X
P (x) = .
k=0 k!
P (x) = x 5 − x 3 − x 2 + 1 P (1) = 0
0 4 2
P (x) = 5x − 3x − 2x P 0 (1) = 0
P 00 (x) = 20x 3 − 6x − 2 P 00 (1) = 12
P 000 (x) = 60x 2 − 6 P 000 (1) = 54
P I V (x) = 120x P I V (1) = 120
P V (x) = 120 P V (1) = 120
© 2019-2020 G. Faccanoni 21
2 Fonctions polynomiales Mercredi 4 septembre 2019
= (x − 1)2 (x 3 + 2x 2 + 2x + 1)
Exercice 2.18
On considère le polynôme
P (x) = x 5 + 2x 4 + x 3 − x 2 − 2x − 1.
1. Écrire la formule de Taylor de P en −1 ; montrer que −1 est une racine de P et déterminer son ordre de multiplicité.
2. Factoriser P sur R en sachant qu’il possède une autre racine réelle évidente (i.e. elle est un diviseur du terme de degré
0).
3. Factoriser P sur C.
Correction
1. Pour un polynôme P de degré 5, la formule de Taylor en −1 s’écrit
k=5 (x + 1)k
P (k) (−1)
X
P (x) = .
k=0 k!
P (x) = x 5 + 2x 4 + x 3 − x 2 − 2x − 1 P (−1) = 0
0 4 3 2
P (x) = 5x + 8x + 3x − 2x − 2 P 0 (−1) = 0
P 00 (x) = 20x 3 + 24x 2 + 6x − 2 P 00 (−1) = −4
P 000 (x) = 60x 2 + 48x + 6 P 000 (−1) = 18
P I V (x) = 120x + 48 P I V (−1) = −72
P V (x) = 120 P V (−1) = 120
= (x + 1)2 (x 3 − 1)
22 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 2 Fonctions polynomiales
Étant donné que x 2 + x + 1 ne se factorise pas sur R, on conclut que la factorisation de P sur R est
1. P (x) = x 3 + 2x 2 − x − 2,
Correction
1. Une racine évidente de P est x = 1. En effectuant la division euclidienne de P par (x − 1) on obtient
x3 +2x 2 −x −2 x −1
−x 3 +x 2 x2 +3x +2
3x 2 −x −2
−3x 2 +3x
2x −2
−2x +2
0
1
2x 3 −(5 + 6i )x 2 +9i x +1 − 3i x 2
−2x 3 +6x 2 x2 +x +(1 − i )
(1 − 6i )x 2 +9i x +1 − 3i
−(1 − 6i )x 2 +3(1 − 6i )x
(3 − 9i )x +1 − 3i
−(3 − 9i )x −(3 − 9i )
0
3. Soit α ∈ R une racine réelle de P , alors elle satisfait l’équation 2α3 − (5 + 6i )α2 + 9i α + 1 − 3i = 0. En comparant les
parties réelle et imaginaire dans cette égalité on trouve
(
2α3 − 5α2 + 1 = 0
−6α2 + 9α − 3 = 0.
La résolution de −6α2 + 9α − 3 = 0 donne deux solutions : α = 1 et α = 1/2. Le réel α = 1 ne vérifiant pas l’équation
2α3 − 5α2 + 1 = 0 et le réel α = 1/2 oui, on a bien P (1/2) = 0. On factorise alors P en effectuant la division euclidienne
© 2019-2020 G. Faccanoni 23
2 Fonctions polynomiales Mercredi 4 septembre 2019
de P par (x − 1/2) :
On obtient P (x) = (x − 2)(x 2 + 3i x − 2). On voit que x = −i est une racine du polynôme x 2 + 3i x − 2 et on trouve
P (x) = (x − 2)(x + i )(x + 2i ).
5. Soit α ∈ R une racine réelle de P , alors elle satisfait l’équation 2α4 + 21 α3 + i α + 12 i = 0. En comparant les parties réelle
et imaginaire dans cette égalité on trouve α = −1/2 et on a bien P (−1/2) = 0. On factorise alors P en effectuant la
division euclidienne de P par (x + 1/2) :
x4 + 12 x 3 +i x + 21 i x +1/2
−x 4 − 12 x 3 x3 +i
+i x + 12 i
−i x − 12 i
0
P (x) = x 5 + 3x 4 + 4x 3 + 4x 2 + 3x + 1 P (−1) = 0
0 4 3 2
P (x) = 5x + 12x + 12x + 8x + 3 P 0 (−1) = 0
P 00 (x) = 20x 3 + 36x 2 + 24x + 8 P 00 (−1) = 0
P 000 (x) = 60x 2 + 72x + 24 P 000 (−1) = 12
On sait alors que P (x) = (x + 1)3Q(x) et on calcule Q en effectuant la division euclidienne de P par (x + 1)3 :
24 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 2 Fonctions polynomiales
7. Comme P (1) = 0 on sait que P est divisible par (x − 1). Avant d’effectuer la division euclidienne, cherchons d’abord la
multiplicité de cette racine :
P (x) = 2x 4 + x 3 − 6x 2 + x + 2 P (1) = 0
0 3 2
P (x) = 8x + 3x − 12x + 1 P 0 (1) = 0
P 00 (x) = 24x 2 + 6x − 12 P 00 (1) = 30
On sait alors que P (x) = (x − 1)2Q(x) et on calcule Q en effectuant la division euclidienne de P par (x − 1)2 :
2x 4 +x 3 −6x 2 +x +2 x 2 − 2x + 1
−2x 4 +4x 3 −2x 2 2x 2 + 5x + 2
5x 3 −8x 2 +x +2
−5x 3 +10x 2 −5x
2x 2 +4x −2
−2x 2 −4x +2
0
Correction
xk
¬ P n (x) = nk=0 (rappel : x 0 = 1 et 0! = 1)
P
k!
α = 0 n’est pas racine de P n (x) car P n (0) = 1
k−1 k= j +1 Pn−1 j
k x k! x
Pn
® P n0 (x) = k=0
= j =0 ( j + 1) ( j +1) j ! = P n−1 (x)
xk xk n n
+ xn! = P n0 (x) + xn!
Pn Pn−1
¯ P n (x) = k=0 k!
= k=0 k!
αn
Si α est racine de P n alors α 6= 0 et P n0 (α) = − n! 6= 0 donc α est une racine simple de P n .
Exercice 2.21
Soient a et b deux nombres complexes et P (x) := x 4 − 2i x 3 + 3(1 + i )x 2 + ax + b.
? Calculer a et b pour que i soit une racine multiple de P .
Dans toute la suite, on suppose que i est une racine multiple de P .
? Quelle est la multiplicité de la racine i de P ?
? Calculer les autres racines de P .
Correction
? Pour que i soit une racine multiple de P il faut que P (i ) = 0 et P 0 (i ) = 0 :
P (x) := x 4 − 2i x 3 + 3(1 + i )x 2 + ax + b P (i ) = −4 + 3i + ai + b
donc P (i ) = 0 si et seulement si b = 4 + 3i − ai ;
P 0 (x) := 4x 3 − 6i x 2 + 6(1 + i )x + a P 0 (i ) = 8i − 6 + a
© 2019-2020 G. Faccanoni 25
2 Fonctions polynomiales Mercredi 4 septembre 2019
? On doit alors calculer deux autres racines de P : soit on effectue la division euclidienne de P par (x − i )2 , soit on
factorise P à l’aide de la formule de Taylor, ce qui donne
(x − i )2 00 (x − i )3 000 (x − i )4 I V
P (x) = P (i ) + (x − i )P 0 (i ) + P (i ) + P (i ) + P (i )
2! 3! 4!
2 3 4
(x − i ) (x − i ) (x − i )
= 6(1 + i ) + 12i + 24
2! 3! µ 4! ¶µ ¶
1 − 3i 1 + 3i
= (x − i )2 (x 2 + 4 + 3i ) = (x − i )2 x − p x− p .
2 2
Exercice 2.22
Montrer que le polynôme P (x) := 2x 4 − 7x 3 + 9x 2 − 5x + 1 admet une racine triple que l’on déterminera, puis factoriser P à
l’aide de la formule de Taylor.
Correction
Comme P (1) = 0 on sait que P est divisible par (x − 1). Cherchons la multiplicité de cette racine :
P (x) = 2x 4 − 7x 3 + 9x 2 − 5x + 1 P (1) = 0
0 3 2
P (x) = 8x − 21x + 18x − 5 P 0 (1) = 0
P 00 (x) = 24x 2 − 42x + 18 P 00 (1) = 0
P 000 (x) = 48x − 42 P 000 (1) = 6
P I V (x) = 48 P I V (1) = 48
Exercice 2.23
Soit x 1 , x 2 , x 3 les trois racines du polynôme P (x) := x 3 − 3x 2 + 2x + 1. Sans calculer les racines de P , donner les valeurs
? σ1 := x 1 + x 2 + x 3 ,
? σ2 := x 1 x 2 + x 2 x 3 + x 3 x 1 ,
? σ3 := x 1 x 2 x 3 ,
? A := x 12 + x 22 + x 32 .
Correction
Soit Q(x) = ni=0 a i x i un polynôme avec a n 6= 0 et x 1 , . . . , x n ses racines. On peut écrire
P
à !
n n a n−k
n i n−i
(−1) σi x avec (−1)k = σk =
Y X X
Q(x) = a n (x − x i ) = a n x + xi 1 xi 2 . . . xi k .
i =1 i =1 an 1≤i 1 <i 2 <···<i k ≤n
En particulier
a n−1 n
= σ1 =
X
− xi = x1 + x2 + · · · + xn
an i =1
a n−2 n
= σ2 =
X
x i = x 1 x 2 + x 1 x 3 + · · · + x n−1 x n
an i<j
a0 n
(−1)n = σn =
X
xi = x1 x2 . . . xn
an i =1
Les quantités σi sont appelées fonctions symétriques élémentaires des racines du polynôme Q.
Pour n = 2, si Q(x) = a 0 + a 1 x + a 2 x 2 = a 2 (x − x 1 )(x − x 2 ), il y a deux fonctions symétriques élémentaires :
a1
σ1 = x 1 + x 2 = ,
a2
26 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 2 Fonctions polynomiales
a0
σ2 = x 1 x 2 = − .
a2
Correction
1 +x 2 +x 4 1 + x + x3
3
−1 −x −x 1 − x + 2x 2 − 3x 3 + 5x 4
−x +x 2 −x 3 +x 4
x +x 2 +x 4
2x 2 −x 3 +2x 4
−2x 2 −2x 3 −2x 5
−3x 3 +2x 4
−2x 5
3x 3 +3x 4 +3x 6
5x 4 −2x 5
+3x 6
−5x 4 −5x 5 −5x 7
−7x 5 +3x 6
−5x 7
d’où la relation
1 + x 2 + x 4 = (1 + x + x 3 )(1 − x + 2x 2 − 3x 3 + 5x 4 ) + x 5 (−7 + 3x − 5x 2 ).
Exercice 2.25
Factoriser P (x) := x 8 + x 4 + 1 sur R et sur C.
On remarquera pour commencer que P (x) = (x 4 + 1)2 − x 4 .
Correction
P (x) = x 8 + x 4 + 1
= (x 4 + 1)2 − x 4
= (x 4 + 1) − x 2 (x 4 + 1) + x 2
£ ¤£ ¤
= x4 − x2 + 1 x4 + x2 + 1
£ ¤£ ¤
= (x 2 − 1)2 − x 2 (x 2 + 1)2 − x 2
£ ¤£ ¤
= (x 2 − 1) − x (x 2 − 1) + x (x 2 + 1) − x (x 2 + 1) + x
£ ¤£ ¤£ ¤£ ¤
© 2019-2020 G. Faccanoni 27
Primitives et Intégrales
3
3.1 Primitives
3.1 Définition (Primitive)
Soit I un intervalle de R. Une fonction f : I → R est intégrable s’il existe une fonction dérivable F : I → R telle que pour tout
x ∈ I , F 0 (x) = f (x). Une telle fonction F est une primitive (ou intégrale indéfinie) de f .
3.3 Propriété
? Si F est une primitive de f alors, pour tout réel c, la fonction F + c est aussi une primitive de f .
? Toute primitive de f est nécessairement de la forme F + c pour une certaine constante c.
R R
Notation L’ensemble des primitives d’une fonction f est noté f ou encore f (x) dx.
Remarque Rx
Si a ∈ I alors F (x) = a f (t ) dt est l’unique primitive qui s’annule en a.
29
3 Primitives et Intégrales Mercredi 4 septembre 2019
Z Z
sin(x) dx = − cos(x) + c =⇒ sin(u(x))u 0 (x) dx = − cos(u(x)) + c
Z Z
cos(x) dx = sin(x) + c =⇒ cos(u(x))u 0 (x) dx = sin(u(x)) + c
Z Z
1 u 0 (x)
cos2 (x)
dx = tan(x) + c =⇒ cos2 (u(x))
dx = tan(u(x)) + c
Z Z
1 1 u 0 (x) 1
sin2 (x)
dx = − tan(x) +c =⇒ sin2 (u(x))
dx = − tan(u(x)) +c
Z Z
p 1 u 0 (x)
dx = arcsin(x) + c = − arccos(x) + c =⇒ p dx = arcsin(u(x)) + c = − arccos(u(x) + c
1−x 2 1−(u(x))2
Z Z
1 u 0 (x)
1+x 2
dx = arctan(x) + c =⇒ 1+(u(x))2
dx = arctan(u(x)) + c
Z Z
cosh(x) dx = sinh(x) + c =⇒ cosh(u(x))u 0 (x) dx = sinh(u(x)) + c
Z Z
sinh(x) dx = cosh(x) + c =⇒ sinh(u(x))u 0 (x) dx = cosh(u(x)) + c
Z Z
1 u 0 (x)
dx = tanh(x) + c =⇒ dx = tanh(u(x)) + c
cosh2 (x) cosh2 (u(x))
Remarque mnémotechnique : pour passer de la colonne de gauche à celle de droite il suffit de remplacer x par u(x) et dx
par u 0 (x) dx.
du
Autrement dit, en posant u = g (x) on obtient dx = g 0 (x), soit encore du = g 0 (x) dx et donc
Z Z
f (g (x))g 0 (x) dx = f (u) du = F (u) + c.
E XEMPLE
ln(x)
Calculer une primitive de x avec les trois méthodes décrites ci-dessus.
Intégration directe :
ln(x) [ f (x)]2 ln(x)
f (x) f 0 (x) dx avec f (x) = ln(x) et comme f (x) f 0 (x) dx =
R R R R
comme x dx = 2 on conclut que x dx =
2
ln (x)
2 + c.
Intégration par changement de variable :
dx
R ln(x) R u2 ln2 (x)
on pose u = ln(x) donc du = x et x dx = u du = 2 +c = 2 + c.
Intégration par parties :
ln(x) ln(x) ln(x)
si on pose g (x) = ln(x) et f 0 (x) = x1 alors g 0 (x) = 1
dx = ln2 (x) −
R R R
x et f (x) = ln(x) donc x x dx, i.e. 2 x dx =
ln2 (x)
ln2 (x) + c et finalement ln(x)
R
x dx = 2 + k.
30 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
3.9 Propriété
Soit N et D deux polynômes de degré respectivement ν et δ et P (x) = N (x)
D(x) une fraction rationnelle impropre (i.e. ν ≥ δ).
Alors, en effectuant la division euclidienne de N par D, on peut réécrire P comme
R(x)
P (x) = Q(x) +
D(x)
R(x)
où Q est un polynôme de degré ν − δ et R un polynôme de degré au plus δ − 1, ainsi D(x) est une fraction rationnelle propre.
On en déduit que
N (x) R(x)
Z Z Z Z
P (x) dx = dx = Q(x) dx + dx.
D(x) D(x)
L’intégration de Q étant triviale, on conclut que la difficulté de l’intégration d’une fraction rationnelle se réduit à l’intégration
d’une fraction rationnelle propre.
3.10 Propriété
R(x)
Si est une fraction rationnelle propre et si D possède
D(x)
? k racines réelles a k chacune de multiplicité m k et
? h couples de racines complexes conjuguées qui sont racines du polynôme x 2 + b h x + d h chacune de multiplicité n h
(ainsi ∆ = b h2 − 4d h < 0 pour tout h),
alors D s’écrit
Pour intégrer une fraction rationnelle il suffit alors de connaître les primitives des quatre fractions simples suivantes :
A A B x +C B x +C
f 1 (x) = , f 2 (x) = , f 3 (x) = , f 4 (x) = .
x −a (x − a)n x 2 + bx + d (x 2 + bx + d )n
© 2019-2020 G. Faccanoni 31
3 Primitives et Intégrales Mercredi 4 septembre 2019
A
2. la primitive de f 2 (x) = (x−a)n est
A A
Z
n
dx = + cnst;
(x − a) (1 − n)(x − a)n−1
B x+C
3. la primitive de f 3 (x) = x 2 +bx+d
est
B x +C B x +C
Z Z
dx = ´2 ³ ´ dx
x 2 + bx + d
³ 2
x + b2 + d − b4
q
2
x + b2 = d − b4 t
µq ¶ q
2
B d − b2
t − b dx = d − b4 dt
1
Z 4 2 +C
=q dt
2 t2 +1
d − b2
B 2t C − B b2 1
Z Z
= 2
dt + dt
1+ t2
q
2 1+t 2
d − b4
B C − B b2
= ln |1 + t 2 | + q arctan(t ) + cnst
2 2
d − b4
¯ ³ ´2 ¯¯
b b b
x
¯
B ¯ ¯ + 2
¯ C − B x +
2
= ln ¯1 + ¯+ q arctan q 2 + cnst;
¯
2
2 ¯¯ d − b4 ¯¯ d− 4b2
d− 4 b2
B x+C
4. la primitive de f 4 (x) = (x 2 +bx+d )n
est
B x +C B 2x + b Bb
µ ¶Z
1
Z Z
dx = dx + C − dx
(x 2 + bx + d )n 2 (x 2 + bx + d )n 2 (x 2 + bx + d )n
| {z } | {z }
I1 I2
avec
B 1
I1 =
2 (1 − n)(x + bx + d )n−1
2
1
Z
I2 = dx q
2
(x + bx + d )n
2
x + b2 = d − b4 t
¶ 1 −n Z q
2
b2 2 dx = d − b4 dt
µ
1
= d− dt
4 (1 + t 2 )n
| {z }
In
t 2n − 3
In = + I n−1 .
2(n − 1)(1 + t 2 )n−1 2(n − 1)
E XEMPLE
x−3
1. On veut intégrer la fraction rationnelle propre f (x) = x 2 −4x+5 . Comme ∆ = 42 − 4 × 5 < 0 il s’agit d’une intégrale du
R B x+C
type x 2 +bx+d dx. Le dénominateur se décompose comme x 2 − 4x + 5 = (x − 2)2 + 1 et l’intégrale s’écrit
x −3 (t + 2) − 3 1 2t 1
Z Z Z Z Z
f (x) dx = dx = dt = dt − dt
(x − 2)2 + 1 2
t +1 2
2 t +1 2
t +1
1 1
= ln(1 + t 2 ) − arctan(t ) + cnst = ln((x + 2)2 + 1) − arctan(x − 2) + cnst
2 2
3x + 1 A1 A2 A3
f (x) = = + +
x(x − 2)(x + 2) x x −2 x +2
32 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
Pour calculer les constantes A i on peut utiliser le principe d’identité des polynômes :
3x + 1 A 1 (x − 2)(x + 2) + A 2 x(x + 2) + A 3 x(x − 2) −4A 1 = 1
= ⇐⇒ 7 = 8A 2
x(x − 2)(x + 2) x(x − 2)(x + 2)
−5 = 8A 3
ainsi
1 1 7 A2 5 A3 1 7 5
Z Z Z Z
f (x) dx = − dx + dx − dx = − ln |x| + ln |x − 2| − ln |x + 2| + c.
4 x 8 x −2 8 x +2 4 8 8
3x 3 +2x−5
3. On veut intégrer la fraction rationnelle impropre f (x) = 3x 2 −5x−2
. On effectue d’abord la division euclidienne
3x 3 +2x −5 3x 2 − 5x − 2
−3x 3 +5x 2 +2x x + 35
5x 2 +4x −5
−5x 2 + 25
3 x + 10
3
37
3 x − 35
37
x− 5 37
x− 5
ainsi f (x) = x + 53 + 3x 32 −5x−2
3
. Maintenant on décompose le terme 3x 32 −5x−2
3
en fractions simples : on a 3x 2 − 5x − 2 =
1
(x + 3 )(x − 2) et on doit chercher les deux constantes A 1 et A 2 telles que
37 5
3 x−3 A1 A2
= 1
+
3x 2 − 5x − 2 x+ 3
x −2
On conclut que
x2 5
¯ ¯
5 52/21 69/7 52 ¯¯ 1 ¯¯ 69
Z Z
f (x) dx = x+ + + dx = + x + ln x + + ln |x − 2| + c.
3 x + 31 x −2 2 3 21 ¯ 3¯ 7
x −4 2
A 1,1 (x + 1) + A 2,1 (x − 3)(x + 1) + A 2,2 (x − 3) A 1,1 = −1/16
= ⇐⇒ A 2,1 = 1/16
x 3 − x 2 − 5x − 3 x 3 − x 2 − 5x − 3
A 2,2 = 5/4
On conclut que
−1/16 1/16 5/4
f (x) = + +
x −3 x + 1 (x + 1)2
et
1 1 5
Z
f (x) dx = − ln |x − 3| + ln |x + 1| + + c.
16 16 4(x + 1)
x +2 2
5. On veut intégrer la fraction rationnelle propre f (x) = (x 2 −2x+5) 2 . Comme ∆ = 4 − 20 < 0, la fonction f se décompose
comme
B 1 x +C 1 B 2 x +C 2
f (x) = 2 + 2 .
x − 2x + 5 (x − 2x + 5)2
© 2019-2020 G. Faccanoni 33
3 Primitives et Intégrales Mercredi 4 septembre 2019
x2 + 2 B 1 x +C 1 B 2 x +C 2
= +
(x 2 − 2x + 5)2 x 2 − 2x + 5 (x 2 − 2x + 5)2
B1 = 0
B 1 x 3 + (C 1 − 2B 1 )x 2 + (5B 1 − 2C 1 + B 2 )x + 5C 1 +C 2
C − 2B = 1
1 1
= ⇐⇒
(x 2 − 2x + 5)2 5B 1 − 2C 1 + B 2 = 0
5C 1 +C 2 = 2
? On calcule I 1 :
1 1
Z Z
dx = dx
x 2 − 2x + 5 (x − 1)2 + 4 x − 1 = 2t
1 1 dx = 2 dt
Z
= dt
2 t2 +1
x −1
µ ¶
1 1 1 1
Z
= dt = arctan(t ) + c = arctan + c;
2 1+ t2 2 2 2
? On calcule I 2 :
2x − 3 2x − 2 1 1 1
Z Z Z Z
dx = dx − dx = − − dx
(x 2 − 2x + 5)2 (x 2 − 2x + 5)2 (x 2 − 2x + 5)2 x 2 − 2x + 5 (x 2 − 2x + 5)2
1 1 1 1 1
Z Z
=− 2 − 2 2
dx = − 2 − dt
x − 2x + 5 ((x − 1) + 4) x − 2x + 5 8 (t + 1)2 2
1 1 t t
µ ¶ µ ¶
1 1 1 1 1
Z
=− 2 − + dt = − − + arctan(t )
x − 2x + 5 8 2 1 + t 2 2 1 + t 2 x 2 − 2x + 5 16 1 + t 2
x −1 x +7 x −1
µ µ ¶¶ µ ¶
1 1 2(x − 1) 1
=− 2 − + arctan + c = − − arctan + c.
x − 2x + 5 16 x 2 − 2x + 5 2 8(x 2 − 2x + 5) 16 2
On conclut que
x −1 x +7
µ ¶
7
Z
f (x) dx = arctan − + c.
16 2 8(x 2 − 2x + 5)
1
6. On veut intégrer la fraction rationnelle propre f (x) = x 3 (x 2 +1)2
qui se décompose comme
A 1 A 1 A 1 B 1 x +C 1 B 2 x +C 2
f (x) = + 2+ 3+ + 2 .
x x x x2 + 1 (x + 1)2
1 1 2x −2x 1 2x
Z Z Z Z Z Z
2
f (x) dx = −2 dx + 3
dx + 2
dx + 2 2
dx = −2 ln |x| − 2
+ ln |x + 1| − dx
x x x +1 (x + 1) 2x (x + 1)2
2
1 2x
= −2 ln |x| − 2 + ln |x 2 + 1| − 2 + c.
2x (x + 1)
34 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
Remarque Rx
Si f est de classe C 1 sur [a; b] alors g (x) = f (x) − f (a) = a f 0 (t ) dt .
y
f (x) = x 2
3 y
4
9 f (x) = (x − 1)2
4
2
1
3
x+
=
x)
g(
2Z 2
1
¸3/2 (g (x) − f (x))dx
3/2 x3 0
·
9
Z
2
x dx = = 1
0 3 0 8
x
E XEMPLE 1 3 2 3 x
2 0 1 2 3
3.15 Propriété
Soit [a; b] un intervalle de R avec a < b, f et g : [a; b] → R deux fonctions intégrables.
Relation de Chasles : ab f (x) dx = ac f (x) dx + cb f (x) dx pour tout c ∈ [a; b].
R R R
Relation d’ordre : si f (x) ≤ g (x) pour toute x ∈ [a; b] alors ab f (x) dx ≤ ab g (x) dx.
R R
Nullité : si f est continue sur [a; b] et f (x) ≥ 0 pour toute x ∈ [a; b] alors ab f (x) dx = 0 si et seulement si f (x) ≡ 0.
R
Ra Ra
ParitéR : si f est paire sur [−a; a] avec a ≥ 0 alors −a f (x) dx = 2 0 f (x) dx ; si f est impaire sur [−a; a] avec a ≥ 0 alors
a
−a f (x) dx = 0.
© 2019-2020 G. Faccanoni 35
3 Primitives et Intégrales Mercredi 4 septembre 2019
Primitives
Exercice 3.1 (Par intégration directe)
Calculer les primitives suivantes :
p p 1
Z Z Z
1. 2x 3 − 3x + 1 dx 2. x+ 3
x dx 3. p dx
x +1
Z p
1
Z p Z
4
4. 2x + 1 dx 5. (x − 1)3 dx 6. p dx
x3x
x x3 + x + 1
Z Z Z
7. (1 + 2x 3 )2 dx 8. dx 9. dx
x +1 x2 + 1
Correction p
1 3 1 p
3 4/3
1. 2 x(x − 3x + 2) + c 2. 12 (8 x + 3)x +c 3. 2 x + 1 + c
p
(2x+1)3 4 7/4 3
4. 3 +c 5. 7 (x − 1) +c 6. − p
3 +c
x
4 7 4 x+1−1 1
R R
7. 7x +x +x +c 8. x+1 dx = 1 − x+1 dx = x − ln(x + 1) + c
R x(x 2 +1)+1 R R 1 x2
9. x 2 +1
dx = x dx + x 2 +1
dx = 2 + arctan(x) + c
Exercice 3.2
Calculer les primitives suivantes :
1 1 1 1
Z Z Z Z
1. dx 2. dx 3. dx 4. dx
1+x 1−x 1 + x2 1 − x2
Correction
u 0 (x)
1. u(x) = 1 + x, u 0 (x) = 1,
R
u(x) dx = ln |1 + x| + c
u 0 (x)
2. u(x) = 1 − x, u 0 (x) = −1, −
R
u(x) dx = − ln |1 − x| + c
3. arctan(x) + c
1 1/2 1/2 1 1 1
dx = 21 (ln |1 + x| − ln |1 − x|) + c = 12 ln
R ¡ 1+x ¢
4. 1−x 2
= 1−x + 1+x , 2 1−x + 1+x 1−x +c
Correction
u 0 (x)
1. u(x) = 1 + e x , u 0 (x) = e x , dx = ln(1 + e x ) + c
R
u(x)
36 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
1+e x −e x u 0 (x)
2. u(x) = 1 + e x , u 0 (x) = e x , dx = x − ln(1 + e x ) + c
R R R
1+e x dx = 1 dx − u(x)
ln4 (x)
3. u(x) = ln(x), u 0 (x) = 1/x, (u(x))3 u 0 (x) dx =
R
+c
4
0 2
R u(x) 0 −1/x
4. u(x) = −1/x, u (x) = 1/x , e u (x) dx = e +c
u 0 (x)
5. u(x) = ln(x), u 0 (x) = 1/x, (u(x))3 dx = − 12 + c
R
2 ln (x)
0
R u 0 (x)
6. u(x) = x + sin(x), u (x) = 1 + cos(x), u(x) dx = ln |x + sin(x)| + c
R u 0 (x)
7. u(x) = x 2 , u 0 (x) = 2x, 1+(u(x)) 2
2 dx = arctan(x ) + c
R 0 (x)
8. u(x) = 1 + sin2 (x), u 0 (x) = 2 sin(x) cos(x) = sin(2x), uu(x) dx = ln(1 + sin2 (x)) + c
4
9. u(x) = sin(x), u 0 (x) = cos(x), (u(x))3 u 0 (x) dx = sin4(x) + c
R
R 0 (x)
10. u(x) = tan(x), u 0 (x) = 1/ cos2 (x), sin(x) cos(x)
dx = uu(x)
R
cos2 (x)
dx = ln | tan(x)| + c
11. u(x) = tan(x), u 0 (x) = 1/ cos2 (x), e u(x) u 0 (x) dx = e tan(x) + c
R
p
12. u(x) = 1 − x 2 , u 0 (x) = −2x, p 1 2 dx + (u(x))−1/2 u 0 (x) dx = arcsin(x) + 2 1 − x 2 + c
R R
1−x
R 0 (x)
13. u(x) = x + 2x + 1, u (x) = 2x + 2, 12 uu(x)
2 0
dx = 21 ln(x 2 + 2x + 1) + c
2 3
14. u(x) = x 2 + 1, u 0 (x) = 2x, 21 (u(x))2 u 0 (x) dx = (x +1)
R
6 +c
2x+1
15. u(x) = 2x + 1, u 0 (x) = 2, 21 e u(x) u 0 (x) dx = 21 e u(x) + c = e 2 + c
R
x2
18. u(x) = x 2 , u 0 (x) = 2x, 1
e u(x) u 0 (x) dx = e 2 + c
R
2
x3
19. u(x) = x 3 , u 0 (x) = 3x 2 , 31 e u(x) u 0 (x) dx = e 3 + c
R
p p
20. u(x) = x, u 0 (x) = 2p1 x , 12 sin (u(x)) u 0 (x) dx = −2 cos x + c
R
1 3p3
21. u(x) = x 2 + 3, u 0 (x) = 2x, (u(x))−1/3 u 0 (x) dx =
R
2 (x 2 + 3)2 + c
4
R u 0 (x) 4)
22. u(x) = 1 + x 4 , u 0 (x) = 4x 3 , 41 dx = ln(1+x
4 +c
u(x)
23. u(x) = 3x, u 0 (x) = 3, 13 sin(u(x))u 0 (x) dx = − 31 cos(3x) + c
R
n+1
24. u(x) = ax + b, u 0 (x) = a, a1 (u(x))n u 0 (x) dx = a1 (ax+b)
R
n+1 +c
−n+1
25. u(x) = ax + b, u 0 (x) = a, a1 (u(x))−n u 0 (x) dx = a1 (ax+b)
R
−n+1 +c
Exercice 3.4
Calculer les primitives suivantes :
cos2 (x) sin(x) + cos(x) 1
Z Z Z
1. dx 2. dx 3. dx
1 − sin(x) sin(x) − cos(x) sin2 (x) cos2 (x)
1 1 1
Z Z Z
4. dx 5. dx 6. dx
cos(x) sin(x) sin(x) cos(x)
Correction
(1−sin(x))(1+sin(x))
1. cos 2 (x) = 1 − sin2 (x) = (1 − sin(x))(1 + sin(x)),
R
dx = x − cos(x) + c
1−sin(x)
R 0 (x)
dx = uu(x)
R sin(x)+cos(x)
2. 1 = cos 2 (x)+sin2 (x), u(x) = sin(x)−cos(x), u 0 (x) = sin(x)+cos(x), sin(x)−cos(x) dx = ln(sin(x)−cos(x))+c
2 2 1
1 = cos 2 (x) + sin2 (x), sin (x)+cos (x)
dx = cos12 (x) dx + sin12 (x) dx = tan(x) −
R R R
3. sin2 (x) cos2 (x)
+c
tan(x)
2 (x)+sin2 (x)
1 = cos 2 (x) + sin2 (x), cos(x)1sin(x) dx = cos
R R R sin(x) cos(x)
4. cos(x) sin(x) dx = cos(x) + sin(x) dx = − ln | cos(x)| + ln | sin(x)| + c
R 1
sin(x) = 2 sin(x/2) cos(x/2), u(x) = x/2, u 0 (x) = 1/2, sin(x) dx = cos(u(x))1sin(u(x)) u 0 (x) dx = − ln | cos(x/2)|+ln | sin(x/2)|+
R
5.
c = ln | tan(x/2)| + c.
cos(x) = sin π − x , u(x) = π − x, u 0 (x) = −1, ¡1 ¢ u 0 (x) dx = ln | cos( π − x )| − ln | sin( π − x )| + c =
¡ ¢ R 1 R
6. 2 2 dx = −
cos(x) 4 2 4 2
sin u(x)
ln(sin(x) − cos(x)) + c
© 2019-2020 G. Faccanoni 37
3 Primitives et Intégrales Mercredi 4 septembre 2019
x3 x5
Z p
1
Z Z Z
(9) dx (10) p dx (11) p dx (12) e x − 1 dx
x(2 + ln2 (x)) 1 − x2 x3 − 1
tan(x)
ln(x) 1 e x
Z Z Z Z
(13) dx (14) x
dx (15) dx (16) p dx
x e + e −x cos2 (x) 1 + x2
1 e 1/x cos(x)
Z p Z Z Z
(17) x a + x 2 dx (18) dx (19) dx (20) dx
x2
p
x 2
1 − ln (x) 1 + sin(x)
x3 x
µ ¶
1 1 1
Z Z Z Z
(21) 2
cos dx (22) p dx (23) dx (24) dx
x x 1 + x2 3 + x2 1 + x4
Correction
1. Pour x > 0, si on pose t = ln(x) alors x1 dx = dt et on obtient − cos(ln(x)) + c
p p p
2. Pour x > 0, si on pose t = x alors dx = 2t dt et on obtient 2( x + e x ) + c
p p
3. Pour x > 0, si on pose t = x alors dx = 2t dt et on obtient 2 ln( x − 1) + c
p p p
4. Pour x > 0, si x > 0. Si on pose t = x alors dx = 2t dt et on obtient 4 1 + x + c
5. Si on pose t = e x alors e x dx = dt et on obtient ln(1 + e x ) + c
6. Pour x > 0, si on pose t = ln(x) alors x1 dx = dt et on obtient ln | ln(x)| + c
q ¡ ¢
7. Pour x > 0, si on pose t = ln x1 alors − x1 dx = 2t dt et on obtient −2 ln x1 + c
¡ ¢
8. Si on pose t = 1 + e x alors e x dx = Rdt ainsi e Rx ln(1 + e x ) dx = ln(t ) dt . Sans utiliser l’intégration par partie, si on
R R
38 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
ln(x) ln(ln(x))
Z Z Z Z
1. dx 2. dx 3. ln(1 + x) dx 4. x 2 e x dx
x2 x
ln(x)
Z Z Z Z
5. p dx 6. x sin(x) dx 7. x ln(x) dx 8. x 2 cos(x) dx
x
sin(x) tan(x)
Z Z Z Z
9. e dx 10. x 3 sin(x 2 ) dx 11. e 2x sin(3x) dx 12. e −3x cos(2x) dx
cos3 (x)
ln(x)
Z Z Z Z
13. x 3 ln(x) dx 14. p4
dx 15. ln2 (x) dx 16. x sin2 (x) dx
x
Correction
1. f (x) = ln(x) ⇒ f 0 (x) = 1
x et g (x) = − x1 ⇐ g 0 (x) = 1
x2
on obtient − 1+ln(x)
x +c
1 1
2. f (x) = ln(ln(x)) ⇒ f 0 (x) = 0
x ln(x) et g (x) = ln(x) ⇐ g (x) = x on obtient (ln(ln(x)) − 1) ln(x) + c
1
3. f (x) = ln(1 + x) ⇒ f 0 (x) = 0
1+x et g (x) = x ⇐ g (x) = 1 on obtient (1 + x) ln(1 + x) − x + c
x x x
4. f (x) = x 2 ⇒ f 0 (x) = 2x et g (x) = e ⇐ g (x) = e on obtient e ((x − 2)x + 2) + c
0
p p
5. f (x) = ln(x) ⇒ f 0 (x) = x1 et g (x) = 2 x ⇐ g 0 (x) = p1x on obtient 2 x(ln(x) − 2) + c
6. f (x) = x ⇒ f 0 (x) = 1 et g (x) = − cos(x) ⇐ g 0 (x) = sin(x) on obtient −x cos(x) + sin(x) + c
7. f (x) = ln(x) ⇒ f 0 (x) = 1/x et g (x) = x 2 /2 ⇐ g 0 (x) = x on obtient 21 x 2 ln(x) − 14 x 2 + c
8. f (x) = x 2 ⇒ f 0 (x) = 2x et g (x) = sin(x) ⇐ g 0 (x) = cos(x) on obtient x 2 sin(x) − 2[−x cos(x) + sin(x)] + c
9. f (x) = sin(x)/cos(x) ⇒ f 0 (x) = 1/cos 2 x et g (x) = e tan(x) ⇐ g 0 (x) = e tan(x)/x 2 on obtient e tan(x) (tan(x) − 1) + c
−x 2 cos(x 2 ) + sin(x 2 )
10. f (x) = sin(x 2 ) ⇒ f 0 (x) = 2x cos(x 2 ) et g (x) = x 4/4 ⇐ g 0 (x) = x 3 on obtient +c
2
2e 2x (sin(3x)− 23 cos(3x))
11. f (x) = sin(3x) ⇒ f 0 (x) = 3 cos(3x) et g (x) = e 2x /2 ⇐ g 0 (x) = e 2x on obtient 13 +c
−e −3x 3e −3x (cos(2x)− 32 sin(2x))
12. f (x) = cos(2x) ⇒ f 0 (x) = −2 sin(2x) et g (x) = 3 ⇐ g 0 (x) = e −3x on obtient − 13 +c
1 x4 1 4
13. f (x) = ln(x) ⇒ f 0 (x) = x et g (x) = 4 ⇐ g 0 (x) = x 3 on obtient 16 x (4 ln(x) − 1) + c
1 4x 3/4 1 4 3/4 ¡
14. f (x) = ln(x) ⇒ f 0 (x) = ⇐ g 0 (x) = ln(x) − 43 + c
¢
x et g (x) = 3
p
4 x
on obtient
x 3
1
15. f (x) = ln(x) ⇒ f 0 (x) = x et g (x) = x ln(x) − x ⇐ g 0 (x) = ln(x) on obtient x(ln2 (x) − 2 ln(x) + 2) + c
x2
16. f (x) = sin2 (x) ⇒ f 0 (x) = 2 sin(x) cos(x) = sin(2x) et g (x) = 2 ⇐ g 0 (x) = x on obtient 41 (x 2 − x sin(2x) − 21 cos(2x)) + c
Exercice 3.7
Calculer la primitive suivante en utilisant un changement de variable. Comparer ensuite au résultat obtenu en utilisant
l’intégration par parties :
x arcsin(x)
Z
p dx
1 − x2
Correction
CV Si on pose t = arcsin(x) alors p 1 dx = dt et x = sin(t ) et on obtient
1−x 2
x arcsin(x)
Z Z
p dx = t sin(t ) dt
1 − x2
On a calculé cette intégrale à l’exercice 3.6(6.) :
Z
t sin(t ) dt = −t cos(t ) + sin(t ) + c = − arcsin(x) cos(arcsin(x)) + x + c
IPP
x arcsin(x)
Z p
p dx = = − 1 − x 2 · arcsin(x) + x + c
1 − x2 f (x) = arcsin(x) ⇒ f 0 (x) = p 1
p 1−x 2
g (x) = − 1 − x 2 ⇐ g 0 (x) = p x
1−x 2
p p
Les deux calculs donnent le même résultat car cos(arcsin(x)) = ± 1 − sin2 (arcsin(x)) = ± 1 − x 2
© 2019-2020 G. Faccanoni 39
3 Primitives et Intégrales Mercredi 4 septembre 2019
Correction
? Si f (x) = k 1 pour tout x ∈ R et g (x) = k 2 pour tout x ∈ R, alors ( f g )0 (x) = (k 1 k 2 )0 = 0 pour tout x ∈ R et f 0 (x)g 0 (x) =
0 × 0 = 0 pour tout x ∈ R.
? Si g = f , on cherche f telle que ( f 2 )0 = ( f 0 )2 , c’est-à-dire 2 f (x) f 0 (x) = ( f 0 (x))2 pour tout x ∈ R. Donc, soit f 0 (x) = 0
pour tout x ∈ R et on trouve à nouveau f (x) = g (x) = k pour tout x ∈ R, soit 2 f (x) = f 0 (x) pour tout x ∈ R et on trouve
f (x) = g (x) = ke 2x pour tout x ∈ R.
? Dire que ( f g )0 = f 0 g 0 revient à dire que f 0 g + f g 0 = f 0 g 0 . En divisant par le produit f g (il est inutile à ce stade de se
préoccuper de la possibilité de diviser par 0, nous cherchons seulement formellement des conditions nécessaires) on
a
f 0 (x) g 0 (x) f 0 (x) g 0 (x)
+ = ·
f (x) g (x) f (x) g (x)
g 0 (x)
f 0 (x) g (x)
c’est-à-dire f (x) = g 0 (x)
, soit encore
1− g (x)
g 0 (x)
[ln( f (x))]0 = .
g 0 (x) − g (x)
g 0 (x)
Si on choisit g , il suffit de poser f = e G où G est une primitive de g 0 (x)−g (x)
.
Voyons quelques exemples :
1 1
g )0 = f 0 g 0 :
R
? si on pose g (x) = x alors G(x) = 1−x dx = − ln(1 − x) et f (x) = 1−x . Vérifions si on a bien ( f
³ x ´0 1
( f g )0 (x) = =
1−x (1 − x)2
1
f 0 (x)g 0 (x) =
(1 − x)2
ax a−1 ax a−1 1
si on pose g (x) = x a alors G(x) =
R R
?
ax a−1 −x a
dx = ax a−1 −x a
dx = −a ln(a − x) et f (x) = (a−x)a . Vérifions si on a
bien ( f g )0 = f 0 g 0 :
¶0
xa
µ
( f g )0 (x) = = a 2 x a−1 (a − x)−a−1
(a − x)a
f 0 (x)g 0 (x) = a(a − x)−a−1 · ax a−1 = a 2 x a−1 (a − x)−a−1
ae ax a
si on pose g (x) = e ax alors G(x) = et f (x) = e bx où b = a/(a − 1). Vérifions si on a bien
R
? ae ax −e ax dx = a−1 x
( f g )0 = f 0 g 0 :
40 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
Correction
1. On pose I n = x n e αx dx. En intégrant par parties ( f (x) = x n et g 0 (x) = e αx ) on trouve
R
e αx n e αx n n−1 e αx n − 1
¶ αx
n n(n − 1) n−2 n n! e
µ ¶ µ
In = xn − I n−1 = x n − x − I n−2 = · · · = x n − x n−1 + x · · · + (−1) +c
α α α α α α α α2 αn α
2. On pose I n = sinn (x) dx. En intégrant par parties ( f (x) = sinn−1 (x) et g 0 (x) = sin(x)) on trouve
R
De la même manière, on pose I n = cosn (x) dx. En intégrant par parties ( f (x) = cosn−1 (x) et g 0 (x) = cos(x)) on trouve
R
3. On pose I n = x n sin(x) dx et J n = x n cos(x) dx. En intégrant par parties ( f (x) = x n et g 0 (x) = sin(x) dans la première
R R
Par conséquence
I n = −x n cos(x) + n x n−1 sin(x) − (n − 1)I n−2 = −x n cos(x) + nx n−1 sin(x) − n(n − 1)I n−2
¡ ¢
J n = x n sin(x) − n −x n−1 cos(x) + (n − 1)J n−2 = x n sin(x) + nx n−1 cos(x) − n(n − 1)J n−2
¡ ¢
4. On pose I n = x α lnn (x) dx. En intégrant par parties ( f (x) = lnn (x) et g 0 (x) = x α ) on trouve
R
x α+1 n
¶ α+1
n n n(n − 1) n−2 n! x
µ
n n−1 n
In = ln (x) − I n−1 = · · · = ln (x) − ln (x) + ln (x) · · · + (−1) + c.
α+1 α+1 α+1 (α + 1)2 (α + 1)n α + 1
a 2x − 1 3x 2 + 2x − 5
Z Z Z
a) dx c) dx e) dx
x −b (x − 1)(x − 2) 3x 2 − 5x − 2
a x2
Z Z
3x + 1
Z
b) dx, n 6= 1 d) dx f) dx
(x − b)n x 3 − 4x (x 2 − 1)2
Correction
Soit c une constante réelle.
R a
a) x−b dx = a ln |x − b| + c
R a a
b) (x−b)n dx = (1−n)(x−b) n−1 + c
2x−1 A B 2x−1 1 1
R R R
c) Comme (x−1)(x−2) = (x−1) + (x−2) ssi A + B = 2 et −2A − B = −1, alors (x−1)(x−2) dx = − (x−1) dx + 3 (x−2) dx = − ln |x −
¯ 3
¯
1| + 3 ln |x − 2| + c = ln ¯ (x−2)
x−1 ¯ + c
¯ ¯
3x+1 −1/4
+ −5/8 7/8 3x+1
dx = − ln4|x| − 5 ln |x+2| + 7 ln |x−2|
R
d) Comme x 3 −4x
= x x+2 + x−2 alors x 3 −4x 8 8 +c
© 2019-2020 G. Faccanoni 41
3 Primitives et Intégrales Mercredi 4 septembre 2019
¯ ¯
3x 2 +2x−5 (3x 2 −5x−2)+(7x−3) R 3x 2 +2x−5
= 1 + 3x 27x−3 16 ¯ 11
e) Comme = alors dx = x + ln + 1 ¯ + 7 ln |x − 2| + c
¯
3x 2 −5x−2 3x 2 −5x−2 −5x−2 2
3x −5x−2 21 ¯3x
2
³ ´
x2
f) Comme (x 2 −1)2 = (x−1)2 (x+1)2 et (x 2x−1)2 = 14 x−1
1 1 1 1
dx = 14 ln |x − 1| − x−1
1 1
R ¡ ¢
+ (x−1) 2 − x+1 + (x+1)2 alors (x 2 −1)2
− ln |x + 1| − x+1
Intégrales
Exercice 3.11
La valeur moyenne de la fonction f (x) = x 3 sur l’intervalle [0; k] est 9. Calculer k.
Correction
1 k 1 k4 p
Z
3
x 3 dx = 9 =⇒ =9 =⇒ k= 36.
k 0 k 4
Correction
? Vitesse maximale : v(t ) = V t (1 − t ), v 0 (t ) = 1 − 2t , v 0 (t ) = 0 ssi t = 12 et v 12 = V4 .
¡ ¢
Z 1
V
? Distance parcourue : v(t ) = x 0 (t ) donc x parcourue = v(t ) dt = .
0 6
Exercice 3.13
Calculer
p
Z
11 Z 1/ 3 1
Z 1 2x − 5
A= dx, B= p dx, C= dx.
−1 1 + x2 −1/ 3 1 + 3x 2 −1 x 2 − 5x + 6
Correction
π −π π
A = [arctan(x)]1−1 = arctan(1) − arctan(−1) = − = ,
4 4 2
A π
Z 1
1 1
B= p dt = p = p ,
3 −1 1 + t 2 3 2 3
Z 1
2x − 5 ¤1
C= dx = ln |x 2 − 5x + 6| −1 = ln(2) − ln(12) = ln(2) − ln(3) − 2 ln(2) = − ln(2) − ln(3).
£
2
−1 x − 5x + 6
x3
a) f (x) = −x 2 + x + 2 et g (x) = x 2 − 3x + 2 b) f (x) = 4 et g (x) = x 2 − x
Correction
a) Comme f (x) = g (x) ssi x ∈ {0, 2} et f (x) ≥ g (x) pour x ∈ [0, 2], l’aire comprise entre le graphe de la fonction f (x) et le
¸2
x3 x2
Z 2 Z 2 ·
2 8
graphe de la fonction g (x) est ( f (x) − g (x)) dx = −2x + 4x dx = −2 + 4 = .
0 0 3 2 0 3
b) Comme f (x) = g (x) ssi x ∈ {0, 2} et f (x) ≥ g (x) pour x ∈ [0, 2], l’aire comprise entre le graphe de la fonction f (x) et le
Z 2 3 · 4 ¸2
x3 x2
Z 2
x x 1
graphe de la fonction g (x) est ( f (x) − g (x)) dx = − x 2 + x dx = − + = .
0 0 4 16 3 2 0 3
42 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
x3
f (x) = −x 2 + x + 2 et g (x) = x 2 − 3x + 2 f (x) = 4 et g (x) = x 2 − x
y y
2 2
−1 −1
g(x)
g(x)
1 2 x f (x) 1 2 x
f (x)
Correction
Remarquons d’abord que
Z Z
cos2 (x)=1−sin2 (x)
cos3 (x) dx = (1 − sin2 (x)) cos(x) dx =
Z Z t =sin x
dt =cos(x) dx
= cos(x) − cos(x) sin2 (x) dx = sin(x) − cos x sin2 (x) dx =
t3 sin3 (x)
Z
= sin(x) − t 2 d t = sin(x) − + k = sin(x) − +k , k ∈ R.
3 3
1 π 5π
Comme cos3 (x) = 8 ssi x = 3 ou x = 3 alors
ÃZ π π
! Ãp p ! p
3
3 3 −0 3 3 π 9 3−π
Aire (A) = 2 cos (x) dx − =2 − − =
0 8 2 8 24 12
et Ãp p p
π π
!
−0 3 3 π 9 3−π
Z
3 3
3
Aire (B ) = cos (x) dx − 2 · =2 − − = .
− π3 8 2 8 24 12
Exercice 3.16 y
Dans la figure ci-contre on a tracé le graphe de la fonction
g : [0; w] → R définie par
Z x
g (x) = f (t ) dt y = g (x)
0 x
0 h k w
avec f : [0; w] → R une fonction continue et dérivable. Le
graphe de g a tangente horizontale en x = 0 et présente un
changement de concavité en x = h et un maximum en x = k.
© 2019-2020 G. Faccanoni 43
3 Primitives et Intégrales Mercredi 4 septembre 2019
Correction
1. g 0 (x) = f (x) pour tout x ∈ [0; w]. Puisque x = 0 et x = k sont des points à tangente horizontale pour le graphe de g ,
alors g 0 (0) = g 0 (k) = 0.
2. f est continue par hypothèse. D’après le théorème de W EIERSTRASS elle admet un maximum et un minimum sur
[0; w].
On a vu au point précédent que f (0) = f (k) = 0 et que g 0 (x) = f (x) pour tout x ∈ [0; w]. g est croissante (g 0 (x) > 0) sur
[0; k] et décroissante (g 0 (x) < 0) sur [k; w], donc f est positive sur [0; h] et négative sur [h; w].
De plus, g 00 (x) = f 0 (x) pour tout x ∈ [0; w]. g est convexe (g 00 (x) > 0) sur [0; h] et concave (g 00 (x) < 0) sur [h; w], donc f
est croissante sur [0; h] et décroissante sur [h; w]. x = h est un maximum absolu pour f et x = w un minimum absolu.
Un graphe plausible de f est donc le suivant :
y
w
0 h k x
y = f (x)
3. g est une fonction polynomiale de degré 3, x = w est un zéro simple et x = 0 est un zéro double (car g 0 (0) = 0), donc
g (x) = a(x − w)(x − 0)2 = ax 2 (x − w) avec a ∈ R∗ un paramètre.
3.1. On a f (x) = g 0 (x) = ax(3x − 2w) : il s’agit d’une parabole. Comme f (0) = f (k) = 0 alors k = 2w/3. De plus,
le sommet de la parabole se trouve en x = k/2 et f 0 (k/2) = 0. Comme x = h est le maximum de f , alors
h = k/2 = w/3.
3.2. Si w = 3 alors g (x) = ax 2 (x − 3) et la condition g (1) = 2/3 implique a = −1/3. On obtient ainsi g (x) = − 13 x 3 + x 2 .
Dans ce tableau, ce qu’on nomme le “taux d’imposition” est en fait le taux marginal : c’est un pourcentage qui ne s’applique
qu’à une partie des revenus, celle de la tranche concernée. On voit que le taux marginal est plus important pour les forts
revenus : c’est ce qu’on appelle la progressivité de l’impôt. On voit aussi que dans chaque tranche, le montant d’impôt à
payer est proportionnel au QF : on dit que la fonction impôts est linéaire par morceaux. Ces tranches d’imposition signifient
la chose suivante :
? Si le revenu est inférieur à 9710 ¤, on ne paye pas d’impôt.
? Si le revenu est compris entre 9710 ¤ et 26 818 ¤, on ne paye pas d’impôt sur les premiers 9710 ¤ de son revenu, et
on paye 14% de la partie qui excède 9710 ¤. Par exemple, si le revenu est 9711 ¤, on payera 14% de 1 ¤, c’est-à-dire
14 centimes.
1. Il s’agit du revenu du contribuable auquel on retire certaines sommes (par exemple au titre des frais professionnels).
44 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
? Si le revenu est compris entre 26 818 ¤ et 71 898 ¤, on ne paye rien sur la première tranche de 9710 ¤, puis 14% sur la
deuxième tranche, allant de 9710 ¤ à 26 818 ¤ (soit 14% de (26818 − 9710) ¤= 2395.12 ¤), et enfin 30% sur la partie
du revenu qui excède 71 898 ¤. Par exemple : si le revenu est de 72 000 ¤, l’impôt sera de 0 + 2395.12 + 30% × (72000 −
71898) = 2425.72 ¤.
? De même on trouve pour les impôts dans chacune des autres tranches : 0.41 × (152260 − 71898) = 32 948.42 ¤, etc.
Prenons l’exemple concret d’un contribuable, célibataire, qui en 2017 déclare 100 000 ¤ de revenu pour l’année 2016. On
veut calculer explicitement les impôts à payer. Ce contribuable est dans une tranche d’imposition marginale de 41% donc
au total a un montant d’impôts de 0 × (9710 − 0) + 14% × (26818 − 9710) + 30% × (71898 − 26818) + 41% × (100000 − 71898) =
27 440.94 ¤ et est ainsi redevable de 27.44% de ses revenus (et non 41%).
Écrire et tracer le graphe des fonctions suivantes en fonction du revenu imposable (ou du QF) :
1. Taux d’imposition marginal,
2. Montant de l’impôt,
3. Taux réel d’imposition,
4. Ce qui reste après avoir payé l’impôt.
En traçant la courbe des impôts payés en fonction du QF on verra tout de suite qu’il n’y a aucun «saut» dans la courbe
lorsqu’on change de tranche. Mathématiquement parlant, l’impôt est une fonction continue du QF pour qu’il n’y ait pas
d’injustices : une petite modification du revenu n’entraîne qu’une petite modification de l’impôt (sauter d’une tranche
n’est pas un drame). Avec des taux marginaux par tranches, on obtient une fonction croissante (plus le revenu est élevé et
plus l’impôt est élevé) et affine par morceaux ; et puisque le taux marginal augmente en fonction du revenu il s’agit d’une
fonction convexe. Cela signifie que la courbe «monte toujours plus vite» : plus le revenu est élevé et plus le taux marginal
augmente. Non seulement l’impôt augmente en fonction du revenu mais il augmente de plus en plus vite.
Correction
1. Commençons par tracer le graphique responsable de la crainte du saut d’une tranche : en horizontale, le revenu
imposable, en verticale, le taux de la tranche correspondante, qu’on appelle le taux marginal.
Taux marginal
45%
41%
30%
14%
0%
Revenu
¤
8¤
8¤
¤
10
60
81
89
97
22
26
71
15
Ce graphique présente effectivement des sauts importants. Mais ces sauts ne concernent pas l’impôt mais la dérivée
de la fonction Impôt I par rapport au revenu R. La fonction Impôt I : R+ → R+ est une fonction continue mais pas sa
dérivée I 0 : R+ → R+ qui est définie par
0 si R < 9710 ¤
0.14 si 9710 ¤ < R < 26 818 ¤
0
I (R) = 0.30 si 26 818 ¤ < R < 71 898 ¤
0.41 si 71 898 ¤ < R < 152 260 ¤
0.45 si R > 152 260 ¤
Ce graphique indique le pourcentage d’imposition qui serait appliqué sur chaque euro supplémentaire qui viendrait
s’ajouter au revenu : si votre revenu est par exemple de 80 000 ¤, alors pour chaque euro supplémentaire que vous
gagnez, 41 centimes seront pour les impôts.
2. Il s’agit de la fonction I : R+ → R+ définie par
Z R
I (R) = I 0 (r ) dr
0
© 2019-2020 G. Faccanoni 45
3 Primitives et Intégrales Mercredi 4 septembre 2019
R
R
0 0 dr si R < 9710 ¤
R 9710 RR
0 dr + 9710 0.14 dr si 9710 ¤ < R < 26 818 ¤
R0
RR
= 09710 0 dr + 9710
R 26818
0.14 dr + 26818 0.30 dr si 26 818 ¤ < R < 71 898 ¤
9710
R R 26818 R 71898 RR
0 0 dr + 9710 0.14 dr + 26818 0.30 dr + 71898 0.41 dr si 71 898 ¤ < R < 152 260 ¤
R 9710 0 dr + R 26818 0.14 dr + R 71898 0.30 dr + R 152260 0.41 dr + R R
si R > 152 260 ¤
0 9710 26818 71898 152260 0.45 dr
0
si R < 9710 ¤
0 + 0.14(R − 9710) si 9710 ¤ < R < 26 818 ¤
= 0 + 0.14(26818 − 9710) + 0.30(R − 26818) si 26 818 ¤ < R < 71 898 ¤
+ 0.14(26818 − 9710) + 0.30(71898 − 26818) + 0.41(R − 71898)
0 si 71 898 ¤ < R < 152 260 ¤
si R > 152 260 ¤
0 + 0.14(26818 − 9710) + 0.30(71898 − 26818) + 0.41(152260 − 71898) + 0.45(R − 152260)
0
si R ≤ 9710 ¤
0.14R − 1359.4 si 9710 ¤ < R ≤ 26 818 ¤
= 0.30R − 5650.28 si 26 818 ¤ < R ≤ 71 898 ¤
0.41R − 13559.06 si 71 898 ¤ < R ≤ 152 260 ¤
si R > 152 260 ¤
0.45R − 19649.46
Voici la représentation graphique de la fonction Impôt en fonction du revenu avec le revenu imposable en horizontale
et l’impôt à payer en verticale.
Impôt à payer
48 867.54 ¤
15 919.12 ¤
2395.12 ¤
Revenu
¤
8¤
8¤
¤
10
60
81
89
97
22
26
71
15
Elle est bien une fonction continue du revenu, croissante, affine par morceaux et convexe.
On peut écrire cette fonction comme suit :
³ ´
I (R) = min max(0, R, 9710) − 0 × 0
³ ´
+ min max(9710, R, 26818) − 9710 × 0.14
³ ´
+ min max(26818, R, 71898) − 26818 × 0.30
³ ´
+ min max(71898, R, 152260) − 71898 × 0.41
³ ´
+ max(R, 152260) − 152260 × 0.45
3. Le taux d’imposition que l’on retient en général est celui de la tranche de revenu dans laquelle on est : le célibataire
qui qui en 2017 déclare 100 000 ¤ de revenu pour l’année 2016 se trouve dans la tranche à 41%. Pourtant, il s’agit là
d’un taux marginal et non pas de l’impôt réellement payé. Ce qui compte, c’est plutôt ce qu’on paye vraiment à la
46 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
fin, au total, sur l’ensemble de son revenu ou QF. Le taux réel d’imposition (ou taux effectif) désigne le pourcentage
d’impôts réellement payé. Dans notre exemple, l’impôt payé est de 0 × (9710 − 0) + 14% × (26818 − 9710) + 30% ×
(71898 − 26818) + 41% × (100000 − 71898) = 27 440.94 ¤ pour un revenu de 100 000 ¤, soit un taux réel d’imposition de
27440.94/100000 ' 27.44% (et non 41%). Depuis quelques années le taux effectif d’imposition est indiqué sur l’avis
d’imposition.
Pour calculer le Taux réel d’imposition, c’est-à-dire le pourcentage qu’il faut appliquer au revenu pour avoir l’impôt,
il s’agit simplement de diviser I par R. Par exemple, si le revenu est de 50 000 ¤, l’impôt est de 9349.72 ¤, si bien que
le taux global est de 9349.72/50000 = 0.1869944, c’est-à-dire ' 18.70%.
Taux réel
32.09%
22.14%
8.93%
Revenu
¤
8¤
¤
10
98
60
81
97
22
26
71
15
Il s’agit de la fonction G : R+ → R+ définie par
0 si R ≤ 9710 ¤
0.14 − 1359.40
si 9710 ¤ < R ≤ 26 818 ¤
R
I (R)
G(R) = = 0.30 − 5650.28
R si 26 818 ¤ < R ≤ 71 898 ¤
R
0.41 − 13559.06 si 71 898 ¤ < R ≤ 152 260 ¤
R
0.45 − 19649.46
si R > 152 260 ¤
R
Chaque changement de tranche se repère ici par un point anguleux dans la courbure de la fonction.
4. Ce dernier graphique montre une autre propriété qui mérite d’être signalée car elle n’est pas évidente pour tout le
monde. Horizontalement, toujours le revenu. Verticalement, on indique «ce qui reste quand on a payé les impôts»,
autrement dit la différence entre le revenu et l’impôt. Eh bien, cette fonction est encore croissante. Qu’est-ce que cela
signifie ? Tout simplement que plus on gagne et plus on est riche, même après avoir déduit les impôts.
Ce qui reste après avoir
payé l’impôt id
103 392.46 ¤
55 978.88 ¤
24 422.88 ¤
9710 ¤
Revenu
0¤
8¤
0¤
18
1
89
26
97
8
26
71
2
15
© 2019-2020 G. Faccanoni 47
3 Primitives et Intégrales Mercredi 4 septembre 2019
Correction
1. Il semble assez évident de délimiter correctement de manière intuitive les zones par deux disques concentriques.
O A B
Attention néanmoins à ne pas répondre que la probabilité de choisir un point dans la zone hachurée est 1/2 (parce
que le rayon du disque hachuré correspond à la moitié du rayon du grand disque). Or, si le petit disque possède un
rayon de r et le grand 2r , on a
2. L’ensemble de points équidistants d’un point et d’une droite est une parabole. La région à considérer est donc
délimitée par quatre paraboles qui ont pour foyer le centre du carré et comme droites directrices les droites qui
supportent les côtés du carré. On s’affaire donc à trouver l’aire de cette région hachurée. On place d’abord le tout
dans un repère cartésien. Les sommets du carré sont (1, 1), (−1, 1), (−1, −1) et (1, −1) (cf. figure à gauche). En vertu des
symétries de la figure, il nous est possible de nous concentrer seulement sur la partie située dans le premier quadrant.
Qui plus est, il est possible de ne s’attarder qu’à la moitié de cette dernière région (cf. figure à droite).
(0, 1) (1, 1)
(−1, 1) (0, 1) (1, 1)
(0, 1
2)
(0, 1 C
(−1, 1) 1 , 1)
(− 2 (0, 1) (1 (1, 1) 2)
2 , 1)
(0, − 1
2)
Ce «croissant de parabole» correspond à la moitié de la région à considérer dans le premier quadrant. On note au
passage que l’aire du carré dans ce premier quadrant est 1. Le «croissant» est à son tour divisé en deux parties : la zone
de forme triangulaire T et la zone A . On cherche l’aire de ces zones. Pour y arriver, on aura besoin des coordonnées
2
de C . L’équation de la parabole qui
p nous intéresse est y = 1 − 2x. On ne s’intéressera qu’à la branche située au dessus
de l’axe des abscisses. Ainsi
p y = 1 − 2x. On cherche ensuite les coordonnées de C p, le point d’intersection entre la
courbe d’équation y = 1 − 2x et la droite d’équation y = x et on p
obtient x = −1 + 2. Il nous est donc déjà possible
3−2 2
de trouver l’aire du triangle, que l’on a identifié comme T = 2 . Il reste à trouver l’aire de la région sous la courbe :
1 p
p ´3 (−1 + 2)3
r
p 1 0 p 1
Z Z
2
³
A= p 1 − 2x dx = − p t dt = 3−2 2 = .
−1+ 2 2 3−2 2 3 3
48 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
h
Correction
1. Calcul de la capacité de la cuve.
Notons L la longueur de la cuve, R son rayon et h la hauteur du mazout. On a 0 ≤ h ≤ 2R. Considérons une coupe
verticale de la cuve :
R
h
Pour calculer l’aire coloriée, on “tourne” le cercle et on le plonge dans un repère orthonormé :
y p
x(2R − x)
h1 h2 x
Z p x=R−y Z q
dx=− dy
x(2R − x) dx = − R 2 − y 2 dy
y=R sin(t ) Z
dy=R cos(t ) dt
= −R 2 cos2 (t ) dt
R2
=− (sin(t ) cos(t ) + t ) + c
2
2µ
R
q ¶
=− sin(t ) 1 − sin2 (t ) + t + c
2
© 2019-2020 G. Faccanoni 49
3 Primitives et Intégrales Mercredi 4 septembre 2019
s
R2 y
¡y¢
t =arcsin R y2 ³y´
+c
= − 1 − 2 + arcsin
2 R R R
y=R−x R −xp R2 ³ x´
= − x(2R − x) − arcsin 1 − + c.
2 2 R
Par conséquent
Z h2 p
A =2 x(2R − x) dx ≈ 1.48 m2 .
h1
? f 0 (h) = 0 ssi h = 0 ou h = 2R ;
h−R
? f 00 (h) = 2000 p ;
h(2R−h)
Correction
On note L(t ) la longueur de l’élastique à l’instant t et `(t ) la distance parcourue par l’escargot à l’instant t . Pour les unités
de mesure, on convient qu’une unité de temps correspond à un jour et les longueurs sont mesurées en mètres. On a
L(0) = 1000, `(0) = 0 et il s’agit de voir si `(t ) = L(t ) pour un certain t .
Pour tout t ≥ 0,
L(t ) = 1000t + 1000
et on peut définir y(t ) la fraction de l’élastique parcourue par l’escargot à l’instant t :
`(t )
y(t ) = ∀t ≥ 0.
L(t )
La vitesse `0 de l’escargot par rapport à l’extrémité fixe de l’élastique est la somme de deux vitesses : la vitesse de l’escargot
sur l’élastique, soit 1 mètre par jour, et la vitesse du point de l’élastique où se trouve l’escargot (on peut faire l’hypothèse
que cette vitesse est proportionnelle à l’abscisse de l’escargot) :
`0 (t ) = 1 + y(t )L 0 (t ) ∀t ≥ 0,
donc
`0 (t ) L 0 (t ) 1 + y(t )L 0 (t ) L 0 (t ) 1
y 0 (t ) = − y(t ) = − y(t ) = ∀t ≥ 0.
L(t ) L(t ) L(t ) L(t ) L(t )
Puisque y(0) = 0, on en conclut que
Z t 1 1
Z t 1 1 1
y(t ) = dτ = dτ = [ln(1 + τ)]0t = ln(1 + t )
0 L(τ) 1000 0 1+τ 1000 1000
50 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 3 Primitives et Intégrales
et donc
1000t + 1000
`(t ) = y(t )L(t ) = ln(1 + t ) = (1 + t ) ln(1 + t ) ∀t ≥ 0.
1000
L’escargot touchera l’extrémité mobile de l’élastique lorsque `(t ) = L(t ), c’est-à-dire à l’instant t f = e 1000 − 1 ' 1.97 × 10434
jours ' 5.397 × 10431 années (ce qui correspond à ' 3.9 × 10421 fois l’age de l’univers).
En étudiant la fonction t 7→ d (t ) = L(t ) − `(t ), on trouve que cette distance est maximale 2 à l’instant t 0 = e 999 − 1 ; après cet
instant l’escargot commence à se rapprocher de l’extrémité de l’élastique pour en arriver au but à l’instant t f = e 1000 − 1. À
l’instant t 0 l’escargot se déplace à une vitesse de 1000 kilomètre par jour et a parcouru y(t 0 ) = 99.9% de l’élastique, elle a
parcouru 99.9% de l’élastique mais elle n’a jamais été aussi loin de son but !
2. d 0 (t ) = 999 − ln(1 + t )
© 2019-2020 G. Faccanoni 51
Équations différentielles ordinaires (EDO)
4
Les équations différentielles décrivent l’évolution de nombreux phénomènes dans des domaines variés. Une équation
différentielle est une équation impliquant une ou plusieurs dérivées d’une fonction inconnue. Si toutes les dérivées sont
prises par rapport à une seule variable, on parle d’équation différentielle ordinaire (EDO). Une équation mettant en jeu des
dérivées partielles est appelée équation aux dérivées partielles (EDP).
Une EDO est une équation exprimée sous la forme d’une relation
F (y(t ), y 0 (t ), y 00 (t ), . . . , y (p) (t )) = g (t )
q 0 (t ) = αq(t ).
La désintégration atomique est un cas de décroissance régi par la même équation mais avec α < 0.
Hypothèse de Verhulst : à chaque instant, la croissance de la population est «proportionnelle» à son effectif, mais
inhibée par des ressources limitées :
q 0 (t ) = αq(t )(m − q(t )).
Il est clair que m sera un point d’équilibre, car la dérivée de q est nulle quand q = m.
Hypothèse de Gompertz : à chaque instant, la croissance de la population est «proportionnelle» à son effectif, mais
inhibée par des ressources limitées :
q 0 (t ) = αq(t )(ln(k) − ln(q(t )).
Il est clair que k sera un point d’équilibre, car la dérivée de q est nulle quand q = k.
Si la fonction g , appelée «second membre» de l’équation, est nulle, on dit que l’équation en question est homogène.
Résoudre une équation différentielle, c’est chercher toutes les fonctions, définies sur un intervalle I ⊂ R, qui satisfont
l’équation (on dit aussi intégrer l’équation différentielle).
E XEMPLE
Résoudre l’équation différentielle y 0 (t ) = −y(t ) signifie chercher toutes les fonctions
y: I ⊂R→R
t 7→ y = f (t )
telles que f 0 (t ) = − f (t ) pour tout t ∈ I . On peut vérifier que y(t ) = ce −t pour tout t ∈ R (où c est constante réelle quelconque)
est une solution de l’EDO (en particulier, pour c = 0 on trouve la solution nulle).
53
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
Une EDO admet généralement une infinité de solutions. Pour choisir, entre les différentes solutions, celle qui décrit le
problème physique, il faut considérer d’autres données qui dépendent de la nature du problème, par exemple la valeur
prise par la solution et/ou éventuellement ses dérivées en un ou plusieurs points de l’intervalle d’intégration.
y(t )
µ ¶
y 0 (t ) = C y(t ) 1 − (4.1)
B
En pratique, se donner une CI revient à se donner le point (t 0 , y 0 ) par lequel doit passer le graphe de la fonction solution et
la valeur de ses dérivées en ce même point.
Le couple EDO-CI porte le nom de problème de C AUCHY ou de problème aux valeurs initiales :
Si ϕ ne dépend pas explicitement de t (i.e. si ϕ(t , y(t )) = ϕ(y(t ))), l’EDO est dite autonome.
L’essentiel de notre analyse concernera le cas où l’on a qu’une seule EDO, c’est-à-dire le cas scalaire. L’extension aux
systèmes sera faite au prochaine semestre.
Résoudre un problème de C AUCHY, c’est chercher toutes les fonctions, définies sur un intervalle I ⊂ R, qui satisfont l’équation
et qui vérifient la condition initiale. On aura donc des questions naturelles telles
? trouver toutes les fonctions solutions de l’EDO,
? parmi toutes ces fonctions, choisir celles qui vérifient la CI (existence ? unicité ?),
? parmi toutes ces fonctions, étudier le domaine de définition (pour chaque fonction trouvée, quel est le plus grande
domaine de définition qui contient le point t 0 ?)
Sa solution, définie sur R, est donnée par y(t ) = (α + 1/3)e −3t + t − 1/3. En effet on a bien
y(0) = (α + 1/3)e 0 + 0 − 1/3 = α, y 0 (t ) = −3(α + 1/3)e −3t + 1 = −3(α + 1/3)e −3t + 1 − 3t + 3t = −3y(t ) + 3t .
54 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
Cet exemple montre le cas où il existe une et une seule solution du problème de C AUCHY définie sur R. Les choses ne se
passent pas toujours si bien. Les exemples ci-dessous montrent que l’étude mathématique de l’existence et de l’unicité des
solutions d’un problème de C AUCHY peut être une affaire délicate.
E XEMPLE (E XISTENCE ET UNICITÉ SUR I ⊂ R ( MAIS NON EXISTENCE SUR R) DE LA SOLUTION D ’ UN PROBLÈME DE C AUCHY )
On se donne ϕ(t , y(t )) = (y(t ))3 et y 0 = 1. On cherche une fonction y : t ∈ R+ 7→ y(t ) ∈ R qui satisfait
(
y 0 (t ) = (y(t ))3 , ∀t > 0,
y(0) = 1.
On vérifie que la solution y est donnée par y(t ) = p 1 qui n’est définie que pour t ∈ [0; 1/2[. Cet exemple montre qu’un
1−2t
problème de C AUCHY n’a pas toujours une solution pour tout t ∈ [0; +∞[ puisqu’ici la solution explose lorsque t tend vars la
valeur 1/2 (en effet, nous avons lim − y(t ) = +∞) : le graphe de la solution a une asymptote verticale en t = 1/2. On parle
t →(1/2)
d’explosion de la solution en temps fini ou encore de barrière.
On verra que ceci est un phénomène général : pour une solution d’une EDO, la seule façon de ne pas être définie sur R est
d’avoir un asymptote verticale.
De façon générale, lorsqu’on se donne une équation différentielle et une condition initiale y(t 0 ) = y 0 , on cherche un
intervalle I , contenant t 0 , sur lequel une solution existe, et qui soit «le plus grand possible» : il n’existe pas d’intervalle
plus grand sur lequel l’équation différentielle ait une solution. Cet intervalle s’appelle intervalle de vie de la solution. Une
solution définie sur cet intervalle le plus grand possible s’appelle solution maximale.
Il y a un résultat qui garantit que, sous certaines hypothèses très générales, deux graphes de fonctions qui sont des solutions
de la même EDO ne se rencontrent jamais. Le théorème garantit aussi l’existence des solutions. Dans ce cours, nous nous
contentons de rappeler un résultat d’existence et d’unicité global, au sens où on peut intégrer le problème de C AUCHY
jusqu’à t = ∞, pour des EDO linéaires.
4.6 Théorème
Soit y = f (t ) une solution maximale définie sur un intervalle de vie I =]a; b[. Si b 6= +∞ alors lim y(t ) = ∞, i.e. le graphe de
t →b−
la solution a une asymptote verticale en t = b. Même chose si a 6= −∞.
On utilise souvent le théorème sous forme contraposée : si les solutions ne peuvent pas «exploser», alors elles sont définies
sur R.
© 2019-2020 G. Faccanoni 55
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
Remarque
D’un point de vue pratique, cet énoncé nous aidera à faire des dessins, en garantissant que les graphes des solutions ne se
rencontrent jamais. On peut en déduire quelques remarques plus subtiles :
? si l’EDO admet comme solution la solution nulle mais y 0 6= 0, alors la solution du problème de C AUCHY est du signe
de y 0 pour tout t ∈ I ;
? si l’EDO admet deux solutions constantes y(t ) = κ1 et y(t ) = κ2 pour tout t ∈ I et y 0 ∈]κ1 ; κ2 [, alors la solution du
problème de C AUCHY vérifie y(t ) ∈]κ1 ; κ2 [ pour tout t ∈ R.
F (y(x)) = G(x) +C où C ∈ R,
y(x) = F −1 (G(x) +C ),
relation qui donne toutes les solutions de l’équation. Cette solution générale dépend de la constante d’intégration C .
E XEMPLE
On veut résoudre l’équation différentielle y 0 (x) = x y(x). Il s’agit d’une EDO du premier ordre à variables séparables :
? Recherche des solutions constantes. Si y(x) = A pour tout x alors y 0 (x) = 0 pour tout x et l’EDO devient 0 = x A pour
tout x. Par conséquent A = 0 : la fonction y(x) = 0 pour tout x est l’unique solution constante de l’EDO.
? Recherche des solutions non constantes. La fonction y(x) = 0 pour tout x étant solution, toute autre solution x 7→ y(x)
sera donc non nulle. On peut alors diviser l’EDO par y et procéder formellement comme suit :
y 0 (x) 1 x2
Z Z
=x =⇒ dy = x dx =⇒ ln |y| = +C avec C ∈ R.
y(x) y 2
56 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
avec
b(x)
? A(x) une primitive de ,
a(x)
g (x) A(x)
? K (x) une primitive de e .
a(x)
On peut montrer que
? y H est la solution générale de l’EDO homogène associée, c’est-à-dire de l’EDO a(x)y 0 (x) + b(x)y(x) = 0 (qui est à
variables séparables) ;
? y P est une solution particulière 1 de l’EDO a(x)y 0 (x) + b(x)y(x) = g (x) ;
? u est la solution générale de l’EDO a(x)u 0 (x)v(x) = g (x) ;
? v est une solution particulière non nulle de l’EDO homogène associée, c’est-à-dire de l’EDO a(x)v 0 (x) + b(x)v(x) = 0.
E XEMPLE
Considérons l’EDO
y 0 (x) − y(x) = x.
On a a(x) = 1, b(x) = −1 et g (x) = e x , donc pour x ∈ R on a
? A(x) = −1 dx = −x,
R
ce qui donne
e x −(1 + x) = C − (1 + x)e −x |{z}
¡ ¢ x
y(x) = C
|{z} e .
| {z } | {z }
yH yP u v
E XEMPLE (L OI DE N EWTON K)
Considérons une tasse de café à la température de 75◦C dans une salle à 25◦C. Après 5 minutes le café est à 50◦C. Si
on suppose que la vitesse de refroidissement du café est proportionnelle à la différence des températures (i.e. que la
température du café suit la loi de Newton), cela signifie qu’il existe une constante γ < 0 telle que la température vérifie
l’EDO du premier ordre
T 0 (t ) = γ(T (t ) − 25)
avec la CI
T (5) = 50,
ayant convenu qu’une unité de temps correspond à une minute et la température est mesuré en degré Celsius.
1. On commence par calculer toutes les solutions de l’EDO. Étant une équation différentielle du premier ordre, la
famille de solutions dépendra d’une constante qu’on fixera en utilisant la CI. Si on réécrit l’EDO sous la forme
T 0 (t ) − γT (t ) = −25γ, on a une EDO linéaire d’ordre 1 avec a(t ) = 1, b(t ) = −γ et g (t ) = −25γ. Donc
? A(t ) = −γ dt = −γt ,
R
© 2019-2020 G. Faccanoni 57
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
ln(2) ln(2)
50 = T (5) = 25 + 50e γt =⇒ γ=− =⇒ T (t ) = 25 + 50e − 5 t
5
T [◦C]
75
50
25
0 t [minutes]
5 10 15
où u, v et w sont des fonctions données, continues sur un intervalle I ⊂ R. Pour la résolution, on se place sur un intervalle
J ⊂ I tel que la fonction u ne s’annule pas sur J et on définit une nouvelle fonction x 7→ z(x) = (y(x))1−α . L’EDO initiale est
alors équivalente à l’EDO linéaire du premier ordre suivante : 2
Par conséquent, pour x ∈ Dv ∩ Dw ∩ { x ∈ Du | u(x) 6= 0 }, toute solution y s’écrit comme y(x) = [z(x)]1/(1−α) avec
−A(x)
? z(x) = C
| e {z } + K (x)e −A(x) = (C + K (x)) e| −A(x) ,
| {z } | {z } {z }
y H (x) y P (x) u(x,C ) v(x)
b(x) v
? A(x) une primitive de i.e. une primitive de (1 − α) ,
a(x) u
g (x) A(x) w
? K (x) une primitive de e i.e. une primitive de (1 − α) e A(x) .
a(x) u
E XEMPLE
On se propose de résoudre l’équation différentielle
1 1
y 0 (x) + y(x) = (x − 1)y 3 (x).
2 2
Il s’agit d’une équation différentielle de B ERNOULLI. Comme u(x) = 1 pour tout x ∈ R, on cherche sa solution générale sur R.
v(x) 1/2
Z Z
? A(x) = (1 − α) dx = −2 dx = −x,
u(x) 1
w(x) A(x) (x − 1)/2 −x
Z Z Z Z
? K (x) = (1 − α) e dx = −2 e dx = (1 − x)e −x dx = −(1 − x)e −x − e −x dx = xe −x ,
u(x) 1
? z(x) = (C + K (x)) e −A(x) = (C + xe −x )e x = C e x + x,
2. Formellement z = y 1−α implique d’une part y = z y α et d’autre part z 0 = (1 − α)y −α y 0 et donc y 0 = (1 − α)z 0 y α .
58 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
1
y(x) = p .
x +C e x
où a, b et c sont des constantes données (a 6= 0) et g est une application continue sur un intervalle I de R.
Toute solution y d’un EDO linéaire du second ordre à coefficients constants dépend de deux constantes arbitraires C 1 et C 2
et est de la forme y H (x) + y P (x) où y P est une solution particulière de l’EDO et y H est la solution générale de l’équation
homogène associée (c’est-à-dire de l’EDO a y 00 (x) + b y 0 (x) + c y(x) = 0). On doit donc résoudre deux problèmes : chercher
d’abord la solution générale de l’équation homogène et ensuite une solution particulière de l’équation complète.
? Résolution de l’équation homogène associée.
On introduit le polynôme caractéristique p(λ) = aλ2 + bλ + c. Soit ∆ = b 2 − 4ac, alors
? si ∆ > 0 on a
p p
λ1 x λ2 x −b − ∆ −b + ∆
y H (x) = C 1 e +C 2 e , C 1 ,C 2 ∈ R, λ1 = , λ2 = ;
2a 2a
? si ∆ = 0 on a
b
y H (x) = (C 1 +C 2 x)e λx , C 1 ,C 2 ∈ R, λ=− ;
2a
? si ∆ < 0 on a
p
σx b |∆|
y H (x) = e (C 1 cos(ωx) +C 2 sin(ωx)), C 1 ,C 2 ∈ R, σ=− , ω= ,
2a 2a
qu’en physique souvent on réécrit comme
q C1 C2
y H (x, A, ϕ) = Ae σx cos(ωx − ϕ), A= C 12 +C 22 , cos(ϕ) = , sin(ϕ) = .
A A
E XEMPLE
Soit m un paramètre qui dépend du polynôme caractéristique.
© 2019-2020 G. Faccanoni 59
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
? Si g (x) = cos(5x) ou g (x) = sin(5x) alors n = 0, µ = 0 et θ = 5 donc on cherchera y P sous la forme y P (x) = x m (A cos(5x)+
B sin(5x)).
? Si g (x) = e 2x sin(5x) alors n = 0, µ = 2 et θ = 5 donc on cherchera y P sous la forme y P (x) = x m e 2x (A cos(5x) +
B sin(5x)).
? Si g (x) = x cos(5x) ou g (x) = x sin(5x) alors n = 1, µ = 0 et θ = 5 donc on cherchera y P sous la forme y P (x) =
x m ((Ax + B ) cos(5x) + (C x + D) sin(5x)).
? Si g (x) = x alors n = 1, µ = 0 et θ = 0 donc on cherchera y P sous la forme y P (x) = x m (Ax + B ).
? Si g (x) = xe 3x alors n = 1, µ = 3 et θ = 0 donc on cherchera y P sous la forme y P (x) = x m e 3x (Ax + B ).
? Si g (x) = e 2x alors n = 0, µ = 2 et θ = 0 donc on cherchera y P sous la forme Ax m e 2x .
E XEMPLE
On veut calculer toutes les solutions de l’EDO
On a alors
y P00 (x) + y P (x) = 3 cos(x) =⇒ (2β − αx) cos(x) − (2α + βx) sin(x) + x(α cos(x) + β sin(x)) = 3 cos(x)
d’où α = 0 et β = 23 .
L’intégrale générale de l’EDO assignée est donc
3
y(x) = y H (x) + y P (x) = c 1 cos(x) + c 2 sin(x) + x cos(x), c 1 , c 2 ∈ R.
2
60 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
p(λ) = λ2 + 2δλ + η2
qui a discriminant
∆ = 4δ2 − 4η2
et racines p p
−2δ − ∆ −2δ + ∆
q q
λ1 = = −δ − δ2 − η2 et λ2 = = −δ + δ2 − η2 .
2 2
Selon le signe de ∆ on a trois comportements différentes :
? Si ∆ > 0, c’est-à-dire si δ > η alors λ1 et λ2 sont réels et différents et la solution de l’EDO est de la forme
x(t ) = C 1 e λ1 t +C 2 e λ2 t , C 1 ,C 2 ∈ R.
Comme λ1 < λ2 < 0 (car δ2 − η2 < δ), x tend vers 0 de façon exponentielle quand t → +∞. Physiquement cela
p
signifie que si la constante de frottement est grande comparée à la constate d’élasticité du ressort alors la masse
n’oscille pas mais va être tirée vers la position d’équilibre.
? Si ∆ = 0, c’est-à-dire si δ = η alors λ1 = λ2 = −δ et la solution de l’EDO est de la forme
x(t ) = (C 1 +C 2 t )e λt , C 1 ,C 2 ∈ R.
Dans ce cas aussi la solution x tend vers 0 de façon exponentielle quand t → +∞.
? Si ∆ < 0, c’est-à-dire si δ < η alors λ1 et λ2 sont deux nombres complexes conjugués et la solution de l’EDO est de la
forme p p
x(t ) = C 1 e −δt cos( −∆t ) +C 2 e −δt sin( −∆t ), C 1 ,C 2 ∈ R.
p q
qui se réécrit x(t ) = r e −δt cos( −∆t + ϕ) avec r = C 12 +C 22 , ϕ = arctan(−C 1 /C 2 ). Dans cette dernière expression on
voit le caractère oscillatoire du mouvement. Dans ce cas le frottement ne suffit pas pour empêcher l’oscillation mais
son effet se traduit par une diminution exponentielle de l’ampleur de l’oscillation : le graphe de x(t ) est compris
entre les courbes d’équation ±r e −δt .
1,0
1,0
1,0
0,5
0,5
0,5
0 0 0
2 4 6 8 10 2 4 6 8 10 2 4 6 8 10
t t t
K0,5
K0,5
K0,5
x(0)=1 x´'(0)=1 x(0)=1 x´'(0)=-1 x(0)=1 x´'(0)=1 x(0)=1 x´'(0)=-1 x(0)=1 x´'(0)=1 x(0)=1 x´'(0)=-1
x(0)=1 x´'(0)=0 x(0)=-1 x´'(0)=1 x(0)=1 x´'(0)=0 x(0)=-1 x´'(0)=1 x(0)=1 x´'(0)=0 x(0)=-1 x´'(0)=1
x(0)=-1 x´'(0)=-1 x(0)=-1 x´'(0)=0 x(0)=-1 x´'(0)=-1 x(0)=-1 x´'(0)=0 x(0)=-1 x´'(0)=-1 x(0)=-1 x´'(0)=0
mx 00 (t ) + kx(t ) = A cos(ϕt ).
© 2019-2020 G. Faccanoni 61
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
? Équation homogène : x 00 (t ) + η2 = 0. Le polynôme caractéristique est p(λ) = λ2 + η2 et a les deux racines complexes
conjugués λ1 = −i η et λ2 = i η. La solution est de la forme
x P (t ) = b cos(ϕt ) + c sin(ϕt ).
On a x P0 (t ) = −bϕ sin(ϕt ) + cϕ cos(ϕt ) et x P00 (t ) = −bϕ2 cos(ϕt ) − cϕ2 sin(ϕt ). En remplaçant dans l’EDO on
obtient
−bϕ2 cos(ϕt ) − cϕ2 sin(ϕt ) + η2 b cos(ϕt ) + η2 c sin(ϕt ) = a cos(ϕt ),
ce qui implique
a
b= , et c = 0.
η2 − ϕ 2
x P (t ) = bt cos(ϕt ) + c t sin(ϕt ).
On a x P0 (t ) = (b + cϕt ) cos(ϕt ) + (c − bϕt ) sin(ϕt ) et x P00 (t ) = (2c − bϕt )ϕ cos(ϕt ) − (2b + cϕt )ϕ sin(ϕt ). En rem-
plaçant dans l’EDO on obtient
qui se réécrit
((−2 − t )b + (1 − t )c)ϕ sin(ϕt ) + ((b − a) + 2c(1 + t )ϕ − bt ϕ2 ) cos(ϕt ) = 0
ce qui implique
a
b = 0, et c= .
2ϕ
62 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
Correction
Il s’agit d’une EDO à variables séparables. La fonction y(x) = 0 pour tout x est solution de l’EDO mais elle ne vérifie pas la
CI. Toute autre solution de l’EDO sera non nulle et se trouve formellement comme suit :
y 0 (x) 1
Z Z
y 0 (x) + 2x y 2 (x) = 0 =⇒ 2
= −2x =⇒ y −2
dy = −2 x dx =⇒ y(x) = 2 , c ∈ R.
y (x) x +c
2
En imposant la CI on obtient 2 = 1/C d’où l’unique solution du problème de C AUCHY : y(x) = 2x 2 +1
.
Exercice 4.2
Résoudre le problème de C AUCHY
(
y 0 (x) − 4x y 2 (x) = 0,
y(0) = 2.
Correction
Il s’agit d’une EDO à variables séparables. La fonction y(x) = 0 pour tout x est solution de l’EDO mais elle ne vérifie pas la
CI. Toute autre solution de l’EDO est non nulle et se trouve formellement comme suit :
y 0 (x) 1
Z Z
y 0 (x) − 4x y 2 (x) = 0 =⇒ 2
= 4x =⇒ y −2
dy = 4 x dx =⇒ y(x) = , c ∈ R.
y (x) −2x 2 + c
2
En imposant la CI on obtient 2 = 1/C d’où l’unique solution du problème de C AUCHY y(x) = 1−4x 2
.
Exercice 4.3
Résoudre le problème de C AUCHY
(
y 0 (t ) = t y 2 (t ),
y(0) = y 0 ,
en fonction de la donnée initiale y 0 .
Correction
Il s’agit d’un problème de C AUCHY avec une CI y(0) = y 0 et une EDO du premier ordre à variables séparable.
On cherche d’abord les solutions constantes, i.e. les fonctions y(x) ≡ A ∈ R qui vérifient l’EDO, c’est-à-dire qui vérifient
0 = t A 2 pour tout y ∈ R ; l’unique solution constante est donc la fonction y(x) ≡ 0.
Comme deux trajectoires ne s’intersectent pas, toutes les autres solution ne s’annulent jamais. Soit donc y(x) 6= 0 ; on peut
alors écrire
y 0 (t ) 1 1 1 t2 1
Z Z
2
= t =⇒ 2
dy = t dt =⇒ 2
dy = t dt =⇒ − = +C =⇒ y(x) = − 2 , pour C ∈ R.
y (t ) y y y 2 t
2 +C
© 2019-2020 G. Faccanoni 63
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
Correction
1. Il s’agit d’une «EDO du premier ordre à variables séparables». Si n(t ) ≡ c est solution alors 0 = −kc d’où c = 0 : l’unique
solution constante est la solution n(t ) = 0 quelque soit t ∈ R+ .
Si n(t ) 6= 0, on peut écrire
n 0 (t )
= −k
n(t )
d’où la famille de solutions
n(t ) = De −kt , D ∈ R+ .
On conclut que, quelque soit la condition initiale n(0) = n 0 ≥ 0, l’unique solution est n(t ) = n 0 e −kt pour tout t ∈ R+ .
2. Puisque n 0 /2 = n(5700) = n 0 e −5700t , on obtient k = ln 2−5700 ≈ 1.216 10−4 .
3. Puisque n 0 /9 = n(t̂ ) = n 0 e −k t̂ , on obtient t̂ = 5700 ln 9
ln 2 ≈ 18000 ans.
Correction
La loi de Newton affirme qu’il existe une constante K < 0 telle que la température du corps suit l’EDO
T 0 (t ) = K (T (t ) − Text ).
On commence par calculer toutes les solutions de l’EDO. Étant une équation différentielle du premier ordre, la famille de
solutions dépendra d’une constante D qu’on fixera en utilisant la CI.
? On cherche d’abord les solutions constantes, i.e. les solutions du type T (t ) ≡ c ∈ R quelque soit t . On a
0 = K (c − Text )
d’où l’unique solution constante T (t ) ≡ Text . Toute autre solution T (t ) vérifiera T (t ) 6= Text pour tout t .
? Soit T (t ) 6= Text pour tout t . Puisqu’il s’agit d’une EDO à variables séparables on peut calculer la solution formellement
comme suit :
T 0 (t ) dT
T 0 (t ) = K (T (t ) − Text ) =⇒ =K =⇒ = K dt =⇒
T (t ) − Text T − Text
1
Z Z
dT = K dt =⇒ ln(T − Text ) = K t + c =⇒ T − Text = De K t =⇒ T (t ) = Text + De K t .
T − Text
T (t ) = −5 + 25e K t .
ln(4/5)
De plus, on sait que 15 = T (30) = −5+25e 30K d’où K = 30 . l’équation
ln(4/5)
La température du corps suit donc la loi 37 = −5 + 25e 30 t
T (t ) = −5 + 25e
ln(4/5)
30 t . d’où t = 30 ln(42/25)
ln(4/5) ∼ −69,7 minutes, c’est-à-dire à 1H10 de
la nuit.
Pour déterminer l’heure du meurtre il faut donc résoudre
64 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
T (◦C)
crime
40◦C
découverte
30◦C
du cadavre après 30
minutes
◦
20 C
10◦C
0 ◦C
-70 -50 -30 -10 10 30 t (minutes)
1H10 1H30 1H50 2H10 2H30 2H50
Correction
La loi de Newton affirme qu’il existe une constante K < 0 telle que la température du corps suit l’EDO
T 0 (t ) = K (T (t ) − Text ).
On commence par calculer toutes les solutions de l’EDO. Étant une équation différentielle du premier ordre, la famille de
solutions dépendra d’une constante D qu’on fixera en utilisant la CI.
? On cherche d’abord les solutions constantes, i.e. les solutions du type T (t ) ≡ c ∈ R quelque soit t . On a
0 = K (c − Text )
T 0 (t ) dT
T 0 (t ) = K (T (t ) − Text ) =⇒ =K =⇒ = K dt =⇒
T (t ) − Text T − Text
1
Z Z
dT = K dt =⇒ ln(T − Text ) = K t + c =⇒ T − Text = De K t =⇒ T (t ) = Text + De K t
T − Text
Ici l’inconnue est Text . On sait que T (t = 0) = 100◦C et ralement, la température du gâteau évolue selon la loi
T (t = 10) = 80◦C et T (t = 20) = 65◦C. Il s’agit donc de ré-
ln(3/4)
soudre le système de trois équations en les trois inconnues T (t ) = 20 + 80e 10 t .
K , D, Text :
100 = Text + D,
80 = Text + De 10K ,
65 = Text + De 20K .
© 2019-2020 G. Faccanoni 65
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
100◦C
80◦C
T (20) = 65◦C
70◦C
60◦C
50◦C
40◦C
10 20 30 t (minutes)
Exercice 4.7
Deux produits chimiques présents dans une cuve avec une concentration de 1g/l à l’instant t = 0 interagissent et produisent
une substance dont la concentration est notée y(t ) à l’instant t ≥ 0. On suppose que y(t ) est régie par l’équation différentielle
Correction
1. Pour montrer qu’une fonction est croissante il suffit de montrer que sa dérivée est de signe positif. Si y est solution de
l’EDO on a
y 0 (t ) = (1 − y(t ))2 ≥ 0 pour tout t ≥ 0
car un carré est toujours positif. y est donc une fonction croissante.
2. On cherche les fonctions constantes solution de l’EDO. Si f (t ) = c est solution de l’EDO alors puisque f 0 (t ) = 0 on
obtient
0 = (1 − c)2
soit c = 1. La seule fonction constante solution de l’EDO est la fonction constante égale à 1.
3. Considérons la solution y telle que y(0) = 0. Tout d’abord on a montré que la fonction y était croissante donc
y(0) ≤ y(t ) pour tout t ≥ 0, par conséquent, puisque 0 ≤ y(0), y(t ) ≥ 0 pour tout t ≥ 0. Supposons qu’il existe un t 0
tel que y(t 0 ) ≥ 1, alors le graphe de y qui relie continument les points (0, y(0)) et (t 0 , y(t 0 )) coupe nécessairement le
graphe de f , i.e. la droite d’équation y = 1. Ceci est impossible, car les graphes de deux solutions distinctes ne se
coupent jamais. Il n’existe donc pas de t 0 tel que y(t 0 ) ≥ 1, c’est-à-dire pour tout t ≥ 0, y(t ) < 1.
4. Considérons la solution y telle que y(0) = 0.
La fonction y est croissante et majorée par 1, elle admet donc une limite pour t → +∞. On note lim y(t ) = `. On
t →+∞
suppose que lim y 0 (t ) = `. En passant à la limite dans l’EDO on obtient :
t →+∞
0 = (1 − `)2
soit ` = 1.
5. ? Si y(0) = 1 on sait que y(t ) = 1 pour tout t > 0.
? Si y(0) = 0 on sait que la fonction y est croissante et lim y 0 (t ) = 1. En effet, il s’agit d’une EDO à variables
t →+∞
séparables et on peut écrire
t +c −1
Z
(1 − y)−2 dy = t , i.e. y(t ) =
t +c
66 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
t
y(t ) = , ∀t ≥ 0.
1+t
? Si y(0) = 2 on sait que la fonction y est croissante mais elle n’existe que pour 0 < t < 1 et on a
t −2
y(t ) = .
t −1
y
0
x
1 2 3
Ce modèle, qui décrit l’évolution d’une population de p individus à l’instant t , suppose que le taux de croissance du nombre
d’individus n’est pas constant mais diminue si la population augment (les ressources se réduisent).
Correction
On doit résoudre l’EDO à variables séparables
k
0 = A(k − h A) ⇐⇒ A = 0 ou A = .
h
On trouve ainsi deux solutions constantes :
k
p(t ) ≡ 0 et p(t ) ≡ .
h
p 0 (t )
= 1;
p(t )(k − hp(t ))
© 2019-2020 G. Faccanoni 67
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
p0 1
D= = k
.
k − hp 0 −h
p0
Remarquons que limt →+∞ p(t ) = hk : une population qui évolue à partir de p 0 individus à l’instant initiale selon la loi
logistique tend à se stabiliser vers un nombre d’individus d’environ k/h, ce qui représente la capacité de l’environnement.
D’autre part, déjà en analysant l’EDO on aurait pu déduire que les solutions sont des fonctions strictement croissantes si
p(t ) ∈]0, k/h[, décroissantes si p(t ) > k/h.
p
k
p(t ) > h
k
h
k
0 < p0 < h
0
t
Correction
On a le problème de C AUCHY
(
I 0 (t ) = k I (t )(5000 − I (t )), (EDO)
I (0) = 160. (CI)
Vu la nature de la question on ne s’intéresse qu’aux solutions positive et que pour t > 0.
1. Tout d’abord on observe qu’il y a deux solutions constantes de l’EDO : la fonction I (t ) ≡ 0 et la fonction I (t ) ≡ 5000.
2. Pour chercher toutes les solutions non constantes on remarque qu’il s’agit d’une EDO à variables séparables donc
formellement on a
I 0 (t )
I 0 (t ) = k I (t )(5000 − I (t )) =⇒ =k =⇒
I (t )(5000 − I (t ))
dI 1
Z Z
= k dt =⇒ d I = k dt =⇒
I (5000 − I ) I (5000 − I )
1 1
Z Z Z
dI − d I = 5000k dt =⇒ ln(I ) + ln(5000 − I ) = 5000kt + c =⇒
I 5000 − I
I I
ln = 5000kt + c =⇒ = De 5000kt =⇒
5000 − I 5000 − I
5000De 5000kt 5000
I (t ) = =⇒ I (t ) =
1 + De 5000kt De −5000kt +1
68 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
4. Il ne reste qu’à établir la valeur numérique de la constante k grâce à l’information sur le nombre d’individus infectés
après 7 jours :
80×5000
5. On cherche t̄ tel que I (t̄ ) = 80%I T = 100 = 4000 :
20000
4000 = t̄
¡ 363 ¢
4 + 121e − 7 ln 38
1
d’où t̄ = 5000 ln(121) ≈ 15 jours.
6. Avec ce modèle lim I (t ) = 5000 mais I ne peut jamais atteindre exactement 100% de la population en un temps fini
t →+∞
(deux solution ne s’intersectent jamais).
I
5000
4000
1200
160
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 t
Exercice 4.10
On note y(t ) le nombre de ménages vivant en France équipés d’un ordinateur (t est exprimé en années et y(t ) en millions
de ménages). Le modèle de VARHULST estime que sur la période 1980 − 2020, y(t ) est solution de l’équation différentielle
Correction
1. On doit résoudre l’EDO à variables séparables
On cherche d’abord les solutions constantes, i.e. des fonctions y(t ) = A pour tout t ∈ R :
0 = 0, 022A(20 − A) ⇐⇒ A = 0 ou A = 20.
y 0 (t )
= 0, 022;
y(t )(20 − y(t ))
© 2019-2020 G. Faccanoni 69
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
et on en déduit
20
y(t ) = pour tout D ∈ R+ .
1 + 20De −0,44t
2. Si t = 0 correspond à l’année 1980 et si y(0) = 0, 01 alors
20
0, 01 = =⇒ D = 1999
1 + 20De −0,44×0
et la fonction qui estime le nombre de ménages en France équipés d’un ordinateur t années après 1980 est
20
y(t ) = .
1 + 1999e −0,44t
Pour prévoir combien de ménages vivant en France seront équipés d’un ordinateur en 2020 il suffit de calculer y(40)
20
y(40) = ≈ 19.99.
1 + 1999e −0,44×40
Calculer toute les solutions de cette équation différentielle pour t > 0 (ne pas oublier les solutions constantes). La population
va-t-elle survivre ?
Correction
1. Il s’agit d’une EDO à variables séparables. On cherche d’abord les solutions constantes, i.e. des fonctions y(t ) = A
pour tout t ∈ R :
0 = A ln(A) ⇐⇒ A = 1.
On trouve ainsi une solution constante :
y(t ) ≡ 1.
Si on suppose que y(t ) 6= 1, l’EDO se réécrit comme
y 0 (t )
= −1;
y(t ) ln(y(t ))
2. Si y(0) > 1 alors y 0 (t ) < 0 (la population décroît) ; si 0 < y(0) < 1 alors y 0 (t ) > 0 (la population croît) ; comme y(t ) = 1
est solution et comme deux solutions ne peuvent pas se croiser, sans faire de calcul on voit que lorsque t tend vers
l’infini, la population tend vers la valeur d’équilibre y(t ) = 1 quelque soit le nombre d’individus à l’instant initial.
1
dx = 1z dz = ln |z| + c = ln | ln(x)| +C
R R
3. x ln(x)
70 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
Correction
On a a(x) = 1, b(x) = 3x 2 + 1 et g (x) = x 2 e −x , donc pour x ∈ R on a
R 2
? A(x) = 3x 1+1 dx = x 3 + x,
R 2 −x 3 3 x3
? K (x) = x e1 e A(x) dx = x 2 e x dx = 31 3x 2 e x dx = 13 e u(x) u 0 (x) dx = 13 e u(x) = e 3 .
R R R
µ ¶
x3 3
Toutes les solutions de l’EDO sont donc les fonctions y(x) = C + e 3 e −x −x pour C ∈ R.
On cherche parmi ces solutions celle qui vérifie y(0) = 1 ; comme y(0) = C + 31 , l’unique solution du problème de C AUCHY
µ ¶
x3 3
donné est la fonction y(x) = 23 + e 3 e −x −x .
Exercice 4.13
Résoudre le problème de C AUCHY
(
y 0 (x) + (3x 2 − 1)y(x) = x 2 e x ,
y(0) = −1.
Correction
On a a(x) = 1, b(x) = 3x 2 − 1 et g (x) = x 2 e x , donc pour x ∈ R on a
R 2
? A(x) = 3x 1−1 dx = x 3 − x,
R 2 x 3 x3
? K (x) = x 1e e A(x) dx = x 2 e x dx = e 3 .
R
µ ¶
x3 3
Toutes les solutions de l’EDO sont donc les fonctions y(x) = C + e 3 e −x +x pour C ∈ R.
On cherche parmi ces solutions celle qui vérifie y(0) = −1 ; comme y(0) = C + 31 , l’unique solution du problème de C AUCHY
µ ¶
x3 3
donné est la fonction y(x) = − 43 + e 3 e −x +x .
Exercice 4.14
Résoudre le problème de C AUCHY
(x−2)2
(
1
y 0 (x) + x−1 y(x) = x−1 ,
y(0) = 1.
Correction
2
On a a(x) = 1, b(x) = x−11
et g (x) = (x−2)
x−1 . b est défini pour x 6= 1 et comme on cherche un solution qui passe par le point
(0, 1), nous allons chercher une solution que pour x < 1. On a
R 1
? A(x) = x−1 dx = ln(1 − x),
R (x−2)2 A(x) 3
? K (x) = dx = − (x − 2)2 dx = (x−2)
R
x−1 e 3 .
3
³ ´
Toutes les solutions de l’EDO pour x < 1 s’écrivent y(x) = C + (x−2) 3
1
x−1 pour C ∈ R. On cherche parmi ces solutions
celle qui vérifie y(0) = 1 ; comme y(0) = −C + 83 , l’unique solution du problème de C AUCHY donné est la fonction y(x) =
(x−2)3
³ ´
5 1
3 + 3 x−1 .
Exercice 4.15
Résoudre le problème de C AUCHY
(
y 0 (x) + (4x 3 + 5)y(x) = x 3 e −5x ,
y(0) = 1.
Correction
On a a(x) = 1, b(x) = 4x 3 + 5 et g (x) = x 3 e −5x . On a
? A(x) = 4x 3 + 5 dx = x 4 + 5x,
R
4 x4
? K (x) = x 3 e −5x e A(x) dx = − x 3 e x dx = e 4 .
R R
© 2019-2020 G. Faccanoni 71
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
µ ¶
−x 4 4
Toutes les solutions de l’EDO sont donc les fonctions y(x) = C − e 4 e −x −5x pour C ∈ R.
On cherche parmi ces solutions celle qui vérifie y(0) = 1 ; comme y(0) = C + 14 , l’unique solution du problème de C AUCHY
µ 4
¶
3 ex 4
donné est la fonction y(x) = 4 + 4 e −x −5x .
Exercice 4.16
Établir s’il existe des solutions de y 0 (x) = −2y(x) + e −2x qui ont dérivée nulle en x = 0.
Correction
On a a(x) = 1, b(x) = 2 et g (x) = e −2x . On a
? A(x) = 2 dx = 2x,
R
Toutes les solutions de l’EDO sont donc les fonctions y(x) = (C + x)e −2x pour C ∈ R.
On cherche si parmi ces solutions il en existe qui vérifient y 0 (0) = 0 ; comme y 0 (x) = (1 − 2C − 2x)e −2x et y 0 (0) = 1 − 2C ,
l’unique solution de l’EDO qui a dérivée nulle en x = 0 est la fonction y(x) = ( 12 + x)e −2x .
Exercice 4.17
Établir s’il existe des solutions de y 0 (x) = −2x y(x) + x.
Correction
On a a(x) = 1, b(x) = 2x et g (x) = x. La solution de cette EDO est du type y(x) = y H (x) + y P (x) où y H (x) est la famille de
solutions de l’EDO homogène y 0 (x) = −2x y(x) et y P (x) est une solution particulière de l’EDO complète y 0 (x) = −2x y(x) + x.
On a y H (x) = C e −A(x) et, par exemple, on cherche y P sous la forme y P (x) = K (x)e −A(x) avec
? A(x) = b(x) 2
R R
a(x) dx = 2x dx = x ,
R g (x) 2 2
? K (x) = a(x) e A(x) dx = xe A(x) dx = xe x dx = 12 e x ,
R R
2
donc toutes les solutions de l’EDO sont les fonctions y(x) = C e −x + 21 pour C ∈ R.
Notons qu’il n’est même pas nécessaire de calculer K (x) ; en effet, il suffit de trouver une solution particulière évidente, par
exemple une solution constante. Si y(x) = B pour tout x est une solution de l’EDO complète, alors 0 = −2xB + x, i.e. B = 1/2.
2
On pose alors y P (x) = 1/2 et on a y(x) = y H (x) + y P (x) = C e −x + 21 .
2
Toutes les solutions de l’EDO sont donc les fonctions y(x) = C e −x + 12 pour C ∈ R.
Exercice 4.18
Résoudre l’équation différentielle
(x + 1)y 0 (x) + y(x) = (x + 1) sin(x)
sur des intervalles à préciser.
Correction
L’équation différentielle est linéaire du premier ordre. On la résout sur un intervalle où le coefficient de y 0 (x) n’est pas nul,
soit sur I 1 =] − ∞; −1[ ou sur I 2 =] − 1; +∞[. Sur chaque intervalle I 1 ou I 2 , l’équation s’écrit
£ ¤0
(x + 1)y(x) = (x + 1) sin(x).
sin(x) +C
y(x) = − cos(x) + avec C ∈ R.
(x + 1)
72 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
Exercice 4.19
Dans un circuit électrique de type résistance-inductance, le courant I évolue avec le temps selon
R V
I 0 (t ) + I (t ) =
L L
où R, L et V sont des constantes associées aux composantes électriques. Résolvez l’équation différentielle. La solution I
tend-elle vers une limite finie ?
Correction
? Résolution de l’équation homogène associée : c’est une équation à variables séparables. Ses solutions sont du type
R
I H (x) = αe − L t
avec α constante arbitraire. Elles comportent la fonction nulle et des fonctions qui ne s’annulent jamais.
R
? Recherche d’une solution particulière. Soit I 1 (t ) = e − L t une solution (non nulle) de l’EDO homogène et introduisons
R
la fonction auxiliaire inconnue K (t ) telle que I (t ) = K (t )I 1 (t ) = K (t )e − L t soit solution de notre EDO ; alors I 0 (t ) =
R R R
K 0 (t )e − L t − RL K (t )e − L t = (K 0 (t ) − RL K (t ))e − L t et si on reporte I 0 (t ) et I (t ) dans notre EDO on obtient
R V
K 0 (t )e − L t = .
L
R R R R
Ainsi K (t ) = VL e L t dt = VR e L t et enfin I P (t ) = VR e L t e − L t = VR .
R
Exercice 4.20
On considère un réservoir de capacité 5000 l rempli d’une solution
Taux : ? kg · min−1
Concentration : 0.03 kg · l−1
Débit : 25 l · min−1
Correction
1. Le taux auquel le sel entre est (0.03 kg)(25 l · min−1 )=0.75 kg · min−1 . Comme le réservoir contient constamment 5000 l
de liquide, la concentration est égale à y(t )/5000 (exprimée en kg · l−1 ). Le débit du mélange qui sort est alors de
y(t ) y(t )
25 l · min−1 , donc le taux auquel le sel sort est ( 5000 kg · l−1 )(25 l · min−1 )= 200 kg · min−1 . L’équation différentielle qui
décrit cette variation s’écrit alors
y(t )
y 0 (t ) = 0.75 −
200
2. On cherche d’abord les solutions constantes, i.e. des fonctions y(t ) = C pour tout t ∈ R :
C
0 = 0.75 − ⇐⇒ C = 150.
200
On trouve ainsi l’unique solution constante y(t ) = 150 pour tout t ∈ R. On sait que toute autre solution ne s’annulera
© 2019-2020 G. Faccanoni 73
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
3. Reste à calculer la quantité de sel après 30 minutes : y(30) = 150 − 130e −3/20 ' 38.1 kg.
Exercice 4.21
L’air d’un garage de 3 m × 5 m × 2 m est initialement chargée de 0.001% de monoxyde de carbone (CO). À l’instant t = 0, on
fait tourner un moteur et des fumées toxiques contenant 5% de CO se dégagent de la pièce à raison de 3 litres par minute.
Heureusement, l’air de la pièce est éliminée à la même vitesse de 3 l · min−1 . On note v(t ) le volume de CO présent dans la
pièce au temps t .
1. En supposant que le mélange se fait instantanément, montrer que cette situation est décrite par le problème de
C AUCHY
v(t )
(
v 0 (t ) = 0.15 − 10000 ,
v(0) = 0.3.
2. Déterminer le volume v(t ) de CO présent dans la pièce au temps t . Calculer vers quelle valeur limite v(t ) tend lorsque
t tend vers l’infini.
3. Le seuil critique pour la santé est de 0.015% de CO. Après combien de temps ce taux est-il atteint ?
Débit : 3 l · min−1
Concentration : 5%
Taux : ? l · min−1
Correction
1. Le taux de CO produit par minute est 0.05 × 3 l · min−1 = 0.15 l · min−1 . Le débit de l’air qui sort est de 3 l · min−1 , donc
v(t ) v(t )
le taux auquel le CO sort est 30 000 × 3 l · min−1 = −1
10 000 l · min . L’équation différentielle qui décrit cette variation
s’écrit alors
v(t )
v 0 (t ) = 0.15 − .
10 000
À l’instant t = 0 le volume de CO présent dans le garage est 0.001% × 30 000 l = 0.3 l.
2. On cherche d’abord les solutions constantes, i.e. des fonctions v(t ) = C pour tout t ∈ R :
C
0 = 0.15 − ⇐⇒ C = 1 500.
10 000
On trouve ainsi l’unique solution constante v(t ) = 1 500 pour tout t ∈ R. On sait que toute autre solution ne s’annulera
pas. L’EDO se réécrit alors comme
v 0 (t ) 1
= ;
1 500 − v(t ) 10 000
74 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
Exercice 4.22
Trouver la solution des problèmes de C AUCHY suivants :
Correction
i πh
b : 0; →R
2
1
t 7→
tan(t )
y 0 (t ) 1 1 cos(t )
Z Z
=− =⇒ dy = − dt =⇒ ln |y| = ln(sin(t )) +C avec C ∈ R.
y(t ) tan(t ) y sin(t )
© 2019-2020 G. Faccanoni 75
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
Nous avons vu que la solution générale de l’équation homogène associée est l’application
i πh
y H : 0; →R
2
κ
t 7→ avec κ ∈ R.
sin(t )
y P (t ) sin(t ) i πh
y P0 (t ) + = cos(t ) + = 2 cos(t ), ∀t ∈ 0; .
tan(t ) tan(t ) 2
Si on n’a pas trouvé directement une solution particulière, on peut utiliser ¯ la méthode ¤ de £variation de la constante :
π
on cherche une solution particulière sous la forme y P (t ) = K (t ) y H (t ) κ=1 où K : 0; 2 → R est une application
¯
inconnue dérivable sur 0; π2 . On a, pour tout t ∈ 0; π2 ,
¤ £ ¤ £
K (t )
y P (t ) =
sin(t )
K 0 (t ) sin(t ) − K (t ) cos(t )
y P0 (t ) =
sin2 (t )
Puisque y P est supposée être une solution particulière de l’EDO, on doit avoir pour tout t ∈ 0; π2 ,
¤ £
y P (t )
y P0 (t ) + = 2 cos(t )
tan(t )
K 0 (t ) sin(t ) − K (t ) cos(t ) K (t )
2
+ = 2 cos(t ).
sin (t ) tan(t ) sin(t )
K 0 (t ) = 2 cos(t ) sin(t )
ce qui impose K (t ) = sin2 (t ) + C , C ∈ R. Puisqu’on est intéressé par une seule solution particulière, on choisit
généralement C = O. On a alors i πh
y P : 0; →R
2
t 7→ sin(t )
On en déduit que la solution générale de l’EDO complète est l’application
i πh
y : 0; →R
2
κ
t 7→ y H (t ) + y P (t ) = + sin(t ) avec κ ∈ R.
sin(t )
Exercice 4.23
Déterminer la solution générale de l’EDO 3y 0 (t ) + 12y(t ) = 4 après avoir indiqué sur quelle intervalle la solution est définie.
76 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
Correction
L’équation 3y 0 (t ) + 12y(t ) = 4 admet pour solution générale l’application
y: R→R
1
t 7→ + κe −4t avec κ ∈ R.
3
En effet :
1. Calcul de la solution générale de l’équation homogène 3y 0 (t ) + 12y(t ) = 0.
Formellement on peut écrire
y 0 (t ) 1
Z Z
= −4 =⇒ dy = − 4 dt =⇒ ln |y| = −4t +C avec C ∈ R.
y(t ) y
yH : R → R
t 7→ κe −4t avec κ ∈ R.
Exercice 4.24
Déterminer la solution générale de l’EDO y 0 (t ) + y(t ) = e 3t après avoir indiqué sur quelle intervalle la solution est définie.
Correction
L’équation y 0 (t ) + y(t ) = e 3t admet pour solution générale l’application
y: R→R
1
t 7→ e 3t + κe −t avec κ ∈ R.
4
En effet :
1. Calcul de la solution générale de l’équation homogène y 0 (t ) + y(t ) = 0.
Toute solution non nulle est de la forme
yH : R → R
t 7→ κe −t avec κ ∈ R.
Exercice 4.25
1
Déterminer la solution générale de l’EDO y 0 (t ) + tan(t )y(t ) = cos(t ) après avoir indiqué sur quelle intervalle la solution est
définie.
Correction
1
L’équation y 0 (t ) + tan(t )y(t ) = cos(t ) admet pour solution générale l’application
i π π h
y : − + κπ; − + κπ → R
2 2
t 7→ sin(t ) +C κ cos(t ) avec C κ ∈ R.
© 2019-2020 G. Faccanoni 77
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
C κ désignant une constante réelle différente sur chacun des intervalles. En effet :
1. L’EDO est déjà écrite sous forme normalisée. Considérons les applications suivantes :
b: R→R g: R→R
1
t 7→ tan(t ) t 7→ cos(t )
π
+ κπ, κ ∈ Z et sont continues sur chacun des intervalles − π2 + κπ; π2 + κπ .
¤ £
Elles ne sont définies que pour t 6= 2
2. Calcul de la solution générale de l’équation homogène y 0 (t ) + tan(t )y(t ) = 0.
Formellement on peut écrire
y 0 (t ) sin(t ) 1 − sin(t )
Z Z
=− =⇒ dy = dt =⇒ ln |y| = ln | cos(t )| +C avec C ∈ R.
y(t ) cos(t ) y cos(t )
1
3. Calcul d’une solution particulière de l’EDO y 0 (t ) + tan(t )y(t ) = cos(t ) .
On vérifie aisément que la fonction
yP : R → R
t 7→ sin(t )
est une solution particulière de l’EDO. Si cette solution ne vous parait pas évidente, il est bien entendu possible de
déterminer une solution particulière en utilisant la méthode de la variation de la constante.
Exercice 4.26 t +2 et
Déterminer la solution générale de l’EDO y 0 (t ) + y(t ) = 2 après avoir indiqué sur quelle intervalle la solution est
t t
définie.
Correction
et
L’équation y 0 (t ) + t +2
t y(t ) = t2
admet pour solution générale l’application
y : R∗ → R
t
e e −t
2 + κ1 2 si t < 0 avec κ1 ∈ R,
t 7→ 2tt t
e e −t
+ κ2 2 avec κ2 ∈ R.
si t > 0
2t 2 t
κ1,2 désignant une constante réelle différente sur chacun des deux intervalles. En effet :
1. L’EDO est déjà écrite sous forme normalisée. Considérons les applications suivantes :
b: R→R g: R→R
t +2 et
t 7→ t 7→
t t2
Elles ne sont définies que pour t 6= 0 et sont continues sur chacun des intervalles ]−∞; 0[ et ]0; +∞[.
2. Calcul de la solution générale de l’équation homogène y 0 (t ) + t +2
t y(t ) = 0.
Formellement on peut écrire
y 0 (t ) 2 1 2
Z Z
= −1 − =⇒ dy = −1 − dt =⇒ ln |y| = −t − 2 ln |t | +C avec C ∈ R.
y(t ) t y t
y H : R∗ → R
e −t
κ1 2 si t < 0 avec κ1 ∈ R,
t 7→ t
e −t
κ2 avec κ2 ∈ R.
si t > 0
t2
78 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
t
3. Calcul d’une solution particulière de l’EDO y 0 (t ) + t +2 e
t y(t ) = t 2 .
On vérifie aisément que la fonction
y P : R∗ → R
e −t
t 7→
2t 2
est une solution particulière de l’EDO. Si cette solution ne vous parait pas évidente, il est bien entendu possible de
déterminer une solution particulière en utilisant la méthode de la variation de la constante.
Correction
1−e −2t
L’équation y 0 (t ) + y(t ) = e t +e −t
admet pour solution générale l’application
y: R→R
t 7→ e −t ln(e t + e −t ) + κe −t avec κ ∈ R.
En effet :
1. Calcul de la solution générale de l’équation homogène y 0 (t ) + y(t ) = 0.
Toute solution non nulle est de la forme
yH : R → R
t 7→ κe −t avec κ ∈ R.
1−e −2t
2. Calcul d’une solution particulière de l’EDO y 0 (t ) + y(t ) = e t +e −t
.
L’application
g: R→R
1 − e −2t
t 7→ = e −t tanh(t )
e t + e −t
est continue sur R. En utilisant la méthode de variation de la constante, on cherche une solution particulière sous la
forme y P (t ) = K (t ) y H (t )¯κ=1 où K : R → R est une application inconnue dérivable sur R. On a, pour tout t ∈ R,
¯
y P (t ) = K (t )e −t
y P0 (t ) = K 0 (t )e −t − K (t )e −t
Puisque y P est supposée être une solution particulière de l’EDO, on doit avoir pour tout t ∈ R,
1 − e −2t
y P0 (t ) + y P (t ) =
e t + e −t
autrement dit la fonction K doit vérifier
e t − e −t sinh(t )
K 0 (t ) = = tanh(t ) =
e t + e −t cosh(t )
ce qui impose K (t ) = ln(cosh(t )) +C = ln(e t + e −t ) − ln(2) +C , C ∈ R. Puisqu’on est intéressé par une seule solution
particulière, on peut choisir C = ln(2). On a alors
yP : R → R
t 7→ e −t ln(e t + e −t )
Exercice 4.28
Déterminer la solution générale de l’EDO t y 0 (t ) + (3t + 1)y(t ) = e −3t après avoir indiqué sur quelle intervalle la solution est
définie.
Correction
e −3t
Étudions d’abord l’EDO normalisée y 0 (t ) + 3tt+1 y(t ) = t qui doit être considérée sur ] − ∞; 0[ et sur ]0; +∞[.
© 2019-2020 G. Faccanoni 79
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
A 1 : ] − ∞; 0[ → R
3t + 1 1
t 7→ = 3+
t t
est continue et strictement positive sur ] − ∞; 0[. Formellement on peut alors écrire
y 0 (t ) 1 1 1
Z Z
= −3 + =⇒ dy = −3 + dt =⇒ ln |y| = −3t + ln(−t ) +C avec C ∈ R.
y(t ) t y t
y H : ] − ∞; 0[ → R
κ1 −3t
t 7→ e avec κ1 ∈ R.
t
−3t
? Calcul d’une solution particulière de l’EDO y 0 (t ) + 3tt+1 y(t ) = e t pour t ∈] − ∞; 0[.
On vérifie aisément que
y P : ] − ∞; 0[ → R
t 7→ e −3t
est une solution particulière de l’EDO puisque
µ ¶
3t + 1 1 1
y P0 (t ) + y P (t ) = −3e −3t + 3 + e −3t = e −3t , ∀t ∈] − ∞; 0[.
t t t
Si cette solution ne vous parait pas évidente, il est bien entendu possible de déterminer une solution particulière
en utilisant la méthode de la variation de la constante.
? Conclusion sur ] − ∞; 0[.
On en déduit que la solution générale de l’EDO sur ] − ∞; 0[ est l’application
y : ] − ∞; 0[ → R
³κ ´
1
t 7→ y H (t ) + y P (t ) = + 1 e −3t avec κ1 ∈ R.
t
? Calcul de la solution générale de l’équation homogène y 0 (t ) + 3tt+1 y(t ) = 0 pour t ∈]0; +∞[.
L’application
A 2 : ]0; +∞[ → R
3t + 1 1
t 7→ = 3+
t t
est continue et strictement positive sur ]0; +∞[. Formellement on peut alors écrire
y 0 (t ) 1 1 1
Z Z
= −3 + =⇒ dy = −3 + dt =⇒ ln |y| = −3t + ln(t ) +C avec C ∈ R.
y(t ) t y t
y H : ]0; +∞[ → R
κ2 −3t
t 7→ e avec κ2 ∈ R.
t
−3t
? Calcul d’une solution particulière de l’EDO y 0 (t ) + 3tt+1 y(t ) = e t pour t ∈]0; +∞[.
On vérifie aisément que
y P : ]0; +∞[ → R
t 7→ e −3t
80 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
Si cette solution ne vous parait pas évidente, il est bien entendu possible de déterminer une solution particulière
en utilisant la méthode de la variation de la constante.
? Conclusion sur ]0; +∞[.
On en déduit que la solution générale de l’EDO sur ]0; +∞[ est l’application
y : ]0; +∞[ → R
³κ ´
2
t 7→ y H (t ) + y P (t ) = + 1 e −3t avec κ2 ∈ R.
t
3. Étude sur R.
La limite de la solution générale pour t → 0− n’est bornée que si κ1 = 0. De même, la limite de la solution générale
pour t → 0+ n’est bornée que si κ2 = 0. Par conséquente, si une solution générale de l’EDO existe pour tout t ∈ R, ce
ne peut être que l’application
ϕ: R → R
t 7→ e −3t .
Cette application est continue et dérivable sur R et
donc ϕ est bien solution de l’EDO sur R. C’est donc l’unique solution de l’EDO sur R.
Exercice 4.29
Déterminer la solution générale de l’EDO 2t (t + 1)y 0 (t ) + (t + 1)y(t ) = 1 après avoir indiqué sur quelle intervalle la solution
est définie.
Correction
1 1
Étudions d’abord l’EDO normalisée y 0 (t ) + 2t y(t ) = 2t (t +1) qui doit être considérée sur ] − ∞; −1[, sur ] − 1; 0[ et sur ]0; +∞[.
1. Étude sur ] − ∞; −1[.
1
? Calcul de la solution générale de l’équation homogène y 0 (t ) + 2t y(t ) = 0 pour t ∈] − ∞; −1[.
L’application
b 1 : ] − ∞; −1[ → R
1
t 7→
2t
est continue et strictement négative sur ] − ∞; −1[. Formellement on peut alors écrire
y 0 (t ) 1 1 1 p
Z Z
=− =⇒ dy = − dt =⇒ ln |y| = − ln( −t ) +C avec C ∈ R.
y(t ) 2t y 2t
y H : ] − ∞; −1[ → R
κ1
t 7→ p avec κ1 ∈ R.
−t
1
? Calcul d’une solution particulière de l’EDO y 0 (t ) + 2t y(t ) = 2t (t1+1) pour t ∈] − ∞; −1[.
L’EDO ne semble pas posséder de solution évidente. On a donc recours à la méthode de la variation de ¯ la
constante pour en déterminer une. On cherche une solution particulière sous la forme y P (t ) = K (t ) y H (t )¯κ1 =1
où K : ] − ∞; −1[→ R est une application inconnue dérivable sur ] − ∞; −1[. On a, pour tout t ∈] − ∞; −1[,
K (t )
y P (t ) = p = (−t )−1/2 K (t )
−t
K 0 (t ) K (t )
y P0 (t ) = p − p
−t 2 (−t )3
© 2019-2020 G. Faccanoni 81
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
Puisque y P est supposée être une solution particulière de l’EDO, on doit avoir pour tout t ∈] − ∞; −1[,
y P (t ) 1
y P0 (t ) + =
2t 2t (t + 1)
K 0 (t ) K (t ) K (t ) 1
p − p + p = .
−t 2 (−t )3 2t −t 2t (t + 1)
1
K 0 (t ) = − p
2 −t (t + 1)
³p ´
ce qui impose K (t ) = 12 ln p−t +1 + C , C ∈ R. Puisqu’on est intéressé par une seule solution particulière, on
−t −1
choisit généralement C = O. On a alors
y P : ] − ∞; −1[ → R
µp ¶
1 −t + 1
t 7→ p ln p
2 −t −t − 1
y : ] − ∞; −1[ → R
µp
κ1
¶
1 −t + 1
t 7→ y H (t ) + y P (t ) = p + p ln p avec κ1 ∈ R.
−t 2 −t −t − 1
b 2 : ] − 1; 0[ → R
1
t 7→
2t
est continue et strictement négative sur ] − 1; 0[. Formellement on peut alors écrire
y 0 (t ) 1 1 1 p
Z Z
=− =⇒ dy = − dt =⇒ ln |y| = − ln( −t ) +C avec C ∈ R.
y(t ) 2t y 2t
y H : ] − 1; 0[ → R
κ2
t 7→ p avec κ2 ∈ R.
−t
1
? Calcul d’une solution particulière de l’EDO y 0 (t ) + 2t y(t ) = 2t (t1+1) pour t ∈] − 1; 0[.
On cherche une solution particulière sous la forme y P (t ) = K (t ) y H (t )¯κ2 =1 où K : ]−1; 0[→ R est une application
¯
K (t )
y P (t ) = p = (−t )−1/2 K (t )
−t
K 0 (t ) K (t )
y P0 (t ) = p − p
−t 2 (−t )3
Puisque y P est supposée être une solution particulière de l’EDO, on doit avoir pour tout t ∈] − 1; 0[,
y P (t ) 1
y P0 (t ) + =
2t 2t (t + 1)
82 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
K 0 (t ) K (t ) K (t ) 1
p − p + p = .
−t 2 (−t ) 3 2t −t 2t (t + 1)
1
K 0 (t ) = − p
2 −t (t + 1)
³p ´
ce qui impose K (t ) = 12 ln p−t +1 + C , C ∈ R. Puisqu’on est intéressé par une seule solution particulière, on
−t −1
choisit généralement C = O. On a alors
y P : ] − 1; 0[ → R
µ p ¶
1 1 + −t
t 7→ p ln p
2 −t 1 − −t
y : ] − 1; 0[ → R
p ¶
κ2
µ
1 1 + −t
t 7→ y H (t ) + y P (t ) = p + p ln p avec κ2 ∈ R.
−t 2 −t 1 − −t
b 3 : ]0; +∞[ → R
1
t 7→
2t
est continue et strictement positive sur ]0; +∞[. Formellement on peut alors écrire
y 0 (t ) 1 1 1 p
Z Z
=− =⇒ dy = − dt =⇒ ln |y| = − ln( t ) +C avec C ∈ R.
y(t ) 2t y 2t
y H : ]0; +∞[ → R
κ3
t 7→ p avec κ3 ∈ R.
t
1
? Calcul d’une solution particulière de l’EDO y 0 (t ) + 2t y(t ) = 2t (t1+1) pour t ∈]0; +∞[.
On cherche une solution particulière sous la forme y P (t ) = K (t ) y H (t )¯κ3 =1 où K : ]0; +∞[→ R est une application
¯
K (t )
y P (t ) = p = (t )−1/2 K (t )
t
0
K (t ) K (t )
y P0 (t ) = p − p
t 2 (t )3
Puisque y P est supposée être une solution particulière de l’EDO, on doit avoir pour tout t ∈]0; +∞[,
y P (t ) 1
y P0 (t ) + =
2t 2t (t + 1)
K 0 (t ) K (t ) K (t ) 1
p − p + p = .
t 2 t 3 2t t 2t (t + 1)
© 2019-2020 G. Faccanoni 83
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
y : ]0; +∞[ → R
p
κ3 arctan( t )
t 7→ y H (t ) + y P (t ) = p + p avec κ3 ∈ R.
t t
4. Étude sur R.
L’EDO donnée et l’EDO normalisée ne sont pas équivalentes puisqu’elles n’ont pas le même domaine de validité. Une
solution de l’EDO donnée, si elle existe, sera solution de l’EDO normalisée sur les trois intervalles ] − ∞; −1[, ] − 1; 0[ et
]0; +∞[. On a déterminé les solutions de l’EDO normalisée sur ces trois intervalles. Se pose maintenant la question
de savoir si à partir des solutions de l’EDO normalisée sur les trois intervalles ] − ∞; −1[, ] − 1; 0[ et ]0; +∞[ on peut
trouver (construire) une solution de l’EDO donnée sur R. Autrement dit : peut-on trouver une application ϕ définie
sur R, dérivable sur R telle que sa restriction sur chacun des trois intervalles est solution de l’EDO normalisée ? Si une
telle solution existe, elle est nécessairement de la forme
³p
κ
´
p 1 + p1 ln p−t +1 si t < −1,
−t 2 −t ³ −t
p−1 ´
κ2 1 1+ −t
ϕ(t ) = p + p ln p si − 1 < t < 0,
−t 2 −t p 1− −t
κ arctan( t )
3
p + p si t > 0,
t t
où κ1 , κ2 , κ3 désignent trois constantes réelles. Pour que cette fonction soit solution, il faut qu’elle soit dérivable sur
R : il nous faut donc regarder si elle est prolongeable par continuité en −1 et en 0 et si le prolongement ainsi défini est
dérivable en −1 et en 0.
4.1. Étude du raccord en −1.
La fonction ϕ ne peut pas être prolongée par continuité en −1 car elle n’est pas bornée en −1 puisque
µp
κ1
¶
1 −t + 1
lim − ϕ(t ) = lim − p + p ln p = +∞.
t →−1 t →−1 −t 2 −t −t − 1
y : ] − ∞; −1[ → R
µp
κ1
¶
1 −t + 1
t 7→ p + p ln p
−t 2 −t −t − 1
84 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
y : ] − 1; 0[ → R
p ¶
κ2
µ
1 1 + −t
t 7→ p + p ln p
−t 2 −t 1 − −t
y : ]0; +∞[ → R
p
κ3 arctan( t )
t 7→ p + p
t t
avec κ1 , κ2 , κ2 ∈ R.
Elle admet une unique solution sur l’intervalle ] − 1; +∞[ qui est l’application
y : ] − 1; +∞[ → R
³ p ´
p1 ln 1+p−t si − 1 < t < 0,
2 −t
1− −t
t 7→ 1 si t = 0,
p
arctan(
t )
p si t > 0.
t
Correction
(a) L’EDO y 0 (t ) − 1t y(t ) = (y(t ))3 sin(t ) est une équation différentielle de B ERNOULLI. Comme u(t ) = 1 pour tout t ∈ R∗ , on
cherche sa solution générale sur ] − ∞; 0[ et sur ]0; +∞[.
v(t ) 1
Z Z
? A(t ) = (1 − α) dt = dt = 2 ln |t |,
u(t ) t
w(t ) A(t ) sin(t ) cos(t )
Z Z
? K (t ) = (1 − α) e dt = − t 2 sin(t ) dt = 2 cos(t ) − 4 −4 2 ,
u(t ) t t
¢ −A(t ) ³ ´ C 1,2 +2 cos(t )−4 sin(t ) cos(t )
t −4 t 2
? z(t ) = C 1,2 + K (t ) e = C 1,2 + 2 cos(t ) − 4 sin(t ) cos(t ) −2 ln |t |
¡
t − 4 t 2 e = t 2 ,
? y(t ) = (z(t ))−1/2 = p1
z(t )
et on conclut que la solution générale de l’EDO de B ERNOULLI assignée est
y : R∗ → R
t
avec C 1 ∈ R+ ,
qC +2 cos(t )−4 sin(t ) −4 cos(t )
si t < 0
1,2 t t2
t 7→ −t
q
sin(t ) cos(t )
si t > 0 avec C 2 ∈ R+ ,
C 1,2 +2 cos(t )−4 −4
t t2
(b) L’EDO y 0 (t ) + t y(t ) = t 3 (y(t ))2 est une équation différentielle de B ERNOULLI. Comme u(t ) = 1 pour tout t ∈ R, on
cherche sa solution générale sur R.
? Solution nulle : la fonction y(t ) = 0 pour tout t ∈ R est solution de l’EDO donnée. Toute autre solution ne s’annule
jamais. Supposons dans la suite que y(t ) 6= 0 pour tout t ∈ R.
? Réduction à une EDO linéaire du premier ordre : si on pose z = 1/y, i.e. y = 1/z, elle se réécrit
z 0 (t ) − t z(t ) = −t 3 .
© 2019-2020 G. Faccanoni 85
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
zP : R → R
t 7→ t 2 + 2
z P0 (t ) − t z P (t ) = 2t − t (t 2 + 2) = 2t − t 3 + 2t = −t 3 pout tout t ∈ R.
Correction
1. Le polynôme caractéristique associé à l’EDO est λ2 − 3λ + 2 = 0 qui a discriminant ∆ = 1 : on a deux solutions réelles
distinctes
3−1 3+1
λ1 = =1 et λ2 = = 2.
2 2
Toutes les solutions sont alors les fonctions y(x) = C 1 e x +C 2 e 2x pour tout C 1 ,C 2 ∈ R.
2. Le polynôme caractéristique associé à l’EDO est λ2 − 4λ + 4 = 0 qui a discriminant ∆ = 0 : on a deux solutions réelles
coïncidentes
4
λ1 = λ2 = = 2.
2
Toutes les solutions sont les fonctions y(x) = (C 1 +C 2 x)e 2x pour tout C 1 ,C 2 ∈ R.
3. Le polynôme caractéristique associé à l’EDO est λ2 − 2λ + 2 = 0 qui a discriminant ∆ = −4 : on a deux solutions
complexes conjuguées
2 − 2i 2 + 2i
λ1 = et λ2 = .
2 2
Comme ℜ(λ2 ) = 1 et ℑ(λ2 ) = 1, toutes les solutions sont les fonctions y(x) = e x (C 1 cos(x) + C 2 sin(x)) pour tout
C 1 ,C 2 ∈ R.
4. Comme l’EDO n’est pas homogène on cherche ses solutions sous la forme y(x) = y H (x) + y P (x) où les y H sont toutes
les solutions de l’EDO homogène associée et y P est une solution particulière de l’EDO complète.
Le polynôme caractéristique associé à l’EDO est λ2 − 1 = 0 qui a discriminant ∆ = 1 : on a deux solutions réelles
distinctes
λ1 = −1 et λ2 = 1.
Toutes les solutions de l’homogène sont les fonctions y H (x) = C 1 e x +C 2 e −x pour tout C 1 ,C 2 ∈ R.
Comme le terme source est g (x) = e 2x , on a n = 0, µ = 2 et ϑ = 0 ; comme ∆ > 0, ϑ = 0 mais µ 6= λ1 et µ 6= λ2 alors
m = 0 : la solution particulière sera de la forme y P (x) = q 1 e 2x . Pour qu’elle soit une solution particulière on doit
imposer qu’elle vérifie l’EDO complète ; comme y P0 (x) = 2q 1 e 2x et y P00 (x) = 4q 1 e 2x il faut que
4q 1 e 2x − q 1 e 2x = e 2x
qui donne q 1 = 31 .
En conclusion toutes les solutions sont les fonctions y(x) = C 1 e x +C 2 e −x + 13 e 2x pour tout C 1 ,C 2 ∈ R.
86 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
5. Toute solution est de la forme y(x) = y H (x)+ y P (x) où y H représente toutes les solutions de l’équation y 00 (x)− y 0 (x) = 0
tandis que y P est une solution particulière de y 00 (x) − y 0 (x) = e x . Le polynôme caractéristique est 2λ2 − 1 = 0 qui
a comme racines les deux réels λ1 = 0 et λ2 = 1. Par conséquent y H (x) = C 1 + C 2 e x . En adoptant la notation du
cours on a n = 0, µ = 1, ϑ = 0, ∆ > 0 et µ = λ2 = 1 par conséquent m = 1 : la solution particulière est alors de la
forme y P (x) = C xe x . Pour déterminer C¢ on 0
¡ impose ࢠy P d’être solution de l’EDO. Comme y P (x) = C (1 + x)e et
x
y P00 (x) = C (2 + x)e x , il faut que C (2 + x)e x − C (1 + x)e x = e x , ce qui donne C = 1. On conclut que toutes les solutions
¡
de l’EDO s’écrivent
y(x) = C 1 +C 2 e x + xe x .
Exercice 4.32
Résoudre le problème de C AUCHY
00 0
y (x) − 2y (x) + 10y(x) = 0,
y(0) = 1,
0
y (0) = 2.
Correction
L’équation différentielle est linéaire du second ordre, à coefficients constants, et sans second membre (i.e. elle est déjà
homogène !). L’équation caractéristique λ2 − 2λ + 10 = 0 a discriminant ∆ < 0. Comme σ = 1 et ω = 3, l’intégrale générale de
l’équation homogène est
y H (x) = e x (c 1 cos(3x) + c 2 sin(3x)) , c 1 , c 2 ∈ R.
Puisque y(0) = 1 alors c 1 = 1. Comme y 0 (0) = 2 et y 0 (x) = e x ((c 1 + 3c 2 ) cos(3x) + (c 1 − 3c 2 ) sin(3x)) alors c 2 = 1/3. On conclut
que la solution du problème de C AUCHY est
µ ¶
1
y(x) = e x cos(3x) + sin(3x) .
3
Exercice 4.33
Calculer toutes les solutions de l’EDO 2y 00 (x) − 5y 0 (x) + 3y(x) = e x .
Correction
Toute solution est de la forme
y(x) = y H (x) + y P (x)
où y H représente toutes les solutions de l’équation 2y 00 (x) − 5y 0 (x) + 3y(x) = 0 tandis que y P est une solution particulière de
2y 00 (x) − 5y 0 (x) + 3y(x) = e x .
? Calcul de y H . Le polynôme caractéristique est 2λ2 − 5λ + 3 = 0 qui a comme racines les deux réels λ1 = 1 et λ = 3/2.
Par conséquent
y H (x) = C 1 e x +C 2 e 3x/2 .
2 C (2 + x)e x − 5 C (1 + x)e x + 3 C xe x = e x ,
¡ ¢ ¡ ¢ ¡ ¢
ce qui donne C = 1.
On conclut que toutes les solutions de l’EDO s’écrivent
y(x) = C 1 e x +C 2 e 3x/2 + xe x .
Exercice 4.34
Calculer toutes les solutions de l’EDO 2y 00 (x) − 7y 0 (x) + 5y(x) = −3e x .
Correction
Toute solution est de la forme
y(x) = y H (x) + y P (x)
où y H représente toutes les solutions de l’équation 2y 00 (x) − 7y 0 (x) + 5y(x) = 0 tandis que y P est une solution particulière de
2y 00 (x) − 7y 0 (x) + 5y(x) = −3e x .
© 2019-2020 G. Faccanoni 87
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
? Calcul de y H . Le polynôme caractéristique est 2λ2 − 7λ + 5 = 0 qui a comme racines les deux réels λ1 = 1 et λ = 5/2.
Par conséquent
y H (x) = C 1 e x +C 2 e 5x/2 .
ce qui donne C = 1.
On conclut que toutes les solutions de l’EDO s’écrivent
y(x) = C 1 e x +C 2 e 5x/2 + xe x .
Exercice 4.35
x
00 0
y (x) − 2y (x) + y(x) = (x + 1)e ,
y(0) = 1,
0
y (0) = 2.
Correction
Comme l’EDO n’est pas homogène on cherche ses solutions sous la forme y(x) = y H (x) + y P (x) où les y H sont toutes les
solutions de l’EDO homogène y 00 (x) − 2y 0 (x) + y(x) = 0 (y H contient deux constantes d’intégration) et y P est une solution
particulière de l’EDO complète y 00 (x) − 2y 0 (x) + y(x) = (x + 1)e x .
Le polynôme caractéristique associé à l’EDO est λ2 − 2λ + 1 = 0 qui a discriminant ∆ = 4 − 4 = 0, l’unique racine réelle est
λ1 = λ2 = 1.
Par conséquent, les solutions de l’EDO homogène sont les fonctions y H (x) = (C 1 +C 2 x)e x pour tout C 1 ,C 2 ∈ R.
Comme le terme source est g (x) = (x + 1)e x , on a n = 1, µ = 1 et ϑ = 0 ; puisque ∆ = 0, ϑ = 0 et µ = λ alors m = 2 : la
solution particulière sera de la forme y P (x) = x 2 e x (α + βx) = (αx 2 + βx 3 )e x . Pour qu’elle soit une solution particulière on
doit imposer qu’elle vérifie l’EDO complète ; comme y P0 (x) = (2αx + 3βx 2 + αx 2 + βx 3 )e x = (2αx + (3β + α)x 2 + βx 3 )e x et
y P00 (x) = (2α + (6β + 4α)x + (6β + α)x 2 + βx 3 )e x , il faut
(2α + (6β + 4α)x + (6β + α)x 2 + βx 3 )e x − 2(2αx + (3β + α)x 2 + βx 3 )e x + (αx 2 + βx 3 )e x = (x + 1)e x
Exercice 4.36
00 0 −x
y (x) + 2y (x) + y(x) = (1 − x)e ,
y(0) = 1,
0
y (0) = 2.
Correction
Comme l’EDO n’est pas homogène on cherche ses solutions sous la forme y(x) = y H (x) + y P (x) où les y H sont toutes les
solutions de l’EDO homogène y 00 (x) + 2y 0 (x) + y(x) = 0 (y H contient deux constantes d’intégration) et y P est une solution
particulière de l’EDO complète y 00 (x) + 2y 0 (x) + y(x) = (1 − x)e −x .
Le polynôme caractéristique associé à l’EDO est λ2 + 2λ + 1 = 0 qui a discriminant ∆ = 4 − 4 = 0, la solution réelle est
λ1 = λ2 = −1.
Par conséquent, les solutions de l’EDO homogène sont les fonctions y H (x) = (C 1 +C 2 x)e −x pour tout C 1 ,C 2 ∈ R.
88 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
Comme le terme source est g (x) = (1 − x)e −x , selon la notation du polycopié on a n = 1, µ = −1 et ϑ = 0 ; comme ∆ = 0,
ϑ = 0 et µ = λ alors m = 2 : la solution particulière sera de la forme y P (x) = x 2 e −x (α + βx) = (αx 2 + βx 3 )e −x . Pour qu’elle
soit une solution particulière on doit imposer qu’elle vérifie l’EDO complète ; comme y P0 (x) = (2αx + (3β − α)x 2 − βx 3 )e −x
et y P00 (x) = (2α + (6β − 4α)x − (6β − α)x 2 + βx 3 )e −x , il faut
(2α + (6β − 4α)x − (6β − α)x 2 + βx 3 )e −x + 2(2αx + (3β − α)x 2 − βx 3 )e −x + (αx 2 + βx 3 )e −x = (1 − x)e −x
c’est-à-dire (2α) + (6β − 4α + 4α)x + (−6β + α + 6β − 2α + α)x 2 + (+β − 2β + β)x 3 = 1 + x, ce qui donne α = 21 et β = 16 .
En conclusion les solutions de l’EDO sont les fonctions y(x) = (C 1 +C 2 x + 12 x 2 + 16 x 3 )e −x pour tout C 1 ,C 2 ∈ R.
On cherche parmi ces solutions celles qui vérifient y(0) = 1 ; comme y(0) = C 1 on obtient les fonctions y(x) = (1 + C 2 x +
1 2 1 3 −x
2 x + 6 x )e pour tout C 2 ∈ R. Parmi ces solutions, on cherche maintenant celle qui vérifie y 0 (0) = 2 ; comme y 0 (0) = 1 +C 2
on conclut que l’unique solution du problème de C AUCHY donné est la fonction y(x) = (1 + x + 12 x 2 + 16 x 3 )e −x .
Exercice 4.37
4x
00 0
y (x) − 5y (x) + 4y(x) = 2xe ,
y(0) = 1,
0
y (0) = 2.
Correction
Comme l’EDO n’est pas homogène on cherche ses solutions sous la forme y(x) = y H (x) + y P (x) où les y H sont toutes les
solutions de l’EDO homogène y 00 (x) − 5y 0 (x) + 4y(x) = 0 (y H contient deux constantes d’intégration) et y P est une solution
particulière de l’EDO complète y 00 (x) − 5y 0 (x) + 4y(x) = 2xe 4x .
Le polynôme caractéristique associé à l’EDO est λ2 − 5λ + 4 = 0 qui a discriminant ∆ = 25 − 16 = 9, les deux solutions réelles
sont λ1 = 1 et λ2 = 4. Par conséquent, les solutions de l’EDO homogène sont les fonctions y H (x) = C 1 e x + C 2 e 4x pour tout
C 1 ,C 2 ∈ R.
Comme le terme source est g (x) = 2xe 4x , selon la notation du polycopié on a n = 1, µ = 4 et ϑ = 0 ; comme ∆ > 0, ϑ = 0
et µ = λ2 alors m = 1 : la solution particulière sera de la forme y P (x) = xe 4x (αx + β) = (αx 2 + βx)e 4x . Pour qu’elle soit
une solution particulière on doit imposer qu’elle vérifie l’EDO complète ; comme y P0 (x) = (β + (2α + 4β)x + 4αx 2 )e 4x et
y P00 (x) = ((2α + 8β) + (16α + 16β)x + 16αx 2 )e 4x , il faut
((2α + 8β) + (16α + 16β)x + 16αx 2 )e 4x − 5(β + (2α + 4β)x + 4αx 2 )e 4x + 4(αx 2 + βx)e 4x = 2xe 4x
c’est-à-dire (2α + 8β − 5β) + (16α + 16β − 10α − 20β + β)x + (16α − 20α + 4α)x 2 = 2x, ce qui donne α = 31 et β = − 29 .
En conclusion les solutions de l’EDO sont les fonctions y(x) = C 1 e x + C 2 − 29 x + 13 x 2 e 4x pour tout C 1 ,C 2 ∈ R. On cherche
¡ ¢
parmi ces solutions celles qui vérifient y(0) = 1 ; comme y(0) = C 1 +C 2 on obtient les fonctions y(x) = C 1 e x + 1 −C 1 − 29 x + 13 x 2 e 4x
¡ ¢
pour tout C 1 ∈ R. Parmi ces solutions, on cherche maintenant celle qui vérifie y 0 (0) = 2 ; comme y 0 (0) = 4 − 29 − 3C 1 on
conclut que l’unique solution du problème de C AUCHY donné est la fonction y(x) = 16 x
¡ 11 2 1 2 4x
¢
27 e + 27 − 9 x + 3 x e .
Exercice 4.38
Résoudre le problème de C AUCHY
00
y (x) + 9y(x) = 108x cos(3x),
y(0) = 0,
0
y (0) = 0.
Correction
Comme l’EDO n’est pas homogène on cherche ses solutions sous la forme y(x) = y H (x) + y P (x) où les y H sont toutes
les solutions de l’EDO homogène y 00 (x) + 9y(x) = 0 (y h contient deux constantes d’intégration) et y P est une solution
particulière de l’EDO complète y 00 (x) + 9y(x) = 108x cos(3x).
Le polynôme caractéristique associé à l’EDO est λ2 +9 = 0 qui a discriminant ∆ = −9, les deux racines complexes conjuguées
sont
λ1 = −3i , λ2 = 3i .
Par conséquent, les solutions de l’EDO homogène sont les fonctions y H (x) = C 1 cos(3x) +C 2 sin(3x) pour tout C 1 ,C 2 ∈ R.
Comme le terme source est g (x) = 108x cos(3x), ¡ on a n = 1, µ = 0 et ϑ = 3 ; puisque
¢ ∆ < 0, ϑ = ω = 3 et µ = 0 alors m = 1 : la
solution particulière sera de la forme y P (x) = x (αx + β) cos(3x) + (γx + δ) sin(3x) = (αx 2 + βx) cos(3x) + (γx 2 + δx) sin(3x).
Pour qu’elle soit une solution particulière on doit imposer qu’elle vérifie l’EDO complète ; comme y P0 (x) = (3γx 2 + (2α +
© 2019-2020 G. Faccanoni 89
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
3δ)x + β) cos(3x) + (−3αx 2 + (2γ − 3β)x + δ) sin(3x) et y P00 (x) = (−9αx 2 + (12γ − 9β)x + 2α + 6δ) cos(3x) + (−9γx 2 + (−12α −
9δ)x + 2γ − 6β) sin(3x), il faut
(−9αx 2 + (12γ − 9β)x + 2α + 6δ) cos(3x) + (−9γx 2 + (−12α − 9δ)x + 2γ − 6β) sin(3x)
+ 9 (αx 2 + βx) cos(3x) + (γx 2 + δx) sin(3x)
¡ ¢
= 108x cos(3x)
c’est-à-dire
ce qui donne
12γ = 108,
2α + 6δ = 0,
=⇒ α = δ = 0, β = 3, γ = 9.
−12α = 0,
2γ − 6β = 0
En conclusion les solutions de l’EDO sont les fonctions y(x) = (C 1 + 3x) cos(3x) + (C 2 + 9x 2 ) sin(3x) pour tout C 1 ,C 2 ∈ R.
On
¡ cherche ¢ parmi ces solutions celles qui vérifient y(0) = 0 ; comme y(0) = C 1 on obtient les fonctions y(x) = 3x cos(3x) +
C 2 + 9x 2 sin(3x) pour tout C 2 ∈ R. Parmi ces solutions, on cherche maintenant celle qui vérifie y 0 (0) = 0 ; comme y 0 (0) =
3 + 3C 2 on conclut que l’unique solution du problème de C AUCHY donné est la fonction y(x) = 3x cos(3x) + (9x 2 − 1) sin(3x).
Exercice 4.39
Trouver la solution des problèmes de C AUCHY suivants :
00 p 0 00 0
00 0
y (t ) − 2 2y (t ) + 2y(t ) = 0
y (t ) − 4y (t ) + 5y(t ) = 0
y (t ) + y (t ) − 12y(t ) = 0
(1) y(0) = 2 (2) y(π/2) = 1 (3) y(0) = 1
0
0
0
y (0) = 3 y (π/2) = 1 y (0) = 4
Correction
p p
(1) L’équation caractéristique λ2 − 2 2λ + 2 = 0 a discriminant ∆ = 0 et l’unique racine double λ = 2. La solution générale
de l’EDO sur R est donc
p
2t
y(t ) = (At + B )e , A, B ∈ R.
Déterminons la valeur des constantes A et B correspondant à l’unique solution vérifiant y(0) = 2 et y 0 (0) = 3. On a
p p
y 0 (t ) = [A + (At + B ) 2]e 2t ,
donc
( ( (
y(0) = 2, B = 2, B = 2,
⇐⇒ p ⇐⇒ p
y 0 (0) = 3 A + 2B = 3 A = 3 − 2 2.
y: R→R
p p
t 7→ [(3 − 2 2)t + 2]e 2t
(2) L’équation caractéristique λ2 − 4λ + 5 = 0 a discriminant ∆ = −4 < 0 et deux racines complexes conjuguées de partie
réelle σ = 2 et partie imaginaire 1 et −1. La solution générale de l’EDO sur R est donc
Déterminons la valeur des constantes A et B correspondant à l’unique solution vérifiant y(π/2) = 1 et y 0 (π/2) = 1. On a
90 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
donc
B e π = 1,
( ( (
y(π/2) = 1, A = e −π ,
⇐⇒ ⇐⇒
y 0 (π/2) = 1, (2B − A)e π = 1, B = e −π .
y: R→R
t 7→ (cos(t ) + sin(t ))e 2t −π
(3) L’équation caractéristique λ2 + λ − 12 = 0 a discriminant ∆ = 49 > 0 et admet deux racines réelles distinctes, λ = −4 et
λ = 3. La solution générale de l’EDO sur R est donc
y(t ) = Ae −4t + B e 3t , A, B ∈ R.
Déterminons la valeur des constantes A et B correspondant à l’unique solution vérifiant y(0) = 1 et y 0 (0) = 4. On a
y 0 (t ) = −4Ae −4t + 3B e 3t ,
donc
( ( (
y(0) = 1, A + B = 1, A = −1/7,
⇐⇒ ⇐⇒
y 0 (0) = 4, −4A + 3B = 4, B = 8/7.
y: R→R
1
t 7→ (−e −4t + 8e 3t )
7
Exercice 4.40
Calculer la solution générale des EDO d’ordre 2 linéaires à coefficients constants suivantes :
Correction
y H (t ) = Ae −2t + B e −t , A, B ∈ R.
y P (t ) = (at + bt 2 + c t 3 )e −t ,
y P0 (t ) = (a + (2b − a)t + (3c − b)t 2 − c t 3 )e −t ,
y P00 (t ) = ((2b − 2a) + (6c − 4b + a)t + (−6c + b)t 2 + c t 3 )e −t ,
d’où
³¡ ´
y P00 (t ) + 3y P0 (t ) + 2y P (t ) = a + 2b + 2b + 6c t + 3c t 2 + e −t = (t 2 + 1)e −t
¢ ¡ ¢ ¡ ¢
© 2019-2020 G. Faccanoni 91
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
Par identification (deux polynômes sont égaux s’ils ont mêmes coefficients) il vient
a + 2b = 1
a = 3
2b + 6c = 0 ⇐⇒ b = −1
3c = 1 c = 1/3
yP : R → R
t 7→ (3t − t 2 + t 3 /3)e −t .
y: R→R
t 7→ Ae −2t + B e −t + (3t − t 2 + t 3 /3)e −t , A, B ∈ R.
y H (t ) = Ae −t + B e 3t , A, B ∈ R.
y P (t ) = At e 3t ,
y P0 (t ) = A(1 + 3t )e 3t ,
y P00 (t ) = A(6 + 9t )e 3t ,
d’où
y P00 (t ) − 2y P0 (t ) − 3y P (t ) = 4Ae 3t = 4e 3t
Par identification (deux polynômes sont égaux s’ils ont mêmes coefficients) il vient A = 1. On en déduit qu’une
solution particulière est donc
yP : R → R
t 7→ t e 3t .
y: R→R
t 7→ Ae −t + (B + t )e 3t , A, B ∈ R.
92 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
d’où
Par identification (deux polynômes sont égaux s’ils ont mêmes coefficients) il vient C = −1/2 et D = 0. On en
déduit qu’une solution particulière est donc
yP : R → R
t 7→ −t cos(t )/2.
y: R→R
t 7→ (A cos(t ) + B sin(t )) − t cos(t )/2, A, B ∈ R.
y P (t ) = C cos(3t ) + D sin(3t ),
y P0 (t ) = −3C sin(3t ) + 3D cos(3t ),
y P00 (t ) = −9C cos(3t ) − 9D sin(3t ),
d’où
Par identification (deux polynômes sont égaux s’ils ont mêmes coefficients) il vient C = 0 et D = −1/8. On en
déduit qu’une solution particulière est donc
yP : R → R
t 7→ − sin(3t )/8.
y: R→R
t 7→ (A cos(3t ) + B sin(3t )) − sin(3t )/8, A, B ∈ R.
yP : R → R
t 7→ −(− sin(3t )/8)/4 + 3(−t cos(t )/2)/4 = sin(3t )/32 − t cos(t )/8.
© 2019-2020 G. Faccanoni 93
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
y: R→R
t 7→ (A cos(3t ) + B sin(3t )) + sin(3t )/32 − t cos(t )/8, A, B ∈ R.
Exercice 4.41
On se propose de déterminer l’évolution au cours du temps de la charge q dans un circuit RLC soumis à une tension
sinusoïdale. La charge est solution de l’équation différentielle
1
Lq 00 (t ) + R q 0 (t ) + q(t ) = ων sin(ωt ) (4.3)
C
où ν est une constante réelle, L et ω sont des constantes réelles positives et où R et C sont des constantes strictement
positives.
F: R→R
t 7→ e αt cos(βt )
Correction
1. Soit f et g deux fonctions dérivables. Alors
Z Z
0
f (x)g (x) dx = f (x)g (x) − f 0 (x)g (x) dx.
e αt β
Z Z
I (t ) = e αt cos(βt ) dt = cos(βt ) + e αt sin(βt ) dt .
α α
De la même manière, on pose J (t ) = e αt sin(βt ) dt . En intégrant par parties ( f (t ) = sin(βt ) et g 0 (x) = e αt ) on trouve
R
e αt β e αt β
Z Z
J (t ) = e αt sin(βt ) dt = sin(βt ) − e αt cos(βt ) dt = sin(βt ) − I (t ).
α α α α
On en déduit
e αt β e αt β e αt β2
I (t ) = cos(βt ) + J (t ) = cos(βt ) + sin(βt ) − 2 I (t )
α α α α α α
d’où
α cos(βt ) + β sin(βt ) αt
Z
I (t ) = e αt cos(βt ) dt = e
α2 + β2
et
e αt β α sin(βt ) − β cos(βt ) αt
Z
J (t ) = e αt sin(βt ) dt = sin(βt ) − I (t ) = e
α α α2 + β2
qH : R → R
t 7→ κe −t /(RC ) avec κ ∈ R.
94 © 2019-2020 G. Faccanoni
Mercredi 4 septembre 2019 4 Équations différentielles ordinaires (EDO)
α=1/(RC ) Ã cos(ωt )
+ ω sin(ωt )
!
0 ων ¡ t /(RC ) ¢ β=ω ων RC t /(RC )
K (t ) = e sin(ωt ) =⇒ K (t ) = 1
e
R R
(RC )2
+ ω2
et donc la fonction
qP : R → R
sin(ωt ) − ωRC cos(ωt ) t /(RC )
µ ¶
t 7→ (ωνC ) e
1 + (ωRC )2
q: R→R
ωνC
t 7→ q H (t ) + q P (t ) = κe −t /(RC ) + (sin(ωt ) − ωRC cos(ωt )) e t /(RC )
1 + (ωRC )2
q P (t ) = α cos(ωt ) + β sin(ωt ),
q P0 (t ) = ω − α sin(ωt ) + β cos(ωt ) ,
¡ ¢
d’où
α´ β
µ ¶
1 ³
2 2
Lq P00 (t ) + R q P0 (t ) + q P (t ) = −Lω α + Rωβ + cos(ωt ) + −Lω β − Rωα + sin(ωt ) = ων sin(ωt )
C C C
ω2 νR
−Lω2 α + Rωβ + Cα = 0 α = − ¡ 1 −Lω2 ¢2 +(Rω)2
( (¡
1 2
− Lω α = −Rωβ
¢
β ⇐⇒ ¡ C1 ⇐⇒ C ¡
ων C1 −Lω2
¢
−Lω2 β − Rωα + C = ων − Lω2 β = ων − Rωα
¢
C β =
¡1 ¢2
2 +(Rω)2
C −Lω
qP : R → R
ων C1 − Lω2
¡ ¢
ω2 νR
t 7→ − ¡ cos(ωt ) + ¡ sin(ωt ).
1 2 2 + (Rω)2 1 2 2 + (Rω)2
¢ ¢
C − Lω C − Lω
© 2019-2020 G. Faccanoni 95
4 Équations différentielles ordinaires (EDO) Mercredi 4 septembre 2019
q: R→R
p p
ων C1 − Lω2
¡ ¢
R2− R 2 −4L/C
t R2+ R 2 −4L/C
t ω2 νR
t 7→ Ae 2L +Be 2L −¡ cos(ωt ) + ¡ sin(ωt ), A, B ∈ R.
1 2 2 + (Rω)2 1 2 2 + (Rω)2
¢ ¢
C − Lω C − Lω
96 © 2019-2020 G. Faccanoni
Bibliographie
[1] Guy AULIAC, Jean AVIGNANT et Elie A ZOULAY : Aide-mémoire de Mathématiques. EdiScience, 2006.
[2] Anne-Emmanuelle B ADEL et François C LAUSSET : Physique tout-en-un - 1re année. Dunod, 2008.
[3] Vivina B ARUTELLO, Monica C ONTI, Davide L. F ERRARIO, Susanna T ERRACINI et Gianmaria V ERZINI : Analisi matematica
con elementi di geometria e calcolo vettoriale, volume 2. Apogeo, 2008.
[4] Vincent B LONDEL : Mathématiques - Analyse. Dunod, 2000.
[5] Xavier B UFF, Josselin G ARNIER, Emmanuel H ALBERTSTADT, Thomas L ACHAND -R OBERT, François M OULIN et Jacques
S AULOY : Mathématiques tout-en-un pour la licence niveau L1. Dunod, 2006.
[6] Xavier B UFF, Josselin G ARNIER, Emmanuel H ALBERTSTADT, François M OULIN, Monique R AMIS et Jacques S AULOY :
Mathématiques tout-en-un pour la licence niveau L2. Dunod, 2007.
[7] Claudio C ANUTO et Anita TABACCO : Analisi matematica II - Teoria ed esercizi con complementi in rete. Springer, 2008.
[8] Alexandre C ASAMAYOU -B OUCAU, Pascal C HAUVIN et Guillaume C ONNAN : Programmation en Python pour les mathé-
matiques. Dunod, 2012.
[9] Yadolah D ODGE : Mathématiques de base pour économistes. Springer, 2007.
[10] Daniel F REDON, Myriam M AUMY-B ERTRAND et Frédéric B ERTRAND : Mathématiques Analyse en 30 fiches. Dunod,
2009.
[11] François G UÉNARD et Patricia H UG : QCM de Mathématiques, volume 1. Dunod, 1993.
[12] Wieslawa J. K ACZOR et Maria T. N OWAK : PROBLÈMES D’ANALYSE I - Nombres réels, suites et séries. EDP Sciences, 2008.
[13] Wieslawa J. K ACZOR et Maria T. N OWAK : PROBLÈMES D’ANALYSE II - Continuité et dérivabilité. EDP Sciences, 2008.
[14] Jean-Pierre L ECOUTRE et Philippe P ILIBOSSIAN : TD Analyse. Dunod, 2008.
[15] François L IRET et Charlotte S CRIBOT : Mini manuel d’Analyse. Dunod, 2010.
[16] Jean-Marie M ONIER : Les méthodes et exercices de Mathématiques PCSI-PTSI. Dunod, 2008.
[17] François M OULIN, Jean François RUAUD, Anne M IQUEL et Jean-Claude S IFRE : Mathématiques tout-en-un - 1re année.
Dunod, 2003.
[18] Bernard M YERS et Dominique S OUDER : Logique et Mathématiques. Dunod, 2009.
[19] James S TEWART : Calculus concepts and contexts. Brooks/Cole, 2010.
[20] James S TEWART : Calculus : early transcendentals. Brooks/Cole Pub Co, 2010.
[21] James S TEWART, Lothar R EDLIN et Saleem WATSON : Precalculus Mathematics for Calculus. Brooks/Cole, 2009.
97