Next Article in Journal
Experimental and Numerical Simulation of a Symmetrical Three-Cylinder Buoy
Previous Article in Journal
A Method for Solving Approximate Partition Boundaries of Spatial Big Data Based on Histogram Bucket Sampling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solving a System of Sylvester-like Quaternion Matrix Equations

1
Department of Mathematics, Shanghai University, Shanghai 200444, China
2
Collaborative Innovation Center for the Marine Artificial Intelligence, Shanghai 200444, China
*
Author to whom correspondence should be addressed.
Submission received: 2 April 2022 / Revised: 4 May 2022 / Accepted: 17 May 2022 / Published: 20 May 2022
(This article belongs to the Section Mathematics)

Abstract

:
Using the ranks and Moore-Penrose inverses of involved matrices, in this paper we establish some necessary and sufficient solvability conditions for a system of Sylvester-type quaternion matrix equations, and give an expression of the general solution to the system when it is solvable. As an application of the system, we consider a special symmetry solution, named the η -Hermitian solution, for a system of quaternion matrix equations. Moreover, we present an algorithm and a numerical example to verify the main results of this paper.

1. Introduction

In 1952, Roth [1] studied the following one-sided generalized Sylvester matrix equation for the first time
A 1 X + Y B 1 = C 1 ,
which is widely used in system and control theory. Since then, many researches have paid attention to Sylvester-type matrix equations (e.g., [2,3,4,5]) because of their wide range of applications, such as in descriptor system control theory [6], neural networks [7], robust, feedback [8], graph theory [9] and other areas. For instance, Baksalary and Kala [10] established a necessary and sufficient condition for Equation (1) to have a solution and gave an expression of its general solution. In [11], Baksalary and Kala give a solvability condition for the equation
A X B + C Y D = E .
Wang investigated Equation (2) over arbitrary regular rings with identity [12].
In 1843, the very famous mathematician Hamilton discovered the quaternion. It is well known that quaternion algebra, denoted by H , is an associative and non-commutative division algebra over the real number field R , where
H = { a 0 + a 1 i + a 2 j + a 3 k | i 2 = j 2 = k 2 = ijk = 1 , a 0 , a 1 , a 2 , a 3 R } .
Since the 1970s, quaternions and the quaternion matrix have been studied a lot (e.g., [13,14,15,16]). The widespread applications of quaternions and the quaternion matrix include theoretical mechanics, optics, computer graphics, flight mechanics and aerospace technology, quantum physics, signal processing and so on (e.g., [17,18,19,20]).
In the last decade, the study of Sylvester-type matrix equations was extended to H (e.g., [21,22,23,24,25,26,27,28]). In 2012, Wang and He [29] presented the necessary and sufficient conditions for the Sylvester-type matrix equation
A 1 X 1 + X 2 B 1 + C 3 X 3 D 3 + C 4 X 4 D 4 = E 1
to be consistent and derived the expression of its general solution, which can be easily generalized to H . For the Sylvester-type matrix equations with multiple variables and multiple equations, Wang [4] gave a solvability condition and the general solution to the system of Sylvester-type matrix equations
A 3 W = B 3 , Z C 3 = D 3 , A 5 W + Z B 5 = D 4 .
Zhang [30] investigated the necessary and sufficient conditions for the solvability of the following system of Sylvester-like matrix equations
A 1 X = B 1 , X C 1 = D 1 , A 2 Y = B 2 , Y C 2 = D 2 , Z C 3 = D 3 , A 4 V = B 4 , A 6 V   + Z B 6 + A 7 X B 7 + A 8 Y B 8 = D 5 ,
and presented a formula of its general solution. We note that Equations (1)–(5) are the special cases of the following Sylvester-type quaternion matrix equations
A 1 X = B 1 , X C 1 = D 1 , A 2 Y = B 2 , Y C 2 = D 2 , A 3 W = B 3 , Z C 3 = D 3 , A 5 W +   Z B 5 = D 4 , A 4 V = B 4 , A 6 V   + Z B 6 + A 7 X B 7 + A 8 Y B 8 = D 5 ,
where A i , B i , C j , D k ( i = 1 , 8 ¯ , j = 1 , 3 ¯ , k = 1 , 5 ¯ ) are given matrices over H ; X , Y , Z , V , W are unknown.
Motivated by the work mentioned above, in this paper we aim to investigate the solvability conditions and the general solutions to a more general system of a Sylvester-type quaternion matrix equation, Equation (6). In 2011, Took et al. [31] defined a special class of symmetric matrices, named η -Hermitian. For η { i , j , k } , a quaternion matrix A is called η -Hermitian if A = A η * , where A η * = η A * η , A * is the conjugate and transpose matrix of A. It is well known that η -Hermitian matrices have some applications in linear modeling (e.g., [32,33,34]) and so on.
As an application of (6), we derive the solvability conditions and an expression of the η -Hermitian solution to the system of matrix equations
A 4 V = B 4 , A 1 X = B 1 , X = X η * , A 2 Y = B 2 , Y = Y η * , A 6 V + ( A 6 V ) η * + A 7 X A 7 η * + A 8 Y A 8 η * = D 5 , D 5 = D 5 η * ,
where A i ( i = 1 , 2 , 4 , 6 , 8 ¯ ) , B 1 , B 2 , B 4 , D 5 are given matrices over H ; X and Y are η -Hermitian matrices over H .
We organize the rest of this article as follows: In Section 2, we introduce the basic knowledge of quaternions and Moore-Penrose inverse of a quaternion matrix, and review some matrix equations. In Section 3, we establish the solvability conditions for the system of (6) in terms of the Moore-Penrose inverses and the ranks of the coefficients’ quaternion matrices in (6). In Section 4, we give an expression of the general solution to the system of (6), and illustrate the main results using a numerical example. In Section 5, we give some solvability conditions and an expression of the η -Hermitian solution to the system (7). Finally, we present a brief conclusion in Section 6 to end this paper.

2. Preliminaries

