A general method for solving linear matrix equations of elliptic biquaternions with applications

1.   Introduction

Throughout the paper, the following notations are used. The set of real numbers, complex numbers, elliptic numbers, elliptic biquaternions are denoted by ${\mathbb{R}}, {\mathbb{C}}, {{{\mathbb{C}}}_{p}}$, $H{{{\mathbb{C}}}_{p}}$, respectively. The set of all matrices on ${\mathbb{R}}({{{\mathbb{C}}}_{p}}\, or H{{{\mathbb{C}}}_{p}})$ are denoted by ${{M}_{m\times n}}({\mathbb{R}})({{M}_{m\times n}}({{{\mathbb{C}}}_{p}})\, or{{M}_{m\times n}}(H{{{\mathbb{C}}}_{p}}))$. For convenience, the set of all square matrices on ${{{\mathbb{C}}}_{p}} (\, or\, H{{{\mathbb{C}}}_{p}})$ are denoted by ${{M}_{s}}({{{\mathbb{C}}}_{p}})(\, or\, {{M}_{s}}(H{{{\mathbb{C}}}_{p}}))$.

In 1843, Hamilton introduced the set of real quaternions ^[1], which can be represented as

where the quaternionic units $\mathbf{i}$, $\mathbf{j}$ and $\mathbf{k}$ satisfy the equalities:

There are many applications of real quaternions in various areas of science. One of these applications is related to matrix theory. In the first half of the 20th century, the real quaternion matrices began to study ^[2]. It is well known that linear matrix equations with their applications have been one of the main topics in matrix theory. For real quaternion matrices, Tian discussed the linear equations in ^[3], and gave a general method to solve them. Also, the real quaternion matrix equation $X-AXF = C$ (Lyapunov equation) is studied by Song et al. ^[4]. On the other hand, in ^[5], Song and Chen studied the real quaternion matrix equation $XF-AX = C$ (Sylvester equation). In the process of preparing this paper, we are motivated by the aforementioned studies ^[3,4,5]. The studies ^[6,7,8,9] can be suggested as some different and qualified studies on quaternion matrix equations.

After the discovery of real quaternion algebra, Hamilton also introduced the complex quaternion algebra ^[10]. The set of complex quaternions is defined by

where $\mathbf{i}, \mathbf{j}$ and $\mathbf{k}$ satisfy the same multiplication rules given in (1.1).

As well as real quaternions, complex quaternions have many applications in many areas of science and they have an important role to explain mathematical and physical events. In these applications of real and complex quaternions, their complex matrix representations have an important place. There can be found some interesting studies on complex matrices in ^{[11,12,13,14]}.

Recently, the elliptic biquaternion algebra, which includes the complex quaternion algebra and real quaternion algebra as special cases, has been introduced. Various studies concerned with elliptic biquaternion algebra have been presented in the literature. We refer the readers to the papers ^{[15,16,17,18,19]}.

This article is organized as follows. In section 2, we review elliptic numbers, elliptic matrices, elliptic biquaternions and elliptic biquaternion matrices. In Section 3, real representations of elliptic biquaternion matrices are obtained. In section 4, in view of these real representations, we develop a general method to study the solutions of linear matrix equations over the elliptic biquaternion algebra $H{{{\mathbb{C}}}_{p}}$. In section 5, we investigate the solutions of the elliptic biquaternion matrix equations $X-AXB = C$ and $AX-XB = C$ by means of this method. In section 6, we provide numerical algorithms for finding the solutions of problems which are discussed in the section 4 and section 5.

2.   Preliminaries

The set of elliptic numbers is represented as

