Next Article in Journal
Generalized Dynamic Inequalities of Copson Type on Time Scales
Previous Article in Journal
An Improved Grey Wolf Optimizer with Multi-Strategies Coverage in Wireless Sensor Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dual Quaternion Matrix Equation AXB = C with Applications

1
Department of Mathematics, Newtouch Center for 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.
Symmetry 2024, 16(3), 287; https://doi.org/10.3390/sym16030287
Submission received: 4 February 2024 / Revised: 26 February 2024 / Accepted: 28 February 2024 / Published: 1 March 2024
(This article belongs to the Section Mathematics)

Abstract

:
Dual quaternions have wide applications in automatic differentiation, computer graphics, mechanics, and others. Due to its application in control theory, matrix equation A X B = C has been extensively studied. However, there is currently limited information on matrix equation A X B = C regarding the dual quaternion algebra. In this paper, we provide the necessary and sufficient conditions for the solvability of dual quaternion matrix equation A X B = C , and present the expression for the general solution when it is solvable. As an application, we derive the ϕ -Hermitian solutions for dual quaternion matrix equation A X A ϕ = C , where the ϕ -Hermitian extends the concepts of Hermiticity and η -Hermiticity. Lastly, we present a numerical example to verify the main research results of this paper.
MSC:
15A03; 15A09; 15A24; 15B33; 15B57

1. Introduction

Let R denote the set of real numbers and H m × n stand for the space of all m × n matrices over quaternions
H = { u 0 + u 1 i + u 2 j + u 3 k | i 2 = j 2 = k 2 = i j k = 1 , u 0 , u 1 , u 2 , u 3 R } .
The symbols r ( A ) , A * , I, and 0 are defined by the rank of a given quaternion matrix A, the conjugate transpose of A, identity matrix, and zero matrix with appropriate sizes, respectively. The Moore–Penrose inverse of A H l × k is denoted as A , which is defined as the solution of A Y A = A , Y A Y = Y , ( A Y ) * = A Y , and ( Y A ) * = Y A . Moreover, let L A = I A A and R A = I A A represent two projectors along A.
Since Hamilton’s discovery of quaternions in 1843, quaternions and quaternion matrices have found a large amount of practical applications in fields such as computer science, statistics, quantum physics, signal and color image processing, flight mechanics, aerospace technology, and so on (see, e.g., [1,2,3,4]). Furthermore, quaternion matrix equations also have significant applications in many fields, such as system and control theory.
Up to this point, matrix equations have witnessed a large number of papers proposing various methods for solving some matrix equations (see, e.g., [5,6,7,8,9,10]). The classical matrix equation
A X B = C ,
has been studied by many authors. Ben-Israel and Greville [11] eatablished the necessary and sufficient conditions for the solvability of matrix Equation (1). In 2003, Liao and Bai [12] investigated the least-squares solution of matrix Equation (1) over symmetric positive semidefinite matrices. Huang et al. [13] provided the skew-symmetric solution and the optimal approximate solution of matrix Equation (1). Peng [14] derived the centro-symmetric solution of matrix Equation (1). Deng et al. [15] studied the general expressions regarding the Hermitian solutions of matrix Equation (1). Xie and Wang [16] considered the reducible solution to matrix Equation (1) when it is solvable. As a special case of matrix Equation (1), the Hermitian solution X to matrix equation
A X A * = B
has attracted extensive attention (see, e.g., [17,18]). Baksalary [19] and Größ [20] studied the nonnegative definite and positive definite solutions to matrix Equation (2), respectively. For η { i , j , k } , a quaternion matrix A is called η -Hermitian and η -skew-Hermitian if A = A η * and A = A η * , respectively, where A η * = η A * η [21]. The utilization of η -Hermitian matrices in linear modeling is extensively recognized [22]. Kyrchei [23] derived the explicit determinantal representation formulas of η -Hermitian and η -skew-Hermitian solutions to the quaternion matrix equation
A X A η * = B .
It is well-known that dual numbers and dual quaternions have wide applications in computer graphics, automatic differentiation, geometry, mechanics, rigid body motions, and robotics (see, e.g., [24,25,26]). For the related definitions of dual numbers and dual quaternions, please see Section 2.
So far, there has been little information on matrix Equation (1) regarding dual quaternion algebra. Motivated by the work mentioned above, in this paper, we aim to investigate the general solution of dual quaternion matrix Equation (1) by using Moore–Penrose inverses and ranks of matrices. Since the ϕ -Hermitian serves as an extended form of both Hermiticity and η -Hermiticity over the quaternions [27], we also provide the definition of ϕ -Hermiticity over the dual quaternions. As an application, we establish the ϕ -Hermitian solution of a special dual quaternion matrix equation
A X A ϕ = C , C ϕ = C .
Further details regarding ϕ -Hermitian matrices will be illustrated in Section 2.
This paper is organized as follows. In Section 2, we provide an overview of essential definitions and lemmas that will be applied in the subsequent sections. In Section 3, we establish some necessary and sufficient conditions for solvability regarding dual quaternion matrix Equation (1) and consider some special cases of dual quaternion matrix Equation (1). As an application, we investigate the ϕ -Hermitian solution of dual quaternion matrix Equation (3) in Section 4. In Section 5, we present a numerical example to illustrate the results of this paper. Finally, a brief conclusion is provided in Section 6.

