ECG appliquées · Chapitre 02 · Premier semestre

Devoir surveillé — Calcul matriciel et résolution de systèmes linéaires

Sujet type, 210 min, barème sur 20 points. À faire en conditions réelles avant de regarder le corrigé.

Sujet type DS — 210 min, barème sur 20 points. Faites-le en conditions réelles avant de regarder le corrigé (PDF).

Exercice 1 (2,5 points) — Cinq affirmations à trancher

Consignes valables pour tout le sujet. Les cinq exercices sont indépendants et peuvent être traités dans l'ordre de votre choix. La calculatrice n'est pas autorisée : tous les calculs numériques du sujet se font à la main. La qualité de la rédaction, la clarté des justifications et le soin apporté aux calculs interviennent dans l'appréciation des copies ; une affirmation non justifiée ne rapporte aucun point.

Chacune des cinq affirmations suivantes est soit vraie, soit fausse. Dire dans chaque cas laquelle des deux réponses est la bonne, puis la justifier : par une démonstration si l'affirmation est vraie, par un contre-exemple explicite si elle est fausse. Un contre-exemple doit être entièrement vérifié par le calcul.

1. (0,5 pt) Si A et B sont deux matrices de M2(R) telles que AB=02, alors A=02 ou B=02.

2. (0,5 pt) Si A est une matrice de M3(R) telle que A3=03, alors A n'est pas inversible.

3. (0,5 pt) Si A est une matrice inversible de M2(R), alors tA est inversible et (tA)1=t(A1).

4. (0,5 pt) Un système linéaire de deux équations à trois inconnues admet toujours une infinité de solutions.

5. (0,5 pt) Pour toute matrice AM3,2(R), le produit AtA est une matrice symétrique d'ordre 3.

Exercice 2 (3,5 points) — Trois tableaux qui s'enchaînent

Un traiteur propose trois menus, notés M1, M2 et M3. Chaque menu est préparé à partir de quatre ingrédients, notés I1, I2, I3 et I4. Le tableau des recettes donne, pour chaque menu, la quantité de chaque ingrédient nécessaire à la préparation de cent couverts, en kilogrammes :

Pour cent couverts I1 I2 I3 I4
Menu M1 12 5 0 8
Menu M2 10 8 6 4
Menu M3 6 4 10 12

Les quatre ingrédients peuvent être achetés chez deux fournisseurs, F1 et F2. Le tableau des prix donne le prix d'un kilogramme de chaque ingrédient chez chaque fournisseur, en euros :

Prix au kilogramme F1 F2
I1 2 5
I2 5 4
I3 4 2
I4 6 5

Le traiteur doit servir deux réceptions, notées E1 et E2. Le tableau des commandes donne, pour chaque réception, le nombre de centaines de couverts de chaque menu :

Centaines de couverts M1 M2 M3
Réception E1 3 2 1
Réception E2 1 2 5

On note R, P et N les trois matrices obtenues en recopiant ces tableaux, dans le même ordre de lignes et de colonnes, sans les en-têtes.

1. (0,5 pt) Donner le format de R, de P et de N. Parmi les cinq produits RP, PR, NR, RN et NP, dire lesquels sont définis et donner le format de ceux qui le sont.

2. (0,75 pt) Calculer la matrice RP. Interpréter en une phrase, dans le contexte du traiteur, le coefficient situé à la ligne 3 et à la colonne 2.

3. (0,75 pt) Calculer la matrice NR. Interpréter en une phrase le coefficient situé à la ligne 2 et à la colonne 3.

4. (0,75 pt) Calculer le produit N×(RP) à partir du résultat de la question 2, puis retrouver son coefficient d'indice (1,1) en calculant à la place le coefficient d'indice (1,1) de (NR)×P. Quelle propriété du produit matriciel garantit que ces deux façons de facturer donnent le même résultat ?

