Analyse Numérique2
Série de TD 2: Dérivation et intégration numérique
ENPConstantine: 2023-2024
�Exercice 1
1. Calculer la valeur exacte de l’intégrale
∫ 1 √
I= x + 1dx
−1
2. Approcher I en utilisant successivement :
• La méthode des Trapèzes avec 6 sous intervalles.
• La méthode de Simpson 1/3 avec 3 intervalles.
• La méthode de Simpson 3/8 avec 2 intervalles.
• La méthode de Gauss à 3 points.
�Exercice 2
1. Calculer la valeur exacte la dérivée de la fonction
f (x) = x2 ex
au point x = 0.
2. Pour h = 0.1, approcher f ′ (0) en utilisant successivement :
f (x + h) − f (x)
• f ′ (x) ≈ .
h
f (x) − f (x − h)
• f ′ (x) ≈ .
h
− f (x + 2h) + 4f (x + h) − 3f (x)
• f ′ (x) ≈ .
2h
f (x + h) − f (x − h)
• f ′ (x) ≈ .
2h
3f (x) − 4f (x − h) + f (x − 2h)
• f ′ (x) ≈ .
2h
�Exercice 3
Soient t0 , t1 ∈ [−1, 1]∫tels que t0 < t1 et c0 , c0 ∈ R. Pour toute fonction f continue sur [−1, 1] et pour
1
approcher l’intégrale −1 f (x)dx, on définit la formule de quadrature J de la façon suivante:
J(f ) = c0 f (x0 ) + c1 f (x1 )
2t1 2t0
1. Montrer que J est exacte d’ordre au moins 1 sur [−1, 1] si et seulement si c0 = et c1 =
t1 − t0 t0 − t1
2. Pour quelles valeurs de c0 , c1 , t0 et t1 , J est-elle au moins exacte d’ordre 3? Quel est alors l’ordre
d’exactitude de la méthode?
3. Déduire des questions précédentes une méthode d’intégration d’ordre 3 sur un intervalle [a, b] quel-
conque.
Page 1 sur 4
�Exercice 4
Soit f une fonction de classe au moins C 3 ([a, b]). On se fixe un nombre h > 0 suffisement petit ainsi qu’un
réel x quelconque de l’intervalle ]a, b[:
1. Montrer que la formulke aux différences finies D , définie par :
3f (x) − 4f (x − h) + f (x − 2h)
D=
2h
approche numériquement une dérivée de f que l’on déterminera
2. Donner l’ordre de précision de cette approximation.
3. Application:
Dans une réaction chimique, on a déterminé l’avancement x de la réaction en fonction du temps t:
t(s) 0 0.1 0.2 0.3 0.4
x(mole) 0 0.5 0.8 0.9 0.91
Déterminer la vitesse d’avancement de la réaction à l’instant t = 0.2(s)
�Exercice 5
Soient les points x0 = 0, x1 = 13 , x2 = 32 , x3 = 1
1. Déterminer P , le polynôme d’interpolation de Lagrange d’une fonction f aux points d’abscisses
(xi )i=0,1,2,3 .
∫
2. J est une formule de quadrature qui permet d’approximer l’intégrale 01 f (x)dx. Utilisez le polynôme
P pour obtenir une expression pour la quadrature J de la forme suivante :
J(f ) = c0 f (x0 ) + c1 f (x1 ) + c2 f (x2 ) + c3 f (x3 ) (1)
c0 , c1 , c2 et c3 sont des réels dont les valeurs doivent être précisées.
3. Cette formule est l’une des formes quadratiques que nous avons abordées dans le cours, veuillez la
spécifier précisément.
4. Application:
Pour
√ un dipôle alimenté à une tension u(t) de 5V et de période T = 0.1s, plus précisément: u(t) =
5 2cos(20πt), et parccouru par un courant i(t) mesuré 4 fois par période T :
t 0.0000 0.0334 0.0667 0.1000
i(t) 1.2000 1.1992 1.1968 1.1928
On définit la puissance active P par l’intégrale:
∫ T
1
P = u(t)i(t)dt
T 0
• Calculer P par la formule (1) .
�Exercice 6
A l’aide d’une certaine méthode d’intégration numérique, on a évalué l’intégrale
∫ π
2
I= sin(x)dx,
0
Page 2 sur 4
en utilisant trois valeurs différente de h. On a obtenu les résultats suivantes:
h I
0.1 1.001325
0.2 1.009872
0.4 1.078979
Compte tenu de la valeur exacte de I déduire l’ordre de convergence de la méthode de quadrature utilisée.
�Exercice 7
Soit f une fonction de classe au moins C 6 ([a, b]). On se fixe un nombre h > 0 suffisement petit ainsi qu’un
réel x quelconque de l’intervalle ]a, b[:
1. Montrer que le rapport
f (x + 2h) − 4f (x + h) + 6f (x) − 4f (x − h) + f (x + 2h)
D=
h4
approche numériquement une dérivée de f que l’on déterminera.
2. Donner l’ordre de précision de cette approximation.
3. En considérant la fonction f (x) = ex au point x0 = 0 et en utilisant successivement h = 0.1 et h = 0.3,
déterminer numériquement l’ordre de cette formule.
�Exercice 8
La fonction de Fresnel est une fonction spéciale définie par l’intégrale:
∫ x
F (x) = sin(t2 )dt
0
Admettant que : F (1) = 0.3102683017233811.
1. Appliquer la quadrature de Gauss Legendre à 2 points sur 4 sous intervalles de [0, 1] pour obtenir une
approximation de F (1).
2. Combien de points devriez-vous utiliser avec la méthode des Trapèzes composite pour obtenir une
précision comparable à celle obtenue en 1). Comparer les deux méthodes.
�Exercice 9
On approche la dérivée première de la fonction f (x) = xex en utilisant la formule de différences finies
suivante:
1
f ′ (x) ≈ Qh = (f (x + h) + f (x) − 2f (x − h))
3h
1. En utilisant ce schéma, approximer f ′ (0) pour h = 0.1 et h = 0.05. Comparer avec la valeur exacte.
2. En vous servant de la solution exacte et des résultats obtenus en 1), montrer que ce schéma est d’ordre
1.
3. A l’aide de la formule, dite de Richardson, donnée par la relation suivante:
( )
h
2n Q − Q(h)
2
f ′ (x) ≈ , où n est l’ordre de l’approximan Qh
2n − 1
calculer une nouvelle approximation de f ′ (0) pour h = 0.1.
4. A l’aide du schéma de différence centrée d’ordre 2, calculer une nouvelle approximation de f ′ (0)
avec h = 0.1.
Page 3 sur 4
�Exercice 10
Soient t0 < t1 < · · · < tm , m + 1 points distincts de l’intervalle [−1, 1] et soit {L0 , L1 , L2 , · · · , Lm } la base
de Lagrange de Pm associée à ces m + 1 points.
On définit la formule de quadrature J par :
∑
m
Jm (g) = ωj g(tj )
j=0
1. Montrer que la formule de quadrature J est exacte de dgré m si et seulement si
∫ 1
ωj = Lj (t)dt
−1
2. Application:
Soit α un nombre réel donné tel que 0 < α < 1, soient t0 = −α, t2 = 0, t3 = α:
• Trouver ω0 , ω1 , ω2 en fonction de α tel que
∫ 1
J2 (P ) = P (t)dx, ∀P ∈ P2
−1
1 ∫
• Existe-t-il α tel que J2 (P ) = −1 P (t)dt, pour tout P ∈ Pn avec n > 3? Si oui, quel est la valeur
maximale de n et que valent alors α et ω0 , ω1 , ω2 .
�Exercice 11
Soit f une fonction suffisement régulière, on définie la formule de différences finies suivante:
− f (x + 2h) + 16f (x + h) − 30f (x) + 16f (x − h) − f (x − 2h)
f ′′ (x) ≈
12h2
1. Montrer que l’ordre de l’approximation proposée est égal à 4.
2. Application:
L’ailette est un dispositif thermique visant à augmenter le transfert de chaleur. Dans un cadre sim-
plifié, l’équation qui régit la température de l’ailette en fonction de la position x est donnée par:
d2 T
− λ2 [T (x) − Tf ] = 0
dx2
• Tf : La température du fluide environnant.
• λ: Paramètre strictement positif caractérisant l’ailette.
Pour une ailette de longueur L = 0.22m, avec Tf = 30(°C)
x(m) 0.02 0.07 0.12 0.17 0.22
T (°C) 75 55 40 35 30
Donner une approximation pour le paramètre λ.
Page 4 sur 4
