penlaur

Parâmetros

Fs

um feixe regular s*E-A

E, A

duas matrizes quadradas de reais

Si,Pi,Di

três matrizes quadradas de reais

order

inteiro

Descrição

penlaur computa os primeiros coeficientes de Laurent de (s*E-A)^-1 no infinito.

(s*E-A)^-1 = ... + Si/s - Pi - s*Di + ... em s = infinito.

order = ordem da singularidade (ordem=índice-1).

O feixe de matrizes Fs=s*E-A deve ser invertível.

Para um feixe de índice 0, Pi, Di,... são zero e Si=inv(E).

Para um feixe de índice 1 (order=0),Di =0.

Para feixes de índices maiores, os termos -s^2 Di(2), -s^3 Di(3),... são dados por:

Di(2)=Di*A*Di, Di(3)=Di*A*Di*A*Di (até Di(order)).

Observação

Versão experimental: há problemas quando se tem mal-condicionamento deso*E-A

Exemplos

F=randpencil([],[1,2],[1,2,3],[]);
F=rand(6,6)*F*rand(6,6);[E,A]=pen2ea(F);
[Si,Pi,Di]=penlaur(F);
[Bfs,Bis,chis]=glever(F);
norm(coeff(Bis,1)-Di,1)

Ver Também

• glever
• pencan
• rowshuff