1.3 Les méthodes directes
1.3.2 Méthode de Gauss méthode LU. Soit A ? Mn(IR) une matrice inversible
Résolution des syst`emes linéaires Méthode de Gauss
Aide : on cherchera d 'abord une relation de récurrence entre Nn et Nn?1. 3. Méthode de Gauss. Transformation de A en une matrice triangulaire supérieure.
Annexe 3 : Inversion de matrices par la méthode du pivot de Gauss
Dans le cas général on utilise la méthode du pivot de Gauss. Pour montrer qu'une matrice M est inversible : On applique les opérations élémentaires : •
Matrices: Gaussian & Gauss-Jordan Elimination
Matrices Handout- Gaussian and Gauss-Jordan. Updated: Fall 2019 Gaussian elimination is a method for solving systems of equations in matrix form.
Étape A : processus délimination de Gauss
Dans chaque cas on écrira les étapes de la méthode sous forme matricielle. 2. (algo) Soit M ? Mn(R) une matrice carrée inversible et soit b ? Rn un
Cours 1: Autour des systèmes linéaires Algorithme du pivot de
Methode plus "automatique" : le pivot de Gauss sur les sytémes linéaires Une méthode pour inverser une matrice : Pivot de Gauss. L'algorithme général.
Méthode de Gauss-Jordan Calcul de linverse dune matrice
Calcul de l'inverse d'une matrice. Méthodes numériques 2003/2004 - D.Pastre licence de mathématiques et licence MASS. 1. Méthode de Gauss-Jordan.
Méthode de Gauss-Jordan Calcul de linverse dune matrice
Calcul de l'inverse d'une matrice. Méthodes numériques 2003/2004 - D.Pastre licence de mathématiques et licence MASS. 1. Méthode de Gauss-Jordan.
Chapitre 3 Méthodes directes de résolution des syst`emes linéaires
Méthodes directes de résolution. 22. 3.2 Méthode d'élimination de Gauss. Le principe de la méthode consiste `a déterminer une matrice M inversible telle que
METHODES NUMERIQUES
5.3 Propriétés des matrices triangulaires unitaires . . . . . . . . . . . . . 28. 6 Factorisation LU. 31. 6.1 Formalisation de l'élimination de Gauss .
Exercice 5
1. Résoudre le système linéaireAx=bpar la méthode d"élimination de Gauss dans les trois cas suivants :
a- A=242 4 6
2 1 1 11 23 5 b=2 441 53
5 b- A=2
42 4 6
24 111 23 5 b=2 412
5 03 5 c- A=2
42 4 6
24 112 23 5 b=2 412
5 13 5 oub=2 42
2 43
5 Dans chaque cas, on écrira les étapes de la méthode sous forme matricielle.
2.(algo)SoitM2 Mn(R)une matrice carrée inversible et soitb2Rnun vecteur (b2 Mn;1(R)). Écrire l"algorithme
d"élimination de Gauss pour résoudre le système linéaireMx=b.Question 1a.Étape A : processus d"élimination de Gauss
On transforme le systèmeAx=ben un système triangulaire supérieurUx=~boùUest une matrice triangulaire supérieure
à diagonale unité en faisant des combinaisons linéaires des équations (ce qui équivaut à faire des combinaisons linéaires des
lignes deA).Étape A1 :
On éliminex1dans les équations2et3(en faisant une combinaison linéaire entre la première et la deuxième équations d"une
part, et entre la première et la troisième équations d"autre part) : L2 L2+L1etL3 L3+12
L1:(1)
On obtient le système
A1x=b1avecA1=2
42 4 6
0 5 70 1 53
5 etb1=2 443 33
5
En terme matriciel, ceci revient à multiplier les systèmeAx=bà gauche par par la matrice triangulaire inférieure
L 1=241 0 0
1 1 00:5 0 13
5L"étape A1 se résume donc à
Ax=b,A1x=b1avecA1=L1A;etb1=L1b:(2)
L1est une matrice triangulaire inférieure inversible à diagonale unité.
Étape A2 :
On éliminex2dans la troisième équation (pour cela, on fait une combinaison linéaire des équations 2 et 3) :
L3 L3 L2=5:(3)
On obtient le système
A2x=b2avecA2=2
41 2 3
0 1 750 0185
3 5 etb2=2 443 185
3 5 (4) 1
En terme matriciel, ceci revient à multiplier les systèmeA1x=b1à gauche par par la matrice triangulaire inférieure
L 2=241 0 0
0 1 001=5 13
5L"étape A2 se résume donc à
A1x=b1,A2x=b2avecA2=L2A1etb2=L2b1:(5)
L2est une matrice triangulaire inférieure inversible à diagonale unité.
Bilan étape A :
La matriceU=A2est une matrice triangulaire supérieure. Ainsi, le systeme(4)(qui peut être réécritUx=b2) est un
système triangulaire supérieur qui va être facile à résoudre à l"étape B. Matriciellement, en combinant(2)et(5)nous pouvons résumer l"étape A comme suit :La matrice
~Lest le produit de deux matrices triangulaires inférieures à diagonales unité. Donc la matrice~Lest une matrice
triangulaire inférieure inversible à diagonale unité. Son inverse, notéeL, est donc aussi triangulaire inférieure à diagonale unité.
On a donc montré queA=LUoùLest une matrice triangulaire inférieure à diagonale unité etUest une matrice triangulaire
supérieure. Autrement dit, la première étape de la méthode du pivot revient à faire de manière implicite la décompositionLU
deA(noter queLy=b, cf. exercice1, question 4).En pratique, lorsqu"on code l"algorithme d"élimination de Gauss, on ne code pas sa version matricielle mais on transforme
progressivement la matriceAet le second membreben faisant combinaisons lineaires des équations (des lignes) (cf.(1)-(3)).
Autrement dit, on n"a pas besoin de calculer pas explicitement la matriceL. Etape B : on résout le système triangulaire supérieurUx=y Ce système est résolu par substitution (cf. exercice 1, question 3). On obtient finalement x=2 412 13 5
Question 1b.
On procède comme dans la question1mais on va voir apparaitre une difficulté supplémentaire due à la présence d"unpivot
nul. On peut remarquer que det2 4 24= 0 si bien que le mineur principal deAd"ordre2est nul.
Étape A : processus d"élimination de Gauss
Étape A1
On éliminex1dans les équations2et3:
L2 L2+L1etL3 L3+12
L1:(7)
On obtient le système
A1x=b1avecA1=2
42 4 6
0 0 70 1 53
5 etb1=2 4127 63
5
En terme matriciel, ceci revient à multiplier les systèmeAx=bà gauche par par la matrice triangulaire inférieure
L 1=241 0 0
1 1 00:5 0 13
5 2L"étape A1 se résume donc à
Ax=b,A1x=b1avecA1=L1A;etb1=L1b:(8)
L1est une matrice triangulaire inférieure inversible à diagonale unité.
Étape A2 : permutation des lignes 2 et 3
Comme le pivotA1(2;2) = 0, on ne peut pas poursuivre directement le processus d"élimination. En particulier, on ne peut
pas éliminerx2dans l"équation3en faisant une combinaison entre la deuxième et la troisième ligne. Cependant, comme
A1(3;2)6= 0, nous allons échanger les lignes2et3.
L2 L3L3 L2:
On obtient alors
A2x=b2avecA2=2
42 4 6
0 1 50 0 73
5 etb2=2 4126 73
5
En terme matriciel, cela revient à multiplier le systèmeA1x=b1par la matrice de permutationP23donnée par
P 23=241 0 0
0 0 10 1 03
5Cette étape A2 se résume donc à
A1x=b1,A2x=b2avecA2=P23A1etb2=P23b1:(9)
Remarque.Comme la matriceA2est de dimension33, il n y a pas besoin de continuer le processus d"élimination (puisque
(A2)3;2= 0). Bien sûr, il faudrait la réaliser dans le cas d"une matrice plus grande. Cela reviendrait alors à multiplier le système
A2x=b2par une matrice triangulaire inférieure.
Bilan étape A
Matriciellement, en combinant(8)-(9), nous pouvons résumer l"étape A comme suit :On remarque queM=P2;3L1=~LP2;3, avec
L=241 0 0
0:5 1 0
1 0 13
5Lest une matrice triangulaire inférieure inversible à diagonale unité. Son inverse, notéeL, est donc aussi triangulaire inférieure
à diagonale unité.
On a donc montré queU=~LP2;3A, c"est à dire que queLU=P2;3A, oùLest une matrice triangulaire inférieure à
diagonale unité,Uest une matrice triangulaire supérieure etP2;3est une matrice de permutation. Autrement dit, la première
étape de la méthode du pivot revient à faire de manière implicite la décompositionPA=LU.
Etape B : on résout le système triangulaire supérieurUx=y Ce système est résolu par substitution (cf. exercice 1, question 3). On obtient finalement x=2 411 13 5
Question 1c.
En reproduisant la technique mise en place dans les questions 1a, 1b, on voit que le systèmeAx=best équivalent au système
suivant : A 1x=b1 3 avec A 1=242 4 6
0 0 70 0 53
5 b1=2 4127 53
5 oub1=2 42
4 53
5
Bien sûr, la matriceA1n"est pas inversible (ce qui signifie d"ailleurs aussi que la matriceAn"est pas inversible) . La deuxième
et la troisième lignes sont liées. Dans le premier cas (b1= (12;7;5)T), le second membre est compatible (i.eb1est orthogonal
à Ker(AT1)). Il y a donc une infinité de solutions données par x=2 411 13 5 +2 42
1 03 5
On remarquera que Ker(A1) =Ker(A) =span(2;1;0)T.
Dans le second cas, (b1= (2;4;5)T), le second membre n"est pas compatible (i.eb1n"est pas orthogonal à Ker(AT1)), il n"y
a pas de solution.Question 2.
4 FonctionElimination(A,b) :U,c[Aest une matrice carrée (inversible),betxsont des vecteurs colonne, Uest un matrice carrée triangulaire supérieure,cest un vecteur colonne.Ce programme transforme un systèmeAx=ben un système équivalentUx=c][Quelques initialisations]" 107[choix de la precision"pour le pivot]n length(b)U A;[On initialise la matriceUpar la matriceA]
c b;[On initialise le vecteurc][On vérifie la compatibilité des dimensions deAetb:]Si(size(A,1)6=nou size(A,2)6=n)AlorsAfficher : "problème dans les tailles des matricesAoub"RetournerU,c;
Fin Si
[Processus d"élimination de Gauss] Pouride1àn1faire[1- On cherche un pivot non nul] k=i;Tant que(jU(k;k)j< "ETkn)fairek=k+ 1 Fait indice=k;[2- On échange la ligneiet la ligneindice][Échange la ligneiet la ligne indice dans la matrice (nécessité d introduire une variable temporaire)]
tempU=U(i;:);U(i;:) =U(indice;:);U(indice;:) =tempU;[Échange la ligne i et la ligne indice dans le second membre]
tempc=c(i);c(i) =c(indice);c(indice) =tempc;[3- On élimine l"inconnue i des equations i+1 à n] pivot=U(i;i);[Test si le pivot est trop petit (matrice possiblement non inversible)] Si(jpivotj< ")Alorsdisp("attention pivot très petit : la matrice est elle inversible?")Fin Si
Pourmdei+1ànfairecoef=U(m;i)=pivot;[Modification de la matrice] U(m;1 :i) = 0;U(m;i+ 1 :n) =U(m;i+ 1 :n)coefU(i;i+ 1 :n);[Modification du second membre] c(m) =c(m)coefc(i);Fin Pour
Fin Pour
Retournerc,U;
FinAlgorithme 1:Algorithme d"élimination de Gauss 5quotesdbs_dbs47.pdfusesText_47[PDF] methode de gauss resolution systeme
[PDF] méthode de gestion du temps pdf
[PDF] methode de horner
[PDF] methode de l'anthropologie
[PDF] méthode de la sécante exercice corrigé
[PDF] méthode de la sécante python
[PDF] methode de la variation de la constant
[PDF] methode de lecture syllabique gratuite
[PDF] méthode de lecture syllabique pour apprendre ? lire pas ? pas pdf
[PDF] methode de maintenance pdf
[PDF] Méthode de Mémoire
[PDF] Méthode de Newton
[PDF] methode de newton analyse numerique exercices corrigés
[PDF] méthode de point fixe exercices corrigés pdf