imagem inversa
[X,dim]=im_inv(A,B [,tol]) [X,dim,Y]=im_inv(A,B, [,tol])
duas matirzes de reais ou complexos com igual número de colunas
matriz quadrada ortogonal ou unitária de ordem igual ao número
de colunas de A
inteiro (dimensão do subespaço)
matriz ortogonal de ordem igual ao número de linhas de
A
e B
.
[X,dim]=im_inv(A,B)
computa(A^-1)(B)
,i.e, vetores cujas imagens através de
A
estão em Im(B
)
As dim
primeiras colunas de X
geram(A^-1)(B)
tol
é um limiar usado para testar a inclusão de
subespaço ; o valor padrão é tol = 100*%eps
. Se
Y
é retornado, então [Y*A*X,Y*B]
é
particionado como segue:
[A11,A12;0,A22]
,[B1;0]
onde B1
tem posto-linha cheio (igual a
posto(B)
) e A22
tem posto-coluna
cheio e tem dim
colunas.