inverse d'un faisceau de matrices
[Bfs,Bis,chis]=glever(E,A [,s])
matrices carrées réelles de même dimensions
chaîne de caractères (indéterminée des polynômes, 's
' par défaut )
deux matrices polynomiales
polynôme
Calcul de
(s*E-A)^-1
par l'algorithme généralisé de Leverrier pour un faisceau de matrices.
chis
= polynôme caractéristique (à une constante multiplicative près).
Bfs
= matrice polynomiale de numérateurs
Bis
= matrice polynomiale ( - développement de (s*E-A)^-1
à l'infini).
Noter le signe - devant Bis
.
Cette fonction utilise cleanp
pour simplifier Bfs,Bis
et chis
.
s=%s;F=[-1,s,0,0;0,-1,0,0;0,0,s-2,0;0,0,0,s-1]; [Bfs,Bis,chis]=glever(F) inv(F)-((Bfs/chis) - Bis) | ![]() | ![]() |