Journal of Applied Mathematics

Volume 2014 (2014), Article ID 857081, 9 pages

http://dx.doi.org/10.1155/2014/857081

## Least Squares Pure Imaginary Solution and Real Solution of the Quaternion Matrix Equation with the Least Norm

School of Mathematics and Computational Science, Wuyi University, Jiangmen 529020, China

Received 5 December 2013; Accepted 3 February 2014; Published 15 April 2014

Academic Editor: Qing-Wen Wang

Copyright © 2014 Shi-Fang Yuan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Using the Kronecker product of matrices, the Moore-Penrose generalized inverse, and the complex representation of quaternion matrices, we derive the expressions of least squares solution with the least norm, least squares pure imaginary solution with the least norm, and least squares real solution with the least norm of the quaternion matrix equation , respectively.

#### 1. Introduction

Quaternions were introduced by Irish mathematician Sir William Rowan Hamilton in 1843. The family of quaternions is a skew field or noncommutative division algebra, since the characteristic property of quaternions is their noncommutativity under multiplication. A quaternion can be uniquely expressed as with real coefficients , , , , , , and can be uniquely expressed as , where and are complex numbers. Thus, every pure imaginary quaternion can be uniquely expressed as .

Throughout this paper, let , , , , and be the skew field of quaternions, the set of all real matrices, the set of all complex matrices, the set of all quaternion matrices, and the set of all pure imaginary quaternion matrices, respectively. For , and denote the real part and the imaginary part of matrix , respectively. For , , , , and denote the conjugate matrix, the transpose matrix, the conjugate transpose matrix, and the Moore-Penrose generalized inverse matrix of matrix , respectively.

For any , can be uniquely expressed as , where . The complex representation matrix of is denoted by Notice that is uniquely determined by . For , , we have (see [1]). Denote the trace of a square matrix by . We define the inner product for all . Then is a Hilbert inner product space and the norm of a matrix generated by this inner product is the quaternion matrix Frobenius norm . The 2-norm of the vector is denoted by .

Various aspects of the solutions of matrix equations such as , , have been investigated. See, for example, [2–35]. For the matrix equation if and are identity matrices, then the matrix equation (2) reduces to the well-known Sylvester equation [5, 36]. If and are identity matrices, then the matrix equation (2) reduces to the well-known Stein equation [37]. There are many important results about the matrix equation (2). For example, Hernández and Gassó [38] obtained the explicit solution of the matrix equation (2). Mansour [11] considered the solvability condition of the matrix equation (2) in the operator algebra. Mitra [39] studied the solvability conditions of matrix equation (2). For the quaternion matrix equation (2). Huang [40] obtained necessary and sufficient conditions for the existence of a solution or a unique solution using the method of complex representation of quaternion matrices.

Note that some authors have investigated the real and pure imaginary solutions to the quaternion matrix equations. For example, Au-Yeung and Cheng [2] considered the pure imaginary quaternionic solutions of the Hurwitz matrix equations. Wang et al. [41] studied the quaternion matrix equation and obtained necessary and sufficient conditions for the existence of a real solution or pure imaginary solution of the quaternion matrix equation . Using the complex representation of quaternion matrices and the Moore-Penrose generalized inverse, Yuan et al. [42] derived the expressions of the least squares solution with the least norm, the least squares pure imaginary solution with the least norm, and the least squares real solution with the least norm for the quaternion matrix equation , respectively. Motivated by the work mentioned above, in this paper, we will consider the related problem of quaternion matrix equation (2).

*Problem 1. *Given , , , , , let

Find such that

*Problem 2. *Given , , , , , let
Find such that

*Problem 3. *Given , , , , , let
Find such that

The solution of Problem 1 is called the least squares solution with the least norm; the solution of Problem 2 is called the least squares pure imaginary solution with the least norm; and the solution of Problem 3 is called the least squares real solution with the least norm for matrix equation (2) over the skew field of quaternions.

This paper is organized as follows. In Section 2, we derive the explicit expression for the solution of Problem 1. In Section 3, we derive the explicit expression for the solution of Problem 2. In Section 4, we derive the explicit expression for the solution of Problem 3. Finally, in Section 5, we report numerical algorithms and numerical examples to illustrate our results.

#### 2. The Solution of Problem 1

To study Problem 1, we begin with the following lemmas.

