Home Mathematics An ℋ-tensor-based criteria for testing the positive definiteness of multivariate homogeneous forms
Article Open Access

An ℋ-tensor-based criteria for testing the positive definiteness of multivariate homogeneous forms

  • Dongjian Bai and Feng Wang EMAIL logo
Published/Copyright: June 25, 2022

Abstract

A positive definite homogeneous multivariate form plays an important role in the field of optimization, and positive definiteness of the form can be identified by a special structured tensor. In this paper, based on the equivalence between the form and the corresponding tensor, and the links of the positive definiteness of a tensor with ℋ-tensor, we propose an ℋ-tensor-based criterion for identifying the positive definiteness of multivariate homogeneous forms. Some numerical examples are provided to illustrate the efficiency and validity of our results.

MSC 2010: 15A69; 15A18; 65F15; 65H17

1 Introduction

Let R ( C ) be the real(complex) field, N = { 1 , 2 , , n } . An m th-order n -dimensional real(complex) tensor A = ( a j 1 j 2 j m ) consists of n m real(complex) entrise:

a j 1 j 2 j m R ( C ) ,

where j i = 1 , 2 , , n , i = 1 , 2 , , m [1,2,3, 4,5,6]. Obviously, a matrix is a 2nd-order tensor. Moreover, tensor A = ( a j 1 j 2 j m ) is called symmetric [7] if

a j 1 j 2 j m = a π ( j 1 j 2 j m ) , π Π m ,

where Π m is the permutation group of m indices. If a j 1 j 2 j m 0 , then tensor A is called a nonnegative tensor. Tensor = ( δ j 1 j 2 j m ) is called the unit tensor [8], where

δ j 1 j 2 j m = 1 , if j 1 = = j m , 0 , otherwise .

Denote an m -th degree multivariate homogeneous form of n variables f ( x ) as follows:

(1) f ( x ) = j 1 , j 2 , , j m N a j 1 j 2 j m x j 1 x j 2 x j m ,

where x R n . While m is even, f ( t ) is labeled positive definite if

f ( t ) > 0 , for any t R n , t 0 .

Function f ( x ) in (1) can be expressed as the tensor product of an m th-order n -dimensional symmetric tensor A and x m defined as follows:

(2) f ( x ) A x m = j 1 , j 2 , , j m N a j 1 j 2 j m x j 1 x j 2 x j m ,

where x R n and x m is an m th-order n -dimensional rank-one tensor with entries x j 1 x j m [5]. The symmetric tensor A is positive definite if f ( x ) is positive definite [9].

The positive definiteness of tensor has received much attention of researchers’ in recent decade [10,11,12]. By the Sturm theorem, the positive definiteness of a multivariate homogeneous form can be identified for n 3 [13]. Ni et al. [9] presented an eigenvalue method for checking positive definiteness of a symmetric tensor. While, all the eigenvalues should be calculated in this method, and it is not practical when tensor order or dimension is large.

Recently, Li et al. [14] provided an -tensor-based method for identifying the positive definiteness of an even-order symmetric tensor. It is well known that an even-order symmetric -tensor with positive diagonal entries is positive definite. Hence, we can check the positive definiteness of a tensor via the aid of -tensor. Subsequently, with the help of generalized diagonal dominance, miscellaneous criterions for -tensors and -tensors are provided [15,16,17, 18,19,20], which depends on the entries of tensors and is efficient to determine -tensor ( -tensor).

Theorem 1

[20] Let A = ( a j 1 j m ) C [ m , n ] . If

a j j j r j > j 2 j 3 j m N 0 m 1 a j j 2 j m + j 2 j 3 j m N 2 m 1 δ j j 2 j m = 0 max t { j 2 , j 3 , , j m } { r t } a j j 2 j m + j 2 j 3 j m N 3 m 1 max t { j 2 , j 3 , , j m } { l t } a j j 2 j m , j N 2

and

a i i i i 2 i 3 i m N 0 m 1 δ i i 2 i m = 0 a i i 2 i m , i N 1 ( o r N 1 = ) ,

then A is an -tensor.

In this paper, we continue to present new criteria for -tensors. These new results improve the corresponding conclusions [20,21,22]. As an application, some sufficient conditions of the positive definiteness for an even-order real symmetric tensor are given. Several numerical experiments demonstrate its efficiency.

Several specially definitions extended from matrices are as follows.

Definition 1

[17] Let A = ( a j 1 j 2 j m ) be a complex tensor of order m dimension n . A is an -tensor if there exists a positive vector x = ( x 1 , x 2 , , x n ) T R n , such that

a j j j x j m 1 > j 2 , , j m N δ j j 2 j m = 0 a j j 2 j m x j 2 x j m , j N .

Definition 2

[8] A complex tensor A = ( a j 1 j 2 j m ) of order m dimension n is called reducible if there is a nonempty subset I N , such that

a j 1 j 2 j m = 0 , j 1 I , j 2 , , j m I .

Otherwise, A is called irreducible.

Definition 3

[18] Let A = ( a j 1 j 2 j m ) be a complex tensor of order m dimension n . For i , j N ( i j ) , if there are indices k 1 , k 2 , , k r , such that

j 2 , , j m N δ k s j 2 j m = 0 , k s + 1 { j 2 , , j m } a k s j 2 j m 0 , s = 0 , 1 , , r .

where k 0 = j , k r + 1 = i , and then there exists a nonzero elements chain from j to i .

The remainder of this paper is as follows. In Section 2, we gives some criteria for the identification of -tensors. In Section 3, some new conditions for the identification of positive definite tensors are presented. Some numerical experiments are also provided to show the effectiveness of new methods.

2 Criteria for -tensors

In this section, three new criteria for identifying -tensors are presented. First, for the convenience of description, some notation and lemmas are given. For a tensor A = ( a j 1 j 2 j m ) of order m dimension n , denote

