An NPDo Approach for Principal Joint Block Diagonalization

Matrix joint block-diagonalization (jbd) frequently arises from diverse applications such as independent component analysis, blind source separation, and common principal component analysis (CPCA), among others. Particularly, CPCA aims at joint diagonalization, i.e., each block size being 1-by-1. This paper is concerned with principal joint block-diagonalization (pjbd), which aim to achieve two goals: 1) partial joint block-diagonalization, and 2) identification of dominant common block-diagonal parts for all involved matrices. This is in contrast to most existing methods, especially the popular ones based on Givens rotation, which focus on full joint diagonalization and quickly become impractical for matrices of even moderate size (300-by-300 or larger). An NPDo approach is proposed and it is built on a nonlinear polar decomposition with orthogonal polar factor dependency that characterizes the solutions of the optimization problem designed to achieve pjbd, and it is shown the associated SCF iteration is globally convergent to a stationary point while the objective function increases monotonically during the iterative process. Numerical experiments are presented to illustrate the effectiveness of the NPDo approach and its superiority to Givens rotation-based methods. MSC Classification: 62H25 , 65F30 , 65K05 , 90C26