Abstract:Quaternions are a non-commutative associative number system that extends complex numbers, first described by Hamilton in 1843. We present algorithms for solving the eigenvalue problem for arrowhead and DPRk (diagonal-plus-rank-k) matrices of quaternions. The algorithms use the Rayleigh Quotient Iteration with double shifts (RQIds), Wielandt's deflation technique and the fact that each eigenvector can be computed in O(n) operations. The algorithms require O(n2) floating-point operations, n being the order of the matrix. The algorithms are backward stable in the standard sense and compare well to the standard QR method in terms of speed and accuracy. The algorithms are elegantly implemented in Julia, using its polymorphism feature.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to efficiently compute the quaternion eigenvalue decomposition of arrow matrices and diagonal plus rank - k matrices (DPRk matrices). Quaternions are a non - commutative associative number system, which extends the concept of complex numbers and was first described by Hamilton in 1843. Due to the non - commutative nature of quaternions, traditional eigenvalue - solving methods need to be modified because quaternion matrices do not have characteristic polynomials. Therefore, the author proposes a method to solve this problem by using the Rayleigh quotient iteration method (RQI) with double - shift techniques, Wielandt's deflation technique, and the fact that each eigenvector can be computed in \(O(n)\) operations. These algorithms are backward - stable in the standard sense and are comparable to the standard QR method in terms of speed and accuracy. Moreover, these algorithms implement elegant polymorphic functionality in the Julia language. ### Specific Problems and Solutions 1. **Problem Background**: - Quaternions are a non - commutative associative number system that extends complex numbers. - The eigenvalue problem of matrices becomes complicated when extended to quaternion matrices because quaternion matrices do not have characteristic polynomials. - Traditional methods need to be modified to adapt to the non - commutativity of quaternions. 2. **Research Objectives**: - Propose efficient algorithms to compute the quaternion eigenvalue decomposition of arrow matrices and DPRk matrices. - Ensure that the algorithms are backward - stable in the standard sense. - The algorithms perform well in terms of speed and accuracy compared to the standard QR method. 3. **Main Methods**: - **Rayleigh Quotient Iteration Method (RQI) with Double - Shift Techniques**: Use the Rayleigh quotient iteration method (RQI) and double - shift techniques to accelerate convergence. - **Wielandt's Deflation Technique**: Gradually reduce the size of the matrix through deflation techniques, thereby simplifying the calculation of eigenvalues and eigenvectors. - **Fast Computation of Inverse Matrices**: Utilize the special structures of arrow matrices and DPRk matrices to quickly compute inverse matrices. 4. **Algorithm Efficiency**: - Each eigenvector can be computed in \(O(n)\) operations. - The entire algorithm is completed in \(O(n^{2})\) floating - point operations. 5. **Implementation**: - The algorithms are implemented in the Julia language, taking advantage of its polymorphic functionality to ensure the elegance and efficiency of the code. ### Conclusion This paper successfully proposes an efficient and stable algorithm for computing the quaternion eigenvalue decomposition of arrow matrices and DPRk matrices. These algorithms not only have good performance theoretically but also perform excellently in practical applications, especially when dealing with large - scale matrices. By using the Rayleigh quotient iteration method, double - shift techniques, and Wielandt's deflation technique, this method can significantly improve computational efficiency while ensuring accuracy.

Fast Eigenvalue Decomposition of Arrowhead and Diagonal-Plus-Rank-k Matrices of Quaternions

On Rayleigh Quotient Iteration for Dual Quaternion Hermitian Eigenvalue Problem

On rank‐revealing QR factorizations of quaternion matrices

Randomized Large-Scale Quaternion Matrix Approximation: Practical Rangefinders and One-Pass Algorithm

Applications of Dual Complex Adjoint Matrix in Eigenvalue Computation of Dual Quaternion Hermitian Matrix

Quaternion matrix decomposition and its theoretical implications

Q-less QR Decomposition in Inner Product Spaces

Block Diagonalization of Quaternion Circulant Matrices with Applications

'Diophantine' and factorization properties of finite orthogonal polynomials in the Askey scheme

2D Eigenvalue Problem II: Rayleigh Quotient Iteration and Applications

Fast and backward stable computation of the eigenvalues and eigenvectors of matrix polynomials

Generalizing Frobenius Inversion to Quaternion Matrices

Real dqds for the nonsymmetric tridiagonal eigenvalue problem

Randomized Low Rank Approximation for Nonnegative Pure Quaternion Matrices

Solutions of the Partially Described Inverse Quadratic Eigenvalue Problem

Adaptive Parallelizable Algorithms for Interpolative Decompositions via Partially Pivoted LU

Multilinear analysis of quaternion arrays: theory and computation

Low-Rank Kernel Matrix Approximation Using Skeletonized Interpolation With Endo- or Exo-Vertices

Fast Eigen Decomposition for Low-Rank Matrix Approximation

Effective Methods of QR-Decompositions of Square Complex Matrices by Fast Discrete Signal-Induced Heap Transforms

The Complexity of Diagonalization