Abstract

We constitute some necessary and sufficient conditions for the system , , , , , to have a solution over the quaternion skew field in this paper. A novel expression of general solution to this system is also established when it has a solution. The least norm of the solution to this system is also researched in this article. Some former consequences can be regarded as particular cases of this article. Finally, we give determinantal representations (analogs of Cramer’s rule) of the least norm solution to the system using row-column noncommutative determinants. An algorithm and numerical examples are given to elaborate our results.

1. Introduction

In the whole article, the notation is reserved for the real number field and stands for the set of all matrices over the quaternion skew field specifies its subset of matrices with a rank . For , let and designate the conjugate transpose, the column right space and the left row space of . dim illustrates the size of and dim =dim by [1], which is known as the rank of denoted by .

Definition 1. The Moore-Penrose inverse of , denoted by , is defined to be the unique solution to the following four matrix equationsMatrices satisfying (1) and (2) are known as reflexive inverses.

Note that the reflexive inverse is denoted most often by but sometimes by (see, e.g., [2]) that is different from the denotation of the Moore-Penrose by . We will use the denotation for the reflexive inverse.

Suppose refers an identity matrix with feasible size. In addition, , represent a pair of orthogonal projectors induced by , respectively, and , , , , and .

Quaternions were invented by Hamilton in 1843. Zhang presented a detail survey on quaternion matrices in [3]. Quaternions provide a concise mathematical method for representing the automorphisms of three- and four-dimensional spaces. The representations by quaternions are more compact and quicker to compute than the representations by matrices [4]. For this reason, an increasing number of applications based on quaternions are found in various fields, such as color imaging, geometry, mechanics, linear adaptive filter, altitude control and computer science, signal processing, in particular as quaternion-valued neural networks, etc. [510].

The research of matrix equations have both applied and theoretical importance. In particular, the Sylvester-type matrix equations have far reaching applications in singular system control [11], system design [12], robust control [13], feedback [14], perturbation theory [15], linear descriptor systems [16], neural networks [17], and theory of orbits [18].

Some recent work on generalized Sylvester matrix equations and their systems can be observed in [1931]. In , Bao [32] examined the least-norm and extremal ranks of the least square solution to the quaternion matrix equationsWang et al. [33] examined the expression of the general solution to the systemAnd, as an application, the -symmetric and -skew-symmetric solution to has been established. Li et al. [34] established a novel expression to the general solution of system (4) and they computed the least-norm of general solution to (4). In , Wang et al. [35] constituted the expression of the general solution toand as an application they explored the -symmetric solution to the system

Some latest findings on the least-norm of matrix equations and -symmetric matrices can be consulted in [3640]. Furthermore, our main system (6) is a special case of the following system:which has been investigated by Zhang in 2014. But the expressions provided for the , and in [41], we are in position to calculate the least-norm of the solutions with its determinantal representations. When some given matrices are zero in (8), then it becomes our system and we will give such kind of expressions in which the least-norm of the solutions can also be computed with its determinantal representations. It is worthy to note that Zhang examined (8) with complex settings and we will consider our system (6) with quaternion settings.

According to our best of knowledge, the least-norm of the general solution to system (6) is not investigated by any one. Motivated by the vast application of quaternion matrices and the latest interest of least-norm of matrix equations, we construct a novel expression of the general solution to system (6) and apply this to investigate the least-norm of the general solution to system (6) over in this paper. Observing that systems (3) and (4) are particular cases of our system (6), solving system (6) will encourage the least-norm to a wide class of problems in the collected work.

Since the general solutions of considered systems are expressed in term of generalized inverses, another goal of the paper is to give determinantal representations of the least-norm of the general solution to system (6) based on determinantal representations of generalized inverses.

