Analyse Numérique et Calcul Scientifique

Corrigé de l’examen final

12/01/2016

Exercice 1. [8 pts] On considère sur l’intervalle [0, ⇡/4] la fonction f (x) = tan(x) + x 1. Soit
↵ > 0. On considère la méthode du point fixe : x0 donné et xn+1 = g↵ (xn ), n = 0, 1, . . ., avec
g↵ (x) = x ↵f (x).
1) [2 pts] Montrer que f admet un unique zéro x⇤ dans l’intervalle [0, ⇡/4].

Corrigé : La fonction f est continue et strictement croissante sur [0, ⇡/4] car f 0 (x) = 2 + tan2 (x)
est strictement positive sur [0, ⇡/4]. De plus, elle change de signe sur [0, ⇡/4] car f (0) = 1 < 0
et f (⇡/4) = ⇡/4 > 0. D’après le théorème des valeurs intermédiaires, f admet un unique zéro x⇤
dans l’intervalle ]0, ⇡/4[.
2) [2 pts] Montrer que si ↵ 2]0, 2/3[, alors la méthode converge vers x⇤ si x0 est suffisamment
proche de x⇤ .

Corrigé : On applique le théorème d’Ostrowski. La fonction g↵ est C 1 (et même C 1 ) au voisinage

de x⇤ et g↵0 (x) = 1 ↵(2 + tan2 (x)). Sur [0, ⇡/4], on a 2 6 2 + tan2 (x) 6 3 et comme ↵ > 0, alors
1 3↵ 6 1 ↵(2 + tan2 (x)) 6 1 2↵. Or ↵ 2]0, 2/3[, donc 1 < 1 3↵ et 1 2↵ < 1, donc
8x 2 [0, ⇡/4], |g↵0 (x)| < 1.
En particulier, on a |g↵0 (x⇤ )| < 1 car x⇤ 2 [0, ⇡/4].
3) [1 pt] Déterminer suivant les valeurs de ↵ l’ordre de convergence de la méthode.

Corrigé : On a g↵0 (x⇤ ) = 1 ↵(2 + tan2 (x⇤ )). Si

1
↵ 6= ,
2 + tan2 (x⇤ )
alors g↵0 (x⇤ ) 6= 0, donc la convergence est linéaire. Dans le cas contraire, la convergence sera au
moins quadratique. On aura alors g↵0 (x⇤ ) = 0 et
2 tan(x⇤ )(1 + tan2 (x⇤ ))
g↵00 (x⇤ ) = .
2 + tan2 (x⇤ )
Comme x⇤ 2]0, ⇡/4[, on a alors g↵00 (x⇤ ) 6= 0, donc si ↵ = 1/(2 + tan2 (x⇤ )), alors la convergence est
quadratique.
4) [1.5 pt] Ecrire l’algorithme de Newton associé à ce problème.

Corrigé : Algorithme de Newton :

(
x0 proche de x⇤ ,
xn+1 = xn ff0(x n)
(xn )
, n = 0, 1, 2, . . .

où f (xn ) = tan(xn ) + xn 1 et f 0 (xn ) = 2 + tan2 (xn ).

1
 5) [1.5 pt] Justifier que la méthode de Newton converge et déterminer son ordre de conver-
gence.

Corrigé : La fonction f est C 2 au voisinage de x⇤ , f (x⇤ ) = 0, f 0 (x⇤ ) = 2 + tan2 (x⇤ ) 6= 0 et

f 00 (x⇤ ) = 2 tan(x⇤ )(1 + tan2 (x⇤ )) 6= 0 (car x⇤ 2]0, ⇡/4[) donc d’après le théorème vu en cours, la
méthode de Newton converge de manière quadratique.

Exercice 2. [3 pts] On souhaite trouver (a, b) 2 IR2 tel que la fonction f :] ⇡/2, ⇡/2[ définie
par f (x) = a cos(x) + b tan(x) passe au plus près des points (0, 1), (⇡/6, 0) et (⇡/4, 2).
1) [2 pts] Montrer que ce problème revient à minimiser la fonction g : IR2 ! IR définie par
✓ ◆ 2
a
g(a, b) = A y
b

où la matrice A et le vecteur y sont à expliciter.

Corrigé : On pose (x1 , y1 ) = (0, 1), (x2 , y2 ) = (⇡/6, 0) et (x3 , y3 ) = (⇡/4, 2). On souhaite
minimiser sur IR2 la fonction g : IR2 ! IR définie par
3
X ✓ ◆ 2
2 a
g(a, b) = (a cos(xi ) + b tan(xi ) yi ) = A y
b
i=1

où on a posé :
0 1 0 1 0 1
cos(0) tan(0) p1 0p 1
A = @cos(⇡/6) tan(⇡/6)A = @ 3/2
p 1/ 3
A et y = 0A .
@
cos(⇡/4) tan(⇡/4) 1/ 2 1 2

2) [1 pt] Déterminer le système linéaire associé à ce problème sans le résoudre .

Corrigé : La matrice A est de rang 2, donc le problème de minimisation de la question 1) admet

une solution unique (a, b) 2 IR2 solution de l’équation normale :
✓ ◆
T a
A A = AT y.
b

On a 0 1
✓ p p ◆ 1 0p ✓ p ◆
1 3/2 1/ 2 p 9/4 (1 + 2)/2
AT A = p @ 3/2 1/ 3A = p
0 1/ 3 1 p (1 + 2)/2 4/3
1/ 2 1
et 0 1
✓ p p ◆ 1 ✓ p ◆
1 3/2
p 1/ 2 @0A = 1 + 2
AT y = .
0 1/ 3 1 2
2
Le système linéaire à résoudre est donc :
( p p
9 1+ 2
a
4 p
+ 2
b=1+ 2
1+ 2
.
2
a + 43 b = 2

2
Supplément : En utilisant les formules de Cramer, on a :
p p
1 + 2 (1 + 2)/2 p
2 4/3 4(1 + 2)
a= p = p
9/4
p (1 + 2)/2 3(9 2 2)
(1 + 2)/2 4/3
et p
9/4
p 1 + 2
p
(1 + 2)/2 2 4(3 2)
b= p = p .
9/4
p (1 + 2)/2 9 2 2
(1 + 2)/2 4/3

Exercice 3. [9 pts]
Soit f 2 C 4 ([ 1, 1]). On considère sur [ 1, 1] la formule de quadrature suivante :
Z 1
f (x)dx ⇡ af (b) + af ( b). (1)
1

où a et b sont des réels strictement positifs.

1) [2.5 pts] Trouver a et b pour que cette formule de quadrature soit d’ordre le plus élevé
possible. Quel est-il ?

Corrigé : On souhaite que la formule de quadrature soit exacte lorsque f (x) = 1, x, x2 , x3 , . . . , xn

avec n le plus grand possible. Si f (x) = 1 alors
Z 1 Z 1
f (x)dx = dx = 2 et af (b) + af ( b) = a + a = 2a.
1 1

D’où 2a = 2 donc a = 1 . Si f (x) = x alors f est impaire donc

Z 1 Z 1
f (x)dx = xdx = 0 et af (b) + af ( b) = ab ab = 0.
1 1

La formule de quadrature est donc vérifiée pour f (x) = x. Si f (x) = x2 , alors

Z 1 Z 1
2
f (x)dx = x2 dx = et af (b) + af ( b) = 2ab2 = 2b2 .
1 1 3
p
D’où b2 = 1/3. Comme b > 0, il en résulte que b = 1/ 3 . Pour f (x) = x3 , on a :
Z 1
f (x)dx = 0 et af (b) + af ( b) = ab3 ab3 = 0.
1

La formule de quadrature est donc vérifiée pour f (x) = x3 . Pour f (x) = x4 , on a :

Z 1 Z 1
2 2
f (x)dx = x4 dx = et af (b) + af ( b) = 2ab4 = .
1 1 5 9

Or 2/5 6= 2/9. Donc la formule n’est pas exacte pour f (x) = x4 .

3
 p
Avec a = 1 et b = 1/ 3, la formule de quadrature
Z 1 ✓ ◆ ✓ ◆
1 1
f (x)dx ⇡ f p +f p
1 3 3
est exacte pour tout polynôme de degré 6 3 : elle est donc d’ordre 3.
2) [1 pt] Quelle formule de quadrature retrouve-t-on ?

Corrigé : Il s’agit de la formule de quadrature de Gauss-Legendre à 2 points.

3) [1 pt] A l’aide d’un changement de variable, montrer que sur l’intervalle [0, h] (h 2]0, 1[),
la formule de quadrature (1) s’écrit :
Z h ⇢ ✓ ◆ ✓ ◆
h h h h h
f (x)dx ⇡ f p +f + p . (2)
0 2 2 2 3 2 2 3

Corrigé : A l’aide du changement de variable

2x
u= 1
h
on a x = h2 (u + 1) donc
Z h Z ✓ ◆
h 1 h
f (x)dx = f (u + 1) du
0 2 1 2
⇢ ✓ ✓ ◆◆ ✓ ✓ ◆◆
h h 1 h 1
⇡ f p +1 +f p +1
2 2 3 2 3
⇢ ✓ ◆ ✓ ◆
h h h h h
⇡ f + p +f p .
2 2 2 3 2 2 3

