glever - inverse d'un faisceau de matrices

Calling Sequence

[Bfs,Bis,chis]=glever(E,A [,s])

Parameters

• E, A : matrices carrées réelles de même dimensions
• s : chaîne de caractères (indéterminée des polynômes, ' s ' par défaut )
• Bfs,Bis : deux matrices polynomiales
• chis : polynôme

Description

Calcul de

(s*E-A)^-1

par l'algorithme généralisé de Leverrier pour un faisceau de matrices.

      (s*E-A)^-1 = (Bfs/chis) - Bis.
      

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 .

Attention

Cette fonction utilise cleanp pour simplifier Bfs,Bis et chis .

Examples

    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)
    

See Also

rowshuff ,   det ,   invr ,   coffg ,   pencan ,   penlaur ,  

Author

F. D. (1988)