This version is available at https://strathprints.strath.ac.uk/57580/ Strathprints is designed to allow users to access the research output of the University of Strathclyde. Unless otherwise explicitly stated on the manuscript, Copyright © and Moral Rights for the papers on this site are retained by the individual authors and/or other copyright owners. Please check the manuscript for details of any other licences that may have been applied. You may not engage in further distribution of the material for any profitmaking activities or any commercial gain. You may freely distribute both the url (https://strathprints.strath.ac.uk/) and the content of this paper for research or private study, educational, or not-for-profit purposes without prior permission or charge.Any correspondence concerning this service should be sent to the Strathprints administrator: strathprints@strath.ac.ukThe Strathprints institutional repository (https://strathprints.strath.ac.uk) is a digital archive of University of Strathclyde research outputs. It has been developed to disseminate open access research outputs, expose data about those outputs, and enable the management and persistent access to Strathclyde's intellectual output. Abstract-Recently a selection of sequential matrix diagonalisation (SMD) algorithms have been introduced which approximate polynomial eigenvalue decomposition of parahermitian matrices. These variants differ only in the search methods that are used to bring energy onto the zero-lag. Here we analyse the search methods in terms of their computational complexities for different sizes of parahermitian matrices which are verified through simulated execution times. Another important factor for these search methods is their ability to transfer energy. Simulations show that the more computationally complex search methods transfer a greater proportion of the off-diagonal energy onto the zero-lag over a selected range of parahermitian matrix sizes. Finally we compare the real time convergence of the search methods as part of their respective SMD algorithms. The real time convergence experiments indicate that despite taking a longer time to compute each iteration the more complex algorithms that transfer more energy converge faster in real time.
I. INTRODUCTIONSequential matrix diagonalisation encompasses a family of iterative algorithms that can factorise a parahermitian matrix into an approximate polynomial matrix eigenvalue decomposition (PEVD). The PEVD extends the wide-ranging utility of the EVD from narrowband to broadband problems, and iterative PEVD algorithms have in the past found use in optimal subband [1] Parahermitian matrices arise e.g. by including an explicit lag τ into the space-time covariance