2. Preliminaries

In this section, we review some definitions of dual numbers, dual quaternions, and related propositions. Moreover, we introduce the definitions of dual quaternion matrix and ϕ -Hermitian matrix, which are fundamental for obtaining the main results.
Definition 1
([28]). Suppose that x 0 , x 1 R ; we say x is a dual number if x has the form
x = x 0 + x 1 ϵ ,
where ϵ is the infinitesimal unit, satisfying ϵ 2 = 0 .
We call x 0 the real part or the standard part of x, x 1 as the dual part or the infinitesimal part of x. The infinitesimal unit ϵ is commutative in multiplication with real numbers, complex numbers, and quaternions. The set of dual numbers is denoted by
D = { x = x 0 + x 1 ϵ | ϵ 2 = 0 , x 0 , x 1 R } .
Assume that x = x 0 + x 1 ϵ , y = y 0 + y 1 ϵ D ; we have x = y if x 0 = y 0 and x 1 = y 1 ; regarding addition and multiplication, there is
x + y = x 0 + y 0 + ( x 1 + y 1 ) ϵ , x y = x 0 y 0 + ( x 0 y 1 + x 1 y 0 ) ϵ .
Definition 2
([28]). Let z 0 , z 1 H . We say z is a dual quaternion if z has the form
z = z 0 + z 1 ϵ ,
where ϵ is the infinitesimal unit, satisfying ϵ 2 = 0 , z 0 , z 1 as the real part and the dual part of z, respectively.
The collection of dual quaternions is denoted by
DQ = { z = z 0 + z 1 ϵ | ϵ 2 = 0 , z 0 , z 1 H } .
Now, we introduce the definition of dual quaternion matrix. Let X 0 , X 1 H m × n . X is said to be a dual quaternion matrix if X has the form X = X 0 + X 1 ϵ ; the set of dual quaternion matrices is denoted by
DQ m × n = { X = X 0 + X 1 ϵ | ϵ 2 = 0 , X 0 , X 1 H m × n } .
The conjugate transpose of X is defined as X * = X 0 * + X 1 * ϵ . For Y = Y 0 + Y 1 ϵ DQ m × n , by analogy, we have X = Y if X 0 = Y 0 and X 1 = Y 1 ; furthermore,
X + Y = X 0 + Y 0 + ( X 1 + Y 1 ) ϵ , X Y = X 0 Y 0 + ( X 0 Y 1 + X 1 Y 0 ) ϵ .
To facilitate our study on the ϕ -Hermitian, we first review the concept of nonstandard involution over quaternions, and then proceed to generalize it to dual quaternions.
Definition 3
([27]). A map ϕ : H H is called an antiendomorphism if ϕ ( p q ) = ϕ ( q ) ϕ ( p ) and ϕ ( p + q ) = ϕ ( p ) + ϕ ( q ) for all p , q H . An antiendomorphism ϕ is called an involution if ϕ ( ϕ ( p ) ) = p for every p H .
Definition 4
([27]). Under the basis 1 , i , j , k , an involution ϕ is called nonstandard if and only if ϕ can be expressed as a real matrix
ϕ = 1 0 0 T ,
where T is a 3 × 3 real orthoganal symmetric matrix with eigenvalues 1 , 1 , 1 .
Proposition 1
([27]). Let z H . Then, every nonstandard involution ϕ of H has the form ϕ ( z ) = β 1 z * β for some β H with β 2 = 1 .
Definition 5
([27]). For a nonstandard involution ϕ, Z H m × n , we denote by Z ϕ the matrix obtained by applying ϕ entrywise to the transposed matrix of Z.
For example, if ϕ is such that ϕ ( i ) = i , ϕ ( j ) = j , ϕ ( k ) = k , then
1 + k j i 3 + k j + k i j ϕ = 1 + k i j + k j 3 + k i + j .
Definition 6
([27]). For a nonstandard involution ϕ, a quaternion matrix Z is said to be ϕ-Hermitian if Z = Z ϕ . In fact, Z ϕ = β 1 Z * β for some β H with β 2 = 1 .
Proposition 2.
Let Z H m × n , β H , and β 2 = 1 . Then,
(1) 
( Z ϕ ) = ( Z ) ϕ ,
(2) 
( L Z ) ϕ = R Z ϕ ,
(3) 
( R Z ) ϕ = L Z ϕ .
Proof. 
Regarding the proof of (1), it is obvious that β 1 = β , β * = β ; by the definition of Moore–Penrose inverse, we obtain
Z ϕ ( Z ) ϕ Z ϕ = ( β Z * β ) [ β ( Z ) * β ] ( β Z * B ) = β Z * ( Z ) * Z * β = β ( Z Z Z ) * β = β Z * β = Z ϕ , ( Z ) ϕ Z ϕ ( Z ) ϕ = [ β ( Z ) * β ] ( β Z * β ) [ β ( Z ) * β ] = β ( Z ) * Z * ( Z ) * β = β ( Z Z Z ) * β = β ( Z ) * β = ( Z ) ϕ , [ Z ϕ ( Z ) ϕ ] * = [ ( β Z * β ) ( β ( Z ) * β ) ] * = [ β Z * ( Z ) * β ] * = ( β Z Z β ) * = β ( Z Z ) * β = ( β Z * β ) [ β ( Z ) * β ] = Z ϕ ( Z ) ϕ , [ ( Z ) ϕ Z ϕ ] * = [ β ( Z ) * β ( β Z * β ) ] * = [ β ( Z ) * Z * β ] * = ( β Z Z β ) * = β ( Z Z ) * β = [ β ( Z ) * β ] ( β Z * β ) = ( Z ) ϕ Z ϕ .
Based on this, we can deduce that ( Z ) ϕ is the Moore–Penrose inverse of Z ϕ .
For (2), we have
( L Z ) ϕ = ( I Z Z ) ϕ = I Z ϕ ( Z ) ϕ = I Z ϕ ( Z ϕ ) = R Z ϕ .
In a similar vein to ( 2 ) , we can offer a demonstration for ( 3 ) ; therefore, we omit it here. □
By analogy, we propose the definition of ϕ -Hermiticity with respect to dual quaternion matrix, where ϕ is a nonstandard involution.
Definition 7.
For X = X 0 + X 1 ϵ DQ m × n , X is called ϕ-Hermitian matrix if X = X ϕ , where
X ϕ : = β 1 X * β = β 1 X 0 * β + β 1 X 1 * β ϵ = X 0 ϕ + X 1 ϕ ϵ ,
with β H and β 2 = 1 .
For example, if ϕ is such that ϕ ( i ) = i , ϕ ( j ) = j , ϕ ( k ) = k ,
X = 1 2 + i 2 + i 3 + k ϕ + 3 + i k k 2 i + 3 k ϕ ϵ = 1 2 + i 2 + i 3 + k + 3 + i k k 2 i + 3 k ϵ ,
then X is a ϕ -Hermitian matrix.
Proposition 3
([27]). Let X , Y H m × n , α , β H . Then,
(1) 
( α X + β Y ) ϕ = X ϕ ϕ ( α ) + Y ϕ ϕ ( β ) ,
(2) 
( X α + Y β ) ϕ = ϕ ( α ) X ϕ + ϕ ( β ) Y ϕ ,
(3) 
( X Y ) ϕ = Y ϕ X ϕ ,
(4) 
( X ϕ ) ϕ = X .
Having outlined the properties of ϕ -Hermitian matrix over quaternions, we now present the corresponding properties of ϕ -Hermitian matrix over the dual quaternion algebra.
Proposition 4.
Let X , Y DQ n × n . Then,
(1) 
( X + Y ) ϕ = X ϕ + Y ϕ ,
(2) 
( X Y ) ϕ = Y ϕ X ϕ ,
(3) 
( X ϕ ) ϕ = X .
Proof. 
By the algebraic properties of ϕ , we have
( X + Y ) ϕ = [ ( X 0 + Y 0 ) + ( X 1 + Y 1 ) ϵ ] ϕ = ( X 0 + Y 0 ) ϕ + ( X 1 + Y 1 ) ϕ ϵ = X 0 ϕ + Y 0 ϕ + X 1 ϕ ϵ + Y 1 ϕ ϵ = X 0 ϕ + X 1 ϕ ϵ + Y 0 ϕ + Y 1 ϕ ϵ = X ϕ + Y ϕ .
In relation to (2), we obtain
( X Y ) ϕ = [ X 0 Y 0 + ( X 0 Y 1 + X 1 Y 0 ) ϵ ] ϕ = ( X 0 Y 0 ) ϕ + ( X 0 Y 1 + X 1 Y 0 ) ϕ ϵ = Y 0 ϕ X 0 ϕ + Y 1 ϕ X 0 ϕ ϵ + Y 0 ϕ X 1 ϕ ϵ = ( Y 0 ϕ + Y 1 ϕ ϵ ) ( X 0 ϕ + X 1 ϕ ϵ ) = Y ϕ X ϕ .
In terms of (3), we have
( X ϕ ) ϕ = [ ( X 0 + X 1 ϵ ) ϕ ] ϕ = ( X 0 ϕ + X 1 ϕ ϵ ) ϕ = X 0 + X 1 ϵ = X .
Now, we provide a few lemmas, which are basic tools for obtaining the key outcomes.
Lemma 1
([29]). Suppose that A, B, and C are provided for matrices with the adequate dimensions over H ; then, quaternion matrix Equation (1) is consistent if and only if
R A C = 0 , C L B = 0 .
In this case, the general solution can be expressed as
X = A C B + L A U + V R B ,
where U , V are any matrices over H with appropriate dimensions.
Lemma 2
([16]). Let A 1 , A 2 , B 1 , B 2 , and C 1 have matrices with appropriate sizes. Set
A = R A 1 C , B = B 1 L B 2 , M = R A 1 A 2 , C 1 = C L B 2 .
Then, the following descriptions are equivalent:
(1) 
The quaternion matrix equation
A 1 X 1 B 1 + A 1 X 2 B 2 + A 2 X 3 B 2 = C
is consistent.
(2) 
R M A = 0 , R A 1 C L B 2 = 0 , C 1 L B = 0 .
(3) 
r A 1 A 2 C = r A 1 A 2 , r B 2 0 C A 1 = r B 2 + r A 1 , r C 1 B 1 B 2 = r B 1 B 2 .
In this case, the general solution to (5) can be expressed as follows:
X 1 = A 1 C 1 B + L A 1 V 1 + V 2 R B , X 2 = A 1 ( C A 1 X 1 B 1 A 2 X 3 B 2 ) B 2 + T 1 R B 2 + L A 1 T 2 , X 3 = M A B 2 + L M U 1 + U 2 R B 2 ,
where U 1 , U 2 , V 1 , V 2 , T 1 , and T 2 are arbitrary matrices over H with appropriate sizes.
The following lemma, originally derived by Marsaglia and Styan [30], can be extended to H .
Lemma 3.
Assume that A H m × n ,   B H m × k , C H l × n , D H j × k , and E H l × i ; then, we have the following rank equality:
(1) 
r R A B = r A B r ( A ) .
(2) 
r C L B = r B C r ( B ) .
(3) 
r R A B L C = r B A C 0 r ( A ) r ( C ) .
(4) 
r R A B C = r B A C 0 r ( A ) .
(5) 
r A L B C = r A C B 0 r ( B ) .
(6) 
r A B L D R E C 0 = r A B 0 C 0 E 0 D 0 r ( D ) r ( E ) .