Let R and H m × n stand for the real number field and the set of all m × n matrix spaces over the quaternion algebra, respectively. The symbols r ( A ) , A * , I and 0 are denoted by the rank of a given quaternion matrix A, the conjugate transpose of A, an identity matrix, and a zero matrix with appropriate sizes, respectively. The Moore-Penrose inverse of A H l × k is defined to be the unique matrix, denoted by A , satisfying
A A A = A , A A A = A , ( A A ) * = A A , ( A A ) * = A A .
Moreover, L A = I A A and R A = I A A represent two projectors. Clearly, ( L A ) η * = R A η * and ( R A ) η * = L A η * of A.
The following lemma was given by Marsaglia and Stynan [35], which is also available over H .
Lemma 1
([35]). Let A H m × n , B H m × k , C H l × n , D H j × k and E H l × i . Then,
r A B L D R E C 0 = r A B 0 C 0 E 0 D 0 r ( D ) r ( E ) .
Lemma 2
([36]). Let A 1 and C 1 be known matrices with feasible dimensions over H . Then, the matrix equation A 1 X = C 1 has a solution if and only if R A 1 C 1 = 0 . In this case, its general solution is expressed as
X = A 1 C 1 + L A 1 T 1 ,
where T 1 is an arbitrary matrix of an appropriate size.
Lemma 3
([36]). Let B 1 and D 1 be known matrices with allowable dimensions over H . Then, the matrix equation Y B 1 = D 1 has a solution if and only if D 1 L B 1 = 0 . In this case, its general solution is
Y = D 1 B 1 + T 2 R B 1 ,
where T 2 is an arbitrary matrix of an appropriate size.
Lemma 4
([37]). Let A 1 , B 1 , C 1 and C 2 be the given matrices. Then, the system of matrix equations
A 1 Y = C 1 , Y B 1 = C 2
is consistent if and only if
R A 1 C 1 = 0 , C 2 L B 1 = 0 , A 1 C 2 = C 1 B 1 .
In this case, its general solution is
Y = A 1 C 1 + L A 1 C 2 B 1 + L A 1 T 3 R B 1 ,
where T 3 is an arbitrary matrix of an appropriate size.
Lemma 5
([10]). Let A , B and C be given over H . Then, the Equation (1) is solvable if and only if R A C L B = 0 . Under this condition, the general solution to Equation (1) can be expressed as
X = A C U 1 B + L A U 2 , Y = R A C B + A U 1 + U 3 R B ,
where U 1 , U 2 and U 3 are arbitrary matrices with appropriate sizes over H .
Lemma 6
([38]). Consider the following matrix equation over H
A 1 X 1 + X 2 B 1 + A 2 Y 1 B 2 + A 3 Y 2 B 3 + A 4 Y 3 B 4 = B ,
where A i , B i ( i = 1 , 4 ¯ ) , B are given and the others are unknown. Let
R A 1 A 2 = A 11 , R A 1 A 3 = A 22 , R A 1 A 4 = A 33 , B 2 L B 1 = B 11 , B 22 L B 11 = N 1 , B 3 L B 1 = B 22 , B 4 L B 1 = B 33 , R A 11 A 22 = M 1 , S 1 = A 22 L M 1 , R A 1 B L B 1 = T 1 , C = R M 1 R A 11 , C 1 = C A 33 , C 2 = R A 11 A 33 , C 3 = R A 22 A 33 , C 4 = A 33 , D = L B 11 L N 1 , D 1 = B 33 , D 2 = B 33 L B 22 , D 3 = B 33 L B 11 , D 4 = B 33 D , E 1 = C T 1 , E 2 = R A 11 T 1 L B 22 , E 3 = R A 22 T 1 L B 11 , E 4 = T 1 D , C 11 = ( L C 2 , L C 4 ) , D 11 = R D 1 R D 3 , C 22 = L C 1 , D 22 = R D 2 , C 33 = L C 3 , D 33 = R D 4 , E 11 = R C 11 C 22 , E 22 = R C 11 C 33 , E 33 = D 22 L D 11 , E 44 = D 33 L D 11 , M = R E 11 E 22 , N = E 44 L E 33 , F = F 2 F 1 , E = R C 11 F L D 11 , S = E 22 L M , F 11 = C 2 L C 1 , G 1 = E 2 C 2 C 1 E 1 D 1 D 2 , F 22 = C 4 L C 3 , G 2 = E 4 C 4 C 3 E 3 D 3 D 4 , F 1 = C 1 E 1 D 1 + L C 1 C 2 E 2 D 2 , F 2 = C 3 E 3 D 3 + L C 3 C 4 E 4 D 4 .
Then, the following statements are equivalent:
( 1 ) Equation (8) is consistent.
( 2 )
R C i E i = 0 , E i L D i = 0 ( i = 1 , 4 ¯ ) , R E 22 E L E 33 = 0 .
( 3 )
r B A 2 A 3 A 4 A 1 B 1 0 0 0 0 = r ( B 1 ) + r ( A 2 A 3 A 4 A 1 ) , r B A 2 A 4 A 1 B 3 0 0 0 B 1 0 0 0 = r ( A 2 A 4 A 1 ) + r B 3 B 1 , r B A 3 A 4 A 1 B 2 0 0 0 B 1 0 0 0 = r ( A 3 A 4 A 1 ) + r B 2 B 1 , r B A 4 A 1 B 2 0 0 B 3 0 0 B 1 0 0 = r B 2 B 3 B 1 + r ( A 4 A 1 ) , r B A 2 A 3 A 1 B 4 0 0 0 B 1 0 0 0 = r ( A 2 A 3 A 1 ) + r B 4 B 1 , r B A 2 A 1 B 3 0 0 B 4 0 0 B 1 0 0 = r B 3 B 4 B 1 + r ( A 2 A 1 ) , r B A 3 A 1 B 2 0 0 B 4 0 0 B 1 0 0 = r B 2 B 4 B 1 + r ( A 3 A 1 ) , r B A 1 B 2 0 B 3 0 B 4 0 B 1 0 = r B 2 B 3 B 4 B 1 + r ( A 1 ) , r B A 2 A 1 0 0 0 A 4 B 3 0 0 0 0 0 0 B 1 0 0 0 0 0 0 0 0 0 B A 3 A 1 A 4 0 0 0 B 2 0 0 0 0 0 0 B 1 0 0 0 B 4 0 0 B 4 0 0 0 = r B 3 0 B 1 0 0 B 2 0 B 1 B 4 B 4 + r A 2 A 1 0 0 A 4 0 0 A 3 A 1 A 4 .
In this case, the general solution to Equation (8) can be expressed as
X 1 = A 1 ( B A 2 Y 1 B 2 A 3 Y 2 B 3 A 4 Y 3 B 4 ) A 1 U 1 B 1 + L A 1 U 2 , X 2 = R A 1 ( B A 2 Y 1 B 2 A 3 Y 2 B 3 A 4 Y 3 B 4 ) B 1 + A 1 A 1 U 1 + U 3 R B 1 , Y 1 = A 11 T B 11 A 11 A 22 M 1 T B 11 A 11 S 1 A 22 T N 1 B 22 B 11 A 11 S 1 U 4 R N 1 B 22 B 11 + L A 11 U 5 + U 6 R B 11 , Y 2 = M 1 T B 22 + S 1 S 1 A 22 T N 1 + L M 1 L S 1 U 7 + U 8 R B 22 + L M 1 U 4 R N 1 , Y 3 = F 1 + L C 2 V 1 + V 2 R D 1 + L C 1 V 3 R D 2 , o r Y 3 = F 2 L C 4 W 1 W 2 R D 3 L C 3 W 3 R D 4 ,
where T = T 1 A 33 Y 3 B 33 ; U i ( i = 1 , 8 ¯ ) represents any matrix with appropriate dimensions over H ,
V 1 = ( I m 0 ) C 11 ( F C 22 V 3 D 22 C 33 W 3 D 33 ) C 11 U 11 D 11 + L C 11 U 12 , W 1 = ( 0 I m ) C 11 ( F C 22 V 3 D 22 C 33 W 3 D 33 ) C 11 U 11 D 11 + L C 11 U 12 , W 2 = R C 11 ( F C 22 V 3 D 22 C 33 W 3 D 33 ) D 11 + C 11 C 11 U 11 + U 21 R D 11 0 I n , V 2 = R C 11 ( F C 22 V 3 D 22 C 33 W 3 D 33 ) D 11 + C 11 C 11 U 11 + U 21 R D 11 I n 0 ,
V 3 = E 11 F E 33 E 11 E 22 M F E 33 E 11 S E 22 F N E 44 E 33 E 11 S U 31 R N E 44 E 33 + L E 11 U 32 + U 33 R E 33 , W 3 = M F E 44 + S S E 22 F N + L M L S U 41 + L M U 31 R N U 42 R E 44 ,
where U 11 , U 12 , U 21 , U 31 , U 32 , U 33 , U 41 and U 42 are any matrix with appropriate dimensions over H .

3. Solvability Conditions to the System (6)

