Table of Contents Author Guidelines Submit a Manuscript
Journal of Applied Mathematics
Volume 2014, Article ID 857081, 9 pages
http://dx.doi.org/10.1155/2014/857081
Research Article

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, [235]. 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

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

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

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

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

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

  1. 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
  2. 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
  3. 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
  4. M. Dehghan and M. Hajarian, “Finite iterative algorithms for the reflexive and anti-reflexive solutions of the matrix equation A1X1B1+A2X2B2=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
  5. J. D. Gardiner, A. J. Laub, J. J. Amato, and C. B. Moler, “Solution of the Sylvester matrix equation AXBT+CXDT=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
  6. 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
  7. 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
  8. N. Li, Q.-W. Wang, and J. Jiang, “An efficient algorithm for the reflexive solution of the quaternion matrix equation AXB+CXHD=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
  9. 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
  10. 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
  11. 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
  12. S. K. Mitra, “Common solutions to a pair of linear matrix equations A1XB1=C1 and A2XB2=C2,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 74, pp. 213–216, 1973. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  13. A. Navarra, P. L. Odell, and D. M. Young, “A representation of the general common solution to the matrix equations A1XB1=C1 and A2XB2=C2 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
  14. 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 A1XB1=C1, A2XB2=C2,” 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
  15. Z.-H. Peng, X.-Y. Hu, and L. Zhang, “An efficient algorithm for the least-squares reflexive solution of the matrix equation A1XB1=C1, A2XB2=C2,” 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
  16. 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
  17. J. W. van der Woude, “On the existence of a common solution X to the matrix equations AiXBj=Cij,(i,j)Γ,” 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
  18. 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
  19. 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
  20. 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
  21. 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
  22. 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
  23. 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
  24. 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
  25. 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
  26. 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
  27. 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
  28. 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
  29. 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
  30. 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
  31. 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
  32. 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
  33. 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
  34. 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
  35. 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
  36. 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
  37. S. F. Yuan and A. P. Liao, “Least squares solution of the quaternion matrix equation X-AX^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
  38. 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
  39. S. K. Mitra, “The matrix equation AXB + 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
  40. 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
  41. 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
  42. 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
  43. A. Ben-Israel and T. N. E. Greville, Generalized Inverses: Theory and Applications, John Wiley & Sons, New York, NY, USA, 1974. View at MathSciNet