An efficient algorithm for solving general coupled matrix equations and its application

The general coupled matrix equations (I){A"1"1X"1B"1"1+A"1"2X"2B"1"2+...+A"1"lX"lB"1"l=C"1,A"2"1X"1B"2"1+A"2"2X"2B"2"2+...+A"2"lX"lB"2"l=C"2,@?A"l"1X"1B"l"1+A"l"2X"2B"l"2+...+A"l"lX"lB"l"l=C"l, (including the generalized coupled Sylvester matrix equations as special cases) have nice applications in various branches of control and system theory. In this paper, by extending the idea of conjugate gradient method, we propose an efficient iterative algorithm to solve the general coupled matrix equations (I). When the matrix equations (I) are consistent, for any initial matrix group, a solution group can be obtained within finite iteration steps in the absence of roundoff errors. The least Frobenius norm solution group of the general coupled matrix equations can be derived when a suitable initial matrix group is chosen. We can use the proposed algorithm to find the optimal approximation solution group to a given matrix group ([email protected]^"1,[email protected]^"2,...,[email protected]^"l) in a Frobenius norm within the solution group set of the matrix equations (I). Also several numerical examples are given to illustrate that the algorithm is effective. Furthermore, the application of the proposed algorithm for solving the system of matrix equations {D"1XE"1=F"1,@?D"pXE"p=F"p, over (R,S)-symmetric and (R,S)-skew symmetric matrices is highlighted.

[1]  Khalide Jbilou,et al.  Block Krylov Subspace Methods for Solving Large Sylvester Equations , 2002, Numerical Algorithms.

[2]  Keith G. Woodgate.  Least-squares solution of F = PG over positive semidefinite symmetric P , 1996 .

[3]  Feng Ding,et al.  Multirate crosstalk identification in xDSL systems , 2006, IEEE Transactions on Communications.

[4]  Qing-Wen Wang,et al.  Ranks and the least-norm of the general solution to a system of quaternion matrix equations , 2009 .

[5]  M. Dehghan,et al.  Efficient iterative method for solving the second-order sylvester matrix equation EVF 2 -AVF-CV=BW , 2009 .

[6]  William F. Trench,et al.  Minimization problems for (R,S)-symmetric and (R,S)-skew symmetric matrices , 2004 .

[7]  Michael A. Epton.  Methods for the solution ofAXD−BXC=E and its application in the numerical solution of implicit ordinary differential equations , 1980 .

[8]  Guang-Ren Duan,et al.  Solutions to a family of matrix equations by using the Kronecker matrix polynomials , 2009, Appl. Math. Comput..

[9]  SI Shu-hong.  Least-squares solutions of inverse problems for anti-symmetric orthogonal anti-symmetric matrices , 2007 .

[10]  Yiqin Lin,et al.  Minimal residual methods augmented with eigenvectors for solving Sylvester equations and generalized Sylvester equations , 2006, Appl. Math. Comput..

[11]  Fikret A. Aliev,et al.  Optimization of linear control systems : analytical methods and computational algorithms , 1998 .

[12]  Tongwen Chen,et al.  Hierarchical least squares identification methods for multivariable systems , 2005, IEEE Transactions on Automatic Control.

[13]  Mehdi Dehghan,et al.  On the reflexive and anti-reflexive solutions of the generalised coupled Sylvester matrix equations , 2010, Int. J. Syst. Sci..

[14]  Thilo Penzl,et al.  A Cyclic Low-Rank Smith Method for Large Sparse Lyapunov Equations , 1998, SIAM J. Sci. Comput..

[15]  Mehdi Dehghan,et al.  An iterative algorithm for the reflexive solutions of the generalized coupled Sylvester matrix equations and its optimal approximation , 2008, Appl. Math. Comput..

[16]  B. Kågström,et al.  Generalized Schur methods with condition estimators for solving the generalized Sylvester equation , 1989 .

[17]  Qingwen Wang,et al.  On solutions to the quaternion matrix equation , 2008 .

[18]  P. Lancaster,et al.  The Algebraic Riccati Equation , 1995 .

[19]  Feng Ding,et al.  Gradient Based Iterative Algorithms for Solving a Class of Matrix Equations , 2005, IEEE Trans. Autom. Control..

[20]  Hsin-Chu Chen.  Generalized Reflexive Matrices: Special Properties and Applications , 1998, SIAM J. Matrix Anal. Appl..

[21]  Peter C. Müller,et al.  Stability of Linear Mechanical Systems With Holonomic Constraints , 1993 .

[22]  M. Dehghan,et al.  The general coupled matrix equations over generalized bisymmetric matrices , 2010 .

[23]  Yang Shi,et al.  Output Feedback Stabilization of Networked Control Systems With Random Delays Modeled by Markov Chains , 2009, IEEE Transactions on Automatic Control.

[24]  M. Sadkane,et al.  Use of near-breakdowns in the block Arnoldi method for solving large Sylvester equations , 2008 .

[25]  Yimin Wei,et al.  A new projection method for solving large Sylvester equations , 2007 .

[26]  Feng Ding,et al.  Iterative least-squares solutions of coupled Sylvester matrix equations , 2005, Syst. Control. Lett..

[27]  Daizhan Cheng,et al.  The polynomial solution to the Sylvester matrix equation , 2006, Appl. Math. Lett..

[28]  Kuan-Yue Wang,et al.  A Matrix Equation , 1991, Econometric Theory.

[29]  Qing-Wen Wang,et al.  Bisymmetric and centrosymmetric solutions to systems of real quaternion matrix equations , 2005 .

[30]  Guang-Jing Song,et al.  A new solvable condition for a pair of generalized Sylvester equations. , 2009 .

[31]  Qing-Wen Wang,et al.  THE REFLEXIVE RE-NONNEGATIVE DEFINITE SOLUTION TO A QUATERNION MATRIX EQUATION ∗ , 2008 .

[32]  J. D. Roberts,et al.  Linear model reduction and solution of the algebraic Riccati equation by use of the sign function , 1980 .

[33]  Mehdi Dehghan,et al.  ON THE REFLEXIVE SOLUTIONS OF THE MATRIX EQUATION AXB + CYD = E , 2009 .

[34]  Huazhen Fang,et al.  Kalman filter‐based adaptive control for networked systems with unknown parameters and randomly missing outputs , 2009 .

[35]  Tongwen Chen,et al.  2-Norm Based Recursive Design of Transmultiplexers with Designable Filter Length , 2006 .

[36]  Guang-Ren Duan,et al.  Weighted least squares solutions to general coupled Sylvester matrix equations , 2009 .

[37]  Danny C. Sorensen,et al.  A Modified Low-Rank Smith Method for Large-Scale Lyapunov Equations , 2004, Numerical Algorithms.

[38]  Feng Ding,et al.  Hierarchical gradient-based identification of multivariable discrete-time systems , 2005, Autom..

[39]  G. Duan,et al.  An explicit solution to the matrix equation AX − XF = BY , 2005 .

[40]  Leiba Rodman,et al.  Algebraic Riccati equations , 1995 .

[41]  Nicholas J. Higham,et al.  Perturbation theory and backward error forAX−XB=C , 1993 .

[42]  J. Reid.  Large Sparse Sets of Linear Equations , 1973 .

[43]  Lothar Reichel,et al.  Application of ADI Iterative Methods to the Restoration of Noisy Images , 1996, SIAM J. Matrix Anal. Appl..

[44]  Mehdi Dehghan,et al.  An iterative method for solving the generalized coupled Sylvester matrix equations over generalized bisymmetric matrices , 2010 .

[45]  Xi-Yan Hu,et al.  An iteration method for the symmetric solutions and the optimal approximation solution of the matrix equation AXB=C , 2005, Appl. Math. Comput..

[46]  Lei Zhang,et al.  An efficient algorithm for the least-squares reflexive solution of the matrix equation A1XB1=C1, A2XB2=C2 , 2006, Appl. Math. Comput..

[47]  E. Wachspress.  Iterative solution of the Lyapunov matrix equation , 1988 .

[48]  R. A. Smith.  Matrix Equation $XA + BX = C$ , 1968 .

[49]  Guang-Ren Duan,et al.  Gradient based iterative algorithm for solving coupled matrix equations , 2009, Syst. Control. Lett..

[50]  Masoud Hajarian,et al.  A lower bound for the product of eigenvalues of solutions to matrix equations , 2009, Appl. Math. Lett..

[51]  Zhi-Bin Yan,et al.  Solutions to right coprime factorizations and generalized Sylvester matrix equations , 2008 .

[52]  Guang-Ren Duan,et al.  Solutions to generalized Sylvester matrix equation by Schur decomposition , 2007, Int. J. Syst. Sci..

[53]  W. Niethammer,et al.  SOR for AX−XB=C , 1991 .

[54]  William F. Trench.  Characterization and properties of .R;S /-commutative matrices , 2012 .

[55]  J. Woude,et al.  A system of real quaternion matrix equations with applications , 2009 .

[56]  A. Laub,et al.  Computation of system balancing transformations and other applications of simultaneous diagonalization algorithms , 1987 .

[57]  Feng Ding,et al.  On Iterative Solutions of General Coupled Matrix Equations , 2006, SIAM J. Control. Optim..

[58]  M. Dehghan,et al.  The Reflexive and Anti-Reflexive Solutions of a Linear Matrix Equation and Systems of Matrix Equations , 2010 .

[59]  Qingwen Wang.  The general solution to a system of real quaternion matrix equations , 2005 .

[60]  Zhen-yun Peng,et al.  An efficient iterative method for solving the matrix equation AXB + CYD = E , 2006, Numer. Linear Algebra Appl..

[61]  Guang-Ren Duan,et al.  On the generalized Sylvester mapping and matrix equations , 2008, Syst. Control. Lett..

[62]  Bo Kågström,et al.  LAPACK-style algorithms and software for solving the generalized Sylvester equation and estimating the separation between regular matrix pairs , 1994, TOMS.

[63]  Mehdi Dehghan,et al.  An iterative algorithm for solving a pair of matrix equations AYB=E, CYD=F over generalized centro-symmetric matrices , 2008, Comput. Math. Appl..

[64]  Feng Ding,et al.  Iterative solutions of the generalized Sylvester matrix equations by using the hierarchical identification principle , 2008, Appl. Math. Comput..

[65]  Wayne H. Enright,et al.  Improving the Efficiency of Matrix Operations in the Numerical Solution of Stiff Ordinary Differential Equations , 1978, TOMS.

[66]  Peter Lancaster,et al.  The theory of matrices , 1969 .

[67]  William F. Trench,et al.  Characterization and Properties of (R, S)-Symmetric, (R, S)-Skew Symmetric, and (R, S)-Conjugate Matrices , 2005, SIAM J. Matrix Anal. Appl..

[68]  Bin Zhou,et al.  A new solution to the generalized Sylvester matrix equation AV-EVF=BW , 2006, Syst. Control. Lett..

[69]  M. R. Osborne,et al.  A Riccati transformation method for solving linear BVPs. I: theoretical aspects , 1988 .

[70]  Mehdi Dehghan,et al.  Finite iterative algorithms for the reflexive and anti-reflexive solutions of the matrix equation A1X1B1+A2X2B2=C , 2009, Math. Comput. Model..

[71]  Qing-Wen Wang,et al.  The common solution to six quaternion matrix equations with applications , 2008, Appl. Math. Comput..