On rappelle que sur l'ensemble des chapitres d'algèbre linéaire, la bibliothèque la plus utilisée est la bibliothèque numpy
et son module linalg
au sein duquel de nombreuses méthodes permettent de répondre aux questions les plus courantes.
On commencera donc par importer cette bibliothèque.
import numpy as np
Remarque : L'ensemble des exercices qui suivent peuvent être résolus avec cette bibliothèque mais il faudra les traiter au préalable avec un papier et un crayon. Ce travail est indispensable car c'est tout ce dont vous disposerez sur les deux épreuves $A$ et $B$, le jour de l'écrit.
L'objectif de ce notebook est essentiellement de vous permette d'acquérir les réflexes qui vous permettront de gagner du temps sur ces calculs. Même s'ils sont parfois répétitifs, il faut les faire !
Remarque : Le cas échéant, si vous n'aboutissez pas, une correction des exercices à partir du paragraphe "2. Espaces vectoriels" est proposée en fin de notebook.
Deux moyens s'offrent à vous pour créer un objet matriciel : A = np.array(L)
ou A = np.matrix(L)
où L
est une liste d'entiers ou de réels.
Nous choisirons ici la deuxième possibilité car elle allège sensiblement l'écriture des produits matriciels (cf. le notebook : "td0BCPST2b-lesPrincipauxModules").
Exercice 1.1 On considère les matrices à coefficients réels :
$A=\begin{pmatrix} 1&2\\3&4\end{pmatrix}$, $B=\begin{pmatrix} -1&0&1\\2&2&3\end{pmatrix}$, $C=\begin{pmatrix} 1&0&1\\-3&-1&2\\5&2&1\end{pmatrix}$, $D=\begin{pmatrix} 1\\-1\\3\end{pmatrix}$, $E=\begin{pmatrix} 1&-1\\-1&2\end{pmatrix}$
B=A.T
ou B=A.transpose()
permettent d'étiqueter B
la transposée de A
).Exemple. Pour le premier calcul, on exécutera :
A = np.matrix([[1,2],[3,4]])
print(A**2)
B = np.matrix([[-1,0,1],[2,2,3]])
print(A*B)
print(B*A)
D = np.matrix([[1],[-1],[3]])
Dt = D.transpose()
print(Dt*D)
Exercice 1.2 : Pour tout $n$ entier naturel, on définit la matrice carrée d'ordre $n$ : $$ A_n=\begin{pmatrix} 1&1&\cdots&1\\0&1&\cdots&1\\ \vdots&\ddots&\ddots&\vdots\\0&\dots&0&1\end{pmatrix} $$
Exercice 1.3 : Trouver toutes les matrices $X$ qui commutent avec la matrice $$A=\begin{pmatrix} 1&2\\0&3\end{pmatrix}$$ Vérifier votre résultat en complétant le code de la fenêtre ci-dessous à l'aide de cas particuliers de matrices $X$ :
A = np.matrix([[1,2],[0,3]])
x,y = 1,20
X = np.matrix([[x,y],[0,x+y]]) # à compléter.
print(A*X == X*A)
[[ True True] [ True True]]
Exercice 1.4 : On considère la matrice $$A=\begin{pmatrix}-1&0&-2\\1&1&1\\1&0&2\end{pmatrix}$$
On cherche à résoudre l'équation d'inconnues réelles : $(xI_3+yA)^2=I_3$
où $I_3$ désigne la matrice identité de $\mathcal{M}_3(\mathbb{R})$.
A = np.matrix([[-1,0,-2],[1,1,1],[1,0,2]]) # à compléter
I3 = np.identity(3)
# question 1. A^2 = ?
x,y = -1,2 # à compléter
B = x*I3+y*A
print(B**2==I3)
[[ True True True] [ True True True] [ True True True]]
Exercice 1.5 : On considère la matrice $$N=\begin{pmatrix}2&-2&1\\2&-3&2\\-1&2&0\end{pmatrix}$$
Préciser notamment la relation qui existe entre $(u_n,v_n)$ et $(u_{n+1},v_{n+1})$.
3. Vérifier que $u_{n+1}+v_{n+1}=u_n+v_n$.
4. En déduire que $u_{n+1}=-3u_n+1$.
5. Exprimer $(u_n,v_n)$ en fonction de $n$ et tester votre réponse en complétant les lignes de code ci-dessous :
N = np.matrix([...]) # à compléter
I3 = np.identity(3)
n = ...
u,v = ... , ...
print(N**(n+1) == u*N+v*I3)
Indication : Vous avez dû obtenir u,v = 1/4+3/4*(-3)**n,3/4-3/4*(-3)**n
Exercice 1.6 : Avec Python il est facile de faire des sommes et des produits de matrices lorsqu'ils sont possibles. Il peut néanmoins arriver (c'est au programme) qu'on demande d'écrire des fonctions Python permettant de faire ces opérations.
Je rappelle ci-dessous deux façons possible d'écrire une fonction permettant d'additionner, quand leur taille le permet, deux matrices $A$ et $B$.
def sommeM1(A,B):
l1,c1 = A.shape
l2,c2 = B.shape
if l1!=l2 or c1!=c2:
print('somme impossible')
return []
else:
return np.matrix([[A[i,j]+B[i,j] for j in range(c1)]for i in range(l1)])
def sommeM2(A,B):
l1,c1 = A.shape
l2,c2 = B.shape
if l1!=l2 or c1!=c2:
print('somme impossible')
return []
else:
C = np.zeros((l1,c1)) # initialisation de la liste qui formera la somme A+B
for i in range(l1):
for j in range(c1):
C[i,j]=A[i,j]+B[i,j]
return C
On vérifie ces deux fonctions à l'aide des matrices de l'exercice 1 :
A = np.matrix([[1,2],[3,4]])
B = np.matrix([[-1,0,1],[2,2,3]])
E = np.matrix([[1,-1],[-1,2]])
C1 = sommeM1(A,B)
C2 = sommeM2(A,E)
print(C2)
print(C2 == A+E)
somme impossible [[ 2. 1.] [ 2. 6.]] [[ True True] [ True True]]
Compléter maintenant les deux fonctions ci-dessous permettant respectivement de faire le produit de deux matrices $A$ et $B$ quand c'est possible et de calculer $A^n$ si $A$ est une matrice carrée et $n\in\mathbb{N}^*$.
On rappelle que si $A=(a_{i,j})\in\mathcal{M}_{n,p}(\mathbb{R})$ et $B=(b_{i,j})\in\mathcal{M}_{p,m}(\mathbb{R})$, alors les coefficients de $C=A*B$ sont :
$$c_{i,j}=\displaystyle\sum_{k=1}^p a_{i,k}\cdot b_{k,j}$$
def produitM(A,B):
la,ca = A.shape
lb,cb = B.shape
if ca != lb : # à compléter
print('produit impossible')
C = []
else:
C = np.zeros((la,cb))
for i in range(la):
for j in range(cb):
for k in range(lb):
C[i,j] = C[i,j]+ A[i,k]*B[k,j]
return C
C = produitM(A,B)
print(C)
print(A*B)
D = produitM(B,E)
def puissanceM(A,n):
# on s'obligera à utiliser la fonction 'produitM()'
la,ca = A.shape
if la != ca:
print('Puissance de A impossible')
return []
else:
pA = A # initialisation de la puissance de A à A (n=1)
for k in ...:
pA = ...
return pA
Application : Utiliser la fonction précédente pour évaluer les puissances successives ($n=1,2,3,\cdots$) pour les matrices suivantes :
$$M=\begin{pmatrix}0&1&0\\0&0&1\\0&0&0\end{pmatrix}\qquad N=\begin{pmatrix}1&1&1\\1&1&1\\1&1&1\end{pmatrix}\qquad J=\begin{pmatrix}0&1&0\\0&0&1\\1&0&0\end{pmatrix}$$Fournir ensuite une preuve algébrique de vos résultats.
M=np.matrix([[0,1,0],[0,0,1],[0,0,0]])
pM = puissanceM(M,1)
print(pM)
Soit $A\in\mathcal{M}_{n,p}(\mathbb{R})$. On appelle rang de la matrice $A$ et on note $\textrm{rg}(A)$ le nombre de pivots (non nuls) à l'issue de l'application de l'algorithme dit "du pivot de Gauss".
A titre d'exemples :
Exemple 1 :$$\textrm{rg}\begin{pmatrix}1&1\\-1&1\end{pmatrix}=\textrm{rg}\begin{pmatrix}1&1\\0&2\end{pmatrix}=2\qquad L_2\leftarrow L_1+L_2$$
Remarque : Dans le cas des matrices d'ordre $2$, il est en fait inutile de passer par un pivot de Gauss car $$\textrm{rg}(A)=2\Leftrightarrow \det(A)\neq 0$$
Ici, $\det\begin{pmatrix}1&1\\-1&1\end{pmatrix}=1+1=2\neq 0$ donc on a directement $\textrm{rg}(A)=2$.
Exemple 2: $$\textrm{rg}\begin{pmatrix}1&1&1\\1&1&1\\1&1&1\end{pmatrix}=\textrm{rg}\begin{pmatrix}1&1&1\\0&0&0\\0&0&0\end{pmatrix}=1\qquad L_2\leftarrow L_2-L_1,\quad L_3\leftarrow L_3-L_1$$
Remarque : Quand c'est possible, on peut utiliser une fonction Python pour déterminer le rang d'une matrice : np.linalg.matrix_rank(A)
A = np.matrix([[1,1],[-1,1]])
B = np.matrix([[1,1,1],[1,1,1],[1,1,1]])
rA = np.linalg.matrix_rank(A)
rB = np.linalg.matrix_rank(B)
print('rang de A : ',rA,'et rang de B : ',rB)
Exercice 7 : Déterminer en utilisant la méthode du pivot de Gauss, le rang des matrices ci-dessous.
Vérifier vos résultats grâce à la fonction np.linalg.matrix_rank()
.
$$A_1=\begin{pmatrix}2&-3&1\\4&-5&2\\5&-7&3\end{pmatrix} ; \quad A_2=\begin{pmatrix}-1&0&-2\\1&1&1\\1&0&2\end{pmatrix} ; \quad A_3=\begin{pmatrix}-1&2&3\\6&-1&3\\0&-1&2\end{pmatrix} ; \quad A_4=\begin{pmatrix}1&0&1\\1&1&0\\2&1&1\end{pmatrix} ;$$
Lien avec l'inversibilité : Une matrice $A$ d'ordre $n$ (c'est-à-dire matrice carrée de $n$ lignes et $n$ colonnes) est inversible si et seulement si $\textrm{rg}(A)=n$.
Dans l'exercice précédent, il est donc immédiat que $A_1$ et $A_3$ sont inversibles tandis que $A_2$ et $A_3$ ne le sont pas.
Exercice 1.8 : Déterminer les valeurs de $a$ réelles pour lesquelles les matrices ci-dessous ne sont pas inversibles.
Attention : Dans la méthode du pivot de Gauss, les pivots ne peuvent pas être nuls.
Vérifier vos réponses en vous assurant que pour ces valeurs de $a$, on obtient $\textrm{rg}(A)<3$.
a = -3# à compléter
A = np.matrix([[-2-a,-2,1],[-2,1-a,-2],[1,-2,-2-a]]) # exemple pour la matrice $A$ ci-dessus.
print(np.linalg.matrix_rank(A))
1
Indication : $rg(B)=2$ si $a=-4$ et pour $C$ on obtient un rang égale à 2 ou 1 selon que a = 3 ou -3
On rappelle que tout système à $p$ inconnues $x_1,\cdots,x_p$ et $n$ équations de la forme :
$$(S)\quad\begin{cases}a_{1,1}x_1+\cdots+a_{1,p}x_p&=b_1\\a_{2,1}x_1+\cdots+a_{2,p}x_p&=b_2\\ \vdots&\vdots\\ a_{n,1}x_1+\cdots+a_{n,p}x_p&=b_n\end{cases}$$peut s'exprimer matriciellement sous la forme : $A\cdot X = B$ où :
$$A=\begin{pmatrix}a_{1,1}&\cdots& a_{1,p}\\ \vdots&&\vdots \\ a_{n,1}&\cdots&a_{n,p}\end{pmatrix} ; \quad X=\begin{pmatrix}x_1\\ \vdots \\x_p\end{pmatrix} ; \quad B=\begin{pmatrix}b_1\\ \vdots \\ b_n\end{pmatrix}$$Sous cette condition, $A$ est dite matrice associée au système et on a : $\textrm{rg}(S)=\textrm{rg}(A)$.
Important ! Cette dernière remarque est essentielle car, si $\textrm{rg}(S)=n=p$ alors $A$ est une matrice inversible et le système $(S)$ admet une unique solution, à savoir : $$X = A^{-1}\cdot B$$
Exercice 1.9: Prendre une à une les matrices de l'exercice 7 et résoudre dans chaque cas le système : $$A\cdot X = B\qquad \text{où }X=\begin{pmatrix}x\\y\\z\end{pmatrix}\text{ et }B=\begin{pmatrix}1\\1\\1\end{pmatrix}$$
Remarque : Dans le cas des matrices $A_1$ et $A_3$ qui sont inversible, vérifier vos réponses grâce à Python en exécutant au choix :
X=np.linalg.solve(A,B)
ou
X=np.linalg.inv(A)*B
# exemple pour la résolution du système : x+y=1 et x-y=0 associé à une matrice A de rang 2 donc inversible :
A = np.matrix([[1,1],[1,-1]])
B = np.matrix([[1],[0]])
X1 = np.linalg.solve(A,B) # première méthode
X2 = np.linalg.inv(A)*B # deuxième méthode
print(X1)
print(X2)
Dans le cas d'une matrice carrée d'ordre $n$ inversible (c'est-à-dire de rang $n$), il est indispensable de savoir produire son inverse rapidement.
Dans le cas d'une matrice d'ordre $2$, c'est immédiat car nous disposons d'une formule de cours à connaître par coeur.
Si $n>2$, nous passerons par la résolution d'un système.
$$A^{-1}=\cfrac{1}{ad-bc}\begin{pmatrix}d&-b\\-c&a\end{pmatrix}$$$A=\begin{pmatrix}a&b\\c&d\end{pmatrix}$ inversible $\Leftrightarrow \det(A)=ad-bc\neq 0$ et dans ce cas :
Exercice 1.10 : Dire parmi ces matrices celles qui sont inversibles et donner leur inverse.
$$A_1=\begin{pmatrix}1&2\\-2&-1\end{pmatrix} ; \quad A_2=\begin{pmatrix}1&2\\-1&-2\end{pmatrix} ; \quad A_1=\begin{pmatrix}1&2\\3&4\end{pmatrix} ; $$Note : Dans le cas des matrices inversibles, vous validerez votre réponse en exécutant : np.linalg.inv(A)
$$A=\begin{pmatrix}a_1&0&\cdots&0\\0&a_2&\cdots&0\\ \vdots&\ddots&\ddots&\vdots\\ 0&\cdots&0&a_n\end{pmatrix},a_i\in\mathbb{R}^*\Rightarrow A^{-1}=\begin{pmatrix}1/a_1&0&\cdots&0\\0&1/a_2&\cdots&0\\ \vdots&\ddots&\ddots&\vdots\\ 0&\cdots&0&1/a_n\end{pmatrix}$$Si $A$ est une matrice diagonale dont les coefficients sont tous non nuls, alors elle est inversible et son inverse s'obtient en inversant chacun des termes de la diagonale. Soit :
Sinon, comme tout système de la forme $A\cdot X = B$ admet une unique solution égale à $X=A^{-1}B$, l'idée est d'exprimer l'unique solution $X$ de ce système en fonction des coordonnées de $B$.
Exemple : Déterminons l'inverse de $A_1=\begin{pmatrix}2&-3&1\\4&-5&2\\5&-7&3\end{pmatrix}$ définie dans l'exercice 7 en considèrant le système $(S)\quad A_1X=B$ avec $X=\begin{pmatrix}x\\y\\z\end{pmatrix}$ et $B=\begin{pmatrix}a\\b\\c\end{pmatrix}$.
On applique la méthode du pivot de Gauss en commençant par $L_2\leftarrow L_2-2L_1$ et $L_3\leftarrow 2L_3-5L_1$ :
Conclusion : $A^{-1}=\begin{pmatrix}-1&2&-1\\-2&1&0\\-3&-1&2\end{pmatrix}$.
# Vérifions avec Python :
A = np.matrix([[2,-3,1],[4,-5,2],[5,-7,3]])
np.linalg.inv(A)
matrix([[ -1.00000000e+00, 2.00000000e+00, -1.00000000e+00], [ -2.00000000e+00, 1.00000000e+00, 5.18104078e-16], [ -3.00000000e+00, -1.00000000e+00, 2.00000000e+00]])
Exercice 1.11 : Justifier l'inversibilité des matrices suivantes et les inverser :
$$B_1=\begin{pmatrix}1&2&-1\\2&4&-1\\-2&-5&3\end{pmatrix} ; \quad B_2=\begin{pmatrix}0&-1&0\\1&1&-1\\-2&-1&0\end{pmatrix} ; \quad B_3=\begin{pmatrix}2&1&3\\1&3&-1\\3&-1&7\end{pmatrix} ; \quad B_4=\begin{pmatrix}-2&-2&2\\-1&2&2\\1&-2&3\end{pmatrix}$$
Remarque importante : Il est une autre méthode pour inverser une matrice $A$.
Elle repose sur la mise en évidence d'un polynôme $P=a_0+a_1X+\cdots+a_nX^n\in\mathbb{R}[X]$ avec $a_0\neq 0$ tel que $P(A)=0$.
En effet :
$$P(A)=a_0I+a_1A+\cdots+a_nA^n=0\Leftrightarrow A\cdot\left(-\cfrac{1}{a_0}(a_1I+\cdots+a_nA^{n-1})\right)=I= \left(-\cfrac{1}{a_0}(a_1I+\cdots+a_nA^{n-1})\right)\cdot A$$
$$\Leftrightarrow A\text{ inversible et }A^{-1}=-\cfrac{1}{a_0}\left(a_1I+\cdots+a_nA^{n-1}\right)$$
Exercice 1.12 : Soit la matrice $A=\begin{pmatrix}0&1&1\\1&0&1\\1&1&0\end{pmatrix}$.
# Vérification grâce à la fonction `np.linalg.inv(A)` de Python :
On supposera dans ce qui suit que $F$ est un espace vectoriel, sous-espace vectoriel de $E=\mathbb{R}^n$ avec $n\in\mathbb{N}^*$.
On rappelle à cette occasion la caractérisation des sous-espaces vectoriels :
Proposition 1 : $F$ est un SEV de $E$ si :
Note : C'est le moment de reprendre votre T.D. sur ce chapitre et de retravailler de telles démonstrations !
Définition 1 : Soit $p\in\mathbb{N}^*$. Une famille $\mathcal{F}=\{u_1,\cdots,u_p\}$ de vecteurs de $F$ est libre (non liée) si :
$$\forall \lambda_1,\cdots,\lambda_p\in\mathbb{R},\quad \lambda_1u_1+\cdots+\lambda_pu_p=0_E\Rightarrow \lambda_1=0=\cdots=\lambda_p$$
Définition 2 : Soit $p\in\mathbb{N}^*$. Une famille $\mathcal{F}=\{u_1,\cdots,u_p\}$ de vecteurs de $F$ est génératrice si :
$$\forall u\in F,\quad \exists \lambda_1,\cdots,\lambda_p\in\mathbb{R}/u=\lambda_1u_1+\cdots+\lambda_pu_p$$
Définition 3 : Toute famille $\mathcal{F}$ à la fois libre et génératrice de $F$ est une base de $F$. On note alors : $$\dim(F)=\textrm{card}(\mathcal{F})$$
Définition 4 : Soit $\mathcal{F}=\{u_1,\cdots,u_p\}$ une famille de vecteurs de $F$. On appelle rang de $\mathcal{F}$ et on note $\textrm{rg}(\mathcal{F})$ la dimension de l'espace vectoriel engendré par les vecteurs de $\mathcal{F}$. Soit : $$\textrm{rg}(\mathcal{F})=\dim\left(\textrm{Vect}\{u_1,\cdots,u_p\}\right)$$
Cette dernière définition est fondamentale car on obtient grâce à elle qu'une famille de $p$ vecteurs est libre si, et seulement si, son rang vaut $p$.
Remarque : Soit $M$ la matrice des coordonnées de la famille $\mathcal{F}=\{u_1,\cdots,u_p\}$ dans une base de $E$. Alors : $$\textrm{rg}(M)=\textrm{rg}(\mathcal{F})$$
Application 1 :
Soit $\mathcal{F}_1=\{u_1,u_2\}$ où $u_1=(1,-1)$ et $u_2=(1,1)$ dans la base canonique de $E=\mathbb{R}^2$.
Cette famille est libre car la matrice des coordonnées de cette famille dans la base canonique de $\mathbb{R}^2$ est : $$M=\begin{pmatrix}1&1\\-1&1\end{pmatrix}\quad\text{! les coordonnées des vecteurs sont en colonne.}$$ et $$\textrm{rg}(M)=\textrm{rg}(\mathcal{F}_1)=2=\textrm{card}(\mathcal{F}_1)$$
Application 2 :
Soit $\mathcal{F}_2=\{u_1,u_2,u_3\}$ où $u_1=(1,1,1)$, $u_2=(1,1,1)$ et $u_3=(1,1,1)$ dans la base canonique de $E=\mathbb{R}^3$.
Cette famille est liée car la matrice des coordonnées de cette famille dans la base canonique de $\mathbb{R}^3$ est : $$M=\begin{pmatrix}1&1&1\\1&1&1\\1&1&1\end{pmatrix}$$ et $$\textrm{rg}(M)=\textrm{rg}(\mathcal{F}_2)=1<\textrm{card}(\mathcal{F}_2)$$
Exercice 2.1 : Pour chacune des familles suivantes de vecteurs de $\mathbb{R}^n$ ($n=2,3$ ou $4$) déterminer si elle est libre ou liée.
Par ailleurs, préciser dans chaque cas la dimension de $F=\textrm{Vect}(\mathcal{F})$.
# Avec l'aide de Python... par exemple pour la quatrième famille de vecteurs (attention à l'écriture en colonnes !)
M = np.matrix([[1,1,-1],[2,0,-2],[-1,1,-3]])
print(np.linalg.matrix_rank(M))
Exercice 2.2: Déterminer si le vecteur $v$ appartient au sous-espace vectoriel $F$ dans chacun des cas suivants et, si c'est le cas, exprimer $v$ dans une base de $F$ :
Indication : On rappelle que si $\{u_1,\cdots,u_p\}$ est une famille libre, $$v\in\mathrm{Vect}\{u_1,\cdots,u_p\}\Leftrightarrow \{v,u_1,\cdots,u_p\}\text{ est une famille liée}$$
# exemple pour répondre avec Python à la question 1. de l'exercice 2 :
A = np.matrix([[0,2,1],[-1,1,1],[3,3,0]])
np.linalg.matrix_rank(A)
Définition 1 : Soit $f$ une application de $E_p=\mathbb{R}^p$ dans $F_n=\mathbb{R}^n$.
$f$ est une application linéaire si, et seulement si :
$$\forall u,v\in E_p,\quad \forall \alpha,\beta\in\mathbb{R},\quad f(\alpha u+\beta v)=\alpha f(u)+\beta f(v)$$
Note : Si $E_p=F_n$ alors on dira que $f$ est un endomorphisme de $E_p$ et on notera $f\in\mathcal{L}(E_p)$
Définition 2 - Proposition 1 : Soit $f\in\mathcal{L}(E_p,F_n)$.
Proposition 2 : $f\in\mathcal{L}(E_p,F_n)$ est injective ssi $\mathrm{Ker}(f)=\{0_E\}$.
Proposition 3 : $f\in\mathcal{L}(E_p,F_n)$ est surjective ssi $\mathrm{Im}(f)=F_n$
Proposition 4 : Formule du rang $$\dim(E_p)=\dim(\mathrm{Ker}(f))+\dim(\mathrm{Im}(f))=\dim(\mathrm{Ker}(f))+\textrm{rg}(f)$$
Proposition 5 : Si $f\in\mathcal{L}(E_p)$ avec $\dim(E_p)$ finie, alors : $$f\text{ injective }\Leftrightarrow f\text{ surjective }\Leftrightarrow f\text{ bijective}$$
Toutes ces propositions sont à connaître. Dans la pratique, on demandera le plus souvent de montrer qu'une application est linéaire (il faut vous entraîner sur la base de votre TD de BCPST 1ère année) puis on demandera son noyau, son image et de préciser si elle est bijective ou non.
Soit $f\in\mathcal{L}(E_p,F_n)$ où $E_p=\mathbb{R}^p$ et $F_n=\mathbb{R}^n$.
On écrira : $u=(x_1,x_2,\cdots,x_p)\in \mathrm{Ker}(f)\Leftrightarrow f(x_1,x_2,\cdots,x_p)=0_F$ et on s'attachera à résoudre le système.
Exemple 1 : $f_1$ définie sur $\mathbb{R}^2$ par $f_1(x,y)=(-x-y,2x+2y)$, $\forall (x,y)\in\mathbb{R}^2$.
$$u=(x,y)\in\mathrm{Ker}(f)\Leftrightarrow f_1(x,y)=(0,0)\Leftrightarrow \begin{cases}-x-y&=0\\2x+2y&=0\end{cases}\Leftrightarrow x+y=0$$
Conclusion 1 : $\mathrm{Ker}(f)=\{(x,y)\in E_2/x+y=0\}=\{(x,-x),x\in\mathbb{R}\}=\mathrm{Vect}\{u_1\}$ où $u_1=(1,-1)$
Exemple 2 : $f_2$ définie sur $\mathbb{R}^2$ par $f_2(x,y)=(4x+y,x-y,2x+3y)$, $\forall (x,y)\in\mathbb{R}^2$.
$$u=(x,y)\in\mathrm{Ker}(f_2)\Leftrightarrow f_2(x,y)=(0,0)\Leftrightarrow \begin{cases}4x+y&=0\\x-y&=0\\2x+3y&=0\end{cases}\Leftrightarrow x=0=y$$
Conclusion 2 : $\mathrm{Ker}(f_2)=\{(x,y)\in E_2/x=0=y\}=\{(0,0)\}$ ; $f_2$ est injective
Exercice 3.1 : Déterminer une base et la dimension du noyau des applications linéaires suivantes. Préciser si ces applications sont injectives, surjectives, bijectives.
Soit $f\in\mathcal{L}(E_p,F_n)$ où $E_p=\mathbb{R}^p$ et $F_n=\mathbb{R}^n$.
On écrira : $v=(a_1,a_2,\cdots,a_n)\in \mathrm{Im}(f)\Leftrightarrow \exists (x_1,x_2,\cdots,x_p)\in E_p\quad/\quad f(x_1,x_2,\cdots,x_p)=v\quad (S)$.
Autrement dit, $v\in\mathrm{Im}(f)$ si le système $(S)$ admet au moins une solution ou encore n'est pas "incompatible".
Reprenons les exemples 1 et 2 du paragraphe précédent consacré à la recherche de noyaux :
Exemple 1 : $f_1$ définie sur $\mathbb{R}^2$ par $f_1(x,y)=(-x-y,2x+2y)$, $\forall (x,y)\in\mathbb{R}^2$.
On a déjà vu que $\dim(\textrm{Ker}(f_1))=1$ donc d'après la formule du rang : $\dim(\textrm{Im}(f_1))=2-1=1$.
$f_1$ n'est pas surjective et ça vaut la peine de chercher à déterminer $\textrm{Im}(f_1)$ !
$$v=(a,b)\in\mathrm{Im}(f_1)\Leftrightarrow \exists (x,y)\in\mathbb{R}^2/f_1(x,y)=(a,b)\Leftrightarrow\exists(x,y)\in\mathbb{R}^2/ \begin{cases}-x-y&=a\\2x+2y&=b\end{cases}\quad (S)$$$$(S)\Leftrightarrow \begin{cases}-x-y&=a\\0&=b+2a\end{cases}\quad L_2\leftarrow L_2+2L_1$$C'est fini ! Les vecteurs $(a,b)$ de $\textrm{Im}(f_1)$ sont les vecteurs de $\mathbb{R}^2$ pour lesquels $2a+b=0$.
Conclusion : $\textrm{Im}(f_1)=\{(a,b)\in\mathbb{R}^2/2a+b=0\}=\{(a,-2a),a\in\mathbb{R}\}=\{a(1,2),a\in\mathbb{R}\}$ ; On retrouve $\dim(\textrm{Im}(f_1))=1$
Exemple 2 : $f_2$ définie sur $\mathbb{R}^2$ par $f_2(x,y)=(4x+y,x-y,2x+3y)$, $\forall (x,y)\in\mathbb{R}^2$.
On a vu que $\dim(\textrm{Ker}(f_2))=0$ donc, d'après la formule du rang : $\dim(\textrm{Im}(f_2))=2-0=2\neq \dim(\mathbb{R}^3)$.
$f_2$ n'est donc pas surjective ; ça vaut encore la peine de chercher à déterminer $\textrm{Im}(f_2)$ !
$$v=(a,b,c)\in\mathrm{Im}(f_2)\Leftrightarrow \exists (x,y)\in\mathbb{R}^2/f_1(x,y)=(a,b,c)\Leftrightarrow\exists(x,y)\in\mathbb{R}^2/ \begin{cases}4x+y&=a\\x-y&=b\\2x+3y&=c\end{cases}\quad (S)$$Résolvons ce système par la méthode du pivot de gauss en échangeant $L_1$ et $L_2$ puis en faisant : $L_2\leftarrow L_2-4L_1$, $L_3\leftarrow L_3-2L_1$ :
$$(S)\Leftrightarrow \begin{cases}x-y&=b\\4x+y&=a\\2x+3y&=c\end{cases}\Leftrightarrow \begin{cases}x-y&=b\\5y&=a-4b\\5y&=c-2b\end{cases}\Leftrightarrow \begin{cases}x-y&=b\\5y&=a-4b\\0&=a-2b-c\end{cases}\quad L_3\leftarrow L_3-L_2$$C'est fini ! Les vecteurs $(a,b,c)$ de $\textrm{Im}(f_2)$ sont les vecteurs de $\mathbb{R}^3$ pour lesquels $a-2b-c=0$.
Conclusion : $\textrm{Im}(f_2)=\{(a,b,c)\in\mathbb{R}^3/a=2b+c\}=\{(2b+c,b,c),b,c\in\mathbb{R}\}=\{b(2,1,0)+c(1,0,1),b,c\in\mathbb{R}\}$
Les vecteurs $(2,1,0)$ et $(1,0,1)$ étant linéairement indépendants, une base de $\textrm{Im}(f_2)$ est $\mathcal{B}_2=(u_1,u_2)$ où $u_1=(2,1,0)$ et $u_2=(1,0,1)$.
Exercice 3.2 : Reprenons les applications linéaires de l'exercice 3.1 mais déterminer cette fois une base et la dimension de l'image de ces applications.
Soit $f\in\mathcal{L}(E_p,F_n)$ où $E_p=\mathbb{R}^p$ et $F_n=\mathbb{R}^n$.
On note $\mathcal{B}_E=(e_1,\cdots,e_p)$ et $\mathcal{B}_F=(c_1,\cdots,c_n)$ deux bases respectivement de $E_p$ et de $F_n$.
On explicite la famille de vecteurs $$f(\mathcal{B}_E)=\{f(e_1),f(e_2),\cdots,f(e_p)\}$$ dans la base $\mathcal{B}_F$.
Pour obtenir la matrice de $f$ dans les bases $\mathcal{B}_E$ et $\mathcal{B}_F$, notée $\mathcal{M}_{\mathcal{B}_E,\mathcal{B}_F}(f)$ il suffit alors d'écrire la matrice des coordonnées de cette famille dans la base $\mathcal{B}_F$.
Exemple 1 : $f_1$ définie sur $\mathbb{R}^2$ par $f_1(x,y)=(-x-y,2x+2y)$, $\forall (x,y)\in\mathbb{R}^2$.
$f_1$ est un endomorphisme de $\mathbb{R}^2$ et on peut prendre $\mathcal{B}_E=\mathcal{B}=\mathcal{B_F}$, base canonique de $\mathbb{R}^2$.
Alors : $$A_1=\mathcal{M}_{\mathcal{B}}(f_1)=\begin{pmatrix}-1&-1\\2&2\end{pmatrix}$$
Exemple 2 : $f_2$ définie sur $\mathbb{R}^2$ par $f_2(x,y)=(4x+y,x-y,2x+3y)$, $\forall (x,y)\in\mathbb{R}^2$.
Notons $\mathcal{B_2}$ la base canonique de $\mathbb{R}^2$ et $\mathcal{B}_3$ la base canonique de $\mathbb{R}^3$ et déterminons $A_2=\mathcal{M}_{\mathcal{B}_2,\mathcal{B}_3}(f_2)$ :
Alors : $$A_2=\mathcal{M}_{\mathcal{B}_2,\mathcal{B}_3}(f_2)=\begin{pmatrix}4&1\\1&-1\\2&3\end{pmatrix}$$
Proposition 1 : Si $A=\mathcal{M}_{\mathcal{B}_E,\mathcal{B}_F}(f)$, alors $$\textrm{rg}(A)=\textrm{rg}(f)=\dim(\textrm{Im}(f))$$
Corollaire : Si $f$ est un endomorphisme de $E_n$ et $A=\mathcal{M}_{\mathcal{B}_E}(f)$, alors :
$$f\text{ bijective}\Leftrightarrow \textrm{rg}(f)=\textrm{rg}(A) = n \Leftrightarrow A\text{ inversible }$$Exercice 3.3 : Ecrire, dans les bases canoniques des espaces vectoriels de départ et d'arrivée, la matrices des applications linéaires définies dans l'exercice 3.1 et que nous rappelons ci-dessous :
Remarque : On précisera dans chaque cas le rang des matrices, les éventuelles matrices inversibles, et on fera le lien avec la dimension du noyau et de l'image.
Réponses pour l'exercice 2.1 :
- Liée car $u_1$ et $u_2$ sont colinéaires. $\dim(F)=2$ et même $F=\mathbb{R}^2$.
- Libre car un seul vecteur non nul. $\dim(F)=1$ (droite vectorielle de $\mathbb{R}^2$.
- Liée car trois vecteurs dans un espace vectoriel de dimension $2$. $\dim(F)=2$ (plan vectoriel de $\mathbb{R}^3$).
- Libre car le rang de la famille vaut $3=\dim(F)$. $F=\mathbb{R}^3$.
- Liée car le rang de la famille vaut $2$. Vérifier que $u_2$ s'exprime facilement en fonction de $u_1$ et $u_3$. $\dim(F)=2$ (plan vectoriel de $\mathbb{R}^3$).
- Liée car quatre vecteurs dans un espace vectoriel de dimension $3$. $\dim(F)=3$. On a donc $F=\mathbb{R}^3$.
- Liée car le rang de la famille vaut $2$ (on note que $u_1+u_3=3u_2$). $\dim(F)=2$ (plan vectoriel).
- Libre car le rang de cette famille vaut $3$. $\dim(F)=3)$ (on parle d'hyperplan de $\mathbb{R}^3$).
- Liée car le rang de cette famille vaut $1$. $\dim(F)=1$ (droite vectorielle de $\mathbb{R}^4$).
Réponse Exercice 2.2 :
- $\{u_1,u_2\}$ est une famille libre de $F$ (vecteurs non colinéaires) donc $\mathcal{B}=(u_1,u_2)$ est une base de $F$.
Par ailleurs $\textrm{rg}\{v,u_1,u_2\}=2$ donc $v\in F$.
On en déduit que : $\exists !(a,b)\in\mathbb{R}^2/v=au_1+bu_2$.
Par identification qui donne lieu à un système de 2 équations à 2 inconnues, on obtient immédiatement $a=1,b=-2$ soit : $$v=u_1-2u_2=(1,-2)_{\mathcal{B}}$$
- $\{u_1,u_2\}$ est une famille libre de $F$ (vecteurs non colinéaires) donc $\mathcal{B}=(u_1,u_2)$ est une base de $F$.
Par ailleurs $\textrm{rg}\{v,u_1,u_2\}=3$ donc $v\notin F$.
- $\{u_1,u_2,u_3\}$ est une famille libre de $F$ (car $\textrm{rg}\{u_1,u_2,u_3\}=3$) donc $\mathcal{B}=(u_1,u_2,u_3)$ est une base de $F$.
Par ailleurs $\textrm{rg}\{v,u_1,u_2,u_3\}=3$ donc $v\in F$.
On en déduit que : $\exists !(a,b,c)\in\mathbb{R}^3/v=au_1+bu_2+cu_3$.
Par identification qui donne lieu à un système de 3 équations à 3 inconnues, on obtient immédiatement $a=1,c=3,c=-2$ soit : $$v=u_1+3u_2-2u_3=(1,3,-2)_{\mathcal{B}}$$
Réponse Exercice 3.1 :
- $\textrm{Ker}(f_1)=\{(x,y)\in\mathbb{R}^2/x=0=y\}=\{0_{\mathbb{R}^2}\}$.
$f_1$ injective et non surjective car d'après la formule du rang : $\dim(\textrm{Im}(f_1))=2-0=2\neq 3$ donc $\textrm{Im}(f_1)\neq \mathbb{R}^3$
- $\textrm{Ker}(f_2)=\{(x,y,z)\in\mathbb{R}^3/y=-2x\text{ et }z=x\}=\{(x,-2x,x),x\in\mathbb{R}\}=\{x(1,-2,1),x\in\mathbb{R}\}$.
$\textrm{Ker}(f_2)$ a pour base $\mathcal{B}_2=(u_1)$ où $u_1=(1,-2,1)$ ; $\dim(\textrm{Ker}(f_2))=1$.
$f_1$ est non injective et surjective car d'après la formule du rang : $$\dim(\textrm{Im}(f_2))=3-1=2=\dim(\mathbb{R}^2)$$ et comme $\textrm{Im}(f_1)\subset \mathbb{R}^2$, on a : $\textrm{Im}(f_1)= \mathbb{R}^2$.
- $\textrm{Ker}(f_3)=\{(x,y,z)\in\mathbb{R}^3/x=y=z\}=\{(x,x,x),x\in\mathbb{R}\}=\{x(1,1,1),x\in\mathbb{R}\}$.
$\textrm{Ker}(f_3)$ a pour base $\mathcal{B}_3=(u_1)$ où $u_1=(1,1,1)$ ; $\dim(\textrm{Ker}(f_3))=1$.
$f_3$ est non injective et non surjective car $f_3$ est un endomorphisme de $\mathbb{R^3}$ et non injective équivaut à non surjective.
- $\textrm{Ker}(f_4)=\{(x,y,z)\in\mathbb{R}^3/x=y=z=0\}=\{0_{\mathbb{R}^3}\}$.
$f_4$ est injective et surjective car $f_4$ est un endomorphisme de $\mathbb{R}^3$.
- $\textrm{Ker}(f_5)=\{(x,y,z)\in\mathbb{R}^3/y=z\}=\{(x,y,y),x,y\in\mathbb{R}\}=\{x(1,0,0)+y(0,1,1),x,y\in\mathbb{R}\}$.
$\textrm{Ker}(f_5)$ a pour base $\mathcal{B}_5=(u_1,u_2)$ où $u_1=(1,0,0)$, $u_2=(0,1,1)$ ; $\dim(\textrm{Ker}(f_5))=2$.
$f_5$ est non injective et non surjective car $f_5$ est un endomorphisme de $\mathbb{R^3}$.
Réponse exercice 3.2 :
$$(a,b,c)\in\textrm{Im}(f_1)\Leftrightarrow \exists (x,y)\in\mathbb{R}^2/\begin{cases}x-y&=a\\x&=b\\x+2y&=c\end{cases}\Leftrightarrow \exists (x,y)\in\mathbb{R}^2/\begin{cases}y&=b-a\\x&=b\\x+2y&=3b-2a=c\end{cases}$$
- $f_1$ est non surjective et $\dim(\textrm{Im}(f_1))=2$ d'après l'ex. 3.1.
Conlusion 1 : $\textrm{Im}(f_1)=\{(a,b,c)\in\mathbb{R}^3/2a-3b+c=0\}=\{(a,b,-2a+3b),a,b\in\mathbb{R}\}=\textrm{Vect}\{(1,0,-2),(0,1,3)\}$.
$\textrm{Im}(f_1)$ a pour base $\mathcal{B'}_1=(v_1,v_2)$ où $v_1=(1,0,-2)$ et $v_2=(0,1,3)$ car ces vecteurs sont non colinéaires.
- $f_2$ est surjective d'après l'exercice 3.1 donc $\textrm{Im}(f_2)=\mathbb{R^2}$.
On peut choisir la base canonique de $\mathbb{R}^2$ comme base de $\textrm{Im}(f_2)$.
$$(a,b,c)\in\textrm{Im}(f_3)\Leftrightarrow \exists (x,y,z)\in\mathbb{R}^3/\begin{cases}x-y&=a\\y-z&=b\\z-x&=c\end{cases}\Leftrightarrow \exists (x,y,z)\in\mathbb{R}^3/\begin{cases}x-y&=a\\x-z&=a+b\\z-x&=c\end{cases}\Leftrightarrow\\ \Leftrightarrow\exists (x,y,z)\in\mathbb{R}^3/\begin{cases}x-y&=a\\x-z&=a+b\\0&=a+b+c\end{cases}$$
- $f_3$ est non surjective et $\dim(\textrm{Im}(f_3))=2$ d'après l'ex. 3.1.
Conlusion 3 : $\textrm{Im}(f_3)=\{(a,b,c)\in\mathbb{R}^3/a+b+c=0\}=\{(a,b,-a-b),a,b\in\mathbb{R}\}=\textrm{Vect}\{(1,0,-1),(0,1,-1)\}$.
$\textrm{Im}(f_3)$ a pour base $\mathcal{B'}_3=(v_1,v_2)$ où $v_1=(1,0,-1)$ et $v_2=(0,1,-1)$ car ces vecteurs sont non colinéaires.
- $f_4$ est surjective d'après l'exercice 3.1 donc $\textrm{Im}(f_4)=\mathbb{R^3}$.
On peut donc choisir la base canonique de $\mathbb{R}^3$ comme base de $\textrm{Im}(f_3)$.
$$(a,b,c)\in\textrm{Im}(f_5)\Leftrightarrow \exists (x,y,z)\in\mathbb{R}^3/\begin{cases}y-z&=a\\-y+z&=b\\-y+z&=c\end{cases}\Leftrightarrow \exists (x,y,z)\in\mathbb{R}^3/\begin{cases}y-z&=a\\0&=a+b\\0&=a+c\end{cases}$$
- $f_5$ est non surjective et $\dim(\textrm{Im}(f_5))=1$ d'après l'ex. 3.1.
Conlusion 5 : $\textrm{Im}(f_5)=\{(a,b,c)\in\mathbb{R}^3/a+b=0\text{ et }a+c=0\}=\{(a,-a,-a),a\in\mathbb{R}\}=\textrm{Vect}\{(1,-1,-1)\}$.
$\textrm{Im}(f_5)$ a pour base $\mathcal{B'}_5=(v_1)$ où $v_1=(1,-1,-1)$.
Réponse exercice 3.3 :
- $A_1=\begin{pmatrix}1&-1\\1&0\\1&2\end{pmatrix}$.
$\textrm{rg}(A)=2=\dim(\textrm{Im}(f_1))$ donc, d'après la formule du rang : $\dim(\textrm{Ker}(f_1))=2-2=0$.
Il découle immédiatement que $f_1$ est non surjective ($\textrm{Im}(f_1)\neq\mathbb{R}^3$) et injective ($\textrm{Ker}(f_1)=\{0_{\mathbb{R}^3}\}$).
- $A_2=\begin{pmatrix}2&1&0\\1&0&-1\end{pmatrix}$.
$\textrm{rg}(A_2)=2=\dim(\textrm{Im}(f_2))$ donc, d'après la formule du rang : $\dim(\textrm{Ker}(f_2))=3-2=1$.
Il découle immédiatement que $f_2$ est surjective ($\textrm{Im}(f_2)=\mathbb{R}^2$) et non injective ($\textrm{Ker}(f_2)\neq\{0_{\mathbb{R}^3}\}$).
- $A_3=\begin{pmatrix}1&-1&0\\0&1&-1\\-1&0&1\end{pmatrix}$.
$\textrm{rg}(A_3)=2=\dim(\textrm{Im}(f_3))$ donc, d'après la formule du rang : $\dim(\textrm{Ker}(f_3))=3-2=1$.
Remarque : Notez en passant que la somme des colonnes de $A_3$ vaut $0$. Autrement dit :
$$f_3(e_1)+f_3(e_2)+f_3(e_3)=f_3(e_1+e_2+e_3)=0$$
D'où $e_1+e_2+e_3\in\textrm{Ker}(f_3)$ et donc $\textrm{Ker}(f_3)$ qui est de dimension $1$ a pour base $(u_1)$ où $u_1=(1,1,1)$.
Il découle par ailleurs immédiatement que $f_3$ est non surjective ($\textrm{Im}(f_3)\neq\mathbb{R}^3$) et non injective ($\textrm{Ker}(f_3)\neq\{0_{\mathbb{R}^3}\}$).
- $A_4=\begin{pmatrix}0&3&2\\-1&0&0\\0&4&3\end{pmatrix}$.
$\textrm{rg}(A_3)=3=\dim(\textrm{Im}(f_3))$ donc $A_4$ est inversible.
Il découle immédiatement que $f_4$ est bijective et donc surjective et injective.
- $A_5=\begin{pmatrix}0&1&-1\\0&-1&1\\0&-1&1\end{pmatrix}$.
$\textrm{rg}(A_5)=1=\dim(\textrm{Im}(f_5))$ donc, d'après la formule du rang : $\dim(\textrm{Ker}(f_5))=3-1=2$.
Il découle immédiatement que $f_3$ est non surjective ($\textrm{Im}(f_5)\neq\mathbb{R}^3$) et non injective ($\textrm{Ker}(f_3)\neq\{0_{\mathbb{R}^3}\}$).
Remarque : Il est aisé de trouver une base de l'image si on rappelle que, si $\mathcal{B}_3=(e_1,e_2,e_3)$ désigne la base canonique de $\mathbb{R}^3$, alors :
$$f_5(\mathcal{B}_3)=(f(e_1),f(e_2),f(e_3))=(0_{\mathbb{R}^3},(1,-1,-1),(-1,1,1))$$ est une famille génératrice de $\textrm{Im}(f_5)$
Sachant que $\dim(\textrm{Im}(f_5))=1$ on déduit immédiatement qu'une base de $\textrm{Im}(f_5)$ est $\mathcal{B'}_5=(v_1)$ où $v_1=(1,-1,-1)$. On retrouve le résultat de l'exercice 3.2...