Determinantal representation of a solution gives a direct method of its finding analogous to the classical Cramer’s rule that has important theoretical and practical significance. Through looking for their more applicable explicit expressions, there are various determinantal representations of generalized inverses even with the complex or real entries, in particular for the Moore-Penrose inverse (see, e.g., [4244]). By virtue of noncommutativity of quaternions, the problem for determinantal representation of generalized quaternion inverses is more complicated, and only now it can be solved due to the theory of column-row determinants introduced in [45, 46]. Within the framework of the theory of noncommutative row-column determinants, determinantal representations of various generalized quaternion inverses and generalized inverse solutions to quaternion matrix equations have been derived by one of our authors (see, e.g.[4754]) and by other researchers (see, e.g. [5557]). Moreover, Song et al. [58] have just recently considered determinantal representations of general solution to the two-sided coupled generalized Sylvester matrix equation over obtained using the theory of row-column determinants as well. But their proposed approach differs from our proposed. In [58], for determinantal representations of the general solution to the equation supplementary matrices have been used that not always easy to get. While, by proposed method only coefficient matrices of the equations are used. More detailed Cramer’s rule to solutions and (skew-)Hermitian solutions of some systems of matrix equations and generalized Sylvester matrix equation over are recently explored in [59, 60] and [61, 62], respectively.

The remainder of our article is directed as follows. In Section 2, we commence with some needed known results about systems of matrix equations and determinantal representations of the Moore-Penrose inverse and of solutions to the quaternion matrix equations. In Section 3, we provide a new expression of the general solution to our system (6) and present an algorithm with an example. We discuss the least-norm of the general solution to (6) over in Section 4. In Section 5, determinantal representations of the general solution to (6) are derived and other example to elaborate obtained Cramer’s Rule to system (6) with data from the example in Section 3 is given. As expected, we get the same solution. Finally, in Section 6, the conclusions are drawn.

2. Preliminaries

We commence with the following lemmas which have crucial function in the construction of the chief outcomes of the following sections.

2.1. The General Solution to System (6)

Lemma 2 (see [63]). Let be given. Then
(1) .
(2)
(3)

Lemma 3 (see [64]). Let , , and be known matrices over with right sizes. Then
(1)
(2)
(3)

Lemma 4 (see [65]). Let be matrices over and Thenwhere , are any fixed reflexive inverses, and stand for the projectors , induced by , , respectively.

Remark 5. Since Moore-Penrose inverses are reflexive inverses, this lemma can be used for Moore-Penrose inverses without any changes. It has taken place in ([64], Lemma 2.4). But for more credibility, we prove this lemma below for the Moore-Penroses inverse as well.

Lemma 6 (see [66]). Suppose thatis consistent linear matrix equation, where , , , and , respectively. Then (1)The general solution of the homogeneous equation, can be expressed by where and are general solution of the following four homogeneous matrix expressions By computing the value of unknowns in the above equations and using them in and , we have where , , , , , and ; the matrices are free to vary over (2) Assume that the matrix expression (11) is solvable, then its general solution can be expressed as where and are any pair of particular solutions to (11).It can also be written as

Lemma 7 (see [67]). Let , , , and be given and to be determined. Then the systemis consistent if and only ifUnder these conditions, the general solution to (18) can be established as where is a free matrix over with accordant dimension.

Lemma 8 (see [35]). Let , , , , , , , ,, , , , , be given and , to be determined. Denote Then the following conditions are tantamount: (1)System (6) is resolvable.(2)The conditions in (19) are met and(3)The equalities in (19) and (22) are satisfied and In these conditions, the general solution to system (6) can be written aswhere are free matrices over with agreeable dimensions.

2.2. Determinantal Representations of Solutions to the Quaternion Matrix Equations

Due to noncommutativity of quaternions there is a problem of a determinant of matrices with noncommutative entries (which are also defined as noncommutative determinants). There are several versions of defining of noncommutative determinants (e.g., see [6870]). But any of the previous noncommutative determinants has not fully retained those properties which it owned for matrices with real settings. Moreover, if functional properties of a noncommutative determinant over a ring are satisfied, then it takes on a value in its commutative subset. This dilemma can be avoided due to the theory of row-column determinants.

For , we define row determinants and column determinants. Suppose is the symmetric group on the set .

Definition 9 (see [45]). The th row determinant of is defined for all by putting where is the left-ordered permutation. It means that its first cycle from the left starts with , other cycles start from the left with the minimal of all the integers which are contained in it,and the order of disjoint cycles (except for the first one) is strictly conditioned by increase from left to right of their first elements, .

