Joint diagonalization of a set of positive (semi)-definite matrices has a wide range of analytical applications, such as estimation of common principal components, estimation of multiple variance components, and blind signal separation. However, when the eigenvectors of involved matrices are not the same, joint diagonalization is a computationally challenging problem. To the best of our knowledge, currently existing methods require at least $O(KN^3)$ time per iteration, when $K$ different $N \times N$ matrices are considered. We reformulate this optimization problem by applying orthogonality constraints and dimensionality reduction techniques. In doing so, we reduce the computational complexity for joint diagonalization to $O(N^3)$ per quasi-Newton iteration. This approach we refer to as JADOC: Joint Approximate Diagonalization under Orthogonality Constraints. We compare our algorithm to two important existing methods and show JADOC has superior runtime while yielding a highly similar degree of diagonalization. The JADOC algorithm is implemented as open-source Python code, available at https://github.com/devlaming/jadoc.
翻译:一组正(semi)- 确定基质的联合二进制(semi)- 确定基质的联合二进制具有广泛的分析应用,例如对共同主要组成部分的估计,对多种差异组成部分的估计,以及盲点信号分离等。然而,当所涉矩阵的分解元体不同时,联合二进制是一个计算上具有挑战性的问题。据我们所知,目前采用的方法要求每转录至少需要O(KN3)3美元的时间,如果考虑的是不同的美元N美元时,则每转录至少需要1美元(KN3)3美元。我们通过应用正方位限制和维度减少技术重新配置这一优化问题。在这样做时,我们将联合对二进制的计算复杂性降低到每准- Newton 迭代号$($3) 。我们称之为 JADOC: 在Othoconomicaldalizalizionalizionalizalization) 。我们将我们的算法比对两种重要的现有方法,并显示JADOC的运行时间优,同时产生高度相似的二进调化程度。JADADODOC算算算法,在 MADDGUDUDI/ ALms/ apps.