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.

An Application of Computer Algebra and Dynamical Systems.

Symbolic matrix factorization for differential-algebraic equations index reduction

A Joint Diagonalisation Approach for Linear Stochastic Systems

Fixed-Point Iteration Schemes to Solve Symmetric Algebraic Riccati Equation $$XBX-XA-A^{T}X-C=0$$

Computing Symmetric Positive Definite Solutions of Three Types of Nonlinear Matrix Equations

Neural Network Approach to Computing Outer Inverses Based on the Full Rank Representation

The Bit Complexity of Dynamic Algebraic Formulas and their Determinants

Algorithms for Robust Inversion of Dynamical Systems

On the Closed Form Expression of Elementary Symmetric Polynomials and the Inverse of Vandermonde Matrix

A rational Krylov methods for large scale linear multidimensional dynamical systems

Jacobi Stability Analysis for Systems of ODEs Using Symbolic Computation

Algebraic Analysis on Asymptotic Stability of Continuous Dynamical Systems

Computer Algebra, Power Series and Summation

Symbolic dynamics and the stable algebra of matrices

On the Solution of Graph Isomorphism by Dynamical Algorithms

Synchronization of dynamical systems: an approach using a Computer Algebra System

Computing the exponential of matrices

A Semi-Algebraic Approach for the Computation of Lyapunov Functions

Symbolic Computation on a (2+1)-Dimensional Generalized Nonlinear Evolution System in Fluid Dynamics, Plasma Physics, Nonlinear Optics and Quantum Mechanics

Structural Preprocessing Method for Nonlinear Differential-Algebraic Equations Using Linear Symbolic Matrices

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