projeção espectral
[Q,M]=psmall(A,thres,flag)
matriz quadrada de reais
número real
string ('c' ou
'd')
matrizes de reais
Projeção sobre auto-subespaço associado com autovalores com parte
real < thres (flag='c') ou com
módulo < thres (flag='d').
A projeção é definda por Q*M,
Q é de posto-coluna cheio, M é de
posto-linha cheio e M*Q=eye.
Se flag='c', os autovalores de
M*A*Q = autovalores de A com parte
real < thres.
Se flag='d', os autovalores de
M*A*Q = autovalores de A com
magnitude < thres.
Se flag='c' e se [Q1,M1] =
fatoração em posto cheio (fullrf) de
eye()-Q*Mentão os autovalores de
M1*A*Q1 = autovalores deA com parte
real >= thres.
Se flag='d' e se [Q1,M1] =
fatoração em posto cheio (fullrf) de
eye()-Q*M então os autovalores de
M1*A*Q1 = autovalores de A com
magnitude >= thres.
Esta função é baseada na forma de Schur ordenada (Função do
Scilabschur).