S m 1 = { j 2 j 3 j m : j i S , i = 2 , 3 , , m } , S N ; N m 1 S m 1 = { j 2 j 3 j m : j 2 j 3 j m N m 1 and j 2 j 3 j m S m 1 } ; N 0 m 1 = N m 1 ( N 2 m 1 N 3 m 1 ) ; Λ j ( A ) = j 2 , , j m N δ j j 2 j m = 0 a j j 2 j m = j 2 , , j m N a j j 2 j m a j j j ; N 1 = { j N : 0 < a j j j = Λ j ( A ) } ; N 2 = { j N : 0 < a j j j < Λ j ( A ) } ; N 3 = { j N : a j j j > Λ j ( A ) } ; w j = Λ j ( A ) Λ j ( A ) + a j j j , w = max { w j } , j N 2 ; r = max j N 3 j 2 j m N 0 m 1 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m a j j j j 2 j m N 3 m 1 δ j j 2 j m = 0 a j j 2 j m ;

P j , r ( A ) = j 2 j m N 0 m 1 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m + r j 2 j m N 3 m 1 δ j j 2 j m = 0 a j j 2 j m , j N 3 ; r 1 = max j N 3 j 2 j m N 0 m 1 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m a j j j j 2 j m N 3 m 1 δ j j 2 i m = 0 max t { j 2 , j 3 , , j m } P t , r ( A ) a t t t a j j 2 j m ; μ = max { w j , r 1 } , j N 2 ;

P j , r 1 ( A ) = j 2 j m N 0 m 1 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m + r 1 j 2 j m N 3 m 1 δ j j 2 j m = 0 max t { j 2 , j 3 , , j m } P t , r ( A ) a t t t a j j 2 j m , j N 3 ; h = max j N 3 μ j 2 j m N 0 m 1 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m w P j , r 1 ( A ) j 2 j m N 3 m 1 δ j j 2 j m = 0 max t { j 2 , j 3 , , j m } P t , r 1 ( A ) a t t t a j j 2 j m .

Throughout this paper, we assume that N 1 and N 2 . Meanwhile, we also assume that tensor A satisfies: a j j j 0 , Λ j ( A ) 0 , j N .

Lemma 1

[17] If tensor A = ( a j 1 j 2 j m ) is strictly diagonally dominant, then A is an -tensor.

Lemma 2

[14] Let A = ( a j 1 j 2 j m ) be an mth-order n-dimensional complex tensor. If there is a positive diagonal matrix X such that A X m 1 is an -tensor, then A is an -tensor.

Lemma 3

[14] Let A = ( a j 1 j 2 j m ) be an mth-order n-dimensional complex tensor. If A is irreducible,

a j j j Λ j ( A ) , j N ,

and N 3 ( A ) , then A is an -tensor.

Lemma 4

[18] Let A = ( ( a j 1 j 2 j m ) ) be an mth-order n-dimensional complex tensor. If

  1. a j j j Λ j ( A ) , j N ;

  2. N 3 = { j N : a j j j > Λ j ( A ) } ;

  3. For j N 1 , there exists a nonzero elements chain from j to i such that i N 1 ,

then A is an -tensor.

Next, we give some new criteria for -tensors.

Theorem 2

Let A = ( a j 1 j 2 j m ) be an mth-order n-dimensional complex tensor. If A satisfies

(3) a j j j > μ w j j 2 j 3 j m N 0 m 1 a j j 2 j m + j 2 j 3 j m N 2 m 1 δ j j 2 j m = 0 a j j 2 j m w w j + h w j j 2 j 3 j m N 3 m 1 max t { j 2 , j 3 , , j m } P t , r 1 ( A ) a t t t a j j 2 j m , j N 2 ,

and there exists j 2 j 3 j m N m 1 N 1 m 1 for any j N 1 such that a j j 2 j 3 j m 0 , then A is an -tensor.

Proof

By 0 r 1 r < 1 , according to the definition of r , P i , r ( A ) , r 1 and P i , r 1 ( A ) , for any i N 3 ,

r 1 a i i i i 2 i m N 0 m 1 a i i 2 i m + i 2 i m N 2 m 1 a i i 2 i m + r 1 i 2 i m N 3 m 1 δ i i 2 i m = 0 max j { i 2 , i 3 , , i m } P j , r ( A ) a j j j a i i 2 i m ,

that is, P i , r 1 ( A ) r 1 a i i i ( i N 3 ) , so

0 < P i , r 1 ( A ) a i i i r 1 r < 1 , i N 3 .

By the definitions of w i and μ , for any i N 3 ,

μ i 2 i m N 0 m 1 a i i 2 i m + i 2 i m N 2 m 1 a i i 2 i m w P i , r 1 ( A ) i 2 i m N 3 m 1 δ i i 2 i m = 0 max j { i 2 , i 3 , , i m } P j , r 1 ( A ) a j j j a i i 2 i m μ P i , r 1 ( A ) r 1 i 2 i m N 3 m 1 δ i i 2 i m = 0 max j { i 2 , i 3 , , i m } P j , r ( A ) a j j j a i i 2 i m P i , r 1 ( A ) i 2 i m N 3 m 1 δ i i 2 i m = 0 max j { i 2 , i 3 , , i m } P j , r 1 ( A ) a j j j a i i 2 i m < 1 .

For i N 3 , from the definition of h , we have 0 < h < 1 , and

(4) h P i , r 1 ( A ) μ i 2 i m N 0 m 1 a i i 2 i m + i 2 i m N 2 m 1 a i i 2 i m w + h i 2 i m N 3 m 1 δ i i 2 i m = 0 max j { i 2 , i 3 , , i m } P j , r 1 ( A ) a j j j a i i 2 i m .

By inequality (3), for all i N 2 , we have

w i a i i i > μ i 2 i 3 i m N 0 m 1 a i i 2 i m + i 2 i 3 i m N 2 m 1 δ i i 2 i m = 0 a i i 2 i m w + h i 2 i 3 i m N 3 m 1 max j { i 2 , i 3 , , i m } P j , r 1 ( A ) a j j j a i i 2 i m .

Let

(5) R j 1 j 2 j 3 j m N 3 m 1 a j j 2 j m w j a j j j μ j 2 j 3 j m N 0 m 1 a j j 2 j m + j 2 j 3 j m N 2 m 1 δ j j 2 j m = 0 a j j 2 j m w + h j 2 j 3 j m N 3 m 1 max t { j 2 , j 3 , , j m } P t , r 1 ( A ) a t t t a j j 2 j m , j N .

