Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Solution-TD1 Méthode Numérique

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 3

Solution Travaux Dirigés n°1:

Les méthodes de résolution directes des systèmes d’équations linéaires

Exercice n° 01

2 1 2 10 2 1 2 10  =
è 01 ∶ 6 4 0 = 26 A 
(0)
 6 4 0 26
8 5 1 35 8 5 1 35 =
..
..

On prend k=1 on annule le terme a21 sous le pivot a11

i=2 j=1 = = 0

i=2 j=2 = = 4 1 = +1

i=2 j=3 = = 0 2 = -6
ici on pose la colonne des bi comme la colonne n+1 de a soit bi= ,

On a alors i=2 j=4 = = 26 10 = -4


On annule ensuite le terme a31 sous le pivot a11
i=3 j=1 = = 0
2 1 2 10 
i=3 j=2 = = 5 1 = +1 A 
(1) 
 0 1  6  4
0  1  7  5
i=3 j=3 = = 1 2 = -7

i=3 j=4 = = 35 10 = -5
2 1 2 10 
On a terminé la 1° transformation sous le pivot 1 on annule maintenant les termes sous
le pivot a22. On prend alors k=2 on aura a32 =0
A 
( 2)
 0 1  6  4
0 0  1  1 
i=k+1=3 j=k+1=3 = = 7 ( 6) = -1

i=3 j=4 = = 5 ( 4) = -1 La solution {x1, x2, x3} = {+3, +2, +1}


2 1 2 +3 10
Alors on trouve à partir de la 3° ligne x3 = (-1)/(-1) = +1 6 4 0 +2 = 26
La 2° ligne nous donne x2= (a24 -a23x3)/a22 = (-4-(-6)*(1))/1 = +2
8 5 1 +1 35
La 1° ligne nous donne x1= (a14 -a13x3- a12x2)/a11 = (10-2*1-2*2)/2 = 3
Système n° 02
3 2 + =2 3 2 1 2
2 + + =7 2 1 1 = 7
4 3 +2 =4 4 3 2 4
On prend k=1 on annule le terme a21 sous le pivot a11
i=2 j=2 = = 1 ( 2) = +7/3

i=2 j=3 = = 1 1 = 1/3


 
ici on pose la colonne des bi comme la colonne n+1 de a soit bi= , 3  2 1 2 
On a alors i=2 j=4 = = 7 2 = +17/3  7 1 17 
On annule ensuite le terme a31 sous le pivot a11
A 
(1)
 0
3 3 3


i=3 j=2 = = 3 ( 2) = -1/3 0  1 2 4 
 3 3 3 
i=3 j=3 = = 2 1 = +2/3

i=3 j=4 = = 4 2 = 4/3

2/4
On a terminé la 1° transformation sous le pivot 1 on annule maintenant les termes sous le
 
pivot a22 3  2 1 2
On prend alors k=2 on aura a32 =0  7 1 17 
i=k+1=3 j=k+1=3 = = + ( )=
A 
( 2)
 0
3 3 3


0 0 5 15 
i=3 j=4 = = + ( )=  7 7 
Alors on trouve
à partir de la 3° ligne x3 = (15/7)/(5/7) = +3
La 2° ligne nous donne x2= (a24 -a23x3)/a22 = ( -( )*(3))/ ( ) = +2
La 1° ligne nous donne x1= (a14 -a13x3- a12x2)/a11 = (2-1*3+2*2)/2 = 1
La solution {x1, x2, x3} = {+1, +2, +3}
3 2 1 +1 2
2 1 1 +2 = 7
4 3 2 +3 4
Exercice n° 02
+3 +3 = 2 1 3 3 2
2 + +5 = 7 2 1 5 = 7
3 + 2 + 6 = 12 3 2 6 12
La stratégie du pivot total consiste à recherche les valeurs maximales de la matrice A est les positionnée en pivot
On peut alors permuter l’équation 1 et 3 pour avoir a11 = 3 > 2 >1
On doit pour permuter les lignes (équations) aussi permuter b3 avec b1 on obtient alors
3 2 6 12
2 1 5 = 7
1 3 3 2
Puis on permute la colonne 1 et 3 pour avoir a11=6 = max(aij i=1,3 et j=1,3)
On doit pour permuter les colonnes aussi permuter l’inconnu x3 avec x1 on obtient alors
6 2 3 12
5 1 2 = 7 Nouveau système equivalent à résoudre
3 3 1 2
On annule la 1° colonne sous le pivot avec k=1 on obtient [A(1)]{X}={B(1)}
6 2 3 12
0 2/3 1/2 = 3
0 2 1/2 8
Puis on cherche la valeur max di pivot a22 sur les lignes 2 et 3. C’est la valeur a32 donc on permute les lignes 2 et 3
6 2 3 12
0 2 1/2 = 8
0 2/3 1/2 3
On annule maintenant la colonne 2 sous le pivot a22 on obtient
6 2 3 12
0 2 1/2 = 8
0 0 2/3 17/3
x1 = 17/2 x2= (-8+1/2*(17/2))/2 = -15/8 x3 = (12-3*(17/2)-2*(-15/8))/6 =-39/24
Vérification
1 3 3 17/2 2
2 1 5 15/8 = 7 Vérifiée donc la solution {x1, x2, x3} = {17/2 ; -15/8 ; -13/8}
3 2 6 13/8 12

3/4
Exercice N° 3
Soient les systèmes linéaires suivants :

3 2 + =2 =
2 + + =7 ( ) ,
4 3 +2 =4 = ∑
,

On calcul les matrices [L] et [U] en faisant varier i de 2 à n et j à chaque fois de i à n ici n=3
Posons d’abord lii = 1 pour i=1 à 3

=
,

=
,

i=1 j=1 u11 = a11/l11 = 3 l11=1


i=1 j=2 u12 = a12= -2 l21=a21/u11=2/3
i=1 j=3 u13 = a13= +1 l31=a31/u11=4/3
i=2 j=2 u22=a22 –l21u12 =1-(2/3)*(-2)=7/3 l22=1
i=2 j=3 u23=a23 –l21u13 =1-2/3*(+1)=+1/3 l32= (a32 –l31u12)/u22= (-3-(4/3)*(-2))/(7/3)=-1/7
i=3 j=3 u33=a33–l31u13 –l32u23 =2-4/3*(+1)-(-1/7)*(1/3)=15/21=5/7 l33=1

      
1 0 0 3  2 1 1 0 0 3  2 1  3  2 1 
   1   7 1
L   2 1 0 U   0 7  A   2 1 0 0  
2 1 1 
 
3   3 3 3  3 3
4  3 2
4 1  0 0 5 4 1  
1 0 0
5
1  3 7   7 
 3 7   7 

  z1 =2/1 = 2
1 0 0  z1  2
     z2 =7-(2/3)*2 =17/3
LZ    2 1 0  z 2   7 
3     z3 = 4-(4/3)*2-(17/3)*(-1/7) =15/7
4  1   z 3  4
1 3  2 1 1 2
 3 7     
AX   2 1 1 2  7 vérifié
4  3 2 3 4
  x3= (15/7) *(7/5) = 3
3  2 1   x1   z1 
 1     x2=(17/3-(1/3)*3) *(3/7)= 2
U X   0 7   x2    z 2 
 3 3     x1= (2+(2)*(2)-1*(3))/3=1
0 0 5   x3   z3 
 7 

4/4

Vous aimerez peut-être aussi