Faculté de la Technologie Année Universitaire 2008/2009

Département 2ème Année ST 2ème Année Sciences Technique

Module : Math VI
Epreuve de Moyenne Durée
Exercice N°1 : (08 points)
On considère la fonction à valeurs réelles f (x ) = 4e x/4 − 6 dans l’intervalle ]0, 4[ .
1-Montrer qu’il existe un zéro α pour la fonction f dans l’intervalle ]0, 4[ et trouver α de
façon analytique. (01 pt)
2-Peut-on appliquer la méthode de dichotomie pour calculer α ? Justifier votre réponse. (01
pt)
3-Pour approcher le zéro α on considère les méthodes de point fixe x (k +1) = φ i x (k ) , i = 1,2,3 , ( )
avec :
e x/4 3
φ1 (x ) = x + 4e x/4 − 6 ; φ 2 (x ) = x − 4 + 6e -x/4 ; φ 3 (x ) = x - +
16 32
Etablir si les trois méthodes sont convergentes (03 pts)
e x/4 3
4-Pour la méthode de fonction φ 3 (x ) = x - + , déterminer le nombre minimal d’itérations
16 32
nécessaires pour avoir une erreur inférieure à 10 −6 , lorsqu’on a choisi x (0 ) tel que x (0 ) − α < 2 . (03
pts)

Exercice N°2 : (06 points)

On considère le système AX = b
⎡ 3 1⎤ ⎧x ⎫ ⎧5⎫
Où A=⎢ ⎥ ; X = ⎨ 1⎬ ; b = ⎨ ⎬
⎣ − 1 4⎦ ⎩x 2 ⎭ ⎩- 6⎭
1-Calculer la matrice G de Gauss-Seidel associée au système. (01 pt)
2- En partant de la relation AX = b , montrer que (D − E )−1 b = X − GX . (01 pt)
et en déduire que : (
X (k +1) − X (k ) = −(I − G ) X (k ) − X ) (01 pt)
3- Déduire alors de ce qui précède que : X (k ) − X = (I − G )−1 X (k +1) − X (k ) (01 pt)
2 2 2

4- Calculer (I − G )−1 et en déduire que X (k ) − X ≤ 2 X (k +1) − X (k ) (01 pt)

2 2 2

5-En déduire une condition suffisante pour avoir X (k ) − X ≤ 10 −3 . (01 pt)

2

Exercice N°2 : (06 points)

Considérons le système linéaire AX = b dont la matrice est tri-diagonale et inversible :
⎡ a1 c1 0 . 0 ⎤
⎢b a2 c2 . 0 ⎥⎥
⎢ 2
A=⎢0 . . . . ⎥
⎢. . . . . ⎥
⎢ ⎥
⎢. . bn −1 a n −1 c n −1 ⎥
⎣⎢ 0 . 0 bn a n ⎦⎥
Les matrices de la factorisation LU de la matrice A sont bi-diagonales de la forme :
⎡α 1 γ 1 0 . 0 ⎤
⎡1 0 0 . 0⎤
⎢β ⎢0 α γ2 0 ⎥⎥
⎢ 2 1 0 . 0⎥⎥ ⎢ 2 .
L=⎢ 0 . . . .⎥ , U =⎢0 . . . . ⎥
⎢ . . . . .⎥ ⎢. . . . . ⎥
⎢ ⎥ ⎢ ⎥
⎢ . . β n −1 1 0⎥ α n −1 γ n −1 ⎥
⎣⎢ 0 . 0 βn 1⎦⎥ ⎢. . 0
⎢⎣ 0 . 0 0 α n ⎥⎦
Donner les expressions des coefficients α i , β i et γ i en fonction des coefficients de A , en utilisant
A = LU .

Page 1/5
Faculté de la Technologie Année Universitaire 2008/2009
Département 2ème Année ST 2ème Année Sciences Technique
Module : Math VI
Corrigé Epreuve de Moyenne Durée
EXERCICE 1- : (08 Pts)
1-on a f (0 ) = −2 < 0 , f (4 ) = 4e − 6 > 0 ; f (x ) étant continue, ∃α ∈ [0,4] tel que f (α ) = 0
f ( x ) est croissante sur l’intervalle [0,4] donc α est unique. (0.5 pt)
3
Son expression analytique est donnée par : f (α ) = 4e α / 4 − 6 = 0 ==> α = 4 ln
2
(0.5 pt)
2- Vu que f (0 ) f (4 ) < 0 et f ( x ) est croissante sur l’intervalle [0,4] ; on peut bien appliquer
la méthode de bissection pour calculer le zéro α . (1 pt)

3. Pour prouver que les méthodes sont convergentes, il faut vérifier que :
1- ∀x ∈ [0,4] , Φ i ( x ) ∈ [0,4]
2- Φ i' ( x) < 1 pour x ∈ ]0,4[
Φ 1 ( x ) = x + 4e x / 4 − 6 , Φ 1' ( x ) = 1 + e x / 4 est toujours positif dans l’intervalle [0,4] , donc
Φ 1 ( x ) = x + 4e x / 4 − 6 est strictement croissante.
5
On a Φ 1' (α ) = > 1 donc la méthode du point fixe x = x + 4e x / 4 − 6 diverge. (1 pt)
2
6
Φ 2 ( x ) = x − 4 + 6e − x / 4 , Φ '2 ( x ) = 1 − e − x / 4 change de signe dans l’intervalle [0,4] . La
4
condition 2 n’est pas vérifiée donc la méthode x = x − 4 + 6e − x / 4 diverge (1 pt)

ex/4 3 ex/4
Φ 3 (x ) = x − + , Φ 3' ( x ) = 1 − est toujours positif dans l’intervalle [0,4] , donc
16 32 64
ex/4 3
Φ 3 (x ) = x − + est strictement croissante.
16 32
125 ex/4 3
On a Φ 3' (α ) = < 1 , la méthode du point fixe x = x − + converge (1 pt)
128 16 32

4- On a :
( ) ( )
x (k +1) − α = Φ 3 x ( k ) − Φ 3 (α ) = Φ 3' (ξ ) x (k ) − α ≤ C x (k ) − α (1 pt)
1
Où ξ ∈ (0,4 ) et C = max ξ ∈[0, 4 ] Φ 3' (ξ ) = 1 − (1 pt)
64
Comme
x (k ) − α ≤ C n x (0 ) − α = 2C n , on veut trouver le plus petit des n tels que 2C n < 10 −6 . Donc
ln 2 + 6 ln 10 ln 2 + 6 ln 10
n> = ==> n ≥ 921 (1 pt)
− ln C ln (64 / 63)

Page 2/5
Faculté de la Technologie Année Universitaire 2008/2009
Département 2ème Année ST 2ème Année Sciences Technique
Module : Math VI
EXERCICE 1 : (06 Pts)
On considère le système : AX = b
⎡ 3 1⎤ ⎧x ⎫ ⎧5⎫
A=⎢ ⎥ ; X = ⎨ 1⎬ ; b = ⎨ ⎬
⎣ − 1 4⎦ ⎩ x2 ⎭ ⎩− 6 ⎭
1-Calcul de la matrice de Gauss-Seidel G associée au système :
On pose A = D − E − F
⎡3 0 ⎤ ⎡0 0 ⎤ ⎡0 − 1⎤
D=⎢ ⎥ , E=⎢ ⎥ et F = ⎢ ⎥
⎣0 4 ⎦ ⎣1 0⎦ ⎣0 0 ⎦
On a G = (D − E ) F
−1

⎡1 ⎤
⎡ 3 0⎤ ⎢ 0⎥
, (D − E ) =⎢3
−1
D−E = ⎢ ⎥
⎣ − 1 4⎦ 1 1⎥
⎢ ⎥
⎣12 4⎦
⎡1 ⎤ ⎡ 1⎤
⎢ 0 ⎥ ⎡0 − 1⎤ ⎢ 0 − ⎥
G = (D − E ) F = ⎢ 3 3
−1
, G=⎢
1 1 ⎥ ⎢⎣0 0 ⎥⎦ 1⎥
⎢ ⎥ ⎢0 − ⎥
⎣12 4 ⎦ ⎣ 12 ⎦
2-On pose A = D − E − F
⎡3 0 ⎤ ⎡0 0 ⎤ ⎡0 − 1⎤
Où D = ⎢ ⎥ , E=⎢ ⎥ ; A=⎢ ⎥ (1 pt)
⎣0 4 ⎦ ⎣1 0⎦ ⎣0 0 ⎦

2-Le système AX = b s’écrit alors (D − E − F )X = b

D’où (D − E )X = FX + b => X = (D − E ) FX + (D − E ) b = GX + (D − E ) b
−1 −1 −1

=> (D − E )−1 b = X − GX (1 pt)

On écrit le processus itératif :

X ( K +1) − GX ( K ) = (D − E ) b => X ( K +1) − X ( K ) − GX ( K ) = (D − E ) b − X ( K )
−1 −1
(a)
X − GX = (D − E ) b => X − X − GX = (D − E ) b − X
(K ) (K )
−1 −1
(b)
‘a)=(b) => X ( K +1) − X ( K ) = X ( K ) − X − GX + GX ( K )
X ( K +1) − X ( K ) = −( I − G )(X ( K ) − X ) (1 pt)