If j 2 j 3 j m N 3 m 1 a j j 2 j m = 0 , we denote R j = + . By Equality (5), we have

R j > 0 ( j N 2 ) , 0 < h P t , r 1 ( A ) a t t t < 1 ( t N 3 ) .

Hence, there exists a positive number ε such that

0 < ε < min j N 2 R j + , max t N 3 h P t , r 1 ( A ) a t t t + ε < 1 .

Let the matrix D = diag ( d 1 , d 2 , , d n ) , denote = A D m 1 = ( b i 1 i 2 i m ) , where

d i = ( μ ) 1 m 1 , i N 1 , ( w i ) 1 m 1 , i N 2 , ε + h P i , r 1 ( A ) a i i i 1 m 1 , i N 3 .

For j N 1 , there exists j 2 j 3 j m N m 1 N 1 m 1 such that a j j 2 j 3 j m 0 , and for t N 3 , let ε > 0 satisfy 0 < ε + h P t , r 1 ( A ) a t t t < μ < 1 , we obtain

Λ j ( ) = μ j 2 j m N 0 m 1 δ j j 2 j m = 0 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m w + j 2 j m N 3 m 1 a j j 2 j m ε + h P j 2 , r 1 ( A ) a j 2 j 2 j 2 1 m 1 ε + h P j m , r 1 ( A ) a j m j m j m 1 m 1 μ j 2 j m N 0 m 1 δ j j 2 j m = 0 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m w + j 2 j m N 3 m 1 a j j 2 j m ε + max t { j 2 , j 3 , , j m } h P t , r 1 ( A ) a t t t < μ j 2 j m N 0 m 1 δ j j 2 j m = 0 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m + j 2 j m N 3 m 1 a j j 2 j m = μ Λ j ( A ) = μ a j j j = b j j j .

For j N 2 , by equality (5), we have

Λ j ( ) = μ j 2 j m N 0 m 1 a j j 2 j m + j 2 j m N 2 m 1 δ j j 2 j m = 0 a j j 2 j m w + j 2 j m N 3 m 1 a j j 2 j m ε + h P j 2 , r 1 ( A ) a j 2 j 2 j 2 1 m 1 ε + h P j m , r 1 ( A ) a j m j m j m 1 m 1 μ j 2 j m N 0 m 1 a j j 2 j m + j 2 j m N 2 m 1 δ j j 2 j m = 0 a j j 2 j m w + j 2 j m N 3 m 1 a j j 2 j m ε + max t { j 2 , j 3 , , j m } h P t , r 1 ( A ) a t t t = ε j 2 j 3 j m N 3 m 1 a j j 2 j m + μ j 2 j 3 j m N 0 m 1 a j j 2 j m + j 2 j 3 j m N 2 m 1 δ j j 2 j m = 0 a j j 2 j m w + h j 2 j 3 j m N 3 m 1 max t { j 2 , j 3 , , j m } P t , r 1 ( A ) a t t t a j j 2 j m < R j j 2 j 3 j m N 3 m 1 a j j 2 j m + μ j 2 j 3 j m N 0 m 1 a j j 2 j m + j 2 j 3 j m N 2 m 1 δ j j 2 j m = 0 a j j 2 j m w + h j 2 j 3 j m N 3 m 1 max t { j 2 , j 3 , , j m } P t , r 1 ( A ) a t t t a j j 2 j m = w j a j j j = b j j j .

Finally, for j N 3 , by inequality (4), we have

Λ j ( ) = μ j 2 j m N 0 m 1 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m w + j 2 j m N 3 m 1 δ j j 2 j m = 0 a j j 2 j m ε + h P j 2 , r 1 ( A ) a j 2 j 2 j 2 1 m 1 ε + h P j m , r 1 ( A ) a j m j m j m 1 m 1 μ j 2 j m N 0 m 1 a j j 2 j m + j 2 j m N 2 m 1 a j j 2 j m w + j 2 j m N 3 m 1 δ j j 2 j m = 0 a j j 2 j m ε + max t { j 2 , j 3 , , j m } h P t , r 1 ( A ) a t t t = ε j 2 j 3 j m N 3 m 1 δ j j 2 j m = 0 a j j 2 j m + μ i 2 i 3 j m N 0 m 1 a j j 2 j m + j 2 j 3 j m N 2 m 1 a j j 2 j m w + h j 2 j 3 j m N 3 m 1 δ j j 2 j m = 0 max t { j 2 , j 3 , , j m } P t , r 1 ( A ) a t t t a j j 2 j m ε j 2 j 3 j m N 3 m 1 δ j j 2 j m = 0 a j j 2 j m + h P j , r 1 ( A ) < ε a j j j + h P j , r 1 ( A ) = b j j j .

Therefore, b j j j > Λ j ( ) ( i N ) . So is an -tensor by Lemma 1. From Lemma 2, A is an -tensor.□

Remark 1

It is hard to theoretically give the comparison between our result and the known ones in [20]. The following numerical example illustrates that our proposed criteria is more effective to theirs in some cases. Moreover, we provide an example that satisfies our conditions in Theorem 2 but not those in Theorem 1 of [20].

Let

K ( A ) = j N 2 : a j j j > μ w j j 2 j 3 j m N 0 m 1 a j j 2 j m + j 2 j 3 j m N 2 m 1 δ j j 2 j m = 0 a j j 2 j m w w j + h w j j 2 j 3 j m N 3 m 1 max t { j 2 , j 3 , , j m } P t , r 1 ( A ) a t t t a j j 2 j m .

Theorem 3

Let A = ( a j 1 j 2 j m ) be an mth-order n-dimensional complex tensor. If A is irreducible,

(6) a j j j μ w j j 2 j 3 j m N 0 m 1 a j j 2 j m + j 2 j 3 j m N 2 m 1 δ j j 2 j m = 0 a j j 2 j m w w j + h w j j 2 j 3 j m N 3 m 1 max t { j 2 , j 3 , , j m } P t , r 1 ( A ) a t t t a j j 2 j m , j N 2 ,

and K ( A ) , then A is an -tensor.

Proof

Since A is irreducible, so

i 2 i 3 i m N m 1 N 3 m 1 a i i 2 i m > 0 , i N 3 .

