Rapport TP3
Exercice 1
1- A=[2,1,-1;0,-1,2;0,0,-2];
( )
21−1
A= 0−1 2
0 0−2
B=[9,5,6;0,-14,8;0,0,3];
( )
956
B= 0−14 8
003
Pour créer les matrices C et D on utilise la fonction transpose :
C=transpose(A);
( )
200
C= 1−1 0
−12−2
D=transpose(B);
( )
900
D= 5−14 0
683
2- E=A+B;
( )
116 5
E= 0−15 10
001
F=A*B;
( )
18−4 17
F= 0 14−2
0 0−6
3- G=C+D;
( )
110 0
G= 0−15 0
5 10 1
H=C*D;
( )
18 0 0
H= 4 14 0
−11−44−6
4- Les matrices A,B,C et D sont triangulaires : A et B triangulaires supérieures, quant à C et D qui sont
triangulaires inférieures
5- Les matrices E,F,G et H sont triangulaires : E et F triangulaires supérieures, quant à G et H qui sont
triangulaires inférieures
6- Les matrices triangulaires conservent la disposition des zéros avec la multiplication et l’addition
donc en sommant/multipliant deux matrices triangulaires supérieures/inférieures , le résultat est une
matrice triangulaire supérieure/inférieure.
Exercice 2
1- Le système est cohérent, et on a :
{ x+ 2 y =3
4 x +5 y=1 } 1 2
, comme ≠ ≠3 et 1 ×5−4 × 2=−3 ≠ 0 donc le système possède une seule solution.
4 5
2- Le système est cohérent, et on a :
{ x −3 y=5
−2 x +6 y =−10 }
, comme
−1 −3 −5
2
= =
6 10
donc le système possède une infinité de solutions.
3- Le système est non cohérent, et on a :
{−2x−3x +6y=5
y =7 }
, comme
−1 −3 5
2
=
6 7
≠ donc le système ne possède aucune solution.
{
3 x 1 +3 x 2−2 x 3=2
3 −3 −2 −1 2
4- x2 +6 x 3=5 on a : = = = ≠ donc le système ne possède aucune
−9 9 6 3 4
−9 x 1−9 x 2+6 x 3=4
solution.
Exercice 3
3 5 7 101
1- La matrice M augmentée est : 2 10 6 134
1 2 3 40
3-
%%Formulation matricielle de la méthode de Gauss
M=[3,5,7,101;2,10,6,134;1,2,3,40];
%étape 1
M(1,:)=M(1,:)-M(2,:)
M(2,:)=M(2,:)-2*M(3,:)
M(1,:)=M(3,:)
M(3,:)=[1,-5,1,-33]
%étape 2
M(1,:)=M(1,:)-1/3*M(2,:)
M(2,:)=1/6*M(2,:)
M(3,:)=M(3,:)-M(1,:)
M(3,:)=M(3,:)+5*M(2,:)
%étape 3
M(3,:)=-1/2*M(3,:)
M(1,:)=M(1,:)-3*M(3,:)
1 0 07
M=0 1 0 9
0 0 15
La solution du système est : (7,9,5)
Exercice 4
1 1 1 x1 1
1- Le système est : 0 1 2 × x 2=1
1 2 4 x3 1
1 1 1 x1 1
Avec A= 0 1 2, x= x 2 et b= 1
1 2 4 x3 1
2-
%%%Méthode de Gauss-Jordan
A=[1,1,1,;0,1,2;1,2,1];
I=[1,0,0;0,1,0;0,0,1];
AI=[A I]
%étape 1
AI(3,:)=AI(3,:)-AI(1,:)
%étape 2
AI(1,:)=AI(1,:)-AI(2,:)
AI(3,:)=AI(3,:)-AI(2,:)
%étape 3
AI(1,:)=AI(1,:)-1/2*AI(3,:)
AI(2,:)=AI(2,:)+AI(3,:)
AI(3,:)=-1/2*AI(3,:)
1 0 0 3/2 −1 /2 −1/2
[ AI ] =0 1 0 −1 0 1
0 0 1 1/2 1 /2 −1/2
3- [ I inv ( A ) ] =[ AI ]
4-
%%%Méthode de Gauss-Jordan
a=[1,1,1,1;0,1,2,1;1,2,4,1]
%étape 1
a(3,:)=a(3,:)-a(1,:)
%étape 2
a(1,:)=a(1,:)-a(2,:)
a(3,:)=a(3,:)-a(2,:)
%étape 2
a(1,:)=a(1,:)+a(3,:)
a(2,:)=a(2,:)-2*a(3,:)
1 0 0 −1
a=0 1 0 3
0 0 1 −1
Donc l’intersection des trois plans est (-1,3,-1-)
Exercice 5
1-
Si a ≠ 1 le système n’a qu’une solution.
Si a=b=1 le système a une infinité de solutions.
Si a=1 et b ≠ 1 le système n’a aucune solution.
2-
Dans le cas où il n’y a qu’une solution :
ab−1 1−b a−ab
S={x 1= , x 2= , x 3= /a ∈ R−{ 1 } et b ∈ R }
1−a 1−a 1−a
3- Dans le cas où le système admet une infinité de solutions :
S={x 1 ∈ R , x 2=x 3 }
Exercice 6
1-
%%%%Décomposition LU
A=[2,3;4,1];
B=[5];
C=[2,1,-1;0,-1,2;-1,0,-2];
M=blkdiag(A,B,C)
I=[1,0,0,0,0,0;0,1,0,0,0,0;0,0,1,0,0,0;0,0,0,1,0,0;0,0,0,0,1,0;0,0,0,0,0,1];
a=[M I]
%étape 1
a(1,:)=1/2*a(1,:)
a(2,:)=a(2,:)-4*a(1,:)
%étape 2
a(1,:)=a(1,:)+3/10*a(2,:)
a(2,:)=-1/5*a(2,:)
%étape 3
a(3,:)=1/5*a(3,:)
%étape 4
a(4,:)=1/2*a(4,:)
a(6,:)=a(6,:)+a(4,:)
%étape 5
a(4,:)=a(4,:)+1/2*a(5,:)
a(5,:)=-1*a(5,:)
a(6,:)=a(6,:)-1/2*a(5,:)
%étape 6
a(4,:)=a(4,:)+1/3*a(6,:)
a(5,:)=a(5,:)-4/3*a(6,:)
a(6,:)=-2/3*a(6,:)
−0 ,1 0 ,3 0 0 0 0
−1
M = 0 , 4 −0 , 2 0 0 0 0
0 0 0,20 0 0
0 0 0 −0,6667 −0,666 7 0,3333
0 0 0 −0,6667 −0,666 7 −1,3333
0 0 0 −0,3333 0,3333 −0,6667
2-
%%%%Décomposition LU
A=[2,3;4,1];
B=[5];
C=[2,1,-1;0,-1,2;-1,0,-2];
M=blkdiag(A,B,C)
M
%ligne 1
M(2,:)=M(2,:)-2*M(1,:)
%ligne 6
M(6,:)=M(6,:)+1/2*M(4,:)
M(6,:)=M(6,:)+1/2*M(5,
L=[1,0,0,0,0,0;2,1,0,0,0,0;0,0,1,0,0,0;0,0,0,1,0,0;0,0,0,0,1,0;0,0,0,-1/2,-1/2,1];
U=M
L*U;
2 3 00 0 0
U =0 −5 0 0 0 0
0 0 50 0 0
0 0 02 1 −1
0 0 0 0 −1 2
0 0 00 0 −1 ,5
1 0 00 0 0
L= 2 2 00 0 0
0 0 50 0 0
0 0 1 1 0 0
0 0 0 0 1 0
0 0 0 −0 ,5 −0 , 5 1
L∗U =M
3- On a L∗U =M et L∗D∗U + M donc D est la matrice identité
