Systems of coupled generalized Sylvester matrix equations

This paper studies some systems of coupled generalized Sylvester matrix equations. We present some necessary and sufficient conditions for the solvability to these systems. We give the expressions of the general solutions to the systems when their solvability conditions are satisfied.     

A note on combined generalized Sylvester matrix equations

The solution of two combined generalized Sylvester matrix equations is studied. It is first shown that the two combined generalized Sylvester matrix equations can be converted into a normal Sylvester matrix equation through extension, and then with the help of a result for solution to normal Sylvester matrix equations, the complete solution to the two combined generalized Sylvester matrix equations is derived. A demonstrative example shows the effect of the proposed approach.     

Efficient iterative solutions to general coupled matrix equations

Linear matrix equations are encountered in many systems and control applications.In this paper,we consider the general coupled matrix equations（including the generalized coupled Sylvester matrix equations as a special case）l t=1EstYtFst = Gs,s = 1,2,···,l over the generalized reflexive matrix group（Y1,Y2,···,Yl）.We derive an efcient gradient-iterative（GI） algorithm for fnding the generalized reflexive solution group of the general coupled matrix equations.Convergence analysis indicates that the algorithm always converges to the generalized reflexive solution group for any initial generalized reflexive matrix group（Y1（1）,Y2（1）,···,Yl（1））.Finally,numerical results are presented to test and illustrate the performance of the algorithm in terms of convergence,accuracy as well as the efciency.     

On global Hessenberg based methods for solving Sylvester matrix equations

In the first part of this paper, we investigate the use of Hessenberg-based methods for solving the Sylvester matrix equation $AX+XB=C$. To achieve this goal, the Sylvester form of the global generalized Hessenberg process is presented. Using this process, different methods based on a Petrov–Galerkin or on a minimal norm condition are derived. In the second part, we focus on the SGl-CMRH method which is based on the Sylvester form of the Hessenberg process with pivoting strategy combined with a minimal norm condition. In order to accelerate the SGl-CMRH method, a preconditioned framework of this method is also considered. It includes both fixed and flexible variants of the SGl-CMRH method. Moreover, the connection between the flexible preconditioned SGl-CMRH method and the fixed one is studied and some upper bounds for the residual norm are obtained. In particular, application of the obtained theoretical results is investigated for the special case of solving linear systems of equations with several right-hand sides. Finally, some numerical experiments are given in order to evaluate the effectiveness of the proposed methods.     

On the generalized Sylvester mapping and matrix equations

General parametric solution to a family of generalized Sylvester matrix equations arising in linear system theory is presented by using the so-called generalized Sylvester mapping which has some elegant properties. The solution consists of some polynomial matrices satisfying certain conditions and a parametric matrix representing the degree of freedom in the solution. The results provide great convenience to the computation and analysis of the solutions to this family of equations, and can perform important functions in many analysis and design problems in linear system theory. It is also expected that this so-called generalized Sylvester mapping tool may have some other applications in control system theory.     

Solutions to the generalized Sylvester matrix equations by a singular value decomposition

In this paper,solutions to the generalized Sylvester matrix equations AX-XF=BY and MXN-X=TY with A,M∈Rn×n,B,T∈Rn×n,F,N∈Rp×p and the matrices N,F being in companion form,are established by a singular value decomposition of a matrix with dimensions n×(n pr).The algorithm proposed in this paper for the euqation AX-XF=BY does not require the controllability of matrix pair(A,B)andthe restriction that A,F do not have common eigenvalues.Since singular value decomposition is adopted,the algorithm is numerically stable and may provide great convenience to the computation of the solution to these equations,and can perform important functions in many design problems in control systems theory.     

Gradient-based maximal convergence rate iterative method for solving linear matrix equations

This paper is concerned with numerical solutions to general linear matrix equations including the well-known Lyapunov matrix equation and Sylvester matrix equation as special cases. Gradient based iterative algorithm is proposed to approximate the exact solution. A necessary and sufficient condition guaranteeing the convergence of the algorithm is presented. A sufficient condition that is easy to compute is also given. The optimal convergence factor such that the convergence rate of the algorithm is maximized is established. The proposed approach not only gives a complete understanding on gradient based iterative algorithm for solving linear matrix equations, but can also be served as a bridge between linear system theory and numerical computing. Numerical example shows the effectiveness of the proposed approach.     

Solutions to the generalized Sylvester matrix equations by a singular value decomposition

In this paper, solutions to the generalized Sylvester matrix equations AX-XF=BY and MXN-X=TY with A,M ∈R ,B,T∈R, F,N∈R and the matrices N, F being in companion form, are established by a singular value decomposition of a matrix with dimensions n×(n + pr). The algorithm proposed in this paper for the euqation AX-XF = BY does not require the controllability of matrix pair (A,B) and the restriction that A,F don’t have common eigenvalues. Since singular value decomposition is adopted, the algorithm is numerically stable and may provide great convenience to the computation of the solution to these equations, and can perform important functions in many design problems in control systems theory.     

A new solution to the generalized Sylvester matrix equation