The goal of this section is to give the necessary and sufficient conditions for the existence of a solution to system (6).
Theorem 1.
Let A i H m i × n i ( i = 1 , 4 ¯ ) , A 5 H m × n 3 , A 6 H m × n 4 , A 7 H m × n 1 , A 8 H m × n 2 , B j H m j × l j ( j = 1 , 2 ¯ ) , B 3 H m 3 × q , B 4 H m 4 × l , B 5 H l 3 × q , B 6 H l 3 × l , B 7 H l 1 × l , B 8 H l 2 × l , C k H l k × p k ( k = 1 , 3 ¯ ) , D j H n j × p j ( j = 1 , 2 ¯ ) , D 3 H m × l 2 , D 4 H m × q and D 5 H m × l . Set
A 11 = A 5 L A 3 , B 11 = R C 3 B 5 , C 11 = D 4 A 5 A 3 B 3 D 3 C 3 B 5 , A 22 = A 6 L A 4 ,
B 22 = R B 11 R C 3 B 6 , A 33 = A 7 L A 1 , B 33 = R C 1 B 7 , A 44 = A 8 L A 2 , B 44 = R C 2 B 8 ,
A 55 = A 11 , B 55 = R C 3 B 6 , M 1 = R A 22 A 33 , M 2 = R A 22 A 44 , M 3 = R A 22 A 55 ,
C 22 = D 5 A 6 A 4 B 4 D 3 C 3 B 6 R A 11 C 11 B 11 R C 3 B 6 A 7 ( A 1 B 1 + L A 1 D 1 C 1 ) B 7 A 8 ( A 2 B 2 + L A 2 D 2 C 2 ) B 8 ,
N 1 = B 33 L B 22 , N 2 = B 44 L B 22 , N 3 = B 55 L B 22 , G 1 = N 2 L N 1 , H 1 = R M 1 M 2 ,
S 1 = M 2 L H 1 , T = R A 22 C 22 L B 22 , P = R H 1 R M 1 , P 1 = P M 3 , P 2 = R M 1 M 3 ,
P 3 = R M 2 M 3 , P 4 = M 3 , Q = L N 1 L G 1 , Q 1 = N 3 , Q 2 = N 3 L N 2 , Q 3 = N 3 L N 1 ,
Q 4 = N 3 Q , E 1 = P T , E 2 = R M 1 T L N 2 , E 3 = R M 2 T L N 1 , E 4 = T Q ,
E 11 = ( L P 2 , L P 4 ) , F 11 = R Q 1 R Q 3 , E 22 = L P 1 , F 22 = R Q 2 , E 33 = L P 3 , F 33 = R Q 4 ,
M 11 = R E 11 E 22 , M 22 = R E 11 E 33 , M 33 = F 22 L F 11 , M 44 = F 33 L F 11 , M = R M 11 M 22 ,
N = M 44 L M 33 , F = F 2 F 1 , E = R E 11 F L F 11 , S = M 22 L M , G 11 = P 2 L P 1 ,
H 11 = E 2 P 2 P 1 E 1 Q 1 Q 2 , G 22 = P 4 L P 3 , H 22 = E 4 P 4 P 3 E 3 Q 3 Q 4 ,
F 2 = P 1 E 1 Q 1 + L P 1 P 2 E 2 Q 2 , F 1 = P 3 E 3 Q 3 + L P 3 P 4 E 4 Q 4 .
Then, the following statements are equivalent:
( 1 ) System (6) has a solution.
( 2 )
A 1 D 1 = B 1 C 1 , A 2 D 2 = B 2 C 2
and
R A 1 B 1 = 0 , D 1 L C 1 = 0 , R A 2 B 2 = 0 , D 2 L C 2 = 0 , R A 3 B 3 = 0 , D 3 L C 3 = 0 , R A 4 B 4 = 0 , R A 11 C 11 L B 11 = 0 , R P i E i = 0 , E i L Q i = 0 ( i = 1 , 4 ¯ ) , R M 22 E L M 33 = 0 .
( 3 ) (22) holds and
r A 1 B 1 = r A 1 , r C 1 D 1 = r C 1 , r A 2 B 2 = r A 2 , r C 2 D 2 = r C 2 ,
r A 3 B 3 = r A 3 , r C 3 D 3 = r C 3 , r A 4 B 4 = r A 4 ,
r D 4 A 5 D 3 B 5 0 C 3 B 3 A 3 0 = r A 5 A 3 + r ( B 5 C 3 ) ,
r D 5 A 7 A 8 A 6 D 4 A 5 D 3 B 6 0 0 0 B 5 0 C 3 B 1 B 7 A 1 0 0 0 0 0 B 2 B 8 0 A 2 0 0 0 0 B 4 0 0 A 4 0 0 0 0 0 0 0 B 3 A 3 0 = r A 7 A 8 A 5 A 6 A 1 0 0 0 0 A 2 0 0 0 0 0 A 4 0 0 A 3 0 + r B 6 B 5 C 3 ,
r D 5 A 7 A 6 A 8 D 2 D 4 A 5 D 3 B 8 0 0 C 2 0 0 0 B 6 0 0 0 B 5 0 C 3 B 1 B 7 A 1 0 0 0 0 0 B 4 0 A 4 0 0 0 0 0 0 0 0 B 3 A 3 0 = r A 7 A 6 A 5 A 1 0 0 0 A 4 0 0 0 A 3 + r B 8 C 2 0 0 B 6 0 B 5 C 3 ,
r D 5 A 8 A 6 A 7 D 1 D 4 A 5 D 3 B 7 0 0 C 1 0 0 0 B 6 0 0 0 B 5 0 C 3 B 2 B 8 A 2 0 0 0 0 0 B 4 0 A 4 0 0 0 0 0 0 0 0 B 3 A 3 0 = r A 8 A 6 A 5 A 2 0 0 0 A 4 0 0 0 A 3 + r B 7 C 1 0 0 B 6 0 B 5 C 3 ,
r D 5 A 6 A 7 D 1 A 8 D 2 D 4 A 5 D 3 B 7 0 C 1 0 0 0 0 B 8 0 0 C 2 0 0 0 B 6 0 0 0 B 5 0 C 3 B 4 A 4 0 0 0 0 0 0 0 0 0 B 3 A 3 0 = r B 7 C 1 0 0 0 B 8 0 C 2 0 0 B 6 0 0 B 5 C 3 + r A 6 A 5 A 4 0 0 A 3 ,
r D 5 A 7 A 8 A 6 D 3 B 6 0 0 0 C 3 B 1 B 7 A 1 0 0 0 B 2 B 8 0 A 2 0 0 B 4 0 0 A 4 0 = r A 7 A 8 A 6 A 1 0 0 0 A 2 0 0 0 A 4 + r B 6 C 3 ,
r D 5 A 7 A 6 A 8 D 2 D 3 B 8 0 0 C 2 0 B 6 0 0 0 C 3 B 1 B 7 A 1 0 0 0 B 4 0 A 4 0 0 = r B 8 C 2 0 B 6 0 C 3 + r A 7 A 6 A 1 0 0 A 4 ,
r D 5 A 8 A 6 A 7 D 1 D 3 B 7 0 0 C 1 0 B 6 0 0 0 C 3 B 2 B 8 A 2 0 0 0 B 4 0 A 4 0 0 = r B 7 C 1 0 B 6 0 C 3 + r A 8 A 6 A 2 0 0 A 4 ,
r D 5 A 6 A 7 D 1 A 8 D 2 D 3 B 7 0 C 1 0 0 B 8 0 0 C 2 0 B 6 0 0 0 C 3 B 4 A 4 0 0 0 = r B 7 C 1 0 0 B 8 0 C 2 0 B 6 0 0 C 3 + r A 6 A 4 ,
r D 5 A 7 A 6 0 0 0 A 8 D 2 D 4 0 0 0 A 5 D 3 0 0 B 8 0 0 0 0 0 C 2 0 0 0 0 0 0 0 0 B 6 0 0 0 0 0 0 B 5 0 0 0 0 C 3 0 0 0 0 0 D 5 A 8 A 6 0 0 A 7 D 1 D 4 0 0 0 A 5 D 3 0 0 0 B 7 0 0 0 0 C 1 0 0 0 0 0 0 0 0 0 B 6 0 0 0 0 0 B 5 0 0 0 0 C 3 B 6 0 0 B 6 0 0 0 0 0 0 C 3 0 0 0 0 B 1 B 7 A 1 0 0 0 0 0 0 0 0 0 0 0 0 0 B 4 0 A 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 B 2 B 8 A 2 0 0 0 0 0 0 0 0 0 0 0 0 0 B 4 0 A 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 B 3 0 0 0 A 3 0 0 0 0 0 0 0 0 A 4 0 0 0 B 3 0 0 0 A 3 0
= r B 8 0 C 2 0 0 0 0 0 0 B 6 0 0 B 5 0 0 0 C 3 0 0 B 7 0 0 C 1 0 0 0 0 0 B 6 0 0 0 B 5 0 0 C 3 B 6 B 6 0 0 0 0 C 3 0 0 + r A 7 A 6 0 0 A 5 0 0 0 A 8 A 6 0 A 5 A 1 0 0 0 0 0 0 A 4 0 0 0 0 0 0 A 2 0 0 0 0 0 0 A 4 0 0 0 0 0 0 A 3 0 0 0 0 0 0 A 3 .
Proof. 
( 1 ) ( 2 )
The proof is divided into three parts:
  • Firstly, we divide the system (6) into the following:
    A 3 W = B 3 , Z C 3 = D 3 , A 4 V = B 4 , A 1 X = B 1 , X C 1 = D 1 , A 2 Y = B 2 , Y C 2 = D 2 ,
    A 5 Z + W B 5 = D 4 ,
    A 6 V + Z B 6 + A 7 X B 7 + A 8 Y B 8 = D 5 ,
    and consider the solvability conditions and the general solution to the system of matrices of Equation (36). For more information, see Step 1.
  • Secondly, substituting the W and Z obtained in the first step into Equation (37) yields
    A 11 T 3 + T 4 B 11 = C 11 ,
    where A 11 , B 11 and C 11 are defined by (9); T 3 and T 4 are unknowns. For more information, see Step 2.
  • Finally, by substituting the X, Y, Z, and V obtained from the above two steps into Equation (38), we obtain a matrix equation with the following form
    A 22 T 5 + U 3 B 22 + A 33 T 1 B 33 + A 44 T 2 B 44 + A 55 U 1 B 55 = C 22 ,
    where A i i , B i i ( i = 2 , 5 ¯ ) and C 22 are given by (9)–(12); T 1 , T 2 , T 5 , U 1 and U 3 are unknowns. For more information, see Step 3.
