Abstract
A new eigenvalue localization set for tensors is given and proved to be tighter than those presented by Qi (J. Symbolic Comput., 2005, 40, 1302-1324) and Li et al. (Numer. Linear Algebra Appl., 2014, 21, 39-50). As an application, a weaker checkable sufficient condition for the positive (semi-)definiteness of an even-order real symmetric tensor is obtained. Meanwhile, an S-type E-eigenvalue localization set for tensors is given and proved to be tighter than that presented by Wang et al. (Discrete Cont. Dyn.-B, 2017, 22(1), 187-198). As an application, an S-type upper bound for the Z-spectral radius of weakly symmetric nonnegative tensors is obtained. Finally, numerical examples are given to verify the theoretical results.
1 Introduction
For a positive integer n, n ≥ 2, N denotes the set {1, 2, …, n}. ℂ(ℝ) denotes the set of all complex (real) numbers. We call 𝓐 = (ai1⋯ im) a complex (real) tensor of order m dimension n, denoted by 𝓐 ∈ ℂ[m,n](ℝ[m,n]), if
where ij ∈ N for j = 1, 2, ⋯, m. A tensor of order m dimension n is called the unit tensor, denoted by 𝓘, if its entries are δi1⋯ im for i1, ⋯, im ∈ N, where
𝓐 is called nonnegative if ai1⋯ im ≥ 0. 𝓐 = (ai1⋯ im) ∈ ℝ[m,n] is called symmetric [1] if
where Πm is the permutation group of m indices. 𝓐 = (ai1⋯ im) ∈ ℝ[m,n] is called weakly symmetric [2] if the associated homogeneous polynomial
satisfies ∇ 𝓐xm = m𝓐xm−1. It is shown in [2] that a symmetric tensor is necessarily weakly symmetric, but the converse is not true in general.
To an n-vector x = (x1, x2⋯, xn)T, real or complex, we define the n-vector:
and
Definition 1.1
([1, 3]). Let 𝓐 = (ai1⋯ im) ∈ ℂ[m,n]. A pair (λ, x) ∈ ℂ × (ℂn ∖ {0}) is called an eigenvalue-eigenvector (or simply eigenpair) of 𝓐 if
(λ, x) is called an H-eigenpair if both of them are real.
Definition 1.2
([1, 3]). Let 𝓐 = (ai1⋯ im) ∈ ℝ[m,n]. A pair (λ, x) ∈ ℂ × (ℂn ∖ {0}) is called an E-eigenpair of 𝓐 if
(λ, x) is called an Z-eigenpair if both of them are real.
We define the Z-spectrum of 𝓐, denoted 𝓩(𝓐) to be the set of all Z-eigenvalues of 𝓐. Assume 𝓩(𝓐) ≠ 0, then the Z-spectral radius [2] of 𝓐, denoted ϱ(𝓐), is defined as
It is shown in [1] that a real even-order symmetric tensor 𝓐 = (ai1⋯ im) is positive definite if and only if all of its H-eigenvalues (Z-eigenvalues) are positive. However, when m and n are very large, it is not easy to compute all H-eigenvalues (Z-eigenvalues) of 𝓐. Then we can try to give a set in the complex which includes all H-eigenvalues (Z-eigenvalues) of 𝓐. If this set is in the right-half complex plane, then we can conclude that all H-eigenvalues (Z-eigenvalues) are positive, consequently, 𝓐 is positive definite; for details, see [1, 4–7].
There are other applications of (E-)eigenvalue inclusion sets, for example we can use them to obtain the lower and upper bounds for the H-eigenvalues (Z-spectral radius) of (nonnegative) tensors and the minimum eigenvalue of 𝓜-tensors; for details, see [8–19].
In 2005, Qi [1] presented the following Geršgorin-type eigenvalue localization set for real symmetric tensors, which can be easily extended to general tensors [4, 20].
Theorem 1.3
([1, Theorem 6]). Let 𝓐 = (ai1⋯ im) ∈ ℂ[m,n]. Then
where σ(𝓐) is the set of all eigenvalues of 𝓐 and
To get a tighter eigenvalue localization set than Γ(𝓐), Li et al. [4] proposed the following Brauer-type eigenvalue localization set for tensors.
Theorem 1.4
([4, Theorem 2.1]). Let 𝓐 = (ai1⋯ im) ∈ ℂ[m,n]. Then
where
To reduce computations, Li et al. [4] gave an S-type eigenvalue localization set by breaking N into disjoint subsets S and S, where S is the complement of S in N.
Theorem 1.5
([4, Theorem 2.2]). Let 𝓐 = (ai1⋯ im) ∈ ℂ[m,n], S be a nonempty proper subset of N. Then
In 2017, Wang et al. established the following Z-eigenvalue localization set for a real tensor 𝓐, which is completely different from eigenvalue localization sets and can be generalized to an E-eigenvalue localization set easily.
Theorem 1.6
([8, Theorem 3.1]). Let 𝓐 = (ai1⋯ im) ∈ ℝ[m,n]. Then
where
The main aim of this paper is to give a new eigenvalue localization set for tensors, which is tighter than those in Theorems 1.3-1.5, and a new E-eigenvalue localization set for tensors, which is tighter than that in Theorem 1.6. As applications, a weaker checkable sufficient condition for the positive (semi-)definiteness of an even-order real symmetric tensor is obtained based on the eigenvalue localization set, and a new upper bound for the Z-spectral radius of weakly symmetric nonnegative tensors is obtained based on the E-eigenvalue localization set.
2 A new eigenvalue localization set for tensors and its applications
In this section, we propose a new eigenvalue localization set for tensors and establish the comparisons between this set with those in Theorems 1.3-1.5. As an application of this set, we give a weaker checkable sufficient condition for the positive (semi-)definiteness of an even-order real symmetric tensor.
Theorem 2.1
Let 𝓐 = (ai1⋯ im) ∈ ℂ[m,n]. Then
where
Proof
Let λ be an eigenvalue of 𝓐 with corresponding eigenvector x = (x1, ⋯, xn)T ∈ ℂn ∖ {0}, i.e.,
Let |xp| = max{|xi| : i ∈ N}. Then, |xp| > 0. From (1), we have
Taking modulus in the above equation and using the triangle inequality give
equivalently,
If |xj| = 0, by |xp| > 0, we have
which implies that λ ∈ 𝓚p,j(𝓐) ⊆ 𝓚∩(𝓐). otherwise, |xj| > 0. Similarly, from (1), we can obtain
Multiplying (2) with (3) and noting that |xp|m−1|xj|m−1 > 0, we have
then λ ∈ 𝓚p,j(𝓐) ⊆ 𝓚∩(𝓐) . From the arbitrariness of j, we have
Next, a comparison theorem is given for Theorems 1.3-1.5 and Theorem 2.1.
Theorem 2.2
Let 𝓐 = (ai1⋯ im) ∈ ℂ[m,n], S be a nonempty proper subset of N. Then
Proof
Let S be the complement of S in N. According to Theorem 2.3 in [4], 𝓚S(𝓐) ⊆ 𝓚(𝓐) ⊆Γ(𝓐) . Hence, we only prove 𝓚∩(𝓐) ⊆ 𝓚S(𝓐) . Let z ∈ 𝓚∩(𝓐), then there exists i0 ∈ N, such that z ∈ 𝓚i0,j(𝓐), ∀ j ∈ N, j ≠ i0. If i0 ∈ S, then for any j ∈S, we have
Remark 2.3
Theorem 2.2 shows that this set in Theorem 2.1 is tighter than those in Theorem 1.3, Theorem 1.4 and Theorem 1.5, that is, 𝓚∩(𝓐) can capture all eigenvalues of 𝓐 more precisely than Γ(𝓐), 𝓚(𝓐) and 𝓚S(𝓐).
As shown in [1, 4–7], an eigenvalue localization set can provide a checkable sufficient condition for the positive (semi-)definiteness of tensors. As an application of Theorem 2.1, we give a checkable sufficient condition for the positive (semi-)definiteness of tensors.
Theorem 2.4
Let 𝓐 = (ai1⋯ im) ∈ ℝ[m,n] be an even-order symmetric tensor with ak⋯ k > 0 for all k ∈ N. If ∀ i ∈ N, ∃ j ∈ N, j ≠ i,
then 𝓐 is positive definite.
Proof
Let λ be an H-eigenvalue of 𝓐. By Theorem 2.1, we have λ ∈ 𝓚∩(𝓐), that is, there is i0 ∈ N, for any j ∈ N, j ≠ i0,
Suppose that λ ≤ 0. Then for i0 ∈ N, ∃ j0, such that ai0⋯ i0 > 0, aj0⋯ j0 > 0, and
This is a contradiction. Hence, λ > 0, and 𝓐 is positive definite. The conclusion follows. □
Similar to the proof of Theorem 2.4, the following sufficient condition is easily obtained.
Theorem 2.5
Let 𝓐 = (ai1⋯ im) ∈ ℝ[m,n] be an even-order symmetric tensor with ak⋯ k ≥ 0 for all k ∈ N. If ∀ i ∈ N, ∃ j ∈ N, j ≠ i,
then 𝓐 is positive semi-definite.
Remark 2.6
When n = 2, Theorem 2.4 is the same as Theorem 4.1 and Theorem 4.2 in [4]. When n ≥ 3, it is easy to see that the conditions of Theorem 2.4 for determining the positive definiteness of tensors are weaker than those in Theorem 4.1 and Theorem 4.2 in [4].
Next, an example is given to verify the fact in Remark 2.6.
Example 2.7
Let 𝓐 = (aijkl) ∈ ℝ[4,3] be a symmetric tensor with elements defined as follows:
By computations, we get that
Let S = {1, 2}, S = {3}. Because (7) holds, we can not use Theorem 4.1 and Theorem 4.2 in [4] to determine the positiveness of 𝓐 under this division. But from (4)-(6) and Theorem 2.4, we can determine that 𝓐 is positive definite. In fact, all the H-eigenvalues of 𝓐 are 2.9074, 3.1633, 3.7705, 4.6282 and 12.4216. By Theorem 5 in [1], 𝓐 is positive definite.
3 A new E-eigenvalue localization set for tensors and its applications
In this section, we give an S-type E-eigenvalue localization set for tensors, and establish the comparison between this set with that in Theorem 1.6. For simplification, we first denote some notations. Given a nonempty proper subset S of N, let
and then
This implies that for a tensor 𝓐 = (ai1⋯ im) ∈ ℝ[m,n], we have that for i ∈ S,
where
Theorem 3.1
Let 𝓐 = (ai1⋯ im) ∈ ℝ[m,n], S be a nonempty proper subset of N, S be the complement of S in N. Then
where σE(𝓐) is the set of all E-eigenvalues of 𝓐 and
Proof
Let λ be an E-eigenvalue of 𝓐 with corresponding eigenvector x = (x1,⋯,xn)T ∈ ℂn ∖ {0}, i.e.,
Let |xp| = max{|xi| : i ∈ S} and |xq| = max{|xj| : j ∈ S}. Then, at least one of |xp| and |xq| is nonzero. We next distinguish two cases to prove.
Case I. Suppose that
Taking modulus in the above equation and using the triangle inequality give
i.e.,
If |xp| = 0, by |xq| > 0, we have
which implies that
Multiplying (9) with (10) and noting that |xp|m−1|xq|m−1 > 0, we have
which leads to
Case II. Suppose that |xp| ≥ |xq|, then
If |xq| = 0, by |xp| > 0, we have
which implies that
Multiplying (11) with (12) and noting that |xp|m−1|xq|m−1 > 0, we have
which leads to
Theorem 3.2
Let 𝓐 = (ai1⋯ im) ∈ ℝ[m,n], S be a nonempty proper subset of N, S be the complement of S in N. Then
Proof
Let λ ∈ ΩS(𝓐) . Then
Without loss of generality, suppose that
Then there are i ∈ S and j ∈S such that
If
Furthermore,
or
which implies that λ ∈ Ri (𝓐)⋃ Rj(𝓐) . □
Remark 3.3
From Theorem 3.2, we known that the set ΩS(𝓐) in Theorem 3.1 localizes all E-eigenvalues of a tensor 𝓐 more precisely than the set
Next, based on Theorem 3.1, we give an S-type upper bound for the Z-spectral radius of a weakly symmetric nonnegative tensor.
Theorem 3.4
Let 𝓐 = (ai1⋯ im) ∈ ℝ[m,n] be a weakly symmetric nonnegative tensor, S be a nonempty proper subset of N, S be the complement of S in N. Then
where
Proof
By Lemma 4.4 in [8], ϱ(𝓐) is the largest Z-eigenvalue of 𝓐. From Theorem 3.1, we know that ϱ(𝓐) ∈ ΩS(𝓐) . Then
We next distinguish two cases to prove.
Case I: If
Then
Furthermore,
Case II: If
The conclusion follows from Cases I and II. □
Theorem 3.5
Let 𝓐 = (ai1⋯ im) ∈ ℝ[m,n] be a weakly symmetric nonnegative tensor, S be a nonempty proper subset of N, S be the complement of S in N. Then
Proof
Here, we only prove that when
Case I: For any i ∈ S,j ∈S, if Ri(𝓐) ≤ Rj(𝓐), then
Case II: For any i ∈ S,j ∈S, if Rj(𝓐) ≥ Ri(𝓐), then 0 ≤
The conclusion follows from Cases I and II. □
Remark 3.6
Theorem 3.5 shows that the upper bound in Theorem 3.4 is better than Corollary 4.5 of [9].
Now, we show that the upper bound in Theorem 3.4 is sharper than those in [8–13] in some cases by the following example.
Example 3.7
Let 𝓐 = (aijk) ∈ ℝ[3,3] with entries defined as follows:
It is not difficult to verify that 𝓐 is a weakly symmetric nonnegative tensor By computation, we obtain (ϱ(𝓐),x) = (7.3450, (0.3908, 0.6421, 0.6596)). By Corollary 4.5 of [9] and Theorem 3.3 of [10], we both have
By Theorem 3.5 of [11], we have
By Theorem 4.6 of [8], we have
By Theorem 4.7 of [8], we have
By Theorem 4.5 of [8] and Theorem 6 of [12], we both have
By Theorem 2.9 of [13], we have
Let S = {1},S = {2,3}. By Theorem 3.4, we obtain
which shows that the upper bound in Theorem 3.4 is sharper.
4 Conclusions
In this paper, we give a new eigenvalue localization set 𝓚∩(𝓐) and prove that 𝓚∩(𝓐) is tighter than those in [1] and [4]. Based on this set, we obtain a weaker checkable sufficient condition to determine the positive (semi-)definiteness for an even-order real symmetric tensor. Meanwhile, we present an S-type E-eigenvalue localization set ΩS(𝓐) and prove that ΩS(𝓐) is tighter than that in [8]. As an application, we obtain an S-type upper bound ΨS(𝓐) for the Z-spectral radius of weakly symmetric nonnegative tensors, and show that ΨS(𝓐) is sharper than those in [8–13] in some cases by a numerical example. Then an interesting problem is how to pick S to make ΨS(𝓐) as small as possible. But this is difficult when n is large. In the future, we will focus on this problem.
Acknowledgement
The authors are very indebted to the reviewers for their valuable comments and corrections, which improved the original manuscript of this paper. This work is supported by National Natural Science Foundation of China (No.11501141), Foundation of Guizhou Science and Technology Department (Grant No.[2015]2073) and Natural Science Programs of Education Department of Guizhou Province (Grant No.[2016]066).
References
[1] Qi L.Q., Eigenvalues of a real supersymmetric tensor, J. Symbolic Comput., 2005, 40, 1302–132410.1016/j.jsc.2005.05.007Search in Google Scholar
[2] Chang K.C., Qi L.Q., Zhang T., A survey on the spectral theory of nonnegative tensors, Numer. Linear Algebra Appl., 2013, 20, 891–91210.1002/nla.1902Search in Google Scholar
[3] Lim L.H., Singular values and eigenvalues of tensors: A variational approach, Proceedings of the IEEE International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP’05), 2015, 1, 129–132Search in Google Scholar
[4] Li C.Q., Li Y.T., Kong X., New eigenvalue inclusion sets for tensors, Numer. Linear Algebra Appl., 2014, 21, 39–5010.1002/nla.1858Search in Google Scholar
[5] Li C.Q., Li Y.T., An eigenvalue localizatiom set for tensor with applications to determine the positive (semi-)definitenss of tensors, Linear Multilinear Algebra, 2016, 64(4), 587–60110.1080/03081087.2015.1049582Search in Google Scholar
[6] Li C.Q., Jiao A.Q., Li Y.T., An S-type eigenvalue location set for tensors, Linear Algebra Appl., 2016, 493, 469–48310.1016/j.laa.2015.12.018Search in Google Scholar
[7] Li C.Q., Chen Z., Li Y.T., A new eigenvalue inclusion set for tensors and its applications, Linear Algebra Appl., 2015, 481, 36–5310.1016/j.laa.2015.04.023Search in Google Scholar
[8] Wang G., Zhou G.L., Caccetta L., Z-eigenvalue inclusion theorems for tensors, Discrete Cont. Dyn.-B, 2017, 22(1), 187–19810.3934/dcdsb.2017009Search in Google Scholar
[9] Song Y.S., Qi L.Q., Spectral properties of positively homogeneous operators induced by higher order tensors, SIAM. J. Matrix Anal. Appl., 2013, 34, 1581–159510.1137/130909135Search in Google Scholar
[10] Li W., Liu D.D., Vong S.W., Z-eigenvalue bounds for an irreducible nonnegative tensor, Linear Algebra Appl., 2015, 483, 182–19910.1016/j.laa.2015.05.033Search in Google Scholar
[11] He J., Bounds for the largest eigenvalue of nonnegative tensors, J. Comput. Anal. Appl., 2016, 20(7), 1290–1301Search in Google Scholar
[12] He J., Liu Y.M., Ke H., Tian J.K., Li X., Bounds for the Z-spectral radius of nonnegative tensors. Springerplus, 2016, 5(1), 172710.1186/s40064-016-3338-3Search in Google Scholar PubMed PubMed Central
[13] Liu Q.L., Li Y.T., Bounds for the Z-eigenpair of general nonnegative tensors, Open Mathematics, 2016, 14, 181–19410.1515/math-2016-0017Search in Google Scholar
[14] Huang Z.G., Wang L.G., Xu Z., Cui J.J., A new S-type eigenvalue inclusion set for tensors and its applications, J. Inequ. Appl., 2016, 2016, 25410.1186/s13660-016-1200-3Search in Google Scholar PubMed PubMed Central
[15] He J., Huang T.Z., Inequalities for 𝓜-tensors, J. Inequal. Appl., 2014, 2014, 11410.1186/1029-242X-2014-114Search in Google Scholar
[16] Zhao J.X., Sang C.L., Two new lower bounds for the minimum eigenvalue of 𝓜-tensors, J. Inequal. Appl., 2016, 2016, 26810.1186/s13660-016-1210-1Search in Google Scholar
[17] Zhou D.M., Chen G.L., Wu G.X., Zhang X.Y., On some new bounds for eigenvalues of the Hadamard product and the Fan product of matrices, Linear Algebra Appl., 2013, 438, 1415–142610.1016/j.laa.2012.09.013Search in Google Scholar
[18] Zhou D.M., Chen G.L., Wu G.X., Zhang X.Y., Some inequalities for the Hadamard product of an M-matrix and an inverse M-matrix, J. Inequal. Appl., 2013, 2013, 1610.1186/1029-242X-2013-16Search in Google Scholar
[19] Jin H.W., Kannan M.R., Bai M.R., Lower and upper bounds for the H-eigenvalues of even order real symmetric tensors, Linear Multilinear Algebra, 2017, 65(7), 1402–141610.1080/03081087.2016.1242112Search in Google Scholar
[20] Yang Y.N., Yang Q.Z., Further results for Perron-Frobenius Theorem for nonnegative tensors, SIAM. J. Matrix Anal. Appl., 2010, 31, 2517–253010.1137/090778766Search in Google Scholar
© 2017 Zhao and Sang
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.