By inequality (6), we obtain

(7) a i i i w i μ i 2 i 3 i m N 0 m 1 a i i 2 i m + i 2 i 3 i m N 2 m 1 δ i i 2 i m = 0 a i i 2 i m w + h i 2 i 3 i m N 3 m 1 max j { i 2 , i 3 , , i m } P j , r 1 ( A ) a j j j a i i 2 i m , i N 2 .

and at least one strict inequality in (7) holds.

Let the matrix D = diag ( d 1 , d 2 , , d n ) , denote = A D m 1 = ( b i 1 i 2 i m ) , where

d i = ( μ ) 1 m 1 , i N 1 , ( w i ) 1 m 1 , i N 2 , h P i , r 1 ( A ) a i i i 1 m 1 , i N 3 .

For i N 1 , by μ > h P i , r 1 ( A ) a i i i ( i N 3 ) , we have

Λ i ( ) = μ i 2 i m N 0 m 1 δ i i 2 i m = 0 a i i 2 i m + i 2 i m N 2 m 1 a i i 2 i m w + i 2 i m N 3 m 1 a i i 2 i m h P i 2 , r 1 ( A ) a i 2 i 2 i 2 1 m 1 h P i m , r 1 ( A ) a i m i m i m 1 m 1 μ i 2 i m N 0 m 1 δ i i 2 i m = 0 a i i 2 i m + i 2 i m N 2 m 1 a i i 2 i m w + i 2 i m N 3 m 1 max j { i 2 , i 3 , , i m } h P j , r 1 ( A ) a j j j a i i 2 i m < μ i 2 i m N 0 m 1 δ i i 2 i m = 0 a i i 2 i m + i 2 i m N 2 m 1 a i i 2 i m + i 2 i m N 3 m 1 a i i 2 i m = μ a i i i = b i i i .

For i N 2 , by inequality (7), we obtain

Λ i ( ) = μ i 2 i m N 0 m 1 a i i 2 i m + i 2 i m N 2 m 1 δ i i 2 i m = 0 a i i 2 i m w + i 2 i m N 3 m 1 a i i 2 i m h P i 2 , r 1 ( A ) a i 2 i 2 i 2 1 m 1 h P i m , r 1 ( A ) a i m i m i m 1 m 1 μ i 2 i m N 0 m 1 a i i 2 i m + i 2 i m N 2 m 1 δ i i 2 i m = 0 a i i 2 i m w + i 2 i m N 3 m 1 max j { i 2 , i 3 , , i m } h P j , r 1 ( A ) a j j j a i i 2 i m w i a i i i = b i i i .

Next, we consider i N 3 , by inequality (4), we obtain

Λ i ( ) = μ i 2 i m N 0 m 1 a i i 2 i m + i 2 i m N 2 m 1 a i i 2 i m w + i 2 i m N 3 m 1 δ i i 2 i m = 0 a i i 2 i m h P i 2 , r 1 ( A ) a i 2 i 2 i 2 1 m 1 h P i m , r 1 ( A ) a i m i m i m 1 m 1 μ i 2 i m N 0 m 1 a i i 2 i m + i 2 i m N 2 m 1 a i i 2 i m w + i 2 i m N 3 m 1 δ i i 2 i m = 0 max j { i 2 , i 3 , , i m } h P j , r 1 ( A ) a j j j a i i 2 i m h P i , r 1 ( A ) = b i i i .

Therefore, b i i i Λ i ( ) for all i N , and at least one strict inequality in (6) holds, that is, there exists an i 0 N 2 such that b i 0 i 0 i 0 > Λ i 0 ( ) .

Notice that A is irreducible and so is . Hence, we obtain that is an -tensor by Lemma 3 and A is an -tensor by Lemma 2.□

Theorem 4

Let A = ( a j 1 j 2 j m ) be an mth-order n-dimensional complex tensor. If A satisfies

(8) a j j j μ w j j 2 j 3 j m N 0 m 1 a j j 2 j m + j 2 j 3 j m N 2 m 1 δ j j 2 j m = 0 a j j 2 j m w w j + h w j j 2 j 3 j m N 3 m 1 max t { j 2 , j 3 , , j m } P t , r 1 ( A ) a t t t a j j 2 j m , j N 2

and if any j N K ( A ) , there is a nonzero elements chain from j to i such that i K ( A ) , then A is an -tensor.

Proof

Let matrix D = diag ( d 1 , d 2 , , d n ) , denote = A D m 1 = ( b j 1 j 2 j m ) , where

d j = ( μ ) 1 m 1 , j N 1 , ( w j ) 1 m 1 , j N 2 , h P j , r 1 ( A ) a j j j 1 m 1 , j N 3 .

By a similar proof to that of Theorem 3, we obtain that b j j j Λ j ( ) ( j N ) , and there exists at least an i 0 N 2 such that b i 0 i 0 i 0 > Λ i 0 ( ) .

On the other hand, if b j j j = Λ j ( ) , then j N K ( A ) , so there is a nonzero elements chain of A from j to i satisfying i K ( A ) . Hence, there is a nonzero elements chain of from j to i satisfying b i i i > Λ i ( ) .

Therefore, we obtain that is an -tensor by Lemma 4 and A is an -tensor by Lemma 2.□

Below, we present a numerical example to show that the criteria for -tensor is obtained by only using the conditions of Theorem 2 but not the ones in [20].

Example 1

Given a 3rd-order three-dimensional tensor A = ( a i j k ) as follows:

A = [ A ( 1 , : , : ) , A ( 2 , : , : ) , A ( 3 , : , : ) ] , A ( 1 , : , : ) = 12 1 0 1 6 0 1 1 15 , A ( 2 , : , : ) = 1 1 0 0 4 0 0 0 6 , A ( 3 , : , : ) = 1 0 0 0 1 0 0 0 16 .

Obviously,

a 111 = 12 , Λ 1 ( A ) = 25 , a 222 = 4 , Λ 2 ( A ) = 8 , a 333 = 16 , Λ 3 ( A ) = 2 ,

so N 1 = , N 2 = { 1 , 2 } , N 3 = { 3 } . By calculations, we have

