Elsevier

Applied Mathematics and Computation

Volume 265, 15 August 2015, Pages 945-957
Applied Mathematics and Computation

Solution to a system of real quaternion matrix equations encompassing η-Hermicity

https://doi.org/10.1016/j.amc.2015.05.104Get rights and content

Abstract

Let Hm×n be the set of all m × n matrices over the real quaternion algebra H={c0+c1i+c2j+c3ki2=j2=k2=ijk=1,c0,c1,c2,c3R}. AHn×n is known to be η-Hermitian if A=Aη*=ηA*η,η{i,j,k} and A* means the conjugate transpose of A. We mention some necessary and sufficient conditions for the existence of the solution to the system of real quaternion matrix equations including η-Hermicity A1X=C1,A2Y=C2,YB2=D2,Y=Yη*,A3Z=C3,ZB3=D3,Z=Zη*,A4X+(A4X)η*+B4YB4η*+C4ZC4η*=D4,and also construct the general solution to the system when it is consistent. The outcome of this paper diversifies some particular results in the literature. Furthermore, we constitute an algorithm and a numerical example to comprehend the approach established in this treatise.

Introduction

In the whole exposition, we symbolize the real number field by R, the complex number field by C and H for the quaternion algebra H={b0+b1i+b2j+b3ki2=j2=k2=ijk=1,b0,b1,b2,b3R}.The set of all matrices of dimension m × n over H is represented by Hm×n. I denotes an identity matrix with right size. For any matrix A over H, the column right space and the row left space of A are denoted by R(A) and N(A), respectively. D[R(A)] denotes the dimension of R(A). By Hungerford [30], we have D[R(A)]=D[N(A)], which is known as rank of A. A* means the conjugate transpose of A. A means the Moore–Penrose inverse of AHm×n, i.e., the exclusive matrix YHn×m obeying AYA=A,YAY=Y,(AY)*=AY,(YA)*=YA.For significant findings on generalized inverses, one can see [1] and [7]. Furthermore, LA=IAA and RA=IAA denote the projectors generated by A, respectively. By the definition of the Moore–Penrose inverse, these projectors are Hermitian and idempotent as well.

The quaternions were first explored by the Irish mathematician Sir William Rowan Hamilton in [2]. Quaternions have massive applications in diverse areas of mathematics like computation, geometry and algebra; see, e.g. [3], [4], [5], [6]. Shoemake [8] introduced them in the field of computer graphics. Nowadays quaternion matrices play a remarkable role in control theory, mechanics, altitude control, quantum physics and signal processing; see, e.g. [9], [10]. In skeletal animation systems, quaternions are mostly applied to interpolate between joint orientations specified with key frames or animation curves [11]. For comprehensive study of quaternions, we refer to [12].

The Hermitian solution to some matrix equations were inspected by many writers in various papers. For instance, Khatri and Mitra [29] furnished some necessary and sufficient conditions for the presence of the Hermitian solution to AY=B,A2Y=C2,YB2=D2,over C respectively, and also presented a general Hermitian solutions to (1.1) with the technique of generalized inverses. The Hermitian solution to AXA*=Bwas considered by Groß in [31].

Linear matrix equations have been one of the principal topics in matrix theory and its applications; see, e.g. [13], [14], [15], [16], [17], [18], [26], [35], [36], [37], [38], [39], [40], [41], [42], [43], [44], [45], [46]. The frequently studied Lyapunov equation BX+(BX)*=Aconstitutes a major function in mathematics, such as optimal control, stability analysis, system theory, model reduction; see, e.g, [19], [20]. The classical matrix equation CXC*+DYD*=Ahas been analyzed by numerous authors with different approaches; see, e.g. [21], [22], [23]. In [24], Yuan and Wang derived the expression of the least squares η-Hermitian solution to the real quaternion matrix expression AXB+CXD=E.

He and Wang [32] gave the general solution of A4X+(A4X)η*+B4YB4η*+C4ZC4η*=D4having η-Hermicity over H. Furthermore, Zhang and Wang [28] obtained the η-Hermitian solution to A1X=C1,A2Y=C2,ZB2=D2,A4X+(A4X)η*+B4YB4η*+C4ZC4η*=D4.Very recently, He and Wang [33] derived the η-Hermitian solution to the system A1Y=C1,YB1=D1,A2Z=C2,ZB2=D2,C3YC3η*+D3ZD3η*=A3.

Observe that the system (1.2)–(1.4) are particular cases of the following system of quaternion matrix equations A1X=C1,A2Y=C2,YB2=D2,Y=Yη*,A3Z=C3,ZB3=D3,Z=Zη*,A4X+(A4X)η*+B4YB4η*+C4ZC4η*=D4.To our best knowledge, there has been little research on the general solution of (1.5). Motivated by the above said findings, we in this paper establish the solvability conditions and the expression of the general solution to the solvable system (1.5).

The rest of this paper is composed as follows. In Section 2, we start with renowned results which are frequently experienced in the sequel. In Section 3, we offer some necessary and sufficient conditions for the presence of a solution (X, Y, Z) to (1.5), where Y and Z are η-Hermitian and present its general solution when the necessary and sufficient conditions are fulfilled. In Section 4, we discuss some particular cases of system (1.5). In Section 5, we present an algorithm and a numerical example to exemplify the notion established in this treatise. In Section 6, we give a conclusion to this paper.

Section snippets

Preliminaries

We begin with the recognized results which are useful in the sequel.

Lemma 2.1

[25]. LetAHs×t,BHs×k andCHl×t be known. Then

  • (1)

    r(A)+r(RAB)=r(B)+r(RBA)=r[AB].

  • (2)

    r(A)+r(CLA)=r(C)+r(ALC)=r[AC].

  • (3)

    r(B)+r(C)+r(RBALC)=r[ABC0].