In this number system, addition and multiplication of any elliptic numbers $\omega = {{x}_{1}}+I{{y}_{1}}, \, \, \varsigma = {{x}_{2}}+I{{y}_{2}}\in {{{\mathbb{C}}}_{p}}$ are defined as $\omega +\varsigma = ({{x}_{1}}+I{{y}_{1}})+({{x}_{2}}+I{{y}_{2}}) = ({{x}_{1}}+{{x}_{2}})+I({{y}_{1}}+{{y}_{2}})$ and $\omega \varsigma = ({{x}_{1}}+I{{y}_{1}})({{x}_{2}}+I{{y}_{2}}) = \left({{x}_{1}}{{x}_{2}}+p{{y}_{1}}{{y}_{2}} \right)+I\left({{x}_{1}}{{y}_{2}}+{{x}_{2}}{{y}_{1}} \right)$, respectively. As it is well known in the literature, ${{{\mathbb{C}}}_{p}}$ is a field under these two operations, ^[20]. The set of matrices, which includes $m\times n$ matrices with elliptic number entries, are investigated in ^[21]. In the set of $m\times n$ elliptic matrices ${{M}_{m\times n}}\left({{{\mathbb{C}}}_{p}} \right)$, ordinary matrix addition and multiplication are defined. Also, the scalar multiplication is defined as $\lambda A = \lambda \left[{{a}_{ij}} \right] = \left[\lambda {{a}_{ij}} \right]\in {{M}_{m\times n}}\left({{{\mathbb{C}}}_{p}} \right)$ where $\lambda \in {{{\mathbb{C}}}_{p}}$ and $A = \left[{{a}_{ij}} \right]\in {{M}_{m\times n}}\left({{{\mathbb{C}}}_{p}} \right)$, ^[21].

The elliptic biquaternion algebra is a four dimensional vector space over the elliptic number field ${{{\mathbb{C}}}_{p}}$. It is expressed by

where $\mathbf{i}, \mathbf{j}$ and $\mathbf{k}$ are the quaternionic units which satisfy (1.1). Let $Q = {{A}_{0}}+{{A}_{1}}\mathbf{i}+{{A}_{2}}\mathbf{j}+{{A}_{3}}\mathbf{k},$ $\, R = {{B}_{0}}+{{B}_{1}}\mathbf{i}+{{B}_{2}}\mathbf{j}+{{B}_{3}}\mathbf{k}\in H{{{\mathbb{C}}}_{p}}$ and $\lambda \in {{{\mathbb{C}}}_{p}}$ be given. Then, the operations of multiplication and addition are expressed as

while the operation of scalar multiplication is expressed as ^[15]

The set of all $m\times n$ type matrices with elliptic biquaternion entries is denoted by ${{M}_{m\times n}}\left(H{{{\mathbb{C}}}_{p}} \right)$. In this set, the ordinary matrix addition and multiplication are defined. Also, the scalar multiplication is defined as

where $A = \left[{{a}_{ij}} \right]\in {{M}_{m\times n}}\left(H{{{\mathbb{C}}}_{p}} \right)$ and $Q\in H{{{\mathbb{C}}}_{p}}$. There is a faithful relation between elliptic matrices and elliptic biquaternion matrices. Where ${{A}_{0}}, \, \, {{A}_{1}}, \, {{A}_{2}}, \, {{A}_{3}}\in {{M}_{m\times n}}\left({{{\mathbb{C}}}_{p}}\right)$, every elliptic biquaternion matrix $A = {{A}_{0}}+{{A}_{1}}\mathbf{i}+{{A}_{2}}\mathbf{j}+{{A}_{3}}\mathbf{k}\in{{M}_{m\times n}}\left(H{{{\mathbb{C}}}_{p}} \right)$ has a $2m\times2n$ elliptic representation

which is determined by means of the following linear isomorphism ^[22]

We aim to obtain the real representations of elliptic biquaternion matrices. In the next section, the representation $\xi \left(A \right)$ will be written in a somewhat different form which is suitable for the purpose of us.

3.   Real representations of elliptic biquaternion matrices

In this section, firstly, we get the real representations of elliptic matrices. Afterwards, we obtain the real representations of elliptic biquaternion matrices which will be useful for investigating the solutions of linear matrix equations over the elliptic biquaternion algebra $H{{{\mathbb{C}}}_{p}}$ in the next section.

