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

Technische Samenvatting: Iteratieve methoden voor het berekenen van de Moore–Penrose-inverse van split-quaternion matrices

Probleemstelling
Het artikel behandelt de berekening van de Moore–Penrose (MP) inverse voor matrices over de algebra van split-quaternionen ($\mathbb{Q}_s$). In tegenstelling tot standaard Hamilton-quaternionen bezitten split-quaternionen een indefinite kwadratische vorm, wat leidt tot het bestaan van niet-triviale nuldelers en nilpotente elementen. Deze algebraïsche kenmerken compliceren de definitie van inversen en de toepassing van de standaard spectrale theorie. Specifiek voorkomt de indefinitie van de kwadratische vorm de directe overdracht van spectrale resultaten van reële of complexe matrices, en vereist de aanwezigheid van nuldelers een zorgvuldige behandeling van singuliere en rang-gebrekkige gevallen. De auteurs streven ernaar een robuust kader vast te stellen voor het definiëren van de MP-inverse in deze setting en efficiënte iteratieve algoritmen te ontwikkelen voor de berekening ervan, met name voor rechthoekige en rang-gebrekkige matrices.

Methodologie
De kern van de methodologische bijdrage is de vestiging van een consistent kader gebaseerd op een $2 \times 2$ reële representatie ($\sigma$) en de $i$-geconjugeerde transponering ($H$).

1. Reële representatie: De auteurs definiëren een bijectieve afbeelding $\sigma: \mathbb{Q}_s \to \mathbb{R}^{2 \times 2}$ die vermenigvuldiging behoudt en de $i$-geconjugeerde transponering afbeeldt op de standaard reële transponering. Voor een matrix $A \in \mathbb{Q}_s^{m \times n}$ maakt de reële representatie $A^\sigma \in \mathbb{R}^{2m \times 2n}$ het mogelijk om de MP-inverse direct te definiëren via de reële MP-inverse: $A^\dagger = \sigma^{-1}((A^\sigma)^\dagger)$. Deze definitie verheldert de behandeling van nuldelers en waarborgt de existentie en uniciteit van de inverse.
2. Newton–Schulz iteratie: Het artikel analyseert de Newton–Schulz iteratie voor het berekenen van de MP-inverse. Twee implementaties worden gepresenteerd:
   • Embedded (ingebed): Het toepassen van de standaard reële Newton–Schulz iteratie op de reële representatie $A^\sigma$.
   • Native (natuurlijk): Een directe iteratie in split-quaternion rekenkunde ($X_{k+1} = X_k(2I - AX_k)$).
     De auteurs bewijzen dat deze twee sequenties exact equivalent zijn onder de representatie $\sigma$.
3. Convergentieanalyse: Met behulp van een dunne Singular Value Decomposition (SVD) van de reële representatie leiden de auteurs exacte formules af voor de iteraties en residuen. Zij stellen vast dat de iteratie convergeert naar de MP-inverse voor matrices van willekeurige rang (inclusief rang-gebrekkige gevallen), mits de initiële schaleringsparameter $\alpha$ voldoet aan $0 < \alpha < 2/\|A^\sigma\|_2^2$. De analyse biedt exacte relaties tussen het residu en de fout, waarbij een kwadratische convergentie van het residu wordt aangetoond.
4. Polynomiale initialisatie: Om de convergentie te versnellen, stellen de auteurs een laag-graads polynomiale initialisatie voor, geïnspireerd door de inverse recursie van Souriau. Dit houdt in dat er een polynoom wordt aangepast aan de inverse van de genormaliseerde Gram-factor ($A^\sigma (A^\sigma)^T$ of $(A^\sigma)^T A^\sigma$) met behulp van een kleinste-kwadraten (LS) benadering. De methode bevat een spectrale acceptatietest om te garanderen dat het initiële residu binnen de convergentieregio ligt, met een veilige terugval naar een geschaalde transponering indien de test faalt.

Belangrijkste Bijdragen