Definition 10 (see [45]). The th column determinant of is defined for all by putting noindent where is the right-ordered permutation. It means that its first cycle from the right starts with , other cycles start from the right with the minimal of all the integers which are contained in it,and the order of disjoint cycles (except for the first one) is strictly conditioned by increase from right to left of their first elements, .
Since [45] for Hermitian we havethe determinant of a Hermitian matrix is defined by puttingIts properties are similar to the properties of an usual (commutative) determinant and they have been completely explored in [46] by using row and column determinants that are so defined only by construction.

For determinantal representations of the Moore-Penrose inverse, we shall use the following notations. Let and be subsets of the order . Let be a submatrix of whose rows are indexed by and the columns indexed by . Similarly, let be a principal submatrix of whose rows and columns indexed by . If is Hermitian, then is the corresponding principal minor of . For , the collection of strictly increasing sequences of integers chosen from is denoted by . For fixed and , let denotes the collection of sequences of row indexes that contain the index , and denotes the collection of sequences of column indexes that contain .

Let be the th column and be the th row of , respectively. Suppose denotes the matrix obtained from by replacing its th column with the column-vector , and denotes the matrix obtained from by replacing its th row with the row-vector . We denote the th row and the th column of by and , respectively.

Lemma 11 (see [47]). If , then the Moore-Penrose inverse have the following determinantal representations,and

Remark 12. For an arbitrary full-rank matrix , a row-vector , and a column-vector , we put

Corollary 13. If , then the projection matrix has the determinantal representation where is the th column of .

Corollary 14. If , then the projection matrix has the determinantal representationwhere is the th row of .

Lemma 15 (see [2]). Let , , be known and be unknown. Then the matrix equationis consistent if and only if . In this case, its general solution can be expressed aswhere are arbitrary matrices over with appropriate dimensions.

In [71], it’s proved that (39) is the least squares solution to (38), and its minimum norm least squares solution is .

Lemma 16 (see [48]). Let , . Then the minimum norm least squares solution to (38) have determinantal representations, or where are the column vector and the row vector, respectively. and are the th row and the th column of .

Corollary 17. Let , be known and be unknown. Then the matrix equation is consistent if and only if . In this case, its general solution can be expressed as , where is an arbitrary matrix over with appropriate dimensions. Its minimum norm least squares solution has the following determinantal representation, where is the th column of .

Corollary 18. Let , be given, and be unknown. Then the equation is solvable if and only if and its general solution is , where is a any matrix with conformable dimension. Moreover, its minimum norm least squares solution has the determinantal representation, where is the th row of .

3. A New Expression of the General Solution to System (6)

First, we show that Lemma 4 is true for the Moore-Penrose inverses.

Lemma 19. Let be matrices over and Thenwhere , are the Moore-Penrose inverses, and , , , , , and are projectors with respect to the corresponding Moore-Penrose inverses.

Proof. In ([65], Lemma 2.4), it is proved that for fixed reflexive inverses and , the reflexive inverse can be expressed as follows, We choose , as the Moore-Penrose inverses, and as the projector with respect to the Moore-Penrose inverse and show that the obtained matrixis the Moore-Penrose inverse of . For this, it is enough to proof that satisfies the conditions (3) and (4) in Definition 1.
Since by Lemma 3, , then can be expressed as So, Since condition (3) is satisfied by components, namely, it follows that satisfies condition (3) as well; i.e., .
Similar, it can be shown that satisfies condition (4). Hence, the Moore-Penrose inverse of can be expressed by (49). From this (46) immediately follow.
The equations (47) can be proved similarly.

Now we demonstrate the principal theorem of this section.

Theorem 20. Assume that , , , , , , , , , , , and system (6) is solvable, then the general solution to our system can be formed aswhere are free matrices over with allowable dimensions.