We can obtain the results from the following steps: First, we consider the solvability conditions and the expression of the general solutions to the system of the matrix Equation (36).
Step 1. It follows from Lemmas 2–4 that system (36) has a solution if and only if (22) holds and
R A 1 B 1 = 0 , D 1 L C 1 = 0 , R A 2 B 2 = 0 , D 2 L C 2 = 0 , R A 3 B 3 = 0 , D 3 L C 3 = 0 , R A 4 B 4 = 0 .
In this case, the general solution to system (36) can be written as
X = A 1 B 1 + L A 1 D 1 C 1 + L A 1 T 1 R C 1 , Y = A 2 B 2 + L A 2 D 2 C 2 + L A 2 T 2 R C 2 , W = A 3 B 3 + L A 3 T 3 , Z = D 3 C 3 + T 4 R C 3 , V = A 4 B 4 + L A 4 T 5 ,
where T i ( i = 1 , 5 ¯ ) are arbitrary matrices over H with appropriate sizes.
Step 2. Substituting W , Z in (42) into (37) yields (39). According to Lemma 5, it follows that Equation (39) has a solution if and only if
R A 11 C 11 L B 11 = 0 .
In this case, the general solution to Equation (39) can be expressed as
T 3 = A 11 C 1 U 1 B 11 + L A 11 U 2 ,
T 4 = R A 11 C 11 B 11 + A 11 U 1 + U 3 R B 11 ,
where U 1 , U 2 and U 3 are any matrix with appropriate sizes over H .
Substituting (45) into Z = D 3 C 3 + T 4 R C 3 yields
Z = D 3 C 3 + R A 11 C 11 B 11 R C 3 + A 11 U 1 R C 3 + U 3 R B 11 R C 3 .
Step 3. By substituting X , Y , V in (42) and Z in (46) into (38), we obtain Equation (40). By using Lemma 6, Equation (40) is consistent if and only if
R P i E i = 0 , E i L Q i = 0 ( i = 1 , 4 ¯ ) , R M 22 E L M 33 = 0 ,
namely,
r C 22 A 33 A 44 A 55 A 22 B 22 0 0 0 0 = r ( B 22 ) + r ( A 33 A 44 A 55 A 22 ) ,
r C 22 A 33 A 55 A 22 B 44 0 0 0 B 22 0 0 0 = r ( A 33 A 55 A 22 ) + r B 44 B 22 ,
r C 22 A 44 A 55 A 22 B 33 0 0 0 B 22 0 0 0 = r ( A 44 A 55 A 22 ) + r B 33 B 22 ,
r C 22 A 55 A 22 B 33 0 0 B 44 0 0 B 22 0 0 = r B 33 B 44 B 22 + r ( A 55 A 22 ) ,
r C 22 A 33 A 44 A 22 B 55 0 0 0 B 22 0 0 0 = r ( A 33 A 44 A 22 ) + r B 55 B 22 ,
r C 22 A 33 A 22 B 44 0 0 B 55 0 0 B 22 0 0 = r B 44 B 55 B 22 + r ( A 33 A 22 ) ,
r C 22 A 44 A 22 B 33 0 0 B 55 0 0 B 22 0 0 = r B 33 B 55 B 22 + r ( A 44 A 22 ) ,
r C 22 A 22 B 33 0 B 44 0 B 55 0 B 22 0 = r B 33 B 44 B 55 B 22 + r ( A 22 ) ,
r C 22 A 33 A 22 0 0 0 A 55 B 44 0 0 0 0 0 0 B 22 0 0 0 0 0 0 0 0 0 C 22 A 44 A 22 A 55 0 0 0 B 33 0 0 0 0 0 0 B 22 0 0 0 B 55 0 0 B 55 0 0 0 = r B 44 0 B 22 0 0 B 33 0 B 22 B 55 B 55 + r A 33 A 22 0 0 A 55 0 0 A 44 A 22 A 55 .
In this case, the general solution to Equation (40) can be expressed as
T 5 = A 22 ( C 22 A 33 T 1 B 33 A 44 T 2 B 44 A 55 U 1 B 55 ) + A 22 V 1 B 22 + L A 22 V 2 , U 3 = R A 22 ( C 22 A 33 T 1 B 33 A 44 T 2 B 44 A 55 U 1 B 55 ) B 22 + A 22 A 22 V 1 + V 3 R B 22 , T 1 = M 1 T 11 N 1 M 1 M 2 H 1 T 11 N 1 M 1 S 1 M 2 T 11 G 1 N 2 N 1 M 1 S 1 V 4 R G 1 N 2 N 1 + L M 1 V 5 + V 6 R N 1 , T 2 = H 1 T 11 N 2 + S 1 S 1 M 2 T 11 G 1 + L H 1 L S 1 V 7 + V 8 R N 2 + L H 1 V 4 R G 1 , U 1 = F 1 + L P 2 W 1 + W 2 R Q 1 + L P 1 W 3 R Q 2 , o r U 1 = F 2 L P 2 W 4 W 5 R Q 3 L P 3 W 6 R Q 4 ,
where T 11 = T M 3 U 1 N 3 , V i ( i = 1 , 8 ¯ ) are any matrix with suitable dimensions over H ,
W 1 = I m 0 E 11 ( F E 22 W 3 F 22 E 33 W 6 F 33 ) E 11 U 11 F 11 + L E 11 U 12 , W 4 = 0 I m E 11 ( F E 22 W 3 F 22 E 33 W 6 F 33 ) E 11 U 11 F 11 + L E 11 U 12 , W 2 = R E 11 ( F E 22 W 3 F 22 E 33 W 6 F 33 ) F 11 + E 11 E 11 U 11 + U 21 R F 11 I n 0 , W 5 = R E 11 ( F E 22 W 3 F 22 E 33 W 6 F 33 ) F 11 + E 11 E 11 U 11 + U 21 R F 11 0 I n , W 3 = M 11 F M 33 M 11 M 22 M F M 33 M 11 S M 22 F N M 44 M 33 M 11 S U 31 R N M 44 M 33 + L M 11 U 32 + U 33 R M 33 , W 6 = M F M 44 + S S M 22 F N + L M L S U 41 + L M U 31 R N U 42 R M 44 ,
where U 11 , U 12 , U 21 , U 31 , U 32 , U 33 , U 41 and U 42 are any matrix with suitable dimensions over H .
To sum up, the system of matrices of Equation (6) has a solution if and only if (41), (43) and (47) hold.
( 2 ) ( 3 ) We divide it into three parts to prove its equivalence.
Part 1. In this part, we prove that (41) holds if and only if (24) and (25) hold. According to Lemma 1, it is easy to show that (41) holds if and only if (24) and (25) hold.
Part 2. In this part, we prove that (43) ⟺ (26). It follows from Lemma 1 and elementary operations that
( 43 ) r C 11 A 11 B 11 0 = r ( A 11 ) + r ( B 11 ) r C 11 A 5 L A 3 R C 3 B 5 0 = r ( A 5 L A 3 ) + r ( R C 3 B 5 ) r D 4 A 5 A 3 B 3 D 3 C 3 B 5 A 5 0 B 5 0 C 3 0 A 3 0 = r A 5 A 3 + r B 5 C 3 r D 4 A 5 D 3 B 5 0 C 3 B 3 A 3 0 = r A 5 A 3 + r B 5 C 3 ( 26 ) .
Part 3. In this part, we show that (47) holds if and only if (27) to (35) hold. By using Lemma 6, (47) holds if and only if (48) to (56) hold. Hence, we only show that (48) to (56) hold if and only if (27) to (35) hold, respectively. We first prove that (48) ⇔ (27).
Note that
X 0 = A 1 B 1 + L A 1 D 1 C 1 , Y 0 = A 2 B 2 + L A 2 D 2 C 2 , Z 0 = D 3 C 3 , V 0 = A 4 B 4 , W 0 = A 3 B 3
are the special solution to the equations
A 1 X = B 1 , X C 1 = D 1 , A 2 Y = B 2 , Y C 2 = D 2 , A 3 W = B 3 , Z C 3 = D 3 , A 4 V = B 4 ,
respectively. Then, we have that
C 11 = D 4 A 5 W 0 Z 0 B 5 ,
C 22 = D 5 A 6 V 0 Z 0 B 6 R A 11 C 11 B 11 R C 3 B 6 A 7 X 0 B 7 A 8 Y 0 B 8 .
It follows from Lemma 1 and elementary operations to (47) that
( 48 ) r C 22 A 7 A 8 A 11 A 6 0 R C 3 B 6 0 0 0 0 B 11 0 A 1 0 0 0 0 0 0 A 2 0 0 0 0 0 0 0 A 4 0 = r ( R C 3 B 6 B 11 ) + r A 7 A 8 A 11 A 6 A 1 0 0 0 0 A 2 0 0 0 0 0 A 4 r D 5 Z 0 B 6 A 7 A 8 A 6 C 11 A 11 R C 3 B 6 0 0 0 B 11 0 B 1 B 7 A 1 0 0 0 0 B 2 B 8 0 A 2 0 0 0 B 4 0 0 A 4 0 0 = r A 7 A 8 A 11 A 6 A 1 0 0 0 0 A 2 0 0 0 0 0 A 4 + r ( R C 3 B 6 B 11 )
r D 5 A 7 A 8 A 6 D 4 A 5 D 3 B 6 0 0 0 B 5 0 C 3 B 1 B 7 A 1 0 0 0 0 0 B 2 B 8 0 A 2 0 0 0 0 B 4 0 0 A 4 0 0 0 0 0 0 0 B 3 A 3 0 = r A 7 A 8 A 5 A 6 A 1 0 0 0 0 A 2 0 0 0 0 0 A 4 0 0 A 3 0 + r ( B 6 B 5 C 3 ) ( 27 ) .
Similarly, we can prove that R P 2 E 2 = 0 ( 28 ) , R P 3 E 3 = 0 ( 29 ) , R P 4 E 4 = 0 ( 30 ) and E i L Q i = 0 ( i = 1 , 4 ¯ ) hold if and only if (31) to (34) hold, respectively. Next, we show that R M 22 E L M 33 = 0 ( 35 ) . According to Lemma 1 and elementary operations, we have that
R M 22 E L M 33 = 0 r E D 22 D 33 0 = r ( D 22 ) + r ( D 33 ) r C 22 A 33 A 22 0 0 0 A 55 B 44 0 0 0 0 0 0 B 22 0 0 0 0 0 0 0 0 0 C 22 A 44 A 22 A 55 0 0 0 B 33 0 0 0 0 0 0 B 22 0 0 0 B 55 0 0 B 55 0 0 0 = r B 44 0 B 22 0 0 B 33 0 B 22 B 55 B 55 + r A 33 A 22 0 0 A 55 0 0 A 44 A 22 A 55 r C 22 A 7 A 6 0 0 0 A 11 0 0 0 0 0 B 8 0 0 0 0 0 0 C 2 0 0 0 0 R C 3 B 6 0 0 0 0 0 0 0 B 11 0 0 0 0 0 0 C 22 A 8 A 6 A 11 0 0 0 0 0 0 0 0 B 7 0 0 0 0 0 C 1 0 0 0 0 0 R C 3 B 6 0 0 0 0 0 0 B 11 0 B 6 0 0 B 6 0 0 0 0 0 0 0 C 3 0 A 1 0 0 0 0 0 0 0 0 0 0 0 0 A 4 0 0 0 0 0 0 0 0 0 0 0 0 0 A 2 0 0 0 0 0 0 0 0 0 0 0 0 A 4 0 0 0 0 0 0 = r B 8 0 C 2 0 0 0 0 R C 3 B 6 0 0 B 11 0 0 0 0 B 7 0 0 C 1 0 0 0 R C 3 B 6 0 0 0 B 11 0 B 6 B 6 0 0 0 0 C 3 + r A 7 A 6 0 0 A 11 0 0 A 8 A 6 A 11 A 1 0 0 0 0 0 A 4 0 0 0 0 0 A 2 0 0 0 0 0 A 4 0 r B 8 0 C 2 0 0 0 0 0 0 B 6 0 0 B 5 0 0 0 C 3 0 0 B 7 0 0 C 1 0 0 0 0 0 B 6 0 0 0 B 5 0 0 C 3 B 6 B 6 0 0 0 0 C 3 0 0 + r A 7 A 6 0 0 A 11 0 0 0 A 8 A 6 0 A 11 A 1 0 0 0 0 0 0 A 4 0 0 0 0 0 0 A 2 0 0 0 0 0 0 A 4 0 0
= r D 5 Z 0 B 6 A 7 A 6 0 0 0 A 8 D 2 C 11 0 0 0 A 11 0 0 0 B 8 0 0 0 0 0 C 2 0 0 0 0 0 0 0 0 B 6 0 0 0 0 0 0 B 5 0 0 0 0 C 3 0 0 0 0 0 Z 0 B 6 D 5 A 8 A 6 0 0 A 7 D 1 C 11 0 0 0 A 11 0 0 0 0 B 7 0 0 0 0 C 1 0 0 0 0 0 0 0 0 0 B 6 0 0 0 0 0 B 5 0 0 0 0 C 3 B 6 0 0 B 6 0 0 0 0 0 0 C 3 0 0 0 0 B 1 B 7 A 1 0 0 0 0 0 0 0 0 0 0 0 0 0 B 4 0 A 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 B 2 B 8 A 2 0 0 0 0 0 0 0 0 0 0 0 0 0 B 4 0 A 4 0 0 0 0 0 0 0 0 0 r B 8 0 C 2 0 0 0 0 0 0 B 6 0 0 B 5 0 0 0 C 3 0 0 B 7 0 0 C 1 0 0 0 0 0 B 6 0 0 0 B 5 0 0 C 3 B 6 B 6 0 0 0 0 C 3 0 0 + r A 7 A 6 0 0 A 5 0 0 0 A 8 A 6 0 A 5 A 1 0 0 0 0 0 0 A 4 0 0 0 0 0 0 A 2 0 0 0 0 0 0 A 4 0 0 0 0 0 0 A 3 0 0 0 0 0 0 A 3 = r D 5 A 7 A 6 0 0 0 A 8 D 2 D 4 0 0 0 A 5 D 3 0 0 B 8 0 0 0 0 0 C 2 0 0 0 0 0 0 0 0 B 6 0 0 0 0 0 0 B 5 0 0 0 0 C 3 0 0 0 0 0 D 5 A 8 A 6 0 0 A 7 D 1 D 4 0 0 0 A 5 D 3 0 0 0 B 7 0 0 0 0 C 1 0 0 0 0 0 0 0 0 0 B 6 0 0 0 0 0 B 5 0 0 0 0 C 3 B 6 0 0 B 6 0 0 0 0 0 0 C 3 0 0 0 0 B 1 B 7 A 1 0 0 0 0 0 0 0 0 0 0 0 0 0 B 4 0 A 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 B 2 B 8 A 2 0 0 0 0 0 0 0 0 0 0 0 0 0 B 4 0 A 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 B 3 0 0 0 A 3 0 0 0 0 0 0 0 0 0 0 0 0 B 3 0 0 0 A 3 0 ( 35 ) .