w 1 = 25 25 + 12 = 25 37 , w 2 = 8 8 + 4 = 2 3 , w = 25 37 , r = 0 + 2 16 0 = 1 8 ,

P 3 , r ( A ) = 0 + 2 + 1 8 × 0 = 2 , r 1 = 0 + 2 16 1 8 × 0 = 1 8 , μ = 25 37 ,

P 3 , r 1 ( A ) = 0 + 2 + 1 8 × 1 8 × 0 = 2 , h = 25 37 × 0 + 25 37 × 2 2 1 8 × 0 = 25 37 .

When i = 1 , we obtain

μ i 2 i 3 N 0 2 a 1 i 2 i 3 + i 2 i 3 N 2 2 δ 1 i 2 i 3 = 0 a 1 i 2 i 3 w + h i 2 i 3 N 3 2 max j { i 2 , i 3 } P j , r 1 ( A ) a j j j a 1 i 2 i 3 = 25 37 × 2 + 25 37 × 8 + 25 37 × 1 8 × 15 = 2375 296 < 300 37 = 12 × 25 37 = a 111 × w 1 .

When i = 2 , we obtain

μ i 2 i 3 N 0 2 a 2 i 2 i 3 + i 2 i 3 N 2 2 δ 2 i 2 i 3 = 0 a 2 i 2 i 3 w + h i 2 i 3 N 3 2 max j { i 2 , i 3 } P j , r 1 ( A ) a j j j a 2 i 2 i 3 = 25 37 × 0 + 25 37 × 2 + 25 37 × 1 8 × 6 = 275 148 < 100 37 = 4 × 25 37 = a 111 × w 1 .

Hence, the conditions of Theorem 2 are satisfied, then A is an -tensor. However,

p 3 = 26 25 , s 3 = 13 200 , M 3 = 1 , F 3 = 26 25 , t 3 = 13 200 .

i 2 i 3 N 0 2 a 1 i 2 i 3 + i 2 i 3 N 2 2 δ 1 i 2 i 3 = 0 max j { i 2 , i 3 } { r j } a 1 i 2 i 3 + i 2 i 3 N 3 2 max j { i 2 , i 3 } { t j } a 1 i 2 i 3 = 2 + 13 25 × 8 + 13 200 × 15 = 1427 200 156 25 = a 111 r 1 .

So, A does not satisfy the conditions of Theorem 1.

3 An application

In this section, we provide some -tensor-based criteria for identifying the positive definite tensor. First, we recall the following lemma.

Lemma 5

[14] Let A = ( a j 1 j 2 j m ) be an mth-order n-dimensional even-order real symmetric tensor with a j j j > 0 for j N . If A is an -tensor, then A is positive definite.

By Theorems 24 and Lemma 5, we can prove easily the following result.

Theorem 5

Let A = ( a j 1 j 2 j m ) be an mth-order n-dimensional even-order real symmetric tensor a j j j > 0 for j N . If one of the following conditions is satisfied:

  1. All the conditions of Theorem 2;

  2. All the conditions of Theorem 3;

  3. All the conditions of Theorem 4,

then A is a positive definite tensor.

Example 2

Consider an 4th-degree homogeneous polynomial:

f ( x ) = A x 4 = 24 x 1 4 + 36 x 2 4 + 25 x 3 4 + 10 x 4 4 8 x 2 3 x 4 + 12 x 1 2 x 2 x 3 12 x 2 x 3 2 x 4 24 x 1 x 2 x 3 x 4 .

We can obtain an 4th-order four-dimensional real symmetric tensor A = ( a j 1 j 2 j 3 j 4 ) with entries

a 1111 = 24 , a 2222 = 36 , a 3333 = 25 , a 4444 = 10 , a 2224 = a 2242 = a 2422 = a 4222 = 2 , a 1123 = a 1132 = a 1213 = a 1312 = a 1231 = a 1321 = 1 , a 2113 = a 2131 = a 2311 = a 3112 = a 3121 = a 3211 = 1 , a 2334 = a 2343 = a 2433 = a 4233 = a 4323 = a 4332 = 1 , a 3234 = a 3243 = a 3324 = a 3342 = a 3423 = a 3432 = 1 , a 1234 = a 1243 = a 1324 = a 1342 = a 1423 = a 1432 = 1 , a 2134 = a 2143 = a 2314 = a 2341 = a 2413 = a 2431 = 1 , a 3124 = a 3142 = a 3214 = a 3241 = a 3412 = a 3421 = 1 , a 4123 = a 4132 = a 4213 = a 4231 = a 4312 = a 4321 = 1 ,

and other a i 1 i 2 i 3 i 4 = 0 . By calculations, we have

a 4444 = 10 < 11 = Λ 4 ( A )

and

a 1111 ( a 4444 Λ 4 ( A ) + a 4111 ) = 24 < 0 = Λ 4 ( A ) a 4111 .

Therefore, tensor A is not strictly diagonally dominate as defined in [21], or quasidoubly strictly diagonally dominant as defined in [22]. However, all the conditions of Theorem 2 can be satisfied.

Λ 1 ( A ) = 12 , Λ 2 ( A ) = 18 , Λ 3 ( A ) = 15 , Λ 4 ( A ) = 11 ,

so N 1 = , N 2 = { 4 } , N 3 = { 1 , 2 , 3 } . By calculations, we have

w 4 = 11 21 = w , r = 6 11 , P 1 , r ( A ) = 102 11 , P 2 , r ( A ) = 183 11 , P 3 , r ( A ) = 150 11 , r 1 = 132 257 , μ = 11 21 , P 1 , r 1 ( A ) = 1974 257 , P 2 , r 1 ( A ) = 4071 257 , P 3 , r 1 ( A ) = 3300 257 , h = 2827 4137 .

When i = 4 , we obtain

μ i 2 i 3 i 4 N 0 3 a 4 i 2 i 3 i 4 + i 2 i 3 i 4 N 2 3 δ 4 i 2 i 3 i 4 = 0 a 4 i 2 i 3 i 4 w + h i 2 i 3 i 4 N 3 3 max j { i 2 , i 3 , i 4 } P j , r 1 ( A ) a j j j j a 4 i 2 i 3 i 4 = 11 21 × 0 + 11 21 × 0 + 2827 4137 × 132 257 × 11 = 5324 1379 < 121 21 = 11 × 11 21 = a 4444 × w 4 .