Proof. Our proof contains three parts. At the first step, we show that the matrices and have the forms ofwhere and are any pair of particular solution to system (6), , , , and are free matrices of able shapes over , are solutions to system (6). At the second step, we display that any couple of solutions and to system (6) can be established as (56) and (57), respectively. At the end, we confirm that and are a couple of particular solutions to system (6).
Now we prove that a couple of matrices and having the shape of (56) and (57), respectively, are solutions to system (6). Observe that It is evident that having the form (56) is a solution of , and and having the form (57) is a solution to . Now we are left to show that is satisfied by and given in (56) and (57). By Lemma 4, we haveandObserve that and by using (61) and (62), we arrive that Conversely, assume that and are any couple of solutions to our system (6). By Lemma 7, we have Observe that producesOn the same lines, we can getIt is manifest that and defined in (66)-(67) are also solution pair ofSince Hence by Lemma 6, and can be written aswhere and are a couple of special solutions to (68) and and are free matrices with agreeable dimensions. Using (70) and (71) in (66) and (67), respectively, we get where and It is evident that and are a couple of solutions to system (6). It is clear that and can be represented by (56) and (57), respectively. Lastly, by putting , and equal to zero in (24) and (25), we conclude that and are special solutions to system (6). Hence the expressions (54) and (55) represent the general solution to system (6) and the theorem is completed.

Remark 21. Due to Lemma 3 and taking into account , we have the following simplification of the solution pair to system (6) that is identical for (24)-(25) and (54)-(55) when , and disappear,

Comment 1. We have established a novel expression of the general solution to system (6) in Theorem 20 which is different from one created in [35]. With the help of this novel expression, we can explore the least-norm of the general solution which can not be studied with the help of the expression given in [35], which is one of the advantage of our new expression.

Now we discuss some special cases of our system.

If and disappear in Theorem 20, then we gain the following conclusion.

Corollary 22. Denote , , , , , and system (4) is solvable, then the general solution to system (4) can be formed as where are the same as in Lemma 6, , , and are free matrices over obeying agreeable dimensions.

Comment 2. The above consequence is a chief result of [64].
If and vanish in our system (6), then we get the following outcome.

Corollary 23. Suppose that and are given. Then the general solution to system (3) is established by where and are arbitrary matrices over with appropriate sizes.

Comment 3. Corollary 23 is the rudimentary result of [32].

Comment 4. When and become zero in (8), then we will get the least-norm of the solution of (8) with the help of Theorem 20 quite smoothly. This is one of the advantage of the our expressions over the expressions given in [41].

An algorithm and numerical example is provided to obtain the general solution of (6) with the help of Theorem 20.

Algorithm 24. (1) Input , , , , , , , , , , , , , with viable dimensions over .
(2) Evaluate and by (54)-(55).

Example 25. For given matrices By these given matrices, the consistency conditions of (6) from Lemma 3 are fulfilled. So, system (6) is resolvable. Now we compute the partial solution to system (6) when , and disappear. Using determinantal representations (33)-(34) for computing Moore-Penrose inverses, we find that Since and , then are zero-matrices. Hence the general solution to our system (6) is

4. The Least-Norm of the General Solution to System (6)

We experience the least-norm to system (6) in this section. We first modify the description of quaternionic inner product space defined in [72] as follows:

A right -vector space is a quaternionic inner product space if there is a mapping such that for all and :(1); ;(2);(3)

It can be achieved by putting for . is referred as the norm of . It is routine to verify that is a quaternionic inner product space under the inner product defined by where . is the matrix norm defined by . The real part of a quaternion is denoted by .

By the definition and [73], we can get the following result easily.

Lemma 26. Let . Then we have
(1) .
(2)

Theorem 27. Assume that system (6) is solvable, then the least-norm of the solution pair and to system (6) can be extracted as follows:

Proof. With the help of Theorem 20 and Remark 21, the general solution to system (6) can be formed as where , and are free matrices over having executable dimensions. By Lemma 26, the norm of can be established aswhereNow we want to show that . Applying Lemmas 3, 4, and 26, we haveBy using (84)-(86) in (83) produces . Since is arbitrary, we get (79) from (82). On the same way, we can prove that (80) hold.

A special cases of our system (6) are given below.

If and become zero matrices in Theorem 27, then again we get the principal result of [30].

