Iterative Methods for Computing the Moore--Penrose Inverse of Split-Quaternion Matrices with Applications

Resumo Técnico: Métodos Iterativos para o Cálculo da Inversa de Moore–Penrose de Matrizes de Split-Quaternions

Enunciado do Problema
O artigo aborda o cálculo da inversa de Moore–Penrose (MP) para matrizes sobre a álgebra dos split-quaternions ($\mathbb{Q}_s$). Diferente dos quaternions de Hamilton padrão, os split-quaternions possuem uma forma quadrática indefinida, o que leva à existência de divisores de zero não triviais e elementos nilpotentes. Essas características algébricas complicam a definição de inversas e a aplicação da teoria espectral padrão. Especificamente, a natureza indefinida da forma quadrática impede a transferência direta de resultados espectrais de matrizes reais ou complexas, e a presença de divisores de zero exige um tratamento cuidadoso de casos singulares e de posto deficiente. Os autores visam estabelecer um framework robusto para definir a inversa MP neste cenário e desenvolver algoritmos iterativos eficientes para seu cálculo, particularmente para matrizes retangulares e de posto deficiente.

Metodologia
A principal contribuição metodológica é o estabelecimento de um framework consistente baseado em uma representação real $2 \times 2$ ($\sigma$) e na transposta conjugada-$i$ ($H$).

1. Representação Real: Os autores definem um mapa bijetivo $\sigma: \mathbb{Q}_s \to \mathbb{R}^{2 \times 2}$ que preserva a multiplicação e mapeia a transposta conjugada-$i$ para a transposta real padrão. Para uma matriz $A \in \mathbb{Q}_s^{m \times n}$, o seu representante real $A^\sigma \in \mathbb{R}^{2m \times 2n}$ permite que a inversa MP seja definida diretamente via a inversa MP real: $A^\dagger = \sigma^{-1}((A^\sigma)^\dagger)$. Esta definição esclarece o tratamento de divisores de zero e garante a existência e unicidade da inversa.
2. Iteração de Newton–Schulz: O artigo analisa a iteração de Newton–Schulz para calcular a inversa MP. Duas implementações são apresentadas:
   • Embebida (Embedded): Aplicando a iteração de Newton–Schulz real padrão ao representante real $A^\sigma$.
   • Nativa (Native): Uma iteração direta na aritmética de split-quaternions ($X_{k+1} = X_k(2I - AX_k)$).
     Os autores provam que estas duas sequências são exatamente equivalentes sob a representação $\sigma$.
3. Análise de Convergência: Utilizando uma Decomposição de Valores Singulares (SVD) "fina" do representante real, os autores derivam fórmulas exatas para as iterações e resíduos. Eles estabelecem que a iteração converge para a inversa MP para matrizes de posto arbitrário (incluindo casos de posto deficiente), desde que o parâmetro de escala inicial $\alpha$ satisfaça $0 < \alpha < 2/\|A^\sigma\|_2^2$. A análise fornece relações exatas entre o resíduo e o erro, mostrando convergência quadrática do resíduo.
4. Inicialização por Polinômio: Para acelerar a convergência, os autores propõem uma inicialização polinomial de baixo grau inspirada na recursão de inversa de Souriau. Isso envolve o ajuste de um polinômio à inversa do fator Gram normalizado ($A^\sigma (A^\sigma)^T$ ou $(A^\sigma)^T A^\sigma$) usando uma abordagem de mínimos quadrados (LS). O método inclui um teste de aceitação espectral para garantir que o resíduo inicial esteja dentro da região de convergência, com um fallback seguro para uma transposta escalonada caso o teste falhe.

Principais Contribuições

• Framework Consistente: O artigo define a inversa MP para matrizes de split-quaternions via representação real e a transposta conjugada-$i$, provando existência, unicidade e propriedades básicas (incluindo uma fórmula escalar para divisores de zero).
• Teoria de Convergência: Uma análise de convergência completa para a iteração de Newton–Schulz é fornecida para matrizes retangulares e de posto deficiente. Os autores derivam recorrências de erro exatas e identidades resíduo-erro, removendo as restrições de posto total encontradas em análises anteriores de matrizes de quaternions de Hamilton.
• Equivalência de Implementações: O artigo demonstra a equivalência exata entre as iterações de Newton–Schulz embebida (real) e nativa (split-quaternion), validando o uso de aritmética nativa sem perda de garantias teóricas.
• Inicialização LS–Gram: Uma inicialização polinomial LS–Gram inovadora é introduzida, que se adapta ao espectro específico da matriz Gram. Demonstra-se que ela reduz o número de iterações necessárias em comparação com a inicialização de transposta escalonada.
• Aproximação CUR: Os autores estendem a teoria de aproximação cross e CUR para matrizes de split-quaternions. Eles provam que $C^\dagger A R^\dagger$ é o fator central Frobenius-ótimo para linhas e colunas amostradas fixas. Além disso, estabelecem uma condição de revelação de posto sob a qual o fator cross mais barato $W^\dagger$ (onde $W$ é a matriz de interseção) produz uma reconstrução exata.

Resultados
Experimentos numéricos conduzidos em MATLAB e Python validam as descobertas teóricas:

• Acurácia: Os métodos iterativos alcançam acurácia comparável ao cálculo direto via SVD (erros relativos em torno de $10^{-11}$ a $10^{-14}$).
• Eficiência: Para matrizes densas de tamanho médio, a iteração de Newton–Schulz é significativamente mais rápida que o cálculo direto por SVD (ex: ~3.7x mais rápida para uma matriz $200 \times 250$).
• Desempenho de Inicialização: A inicialização polinomial LS–Gram consistentemente reduz o número de iterações de Newton–Schulz em comparação com a transposta escalonada. Em experimentos com números de condição até $10^4$, o método LS–Gram reduziu as iterações em aproximadamente 20% em relação à transposta escalonada e permaneceu competitivo com uma inicialização de polinômio de Chebyshev, que exigia limites espectrais exatos.
• Aproximação CUR: Os experimentos confirmam que a aproximação cross usando $W^\dagger$ pode alcançar reconstrução exata para matrizes de baixo posto quando a condição de revelação de posto é atendida, sendo que o cálculo iterativo da inversa do núcleo é substancialmente mais rápido que os métodos diretos.

Significância e Alegações
O artigo afirma fornecer o primeiro framework consistente para a inversa MP de matrizes de split-quaternions que lida explicitamente com divisores de zero e deficiência de posto através de uma representação real. Ao transferir a análise das iterações de Newton–Schulz de matrizes reais para o cenário de split-quaternions, os autores estabelecem garantias de convergência rigorosas para matrizes de posto arbitrário. A proposta de inicialização LS–Gram oferece uma melhoria prática sobre as estratégias de inicialização padrão, reduzindo o custo computacional sem exigir conhecimento prévio dos limites espectrais (ao contrário dos métodos de Chebyshev). Finalmente, a extensão dos resultados de aproximação CUR fornece uma base teórica para aproximações de baixo posto em processamento de sinais e aplicações de física com split-quaternions, caracterizando a otimalidade de fatores específicos e condições para reconstrução exata. Os autores observam que, embora os métodos sejam eficientes, a principal contribuição é a unificação teórica e o desenvolvimento de solvers iterativos robustos, com trabalhos futuros sugeridos para matrizes esparsas e a inversa de Drazin.