A structured method for the real quadratic eigenvalue problem for specific gyroscopic systems

This study examines a specific numerical approach that computes the eigenvalues (normal modes) of a Quadratic Eigenvalue Problem (QEP) of the form (lambda 2I + lambdaB + C) x = 0 where B is constrained to a real skew-symmetric matrix and C is constrained to a real symmetric positive definite matrix. A widely used linearization of this QEP is the companion matrix A = -B-C IO . The goal is to find an algorithm method which diagonalizes matrix A without contaminating the (2,2) zero block. Once this algorithm is developed, the study measures the eigenvalue error bounds and compare its efficiency to the standard symmetric QR workhorse. Also, this approach preserves the structure of the error matrix in the same form as the QEP. In ensuring that the error matrix structure is a QEP, this algorithm provides fertile ground for future analysis in sensitivity and perturbation errors in the algorithm's eigenvalues. This study concludes that the algorithm appears to have a reasonable error bound; and it is more cost efficient in finding the eigenvalues then the symmetric QR algorithm.