So, A is positive definite by Theorem 5, that is, f ( x ) is positive definite.

4 Conclusions

In this paper, we introduced some new criterions for identifying -tensors, which only depended on the elements of tensor. They were well used to check the positive definiteness of an even-order homogeneous polynomial form f ( x ) A x m . We also verified the effectiveness of new criteria by several numerical examples.

Acknowledgements

The authors would like to express their cordial gratitude to the referee for valuable comments which improved the paper.

  1. Funding information: This work was supported by the Foundation of Science and Technology Department of Guizhou Province (20191161, 20181079), the Foundation of Education Department of Guizhou Province (2018143) and the Research Foundation of Guizhou Minzu University (2019YB08).

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

References

[1] K. C. Chang, K. Pearson, and T. Zhang, Perron-Frobenius theorem for nonnegative tensors, Commun. Math. Sci. 6 (2008), no. 2, 507–520. 10.4310/CMS.2008.v6.n2.a12Search in Google Scholar

[2] L. D. Lathauwer, B. D. Moor, and J. Vandewalle, On the best rank-1 and rank-(R1, R2, RN) approximation of higer-order tensors, SIAM J. Matrix Anal. Appl. 21 (2000), no. 4, 1324–1342. 10.1137/S0895479898346995Search in Google Scholar

[3] Y. Liu, G. L. Zhou, and N. F. Ibrahim, An always convergent algorithm for the largest eigenvalue of an irreducible nonnegative tensor, J. Comput. Appl. Math. 235 (2010), no. 1, 286–292. 10.1016/j.cam.2010.06.002Search in Google Scholar

[4] M. Ng, L. Q. Qi, and G. L. Zhou, Finding the largest eigenvalue of a nonnegative tensor, SIAM J. Matrix Anal. Appl. 31 (2009), 1090–1099. 10.1137/09074838XSearch in Google Scholar

[5] L. Q. Qi, Eigenvalues of a real supersymetric tensor, J. Symbolic Comput. 40 (2005), no. 6, 1302–1324. 10.1016/j.jsc.2005.05.007Search in Google Scholar

[6] L. H. Lim, Singular values and eigenvalues of tensors: a variational approach, in: CAMSAP05: Proceeding of the IEEE International Workshop on Computational Advances in Multi Sensor Adaptive Processing, 2005, p. 129–132. Search in Google Scholar

[7] E. Kofidis and P. A. Regalia, On the best rank-1 approximation of higher-order supersymmetric tensors, SIAM J. Matrix Anal. Appl. 23 (2002), no. 3, 863–884. 10.1137/S0895479801387413Search in Google Scholar

[8] Y. N. Yang and Q. Z. Yang, Further results for Perron-Frobenius Theorem for nonnegative tensors, SIAM J. Matrix Anal. Appl. 31 (2010), 2517–2530. 10.1137/090778766Search in Google Scholar

[9] Q. Ni, L. Q. Qi, and F. Wang, An eigenvalue method for testing positive definiteness of a multivariate form, IEEE Trans. Automat. Control 53 (2008), 1096–1107. 10.1109/TAC.2008.923679Search in Google Scholar

[10] H. B. Chen and L. Q. Qi, Positive deniteness and semi-deniteness of even order symmetric Cauchy tensors, J. Ind. Manag. Optim. 11 (2015), 1263–1274. 10.3934/jimo.2015.11.1263Search in Google Scholar

[11] Y. S. Song and L. Q. Qi, Necessary and sufficient conditions for copositive tensors, Linear Multilinear Algebra 63 (2015), no. 1, 120–131. 10.1080/03081087.2013.851198Search in Google Scholar

[12] K. L. Zhang and Y. J. Wang, An ℋ-tensor-based iterative scheme for identifying the positive deniteness of multivariate homogeneous forms, J. Comput. Appl. Math. 305 (2016), 1–10. 10.1016/j.cam.2016.03.025Search in Google Scholar

[13] N. K. Bose and A. R. Modaress, General procedure for multivariable polynomial positivity with control applications, IEEE Trans. Automat. Control 21 (1976), no. 5, 696–701. 10.1109/TAC.1976.1101356Search in Google Scholar

[14] C. Q. Li, F. Wang, J. X. Zhao, Y. Zhu, and Y. Li, Criterions for the positive definiteness of real supersymmetric tensors, J. Comput. Appl. Math. 255 (2014), 1–14. 10.1016/j.cam.2013.04.022Search in Google Scholar

[15] L. P. Zhang, L. Q. Qi, and G. L. Zhou, ℳ-tensors and some applications, SIAM J. Matrix Anal. Appl. 35 (2014), 437–542. 10.1137/130915339Search in Google Scholar

[16] M. R. Kannana, N. S. Mondererb, and A. Bermana, Some properties of strong ℋ-tensors and general ℋ-tensors, Linear Algebra Appl. 476 (2015), 42–55. 10.1016/j.laa.2015.02.034Search in Google Scholar

[17] W. Y. Ding, L. Q. Qi, and Y. M. Wei, ℳ-tensors and nonsingular ℳ-tensors, Linear Algebra Appl. 439 (2013), no. 10, 3264–3278. 10.1016/j.laa.2013.08.038Search in Google Scholar

[18] F. Wang and D. S. Sun, New criteria for ℋ-tensors and an application, J. Inequal. Appl. 96 (2016), 1–12. 10.1186/s13660-016-1041-0Search in Google Scholar

[19] F. Wang, D. S. Sun, J. X. Zhao, and C. Q. Li, New practical criteria for ℋ-tensors and its application, Linear Multilinear Algebra 65 (2017), no. 2, 269–283. 10.1080/03081087.2016.1183558Search in Google Scholar

[20] F. Wang, D. S. Sun, and Y. M. Xu, Some criteria for identifying ℋ-tensors and its applications, Calcolo 56 (2019), 19. 10.1007/s10092-019-0316-5Search in Google Scholar

[21] L. Q. Qi and Y. S. Song, An even-order symmetric ℬ-tensor is positive definite, Linear Algebra Appl. 457 (2014), 303–312. 10.1016/j.laa.2014.05.026Search in Google Scholar

