pencan

Parâmetros

Fs

um feixe regular s*E-A

E,A

duas matrizes quadradas de reais

Q,M

duas matrizes não-singulares de reais

i1

inteiro

Descrição

Dado o feixe regular Fs=s*E-A, pencan retorna as matrizes Q e M tais que M*(s*E-A)*Q está na forma "canônica".

M*E*Q é uma matriz de blocos

[I,0;
 0,N]

com N nilpotente e i1 = tamanho da matriz acima I.

M*A*Q é uma matriz de blocos:

[Ar,0;
 0,I]

Exemplos

F=randpencil([],[1,2],[1,2,3],[]);
F=rand(6,6)*F*rand(6,6);
[Q,M,i1]=pencan(F);
W=clean(M*F*Q)
roots(det(W(1:i1,1:i1)))
det(W($-2:$,$-2:$))

Ver Também

• glever
• penlaur
• rowshuff