Lemma 4 (see [43]). *The matrix equation , with and , has a solution if and only if
**
in this case it has the general solution
**
where is an arbitrary vector.*

Lemma 5 (see [43]). *The least squares solutions of the matrix equation , with and , can be represented as
**
where is an arbitrary vector, and the least squares solution of the matrix equation with the least norm is .*

We identify with a complex vector and denote such an identification by the symbol *≅*, that is;
For , we have and
We denote ,
Notice that . In particular, for with , we have , and

Addition of two quaternion matrices and is defined by whereas multiplication is defined as So , ; moreover, can be expressed as

Lemma 6 (see [34]). *Let , , and be given. Then
*

*Lemma 7. For , let
Then
*

*Proof. *For , we have

*By Lemmas 6 and 7, we have the following.*

*Lemma 8. If , , and , then
*

*Based on our earlier discussions, we now turn our attention to Problem 1. The following notations are necessary for deriving the solutions of Problem 1. For , , , , , set
From the results in [33], one has
*

*Theorem 9. Let , , , , and , and let , , be as in (24). Then
where is an arbitrary vector of appropriate order.*

*Proof. *By Lemmas 5 and 8, we can get
By Lemma 5, it follows that
thus
The proof is completed.

*By Lemma 4 and Theorem 9, we get the following conclusion.*

*Corollary 10. The quaternion matrix equation (2) has a solution if and only if
In this case, denote by the solution set of (2). Then
where is an arbitrary vector of appropriate order.*

Furthermore, if (31) holds, then the quaternion matrix equation (2) has a unique solution if and only if In this case,

*Theorem 11. Problem 1 has a unique solution . This solution satisfies
*

*Proof. *From (27), it is easy to verify that the solution set is nonempty and is a closed convex set. Hence, Problem 1 has a unique solution .

We now prove that the solution can be expressed as (35).

From (27), we have
by Lemma 5 and (27),
Thus,
Thus we have completed the proof.

*Corollary 12. The least norm problem
has a unique solution and can be expressed as (35).*

*3. The Solution of Problem 2*

*3. The Solution of Problem 2*

*We now discuss the solution of Problem 2. For , we have .*

*For , , , , , set
Thus we have
*

*We now study Problem 2. Since the methods are the same as in Section 2, we only describe the following results using Lemmas 4 and 5 and Theorem 9 and omit their detailed proofs.*

*Theorem 13. Let , , , , and ; let , , be as in (41). Then the set of Problem 2 can be expressed as
where is an arbitrary vector of appropriate order.*

*Corollary 14. The quaternion matrix equation (2) has a solution if and only if
In this case, denote by the pure imaginary solution set of (2). Then
where is an arbitrary vector of appropriate order.*

Furthermore, if (45) holds, then the quaternion matrix equation (2) has a unique solution if and only if In this case,

*Theorem 15. Problem 2 has a unique solution . This solution satisfies
*

*Corollary 16. The least norm problem
has a unique solution and can be expressed as (49).*

*4. The Solution of Problem 3*

*4. The Solution of Problem 3*

*We now discuss the solution of Problem 3. For , we have , and . Thus, we have the following lemmas.*

*For , , , , , set
We have
*

*By Lemma 5, we can easily get the following results for Problem 3.*

*Theorem 17. Let , , , , and ; let , , be as in (52). Then the set of Problem 3 can be expressed as
where is an arbitrary vector of appropriate order.*

*Corollary 18. The quaternion matrix equation (2) has a solution if and only if
In this case, denote by the real solution set of (2). Then
where is an arbitrary vector of appropriate order.*

Furthermore, if (56) holds, then the quaternion matrix equation (2) has a unique solution if and only if In this case,

*Theorem 19. Problem 3 has a unique solution . This solution satisfies
*

*Corollary 20. The least norm problem
has a unique solution and can be expressed as (60).*

*5. Numerical Verification*

*5. Numerical Verification*

*Based on the discussions in Sections 2, 3, and 4, we report numerical tests in this section. We give three numerical algorithms and four numerical examples to find the solutions of Problems 1, 2, and 3.*