5. (0,25 pt) Le traiteur commande la totalité des ingrédients d'une réception chez un seul fournisseur, mais il peut changer de fournisseur d'une réception à l'autre. Quel fournisseur doit-il choisir pour chacune des deux réceptions, et quelle est alors la dépense totale ?

6. (0,5 pt) Le service comptable range ses tableaux dans l'autre sens : il écrit les ingrédients en lignes dans le tableau des recettes, et les fournisseurs en lignes dans le tableau des prix. Il travaille donc avec tR, tP et tN. Exprimer t(NRP) à l'aide de ces trois matrices, en justifiant l'ordre des facteurs, puis vérifier la formule obtenue sur le coefficient d'indice (2,1).

Exercice 3 (3 points) — Le déterminant d'un produit

Dans tout l'exercice, les matrices considérées appartiennent à M2(R). On rappelle que le déterminant de M=(abcd) est le réel det(M)=adbc, et que M est inversible si et seulement si det(M)0.

1. (1 pt) Soient A=(abcd) et B=(abcd) deux matrices quelconques de M2(R). Calculer le produit AB, puis démontrer par le calcul que

det(AB)=det(A)×det(B).

2. (0,5 pt) En déduire que si A est inversible, alors det(A1)=1det(A).

3. (0,5 pt) Démontrer par récurrence que det(An)=(det(A))n pour tout entier n1.

4. (0,5 pt) On pose Q=(3152) et S=(2346). Sans calculer aucune puissance, dire si Q7 est inversible, et démontrer qu'aucune puissance Sn, avec n1, ne l'est.

5. (0,5 pt) Soient A et B deux matrices de M2(R) telles que le produit AB soit inversible. Démontrer que A et B sont toutes les deux inversibles.

Exercice 4 (4 points) — Le chargement d'un camion

Un producteur de granulés de bois conditionne sa production dans des sacs de trois formats : des sacs de 200 kilogrammes, des sacs de 300 kilogrammes et des sacs de 500 kilogrammes. Il charge un camion et doit respecter deux contraintes en même temps :

  • le camion emporte exactement 30 sacs, quel que soit leur format, car il ne dispose que de 30 emplacements ;
  • la masse totale chargée doit être exactement de 11 tonnes, soit 11000 kilogrammes.

On note x, y et z les nombres de sacs de 200, 300 et 500 kilogrammes chargés.

1. (0,5 pt) Traduire les deux contraintes par un système linéaire (S) d'inconnue (x,y,z), puis écrire (S) sous la forme matricielle AX=B en précisant A, X, B et leurs formats. Expliquer en une phrase pourquoi la question « la matrice A est-elle inversible ? » n'a ici aucun sens.

2. (1 pt) Résoudre (S) dans R3 par la méthode du pivot de Gauss, en écrivant l'opération élémentaire utilisée. Décrire l'ensemble S des solutions à l'aide d'un paramètre réel, puis vérifier le résultat dans les deux équations de départ.

3. (0,75 pt) Un nombre de sacs est un entier naturel. Déterminer tous les chargements possibles et dire combien il y en a.

4. (0,75 pt) Le client impose une contrainte supplémentaire : le camion doit contenir au moins deux fois plus de sacs de 500 kilogrammes que de sacs de 200 kilogrammes. Combien de chargements restent-ils possibles ? Les donner tous.

5. (0,5 pt) La marge du producteur est de 400 euros par sac de 200 kilogrammes, 250 euros par sac de 300 kilogrammes et 600 euros par sac de 500 kilogrammes. Exprimer la marge totale en fonction du seul paramètre choisi à la question 2, puis déterminer, parmi les chargements de la question 4, celui qui rapporte le plus.

6. (0,5 pt) Un commercial affirme : « le système (S) a une infinité de solutions, donc le producteur dispose d'une infinité de chargements possibles ». Expliquer précisément où est l'erreur.

Exercice 5 (7 points) — Problème : une matrice qui se coupe en deux