Corollary 28. Assume that system (4) is solvable, then the least-norm of the solution pair and to system (4) can be furnished as

If , , , , , and vanish in our system, then we get the next consequence.

Corollary 29. Suppose that , and are given. Then the least-norm of the least square solution to system (3) is launched by

Comment 5. Corollary 29 is the key result of [32].

5. Determinantal Representations of the Least-Norm Solution to System (6)

In this section, we give determinantal representations of the least-norm solution to system (6). Let , , , , , , , , , , , , , and .

First, consider each term of (79) separately.

(i) Denote . Due to Corollary 17 for the first term of (79), , we havewhere is the th column of .

(ii) For the second term of (79) we have, . So, due to Corollaries 18 and 13,where is the th row of and is the th column of .

Construct the matrix , where

and denote . Then, from (90), it follows that where is the th row of the matrix .

If we construct the matrix , where and denote , and then, from (90), we obtainwhere is the th column of the matrix .

(iii) Due to Theorem 2.15 for the third term , we obtainorwhereare the column vector and the row vector, respectively. and are the th row and the th column of .

(iv) For the fourth term of (79), , using the determinantal representation (33) for and by Theorem 2.15, we havewhereand are the column vector and the row vector, respectively. is the th column of , , and are the th row and the th column of , respectively. Construct the matrix , where is given by (99) and denote . Then, from (98), we get the following final determinantal representation of the fourth term of (79),where is the th column of .

(v) For the fifth term of (79), , due to Corollary 17 to , by Theorem 2.15 to , and Corollary 18 to , we obtain where is the th column of , is the th row of ,andare the column vector and the row vector, respectively. and are the th row and the th column of . Construct the matrix , where is determined by (103), and denote . Then, from (102), it follows that where and are the unit row-vector and the unit column-vector whose components are 0 except the th or th components which are 1, respectively.

If we denotethen, from (102), it follows the determinantal representationwhere is the th row of the matrix that is determined by (106).

If we denotethen, from (102), it follows the determinantal representationwhere is the th column of the matrix that is determined by (108).

Similarly, consider each term of (80) separately.

(i) Denote . Due to Corollary 17 for the first term of (80), , we havewhere is the th column of .

(ii) For the second term of (80) we have, . So, due to Corollaries 18 and 13,where is the th row of and is the th column of .

Construct the matrix , whereand denote . Then, from (111), it follows thatwhere is the th row of .

If we construct the matrix , whereand denote , then, from (111), we obtainwhere is the th row of .

(iii) Due to Theorem 2.15 for the third term , we obtain

orwhere are the column vector and the row vector, respectively. and are the th row and the th column of .

(iv) Using Corollary 13 to and by Theorem 2.15 to , we obtain the the following representation of the fourth term, , of (80)where is determined by (103). Construct the matrix and denote . Then, from (119) finally, we havewhere is the th column of .

So, we prove the following theorem.

Theorem 30. Let , , , , , , , , , , , , , and . The least-norm solution (79)-(80) to system (8), , , by componentshas determinantal representations, where the term is (89), is (92) or (94), is (95) or (96), is (101), and is (107) or (109); similarly, is (110), is (113) or (115), is (116) or (117), and is (120).

A numerical example is provided to obtain the least norm of the general solution of (6) with the help of Theorem 30.

Example 31. We use the given matrices from the Example 25. Since andand then, by (89), Now, by (92), we find for all . So,Similarly, by (91), . SoSince , then by (92), Since andand then, by (95), and Further, due to Example 25, for all . So, Since, andand, due to Example 25, for all ; then by (110) and (121) Hence, the least norm solution of (6) obtained by Cramer’s Rule and the matrix method in Example 25 are the same as expected.
Note that we used Maple with the package CLIFFORD in the calculations.

6. Conclusion

We have constructed a novel expression of the general solution to system (6) over and used this result to explore the least-norm of the general solution to this system when it is solvable. Some particular cases of our system are also discussed. Our results carry the principal results of [32, 64]. Finally, we give determinantal representations (analogous of Cramer’s Rule) of the least norm solutions to the systems using row-column noncommutative determinants. Numerical examples are also provided to interpret the results.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.