[22] C. Q. Li and Y. T. Li, Double ℬ-tensors and quasi-double ℬ-tensors, Linear Algebra Appl. 466 (2015), 343–356. 10.1016/j.laa.2014.10.027Search in Google Scholar

Received: 2021-06-26
Revised: 2022-01-24
Accepted: 2022-04-08
Published Online: 2022-06-25

© 2022 Dongjian Bai and Feng Wang, published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. A random von Neumann theorem for uniformly distributed sequences of partitions
  3. Note on structural properties of graphs
  4. Mean-field formulation for mean-variance asset-liability management with cash flow under an uncertain exit time
  5. The family of random attractors for nonautonomous stochastic higher-order Kirchhoff equations with variable coefficients
  6. The intersection graph of graded submodules of a graded module
  7. Isoperimetric and Brunn-Minkowski inequalities for the (p, q)-mixed geominimal surface areas
  8. On second-order fuzzy discrete population model
  9. On certain functional equation in prime rings
  10. General complex Lp projection bodies and complex Lp mixed projection bodies
  11. Some results on the total proper k-connection number
  12. The stability with general decay rate of hybrid stochastic fractional differential equations driven by Lévy noise with impulsive effects
  13. Well posedness of magnetohydrodynamic equations in 3D mixed-norm Lebesgue space
  14. Strong convergence of a self-adaptive inertial Tseng's extragradient method for pseudomonotone variational inequalities and fixed point problems
  15. Generic uniqueness of saddle point for two-person zero-sum differential games
  16. Relational representations of algebraic lattices and their applications
  17. Explicit construction of mock modular forms from weakly holomorphic Hecke eigenforms
  18. The equivalent condition of G-asymptotic tracking property and G-Lipschitz tracking property
  19. Arithmetic convolution sums derived from eta quotients related to divisors of 6
  20. Dynamical behaviors of a k-order fuzzy difference equation
  21. The transfer ideal under the action of orthogonal group in modular case
  22. The multinomial convolution sum of a generalized divisor function
  23. Extensions of Gronwall-Bellman type integral inequalities with two independent variables
  24. Unicity of meromorphic functions concerning differences and small functions
  25. Solutions to problems about potentially Ks,t-bigraphic pair
  26. Monotonicity of solutions for fractional p-equations with a gradient term
  27. Data smoothing with applications to edge detection
  28. An ℋ-tensor-based criteria for testing the positive definiteness of multivariate homogeneous forms
  29. Characterizations of *-antiderivable mappings on operator algebras
  30. Initial-boundary value problem of fifth-order Korteweg-de Vries equation posed on half line with nonlinear boundary values
  31. On a more accurate half-discrete Hilbert-type inequality involving hyperbolic functions
  32. On split twisted inner derivation triple systems with no restrictions on their 0-root spaces
  33. Geometry of conformal η-Ricci solitons and conformal η-Ricci almost solitons on paracontact geometry
  34. Bifurcation and chaos in a discrete predator-prey system of Leslie type with Michaelis-Menten prey harvesting
  35. A posteriori error estimates of characteristic mixed finite elements for convection-diffusion control problems
  36. Dynamical analysis of a Lotka Volterra commensalism model with additive Allee effect
  37. An efficient finite element method based on dimension reduction scheme for a fourth-order Steklov eigenvalue problem
  38. Connectivity with respect to α-discrete closure operators
  39. Khasminskii-type theorem for a class of stochastic functional differential equations
  40. On some new Hermite-Hadamard and Ostrowski type inequalities for s-convex functions in (p, q)-calculus with applications
  41. New properties for the Ramanujan R-function
  42. Shooting method in the application of boundary value problems for differential equations with sign-changing weight function
  43. Ground state solution for some new Kirchhoff-type equations with Hartree-type nonlinearities and critical or supercritical growth
  44. Existence and uniqueness of solutions for the stochastic Volterra-Levin equation with variable delays
  45. Ambrosetti-Prodi-type results for a class of difference equations with nonlinearities indefinite in sign
  46. Research of cooperation strategy of government-enterprise digital transformation based on differential game
  47. Malmquist-type theorems on some complex differential-difference equations
  48. Disjoint diskcyclicity of weighted shifts
  49. Construction of special soliton solutions to the stochastic Riccati equation
  50. Remarks on the generalized interpolative contractions and some fixed-point theorems with application
  51. Analysis of a deteriorating system with delayed repair and unreliable repair equipment
  52. On the critical fractional Schrödinger-Kirchhoff-Poisson equations with electromagnetic fields
  53. The exact solutions of generalized Davey-Stewartson equations with arbitrary power nonlinearities using the dynamical system and the first integral methods
  54. Regularity of models associated with Markov jump processes
  55. Multiplicity solutions for a class of p-Laplacian fractional differential equations via variational methods
  56. Minimal period problem for second-order Hamiltonian systems with asymptotically linear nonlinearities
  57. Convergence rate of the modified Levenberg-Marquardt method under Hölderian local error bound
  58. Non-binary quantum codes from constacyclic codes over 𝔽q[u1, u2,…,uk]/⟨ui3 = ui, uiuj = ujui
  59. On the general position number of two classes of graphs
  60. A posteriori regularization method for the two-dimensional inverse heat conduction problem
  61. Orbital stability and Zhukovskiǐ quasi-stability in impulsive dynamical systems
  62. Approximations related to the complete p-elliptic integrals
  63. A note on commutators of strongly singular Calderón-Zygmund operators
  64. Generalized Munn rings
  65. Double domination in maximal outerplanar graphs
  66. Existence and uniqueness of solutions to the norm minimum problem on digraphs
  67. On the p-integrable trajectories of the nonlinear control system described by the Urysohn-type integral equation
  68. Robust estimation for varying coefficient partially functional linear regression models based on exponential squared loss function
  69. Hessian equations of Krylov type on compact Hermitian manifolds
  70. Class fields generated by coordinates of elliptic curves
  71. The lattice of (2, 1)-congruences on a left restriction semigroup
  72. A numerical solution of problem for essentially loaded differential equations with an integro-multipoint condition
  73. On stochastic accelerated gradient with convergence rate
  74. Displacement structure of the DMP inverse
  75. Dependence of eigenvalues of Sturm-Liouville problems on time scales with eigenparameter-dependent boundary conditions
  76. Existence of positive solutions of discrete third-order three-point BVP with sign-changing Green's function
  77. Some new fixed point theorems for nonexpansive-type mappings in geodesic spaces
  78. Generalized 4-connectivity of hierarchical star networks
  79. Spectra and reticulation of semihoops
  80. Stein-Weiss inequality for local mixed radial-angular Morrey spaces
  81. Eigenvalues of transition weight matrix for a family of weighted networks
  82. A modified Tikhonov regularization for unknown source in space fractional diffusion equation
  83. Modular forms of half-integral weight on Γ0(4) with few nonvanishing coefficients modulo
  84. Some estimates for commutators of bilinear pseudo-differential operators
  85. Extension of isometries in real Hilbert spaces
  86. Existence of positive periodic solutions for first-order nonlinear differential equations with multiple time-varying delays
  87. B-Fredholm elements in primitive C*-algebras
  88. Unique solvability for an inverse problem of a nonlinear parabolic PDE with nonlocal integral overdetermination condition
  89. An algebraic semigroup method for discovering maximal frequent itemsets
  90. Class-preserving Coleman automorphisms of some classes of finite groups
  91. Exponential stability of traveling waves for a nonlocal dispersal SIR model with delay
  92. Existence and multiplicity of solutions for second-order Dirichlet problems with nonlinear impulses
  93. The transitivity of primary conjugacy in regular ω-semigroups
  94. Stability estimation of some Markov controlled processes
  95. On nonnil-coherent modules and nonnil-Noetherian modules
  96. N-Tuples of weighted noncommutative Orlicz space and some geometrical properties
  97. The dimension-free estimate for the truncated maximal operator
  98. A human error risk priority number calculation methodology using fuzzy and TOPSIS grey
  99. Compact mappings and s-mappings at subsets
  100. The structural properties of the Gompertz-two-parameter-Lindley distribution and associated inference
  101. A monotone iteration for a nonlinear Euler-Bernoulli beam equation with indefinite weight and Neumann boundary conditions
  102. Delta waves of the isentropic relativistic Euler system coupled with an advection equation for Chaplygin gas
  103. Multiplicity and minimality of periodic solutions to fourth-order super-quadratic difference systems
  104. On the reciprocal sum of the fourth power of Fibonacci numbers
  105. Averaging principle for two-time-scale stochastic differential equations with correlated noise
  106. Phragmén-Lindelöf alternative results and structural stability for Brinkman fluid in porous media in a semi-infinite cylinder
  107. Study on r-truncated degenerate Stirling numbers of the second kind
  108. On 7-valent symmetric graphs of order 2pq and 11-valent symmetric graphs of order 4pq
  109. Some new characterizations of finite p-nilpotent groups
  110. A Billingsley type theorem for Bowen topological entropy of nonautonomous dynamical systems
  111. F4 and PSp (8, ℂ)-Higgs pairs understood as fixed points of the moduli space of E6-Higgs bundles over a compact Riemann surface
  112. On modules related to McCoy modules
  113. On generalized extragradient implicit method for systems of variational inequalities with constraints of variational inclusion and fixed point problems
  114. Solvability for a nonlocal dispersal model governed by time and space integrals
  115. Finite groups whose maximal subgroups of even order are MSN-groups
  116. Symmetric results of a Hénon-type elliptic system with coupled linear part
  117. On the connection between Sp-almost periodic functions defined on time scales and ℝ
  118. On a class of Harada rings
  119. On regular subgroup functors of finite groups
  120. Fast iterative solutions of Riccati and Lyapunov equations
  121. Weak measure expansivity of C2 dynamics
  122. Admissible congruences on type B semigroups
  123. Generalized fractional Hermite-Hadamard type inclusions for co-ordinated convex interval-valued functions
  124. Inverse eigenvalue problems for rank one perturbations of the Sturm-Liouville operator
  125. Data transmission mechanism of vehicle networking based on fuzzy comprehensive evaluation
  126. Dual uniformities in function spaces over uniform continuity
  127. Review Article
  128. On Hahn-Banach theorem and some of its applications
  129. Rapid Communication
  130. Discussion of foundation of mathematics and quantum theory
  131. Special Issue on Boundary Value Problems and their Applications on Biosciences and Engineering (Part II)
  132. A study of minimax shrinkage estimators dominating the James-Stein estimator under the balanced loss function
  133. Representations by degenerate Daehee polynomials
  134. Multilevel MC method for weak approximation of stochastic differential equation with the exact coupling scheme
  135. Multiple periodic solutions for discrete boundary value problem involving the mean curvature operator
  136. Special Issue on Evolution Equations, Theory and Applications (Part II)
  137. Coupled measure of noncompactness and functional integral equations
  138. Existence results for neutral evolution equations with nonlocal conditions and delay via fractional operator
  139. Global weak solution of 3D-NSE with exponential damping
  140. Special Issue on Fractional Problems with Variable-Order or Variable Exponents (Part I)
  141. Ground state solutions of nonlinear Schrödinger equations involving the fractional p-Laplacian and potential wells
  142. A class of p1(x, ⋅) & p2(x, ⋅)-fractional Kirchhoff-type problem with variable s(x, ⋅)-order and without the Ambrosetti-Rabinowitz condition in ℝN
  143. Jensen-type inequalities for m-convex functions
  144. Special Issue on Problems, Methods and Applications of Nonlinear Analysis (Part III)
  145. The influence of the noise on the exact solutions of a Kuramoto-Sivashinsky equation
  146. Basic inequalities for statistical submanifolds in Golden-like statistical manifolds
  147. Global existence and blow up of the solution for nonlinear Klein-Gordon equation with variable coefficient nonlinear source term
  148. Hopf bifurcation and Turing instability in a diffusive predator-prey model with hunting cooperation
  149. Efficient fixed-point iteration for generalized nonexpansive mappings and its stability in Banach spaces
Downloaded on 30.1.2026 from https://www.degruyterbrill.com/document/doi/10.1515/math-2022-0043/html?lang=en
Scroll to top button