Le pivot de Gauss transforme une matrice carrée en une matrice triangulaire supérieure. Ce problème montre que cette transformation garde la trace de tout ce qu'elle a fait, sous la forme d'une seconde matrice triangulaire, et que le produit des deux redonne la matrice de départ. On obtient ainsi une écriture A=LU, où L est triangulaire inférieure et U triangulaire supérieure, qui permet ensuite de résoudre très vite le système AX=B pour autant de seconds membres B que l'on veut.

Pour tout triplet de réels (a,b,c), on note

L(a,b,c)=(100a10bc1)M3(R).

Ce sont exactement les matrices triangulaires inférieures d'ordre 3 dont les trois coefficients diagonaux valent 1. On remarquera en particulier que L(0,0,0)=I3.

Les trois parties peuvent être traitées indépendamment les unes des autres, à condition d'admettre les résultats des parties précédentes.

Partie A — les matrices L(a,b,c).

1. (0,75 pt) Soient (a,b,c) et (a,b,c) deux triplets de réels. Calculer le produit L(a,b,c)×L(a,b,c) et vérifier qu'il est encore de la forme L(a,b,c), en donnant a, b et c en fonction de a, b, c, a, b et c.

2. (0,75 pt) En déduire que L(a,b,c) est inversible et que son inverse est encore de cette forme, en précisant lequel. On vérifiera les deux produits.

3. (0,5 pt) Démontrer que la seule matrice qui soit à la fois de la forme L(a,b,c) et triangulaire supérieure est I3.

Partie B — la décomposition d'une matrice. On considère désormais

A=(2114502811)M3(R).

4. (1 pt) Échelonner A par la méthode du pivot de Gauss, sans jamais échanger deux lignes, en écrivant chaque opération sous la forme LiLimLj. On note U la matrice triangulaire supérieure obtenue, et m2,1, m3,1, m3,2 les trois multiplicateurs utilisés, dans cet ordre.

5. (0,75 pt) Calculer le produit L ⁣(m2,1,m3,1,m3,2)×U et constater le résultat.

6. (0,75 pt) Retrouver ce résultat par identification, sans utiliser le pivot : on cherche des réels a, b, c et une matrice triangulaire supérieure V=(pqr0su00w) tels que L(a,b,c)V=A. Écrire les neuf équations obtenues et les résoudre dans un ordre bien choisi, en expliquant pourquoi cet ordre s'impose.

Partie C — à quoi cette écriture sert. On pose L=L ⁣(m2,1,m3,1,m3,2), de sorte que A=LU.

7. (0,75 pt) On veut résoudre AX=B1, où B1=(11447). Justifier que ce système équivaut à la résolution successive de LY=B1, d'inconnue la colonne Y, puis de UX=Y, d'inconnue X. Effectuer ces deux résolutions et vérifier la solution obtenue dans le système de départ.

8. (0,5 pt) Reprendre la question 7 avec le second membre B2=(231).

9. (0,25 pt) Expliquer en deux ou trois phrases l'intérêt de cette méthode lorsqu'on doit résoudre AX=B pour un grand nombre de seconds membres B différents.

10. (0,75 pt) Démontrer l'unicité de l'écriture : si A=L(a,b,c)U1=L(a,b,c)U2U1 et U2 sont deux matrices triangulaires supérieures inversibles, alors L(a,b,c)=L(a,b,c) et U1=U2. On admettra les deux résultats du cours suivants : le produit de deux matrices triangulaires supérieures est triangulaire supérieure, et l'inverse d'une matrice triangulaire supérieure inversible est triangulaire supérieure.

11. (0,25 pt) Une telle écriture n'existe pas toujours. En raisonnant sur les coefficients, démontrer qu'il n'existe aucun réel a ni aucune matrice triangulaire supérieure V de M2(R) tels que

(0110)=(10a1)V.

Bloqué sur « Calcul matriciel et résolution de systèmes linéaires » ?

On peut le travailler ensemble dès cette semaine. La première heure est offerte — on fait le point honnêtement, et vous repartez au minimum avec une méthode.