1. Vocabulaire
↩ Revoir le cours
Exercice 1 — ⭐
Écris la matrice A A A de taille 2 × 3 2 \times 3 2 × 3 définie par a i j = i + 2 j a_{ij} = i + 2j a ij = i + 2 j .
Écris la matrice carrée B B B d’ordre 3 3 3 définie par b i j = 1 b_{ij} = 1 b ij = 1 si i ⩽ j i \leqslant j i ⩽ j et b i j = 0 b_{ij} = 0 b ij = 0 sinon.
Détermine les réels x x x et y y y tels que ( x + y 2 1 x − y ) = ( 5 2 1 1 ) \begin{pmatrix} x + y & 2 \\ 1 & x - y \end{pmatrix} = \begin{pmatrix} 5 & 2 \\ 1 & 1 \end{pmatrix} ( x + y 1 2 x − y ) = ( 5 1 2 1 ) .
Voir le corrigé
A = ( 3 5 7 4 6 8 ) A = \begin{pmatrix} 3 & 5 & 7 \\ 4 & 6 & 8 \end{pmatrix} A = ( 3 4 5 6 7 8 ) .
B = ( 1 1 1 0 1 1 0 0 1 ) B = \begin{pmatrix} 1 & 1 & 1 \\ 0 & 1 & 1 \\ 0 & 0 & 1 \end{pmatrix} B = 1 0 0 1 1 0 1 1 1 (matrice « triangulaire supérieure »).
x + y = 5 x + y = 5 x + y = 5 et x − y = 1 x - y = 1 x − y = 1 : x = 3 x = 3 x = 3 , y = 2 y = 2 y = 2 .
2. Opérations sur les matrices
↩ Revoir le cours
Exercice 2 — ⭐⭐
On pose A = ( 1 − 2 3 0 ) A = \begin{pmatrix} 1 & -2 \\ 3 & 0 \end{pmatrix} A = ( 1 3 − 2 0 ) , B = ( 2 1 − 1 4 ) B = \begin{pmatrix} 2 & 1 \\ -1 & 4 \end{pmatrix} B = ( 2 − 1 1 4 ) et C = ( 1 0 2 − 1 3 1 ) C = \begin{pmatrix} 1 & 0 & 2 \\ -1 & 3 & 1 \end{pmatrix} C = ( 1 − 1 0 3 2 1 ) .
Calcule A + B A + B A + B , 3 A − 2 B 3A - 2B 3 A − 2 B , A B AB A B et B A BA B A . Que remarques-tu ?
Le produit A C AC A C est-il défini ? Et C A CA C A ? Calcule celui qui existe.
Voir le corrigé
A + B = ( 3 − 1 2 4 ) A + B = \begin{pmatrix} 3 & -1 \\ 2 & 4 \end{pmatrix} A + B = ( 3 2 − 1 4 ) ; 3 A − 2 B = ( − 1 − 8 11 − 8 ) 3A - 2B = \begin{pmatrix} -1 & -8 \\ 11 & -8 \end{pmatrix} 3 A − 2 B = ( − 1 11 − 8 − 8 ) ; A B = ( 4 − 7 6 3 ) AB = \begin{pmatrix} 4 & -7 \\ 6 & 3 \end{pmatrix} A B = ( 4 6 − 7 3 ) ; B A = ( 5 − 4 11 2 ) BA = \begin{pmatrix} 5 & -4 \\ 11 & 2 \end{pmatrix} B A = ( 5 11 − 4 2 ) .
A B ≠ B A AB \neq BA A B = B A : le produit n’est pas commutatif.
A A A est 2 × 2 2 \times 2 2 × 2 , C C C est 2 × 3 2 \times 3 2 × 3 : A C AC A C existe (2 × 3 2 \times 3 2 × 3 ), C A CA C A non (3 colonnes pour C C C , 2 lignes pour A A A ).
A C = ( 3 − 6 0 3 0 6 ) AC = \begin{pmatrix} 3 & -6 & 0 \\ 3 & 0 & 6 \end{pmatrix} A C = ( 3 3 − 6 0 0 6 ) .
3. Puissances d’une matrice carrée
↩ Revoir le cours
Exercice 3 — ⭐⭐
Calcule D 5 D^5 D 5 pour D = ( 2 0 0 − 1 ) D = \begin{pmatrix} 2 & 0 \\ 0 & -1 \end{pmatrix} D = ( 2 0 0 − 1 ) .
Soit N = ( 0 1 1 0 0 1 0 0 0 ) N = \begin{pmatrix} 0 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix} N = 0 0 0 1 0 0 1 1 0 . Calcule N 2 N^2 N 2 et N 3 N^3 N 3 . Que vaut N n N^n N n pour n ⩾ 3 n \geqslant 3 n ⩾ 3 ?
Soit A = ( 3 0 1 1 ) A = \begin{pmatrix} 3 & 0 \\ 1 & 1 \end{pmatrix} A = ( 3 1 0 1 ) . Montre par récurrence que A n = ( 3 n 0 3 n − 1 2 1 ) A^n = \begin{pmatrix} 3^n & 0 \\ \frac{3^n - 1}{2} & 1 \end{pmatrix} A n = ( 3 n 2 3 n − 1 0 1 ) .
Voir le corrigé
D 5 = ( 32 0 0 − 1 ) D^5 = \begin{pmatrix} 32 & 0 \\ 0 & -1 \end{pmatrix} D 5 = ( 32 0 0 − 1 ) .
N 2 = ( 0 0 1 0 0 0 0 0 0 ) N^2 = \begin{pmatrix} 0 & 0 & 1 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix} N 2 = 0 0 0 0 0 0 1 0 0 , N 3 = 0 N^3 = 0 N 3 = 0 , donc N n = N 3 N n − 3 = 0 N^n = N^3N^{n-3} = 0 N n = N 3 N n − 3 = 0 pour n ⩾ 3 n \geqslant 3 n ⩾ 3 .
Pour n = 0 n = 0 n = 0 : ( 1 0 0 1 ) = I \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} = I ( 1 0 0 1 ) = I ✓. Hérédité : A n + 1 = A n A = ( 3 n 0 3 n − 1 2 1 ) ( 3 0 1 1 ) = ( 3 n + 1 0 3 ( 3 n − 1 ) 2 + 1 1 ) A^{n+1} = A^nA = \begin{pmatrix} 3^n & 0 \\ \frac{3^n - 1}{2} & 1 \end{pmatrix}\begin{pmatrix} 3 & 0 \\ 1 & 1 \end{pmatrix} = \begin{pmatrix} 3^{n+1} & 0 \\ \frac{3(3^n - 1)}{2} + 1 & 1 \end{pmatrix} A n + 1 = A n A = ( 3 n 2 3 n − 1 0 1 ) ( 3 1 0 1 ) = ( 3 n + 1 2 3 ( 3 n − 1 ) + 1 0 1 )
et 3 n + 1 − 3 2 + 1 = 3 n + 1 − 1 2 \dfrac{3^{n+1} - 3}{2} + 1 = \dfrac{3^{n+1} - 1}{2} 2 3 n + 1 − 3 + 1 = 2 3 n + 1 − 1 : la formule est vraie au rang n + 1 n + 1 n + 1 .
Exercice 4 — ⭐⭐⭐
On pose A = ( 4 − 2 1 1 ) A = \begin{pmatrix} 4 & -2 \\ 1 & 1 \end{pmatrix} A = ( 4 1 − 2 1 ) , P = ( 1 2 1 1 ) P = \begin{pmatrix} 1 & 2 \\ 1 & 1 \end{pmatrix} P = ( 1 1 2 1 ) et D = ( 2 0 0 3 ) D = \begin{pmatrix} 2 & 0 \\ 0 & 3 \end{pmatrix} D = ( 2 0 0 3 ) .
Calcule P − 1 P^{-1} P − 1 et vérifie que A = P D P − 1 A = PDP^{-1} A = P D P − 1 .
Montre par récurrence que A n = P D n P − 1 A^n = PD^nP^{-1} A n = P D n P − 1 pour tout entier n ⩾ 1 n \geqslant 1 n ⩾ 1 .
En déduire les coefficients de A n A^n A n en fonction de n n n .
Voir le corrigé
det P = 1 − 2 = − 1 \det P = 1 - 2 = -1 det P = 1 − 2 = − 1 , P − 1 = − ( 1 − 2 − 1 1 ) = ( − 1 2 1 − 1 ) P^{-1} = -\begin{pmatrix} 1 & -2 \\ -1 & 1 \end{pmatrix} = \begin{pmatrix} -1 & 2 \\ 1 & -1 \end{pmatrix} P − 1 = − ( 1 − 1 − 2 1 ) = ( − 1 1 2 − 1 ) . P D = ( 2 6 2 3 ) PD = \begin{pmatrix} 2 & 6 \\ 2 & 3 \end{pmatrix} P D = ( 2 2 6 3 ) et
P D P − 1 = ( − 2 + 6 4 − 6 − 2 + 3 4 − 3 ) = ( 4 − 2 1 1 ) = A PDP^{-1} = \begin{pmatrix} -2 + 6 & 4 - 6 \\ -2 + 3 & 4 - 3 \end{pmatrix} = \begin{pmatrix} 4 & -2 \\ 1 & 1 \end{pmatrix} = A P D P − 1 = ( − 2 + 6 − 2 + 3 4 − 6 4 − 3 ) = ( 4 1 − 2 1 ) = A .
Vrai pour n = 1 n = 1 n = 1 . Si A n = P D n P − 1 A^n = PD^nP^{-1} A n = P D n P − 1 , alors A n + 1 = A n A = P D n P − 1 P D P − 1 = P D n + 1 P − 1 A^{n+1} = A^nA = PD^nP^{-1}PDP^{-1} = PD^{n+1}P^{-1} A n + 1 = A n A = P D n P − 1 P D P − 1 = P D n + 1 P − 1 .
P D n = ( 2 n 2 × 3 n 2 n 3 n ) PD^n = \begin{pmatrix} 2^n & 2 \times 3^n \\ 2^n & 3^n \end{pmatrix} P D n = ( 2 n 2 n 2 × 3 n 3 n ) , puis A n = ( − 2 n + 2 × 3 n 2 n + 1 − 2 × 3 n − 2 n + 3 n 2 n + 1 − 3 n ) A^n = \begin{pmatrix} -2^n + 2 \times 3^n & 2^{n+1} - 2 \times 3^n \\ -2^n + 3^n & 2^{n+1} - 3^n \end{pmatrix} A n = ( − 2 n + 2 × 3 n − 2 n + 3 n 2 n + 1 − 2 × 3 n 2 n + 1 − 3 n ) .
4. Inverse d’une matrice et systèmes linéaires
↩ Revoir le cours
Exercice 5 — ⭐⭐
Les matrices ( 3 2 6 4 ) \begin{pmatrix} 3 & 2 \\ 6 & 4 \end{pmatrix} ( 3 6 2 4 ) et ( 5 3 3 2 ) \begin{pmatrix} 5 & 3 \\ 3 & 2 \end{pmatrix} ( 5 3 3 2 ) sont-elles inversibles ? Calcule l’inverse quand il existe.
Résous le système { 5 x + 3 y = 1 3 x + 2 y = 2 \begin{cases} 5x + 3y = 1 \\ 3x + 2y = 2 \end{cases} { 5 x + 3 y = 1 3 x + 2 y = 2 par le calcul matriciel.
Interpolation. On cherche une parabole y = a x 2 + b x + c y = ax^2 + bx + c y = a x 2 + b x + c passant par ( 1 ; 2 ) (1\,;\,2) ( 1 ; 2 ) , ( 2 ; 3 ) (2\,;\,3) ( 2 ; 3 ) et ( 3 ; 6 ) (3\,;\,6) ( 3 ; 6 ) . Écris le système sous la
forme M X = Y MX = Y M X = Y , puis résous-le à la calculatrice ou avec numpy.linalg.solve.
Voir le corrigé
3 × 4 − 2 × 6 = 0 3 \times 4 - 2 \times 6 = 0 3 × 4 − 2 × 6 = 0 : la première n’est pas inversible. 5 × 2 − 3 × 3 = 1 5 \times 2 - 3 \times 3 = 1 5 × 2 − 3 × 3 = 1 : la seconde l’est, d’inverse ( 2 − 3 − 3 5 ) \begin{pmatrix} 2 & -3 \\ -3 & 5 \end{pmatrix} ( 2 − 3 − 3 5 ) .
X = ( 2 − 3 − 3 5 ) ( 1 2 ) = ( − 4 7 ) X = \begin{pmatrix} 2 & -3 \\ -3 & 5 \end{pmatrix}\begin{pmatrix} 1 \\ 2 \end{pmatrix} = \begin{pmatrix} -4 \\ 7 \end{pmatrix} X = ( 2 − 3 − 3 5 ) ( 1 2 ) = ( − 4 7 ) : x = − 4 x = -4 x = − 4 , y = 7 y = 7 y = 7 .
M = ( 1 1 1 4 2 1 9 3 1 ) M = \begin{pmatrix} 1 & 1 & 1 \\ 4 & 2 & 1 \\ 9 & 3 & 1 \end{pmatrix} M = 1 4 9 1 2 3 1 1 1 , X = ( a b c ) X = \begin{pmatrix} a \\ b \\ c \end{pmatrix} X = a b c , Y = ( 2 3 6 ) Y = \begin{pmatrix} 2 \\ 3 \\ 6 \end{pmatrix} Y = 2 3 6 . On trouve a = 1 a = 1 a = 1 , b = − 2 b = -2 b = − 2 , c = 3 c = 3 c = 3 : y = x 2 − 2 x + 3 y = x^2 - 2x + 3 y = x 2 − 2 x + 3 .
↩ Revoir le cours
Exercice 6 — ⭐⭐
On note R R R la matrice de la rotation de centre O O O d’angle π 2 \dfrac\pi2 2 π et S S S celle de la symétrie d’axe ( O x ) (Ox) ( O x ) .
Écris R R R et S S S . Calcule l’image du point M ( 3 ; 1 ) M(3\,;\,1) M ( 3 ; 1 ) par la rotation.
Calcule R 2 R^2 R 2 et R 4 R^4 R 4 . Interprète.
Calcule S R SR S R et R S RS R S . L’ordre dans lequel on enchaîne les transformations compte-t-il ?
Voir le corrigé
R = ( 0 − 1 1 0 ) R = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} R = ( 0 1 − 1 0 ) , S = ( 1 0 0 − 1 ) S = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} S = ( 1 0 0 − 1 ) . R M = ( − 1 3 ) RM = \begin{pmatrix} -1 \\ 3 \end{pmatrix} R M = ( − 1 3 ) : M ′ ( − 1 ; 3 ) M'(-1\,;\,3) M ′ ( − 1 ; 3 ) .
R 2 = ( − 1 0 0 − 1 ) = − I R^2 = \begin{pmatrix} -1 & 0 \\ 0 & -1 \end{pmatrix} = -I R 2 = ( − 1 0 0 − 1 ) = − I (rotation d’angle π \pi π , symétrie de centre O O O ) ; R 4 = I R^4 = I R 4 = I (quatre quarts de tour).
S R = ( 0 − 1 − 1 0 ) SR = \begin{pmatrix} 0 & -1 \\ -1 & 0 \end{pmatrix} S R = ( 0 − 1 − 1 0 ) et R S = ( 0 1 1 0 ) RS = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} R S = ( 0 1 1 0 ) : différents, l’ordre compte.
6. Suites de matrices colonnes
↩ Revoir le cours
Exercice 7 — ⭐⭐⭐
Deux suites sont définies par u 0 = 1 u_0 = 1 u 0 = 1 , v 0 = 0 v_0 = 0 v 0 = 0 et, pour tout n n n : u n + 1 = 0,5 u n + 0,25 v n + 1 u_{n+1} = 0{,}5u_n + 0{,}25v_n + 1 u n + 1 = 0 , 5 u n + 0 , 25 v n + 1 et v n + 1 = 0,25 u n + 0,5 v n + 1 v_{n+1} = 0{,}25u_n + 0{,}5v_n + 1 v n + 1 = 0 , 25 u n + 0 , 5 v n + 1 .
Écris la relation sous la forme U n + 1 = A U n + C U_{n+1} = AU_n + C U n + 1 = A U n + C .
Détermine la matrice colonne S S S telle que S = A S + C S = AS + C S = A S + C .
On pose V n = U n − S V_n = U_n - S V n = U n − S . Montre que V n + 1 = A V n V_{n+1} = AV_n V n + 1 = A V n et en déduire U n U_n U n en fonction de A n A^n A n .
À l’aide de Python, calcule U 30 U_{30} U 30 et conjecture la limite de ( u n ) (u_n) ( u n ) et de ( v n ) (v_n) ( v n ) .
Voir le corrigé
A = ( 0,5 0,25 0,25 0,5 ) A = \begin{pmatrix} 0{,}5 & 0{,}25 \\ 0{,}25 & 0{,}5 \end{pmatrix} A = ( 0 , 5 0 , 25 0 , 25 0 , 5 ) , C = ( 1 1 ) C = \begin{pmatrix} 1 \\ 1 \end{pmatrix} C = ( 1 1 ) .
S = ( s t ) S = \begin{pmatrix} s \\ t \end{pmatrix} S = ( s t ) : s = 0,5 s + 0,25 t + 1 s = 0{,}5s + 0{,}25t + 1 s = 0 , 5 s + 0 , 25 t + 1 et t = 0,25 s + 0,5 t + 1 t = 0{,}25s + 0{,}5t + 1 t = 0 , 25 s + 0 , 5 t + 1 , soit 0,5 s − 0,25 t = 1 0{,}5s - 0{,}25t = 1 0 , 5 s − 0 , 25 t = 1 et − 0,25 s + 0,5 t = 1 -0{,}25s + 0{,}5t = 1 − 0 , 25 s + 0 , 5 t = 1 ; par symétrie
s = t s = t s = t , puis 0,25 s = 1 0{,}25s = 1 0 , 25 s = 1 : S = ( 4 4 ) S = \begin{pmatrix} 4 \\ 4 \end{pmatrix} S = ( 4 4 ) .
V n + 1 = U n + 1 − S = A U n + C − ( A S + C ) = A ( U n − S ) = A V n V_{n+1} = U_{n+1} - S = AU_n + C - (AS + C) = A(U_n - S) = AV_n V n + 1 = U n + 1 − S = A U n + C − ( A S + C ) = A ( U n − S ) = A V n , donc V n = A n V 0 V_n = A^nV_0 V n = A n V 0 et U n = A n ( U 0 − S ) + S U_n = A^n(U_0 - S) + S U n = A n ( U 0 − S ) + S .
import numpy as np
A = np.array([[ 0.5 , 0.25 ], [ 0.25 , 0.5 ]]); C = np.array([ 1 , 1 ]); U = np.array([ 1 , 0 ])
for n in range ( 30 ):
U = A @ U + C
print (U) # [3.99937... 3.99937...]
Les deux suites semblent tendre vers 4 4 4 (les coefficients de A n A^n A n tendent vers 0 0 0 ).