This note deals with the problem of solving the generalized Sylvester matrix equation AV-EVF=BW, with F being an arbitrary matrix, and provides complete general parametric expressions for the matrices V and W satisfying this equation. The primary feature of this solution is that the matrix F does not need to be in any canonical form, and may be even unknown a priori. The results provide great convenience to the computation and analysis of the solutions to this class of equations, and can perform important functions in many analysis and design problems in control systems theory.     

A symmetric preserving iterative method for generalized Sylvester equation

The generalized Sylvester matrix equation AX + YB = C is encountered in many systems and control applications, and also has several applications relating to the problem of image restoration, and the numerical solution of implicit ordinary differential equations. In this paper, we construct a symmetric preserving iterative method, basing on the classic Conjugate Gradient Least Squares (CGLS) method, for AX + YB = C with the unknown matrices X, Y having symmetric structures. With this method, for any arbitrary initial symmetric matrix pair, a desired solution can be obtained within finitely iterate steps. The unique optimal (least norm) solution can also be obtained by choosing a special kind of initial matrix. We also consider the matrix nearness problem. Some numerical results confirm the efficiency of these algorithms. It is more important that some numerical stability analysis on the matrix nearness problem is given combined with numerical examples, which is not given in the earlier papers. Copyright © 2010 John Wiley and Sons Asia Pte Ltd and Chinese Automatic Control Society     

An iterative algorithm for solving a pair of matrix equations over generalized centro-symmetric matrices

A matrix is said to be a symmetric orthogonal matrix if . A matrix is said to be generalized centro-symmetric (generalized central anti-symmetric) with respect to P, if A=PAP (A=−PAP). The generalized centro-symmetric matrices have wide applications in information theory, linear estimate theory and numerical analysis. In this paper, we propose a new iterative algorithm to compute a generalized centro-symmetric solution of the linear matrix equations . We show, when the matrix equations are consistent over generalized centro-symmetric matrix Y, for any initial generalized centro-symmetric matrix Y₁, the sequence {Y_k} generated by the introduced algorithm converges to a generalized centro-symmetric solution of matrix equations . The least Frobenius norm generalized centro-symmetric solution can be derived when a special initial generalized centro-symmetric matrix is chosen. Furthermore, the optimal approximation generalized centro-symmetric solution to a given generalized centro-symmetric matrix can be derived. Several numerical examples are given to show the efficiency of the presented method.     

An iterative algorithm for solving a class of matrix equations

In this paper, an iterative algorithm is presented to solve the Sylvester and Lyapunov matrix equations. By this iterative algorithm, for any initial matrix X1, a solution X＊ can be obtained within finite iteration steps in the absence of roundoff errors. Some examples illustrate that this algorithm is very efficient and better than that of [ 1 ] and [2].     

Gradient-based iterative algorithms for generalized coupled Sylvester-conjugate matrix equations

By applying the hierarchical identification principle, the gradient-based iterative algorithm is suggested to solve a class of complex matrix equations. With the real representation of a complex matrix as a tool, the sufficient and necessary conditions for the convergence factor are determined to guarantee that the iterative solutions given by the proposed algorithm converge to the exact solution for any initial matrices. Also, we solve the problem which is proposed by Wu et al. (2010). Finally, some numerical examples are provided to illustrate the effectiveness of the proposed algorithms and testify the conclusions suggested in this paper.     

预条件的平方Smith法求解大型Sylvester矩阵方程

提出了一种预条件的平方Smith算法求解大型连续Sylvester矩阵方程,该算法利用交替方向隐式迭代(ADI)来构造预条件算子,将原方程转换为非对称Stein方程,并在Krylov子空间中应用平方Smith法迭代产生低秩逼近解。数值实验表明,与已知的Jacobi迭代法等算法相比,该算法有更好的迭代效率和收敛精度。     

Finite iterative algorithm for solving coupled Lyapunov equations appearing in continuous-time Markov jump linear systems

An iterative algorithm for solving coupled algebraic Lyapunov equations appearing in continuous-time linear systems with Markovian transitions is established. The algorithm is computationally efficient since it can obtain the solution within finite steps in the absence of round-off errors.     

Speeding up solving of differential matrix Riccati equations using GPGPU computing and MATLAB

In this work, we developed a parallel algorithm to speed up the resolution of differential matrix Riccati equations using a backward differentiation formula algorithm based on a fixed‐point method. The role and use of differential matrix Riccati equations is especially important in several applications such as optimal control, filtering, and estimation. In some cases, the problem could be large, and it is interesting to speed it up as much as possible. Recently, modern graphic processing units (GPUs) have been used as a way to improve performance. In this paper, we used an approach based on general‐purpose computing on graphics processing units. We used NVIDIA © GPUs with unified architecture. To do this, a special version of basic linear algebra subprograms for GPUs, called CUBLAS, and a package (three different packages were studied) to solve linear systems using GPUs have been used. Moreover, we developed a MATLAB © toolkit to use our implementation from MATLAB in such a way that if the user has a graphic card, the performance of the implementation is improved. If the user does not have such a card, the algorithm can also be run using the machine CPU. Experimental results on a NVIDIA Quadro FX 5800 are shown. Copyright © 2011 John Wiley & Sons, Ltd.