*Algorithms 21, 22, and 23 provide the methods to find the solutions of Problems 1, 2, and 3. If the consistent conditions for matrix equation (2) hold, Examples 24 and 25 consider the numerical solutions of Problem 1 for . In Examples 26 and 27, if the consistent conditions for matrix equation (2) are not satisfied, we can compute the least squares solution with the least norm in Problems 2 and 3 by Algorithms 22 and 23, respectively. For demonstration purpose and avoiding the matrices with large norm to interrupt the solutions of Problems 1 and 2, we only consider the coefficient matrices of small sizes in numerical experiments.*

*Algorithm 21 (for Problem 1). *We have the following.(1)Input , , , , and (, , , , and ).(2)Compute , , , , , , , , .(3)If (31) and (33) hold, then calculate according to (34).(4)If (31) holds, then calculate according to (35). Otherwise go to next step.(5)Calculate according to (35).

*Algorithm 22 (for Problem 2). *We have the following.(1)Input , , , , and (, , , , and ).(2)Compute , , , , , , , , .(3)If (45) and (47) hold, then calculate according to (48).(4)If (45) holds, then calculate according to (49). Otherwise go to next step.(5)Calculate according to (49).

*Algorithm 23 (for Problem 3). *We have the following.(1)Input , , , , and (, , , , and ).(2)Compute , , , , , , , , .(3)If (56) and (58) hold, then calculate according to (59).(4)If (56) holds, then calculate according to (60). Otherwise go to next step.(5)Calculate according to (60).

*Example 24. *Let , , , ,
, where

Let
By using matlab 7.7 and Algorithm 21, we obtain
According to Algorithm 21 , we can see the matrix equation has a unique solution which is a unique solution with the least norm . We can get .

*Example 25. *Let , , , ,
, where

Let
By using matlab 7.7 and Algorithm 21, we obtain
According to Algorithm 21 , we can see the matrix equation has infinite solution and a unique solution with the least norm , and we can get .

*Example 26. *Suppose , , , , , , , , , are the same as in Example 25; by using matlab 7.7 and Algorithm 22, we obtain
According to Algorithm 22 , we can see the matrix equation has infinite pure imaginary least squares solutions and a unique pure imaginary solution with the least norm for Problem 2 and we can get , and , where

*Example 27. *Suppose , , , , , , , , , are the same as in Example 25. By using matlab 7.7 and Algorithm 22, we obtain
According to Algorithm 23 , we can see the matrix equation has infinite least squares real solutions and a unique least squares real solution with the least norm for Problem 3 and we can get , and
Examples 24, 25, 26, and 27 are used to show the feasibility of Algorithms 21, 22, and 23.

*Conflict of Interests*

*Conflict of Interests*

*The author declares that there is no conflict of interests regarding the publication of this paper.*

*Acknowledgments*

*Acknowledgments*

*This work is supported by Natural Science Foundation of China (no. 11301397), Guangdong Natural Science Fund of China (no. 10452902001005845), and Science and Technology Project of Jiangmen City, China.*

*References*

*References*