By means of the study ^[23] which was presented by Yaglom in 1968, we know that in the case ${{I}^{2}} = -q-rI({{r}^{2}}-4q < 0)$, the transformation that is obtained by making the generalized complex number ${{c}_{1}}+I{{d}_{1}}$ correspond to the ordinary complex number $c+id$, where $c = {{c}_{1}}-\frac{r}{2}{{d}_{1}}$ and $d = \frac{{{d}_{1}}}{2}\sqrt{4q-{{r}^{2}}}$, is an isomorphism. As it is well known, ${{I}^{2}} = p, \, p < 0$ for elliptic numbers. By taking into consideration this case, the restriction of this isomorphism is obtained as follows:

On the other hand, we know that an ordinary complex number $z = a+ib$ has a faithful real matrix representation

which is determined by means of the following linear isomorphism

where $M_{2}^{\Omega }({\mathbb{R}}) = \left\{ \left[\begin{array}{cc} x & -y \\ y & x \\ \end{array} \right]:\, \, x, y\in {\mathbb{R}} \right\}$, ^[24].

Now, we take into consideration the compound function $\delta_{p} = \alpha \circ \varepsilon$ which is given as

Since the functions $\varepsilon$ and $\alpha$ are linear isomorphisms, there is no doubt that $\delta_{p}$ is a linear isomorphism. Thus, we have a faithful real matrix representation of an elliptic number $\omega = x+Iy\in {{{\mathbb{C}}}_{p}}$ as

As a natural consequence of the linear isomorphism ${{\delta }_{p}}$, the following function

that is expected to be a linear isomorphism can be immediately defined where $M_{2m\times 2n}^{\Omega }({\mathbb{R}}) = \left\{\left[\begin{array}{cc} G & -\sqrt{\left| p \right|}H \\ \sqrt{\left| p \right|}H & G \\ \end{array} \right]:G, H\in {{M}_{m\times n}}({\mathbb{R}}) \right\}$. One can see easily that this function is bijection and satisfies the following equalities

for any elliptic matrices $A, \, B, \, C$ of appropriate sizes. Thus, $\gamma _{p}$ is a linear isomorphism as anticipated and we have a faithful real matrix representation of an elliptic matrix $A = {{A}_{1}}+I{{A}_{2}}\in{{M}_{m\times n}}({{{\mathbb{C}}}_{p}})$ as

On the other hand, the elliptic matrix representation of an elliptic biquaternion matrix $A = {{A}_{0}}+{{A}_{1}}\mathbf{i}+{{A}_{2}}\mathbf{j}+{{A}_{3}}\mathbf{k}\in {{M}_{m\times n}}(H{{{\mathbb{C}}}_{p}})$ given in (2.1) can be written in a somewhat different form as

