scholarly journals Biconjugate residual algorithm for solving general Sylvester-transpose matrix equations

Filomat ◽  
2018 ◽  
Vol 32 (15) ◽  
pp. 5307-5318 ◽  
Author(s):  
Masoud Hajarian

The present paper is concerned with the solution of the coupled generalized Sylvester-transpose matrix equations {A1XB1 + C1XD1 + E1XTF1 = M1, A2XB2 + C2XD2 + E2XTF2 = M2, including the well-known Lyapunov and Sylvester matrix equations. Based on a variant of biconjugate residual (BCR) algorithm, we construct and analyze an efficient algorithm to find the (least Frobenius norm) solution of the general Sylvester-transpose matrix equations within a finite number of iterations in the absence of round-off errors. Two numerical examples are given to examine the performance of the constructed algorithm.

2016 ◽  
Vol 40 (2) ◽  
pp. 647-656 ◽  
Author(s):  
Masoud Hajarian

The periodic discrete-time matrix equations have wide applications in stability theory, control theory and perturbation analysis. In this work, the biconjugate residual algorithm is generalized to construct a matrix iterative method to solve the periodic discrete-time generalized coupled Sylvester matrix equations [Formula: see text] The constructed method is shown to be convergent in a finite number of iterations in the absence of round-off errors. By comparing with other similar methods in practical computation, we give numerical results to demonstrate the accuracy and the numerical superiority of the constructed method.


2019 ◽  
Vol 42 (3) ◽  
pp. 503-517
Author(s):  
Masoud Hajarian

The study of linear matrix equations is extremely important in many scientific fields such as control systems and stability analysis. In this work, we aim to design the Hestenes-Stiefel (HS) version of biconjugate residual (Bi-CR) algorithm for computing the (least Frobenius norm) partially doubly symmetric solution [Formula: see text] of the general Sylvester matrix equations [Formula: see text] for [Formula: see text]. We show that the proposed algorithm converges in a finite number of iterations. Finally, numerical results compare the proposed algorithm to alternative algorithms.


2017 ◽  
Vol 34 (5) ◽  
pp. 1446-1467 ◽  
Author(s):  
Masoud Hajarian

Purpose The purpose of this paper is to obtain an iterative algorithm to find the solution of the coupled Sylvester-like matrix equations. Design/methodology/approach In this work, the matrix form of the conjugate direction (CD) algorithm to find the solution X of the coupled Sylvester-like matrix equations: {A1XB1+M1f1(X)N1=F1,A2XB2+M2f2(X)N2=F2,with fi(X) = X, fi(X) = X¯, fi(X) = XT and fi(X) = XH for i = 1; 2 has been established. Findings It is proven that the algorithm converges to the solution within a finite number of iterations in the absence of round-off errors. Finally, four numerical examples were used to test the proficiency and convergence of the established algorithm. Originality/value The numerical examples have led the author to believe that the generalized CD (GCD) algorithm is efficient and it converges more rapidly in comparison with the CGNR and CGNE algorithms.


2012 ◽  
Vol 2012 ◽  
pp. 1-6
Author(s):  
Xuefeng Duan ◽  
Chunmei Li

Based on the alternating projection algorithm, which was proposed by Von Neumann to treat the problem of finding the projection of a given point onto the intersection of two closed subspaces, we propose a new iterative algorithm to solve the matrix nearness problem associated with the matrix equations AXB=E, CXD=F, which arises frequently in experimental design. If we choose the initial iterative matrix X0=0, the least Frobenius norm solution of these matrix equations is obtained. Numerical examples show that the new algorithm is feasible and effective.


2013 ◽  
Vol 2013 ◽  
pp. 1-15
Author(s):  
Zhongli Zhou ◽  
Guangxin Huang

The general coupled matrix equations (including the generalized coupled Sylvester matrix equations as special cases) have numerous applications in control and system theory. In this paper, an iterative algorithm is constructed to solve the general coupled matrix equations over reflexive matrix solution. When the general coupled matrix equations are consistent over reflexive matrices, the reflexive solution can be determined automatically by the iterative algorithm within finite iterative steps in the absence of round-off errors. The least Frobenius norm reflexive solution of the general coupled matrix equations can be derived when an appropriate initial matrix is chosen. Furthermore, the unique optimal approximation reflexive solution to a given matrix group in Frobenius norm can be derived by finding the least-norm reflexive solution of the corresponding general coupled matrix equations. A numerical example is given to illustrate the effectiveness of the proposed iterative algorithm.


2012 ◽  
Vol 2012 ◽  
pp. 1-28 ◽  
Author(s):  
Feng Yin ◽  
Guang-Xin Huang

An iterative algorithm is constructed to solve the generalized coupled Sylvester matrix equations(AXB-CYD,EXF-GYH)=(M,N), which includes Sylvester and Lyapunov matrix equations as special cases, over generalized reflexive matricesXandY. When the matrix equations are consistent, for any initial generalized reflexive matrix pair[X1,Y1], the generalized reflexive solutions can be obtained by the iterative algorithm within finite iterative steps in the absence of round-off errors, and the least Frobenius norm generalized reflexive solutions can be obtained by choosing a special kind of initial matrix pair. The unique optimal approximation generalized reflexive solution pair[X̂,Ŷ]to a given matrix pair[X0,Y0]in Frobenius norm can be derived by finding the least-norm generalized reflexive solution pair[X̃*,Ỹ*]of a new corresponding generalized coupled Sylvester matrix equation pair(AX̃B-CỸD,EX̃F-GỸH)=(M̃,Ñ), whereM̃=M-AX0B+CY0D,Ñ=N-EX0F+GY0H. Several numerical examples are given to show the effectiveness of the presented iterative algorithm.


Filomat ◽  
2016 ◽  
Vol 30 (9) ◽  
pp. 2503-2520 ◽  
Author(s):  
Masoud Hajarian

The periodic matrix equations are strongly related to analysis of periodic control systems for various engineering and mechanical problems. In this work, a matrix form of the conjugate gradient for least squares (MCGLS) method is constructed for obtaining the least squares solutions of the general discrete-time periodic matrix equations ?t,j=1 (Ai,jXi,jBi,j + Ci,jXi+1,jDi,j)=Mi, i=1,2,.... It is shown that the MCGLS method converges smoothly in a finite number of steps in the absence of round-off errors. Finally two numerical examples show that the MCGLS method is efficient.


2020 ◽  
Vol 2020 ◽  
pp. 1-5
Author(s):  
Ehab A. El-Sayed ◽  
Eid E. El Behady

This paper considers a new method to solve the first-order and second-order nonhomogeneous generalized Sylvester matrix equations AV+BW= EVF+R and MVF2+DV F+KV=BW+R, respectively, where A,E,M,D,K,B, and F are the arbitrary real known matrices and V and W are the matrices to be determined. An explicit solution for these equations is proposed, based on the orthogonal reduction of the matrix F to an upper Hessenberg form H. The technique is very simple and does not require the eigenvalues of matrix F to be known. The proposed method is illustrated by numerical examples.


2018 ◽  
Vol 41 (10) ◽  
pp. 2708-2725 ◽  
Author(s):  
Masoud Hajarian

As is well known, periodic matrix equations have wide applications in many areas of control and system theory. This paper is devoted to a study of the numerical solutions of a general type of periodic matrix equations. We present three types of biconjugate residual (BCR) method to find the generalized bisymmetric periodic solutions [Formula: see text] of general periodic matrix equations [Formula: see text] The main theorems of this paper show that the presented methods can compute the generalized bisymmetric periodic solutions in a finite number of steps in the absence of round-off errors. We give two numerical examples to illustrate and interpret the theoretical results.


2012 ◽  
Vol 2012 ◽  
pp. 1-18 ◽  
Author(s):  
Feng Yin ◽  
Guang-Xin Huang

The generalized coupled Sylvester systems play a fundamental role in wide applications in several areas, such as stability theory, control theory, perturbation analysis, and some other fields of pure and applied mathematics. The iterative method is an important way to solve the generalized coupled Sylvester systems. In this paper, an iterative algorithm is constructed to solve the minimum Frobenius norm residual problem: min over generalized reflexive matrix . For any initial generalized reflexive matrix , by the iterative algorithm, the generalized reflexive solution can be obtained within finite iterative steps in the absence of round-off errors, and the unique least-norm generalized reflexive solution can also be derived when an appropriate initial iterative matrix is chosen. Furthermore, the unique optimal approximate solution to a given matrix in Frobenius norm can be derived by finding the least-norm generalized reflexive solution of a new corresponding minimum Frobenius norm residual problem: with , . Finally, several numerical examples are given to illustrate that our iterative algorithm is effective.


Sign in / Sign up

Export Citation Format

Share Document