Relative Perturbation Bounds for the Unitary Polar Factor
Abstract: Let B be an m x n (m >= n) complex matrix. It is known that there is a unique polar decomposition B = QH, where Q*Q = I, the n x n identity matrix, and H is positive definite, provided B has full column rank. This paper addresses the following question: how much may Q change if B is perturbed to ~B = D1BD2? Here D1 and D2 are two nonsingular matrices and close to the identities of suitable dimensions.
Known perturbation bounds for complex matrices indicate that in the worst case, the change in Q is proportional to the reciprocal of the smallest singular value of B. In this paper, we will prove that for the above mentioned perturbations to B, the change in Q is bounded only by the distances from D1 and D2 to identities!
As an application, we will consider perturbations for one-side scaling, i.e., the case when G = D * B is perturbed to ~G = D * ~B, where D is usually a nonsingular diagonal scaling matrix but for our purpose we do not have to assume this, and B and ~B are nonsingular.