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) | ![]() | ![]() |