glever

Parâmetros

E, A

duas matrizes de reais quadradas de igual dimensão

s

string (o valor padrão é 's')

Bfs,Bis

duas matrizes de polinômios

chis

polinômio

Descrição

Computação de

(s*E-A)^-1

pelo algoritmo generalizado de Leverrier para um feixe de matrizes.

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

chis = polinômio característico (até uma constante multiplicativa).

Bfs = nmatriz de polinômios como numerador.

Bis = matriz de polinômios ( - expansão de (s*E-A)^-1 ao infinito).

Note o sinal - antes de Bis.

Cuidado

Esta função usa cleanp para simplificar Bfs,Bis e chis.

Exemplos

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)

Ver Também

• rowshuff
• det
• invr
• coffg
• pencan
• penlaur