An iterative algorithm converging to the principal subspace projection operator with applications

Abstract

This paper analyses a simple method for iteratively computing the projection operator whose range space is the principal subspace of a given matrix. The iteration appears to be numerically very stable and is proven to converge at a quadratic rate. Several applications are briefly explored, as are extensions to Jordan algebras and C*-algebras. © 2004 IEEE.