Lemma 2.2

[32]. LetAHm×n be given. Then

  • (1)

    (Aη)=(A)η,(Aη*)=(A)η*.

  • (2)

    r(A)=r(Aη*)=r(Aη)=r(AηAη*)=r(Aη*Aη).

  • (3)

    (AA)η*=Aη*(A)η*=(AA)η=(A)ηAη.

  • (4)

    (AA)η*=(A)η*Aη*=(AA)η=Aη(A)η.

  • (5)

    (LA)η*=η(LA)η=(LA)η=LAη=RAη*.

  • (6)

    (RA)η*=η(RA)η=(RA)η=LAη*=RAη.

Lemma 2.3

[34]. Let A, B and C be given matrices with right sizes overH . Then

  • (1)

    A=(A*A

Some solvability conditions and the general solution to (1.5)

In this section, we provide some necessary and sufficient conditions for the system (1.5) to have a solution (X, Y, Z), where Y and Z are η-Hermitian matrices. Furthermore, we give its general solution when the solvability conditions are satisfied.

Theorem 3.1

Let A1, C1, A2, B2, C2, D2, A3, B3, C3, D3, A4, B4, C4 and D4=D4η*,(η{i,j,k}), be coefficient matrices in (1.5) overH with feasible sizes. Denote E=[A2B2η*],F=[C2D2η*],G=[A3B3η*],H=[C3D3η*],A5=A4LA1,Y0=EF+(EF)η*EE(EF)η*,Z0=GH+(GH)η*GG(GH)η*

Solutions to particular cases of system (1.5)

Some particular cases of system (1.5) are mentioned in this section.

If matrices A1, A2, A3, B2, B3, C1, C2, C3, D2 and D3 are equal to zero in our system, then we acquire the following outcome:

Corollary 4.1

Let A4, B4, C4 andD4=D4η*,(η{i,j,k}), be given coefficient matrices in (1.2) overH with agreeable sizes. Assume A=RA4B4,B=RA4C4,C=RA4D4(RA4)η*,M=RAB,S=BLM.Then the conditions expressed below are alike:

  • (1)

    The system (1.2) has a solution (X, Y, Z), where Y and Z are η-Hermitian matrices.

  • (2)

    The coefficient

An algorithm and a numerical example

By using Theorem 3.1, we establish an algorithm to determine the η-Hermitian solution to the system (1.5) in this section. A numerical example is also mentioned to justify our result.

Algorithm 5.1

  • (1)

    Input A1,C1,A2,B2,C2,D2,A3,B3,C3,D3,A4,B4,C4,D4=D4η* with viable dimensions over over H and η ∈ {i, j, k}.

  • (2)

    Compute the matrices E, F, G, H, A5, B5, C5, D5, A, B, C, M and S by (3.1).

  • (3)

    Check whether (3.2) and (3.3) are all agree or not. If yes, then go into the following.

  • (4)

    Evaluate X, Y and Z by (3.4)–(3.6) with U, V

Conclusion

We have furnished some necessary and sufficient conditions for the presence of the η-Hermitian solution to system (1.5) and also established an expression of the general η-Hermitian solution to (1.5) when it has a solution. The outcome of our exposition extends the noteworthy result of [28], [32], [33] and [46]. Furthermore, we have given an algorithm and a numerical example to interpret our rudimentary consequence.

Acknowledgements

I am highly grateful to the anonymous reviewers because their precious comments on the earlier version of paper really improved our paper a lot. I will also appreciate Laiyuan Gao to help me in the proofreading of my paper.

References (46)

  • Q.-W. Wang et al.

    Ranks and the least-norm of the general solution to a system of quaternion matrix equations

    Linear Algebra Appl.

    (2009)
  • A.B. Israel et al.

    Generalized Inverses: Theory and Applications

    (1974)
  • W.R. Hamilton

    On quaternions, or on a new system of imaginaries in algebra

    Philos. Mag.

    (1844)
  • J.H. Conway et al.

    On Quaternions and Octonions: Their Geometry, Arithmetic and Symmetry

    (2002)
  • G. Kamberov et al.

    Quaternions Spinors and Surfaces, Contemporary Mathematics

    (2002)
  • G. Nebe

    Finite quaternionic matrix groups

    Represent. Theory

    (1998)
  • J.P. Ward

    Quaternions and Cayley Numbers, Mathematics and Its Applications

    (1997)
  • S.R. Caradus

    Generalized Inverses and Operator Theory, Queen’s Paper in Pure and Applied Mathematics

    (1978)
  • K. Shoemake

    Animating rotation with quaternion curves

    Comput. Graph.

    (1985)
  • J.B. Kuipers

    Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality

    (2002)
  • S.L. Adler

    Quaternionic Quantum Mechanics and Quantum Fields

    (1995)
  • E.B. Dam et al.

    Quaternions, Interpolation and Animation, Department of Computer Science

    Technical report DIKU-TR-98/5

    (July 1998)
  • A.B. Özgöler et al.

    A common solution to a pair of linear matrix equations over a principle domain

    Linear Algebra Appl.

    (1991)
  • Cited by (33)

    • Weighted singular value decomposition and determinantal representations of the quaternion weighted Moore–Penrose inverse

      2017, Applied Mathematics and Computation
      Citation Excerpt :

      The definitions of the generalized inverse matrices can be extended to quaternion matrices. Generalized inverses and their using to solutions of matrix equations over the quaternion skew field (or arbitrary non-division ring) are subjects of current research (see, e.g., [3–12]). It is known various representations of the weighted Moore–Penrose.

    View all citing articles on Scopus
    View full text