- F. Z. Zhang, “Quaternions and matrices of quaternions,”
*Linear Algebra and Its Applications*, vol. 251, pp. 21–57, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y.-H. Au-Yeung and C.-M. Cheng, “On the pure imaginary quaternionic solutions of the Hurwitz matrix equations,”
*Linear Algebra and Its Applications*, vol. 419, no. 2-3, pp. 630–642, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K.-W. E. Chu, “Singular value and generalized singular value decompositions and the solution of linear matrix equations,”
*Linear Algebra and Its Applications*, vol. 88-89, pp. 83–98, 1987. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Dehghan and M. Hajarian, “Finite iterative algorithms for the reflexive and anti-reflexive solutions of the matrix equation ${A}_{1}{X}_{1}{B}_{1}+{A}_{2}{X}_{2}{B}_{2}=C$,”
*Mathematical and Computer Modelling*, vol. 49, no. 9-10, pp. 1937–1959, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. D. Gardiner, A. J. Laub, J. J. Amato, and C. B. Moler, “Solution of the Sylvester matrix equation $AX{B}^{T}+CX{D}^{T}=E$,”
*ACM Transactions on Mathematical Software*, vol. 18, no. 2, pp. 223–231, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z.-H. He and Q.-W. Wang, “A real quaternion matrix equation with applications,”
*Linear and Multilinear Algebra*, vol. 61, no. 6, pp. 725–740, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. Li and Q.-W. Wang, “Iterative algorithm for solving a class of quaternion matrix equation over the generalized $(P,Q)$-reflexive matrices,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 831656, 15 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - N. Li, Q.-W. Wang, and J. Jiang, “An efficient algorithm for the reflexive solution of the quaternion matrix equation $AXB+C{X}^{H}D=F$,”
*Journal of Applied Mathematics*, vol. 2013, Article ID 217540, 14 pages, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y.-T. Li and W.-J. Wu, “Symmetric and skew-antisymmetric solutions to systems of real quaternion matrix equations,”
*Computers & Mathematics with Applications*, vol. 55, no. 6, pp. 1142–1147, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A.-P. Liao and Y. Lei, “Least-squares solution with the minimum-norm for the matrix equation $(AXB,GXH)=(C,D)$,”
*Computers & Mathematics with Applications*, vol. 50, no. 3-4, pp. 539–549, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Mansour, “Solvability of $AXB-CXD=E$ in the operators algebra $B(H)$,”
*Lobachevskii Journal of Mathematics*, vol. 31, no. 3, pp. 257–261, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. K. Mitra, “Common solutions to a pair of linear matrix equations ${A}_{1}X{B}_{1}={C}_{1}$ and ${A}_{2}X{B}_{2}={C}_{2}$,”
*Mathematical Proceedings of the Cambridge Philosophical Society*, vol. 74, pp. 213–216, 1973. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Navarra, P. L. Odell, and D. M. Young, “A representation of the general common solution to the matrix equations ${A}_{1}X{B}_{1}={C}_{1}$ and ${A}_{2}X{B}_{2}={C}_{2}$ with applications,”
*Computers & Mathematics with Applications*, vol. 41, no. 7-8, pp. 929–935, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y.-X. Peng, X.-Y. Hu, and L. Zhang, “An iterative method for symmetric solutions and optimal approximation solution of the system of matrix equations ${A}_{1}X{B}_{1}={C}_{1}$, ${A}_{2}X{B}_{2}={C}_{2}$,”
*Applied Mathematics and Computation*, vol. 183, no. 2, pp. 1127–1137, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z.-H. Peng, X.-Y. Hu, and L. Zhang, “An efficient algorithm for the least-squares reflexive solution of the matrix equation ${A}_{1}X{B}_{1}={C}_{1}$, ${A}_{2}X{B}_{2}={C}_{2}$,”
*Applied Mathematics and Computation*, vol. 181, no. 2, pp. 988–999, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z.-Y. Peng, “Solutions of symmetry-constrained least-squares problems,”
*Numerical Linear Algebra with Applications*, vol. 15, no. 4, pp. 373–389, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. W. van der Woude, “On the existence of a common solution $X$ to the matrix equations ${A}_{i}X{B}_{j}={C}_{ij},(i,j)\in \Gamma $,”
*Linear Algebra and Its Applications*, vol. 375, pp. 135–145, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. X. Yuan, “Two classes of best approximation problems of matrices,”
*Mathematica Numerica Sinica*, vol. 23, no. 4, pp. 429–436, 2001. View at Google Scholar · View at MathSciNet - Y. X. Yuan, “The minimum norm solutions of two classes of matrix equations,”
*Numerical Mathematics*, vol. 24, no. 2, pp. 127–134, 2002. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. X. Yuan, “The optimal solution of linear matrix equations by using matrix decompositions,”
*Mathematica Numerica Sinica*, vol. 24, no. 2, pp. 165–176, 2002 (Chinese). View at Google Scholar · View at MathSciNet - T. S. Jiang, Y. H. Liu, and M. S. Wei, “Quaternion generalized singular value decomposition and its applications,”
*Applied Mathematics*, vol. 21, no. 1, pp. 113–118, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. S. Jiang and M. S. Wei, “Real representations of quaternion matrices and quaternion matrix equations,”
*Acta Mathematica Scientia A*, vol. 26, no. 4, pp. 578–584, 2006 (Chinese). View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q.-W. Wang, “The general solution to a system of real quaternion matrix equations,”
*Computers & Mathematics with Applications*, vol. 49, no. 5-6, pp. 665–675, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q.-W. Wang, H.-X. Chang, and C.-Y. Lin, “$P$-(skew)symmetric common solutions to a pair of quaternion matrix equations,”
*Applied Mathematics and Computation*, vol. 195, no. 2, pp. 721–732, 2008. View at Publisher · View at Google Scholar · View at MathSciNet - Q.-W. Wang and Z.-H. He, “Solvability conditions and general solution for mixed Sylvester equations,”
*Automatica*, vol. 49, no. 9, pp. 2713–2719, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - Q. Wang and Z. He, “A system of matrix equations and its applications,”
*Science China Mathematics*, vol. 56, no. 9, pp. 1795–1820, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q.-W. Wang, J. W. van der Woude, and H.-X. Chang, “A system of real quaternion matrix equations with applications,”
*Linear Algebra and Its Applications*, vol. 431, no. 12, pp. 2291–2303, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q. Wang, J. W. van der Woude, and S. W. Yu, “An equivalence canonical form of a matrix triplet over an arbitrary division ring with applications,”
*Science China Mathematics*, vol. 54, no. 5, pp. 907–924, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q.-W. Wang, S.-W. Yu, and C.-Y. Lin, “Extreme ranks of a linear quaternion matrix expression subject to triple quaternion matrix equations with applications,”
*Applied Mathematics and Computation*, vol. 195, no. 2, pp. 733–744, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q.-W. Wang and F. Zhang, “The reflexive re-nonnegative definite solution to a quaternion matrix equation,”
*Electronic Journal of Linear Algebra*, vol. 17, pp. 88–101, 2008. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q. W. Wang, H.-S. Zhang, and G.-J. Song, “A new solvable condition for a pair of generalized Sylvester equations,”
*Electronic Journal of Linear Algebra*, vol. 18, pp. 289–301, 2009. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q.-W. Wang, X. Zhang, and J. W. van der Woude, “A new simultaneous decomposition of a matrix quaternity over an arbitrary division ring with applications,”
*Communications in Algebra*, vol. 40, no. 7, pp. 2309–2342, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. R. Magnus, “$L$-structured matrices and linear matrix equations,”
*Linear and Multilinear Algebra*, vol. 14, no. 1, pp. 67–88, 1983. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. F. Yuan, A. P. Liao, and Y. Lei, “Least squares Hermitian solution of the matrix equation $(AXB,CXD)=(E,F)$ with the least norm over the skew field of quaternions,”
*Mathematical and Computer Modelling*, vol. 48, no. 1-2, pp. 91–100, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S.-F. Yuan, Q.-W. Wang, and X. Zhang, “Least-squares problem for the quaternion matrix equation $AXB+CYD=E$ over different constrained matrices,”
*International Journal of Computer Mathematics*, vol. 90, no. 3, pp. 565–576, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A.-G. Wu, F. Zhu, G.-R. Duan, and Y. Zhang, “Solving the generalized Sylvester matrix equation $AV+BW=EVF$ via a Kronecker map,”
*Applied Mathematics Letters*, vol. 21, no. 10, pp. 1069–1073, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. F. Yuan and A. P. Liao, “Least squares solution of the quaternion matrix equation $X-A\widehat{X}B=C$ with the least norm,”
*Linear and Multilinear Algebra*, vol. 59, no. 9, pp. 985–998, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. Hernández and M. Gassó, “Explicit solution of the matrix equation $AXB-CXD=E$,”
*Linear Algebra and Its Applications*, vol. 121, pp. 333–344, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. K. Mitra, “The matrix equation $\mathrm{AXB}+\mathrm{CXD}=E$,”
*SIAM Journal on Applied Mathematics*, vol. 32, no. 4, pp. 823–825, 1977. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. P. Huang, “The matrix equation $AXB-GXD=E$ over the quaternion field,”
*Linear Algebra and Its Applications*, vol. 234, pp. 197–208, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q. W. Wang, S. W. Yu, and W. Xie, “Extreme ranks of real matrices in solution of the quaternion matrix equation $AXB=C$ with applications,”
*Algebra Colloquium*, vol. 17, no. 2, pp. 345–360, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S.-F. Yuan, Q.-W. Wang, and X.-F. Duan, “On solutions of the quaternion matrix equation $AX=B$ and their applications in color image restoration,”
*Applied Mathematics and Computation*, vol. 221, pp. 10–20, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - A. Ben-Israel and T. N. E. Greville,
*Generalized Inverses: Theory and Applications*, John Wiley & Sons, New York, NY, USA, 1974. View at MathSciNet

*
*