351
Views
0
CrossRef citations to date
0
Altmetric
Research Article

Computation of the zeros of a quaternionic polynomial using matrix methods

, , , &
Pages 65-70 | Received 10 Jun 2023, Accepted 23 Dec 2023, Published online: 08 Jan 2024

Abstract

In a recent paper, Ishfaq Dar (Citation2024), worked on the problem of locating the zeros of quaternion polynomials by introducing various matrix techniques. In this paper, we use those newly developed matrix methods to locate the left eigenvalues of the quaternion companion matrix of the polynomial, which in turn yield various results concerning the location of the zeros of quaternionic polynomials including extensions of the results of Rubinstein and A. Aziz to the quaternionic setting as well.

Mathematics Subject Classification (2020):

1. Introduction

In an effort to expand complex numbers to greater spatial dimensions, the Irish mathematician Sir William Rowan Hamilton (1805-1865) invented quaternions in 1843. Hamilton became obsessed with quaternions and their uses (Hankins, Citation1980) after inventing them, and he did so for the remainder of his life. However, he probably never imagined that his invention, quaternions, would one day be used to programme video games and steer spacecraft (Turner, Citation2006). Quaternions are used in many areas of science and engineering, including computer graphics, robotics, and aerospace. In computer graphics they are particularly used in computer vision and machine learning, where they are used to represent and manipulate 3D data. In engineering they are used in the field of augmented reality, where they are used to track the position and orientation of objects in 3D space. Similarly they are used to solve problems in robotics such as inverse kinematics and path planning. Quaternions are denoted by H={a=a0+a1i+a2j+a3k:a0,a1,a2,a3R and i, j, k satisfy i2=j2=k2=ijk=1,ij=ji=k,jk=kj=i,ki=ik=j.. The set of quaternions is a skew-field and because of non-commutative nature, they differ from complex numbers C and real numbers R. A number in the quaternions is denoted by q where q=α+βi+γj+δkH, contains one real part α and three imaginary parts β, γ and δ. The conjugate of q, denoted by q¯ is a quaternion q=αβiγjδk and the norm of q is |q|=qq¯=α2+β2+γ2+δ2. The inverse of each non zero element q of H is given by q1=|q|2q¯.

Depending upon the position of the coefficients, the quaternion polynomial of degree n in indeterminate q is defined as f(q)=qn+qn1a1++qan1+an or g(q)=qn+a1qn1++an1q+an.

The quaternion Companion Matrix: The n×n companion matrix of a monic quaternion polynomial of the form f(q)=qn+qn1a1++qan1+an, is given by Cf=[0000an1000an10100an20010an30001a1], whereas, the n×n companion matrix for a monic quaternion polynomial of the form g(q)=qn+a1qn1++an1q+an, is given by Cg=[010000100001anan1an2a1]

Right Eigenvalue: Given an n×n matrix A=[aμν] of quaternions, λH is called the right eigenvalue of A, if Ax=xλ for some non zero eigenvector x=[x1,x2,xn]T of quaternions.

Left Eigenvalue: Given an n×n matrix A=[aμν] of quaternions, λH is called the left eigenvalue of A, if Ax=λx for some non zero eigenvector x=[x1,x2,,xn]T of quaternions.

For complex case, concerning the location of the eigenvalues, the famous Geršgorin theorem (Marden, Citation1949) can be stated as;

Theorem A.

All the eigenvalues of a n×n complex matrix A=(aμν) are contained in the union of n Geršgorin discs defined by Dμ={zC:|zaμμ|ν=1νμn|aμν|}.

Recently, Ishfaq Dar (Citation2024) proved the following quaternion version of Geršgorin theorem.

Theorem B.

All the left eigenvalues of a n×n matrix A=(aμν) of quaternions lie in the union of the n Geršgorin balls defined by Bμ={qH:|qaμμ|ρμ(A)}, where ρμ(A)=ν=1νμn|aμν|.

In the same paper, they considered the quaternion polynomial with coefficients on left side and gave connection between its zeros and the left eigen values of its corresponding companion matrix by proving the following result

Theorem C.

Let P(q)=qn+an1qn1++a1q+a0 be a quaternion polynomial with quaternionic coefficients and q be quaternionic variable, then for any diagonal matrix D=diag(d1,d2,,dn1,dn), where d1,d2,,dn are positive real numbers, the left eigenvalues of D1CPD and the zeros of P(q) are same.

Since variable and the coefficients are quaternions, the question arises, does above theorem holds for polynomials with coefficients on right side, the answer is provided by the following theorem.

Theorem 1.

Let f(q)=qn+qn1a1++qan1+an be a quaternion polynomial with quaternionic coefficients and q be quaternionic variable, then for any diagonal matrix T=diag(d1,d2,,dn1,dn), where d1,d2,,dn are positive real numbers, the left eigenvalues of T1CfT and the zeros of f(q) are same.

From this theorem, we conclude that the zeros of quaternion polynomial and left eigen values of its corresponding companion matrix are same irrespective of the position of its coefficients.

Remark 1.

It is easy to see that, if we take d1=d2=d3==dn=1, then both Theorem C and Theorem 1 reduce to the following result.

Corollary 1.

If λ is the left eigenvalue of the companion matrix CP associated with the quaternionic polynomial P(q), then λ is a zero of P(q).

In last two years various results were proved by different authors regarding the location of zeros of quaternion polynomials. Recently, Carney, Gardner, Keaton, and Powers (Citation2020) extended Eneström-Kakeya theorem to quaternion settings by proving following result.

Theorem D.

If p(q)=qnan+qn1an1+qn2an2++qa1+a0 is a polynomial of degree n (where q is a quaternionic variable) with real coefficients satisfying 0a0a1an, then all the zeros of p lie in |q|1.

In the same paper, they generalized of Theorem D to the polynomials whose coefficients are monotonic but not necessarily non-negative by establishing the following result.

Theorem E.

If p(q)=qnan+qn1an1+qn2an2++qa1+a0 is a polynomial of degree n (where q is a quaternionic variable) with real coefficients satisfying a0a1an, then all the zeros of p lie in |q|an+|a0|a0|an|.

Milovanović, Mir, and Ahmad (Citation2022) generalized Theorem D and Theorem E by proving the following result.

Theorem F.

If p(q)=qnan+qn1an1+qn2an2++qa1+a0 is a polynomial of degree n (where q is a quaternionic variable) with real coefficients satisfying a0a1aλ1aλan1an, where 0λn, then all the zeros of p lie in |q|2aλan+|a0|a0|an|.

Because of the restriction on the coefficients that they should be real and monotonic, the results discussed above are applicable to a small class of polynomials, so its interesting to look for the results without any restriction on the coefficients and applicable to every quaternionic polynomial with quaternion/complex or real coefficients. In this direction, Ishfaq Dar (Citation2024) proved various results concerning the location of the zeros of quaterionic polynomials with quaternionic coefficients without any restriction on the coefficients and besides this extended Cauchy’s theorem to quaternion settings by proving the following result.

Theorem G.

If f(q)=qn+qn1a1++qan1+an is a quaternion polynomial with quaternion coefficients and q is quaternionic variable, then all the zeros of f(q) lie inside the ball |q|1+max1νn|aν|.

Now in view of Theorem B and the fact that the zeros of a quaternion polynomial and the left eigenvalues of corresponding companion matrix are same, here we prove some results concerning the location of the zeros of quaternion polynomials, which extend various well known results to quaternion settings and show by examples that obtained results give better estimation of the zeros as compared to the results already available in the literature. We begin by proving the following refinement of Theorem G.

Theorem 2.

Let f(q)=qn+qn1a1+qn2a2++qan1+an be a monic quaternion polynomial of degree n with quaternion coefficients and q be a quaternion variable. If α2 α3 αn are ordered positive numbers, (1) αν=|aν|rν,ν=2,3,,n(1) where r is a positive real number. Then all the zeros of f(q) lie in the union of balls {qH:|q|r(1+α2)}and{qH:|q+a1|r}.

Next as an application of Theorem 2, we prove the following result concerning the location of the zeros of the quaternion polynomials with quaternion coefficients.

Theorem 3.

Let f(q)=qn+qn1a1+qn2a2++qan1+an be a monic quaternion polynomial with quaternion coefficients and q be a quaternion variables. If α2 α3 αn are ordered positive numbers with αν = |aν|rν,ν=2,3,,n, where r is a positive real number, then all the zeros of f(q) lie in the union of balls (2) {qH:|q|r(1+β)}and{qH:|q+a1|r},(2) where β=α2δ21+α2δ3(1+α2)2δn(1+α2)n1, δν=αναν+1,ν=2,3,,n and αn+1=0.

Remark 2.

If we assume that the variable q and coefficients a1,a2,,an in Theorem 3 are complex, then Theorem 3 reduces to a result of Zalman Rubinstein [Theorem 1, 6].

Lastly, we prove the following result.

Theorem 4.

Let f(q)=qn+qn1an1+qn2an2++qa1+a0 be a monic quaternion polynomial with quaternion coefficients and q be a quaternion variable. If λ={max0jn1|aj|}1n, then all the zeros of f(q) lie in the ball (3) {qH:|q|λ+max(λ2+λ3++λn)}.(3)

2. Computations and analysis

In this section, we present some examples for which existing Eneström-Kakeya type results are not applicable but the obtained results are applicable and give better information about the location of the zeros than existing results present in the literature. It is worth mentioning that all existing Eneström-Kakeya type results are applicable to a small class of polynomials with real coefficients satisfying monotonicity condition, whereas the results proved in this paper are applicable to a larger class of polynomial with quaternion/complex or real coefficients.

Example 1:

Let p(q)=q3+3q2+q1. It is easy to see that Theorem D and Theorem E are not applicable and on using Theorem F (with λ=2,n=3), it follows that all the zeros of p(q) lie in the ball |q+2|7. On using Theorem G, it follows that all the zeros of p(q) lie in the ball |q|4. Whereas, if we use Theorem 2 (with r=1), it follows that all the zeros of p(q) lie in the union of the balls {|q|2} and {|q+3|1}. Thus, Theorem 2 gives better bound with a significant improvement.

Example 2:

Let p(q)=q3+3q2+2q1. It is easy to see that Theorem D and Theorem E are not applicable and on using Theorem F (with λ=2,n=3), it follows that all the zeros of p(q) lie in the ball |q+2|7. On using Theorem G, it follows that all the zeros of p(q) lie in the ball |q|4. Whereas, if we use Theorem 3 (with r=1), it follows that all the zeros of p(q) lie in the union of the balls {|q|2.66} and {|q+3|1}. Thus, Theorem 3 gives better bound with a significant improvement.

3. Proof of the main theorems

Proof of Theorem 1.

Let Cf be the companion matrix of the polynomial f(q)=qn+qn1a1++qan1+an and T=diag(d1,d2,,dn1,dn) be diagonal matrix, where d1,d2,,dn are positive real numbers, then one can easily see that T1CfT=[0000andnd1d1d2000an1dnd20d2d300an2dnd3000dn1dna1].

Now if λ is the left eigenvalue of T1CfT corresponding to the non-zero eigenvector v=[v1,v2,,vn]T, then by definition of left eigenvalues (T1CfT)v=λv, implies [0000andnd1d1d2000an1dnd20d2d300an2dnd3000dn1dna1][v1v2v3vn]=λ[v1v2v3vn]

Which on simplification yields the following system of linear equations andnd1vn=λv1d1d2v1an1dnd2vn=λv2dn2dn1vn2a2dndn1vn=λvn1dn1dnvn1a1vn=λvn

Since v=[v1,v2,,vn]T being a non-zero vector, implies that vn0, hence from above equations we get an=λd1dnv1vn1an1=λd2dnv2vn1+d1dnv1vn1a2=λdn1dnvn1vn1+dn2dnvn2vn1a1=dn1dnvn1vn1λ

Multiplying above equations respectively by 1,λ,λ2,,λn1 and then on adding we get λn1a1+λn2a2++λan1+an=λn that is, λn+λn1a1+λn2a2++λan1+an=0

This shows that λ is a zero of f(q). Since λ was chosen arbitrarily, it follows that the left eigenvalues of T1CfT are the zeros of f(q).

Proof of Theorem 2.

By definition the companion matrix of the quaternion polynomial f(q)=qn+qn1a1+qn2a2++qan1+an is given by Cf=[0000an1000an10100an20001a1].

We take a diagonal matrix T = diag(rn1,rn2,,r,1), where r is a positive real number, then T1CfT=[0000anrn1r000an1rn20r00an2rn3000ra1].

Applying Theorem B to the matrix T1CfT, it follows that all the left eigenvalues of T1CfT lie in the union of balls |q||anrn1|=r|an|rn<r+r|an|rn,|q|r+|an1rn2|<r+r|an1|rn1,|q|r+|a2r|<r+r|a2|r2 and |q+a1|r.

Since by (1), αν=|aνrν|,ν=2,3,4,..,n. Therefore all the left eigenvalues of T1CfT lie in the union of balls |q|<r(1+αn),|q|r(1+αn1),|q|r(1+α2), and |q+a1|r.

Using the fact that α2α3α4αn, it follows that all left eigenvalues of T1CfT lie in the union of balls (4) {qH:|q|r(1+α2)}and{qH:|q+a1|r}.(4)

Since T is a diagonal matrix with real positive entries, by Theorem 1, it follows that all the zeros of f(q) lie in the union of the balls given by (Ishfaq Dar, Citation2024).

This completes the proof of Theorem 2.

Proof of Theorem 3.

Let w be a zero of quaternion polynomial f(q)=qn+qn1a1+qn2a2++qan1+an. Since by Theorem 2, all the zeros of quaternion polynomial f(q) lie in the union of balls (5) {qH:|q|r(1+α2)} and {qH:|q+a1|r},(5) it follows that |w|r(1+α2) and |w+a1|r. If w lies in the ball {qH:|q+a1|r}, then clearly w lies in the union of balls defined by (Milovanović et al., Citation2022). Hence we assume that the zero w does not lie in the ball {|q+a1|r}, then clearly from (Marden, Citation1949), it follows that |w|r(1+α2). Since (1+α2)(1+α2)ν,ν=1,2,3,,n1,therefore δ2(1+α2)+δ3(1+α2)2++δn(1+α2)n1δ2+δ3++δn(1+α2)=(α2α3)+(α3α4)++(αnαn+1)(1+α2)=α2(1+α2)  where  αn+1=0.

Thus β=α2{δ2(1+α2)+δ3(1+α2)2++δn(1+α2)n1}α2α2(1+α2)=α22(1+α2)>0.

In view of this, if |w|r, then |w|r(1+β), therefore we assume that |w|>r, that is, r<|w|r(1+α2). So that (6) |w|=rtwhere1<t(1+α2).(6)

Now q(w)=0 implies that |wn+wn1a1|=|(wn2a2++wan1+an)||w|n2|a2|++|w||an1|+|an|, that is, |w|n1|w+a1||w|n2|a2|++|w||an1|+|an|.

Hence by using (Rubinstein, Citation1963) and the hypothesis, we get (rt)n1|w+a1|(rt)n2|a2|++|rt||an1|+|an|=rn{tn2α2++tαn1+αn}.

This leads to tn1|w+a1|r{tn2α2++tαn1+αn}.

Multiplying both sides by (t1)tn1 we get (t1)|w+a1|r(t1)α2t+r(t1)α3t2++r(t1)αn1tn2+r(t1)αn1tn1=rα2rt(α2α3)rt2(α3α4)rtn2(αn1αn)rtn1αn=rα2rtδ2rt2δ3rtn2δn1rtn1δn.

Since t(1+α2), (t1)|w+a1|<rα2rδ2(1+α2)rδ3(1+α2)2rδn1(1+α2)n2rδn(1+α2)n1 which implies (7) (t1)|w+a1|r{α2δ2(1+α2)δ3(1+α2)2δn(1+α2)n1}.(7)

Using the fact that |w+a1|>r, and t>1, implies from inequality (Ell, Bihan, & Sangwine, Citation2014), we obtain (t1)<α2δ2(1+α2)δ3(1+α)2δn(1+α2)n1 equivalently t<1+α2δ2(1+α2)δ3(1+α)2δn(1+α2)n1=1+β.

Using this in (Rubinstein, Citation1963), we get |w|=rt<r(1+β)whereβ=α2δ2(1+α2)δ3(1+α2)2δn(1+α2)n1

Hence, it follows that all the zeros of f(q) which are outside the ball {qH:|q+a1|r} lie in the ball {qH:|q|<r(1+β)}. Therefore, it follows that all the zeros of f(q) lie in the union of balls {qH:|q|r(1+β)} and {qH:|q+a1|r}.

This completes the proof of Theorem 3.

Proof of Theorem 4.

By definition the companion matrix of the polynomial f(q)=qn+qn1an1++qa1+a0 is given by Cf=[0000a01000a10100a20001an1].

Now for any positive real number λ define a diagonal matrix T=diag(λn1,λn2,,λ,1), then T1CfT=[0000a0λn1λ000a1λn20λ00a2λn3000λan1].

Applying Theorem B to matrix T1CfT, it follows that all the left eigenvalues of T1CfT lie in the union of the balls, |q||a0|λn1,|q|λ+|a1|λn2,|q|λ+|an1|. that is, all the left eigenvalues of T1CfT lie in the ball |q|max1νn1{|a0|λn1,λ+|aν|λnν1}.

Since by hypothesis λn=max0jn1|aj|,λ0 we have |aj|λn,j=0,1,2,,n1,

This implies |q|max1νn1{λ,λ+λν+1}=max{λ,λ+λ2,λ+λ3,,λ+λn}. that is, all the left eigenvalues of the matrix T1CfT lie in the ball (8) |q|λ+max{λ2,λ3,,λn}.(8)

Since T is a diagonal matrix with real positive entries, by Theorem 1, it follows that all the zeros of f(q) lie in the ball given by (Turner, Citation2006).

This completes the proof of Theorem 4.▪

Author’s contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Acknowledgement

The author’s are highly grateful to the referees for their valuable suggestions and comments which has improved the quality of this work.

Disclosure statement

The authors declare that they have no competing interests.

References

  • Carney, N., Gardner, R., Keaton, R., & Powers, A. (2020). The Eneström-Kakeya theorem for polynomials of a quaternionic variable. Journal of Approximation Theory, 250(105325), 105325. doi:10.1016/j.jat.2019.105325
  • Ell, T. A., Bihan, N. L., & Sangwine, S. J. (2014). Quaternion Fourier transforms for signal and image processing. Willey-ISTE.
  • Hankins, T. L. (1980). Sir William Rowan Hamilton, Johns Hopkins University Press, Baltimore
  • Ishfaq Dar, N. A. R. (2024). Irfan Faiq, bounds on the zeros of a quaternionic polynomials using matrix methods. Filomat, 38(9).
  • Marden, M. (1949). Geometry of polynomials, math surveys, no. 3. Providence: American Mathematical Society.
  • Milovanović, G. V., Mir, A., & Ahmad, A. (2022). On the zeros of a quaternionic polynomial with restricted coefficients. Linear Algebra and Its Applications, 653, 231–245. doi:10.1016/j.laa.2022.08.010
  • Rubinstein, Z. (1963). Analytic method in the study of zeros of polynomials. Pacific Journal of Mathematics, 13(1), 237–249. doi:10.2140/pjm.1963.13.237
  • Turner, J. (2006, January 16). Private communication.