4. The General Solution to the System (6)

In this section, we give an expression for the general solution of Equation (6) by using the Moore-Penrose inverse. According to the proof of Theorem 1, we obtain the following theorem:
Theorem 2.
The general solution to system (6) can be expressed as follows when the solvability conditions are met:
X = A 1 B 1 + L A 1 D 1 C 1 + L A 1 T 1 R C 1 , Y = A 2 B 2 + L A 2 D 2 C 2 + L A 2 T 2 R C 2 , Z = D 3 C 3 + R A 11 C 11 B 11 R C 3 + A 11 U 1 R C 3 + U 3 R B 11 R C 3 , W = A 3 B 3 + L A 3 A 11 C 1 L A 3 A 11 U 1 B 11 + L A 3 L A 11 U 2 ,
V = A 4 B 4 + L A 4 A 22 ( C 22 A 33 T 1 B 33 A 44 T 2 B 44 A 55 U 1 B 55 ) + L A 4 A 22 V 1 B 22 + L A 4 L A 22 V 2 ,
where T 11 = T M 3 U 1 N 3 , V i ( i = 1 , 8 ¯ ) are arbitrary matrices with appropriate sizes.
T 1 = M 1 T 11 N 1 M 1 M 2 H 1 T 11 N 1 M 1 S 1 M 2 T 11 G 1 N 2 N 1 M 1 S 1 V 4 R G 1 N 2 N 1 + L M 1 V 5 + V 6 R N 1 , T 2 = H 1 T 11 N 2 + S 1 S 1 M 2 T 11 G 1 + L H 1 L S 1 V 7 + V 8 R N 2 + L H 1 V 4 R G 1 ,
U 3 = R A 22 ( C 22 A 33 T 1 B 33 A 44 T 2 B 44 A 55 U 1 B 55 ) B 22 + A 22 A 22 V 1 + V 3 R B 22 , U 1 = F 1 + L P 2 W 1 + W 2 R Q 1 + L P 1 W 3 R Q 2 , o r U 1 = F 2 L P 2 W 4 W 5 R Q 3 L P 3 W 6 R Q 4 , W 1 = I m 0 E 11 ( F E 22 W 3 F 22 E 33 W 6 F 33 ) E 11 U 11 F 11 + L E 11 U 12 , W 4 = 0 I m E 11 ( F E 22 W 3 F 22 E 33 W 6 F 33 ) E 11 U 11 F 11 + L E 11 U 12 , W 2 = R E 11 ( F E 22 W 3 F 22 E 33 W 6 F 33 ) F 11 + E 11 E 11 U 11 + U 21 R F 11 I n 0 , W 5 = R E 11 ( F E 22 W 3 F 22 E 33 W 6 F 33 ) F 11 + E 11 E 11 U 11 + U 21 R F 11 0 I n , W 3 = M 11 F M 33 M 11 M 22 M F M 33 M 11 S M 22 F N M 44 M 33 M 11 S U 31 R N M 44 M 33 + L M 11 U 32 + U 33 R M 33 , W 6 = M F M 44 + S S M 22 F N + L M L S U 41 + L M U 31 R N U 42 R M 44 ,
where U 11 , U 12 , U 21 , U 31 , U 32 , U 33 , U 41 and U 42 are arbitrary matrices over H of appropriate sizes.
Next, we discuss the special cases of the system of matrices of Equation (6). Letting A 3 , B 3 , A 5 , B 5 and D 4 vanish yields the following:
Corollary 1.
Suppose that A i , B i , C j , D j ( i = 1 , 4 ¯ , j = 1 , 5 ¯ ) and E 1 are given, denote
A 6 = A 4 L A 1 , B 6 = R B 1 B 4 , C 6 = C 4 L A 2 , D 6 = R B 2 D 4 , C 7 = C 5 L A 3 , D 7 = R B 3 D 5 , E 6 = E 1 A 4 A 1 C 1 D 1 B 1 B 4 C 4 A 2 C 2 + L A 2 D 2 B 2 D 4 C 5 A 3 C 3 + L A 3 D 3 B 3 D 5 , A = R A 6 C 6 , B = D 6 L B 6 , C = R A 6 C 7 , D = D 7 L B 6 ,
E = R A 6 E 6 L B 6 , M = R A C , N = D L B , S = C L M .
Then, the following statements are equivalent:
( 1 ) System (5) is consistent.
( 2 )
R A i C i = 0 , D i L B i = 0 ( i = 1 , 2 , 3 ) , A 2 D 2 = C 2 B 2 , A 3 D 3 = C 3 B 3 , R A E = M M E , E L B = E N N , R A E L D = 0 , R C E L B = 0 .
( 3 )
r A i C i = r ( A i ) , r B i D i = r ( B i ) ( i = 1 , 2 , 3 ) , A 2 D 2 = C 2 B 2 , A 3 D 3 = C 3 B 3 , r E 1 A 4 D 1 C 4 D 2 C 5 D 3 B 4 0 B 1 0 0 D 4 0 0 B 2 0 D 5 0 0 0 B 3 C 1 A 1 0 0 0 = r A 1 A 4 + r B 4 B 1 0 0 D 4 0 B 2 0 D 5 0 0 B 2 ,
r E 1 A 4 C 4 C 5 D 1 B 4 0 0 0 B 1 C 1 A 1 0 0 0 C 2 D 4 0 A 2 0 0 C 3 D 5 0 0 A 3 0 = r A 4 C 4 C 5 A 1 0 0 0 A 2 0 0 0 A 3 + r B 4 B 1 , r E 1 A 4 C 4 D 1 C 5 D 3 B 4 0 0 B 1 0 D 5 0 0 0 B 3 C 1 A 1 0 0 0 C 2 D 4 0 A 2 0 0 = r A 3 C 4 A 1 0 0 A 2 + r B 4 B 1 0 D 5 0 B 3 , r E 1 A 4 C 5 D 1 C 4 D 2 B 4 0 0 B 1 0 D 4 0 0 0 B 2 C 1 A 1 0 0 0 C 3 D 5 0 A 3 0 0 = r A 4 C 5 A 1 0 0 A 3 + r B 4 B 1 0 D 4 0 B 2 .
In this case, the general solution to system (5) can be expressed as
X 1 = A 1 C 1 + L A 1 U 1 , X 2 = D 1 B 1 + U 2 R B 1 , X 3 = A 2 C 2 + L A 2 D 2 B 2 + L A 2 U 3 R B 2 , X 4 = A 3 C 3 + L A 3 D 3 B 3 + L A 3 U 4 R B 3 , U 1 = A 6 E 6 C 6 U 3 D 6 C 7 U 4 D 7 A 6 W 2 B 6 + L A 6 W 1 , U 2 = R A 6 E 6 C 6 U 3 D 6 C 7 U 4 D 7 B 6 + A 6 A 6 W 2 + W 3 R B 6 , U 3 = A E B A C M E B A S C E N D B A S V 4 R N D B + L A V 1 + V 2 R B , U 4 = M E D + S S C E N + L M L S V 3 + L M V 4 R N + V 5 R D ,
where V i , W j ( i = 1 , 5 ¯ , j = 1 , 3 ¯ ) are arbitrary matrices over H with appropriate sizes.
Remark 1.
The above corollary is from the important findings of [30].
Letting A i , B i , C j , D j ( i = 1 , 2 , 4 , 6 , 7 , 8 , j = 1 , 2 ) and D 5 vanish, we have the following:
Corollary 2.
Given A 3 , B 3 , C 3 , D 3 , A 5 , B 5 and D 4 of feasible dimensions over H . Set A 11 = A 5 L A 3 , B 11 = R C 3 C and E 11 = D 4 A 5 A 3 B 3 D 3 C 3 B 5 . Then, the following statements are equivalent:
( 1 ) System (4) is consistent.
( 2 ) r ( A 3 B 3 ) = r A 3 , r C 3 D 3 = r C 3 , r D 4 A 5 D 3 B 5 0 C 3 B 3 A 3 0 = r A 5 A 3 + r ( B 5 C 3 ) .
In this case, the general solution to system (4) can be expressed as
W = A 3 B 3 + L A 3 A 11 E 11 A 3 W 2 B 11 + L A 11 W 1 , Z = D 3 C 3 + R A 11 E 11 B 11 + A 11 A 11 W 2 + W 3 R B 11 R C 3 ,
where W 1 , W 2 and W 3 are arbitrary matrices over H of appropriate sizes.
Remark 2.
The above corollary is from the vital investigation of [4].
Finally, we give Algorithm 1 and an example to illustrate the main results of this paper.
Algorithm 1: Algorithm for solving Equation (6)
(1) Feed the values of A i , B i , C j , D k ( i = 1 , 8 ¯ , j = 1 , 3 ¯ , k = 1 , 5 ¯ ) with conformable shapes over H .
(2) Compute the symbols in (9) to (21).
(3) Check (22), (23) or rank equalities in (24) to (35) hold or not. If no, then return “inconsistent”.
(4) Otherwise, compute X , Y , Z , V , W .
Example 1.
Consider the matrix of Equation (6). Assume
A 1 = 1 j i j i 1 1 i j , A 2 = i + j 1 + k i + k k 1 + i 1 + i + k , A 3 = j + k 1 + j i + k i i + j j + k , A 4 = 1 + j 2 + k 1 i + k 2 j + 2 k 1 + i + j + k 1 + i + j , A 5 = j j + k 1 j + k k 1 j i + k , A 6 = 1 + i + j + k 1 0 1 + i 1 + j 1 + i , A 7 = i + j 1 + i 1 + i + k 1 + j 1 + i + j 1 + j , A 8 = 1 + i + j k i 1 + j k i j 1 + j + k 1 + i , B 1 = i 5 j 1 + i 2 + i 1 2 i + j 1 2 i j k 2 + i j + k , B 2 = 1 + i + j k 1 + 3 i + j + k 1 + 2 i + j + k 1 + i + 2 j 1 + 3 i j 3 + 3 i , B 3 = 2 2 i + j k 3 + 2 j 3 k 3 2 j + k 1 k 2 i + 4 j 1 i + k 1 + 3 k 3 j 2 k 1 2 i j , B 4 = 2 + i j + 3 k 2 + i + 5 j 2 k 2 j 1 + 3 i + 2 j + 3 k 2 3 i + 2 j k 1 2 i + j 3 + 3 i + k 3 i + 5 j 2 j + 2 k , B 5 = i 1 + j + k 1 + i + k 0 i 1 + i + j , B 6 = 2 + 2 i + j + k 2 + 2 j + k 1 + i + j 1 + 2 i + j + 2 k 1 + 2 i + k 2 + 2 i + 2 j + k , B 7 = k i + j i + k i + k j 1 , B 8 = j j 1 + i + j + k 1 1 + j + k 1 + j , C 1 = 1 + i + k 1 + i + k 1 + i + k 1 + i k 0 , C 2 = i 1 0 1 i + j k i j i + k , C 3 = 0 k 0 1 + j + k i + j + k 1 + i + k , D 1 = i + j 1 + 2 i + j k 1 + 2 i + j 1 + 4 i + 4 k 1 + i + 4 k 1 + 2 i + j + 3 k , D 2 = i j + k 1 i + j 0 4 + 2 i + 3 k i j 3 k 3 i j + k , D 3 = i + j + 2 k 1 + j + 2 k i + j + 2 k 1 j + 2 k 3 + 2 k 2 + i + k 0 1 j 0 , D 4 = i + 2 k 3 + j + 5 k 4 + k 1 + j 2 6 i + 5 j + 4 k i + j + 3 k 1 + i + 3 j + k i + 7 k i + 3 k , D 5 = 7 i 2 j + 9 k 2 + 10 i + 13 j + 10 k 4 i + 6 j 7 + 8 i j + 7 k 12 + 2 i + 10 j k 5 + 5 i + j k 7 + 4 i 4 j + 6 k 8 + 2 i + 6 j 4 k 11 + 9 i + j 8 k .
Computing directly yields
r A i B i = r A i = 2 , r C j D j = r C j = 2 ( i = 1 , 4 ¯ , j = 1 , 3 ¯ ) , ( 26 ) = 4 , ( 27 ) = 10 , ( 28 ) = 10 , ( 29 ) = 10 , ( 30 ) = 10 , ( 31 ) = 8 , ( 32 ) = 8 , ( 33 ) = 8 , ( 34 ) = 8 , ( 35 ) = 24 .
All rank equations hold. Thus, according to Theorem 1, the system of matrix equations has a solution, and the general solution to the system can be expressed as
X = 1 + j 1 2 + k 1 + i j , Y = 1 + i j 0 1 2 i j + k , Z = 1 + i 1 + k 1 + j + k i + k i + k 0 , V = 1 2 + i k 1 + j 1 j + k i + j 0 , W = j 3 i + j + k j 1 + i i + j i + 2 j ,

