Exercices d’algèbre linéaire tombés à Centrale PC en 2020. Cherchez d’abord à partir de l’énoncé, puis ouvrez une indication ou le corrigé lorsque vous en avez besoin.
Les polynômes P1,…,Pp+2sont p+2 vecteurs de l’espaceRp[X], qui est de dimension p+1. La famille est donc liée : il existe des réelsα1,…,αp+2, non tous nuls, tels que
α1P1+⋯+αp+2Pp+2=0.
La même combinaison linéaire des lignes deBp est donc la ligne nulle. Les lignes sont liées, d’où
Pour C, on inverse l’ordre des lignes. On obtient alors un Vandermonde évalué en0,1,…,p. Le renversement dep+1 lignes apporte le signe(−1)⌊(p+1)/2⌋, donc
det(C)=(−1)⌊(p+1)/2⌋j=1∏pj!.
Finalement les produits de factorielles se simplifient :
det(Ap)=(−1)⌊(p+1)/2⌋(p!)p+1
Exercice 2—Matrices à diagonale dominante et déterminant
Exercice tombé en 2020 · RMS 1061 · Python
DifficultéÀ définir
Énoncé
Soit A=(ai,j)1≤i,j≤n∈Mn(C). Pour i∈[[1,n]], on note
ρi(A)=∣ai,i∣−j=i∑∣ai,j∣.
Une matrice est dite à diagonale dominante siρi(A)>0 pour touti.
Montrer qu’une matrice à diagonale dominante est inversible.
Étudier la réciproque.
Programmer en Python une fonction d’argument une matrice à diagonale dominante et renvoyant∏k=1nρk(A).
Comparer sur des exemples∏k=1nρk(A) et∣det(A)∣. Conjecture ?
Prouver la conjecture.
Indication›
Pour l’inversibilité, partir de AX=0et choisir une coordonnée de Xde module maximal.
Pour la dernière question, normaliser chaque ligne en la divisant parρi(A) et utiliser le fait qu’une matrice non inversible ne peut pas être à diagonale dominante.
Corrigé›
1. Théorème d’Hadamard
Supposons AX=0 avecX=(x1,…,xn)T. Choisissonsi0 tel que
∣xi0∣=1≤i≤nmax∣xi∣.
La ligne i0 de l’égalité donne
ai0,i0xi0=−j=i0∑ai0,jxj.
Donc
∣ai0,i0∣∣xi0∣≤∣xi0∣j=i0∑∣ai0,j∣.
Autrement dit,
ρi0(A)∣xi0∣≤0.
Comme ρi0(A)>0, on axi0=0, puis toutes les coordonnées de X sont nulles. Ainsi le noyau est réduit à zéro et
A est inversible.
2. Réciproque
La réciproque est fausse. Par exemple
A=(1021)
est inversible, mais la première ligne ne vérifie pas la domination diagonale puisque 1−2<0.
3. Fonction Python
import numpy as np
import numpy.linalg as alg
def rho(A):
n, res = len(A[0]), 1
for k in range(n):
s = abs(A[k, k])
for j in range(n):
if j != k:
s -= abs(A[k, j])
res *= s
return res
4. Conjecture
Les essais conduisent à conjecturer
k=1∏nρk(A)≤∣det(A)∣
5. Preuve
Commençons par une remarque. Siλ∈Sp(A), alors A−λIn n’est pas inversible. D’après la première question, elle n’est donc pas à diagonale dominante. Il existe donc un indicei0 tel que
∣ai0,i0−λ∣≤j=i0∑∣ai0,j∣.
Par l’inégalité triangulaire renversée,
∣ai0,i0∣−∣λ∣≤∣ai0,i0−λ∣≤j=i0∑∣ai0,j∣,
d’où
ρi0(A)≤∣λ∣.
Posons maintenant
A′=(ρi(A)ai,j)1≤i,j≤n.
On a, par construction,ρi(A′)=1 pour touti. Par la remarque précédente, toute valeur propre λ deA′ vérifie∣λ∣≥1. Par conséquent,
∣det(A′)∣=λ∈Sp(A′)∏∣λ∣≥1
en comptant les valeurs propres avec leurs multiplicités. Enfin, puisque chaque ligne de A a été divisée par ρi(A),