(
3- On a X ( K +1) − X ( K ) = −( I − G ) X ( K ) − X )
=> ( ) (
X ( K ) − X = −( I − G ) −1 X ( K +1) − X ( K ) )
=> X (K ) − X
2
= ( I − G ) −1 (X ( K +1)
− X (K )
) 2
= ( I − G ) −1
2
X ( K +1) − X ( K )
2
(1 pt)

Page 3/5
Faculté de la Technologie Année Universitaire 2008/2009
Département 2ème Année ST 2ème Année Sciences Technique
Module : Math VI
4- Calcule de ( I − G ) −1
2

⎡ 1⎤ ⎡ 1⎤ ⎡ 4⎤
0 − 1 ⎢ 1 −
⎢ 3 ⎥, I −G = ⎢ 3⎥ 13 ⎥
⎢ 13 ⎥ (I − G ) = ⎢
−1
On a G = ⎢
1⎥ 12 ⎥
⎢0 − ⎥ ⎢0 ⎥ ⎢0 ⎥
⎣ 12 ⎦ ⎣ 12 ⎦ ⎣ 13 ⎦

1
⎡n ⎛ n ⎞⎤ 2
On a la définition A 2 = ⎢∑ ⎜⎜ ∑ aij
2
⎟⎥
⎟
⎣⎢ i =1 ⎝ j =1 ⎠⎦⎥