5. The General Solution to the System (7) with η -Hermicity

As an application of the results of system (6), we study the necessary and sufficient conditions for system (7) to have a solution involving η -Hermicity and derive a formula of its general solution, where X , Y are η -Hermitian matrices.
Theorem 3.
Given A i , B j ( i = 1 , 2 , 7 , 8 , j = 1 , 2 , 5 , 8 ¯ ) , C 3 , D 3 , D 4 of appropriate dimensions over H . Set
A 22 = A 6 L A 4 , A 33 = A 7 L A 1 , A 44 = A 8 L A 2 ,
C 22 = D 5 A 6 A 4 B 4 A 7 A 1 ( A 1 ) η * + L A 1 C 1 η * ( C 1 ) η * A 7 η * A 8 A 2 ( A 2 ) η * + L A 2 C 2 η * ( C 2 ) η * ) A 8 η * ,
M 1 = R A 22 A 33 , M 2 = R A 22 A 44 , T = R A 22 C 22 R A 22 η * , M = R M 1 M 2 , S = M 2 L M .
Then, the following statements are equivalent:
( 1 ) System (7) has a solution.
( 2 )
R A 1 B 1 = 0 , R A 2 B 2 = 0 , R A 4 B 4 = 0 , R M 1 R M T = 0 , R A 22 T ( R A 44 ) η * = 0 .
( 3 )
r A 1 B 1 = r A 1 , r A 2 B 2 = r A 2 , r A 4 B 4 = r A 4 , r D 5 A 6 B 4 η * A 7 B 1 η * A 8 B 2 η * A 6 η * 0 A 4 η * 0 0 A 7 η * 0 0 A 1 η * 0 A 8 η * 0 0 0 A 2 η * B 4 A 4 0 0 0 = r A 6 η * A 4 η * 0 0 A 7 η * 0 A 1 η * 0 A 8 η * 0 0 A 1 η * + r A 4 A 6 , r D 5 A 6 A 7 B 4 η * A 8 B 2 η * A 6 η * 0 0 A 4 η * 0 A 8 η * 0 0 0 A 2 η * B 4 A 4 0 0 0 B 1 A 7 η * 0 A 1 0 0 = r A 6 η * A 4 η * 0 A 8 η * 0 A 2 η * + r A 6 A 7 A 4 0 0 A 1 .
Under these conditions, the general solution with η-Hermicity to the system (7) can be stated as
V = A 4 B 4 + L A 4 U 1 , X = A 1 B 1 + L A 1 B 1 η * ( A 1 ) η * + L A 1 U 2 L A 1 η * , Y = A 2 B 2 + L A 2 B 2 η * ( A 2 ) η * + L A 2 U 3 L A 2 η * , U 1 = A 22 ( C 22 A 33 U 2 A 33 η * A 44 U 3 A 44 η * ) A 22 W 2 A 22 η * + L A 22 W 1 , U 2 = M 1 T M 1 η * M 1 M 2 M T M 1 η * M 1 S M 2 T ( M η * ) M 2 η * M 1 η * M 1 S V 4 ( L M ) η * M 2 η * M 1 η * + L M 1 V 1 + V 2 ( L M ) η * , U 3 = M T M 2 η * + S S M 2 T M η * + L M L S V 3 + L M V 4 ( L M ) η * + V 5 ( L M 2 ) η * ,
where V i ( i = 1 , 5 ¯ ) and W j ( j = 1 , 3 ¯ ) are arbitrary matrices with appropriate sizes over H .
Proof. 
Since the solvability of the system (7) is equivalent to the system
A 4 V 1 = B 4 , V 2 ( A 4 ) η * = ( B 4 ) η * , V 2 = ( V 1 ) η * , A 1 X 1 = B 1 , X 1 A 1 η * = B 1 η * , X 1 = X 1 η * , A 2 Y 1 = B 2 , Y 1 A 2 η * = B 2 η * , Y 1 = Y 1 η * , A 6 V 1 + V 2 A 6 η * + A 7 X 1 A 7 η * + A 8 Y 1 A 8 η * = D 5 , D 5 = D 5 η * .
If system (7) has a solution, say, ( V , X , Y ) , then system (59) has a solution, ( V 1 , V 2 , X 1 , Y 1 ) = ( V , V η * , X , Y ) . Conversely, if system (59) has a solution ( V 1 , V 2 , X 1 , Y 1 ) , then
( V , X , Y ) = ( V 1 + V 2 η * 2 , X 1 + X 1 η * 2 , X 2 + X 2 η * 2 )
is the solution of (7). It follows from Corollary 1 that this proof can be completed. □