3. The Solution of Matrix Equation (1)

In this section, we establish the necessary and sufficient conditions for the solvability of dual quaternion matrix Equation (1) and provide the expressions of its general solution. Additionally, we investigate some special cases of dual quaternion matrix Equation (1).
Theorem 1.
Let A = A 0 + A 1 ϵ DQ m × n , B = B 0 + B 1 ϵ DQ k × l , C = C 0 + C 1 ϵ DQ m × l be known. Put
A 2 = A 1 L A 0 , B 2 = R B 0 B 1 , C 11 = A 0 A 0 C 0 B 0 B 1 ,
C 22 = A 1 A 0 C 0 B 0 B 0 , C 2 = C 1 C 11 C 22 ,
M = R A 0 A 2 , N = R A 0 C 2 , E = B 2 L B 0 , F = C 2 L B 0 .
Then, the following statements are equivalent:
(1) 
Dual quaternion matrix Equation (1) is consistent.
(2) 
R A 0 C 0 = 0 , C 0 L B 0 = 0 ,
R M N = 0 , R A 0 C 2 L B 0 = 0 , F L E = 0 .
(3) 
r A 0 C 0 = r A 0 , r B 0 C 0 = r B 0 ,
r A 1 A 0 C 1 A 0 0 C 0 = r A 1 A 0 A 0 0 ,
r C 1 A 0 B 0 0 = r A 0 + r B 0 ,
r B 1 B 0 B 0 0 C 1 C 0 = r B 1 B 0 B 0 0 .
In this case, the general solution X of dual quaternion matrix Equation (1) can be expressed as X = X 0 + X 1 ϵ , where
X 0 = A 0 C 0 B 0 + L A 0 U + V R B 0 , X 1 = A 0 ( C 2 A 0 V B 2 A 2 U B 0 ) B 0 + W 1 R B 0 + L A 0 W 2 , U = M N B 0 + L M Q 1 + Q 2 R B 0 , V = A 0 F E + L A 0 Q 3 + Q 4 R E ,
and Q i ( i = 1 , 4 ¯ ) , W i ( i = 1 , 2 ¯ ) are arbitrary matrices over H with appropriate dimensions.
Proof. 
( 1 ) ( 2 ) : Suppose that dual quaternion matrix Equation (1) is solvable and its solution is X DQ m × n , which can be expressed as
X = X 0 + X 1 ϵ ,
substituting (16) into (1), by the definition of equality of dual quaternion matrices, we can obtain that dual quaternion matrix Equation (1) is equivalent to the system of quaternion matrix equations
{ A 0 X 0 B 0 = C 0 , A 0 X 0 B 1 + A 0 X 1 B 0 + A 1 X 0 B 0 = C 1 .
The structure of the proof goes as follows. We first prove that ( 1 ) ( 2 ) and illustrate the general solution of (17) with the form of (15), and then we prove ( 2 ) ( 3 ) .
We divide the system (17) into the following:
A 0 X 0 B 0 = C 0 ,
and
A 0 X 0 B 1 + A 0 X 1 B 0 + A 1 X 0 B 0 = C 1 .
By Lemma 1, we obtain that (18) is consistent if and only if
R A 0 C 0 , C 0 L B 0 = 0 .
In this case, the general solution of (18) can be written as
X 0 = A 0 C 0 B 0 + L A 0 U + V R B 0 ,
where U , V are any matrices over H with appropriate sizes.
Substituting Equation (20) into Equation (19) provides
A 0 V B 2 + A 0 X 1 B 0 + A 2 U B 0 = C 2 ,
where A 2 , B 2 , and C 2 are defined by (6) and (7). Using Lemmas 2 to (21), we know that matrix Equation (21) is solvable if and only if
R M N = 0 , R A 0 C 2 L B 0 = 0 , F L E = 0 ,
where M , N , E , and F are provided by (8). In this case, the general solution of (21) can be expressed as
X 1 = A 0 ( C 2 A 0 V B 2 A 2 U B 0 ) B 0 + W 1 R B 0 + L A 0 W 2 ,
U = M N B 0 + L M Q 1 + Q 2 R B 0 ,
V = A 0 F E + L A 0 Q 3 + Q 4 R E ,
where Q 1 , Q 2 , Q 3 , Q 4 , W 1 , and W 2 are any matrices with the suitable dimensions over H . To sum up, we have shown that matrix Equation (1) has a dual solution X DQ m × n if and only if ( 2 ) holds.
( 2 ) ( 3 ) : We divide it into two parts to prove its equivalence.
  • Firstly, we prove that (9) holds if and only if (11) holds. According to Lemma 3, it is easy to verify that ( 9 ) ( 11 ) .
  • Now, we turn to prove that ( 10 ) ( 12 ) ( 14 ) . Let X 0 = A 0 C 0 B 0 . Then, it is easy to verify that X 0 is a particular solution to the matrix equation A 0 X 0 B 0 = C 0 . By Lemma 3 and block elementary operations, we obtain
    R M N = 0 r ( R M N ) = 0 r R A 0 A 2 R A 0 C 2 = r ( R A 0 A 2 ) , r A 0 A 2 C 2 = r A 0 A 2 , r A 0 A 1 L A 0 C 2 = r A 0 A 1 L A 0 , r A 1 A 0 C 2 A 0 0 0 = r A 1 A 0 A 0 0 , r A 1 A 0 C 1 A 0 A 0 C 0 B 0 B 1 A 1 A 0 C 0 B 0 B 0 A 0 0 0 = r A 1 A 0 A 0 0 , r A 1 A 0 C 1 A 0 0 A 0 A 0 C 0 B 0 B 0 = r A 1 A 0 A 0 0 , r A 1 A 0 C 1 A 0 0 C 0 = r A 1 A 0 A 0 0 ,
    R A 0 C 2 L B 0 = 0 r ( R A 0 C 2 L B 0 ) = 0 r C 2 A 0 B 0 0 = r ( A 0 ) + r ( B 0 ) , r C 1 C 11 C 22 A 0 B 0 0 = r ( A 0 ) + r ( B 0 ) , r C 1 A 0 A 0 C 0 B 0 B 1 A 1 A 0 C 0 B 0 B 0 A 0 B 0 0 = r ( A 0 ) + r ( B 0 ) , r C 1 A 0 B 0 0 = r ( A 0 ) + r ( B 0 ) , F L E = 0 r ( F L E ) = 0 r E F = r ( E ) r B 2 L B 0 C 2 L B 0 = r ( B 2 L B 0 ) , r B 0 B 2 C 2 = r B 0 B 2 r B 0 R B 0 B 1 C 2 = r B 0 R B 0 B 1 , r B 1 B 0 B 0 0 C 2 0 = r B 1 B 0 B 0 0 , r B 1 B 0 B 0 0 C 1 A 0 A 0 C 0 B 0 B 1 A 1 A 0 C 0 B 0 B 0 0 = r B 1 B 0 B 0 0 , r B 1 B 0 B 0 0 C 1 A 0 A 0 C 0 B 0 B 0 = r B 1 B 0 B 0 0 , r B 1 B 0 B 0 0 C 1 C 0 = r B 1 B 0 B 0 0 .
Now, we consider some special cases of dual quarernion matrix Equation (1).
Corollary 1
([31]). Assume that A = A 0 + A 1 ϵ DQ m × n , C = C 0 + C 1 ϵ DQ m × l are given. Put
A 2 = A 1 L A 0 , C 22 = A 1 A 0 C 0 , C 2 = C 1 C 22 , M = R A 0 A 2 , N = R A 0 C 2 .
Then, the following statements are equivalent:
(1) 
Dual quaternion matrix equation A X = C is consistent.
(2) 
R A 0 C 0 = 0 , R M N = 0 .
(3) 
r A 0 C 0 = r A 0 , r A 1 A 0 C 1 A 0 0 C 0 = r A 1 A 0 A 0 0 .
In this case, the general solution X of dual quaternion matrix equation A X = C can be expressed as X = X 0 + X 1 ϵ , where
X 0 = A 0 C 0 + L A 0 U , X 1 = A 0 ( C 2 A 2 U ) + L A 0 W 1 , U = M N + L M W 2 ,
and W 1 , W 2 are arbitrary matrices over H with appropriate dimensions.
Corollary 2
([31]). Let B = B 0 + B 1 ϵ DQ k × l , C = C 0 + C 1 ϵ DQ m × l be known. Denote
B 2 = R B 0 B 1 , C 11 = C 0 B 0 B 1 , C 2 = C 1 C 11 , E = B 2 L B 0 , F = C 2 L B 0 .
Then, the following statements are equivalent:
(1) 
Dual quaternion matrix equation X B = C is consistent.
(2) 
C 0 L B 0 = 0 , F L E = 0 .
(3) 
r B 0 C 0 = r B 0 , r B 1 B 0 B 0 0 C 1 C 0 = r B 1 B 0 B 0 0 .
In this case, the general solution can be expressed as X = X 0 + X 1 ϵ , where
X 0 = C 0 B 0 + V R B 0 , X 1 = ( C 2 V B 2 ) B 0 + W 1 R B 0 , V = F E + W 2 R E ,
and W 1 , W 2 are arbitrary matrices over H with appropriate dimensions.
Remark 1.
Matrix equations A X = B and X C = D have significant applications in eigenvalue problems, image processing, and solving linear systems. However, the matrix equation A X B = C is a general case that encompasses either matrix equation A X = B or X C = D . Therefore, the applications regarding matrix equations A X = B and X C = D are applicable to matrix equation A X B = C .

4. Applications

As an application of Theorem 1, we can investigate dual quaternion matrix Equation (3).
Theorem 2.
Suppose that A = A 0 + A 1 ϵ DQ m × n , C = C 0 + C 1 ϵ = C ϕ DQ m × m are provided; denote
B 2 = R A 0 ϕ A 1 ϕ , C 11 = A 0 A 0 C 0 ( A 0 ϕ ) A 1 ϕ , C 22 = A 1 A 0 C 0 ( A 0 ϕ ) A 0 ϕ , C 2 = C 1 C 11 C 22 , M = R A 0 B 2 ϕ , N = R A 0 C 2 .
Then, the following statements are equivalent:
(1) 
Dual quaternion matrix Equation (3) is consistent.
(2) 
The following equalities are satisfied:
R A 0 C 0 = 0 , R M N = 0 , R A 0 C 2 L A 0 ϕ = 0 .
(3) 
The following rank equalities hold:
r A 0 C 0 = r A 0 , r A 1 A 0 C 1 A 0 0 C 0 = r A 1 A 0 A 0 0 , r C 1 A 0 A 0 ϕ 0 = r A 0 + r A 0 ϕ = 2 r ( A 0 ) .
In this case, the general solution X of (3) can be expressed as X = X 0 + X 1 ϵ , where
X 0 = X 0 ˜ + X 0 ˜ ϕ 2 , X 1 = X 1 ˜ + X 1 ˜ ϕ 2
and
X 0 = A 0 C 0 ( A 0 ϕ ) + L A 0 U + V R A 0 ϕ , X 1 = A 0 ( C 2 A 0 V B 2 B 2 ϕ U A 0 ϕ ) ( A 0 ϕ ) + W 1 R A 0 ϕ + L A 0 W 2 , U = M N ( A 0 ϕ ) + L M Q 1 + Q 2 R A 0 ϕ , V = A 0 N ϕ ( M ϕ ) + L A 0 Q 3 + Q 4 R M ϕ ,
Q i ( i = 1 , 4 ¯ ) , W i ( i = 1 , 2 ¯ ) are any matrices with appropriate dimensions over H .
Proof. 
By using the definitions of equality of dual quaternion matrices and dual quaternion matrix multiplication, we can conclude that the consistency of dual quaternion matrix Equation (3) is contingent on the existence of the solutions to the system of quaternion matrix equations
A 0 X 0 ˜ A 0 ϕ = C 0 , A 0 X 0 ˜ A 1 ϕ + A 0 X 1 ˜ A 0 ϕ + A 1 X 0 ˜ A 0 ϕ = C 1 .
In fact, if matrix Equation (3) has a ϕ -Hermitian solution X = X 0 + X 1 ϵ , it is obvious that X 0 and X 1 must be solutions to (33). Conversely, if the system (33) has solutions X 0 ˜ and X 1 ˜ , then matrix Equation (3) has solution X = X 0 + X 1 ϵ , where
X 0 = X 0 ˜ + X 0 ˜ ϕ 2 , X 1 = X 1 ˜ + X 1 ˜ ϕ 2 .
According to Theorem 1, we can present the necessary and sufficient conditions for the solvability of (33), along with the general expression for its solutions. □

5. Numerical Example

Now, we provide a numerical example to illustrate the main results of this paper.
Example 1.
Given the dual quaternion matrices:
A = A 0 + A 1 ϵ = 2 i + k 3 i + j j 0 3 j 4 k i + k + 2 3 i + k i 3 k i j 0 4 i + j ϵ , B = B 0 + B 1 ϵ = 1 + i + j j 3 i + k 0 k + j 0 + i 2 k j + 3 k 0 i k i + k ϵ , C = C 0 + C 1 ϵ = 6 3 j + 3 k 9 + 3 i + 4 j 6 k 9 + 9 i 3 j 3 k 1 + 2 i j i k 4 + 2 j 4 + 10 i 2 j + 6 k 7 i j 3 k 14 16 i + 2 j 8 k + 2 i + 9 j + 5 k 16 + 28 i 16 j 18 k 2 + 14 i 6 j 10 k 2 + 5 i + j + 7 k 3 + i 6 j + k 9 10 i 3 j 10 k 11 + 24 i + 17 j 2 k 11 6 i 11 j + 4 k 42 6 i 12 j + 10 k ϵ .
Computing directly yields
r A 0 C 0 = r ( A 0 ) = 2 , r B 0 C 0 = r ( B 0 ) = 2 , r A 1 A 0 C 1 A 0 0 C 0 = r A 1 A 0 A 0 0 = 4 , r C 1 A 0 B 0 0 = r A 0 + r B 0 = 4 , r B 1 B 0 B 0 0 C 1 C 0 = r B 1 B 0 B 0 0 = 4 .
All rank equations are satisfied and a solution of dual quaternion matrix Equation (1) can be expressed as
X = X 0 + X 1 ϵ = i + k 0 1 j 2 k + j + 2 k 1 i 2 + 3 i 4 k ϵ .

6. Conclusions

Matrix equations A X = B and X C = D have specific applications in areas such as eigenvalue problems, image processing, and linear system solving. On the other hand, A X B = C is a more general matrix equation that has broader use. In this paper, we have established the solvability conditions for dual quaternion matrix Equation (1) by using Moore–Penrose inverses and ranks of matrices; we have also derived the expressions of its general solution to (1) when the solvability conditions are met. As special cases, some dual quaternion matrix equations have also been discussed. Moreover, we have investigated the ϕ -Hermitian matrix over dual quaternion algebra and provided its related properties. As an application of the aforementioned research, we have considered a special case of (1) and provided the ϕ -Hermitian solutions to (3). Finally, we have presented an example to illustrate the main results. In the future, we will focus on researching more complex matrix and tensor equations over the dual quaternion algebra.

Author Contributions

Methodology, Y.C. and Q.-W.W.; software, Y.C. and L.-M.X.; writing—original draft preparation, Q.-W.W. and Y.C.; writing—review and editing, Q.-W.W., Y.C. and L.-M.X.; 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 is supported by the National Natural Science Foundation of China (No. 12371023).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. 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]
  2. Cyrus, J.; Clive, C.T.; Danilo, P.M. A class of quaternion valued affine projection algorithms. Signal Process. 2013, 93, 1712–1723. [Google Scholar] [CrossRef]
  3. Bihan, N.L.E.; Mars, J. Singular value decomposition of matrices of quaternions: A new tool for vector-sensor signal processing. IEEE Trans. Signal Process. 2004, 84, 1177–1199. [Google Scholar] [CrossRef]
  4. Qi, L.Q.; Luo, Z.Y.; Wang, Q.W.; Zhang, X.Z. Quaternion matrix optimization: Motivation and analysis. J. Optim. Theory Appl. 2022, 193, 621–648. [Google Scholar] [CrossRef]
  5. Xu, X.; Wang, Q.W. The consistency and the general common solution to some quaternion matrix equations. Ann. Funct. Anal. 2023, 14, 53. [Google Scholar] [CrossRef]
  6. Dehghan, M.; Hajarian, M. On the generalized reflexive and anti-reflexive solutions to a system of matrix equations. Linear Algebra Appl. 2012, 437, 2793–2812. [Google Scholar] [CrossRef]
  7. 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]
  8. Dmytryshyn, A.; Kagstrom, B. Coupled Sylvester-type matrix equations and block diagonalization. SIAM J. Matrix Anal. Appl. 2015, 36, 580–593. [Google Scholar] [CrossRef]
  9. Zhang, H.; Yin, H. Conjugate gradient least squares algorithm for solving the generalized coupled sylvester matrix equations. Comput. Math. Appl. 2017, 12, 2529–2547. [Google Scholar] [CrossRef]
  10. Liu, L.S.; Wang, Q.W.; Chen, J.F. An exact solution to a quaternion matrix equation with an application. Symmetry 2022, 14, 375. [Google Scholar] [CrossRef]
  11. Ben-Israel, A.; Greville, T.N.E. Generalized Inverses: Theory and Application; John Wiley and Sons: New York, NY, USA, 1974. [Google Scholar] [CrossRef]
  12. Liao, A.P.; Bai, Z.Z. The constrained solutions of two matrix equations. Acta Math. Appl. Sin. Engl. Ser. 2002, 18, 671–678. [Google Scholar] [CrossRef]
  13. Huang, G.X.; Yin, F.; Guo, K. An iterative method for the skew-symmetric solution and the optimal approximate solution of the matrix equation AXB=C. J. Comput. Appl. Math. 2008, 212, 231–244. [Google Scholar] [CrossRef]
  14. Peng, Z.Y. The centro-symmetric solutions of linear matrix equation AXB = C and its optimal approximation. Chin. J. Eng. Math. 2003, 20, 60–64. [Google Scholar] [CrossRef]
  15. Deng, Y.B.; Bai, Z.Z.; Gao, Y.H. Iterative orthogonal direction methods for Hermitian minimum norm solutions of two consistent matrix equations. Numer. Linear Algebra Appl. 2006, 13, 801–823. [Google Scholar] [CrossRef]
  16. Xie, M.Y.; Wang, Q.W. The reducible solution to a quaternion tensor equation. Front. Math. China 2020, 15, 1047–1070. [Google Scholar] [CrossRef]
  17. Cvetković-IIić, D.S.; Dajić, A.; Koliha, J.J. Positive and real-positive solutions to the equation axa* = c in C*-algebras. Linear Multilinear Algebra 2007, 55, 535–543. [Google Scholar] [CrossRef]
  18. Größ, J. A note on the general Hermitian solution to AXA* = B. Bull. Malays. Math. Sci. Soc. 1998, 21, 57–62. [Google Scholar]
  19. Baksalary, J.K. Nonnegative definite and positive definite solutions to the matrix equation AXA* = B. Linear Multilinear Algebra 1984, 16, 133–139. [Google Scholar] [CrossRef]
  20. Größ, J. Nonnegative-definite and positive-definite solutions to the matrix equation AXA* = B revisited. Linear Algebra Appl. 2000, 321, 123–129. [Google Scholar] [CrossRef]
  21. Took, C.C.; Mandic, D.P. On the unitary diagonalisation of a special class of quaternion matrices. Appl. Math. Lett. 2011, 24, 1806–1809. [Google Scholar] [CrossRef]
  22. Took, C.C.; Mandic, D.P. Augmented second-order statistics of quaternion random signals. Signal Process. 2011, 91, 214–224. [Google Scholar] [CrossRef]
  23. Kyrchei, I. Cramer’s rules of η-(skew-)Hermitian solutions to the quaternion Sylvester-type matrix equations. Adv. Appl. Clifford Algebr. 2019, 29, 56. [Google Scholar] [CrossRef]
  24. Brambley, G.; Kim, J. Unit dual quaternion-based pose optimization for visual runway observations. IET Cyber Syst. Robot. 2020, 2, 181–189. [Google Scholar] [CrossRef]
  25. Cheng, J.; Kim, J.; Jiang, Z.; Che, W. Dual quaternion-based graph SLAM. Robot. Auton. Syst. 2016, 77, 15–24. [Google Scholar] [CrossRef]
  26. Wang, X.; Yu, C.; Lin, Z. A dual quaternion solution to attitude and position control for rigid body coordination. IEEE Trans. Rob. 2012, 28, 1162–1170. [Google Scholar] [CrossRef]
  27. Rodman, L. Topics in Quaternion Linear Algebra; Princeton University Press: Princeton, NJ, USA, 2014. [Google Scholar]
  28. Qi, L.Q.; Ling, C.; Yan, H. Dual quaternions and dual quaternion vectors. Commun. Appl. Math. Comput. 2022, 4, 1494–1508. [Google Scholar] [CrossRef]
  29. He, Z.H.; Wang, Q.W. The general solutions to some systems of matrix equations. Linear Multilinear Algebra 2015, 63, 2017–2032. [Google Scholar] [CrossRef]
  30. Marsaglia, G.; Styan, G.P. Equalities and inequalities for ranks of matrices. Linear Multilinear Algebra 1974, 2, 269–292. [Google Scholar] [CrossRef]
  31. Xie, L.M.; Wang, Q.W. A system of dual quaternion matrix equations with its applications. arXiv 2023, arXiv:2312.10037. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, Y.; Wang, Q.-W.; Xie, L.-M. Dual Quaternion Matrix Equation AXB = C with Applications. Symmetry 2024, 16, 287. https://doi.org/10.3390/sym16030287

AMA Style

Chen Y, Wang Q-W, Xie L-M. Dual Quaternion Matrix Equation AXB = C with Applications. Symmetry. 2024; 16(3):287. https://doi.org/10.3390/sym16030287

Chicago/Turabian Style

Chen, Yan, Qing-Wen Wang, and Lv-Ming Xie. 2024. "Dual Quaternion Matrix Equation AXB = C with Applications" Symmetry 16, no. 3: 287. https://doi.org/10.3390/sym16030287

APA Style

Chen, Y., Wang, Q. -W., & Xie, L. -M. (2024). Dual Quaternion Matrix Equation AXB = C with Applications. Symmetry, 16(3), 287. https://doi.org/10.3390/sym16030287

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