Algorithms on the Stiefel manifold for joint diagonalisation
Algorithms on the Stiefel manifold for joint diagonalisation
复制标题
用于联合对角化的 Stiefel 流形算法
DOI:
10.1109/icassp.2002.5744893
复制
发表时间:
2002
期刊:
影响因子:
--
通讯作者:
G. Hori
中科院分区:
文献类型:
--
作者:
M. Nikpour;J. Manton;G. Hori
This paper presents a novel approach to the unitary joint diagonalisation problem. We approach the problem as an optimisation problem over the manifold of unitary matrices and obtain steepest descent and Newton algorithms. Both algorithms converge to the same solution as the well known Joint Approximate Diagonalisation of Eigen-matrices (JADE) algorithm. However, our Newton algorithm achieves a quadratic convergence rate almost immediately if suitably initialised, whereas JADE can (in some cases) take up to 30 sweeps before it starts converging to the solution at a quadratic rate. Simulations illustrate that our steepest descent algorithm can be used to find a suitable starting point for our Newton algorithm which then refines our estimate at a quadratic rate. This two stage algorithm requires less iterations overall than JADE, with the cost per iteration being the same as JADE. Our methods, therefore, offer a novel alternative to JADE.