6. Conclusions

We established the solvability conditions for system (6) by using the Moore-Penrose inverses and ranks of the coefficient quaternion matrices in (6), and derived a formula of its general solution when it is solvable. In terms of applications, we derived the necessary and sufficient conditions for system (7) to have an η -Hermitian solution as well as the expression of the general solution. In addition, we used an algorithm and a numerical example to verify the main results of this paper. It is worth noting that the main results of (6) are available not only for R and C , but also any division ring. Moreover, inspired by [39], we can investigate the system (6) tensor equations over the quaternion algebra.

Author Contributions

Methodology, R.-N.W. and Q.-W.W.; software, R.-N.W. and L.-S.L.; writing—original draft preparation, Q.-W.W. and R.-N.W.; writing—review and editing, Q.-W.W., R.-N.W. and L.-S.L.; supervision, Q.-W.W.; project administration, Q.-W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by National Natural Science Foundation of China (11971294) and (12171369).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors have declared that there is no conflict of interest.

References

  1. Roth, W.E. The equations AXYB = C and AXXB = C in matrices. Proc. Am. Math. Soc. 1952, 3, 392–396. [Google Scholar] [CrossRef]
  2. Cvetković-Ilić, S.S.; Radenkovixcx, J.N.; Wang, Q.W. Algebraic conditions for the solvability to some systems of matrix equations. Linear Multilinear Algebra 2021, 69, 1579–1609. [Google Scholar] [CrossRef]
  3. Sasane, A. The Sylvester equation in Banach algebras. Linear Algebra Appl. 2021, 631, 1–9. [Google Scholar] [CrossRef]
  4. Wang, Q.W.; Rehman, A.; He, Z.H.; Zhang, Y. Constraint generalized Sylvester matrix equations. Automatica 2016, 69, 60–64. [Google Scholar] [CrossRef]
  5. Wang, L.; Wang, Q.W.; He, Z.H. The common solution of some matrix equations. Algebra Colloq. 2016, 23, 71–81. [Google Scholar] [CrossRef]
  6. Castelan, E.B.; da Silva, V.G. On the solution of a Sylvester equation appearing in descriptor systems control theory. Syst. Control Lett. 2005, 54, 109–117. [Google Scholar] [CrossRef]
  7. Zhang, Y.N.; Jiang, D.C.; Wang, J. A recurrent neural network for solving Sylvester equation with time-varying coefficients. IEEE Trans. Neural Netw. 2002, 13, 1053–1063. [Google Scholar] [CrossRef]
  8. Villareal, E.R.L.; Vargas, J.A.R.; Hemerly, E.M. Static output feedback stabilization using invariant subspaces and Sylvester equations. TEMA Tend. Mat. Apl. Comput. 2009, 10, 99–110. [Google Scholar] [CrossRef]
  9. Dmytryshyn, A.; Kågström, B. Coupled Sylvester-type matrix equations and block diagonalization. SIAM J. Matrix Anal. Appl. 2015, 36, 580–593. [Google Scholar] [CrossRef]
  10. Baksalary, J.K.; Kala, R. The matrix equation AXYB = C. Linear Algebra Appl. 1979, 25, 41–43. [Google Scholar] [CrossRef] [Green Version]
  11. Baksalary, J.K.; Kala, R. The matrix equation AXBCYD = E. Linear Algebra Appl. 1980, 30, 141–147. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, Q.W. A system of matrix equations and a linear matrix equation over arbitrary regular rings with identity. Linear Algebra Appl. 2004, 384, 43–54. [Google Scholar] [CrossRef]
  13. Kyrchei, I. Explicit representation formulas for the minimum norm least squares solutions of some quaternion matrix equations. Linear Algebra Appl. 2018, 438, 136–152. [Google Scholar] [CrossRef] [Green Version]
  14. Kyrchei, I. Determinantal representations of the Drazin inverse over the quaternion skew field with applications to some matrix equations. Appl. Math. Comput. 2014, 238, 193–207. [Google Scholar] [CrossRef]
  15. Liu, X. The η-anti-Hermitian solution to some classic matrix equations. Appl. Math. Comput. 2018, 320, 264–270. [Google Scholar] [CrossRef]
  16. Liu, X.; Zhang, Y. Consistency of split quaternion matrix equations AX* − XB = CY + D and XAX*B = CY + D. Adv. Appl. Clifford Algebr. 2019, 64, 1–20. [Google Scholar] [CrossRef]
  17. Assefa, D.; Mansinha, L.; Tiampo, K.F.; Rasmussen, H.; Abdella, K. Local quaternion Fourier transform and color image texture analysis. Signal Process. 2010, 90, 1825–1835. [Google Scholar] [CrossRef]
  18. Bihan, N.L.; Mars, J. Singular value decomposition of quaternion matrices: A new tool for vectorsensor signal processing. Signal Process. 2004, 84, 1177–1199. [Google Scholar] [CrossRef]
  19. He, Z.H.; Chen, C.; Wang, X.X. A simultaneous decomposition for three quaternion tensors with applications in color video signal processing. Anal. Appl. 2021, 19, 529–549. [Google Scholar] [CrossRef]
  20. Qi, L.; Luo, Z.Y.; Wang, Q.W.; Zhang, X.Z. Quaternion matrix optimization: Motivation and analysis. J. Optim. Theory Appl. 2021, 1–28. [Google Scholar] [CrossRef]
  21. Xu, Y.F.; Wang, Q.W.; Liu, L.S. A constrained system of matrix equations. Comput. Appl. Math. 2022, 41, 166. [Google Scholar] [CrossRef]
  22. Liu, L.S.; Wang, Q.W.; Cheng, J.F.; Xie, Y.Z. An exact solution to a quaternion matrix equation with an application. Symmetry 2022, 14, 375. [Google Scholar] [CrossRef]
  23. Liu, X.; Song, G.J.; Zhang, Y. Determinantal representations of the solutions to systems of generalized sylvester equations. Adv. Appl. Clifford Algebr. 2019, 30, 12. [Google Scholar] [CrossRef]
  24. Mehany, M.S.; Wang, Q.W. Three symmetrical systems of coupled Sylvester-like quaternion matrix equations. Symmetry 2022, 14, 550. [Google Scholar] [CrossRef]
  25. Rehman, A.; Wang, Q.W.; Ali, I.; Akram, M.; Ahmad, M.O. A constraint system of generalized Sylvester quaternion matrix equations. Adv. Appl. Clifford Algebr. 2017, 27, 3183–3196. [Google Scholar] [CrossRef]
  26. Wang, Q.W.; He, Z.H.; Zhang, Y. Constrained two-sided coupled Sylvester-type quaternion matrix equations. Automatica 2019, 101, 207–213. [Google Scholar] [CrossRef]
  27. Yu, S.W.; He, Z.H.; Qi, T.C.; Wang, X.X. The equivalence canonical form of five quaternion matrices with applications to imaging and Sylvester-type equations. J. Comput. Appl. Math. 2021, 393, 113494. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Wang, R.H. The exact solution of a system of quaternion matrix equations involving η-Hermicity. Appl. Math. Comput. 2013, 222, 201–209. [Google Scholar] [CrossRef]
  29. Wang, Q.W.; He, Z.H. Some matrix equations with applications. Linear Multilinear Algebra 2012, 60, 1327–1353. [Google Scholar] [CrossRef]
  30. Zhang, X. Characterization for the general solution to a system of matrix equations with quadruple variables. Appl. Math. Comput. 2014, 226, 274–287. [Google Scholar] [CrossRef]
  31. Took, C.C.; Mandic, D.P.; Zhang, F.Z. On the unitary diagonalisation of a special class of quaternion matrices. Appl. Math. Lett. 2011, 24, 1806–1809. [Google Scholar] [CrossRef] [Green Version]
  32. Took, C.C.; Mandic, D.P. Augmented second-order statistics of quaternion random signals. Signal Process. 2011, 91, 214–224. [Google Scholar] [CrossRef] [Green Version]
  33. He, Z.H.; Wang, Q.W. The η-bihermitian solution to a system of real quaternion matrix equations. Linear Multilinear Algebra 2014, 62, 1509–1528. [Google Scholar] [CrossRef]
  34. He, Z.H.; Wang, Q.W.; Zhang, Y. Simultaneous decomposition of quaternion matrices involving η-Hermicity with applications. Appl. Math. Comput. 2017, 298, 13–35. [Google Scholar] [CrossRef]
  35. Marsaglia, G.; Styan, G.P. Equalities and inequalities for ranks of matrices. Linear Multilinear Algebra 1974, 2, 269–292. [Google Scholar] [CrossRef]
  36. Buxton, J.N.; Churchouse, R.F. Tayler, A.B. Matrices Methods and Applications; Clarendon Press: Oxford, UK, 1990. [Google Scholar]
  37. Li, F.L.; Hu, X.Y.; Zhang, L. The generalized reflexive solution for a class of matrix equations (AX = B, XC = D). Acta Math. Sci. 2008, 28, 185–193. [Google Scholar]
  38. Liu, L.S.; Wang, Q.W.; Mahmoud, S.M. A Sylvester-type Hamilton quaternion matrix equation with an application. arXiv 2021, arXiv:2109.10045v2. [Google Scholar]
  39. Li, T.; Wang, Q.W.; Zhang, X.F. A Modified Conjugate Residual Method and Nearest Kronecker Product Preconditioner for the Generalized Coupled Sylvester Tensor Equations. Mathematics 2022, 10, 1730. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, R.-N.; Wang, Q.-W.; Liu, L.-S. Solving a System of Sylvester-like Quaternion Matrix Equations. Symmetry 2022, 14, 1056. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14051056

AMA Style

Wang R-N, Wang Q-W, Liu L-S. Solving a System of Sylvester-like Quaternion Matrix Equations. Symmetry. 2022; 14(5):1056. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14051056

Chicago/Turabian Style

Wang, Ruo-Nan, Qing-Wen Wang, and Long-Sheng Liu. 2022. "Solving a System of Sylvester-like Quaternion Matrix Equations" Symmetry 14, no. 5: 1056. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14051056

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop