Abstract:PurposeThe purpose of this paper is to present a novel gradient‐based iterative algorithm for the joint diagonalization of a set of real symmetric matrices. The approximate joint diagonalization of a set of matrices is an important tool for solving stochastic linear equations. As an application, reliability analysis of structures by using the stochastic finite element analysis based on the joint diagonalization approach is also introduced in this paper, and it provides useful references to practical engineers.Design/methodology/approachBy starting with a least squares (LS) criterion, the authors obtain a classical nonlinear cost‐function and transfer the joint diagonalization problem into a least squares like minimization problem. A gradient method for minimizing such a cost function is derived and tested against other techniques in engineering applications.FindingsA novel approach is presented for joint diagonalization for a set of real symmetric matrices. The new algorithm works on the numerical gradient base, and solves the problem with iterations. Demonstrated by examples, the new algorithm shows the merits of simplicity, effectiveness, and computational efficiency.Originality/valueA novel algorithm for joint diagonalization of real symmetric matrices is presented in this paper. The new algorithm is based on the least squares criterion, and it iteratively searches for the optimal transformation matrix based on the gradient of the cost function, which can be computed in a closed form. Numerical examples show that the new algorithm is efficient and robust. The new algorithm is applied in conjunction with stochastic finite element methods, and very promising results are observed which match very well with the Monte Carlo method, but with higher computational efficiency. The new method is also tested in the context of structural reliability analysis. The reliability index obtained with the joint diagonalization approach is compared with the conventional Hasofer Lind algorithm, and again good agreement is achieved.

A shift‐splitting Jacobi‐gradient iterative algorithm for solving the matrix equation AV−V ̅B=C

A modified gradient‐based iterative algorithm for solving the complex conjugate and transpose matrix equations

Gradient-based iterative approach for solving constrained systems of linear matrix equations

Stationary splitting iterative methods for the matrix equation A X B = C

The relaxed gradient based iterative algorithm for solving the generalized coupled complex conjugate and transpose Sylvester matrix equations

A Modified Newton-based Matrix Splitting Iteration Method for Generalized Absolute Value Equations

Modified and accelerated relaxed gradient-based iterative algorithms for the complex conjugate and transpose matrix equations

The Shift-Invert Residual Arnoldi Method and the Jacobi-Davidson Method: Theory and Algorithms ∗

A general alternating direction implicit iteration method for solving complex symmetric linear systems

A general class of iterative splitting methods for solving linear systems

Numerical Solution of the Absolute Value Equations Using Two Matrix Splitting Fixed Point Iteration Methods

Inner Iterations in the Shift-Invert Residual Arnoldi Method and the Jacobi-Davidson Method

A Novel Joint Diagonalization Approach for Linear Stochastic Systems and Reliability Analysis

N A ] 2 2 N ov 2 01 1 The Shift-Invert Residual Arnoldi Method and the Jacobi – Davidson Method : Theory and Algorithms ∗

An Effective Iterative Solution for Independent Vector Analysis with Convergence Guarantees

Iterative algorithm to the generalized coupled Sylvester‐transpose matrix equations with application in robust and minimum norm observer design of linear systems

The Joint Bidiagonalization of a Matrix Pair with Inaccurate Inner Iterations

Low-rank generalized alternating direction implicit iteration method for solving matrix equations

Enhanced shifted Jacobi operational matrices of integrals: spectral algorithm for solving some types of ordinary and fractional differential equations

A sketch-and-project method for solving the matrix equation AXB = C