• Consistent Kader: Het artikel definieert de MP-inverse voor split-quaternion matrices via de reële representatie en de $i$-geconjugeerde transponering, waarbij existentie, uniciteit en basiseigenschappen (inclusclusief een scalaire formule voor nuldelers) worden bewezen.
• Convergentietheorie: Een volledige convergentieanalyse voor de Newton–Schulz iteratie wordt geboden voor rechthoekige en rang-gebrekkige matrices. De auteurs leiden exacte foutrecurrenties en residu-fout-identiteiten af, waardoor de restricties op volledige rang die in eerdere analyses van Hamilton-quaternion matrices werden gevonden, worden opgeheven.
• Equivalentie van Implementaties: Het artikel demonstreert de exacte equivalentie tussen de embedded (reële) en native (split-quaternion) Newton–Schulz iteraties, wat het gebruik van native rekenkunde zonder verlies van theoretische garanties valideert.
• LS–Gram Initialisatie: Een nieuwe polynomiale initialisatie wordt geïntroduceerd die zich aanpast aan het specifieke spectrum van de Gram-matrix. Er wordt aangetoond dat dit het aantal benodigde Newton–Schulz iteraties vermindert vergeleken met een geschaalde transponering.
• CUR Approximatie: De auteurs breiden de cross- en CUR-approximatietheorie uit naar split-quaternion matrices. Zij bewijzen dat $C^\dagger A R^\dagger$ de Frobenius-optimale middelste factor is voor gefixeerde gesamplede rijen en kolommen. Bovendien stellen zij een rang-onthullende conditie vast waaronder de goedkopere cross-factor $W^\dagger$ (waarbij $W$ de intersectiematrix is) een exacte reconstructie oplevert.

Resultaten
Numerieke experimenten uitgevoerd in MATLAB en Python valideren de theoretische bevindingen:

• Nauwkeurigheid: De iteratieve methoden bereiken een nauwkeurigheid die vergelijkbaar is met directe SVD-gebaseerde berekeningen (relatieve fouten rond $10^{-11}$ tot $10^{-14}$).
• Efficiëntie: Voor matrices van gemiddelde grootte en dichtheid is de Newton–Schulz iteratie aanzienlijk sneller dan directe SVD-berekening (bijv. ~3,7x sneller voor een $200 \times 250$ matrix).
• Initialisatieprestaties: De LS–Gram polynomiale initialisatie vermindert consequent het aantal Newton–Schulz iteraties vergeleken met een geschaalde transponering. In experimenten met condities van tot $10^4$ verminderde de LS–Gram methode het aantal iteraties met ongeveer 20% ten opzichte van de geschaalde transponering en bleef deze competitief met een Chebyshev-polynomiale initialisatie die exacte spectrale grenzen vereiste.
• CUR Approximatie: De experimenten bevestigen dat de cross-approximatie met behulp van $W^\dagger$ een exacte reconstructie kan bereiken voor low-rank matrices wanneer de rang-onthullende conditie wordt vervuld, waarbij de iteratieve berekening van de kern-inverse aanzienlijk sneller is dan directe methoden.

Betekenis en Claims
Het artikel claimt het eerste consistente kader te bieden voor de MP-inverse van split-quaternion matrices dat expliciet nuldelers en rang-gebrek behandelt via een reële representatie. Door de analyse van Newton–Schulz iteraties van reële matrices over te dragen naar de split-quaternion setting, vestigen de auteurs rigoureuze convergentiegaranties voor matrices van willekeurige rang. De voorgestelde LS–Gram initialisatie biedt een praktische verbetering ten opzichte van standaard initialisatiestrategieën, waarbij de computationele kosten worden verlaagd zonder voorafgaande kennis van de spectrale grenzen te vereisen (in tegen tegenstelling tot Chebyshev-methoden). Ten slotte biedt de uitbreiding van de CUR-approximatie resultaten een theoretische basis voor low-rank benaderingen in split-quaternion signaalverwerking en fysica-toepassingen, waarbij de optimaliteit van specifieke factoren en condities voor exacte reconstructie wordt gekarakteriseerd. De auteurs merken op dat hoewel de methoden efficiënt zijn, de primaire bijdrage de theoretische unificatie en de ontwikkeling van robuuste iteratieve solvers is, waarbij toekomstig werk voor ijle (sparse) matrices en de Drazin-inverse wordt gesuggereerd.