0% ont trouvé ce document utile (0 vote)
121 vues3 pages

Résolution de Systèmes Linéaires

Transféré par

mauricelouba9
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
121 vues3 pages

Résolution de Systèmes Linéaires

Transféré par

mauricelouba9
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Année universitaire: 2024 − 2025

École Nationale Supérieure des Mines De


Rabat Analyse Numérique: Tronc commun

TD1: Résolution des Systèmes Linéaires et Non-linéaires

Exercice 1 :
Soit le système linéaire Ax = b suivant :

    
1 2 3 4 x1 10
2 3 4 1 x  10
 2  
 = 
 
3 4 1 2 x3  10
 

4 1 2 3 x4 10

1. Résoudre le systèmes linéaire par la méthode du pivot de Gauss.

2. Déduire le déterminant de A.

3. Donner les matrices de passage Mk vérifiant A(k+1) = Mk A(k) où A(k) est la matrice à
l’étape k avec A(1) = A.

4. Déduire les matrices L et U avec L une matrice triangulaire inférieur à diagonale unitaire
et A = LU .

5. Utiliser maintenant l’algorithme de décomposition LU pour retrouver les matrices L et U .

Exercice 2 :
Considérons la matrice symétrique défnie positive A et le vecteur b suivants :
   
9 −3 3 3
A = −3 10 −4 b = −4
   
3 −4 18 18

1. Calculer la factorisation la plus adéquate de A. Justifier votre choix.

2. Utiliser cette factorisation pour calculer la solution du système AX = b.

3. Les méthodes de Gauss-Seidel et Jacobi appliquées au système Ax = b seront-t-elles con-


vergentes si le point de départ x0 est bien choisi ?

On donne X (0) = [1, 1, 0]T .

4. Effectuer deux itérations de la méthode de Jacobi pour estimer la solution X.

ENSMR— RABAT 1ere année Tronc Commun I. MEDARHRI & [Link] 1


5. Effectuer deux itérations de la méthode de Gauss-Siedel pour estimer la solution X.

6. A chaque itération et pour chaque méthode, calculer l’écart entre la solution exacte X et
son approximée. Commenter le résultat..

Exercice 3 :
On considère le système linéaire A⃗x = ⃗b suivant:
" #" # " #
0, 780 0, 563 x1 0, 217
=
0, 913 0, 659 x2 0, 254

dont la solution exacte est ⃗x = [1, −1]T . En résolvant le système à l’aide d’une méthode
numérique, on a trouvé une solution approximative ⃗x∗ dont le résidu est donné par
" #
⃗b − A⃗x∗ = 0, 000001
0, 0

1. Sachant que A−1 ∞ = 1, 693 × 106 , est-ce que la solution approximative ⃗x∗ est accept-
able? Justifier votre réponse.

Exercice 4 :
Considérons le calcul numérique de la racine positive de la fonction :

f (x) = x3 − 3x − 2

1. Donner la formule itérative de Newton.

2. Calculer les solutions approchées xk , k = 1; 2; 3; 4, pour x0 = 2, 1.

Exercice 5 :
Soit la fonction f définie par f (x) = x3 − x2 + 8x − 8. Pour calculer ses racines on utilise 4
méthodes de point fixe différentes décrites par les fonctions d’itération suivantes:

8 − x3 1 3 1 1 4 2x3 − x2 + 8
φ1 (x) = −x3 + x2 − 7x + 8 φ2 (x) = φ3 (x) = − x + x2 + x + φ4 (x) =
8−x 10 10 5 5 3x2 − 2x + 8
Le tableau ci-dessous présente les suites des itérées obtenues par ces quatre méthodes.
Méth Méthode B Méthode C Méthode D
x0 0.5000000000000000 0.500000000000000 0.5000000000000000 0.5000000000000000
x1 0.9125000000000001 1.032258064516129 4.625000000000000 1.050000000000000
x2 0.9897857421875000 1.000235245684712 −101.9160156250000 0.9845143884892086
x3 0.9989578145726552 1.000000012299503 1.069697123778202 × 106 1.004312677086027
x4 0.9998955643403695 1.000000000000000 −1.224001861234915 × 1018 0.9987590594698483
x5 0.9999895542527895 1.000000000000000 1.833775789385161 × 1054 1.000353832012369
x6 0.9999989554034564 1.000000000000000 −6.166499545700052 × 10162 0.9998988463640411

1. Montrer que ℓ = 1 est l’unique racine réelle de f .

ENSMR— RABAT 1ere année Tronc Commun I. MEDARHRI & [Link] 2


2. Associer chaque méthode à sa fonction d’itération (justifier chaque réponse). Commenter
le résultat.

Exercice 6 :
Pour chaque valeur du paramètre s, on a une fonction :

gs (x) = s3 + x(1 − s) + (x − s2 )3

1. Montrer que pour chaque valeur de s, x = s2 est un point fixe de gs .

2. La méthode des points fixes converge vers le point fixe x = s2 pour quelles valeurs de
s? si la convergence est linéaire, donner le taux de convergence. Sinon donner l’ordre de
convergence.

Exercice 7 :

1. Obtenir tous les points fixes de la fonction :

g (x) = λx(1 − x)

où λ est un paramètre (λ ̸= 0).

2. Déterminer pour chaque point fixe trouvé en 1) les valeurs de λ pour lesquelles ces points
fixes sont attractifs.

3. Déterminer pour chaque point fixe trouvé en 1) la valeur de λ pour laquelle la convergence
de la méthode des points fixes sera quadratique.

Exercice 8 :
On considère la matrice:  
0 3 1
A = 4 1 1
 
2 2 4

1. Montrer que la décomposition LU de la matrice obtenue en permutant les lignes 1 et 2 de


la matrice A s’écrit P A = LU , où P est une matrice de permutation. Donner P , L et U .

2. Résoudre Ax = b, où b = (1, 5, 6)T à l’aide de la décomposition P A = LU .

ENSMR— RABAT 1ere année Tronc Commun I. MEDARHRI & [Link] 3

Vous aimerez peut-être aussi