1
⎡ 16 144 ⎤ 2 329
(I − G )−1 = ⎢1 + + ⎥ = = 1,9467 = 1,3952
2
⎣ 169 169 ⎦ 169

D’où X ( K ) − X = 1,3952 X ( K +1) − X ( K ) ≤ 2 X ( K +1) − X ( K ) (01 pt)

2 2 2

La condition suffisante pour avoir X (k ) − X ≤ 10 −3

2

On a X ( K ) − X ≤ 2 X ( K +1) − X ( K )
2 2

X (k ) − X ≤ 10 −3 ==> 1,3952 X ( K +1) − X ( K ) ≤ 10 −3 ==> X ( K +1) − X ( K ) ≤ 0.716.10 −3

2 2 2

(01 pt)

Page 4/5
Faculté de la Technologie Année Universitaire 2008/2009
Département 2ème Année ST 2ème Année Sciences Technique
Module : Math VI
EXERCICE 3 : (06 Pts)
On a
⎡ a1 c1 0 . . . 0 ⎤
⎢b a c . . . 0 ⎥⎥
⎢ 2 2 2

⎢. ⎥
⎢ ⎥
A=⎢ . ⎥
⎢. ⎥
⎢ ⎥
⎢0 . . . bn −1 a n −1 c n −1 ⎥
⎢0 . . . 0 c n ⎥⎦
⎣ bn
⎡1 0 0 . . . 0⎤ ⎡α 1 γ 1 0 . . . 0 ⎤
⎢β 1 0 . . . ⎥
0⎥ ⎢0 α γ2 . . . 0 ⎥⎥
⎢ 2 ⎢ 2

⎢ . ⎥ ⎢. ⎥
⎢ ⎥ ⎢ ⎥
L=⎢ . ⎥ U =⎢ . ⎥
⎢ . ⎥ ⎢. ⎥
⎢ ⎥ ⎢ ⎥
⎢ 0 . . . β n −1 1 0⎥ ⎢0 . . . 0 α n −1 γ n −1 ⎥
⎢0 . . . 0 βn ⎥
1⎦ ⎢ α n ⎥⎦
⎣ ⎣0 . . . 0 0
Le produit LU donne

⎡ α1 γ1 0 . . . 0 ⎤
⎢α β β 2γ 1 + α 2 γ 2 . . . 0 ⎥
⎢ 1 2 ⎥
⎢ . ⎥
⎢ ⎥
LU = ⎢ . ⎥
⎢ . ⎥
⎢ ⎥
⎢ 0 . . . α n − 2 β n −1 β n −1γ n − 2 + α n −1 γ n −1 ⎥
⎢ 0 . . . 0 α n −1 β n β n γ n −1 + α n ⎥⎦
⎣
⎡ a1 c1 0 . . . 0 ⎤
⎢b a c . . . 0 ⎥⎥
⎢ 2 2 2

⎢. ⎥
⎢ ⎥
=⎢. ⎥=A
⎢. ⎥
⎢ ⎥
⎢0 . . . bn −1 a n −1 c n −1 ⎥
⎢0 . . . 0 c n ⎥⎦
⎣ bn
D’où
b
β i = i pour 2 < i < n (02 pts)
ai −1
γ i = ci pour 1 < i < n − 1 (02 pts)
b
α 1 = a1 , α i = ai − β i γ i −1 = ai − i ci −1 pour 2 < i < n (02 pts)
ai −1

Page 5/5