4) [2.5 pts] Montrer qu’il existe une constante C > 0 indépendante de h telle que :
Z h ⇢ ✓ ◆ ✓ ◆
h h h h h
f (x)dx f p +f + p 6 C max |f (4) (x)|h5 .
0 2 2 2 3 2 2 3 x2[0,h]

Corrigé : On pose
h h h h
+ p
x+ = et x = p .
2 2 3 2 2 3
Soit L le polynôme de Lagrange de degré 6 3 qui interpole f en x+ , x , 0 et h. On a en
particulier L(x ) = f (x ) et L(x+ ) = f (x+ ). La formule de quadrature (2) étant d’ordre 3 (car
le changement de variable est affine), on a
Z h
h h
L(x)dx = {L(x ) + L(x+ )} = {f (x ) + f (x+ )} .
0 2 2
D’après le cours, on sait que pour tout x 2 [0, h], il existe ⇠ 2]0, h[ tel que

f (4) (⇠)
f (x) L(x) = (x 0)(x h)(x x )(x x+ ).
4!

4
Par conséquent, pour tout x 2 [0, h],
1
|f (x) L(x)| 6 max |f (4) (x)|h4 .
4! x2[0,h]

Donc
Z h ⇢ ✓ ◆ ✓ ◆ Z h
h h h h h h
f (x)dx f p +f + p = f (x)dx {f (x ) + f (x+ )}
0 2 2 2 3 2 2 3 0 2
Z h Z h
= f (x)dx L(x)dx
0 0
Z h
= (f (x) L(x))dx
0
Z h
6 |f (x)
L(x)| dx
0
Z h
1
6 (4)
max |f (x)|h 4
dx
4! x2[0,h] 0
1
6 max |f (4) (x)|h5 .
4! x2[0,h]

5) Soit g 2 C 4 ([0, 1]). On considère une subdivision de [0, 1] en n intervalles de longueur h.

On pose ti = ih pour i = 0, 1, . . . n.

a) (1 pt) Montrer que la formule de quadrature (1) répétée sur [0, 1] est donnée par
Z n 1⇢ ✓ ◆ ✓ ◆
1
hX h h h h
g(x)dx ⇡ g p + ih + g + p + ih .
0 2 i=0 2 2 3 2 2 3

Corrigé : On a
Z 1 n 1Z
X ti+1
g(x)dx = g(x)dx.
0 i=0 ti

Avec le changement de variable t = x ti , obtient

Z ti+1 Z h Z h
g(x)dx = g(t + ti )dt = g(t + ih)dt.
ti 0 0

Ainsi
Z 1 n 1Z
X h
g(x)dx = g(t + ih)dt.
0 i=0 0

En appliquant la formule de quadrature (2), on obtient

Z h ⇢ ✓ ◆ ✓ ◆
h h h h h
g(t + ih)dt ⇡ g p + ih + g + p + ih .
0 2 2 2 3 2 2 3

5
 b) (1 pt) Montrer qu’il existe une constante C > 0 indépendante de h telle que
Z n 1⇢ ✓ ◆ ✓ ◆
1
hX h h h h
g(x)dx g p + ih + g + p + ih 6 CM4 h4
0 2 i=0 2 2 3 2 2 3

où M4 = maxx2[0,1] |g (4) (x)|.

Corrigé : On a
Z n 1⇢ ✓ ◆ ✓ ◆
1
hX h h h h
g(x)dx g p + ih + g + p + ih
0 2 i=0 2 2 3 2 2 3
n 1Z
X n 1⇢ ✓ ◆ ✓ ◆
h
hX h h h h
= g(t + ih)dt g p + ih + g + p + ih
i=0 0 2 i=0 2 2 3 2 2 3
X1
n Z h ⇢ ✓ ◆ ✓ ◆
h h h h h
6 g(t + ih)dt g p + ih + g + p + ih .
i=0 0 2 2 2 3 2 2 3

En appliquant le résultat de 4) à la fonction t 7! g(t + ih), on obtient

Z h ⇢ ✓ ◆ ✓ ◆
h h h h h
g(t + ih)dt g p + ih + g + p + ih
0 2 2 2 3 2 2 3
1 1 1
6 max |g (4) (t + ih)|h5 6 max |g (4) (t)|h5 6 max |g (4) (t)|h5 .
4! t2[0,h] 4! t2[ti ,ti+1 ] 4! t2[0,1]

Il en résulte que
Z n 1⇢ ✓ ◆ ✓ ◆
1
hX h h h h
g(x)dx g p + ih + g + p + ih
0 2 i=0 2 2 3 2 2 3
n 1
X 1 1 1
6 max |g (4) (t)|h5 6 max |g (4) (t)|nh5 6 max |g (4) (t)|h4
i=0
4! t2[0,1] 4! t2[0,1] 4! t2[0,1]

car nh = 1.