algoritmo de embaralhamento
[Ws,Fs1]=rowshuff(Fs, [alfa])
feixe quadrado de reais Fs = s*E-A
matriz de polinômios
feixe quadrado de reais F1s = s*E1 -A1
com
E1
não-singular
número real (alfa = 0
é o valor
padrão)
Algoritmo de embaralhamento: dado o feixe
Fs=s*E-A
, retorna Ws=W(s) (matriz quadrada de
polinômios) tal que:
Fs1 = s*E1-A1 = W(s)*(s*E-A)
é um feixe com
matriz E1
não-singular.
Isto é possível se, e só se, o feixe Fs = s*E-A
é
regular (i.e., invertível). O grau de Ws
é igual ao
índice do feixe.
Os pólos no infinito deFs
asão colocados para
alfa
e os zeros de Ws
estão em
alfa
.
Note que (s*E-A)^-1 = (s*E1-A1)^-1 * W(s) =
(W(s)*(s*E-A))^-1 *W(s)