inverso do feixe de matrizes
[Bfs,Bis,chis]=glever(E,A [,s])
duas matrizes de reais quadradas de igual dimensão
string (o valor padrão é 's')
duas matrizes de polinômios
polinômio
Computação de
(s*E-A)^-1
pelo algoritmo generalizado de Leverrier para um feixe de matrizes.
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.
Esta função usa cleanp para simplificar
Bfs,Bis e 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) | ![]() | ![]() |