where ${{A}_{i}} = {{A}_{i}}^{\#}+I\, {{A}_{i}}^{\prime }\in {{M}_{m\times n}}({{{\mathbb{C}}}_{p}})\,$, ${{A}_{i}}^{\#}, \, {{A}_{i}}^{\prime }\in {{M}_{m\times n}}({\mathbb{R}})\,$, $0\le i\le 3$. Then, applying (3.4) to (3.5), we get $4m\times4n$ real matrix representation of the elliptic matrix $\xi \left(A \right)$ (in other words the real representation of the elliptic biquaternion matrix $A$) as follows:

For convenience, let us denote ${{\gamma }_{p}}(\xi(A))$ by $({A})_{\gamma p}$ where $A$ is any elliptic biquaternion matrix. One can immediately see that the unit matrix ${{I}_{n}}$ satisfies the equation

Some properties which are satisfied by the real representation are given below.

Proposition 3.1. Let $A, \, B\in{{M}_{m\times n}}(H{{{\mathbb{C}}}_{p}})$, $C\in{{M}_{n\times l}}(H{{{\mathbb{C}}}_{p}})$, $D\in{{M}_{s}}(H {{{\mathbb{C}}}_{p}})$ be arbitrary elliptic biquaternion matrices. In that case

1. $A = B\Leftrightarrow {({A})_{\gamma p}} = {({B})_{\gamma p}}$,

2. ${{(A+B)}_{\gamma p}} = {({A})_{\gamma p}}+{({B})_{\gamma p}}, \, \, {{(AC)}_{\gamma p}} = {({A})_{\gamma p}}{({C})_{\gamma p}}$,

3. If $D$ is invertible, then ${({D})_{\gamma p}}$ is invertible and ${{({{D}^{-1}})}_{\gamma p}} = {{({{(D)}_{\gamma p}})}^{-1}}$,

4. ${({A})_{\gamma p}} = {{S}_{4m}}^{-1}{({A})_{\gamma p}}{{S}_{4n}}$ where ${{S}_{4t}} = \left[\begin{array}{cccc} 0 & 0 & -{{I}_{t}} & 0 \\ 0 & 0 & 0 & -{{I}_{t}} \\ {{I}_{t}} & 0 & 0 & 0 \\ 0 & {{I}_{t}} & 0 & 0 \\ \end{array} \right], \, \, \, t = m, \, n.$

Proof. Since ${{\gamma }_{p}}\circ \varepsilon$ is a linear isomorphism, 1 and 2 are obvious. Also, the proof of 4 can be easily completed by direct calculation. Now we will prove 3.

$3.$ From the inverse property, we can write

Then, we get the equalities

and

by means of (3.7) and first two properties in this proposition. It means that ${{({{D}^{-1}})}_{\gamma p}} = {{({{(D)}_{\gamma p}})}^{-1}}$.

4.   On the solutions of linear matrix equations over $H{{{\mathbb{C}}}_{p}}$

In this section, we give a general method on the solutions of linear matrix equations over the elliptic biquaternion algebra $H{{{\mathbb{C}}}_{p}}$ with the aid of the real representations. To do so, we take into consideration the elliptic biquaternion matrix equation:

where ${{A}_{1}}, ..., {{A}_{k}}\in {{M}_{m\times n}}(H{{{\mathbb{C}}}_{p}}), \, {{B}_{1}}, ..., {{B}_{k}}\in {{M}_{u\times v}}(H{{{\mathbb{C}}}_{p}}), \, C\in {{M}_{m\times v}}(H{{{\mathbb{C}}}_{p}})$ and $X\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$.

Let us define the real representation of the elliptic biquaternion matrix equation (4.1) as in the following:

Thanks to the first two properties given in Proposition 3.1, the elliptic biquaternion matrix equation (4.1) is equivalent to the following real matrix equation

Hence, we have the following proposition.

Proposition 4.1. The elliptic biquaternion matrix equation $(4.1)$ has an elliptic biquaternion matrix solution $X\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$ if and only if the real matrix equation $(4.2)$ has a real matrix solution $Y = {({X})_{\gamma p}}\in {{M}_{4n\times 4u}}({\mathbb{R}})$.

Theorem 4.2. Let ${{A}_{1}}, {{A}_{2}}, ..., {{A}_{k}}\in {{M}_{m\times n}}(H{{{\mathbb{C}}}_{p}}), \, {{B}_{1}}, {{B}_{2}}, ..., {{B}_{k}}\in {{M}_{u\times v}}(H{{{\mathbb{C}}}_{p}})$ and $C\in {{M}_{m\times v}}(H{{{\mathbb{C}}}_{p}})$ be given. In this case the elliptic biquaternion matrix equation $(4.1)$ has a solution $X\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$ if and only if its real representation equation $(4.2)$ has a solution $Y\in {{M}_{4n\times 4u}}({\mathbb{R}})$. In that case, if the block matrix

is a solution of $(4.2)$, then $n\times u$ elliptic biquaternion matrix

is a solution of $(4.1)$ where

Proof. In view of the Proposition 4.1, the proof remains to show that if the block real matrix

is a solution of (4.2), in that case the elliptic biquaternion matrix, which is given in (4.4), is a solution of (4.1). When $Y$ is a solution of (4.2), in view of the fourth property in Proposition 3.1, we have the equation

This last equation shows that ${{S}_{4n}}^{-1}Y{{S}_{4u}}$ is also a solution of (4.2). Then, according to the matrix theory, the following matrix

satisfies the real matrix equation (4.2), that is, ${Y}'$ is another solution of (4.2). If (4.6) is substituted in (4.8), after some calculations, the equality

is obtained where

By taking into consideration (4.9) and (3.6), we construct the elliptic biquaternion matrix

where ${{X}^{\#}_{i}}, \, {{X}^{\prime }_{i}}, \, 0\le i\le 3$ are as in $(4.5)$.

Obviously, ${({X})_{\gamma p}} = {Y}'$. Then, according to Proposition 4.1, the elliptic biquaternion matrix $X$ is a solution of the equation (4.1).

5.   Some results

In this section, we investigate the solutions of the elliptic biquaternion matrix equations $X-AXB = C$ and $AX-XB = C$ by means of Theorem 4.2.

For $k = 2$, the special case of (4.1) is given by

where ${{A}_{1}}, {{A}_{2}}\in {{M}_{m\times n}}(H{{{\mathbb{C}}}_{p}}), \, \, {{B}_{1}}, {{B}_{2}}\in {{M}_{u\times v}}(H{{{\mathbb{C}}}_{p}}), \, \, C\in {{M}_{m\times v}}(H{{{\mathbb{C}}}_{p}})$ and $X\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$. If ${{B}_{1}} = {{I}_{u}}, \, \, {{A}_{2}} = -{{I}_{n}}, \, \, m = n, \, \, u = v$ are taken in (5.1) and also some notation changes are made as follows: ${{A}_{1}} = A, \, \, {{B}_{2}} = B$,

is obtained where $A\in {{M}_{n}}(H{{{\mathbb{C}}}_{p}}), \, \, B\in {{M}_{u}}(H{{{\mathbb{C}}}_{p}})$ and $C\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$. It is not difficult to see that the real representation of the last equation is

In view of the above derivation and Theorem 4.2, we have the following corollary:

Corollary 5.1. Let $A\in {{M}_{n}}(H{{{\mathbb{C}}}_{p}}), \, \, B\in {{M}_{u}}(H{{{\mathbb{C}}}_{p}})$ and $C\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$. In this case the elliptic biquaternion matrix equation

has a solution $X\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$ if and only if the real matrix equation

has a solution $Y\in {{M}_{4n\times 4u}}({\mathbb{R}})$, in which case, if

is a solution of $(5.3)$, then $n\times u$ elliptic biquaternion matrix

is a solution of $(5.2)$ where ${{X}^{\#}_{i}}, \, {{X}^{\prime }_{i}}, \, 0\le i\le 3$ are calculated as in $(4.5)$.

Similarly above, if ${{A}_{1}} = {{I}_{n}}, \, \, {{B}_{1}} = {{I}_{u}}, \, \, m = n, \, \, u = v$ are taken in (5.1) and also some notation changes are made as follows: ${{A}_{2}} = -A, \, \, {{B}_{2}} = B$,

is obtained where $A\in {{M}_{n}}(H{{{\mathbb{C}}}_{p}}), \, \, B\in {{M}_{u}}(H{{{\mathbb{C}}}_{p}})$ and $C\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$. It is easy to see that the real representation of the last equation is

Considering the above derivation and Theorem 4.2, we have the following corollary:

Corollary 5.2. Let $A\in {{M}_{n}}(H{{{\mathbb{C}}}_{p}}), \, \, B\in {{M}_{u}}(H{{{\mathbb{C}}}_{p}})$ and $C\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$. In this case the elliptic biquaternion matrix equation

has a solution $X\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$ if and only if the real matrix equation

has a solution $Y\in {{M}_{4n\times 4u}}({\mathbb{R}})$, in which case, if

is a solution of $(5.5)$, then $n\times u$ elliptic biquaternion matrix

is a solution of $(5.4)$ where ${{X}^{\#}_{i}}, \, {{X}^{\prime }_{i}}, \, 0\le i\le 3$ are calculated as in $(4.5)$.

6.   Numerical algorithms and examples

Based on the discussions in Section 4 and Section 5, in this section we provide numerical algorithms for finding the solutions of problems which are related to Corollary 5.1, Corollary 5.2 and Theorem 4.2.

Note that all computations in the rest of the paper are performed on an Intel i7-3630QM@2.40 Ghz/16GB computer using MATLAB R2016a software. Another thing that can be of importance is that we use the standard MATLAB package procedures.

Firstly, we give an example related to Corollary 5.1.

Example 6.1. Solve elliptic biquaternion matrix equation

over the elliptic biquaternion algebra $H{{{\mathbb{C}}}_{-9}}$ with the aid of its real representation.

By taking into consideration the Corollary 5.1, real representation of given matrix equation can be written as in the following:

If we solve this equation, we have

It means that

Consequently, it is concluded that

by means of the equations (4.4) and (4.5).

We can generate an algorithm for problems related to Corollary 5.1 as in the following:

Algorithm 1

1. Input $A, \, B, \, C\, (A\in {{M}_{n}}(H{{{\mathbb{C}}}_{p}}), \, B\in {{M}_{u}}(H{{{\mathbb{C}}}_{p}})$ and $C\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$).

2. Form ${({A})_{\gamma p}}, \, {({B})_{\gamma p}}, \, {({C})_{\gamma p}}$.

3. Compute $Y = \left[{{Y}_{ij}} \right]_{i, j = 1}^{4}$ satisfying (5.3) (${{Y}_{ij}}\in {{M}_{n\times u}}({\mathbb{R}})$).

4. Calculate ${{X}^{\#}_{i}}, \, {{X}^{\prime }_{i}}$ according to (4.5)($0\le i\le 3$).

5. Output $X = ({{X}^{\#}_{0}}+I{{X}^{\prime }_{0}})+({{X}^{\#}_{1}}+I{{X}^{\prime }_{1}})\mathbf{i}+({{X}^{\#}_{2}}+I{{X}^{\prime }_{2}})\mathbf{j}+({{X}^{\#}_{3}}+I\, {{X}^{\prime }_{3}})\mathbf{k}$.

Now, we give an example related to Corollary 5.2.

Example 6.2. Solve elliptic biquaternion matrix equation

over the elliptic biquaternion algebra $H{{{\mathbb{C}}}_{-2}}$ by using its real representation.

If the Corollary 5.2 is taken into consideration, the real representation equation

can be written. By solving this equation, we have

It means that

Consequently, we obtain

by means of the equations (4.4) and (4.5).

We can generate an algorithm for problems related to Corollary 5.2 as follows:

Algorithm 2

1. Input $A, \, B, \, C\, (A\in {{M}_{n}}(H{{{\mathbb{C}}}_{p}}), \, B\in {{M}_{u}}(H{{{\mathbb{C}}}_{p}})$ and $C\in {{M}_{n\times u}}(H{{{\mathbb{C}}}_{p}})$).

2. Form ${({A})_{\gamma p}}, \, {({B})_{\gamma p}}, \, {({C})_{\gamma p}}$.

3. Compute $Y = \left[{{Y}_{ij}} \right]_{i, j = 1}^{4}$ satisfying (5.5) (${{Y}_{ij}}\in {{M}_{n\times u}}({\mathbb{R}})$).

4. Calculate ${{X}^{\#}_{i}}, \, {{X}^{\prime }_{i}}$ according to (4.5)($0\le i\le 3$).

5. Output $X = ({{X}^{\#}_{0}}+I{{X}^{\prime }_{0}})+({{X}^{\#}_{1}}+I{{X}^{\prime }_{1}})\mathbf{i}+({{X}^{\#}_{2}}+I{{X}^{\prime }_{2}})\mathbf{j}+({{X}^{\#}_{3}}+I\, {{X}^{\prime }_{3}})\mathbf{k}$.

As a result of using the standard MATLAB package procedures, when our calculations include the rational numbers, root numbers, exponential expressions, logarithmic expressions etc., our method gives an approximate solution of the desired equation just like in the case of Example 6.2. Otherwise, our method gives the exact solution of the desired equation just like in the case of Example 6.1. To ensure the exact solution of the desired equation in Example 6.2, one can use the MATLAB package Symbolic Math Toolbox. If the same steps are followed by using this package, the exact solution

of the aforementioned elliptic biquaternion matrix equation is immediately found. It must be noted that using this package always provides an advantage in terms of the exact solution, however using it sometimes causes a disadvantage in terms of the length of the solution.

Finally, we can give an algorithm for the most general case, that is, for the problems related to Theorem 4.2.

Algorithm 3

1. Input ${{A}_{i}}, \, {{B}_{i}}, \, C\, ({{A}_{i}}\in {{M}_{m\times n}}(H{{{\mathbb{C}}}_{p}}), \, {{B}_{i}}\in {{M}_{u\times v}}(H{{{\mathbb{C}}}_{p}}), \, 1\le i\le k$ and $C\in {{M}_{m\times v}}(H{{{\mathbb{C}}}_{p}})$).

2. Form ${{({{A}_{i}})}_{\gamma p}}, \, \, {{({{B}_{i}})}_{\gamma p}}, \, \, {({C})_{\gamma p}}\, (1\le i\le k)$.

3. Compute $Y = \left[{{Y}_{ij}} \right]_{i, j = 1}^{4}$ satisfying (4.2) (${{Y}_{ij}}\in {{M}_{n\times u}}({\mathbb{R}})$).

4. Calculate ${{X}^{\#}_{i}}, \, {{X}^{\prime }_{i}}$ according to (4.5)($0\le i\le 3$).

5. Output $X = ({{X}^{\#}_{0}}+I{{X}^{\prime }_{0}})+({{X}^{\#}_{1}}+I{{X}^{\prime }_{1}})\mathbf{i}+({{X}^{\#}_{2}}+I{{X}^{\prime }_{2}})\mathbf{j}+({{X}^{\#}_{3}}+I\, {{X}^{\prime }_{3}})\mathbf{k}$.

7.   Conclusion

In this paper, real representations of elliptic biquaternion matrices, which may be needed to investigate various topics on elliptic biquaternion matrices in the future, are obtained. By means of these representations, a general method on the solutions of linear matrix equations over the elliptic biquaternion algebra $H{{{\mathbb{C}}}_{p}}$ is developed. Also, some problems are considered as applications of this method. Lastly, the numerical algorithms for finding the solutions of these problems are provided.

When $p = -1$ the number system ${{{\mathbb{C}}}_{p}}$ and the set of elliptic biquaternions $H{{{\mathbb{C}}}_{p}}$ correspond to the complex number system $\mathbb{C}$ and the set of complex quaternions ${{H}_{\mathbb{C}}}$, respectively. As a natural consequence of this case, the set of elliptic biquaternion matrices ${{M}_{m\times n}}(H{{{\mathbb{C}}}_{p}})$ is reduced to set of complex quaternion matrices ${{M}_{m\times n}}({{H}_{\mathbb{C}}})$ when $p = -1$. Therefore, our method solves the linear equations of complex quaternion matrices as well.

Real or complex quaternion matrices have an important role in many areas of science. Since elliptic biquaternion matrices are generalized form of complex quaternion matrices and so real quaternion matrices, it is expected that the results obtained here will be used as a valuable tool in many areas of science.

Acknowledgements

The author would like to thank to Professor Murat Tosun and Assistant Professor Hidayet Hüda Kösal for their help and useful discussions.

Conflict of interest

The author declares that there is no conflict of interest.