Research article

Some new bounds on the spectral radius of nonnegative matrices

  • Received: 01 August 2019 Accepted: 02 December 2019 Published: 19 December 2019
  • MSC : 15A18, 15A42, 05C50

  • In this paper, we determine some new bounds for the spectral radius of a nonnegative matrix with respect to a new defined quantity, which can be considered as an average of average 2-row sums. The new formulas extend previous results using the row sums and the average 2-row sums of a nonnegative matrix. We also characterize the equality cases of the bounds if the matrix is irreducible and we provide illustrative examples comparing with the existing bounds.

    Citation: Maria Adam, Dimitra Aggeli, Aikaterini Aretaki. Some new bounds on the spectral radius of nonnegative matrices[J]. AIMS Mathematics, 2020, 5(1): 701-716. doi: 10.3934/math.2020047

    Related Papers:

    [1] Mengting Gan . Global error bounds for the extended vertical LCP of $ \sum-SDD $ matrices. AIMS Mathematics, 2024, 9(9): 24326-24335. doi: 10.3934/math.20241183
    [2] Sumaira Hafeez, Rashid Farooq . On generalized inverse sum indeg index and energy of graphs. AIMS Mathematics, 2020, 5(3): 2388-2411. doi: 10.3934/math.2020158
    [3] Zhen Lin . On the sum of the largest $ A_{\alpha} $-eigenvalues of graphs. AIMS Mathematics, 2022, 7(8): 15064-15074. doi: 10.3934/math.2022825
    [4] Qin Zhong, Na Li . Generalized Perron complements in diagonally dominant matrices. AIMS Mathematics, 2024, 9(12): 33879-33890. doi: 10.3934/math.20241616
    [5] Qin Zhong, Ling Li . Notes on the generalized Perron complements involving inverse $ {{N}_{0}} $-matrices. AIMS Mathematics, 2024, 9(8): 22130-22145. doi: 10.3934/math.20241076
    [6] Wafaa Fakieh, Zakeiah Alkhamisi, Hanaa Alashwali . On the $ A_{\alpha^-} $-spectra of graphs and the relation between $ A_{\alpha} $- and $ A_{\alpha^-} $-spectra. AIMS Mathematics, 2024, 9(2): 4587-4603. doi: 10.3934/math.2024221
    [7] Chuang Lv, Lihua You, Yufei Huang . A general result on the spectral radii of nonnegative k-uniform tensors. AIMS Mathematics, 2020, 5(3): 1799-1819. doi: 10.3934/math.2020121
    [8] Qin Zhong . Some new inequalities for nonnegative matrices involving Schur product. AIMS Mathematics, 2023, 8(12): 29667-29680. doi: 10.3934/math.20231518
    [9] Huan Xu, Tao Yu, Fawaz E. Alsaadi, Madini Obad Alassafi, Guidong Yu, Jinde Cao . Some spectral sufficient conditions for a graph being pancyclic. AIMS Mathematics, 2020, 5(6): 5389-5401. doi: 10.3934/math.2020346
    [10] Panpan Liu, Haifeng Sang, Min Li, Guorui Huang, He Niu . New criteria for nonsingular $ H $-matrices. AIMS Mathematics, 2023, 8(8): 17484-17502. doi: 10.3934/math.2023893
  • In this paper, we determine some new bounds for the spectral radius of a nonnegative matrix with respect to a new defined quantity, which can be considered as an average of average 2-row sums. The new formulas extend previous results using the row sums and the average 2-row sums of a nonnegative matrix. We also characterize the equality cases of the bounds if the matrix is irreducible and we provide illustrative examples comparing with the existing bounds.


    Let Mn(R) be the algebra of n×n real matrices. We refer to A=(aij)ni,j=1Mn(R) as a nonnegative or a positive matrix, when each aij0 or aij>0, respectively, denoted by writing A0 or A>0. The matrix AMn(R) is called irreducible if and only if (I+A)n1>0. We also define the spectral radius of A by

    ρ(A)=max{|λ|:λσ(A)},

    where σ(A) denotes the spectrum of A, that is, the set of eigenvalues of A.

    The spectral radius of a nonnegative matrix has been studied extensively in the mathematical fields of dynamical systems and graph theory (see, e.g. [1,2,3,4,6,9,10,11,12,13]). For example, the topological entropy, one of the main invariants of a topological dynamical system telling us how chaotic the system is, can often be computed as a logarithm of the spectral radius of a certain nonnegative matrix [13]. Also, from the perspective of graph theory, several combinatorial properties of simple undirected graphs have been interpreted via the spectral radius of the signless Laplacian matrix derived accordingly from adjacency and degree matrices of a graph (see e.g. [4,6,9,11] and the references therein).

    The problem of bounding the largest eigenvalue in modulus of a nonnegative matrix has attracted the interest of many researchers [2,4,5,6,7,8,9,10,11]. The celebrated Perron-Frobenius theory investigated the existence of positive eigenvalues of nonnegative matrices [5,8]. In particular, Perron proved that ρ(A) is a positive and simple eigenvalue of A>0 [8,Theorem 8.2.11], and Frobenius generalized Perron's statement to nonnegative and irreducible matrices [8,Theorem 8.4.4]. Moreover, apart from the familiar power method, other numerical algorithms have been also constructed and implemented for locating the spectrum of matrices as in [3,12]. However, the proposed methods are only valid on the limited class of diagonalizable matrices or on irreducible nonnegative matrices, whose spectral radius has been proved to be a simple eigenvalue.

    In this paper, we prove some new formulas bounding the spectral radius of any nonnegative matrix. They are expressed in terms of the elements of the matrix and they extend previous results, [4,11]. In the remainder, we give the necessary notation required for our results. For 1in, the quantities

    ri(A)=nj=1aijandMi(A)=nj=1aijrj(A)

    are known as the i-th row and the i-th 2-row sums of A, respectively, and for ri(A)>0 their ratio

    mi(A)=Mi(A)ri(A)

    is called the i-th average 2-row sum of A, (see [9,11]). Motivated by the expression of mi(A), we define a new quantity as the ratio

    wi(A)=nj=1aijmj(A)mi(A), (1.1)

    which can be seen as a further development of the quantity mi(A) and it will be used to compute new bounds for the spectral radius of a nonnegative matrix. Moreover, the following specific entries of matrix A will be used: its largest diagonal and off-diagonal elements,

    M=max1in{aii}andN=max1i,jnij{aij},

    respectively, as well as, its smallest diagonal and off-diagonal elements,

    S=min1in{aii}andT=min1i,jnij{aij},

    respectively.

    In this article we take into account the new quantities wi(A), i=1,,n, defined in (1.1), which can be considered as an average of averages 2-row, to present some new bounds for the spectral radius of a nonnegative matrix. Analytically, in section 2 we prove sharp upper bounds of ρ(A) with respect to wi(A), while in section 3 we turn our attention to a sharp lower bound of ρ(A). In both sections we characterize the equality cases of the bounds if the matrix is irreducible. Illustrative numerical examples are also provided testing our results and comparing with known bounds.

    In this section, we investigate some upper bounds for the spectral radius ρ(A) of a nonnegative matrix extending established results. Motivated by [4,9,11] and adopting the techniques used therein, we obtain a new expression for a sharp upper bound of ρ(A). In addition, we compare our findings with the ones presented in [4,11] providing illustrative examples.

    The next lemmas demonstrate well-established bounds for ρ(A), and since they are used in our arguments, they are stated here for the sake of completeness.

    Lemma 1. ([5]) Let AMn(R), A0 with i-th row sum ri(A), i=1,,n and largest diagonal element M. Then

    ρ(A)M,

    and

    min1in{ri(A)}ρ(A)max1in{ri(A)}.

    Moreover, if A is irreducible, then either equality holds if and only if r1(A)==rn(A).

    Lemma 2. ([8,Theorem 8.1.26,Corollary 8.1.31]) Let AMn(R), A0. Then for any positive vector xRn we have

    min1in{1xinj=1aijxj}ρ(A)max1in{1xinj=1aijxj}.

    If A is also irreducible, then

    ρ(A)=maxx>0min1in{1xinj=1aijxj}=minx>0max1in{1xinj=1aijxj}.

    Combining the definition of wi(A) in (1.1) with Lemma 2 by taking a vector x with components mi(A)>0, 1in, we immediately obtain the following proposition, which locates the spectral radius of a nonnegative matrix among the maximum and minimum values of wi(A).

    Proposition 3. Let AMn(R), A0. Then

    min1in{wi(A)}ρ(A)max1in{wi(A)}.

    If A is also irreducible, then either equality holds if and only if w1(A)==wn(A).

    Next we state and prove the main result in this section.

    Theorem 4. Let AMn(R), A0 with row sums ri(A)>0 and average 2-row sums mi(A)>0 for all i=1,,n, and wi(A) as defined in (1.1) such that w1(A)w2(A)wn(A)>0. Let M,N be the largest diagonal and off-diagonal elements of A, respectively, with N>0. Denoting by bm=max{mj(A)mi(A):1i,jn,ij}, consider

    Ψ=12(w(A)+MNbm+Δ),=1,,n, (2.1)

    where

    Δ=(w(A)M+Nbm)2+4Nbm1j=1(wj(A)w(A)). (2.2)

    Then

    ρ(A)min{Ψ:1n}. (2.3)

    Proof. To simplify the exposition of our calculations, we let mi=mi(A), and wi=wi(A) for 1in.

    Consider =1. Due to Lemma 1, Proposition 3 and our assumptions,

    Mρ(A)max1in{wi}w1,

    which yield that

    w1M0w1M+NbmNbm>0.

    Then the quantities in (2.1), (2.2) give

    Ψ1=12(w1+MNbm+(w1M+Nbm)2)=12(w1+MNbm+w1M+Nbm)=w1.

    Then ρ(A)Ψ1.

    Consider 2n. Let U=diag(m1x1,,m1x1,m,,mn) be an n×n diagonal matrix, where xj1 is a variable to be determined later for 1j1 and let B=U1AU. Due to similarity, A and B have the same eigenvalues, hence ρ(A)=ρ(B).

    For 1i1, we derive

    ri(B)=ri(U1AU)=1xi(1j=1aijmjmixj+nj=aijmjmi)=1xi(1j=1jiaijmjmixj+aiixi+nj=1aijmjmi1j=1aijmjmi)=1xi(nj=1aijmjmi+1j=1jiaijmjmixj1j=1jiaijmjmi+aiixiaii)=1xi(nj=1aijmjmi+1j=1jiaijmjmi(xj1)+aii(xi1)). (2.4)

    Obviously, aiiM, aijN, and mjmibm for 1j1 and ji. Combining these inequalities with the definition of wi in (1.1), the equality (2.4) is formulated as

    ri(B)1xi(wi+Nbm1j=1ji(xj1)+M(xi1)),i=1,,1. (2.5)

    For in, we derive

    ri(B)=ri(U1AU)=1j=1aijmjmixj+nj=aijmjmi=1j=1aijmjmixj+nj=1aijmjmi1j=1aijmjmi=nj=1aijmjmi+1j=1aijmjmi(xj1) (2.6)
    w+Nbm1j=1(xj1). (2.7)

    At this point, in order to construct the variable xj for j=1,,1 and =2,,n, we consider the real roots of the quadratic equations

    Ψ2(w+MNbm)Ψ+w(MNbm)Nbm1j=1(wjw)=0. (2.8)

    In particular, the trinomials in (2.8) have discriminant

    Δ(w+MNbm)24(w(MNbm)Nbm1j=1(wjw))=(wM+Nbm)2+4Nbm1j=1(wjw).

    Due to the hypotheses that N>0, bm>0 and {wi}ni=1 is a decreasing sequence of positive numbers, the discriminant Δ is a positive number for all =2,,n, which yields that the quadratic equations in (2.8) have two distinct real roots with a positive one

    Ψ=12(w+MNbm+Δ),=2,,n. (2.9)

    Now, for 1j1, we consider

    xj=1+wjwΨM+Nbmxj1=wjwΨM+Nbm, (2.10)

    where Ψ is given by (2.9). If 1j=1(wjw)>0, it is clear by relation (2.9) that

    Ψ>12(w+MNbm+|wM+Nbm|)12(w+MNbm(wM+Nbm))=MNbm,

    otherwise, w1==wMwM+Nbm>0 and (2.9) yields

    Ψ=12(w+MNbm+|wM+Nbm|)>12(w+MNbm(wM+Nbm))=MNbm.

    Both cases ensure xj10 in (2.10).

    Additionally, by the expression (2.8), we may write

    Nbm1j=1(wjw)=Ψ2(w+MNbm)Ψ+w(MNbm)=Ψ(ΨwM+Nbm)+w(MNbm)=Ψ(ΨM+Nbm)Ψw+w(MNbm)=Ψ(ΨM+Nbm)w(ΨM+Nbm)=(ΨM+Nbm)(Ψw). (2.11)

    Overall, we take xj10, j=1,,1 from (2.10) and Nbm1j=1(wjw) from (2.11) and substitute them into the inequality (2.5). Hence, for 1i1, we obtain

    ri(B)1xi(wi+Nbm1j=1ji(xj1)+M(xi1))=1xi(wi+Nbm1j=1(xj1)+M(xi1)Nbm(xi1))=1xi(wi+Nbm1j=1(xj1)+(MNbm)(xi1))=1xi(wi+Nbm1j=1wjwΨM+Nbm+(MNbm)wiwΨM+Nbm)=1xi(wi+(ΨM+Nbm)(Ψw)ΨM+Nbm+(MNbm)wiwΨM+Nbm)=wi(ΨM+Nbm)+(ΨM+Nbm)(Ψw)+(MNbm)(wiw)xi(ΨM+Nbm)=(ΨM+Nbm)Ψ+(wiw)(ΨM+Nbm+MNbm)xi(ΨM+Nbm)=(ΨM+Nbm)Ψ+(wiw)ΨΨM+Nbm+wiwΨM+Nbm(ΨM+Nbm)=(ΨM+Nbm+wiw)ΨΨM+Nbm+wiw=Ψ, (2.12)

    where =2,,n.

    Similarly, for in and 1j1 we substitute the relations (2.10) and (2.11) into the inequality (2.7), which can be written as

    ri(B)wi+Nbm1j=1(xj1)=wi+NbmΨM+Nbm1j=1(wjw)=wi+(ΨM+Nbm)(Ψw)ΨM+Nbm=wi+Ψww+Ψw=Ψ. (2.13)

    Thus, for all 2n and 1in the inequalities (2.12) and (2.13) verify

    ri(B)Ψ,

    and by Lemma 1,

    ρ(A)=ρ(B)max1in{ri(B)}Ψ. (2.14)

    The validity of (2.3) follows readily from (2.14).

    Remark 1. Obviously, inequalities (2.5) and (2.7) are stated as equalities for the special values of M,N,bm,xi, as the following cases indicate.

    (ⅰ) For 1i1, (2.5) is given by equality if and only if (a) and (b) hold:

    (a) xi=1 or aii=M, when xi>1,

    (b) xj=1 or aij=N and mjmi=bm, when xj>1 with 1j1 and ji.

    (ⅱ) For in, (2.7) is given by equality if and only if (c) and (d) hold:

    (c) xj=1 or aij=N and mjmi=bm, when xj>1 with 1j1,

    (d) wi=w.

    Using Remark 1 and Proposition 3 for a nonnegative and irreducible matrix, a special formulation of the spectral radius of the matrix is derived shown in the following proposition.

    Proposition 5. Let AMn(R) be nonnegative and irreducible and let the quantities M,N,bm,Ψ, mi(A),wi(A), i=1,,n satisfy the notations and assumptions of Theorem 4. Then ρ(A)=Ψ holds for some =1,,n if and only if w1(A)==wn(A)>0 or for some t=2,,, A satisfies the following conditions:

    (i) aii=M, for 1it1,

    (ii) aij=N and mj(A)mi(A)=bm, for 1in, 1jt1 and ji,

    (iii) wt(A)==wn(A).

    Proof. Suppose that A is nonnegative and irreducible with ρ(A)=Ψ for some 2n. Consider B=U1AU as constructed in the proof of Theorem 4, then B is also nonnegative and irreducible with ρ(B)=max1inri(B)=ρ(A)=Ψ. By Lemma 1, r1(B)==rn(B)=Ψ, and thus, for 1i1 cases (a) and (b) and for in cases (c) and (d) in Remark 1 hold. Let us firstly assume that w1(A)=w(A), then case (d) implies w1(A)==wn(A). On the other hand, if w1(A)>w(A), consider the smallest integer 2t such that wt(A)=w(A). Clearly, wi(A)>w(A) for integers 1it1, which implies that xi>1. Therefore, conditions (ⅰ) and (ⅱ) follow from the corresponding cases (a) and (b) of Remark 1 for 1i1 and the case (c) for in. Condition (ⅲ) is verified by case (d), since we have wt(A)==w(A)==wn(A).

    For the converse statement, if w1(A)==wn(A), then for 1n we obtain Ψ=w1(A) and by Proposition 3, ρ(A)=w1(A)=Ψ. If conditions (ⅰ)–(ⅲ) hold for some 2tn, then (a) and (b) hold for 1i1, and (c) and (d) hold for in, implying that ri(B)=Ψ for 1in, and thus by Lemma 1, ρ(A)=ρ(B)=Ψ.

    Next we present a sharp upper bound for the spectral radius of a nonnegative matrix proved by Duan, Zhou in [4] and Xing, Zhou in [11]. We outline the corresponding theorems, since their expressions will be compared to our results proved in Theorem 4.

    Theorem 6. ([4,Theorem 2.1]) Let AMn(R), A0 with row sums r1(A)r2(A)rn(A). Let M, N be the largest diagonal and non-diagonal elements of A, respectively. Suppose that N>0. For 1n, let

    ˆΦ=r(A)+MN2+(r(A)M+N2)2+N1i=1(ri(A)r(A)). (2.15)

    Then, ρ(A)ˆΦ for 1n.

    Theorem 7. ([11,Theorem 2.1]) Let AMn(R), A0 with row sums ri(A)>0, i=1,,n and average 2-row sums m1(A)m2(A)mn(A). Let M, N be the largest diagonal and non-diagonal elements of A, respectively. Suppose that N>0. Let b=max{rj(A)ri(A):1i,jn} and

    ˜Φ=m(A)+MNb2+(m(A)M+Nb2)2+Nb1i=1(mi(A)m(A)),

    for 1n. Then, ρ(A)˜Φ for 1n.

    At the following example we test the results of Theorem 4 in comparison to the ones of Theorems 6 and 7.

    Example 8. Consider the matrix

    A=(3110031122000001).

    The spectrum of A is σ(A)={0.8951,1,2.6027,4.2924}, which means that ρ(A)=4.2924. Clearly, its row sums are r1(A)=r2(A)=5,r3(A)=4,r4(A)=1, its average 2-row sums are m1(A)=4.8,m2(A)=4,m3(A)=5,m4(A)=1 and w1(A)=4.875,w2(A)=4.5,w3(A)=3.52,w4(A)=1, M=3, N=2, and bm=5. The assumptions of Theorem 4 hold and the quantities Ψ, 14, given by (2.1), are Ψ1=4.875,Ψ2=4.8173,Ψ3=5.4027, and Ψ4=7.7215. Hence, the inequality (2.3) yields ρ(A)4.8173.

    The assumptions of Theorem 6 are also satisfied, since r1(A)=r2(A)>r3(A)>r4(A) and the quantities ˆΦ, given by (2.15), are ˆΦ1=ˆΦ2=ˆΦ3=5, and ˆΦ4=5.6904. Then, Theorem 6 yields ρ(A)5.

    In order to apply Theorem 7, we use the permutation matrix

    P=(0010100001000001),

    such that A is similar to

    B=PAPT=(0220131010310001).

    Then the assumptions of Theorem 7 hold for the matrix B, since r1(B)=4, r2(B)=r3(B)=5,r4(B)=1, m1(B)=5,m2(B)=4.8, m3(B)=4,m4(B)=1 and M=3, N=2, b=5; the quantities ˜Φ are computed to be equal to ˜Φ1=5,˜Φ2=4.9671,˜Φ3=5.4462, and ˜Φ4=8.1355. Due to the similarity of matrices A,B and Theorem 7, ρ(A)=ρ(B)4.9671.

    Overall, Theorem 4 appears to be a refinement, since the upper bound of the spectral radius ρ(A) computed by the expressions of Theorem 4 is sharper than the ones computed by Theorems 6 and 7.

    At this point, let us consider some well-known nonnegative matrices representing connectivity properties of finite graphs. Let G be a simple undirected graph on vertex set V(G)={v1,,vn} and edge set E(G) with no isolated vertices. Recall that the adjacency matrix A(G)Mn(R) of G is a symmetric (0,1)-matrix, whose entries depend on whether the corresponding vertices are adjacent and the degree matrix D(G)=diag(d1,,dn)Mn(R) is the diagonal matrix with vertex degrees of G. Both of these matrices are linked to the signless Laplacian matrix of the graph defined as Q(G)=D(G)+A(G). In particular, the entries qij of Q(G) are given by

    qij={di,ifi=j,1,if ij and vivjE(G)0,otherwise.

    Apparently, Q(G) is a symmetric positive semi-definite matrix whose spectral radius ρ(Q(G)) is known as the signless Laplacian spectral radius of G. More information about the graph can be provided by the average 2-degree of vertex viV(G), mi=d1ij:vivjE(G)dj. The sequence {(di,mi)}ni=1 of pairs is called the sequence of degree pairs of G.

    The next example illustrates the formulas proved in Theorem 4 applied on the signless Laplacian matrix of an undirected graph as defined above.

    Example 9. Let an undirected graph G with signless Laplacian matrix

    Q(G)=(3111131111201102).

    The spectrum of Q(G) is σ(Q(G))={0.7639,2,5.2361}, which means that ρ(Q(G))=5.2361. Then r1(Q(G))=r2(Q(G))=6, r3(Q(G))=r4(Q(G))=4, m1(Q(G))=m2(Q(G))=5.3333, m3(Q(G))=m4(Q(G))=5 and w1(Q(G))=w2(Q(G))=5.8750, w3(Q(G))=w4(Q(G))=4.1333, M=3, N=1, and bm=1.0667. Apparently, the assumptions of Theorem 4 are satisfied and the quantities Ψ, 14, given by (2.1), are Ψ1=Ψ2=5.8750 and Ψ3=Ψ4=5.2527. Hence, the inequality (2.3) yields ρ(Q(G))5.2527.

    The assumptions of Theorem 6 are also satisfied, since r1(Q(G))=r2(Q(G))>r3(Q(G))=r4(Q(G)), and the quantities ˆΦ, given by (2.15), are ˆΦ1=ˆΦ2=6, ˆΦ3=ˆΦ4=5.2361. Thus, Theorem 6 yields ρ(A2)5.2361.

    The assumptions of Theorem 7 also hold, since m1(Q(G))=m2(Q(G))>m3(Q(G))=m4(Q(G)), and the quantities ˜Φ are equal to ˜Φ1=˜Φ2=5.3333 and ˜Φ3=˜Φ4=5.2656. Evidently, ρ(Q(G))5.2656.

    Concluding, we notice that the upper bound of ρ(Q(G)) calculated by Theorem 6 coincides with the exact value of ρ(Q(G)), which reveals that in this case Theorem 6 provides a sharper estimate than the corresponding ones computed by Theorems 4 and 7.

    In this section, we obtain a new result on the lower bound of the spectral radius of nonnegative matrices, and we compare these with the bounds investigated in [4,11].

    Theorem 10. Let AMn(R), A0 with row sums ri(A)>0 and average 2-row sums mi(A)>0 for all i=1,,n, and wi(A) as defined in (1.1) such that w1(A)w2(A)wn(A)>0. Let S,T be the smallest diagonal and off-diagonal elements of A, respectively. Denoting by cm=min{mj(A)mi(A):1i,jn,ij}, consider

    ψn=12(wn(A)+STcm+Δn), (3.1)

    where

    Δn=(wn(A)S+Tcm)2+4Tcmn1j=1(wj(A)wn(A)). (3.2)

    Then

    ρ(A)ψn. (3.3)

    Proof. To simplify the exposition of our calculations, we let mi=mi(A), and wi=wi(A) for 1in.

    If T=0, equality (3.1) degenerates to ψn=wn due to the fact that wnannS and Theorem 10 is apparent from Proposition 3. Consequently, in what follows we assume T>0.

    Let U=diag(m1x1,m2x2,,mn1xn1,mn) be an n×n diagonal matrix, where xj1 for 1jn1 is a variable to be determined later and let B=U1AU. Due to similarity, A and B have the same eigenvalues, hence, ρ(A)=ρ(B).

    For 1in1 we refer to the equality (2.4) with =n, that is,

    ri(B)=ri(U1AU)=1xi(nj=1aijmjmi+n1j=1jiaijmjmi(xj1)+aii(xi1)).

    Moreover, aiiS, aijT, and mjmicm, for 1jn1 and ji. Combining these inequalities with the definition of wi in (1.1), the latter equality is formulated as

    ri(B)1xi(wi+Tcmn1j=1ji(xj1)+S(xi1)). (3.4)

    Furthermore, for the special case i=n the equality (2.6) can be written as

    rn(B)=rn(U1AU)=n1j=1anjmjmn(xj1)+nj=1anjmjmn=n1j=1anjmjmn(xj1)+wnTcmn1j=1(xj1)+wn. (3.5)

    Now, the variable xj will be formed by the roots of the equation

    ψ2n(wn+STcm)ψn+wn(STcm)Tcmn1j=1(wjwn)=0. (3.6)

    The quadratic equation in (3.6) has real roots, since its discriminant

    Δn(wn+STcm)24(wn(STcm)Tcmn1j=1(wjwn))=(wnS+Tcm)2+4Tcmn1j=1(wjwn)

    is a positive number, due to T>0, cm>0 and the monotonicity of the sequence {wi}ni=1 of positive numbers. Hence, (3.6) has a positive real root

    ψn=12(wn+STcm+Δn), (3.7)

    which is used in the construction of

    xj=1+wjwnψnS+Tcmxj1=wjwnψnS+Tcm, (3.8)

    for 1jn1. If n1j=1(wjwn)>0, it is clear by relation (3.7) that

    ψn>12(wn+STcm+|wnS+Tcm|)12(wn+STcm(wnS+Tcm))=STcm,

    otherwise, w1==wnSwnS+Tcm>0 and (3.7) yields

    ψn=12(wn+STcm+|wnS+Tcm|)>12(wn+STcm(wnS+Tcm))=STcm.

    Both cases ensure xj10 in (3.8).

    Moreover, from (3.6) we derive

    Tcmn1j=1(wjwn)=ψ2n(wn+STcm)ψn+wn(STcm)=ψn(ψnwnS+Tcm)+wn(STcm)=ψn(ψnS+Tcm)ψnwn+wn(STcm)=ψn(ψnS+Tcm)wn(ψnS+Tcm)=(ψnS+Tcm)(ψnwn). (3.9)

    Overall, we substitute xj10, j=1,,n1, from (3.8) and Tcmn1j=1(wjwn) from (3.9) into the inequality (3.4). Hence, for 1in1, we have

    ri(B)1xi(wi+Tcmn1j=1ji(xj1)+S(xi1))=1xi(wi+Tcmn1j=1(xj1)+S(xi1)Tcm(xi1))=1xi(wi+Tcmn1j=1(xj1)+(STcm)(xi1))=1xi(wi+Tcmn1j=1wjwnψnS+Tcm+(STcm)wiwnψnS+Tcm)=1xi(wi+(ψnS+Tcm)(ψnwn)ψnS+Tcm+(STcm)(wiwn)ψnS+Tcm)=wi(ψnS+Tcm)+(ψnS+Tcm)(ψnwn)+(STcm)(wiwn)xi(ψnS+Tcm)=(ψnS+Tcm)(ψn+wiwn)+(STcm)(wiwn)xi(ψnS+Tcm)=ψn(ψnS+Tcm)+ψn(wiwn)xi(ψnS+Tcm)=ψn(ψnS+Tcm+wiwn)xi(ψnS+Tcm)=ψn(ψnS+Tcm+wiwn)ψnS+Tcm+wiwnψnS+Tcm(ψnS+Tcm)=ψn. (3.10)

    Also, by substituting xj10 and Tcmn1j=1(wjwn) from (3.8) and (3.9), respectively, into the inequality (3.5), we may write

    rn(B)wn+Tcmn1j=1wjwnψnS+Tcm=wn+TcmψnS+Tcmn1j=1(wjwn)=wn+(ψnS+Tcm)(ψnwn)ψnS+Tcm=ψn. (3.11)

    Hence, for all 1in the inequalities (3.10) and (3.11) confirm

    ri(B)ψn.

    By Lemma 1,

    ρ(A)=ρ(B)min1in{ri(B)}ψn,

    verifying the validity of (3.3).

    Analogously to equality cases stated in Remark 1 for the upper bounds of the spectral radius of a nonnegative matrix, we may have corresponding equality cases for the lower bounds as stated in the following arguments.

    Remark 2. Inequalities (3.4) and (3.5) are reduced to equalities for the special values of S,T,cm,xi, as the following cases indicate.

    (ⅰ) For 1in1, (3.4) is given by equality if and only if (a) and (b) hold:

    (a) xi=1 or aii=S, when xi>1,

    (b) xj=1 or aij=T and mjmi=cm, when xj>1 with 1jn1 and ji.

    (ⅱ) Inequality (3.5) is given by equality if and only if (c) holds:

    (c) xj=1 or anj=T and mjmn=cm, when xj>1 with 1jn1.

    Using Remark 2 and Proposition 3 for a nonnegative and irreducible matrix, a special formulation of the spectral radius of the matrix is derived shown in the following proposition.

    Proposition 11. Let AMn(R) be nonnegative and irreducible and let the quantities S,T,cm,ψn, mi(A),wi(A), i=1,,n satisfy the notations and assumptions of Theorem 10. Then ρ(A)=ψn holds if and only if w1(A)==wn(A)>0 or T>0 and for some t=2,,n, A satisfies the following conditions:

    (i) aii=S, for 1it1,

    (ii) aij=T and mj(A)mi(A)=cm, for 1in, 1jt1 and ji,

    (iii) wt(A)==wn(A),

    Proof. Suppose that A is nonnegative and irreducible with ρ(A)=ψn. Consider B=U1AU as constructed in the proof of Theorem 10, then B is also nonnegative and irreducible with ρ(B)=max1inri(B)=ρ(A)=Ψ. By Lemma 1, r1(B)==rn(B)=ψn, and thus, for 1in1 cases (a) and (b) and for i=n case (c) in Remark 2 hold. If w1(A)>wn(A), consider the smallest integer 2tn such that wt(A)=wn(A). Clearly, wi(A)>wn(A) for 1it1, which implies that xi>1. Therefore, conditions (ⅰ)-(ⅲ) of Theorem 10 follow from cases (a) and (b) for 1in1 and (c) of Remark 2.

    Conversely, if w1(A)==wn(A), then by Proposition 3, ρ(A)=wn(A)=ψn. If conditions (ⅰ)–(ⅲ) hold for some 2tn, then (a) and (b) for 1in1, and (c) hold, implying that ri(B)=ψn for 1in, and thus by Lemma 1, ρ(A)=ρ(B)=ψn.

    The following statements concern lower bounds for the spectral radius of nonnegative matrices proved by Duan, Zhou [4] and Xing, Zhou [11]. They are presented here for reasons of comparison.

    Theorem 12. ([4,Theorem 2.2]) Let AMn(R), A0 with row sums r1(A)r2(A)rn(A). Let S, T be the smallest diagonal and non-diagonal elements of A, respectively. Let

    ˆϕn=rn(A)+ST2+(rn(A)S+T2)2+Tn1i=1(ri(A)rn(A)). (3.12)

    Then, ρ(A)ˆϕn.

    Theorem 13. ([11,Theorem 2.3]) Let AMn(R), A0 with all row sums positive and average 2-row sums such that m1(A)m2(A)mn(A). Let S, T be the smallest diagonal and non-diagonal elements of A, respectively. Denote by c=min{rj(A)ri(A):1i,jn} and by

    ˜ϕn=mn(A)+STc2+(mn(A)S+Tc2)2+Tcn1i=1(mi(A)mn(A)). (3.13)

    Then, ρ(A)˜ϕn.

    We provide the following illustrative example, comparing the results among Theorems 10, 12 and 13.

    Example 14. I. Consider the matrix

    A=(6222022122022220).

    The spectrum of A is σ(A)={2,0.7321,2.7321,8}, which means that ρ(A)=8. With respect to the notations of Theorem 10, r1(A)=12,r2(A)=5,r3(A)=r4(A)=6, m1(A)=8.8333,m2(A)=5.6,m3(A)=m4(A)=7.6667 and w1(A)=10.7396,w2(A)=6.1071,w3(A)=w4(A)=5.7652, and cm=0.6340. The equality S=T=0 in (3.1) and (3.2) yields ρ(A)w4(A)=5.7652.

    Using the permutation matrix P=(1000000100100100), we obtain the similar matrix

    B=PAPT=(6222202222020122),

    for which r1(B)=12,r2(B)=r3(B)=6,r4(B)=5, and m1(B)=8.8333,m2(B)=m3(B)=7.6667,m4(B)=5.6, and S=T=0.

    Clearly, the assumptions of Theorem 12 are satisfied and by (3.12), ˆϕ4=r4(B)=5. Then, the similarity of A,B implies ρ(A)=ρ(B)5.

    The assumptions of Theorem 13 are also satisfied and by (3.13), ˜ϕ4=m4(B)=5.6. Then ρ(A)=ρ(B)5.6.

    Overall, the lower bound of Theorem 10 appears to be sharper than the ones evaluated by Theorems 12 and 13.

    II. Let again the signless Laplacian matrix Q(G) be presented at Example 9 with ρ(Q(G))=5.2361. Recall the orderings

    r1(Q(G))=r2(Q(G))>r3(Q(G))=r4(Q(G))m1(Q(G))=m2(Q(G))>m3(Q(G))=m4(Q(G))w1(Q(G))=w2(Q(G)))>w3(Q(G))=w4(Q(G)),

    each of them satisfying the assumption of Theorems 12, 13 and 10, respectively. Morover, S=2, T=0, and cm=0.9375. Then, Theorem 10 provides the estimate ρ(Q(G))ψ4=4.1333, Theorem 12 yields ρ(Q(G))ˆϕ4=4 and lastly, Theorem 13 implies ρ(Q(G))˜ϕ4=5. In this case, the lower bound of Theorem 13 appears to be sharper than the ones from Theorems 10 and 12.

    The authors declare no conflict of interest.



    [1] M. Adam, Aik. Aretaki, Sharp bounds for eigenvalues of the generalized k, m-step Fibonacci matrices, Proceedings of the 3rd International Conference on Numerical Analysis and Scientific Computation with Applications (NASCA18), Kalamata, Greece, (2018). Available from: http://nasca18.math.uoa.gr/participants-nbsp.html.
    [2] A. Brauer, I.C. Gentry, Bounds for the greatest characteristic root of an irreducible nonnegative matrix, Linear Algebra its Appl., 8 (1974), 105-107. doi: 10.1016/0024-3795(74)90048-2
    [3] F. Duan, K. Zhang, An algorithm of diagonal transformation for Perron root of nonnegative irreducible matrices, Appl. Math. Comput., 175 (2006), 762-772.
    [4] X. Duan, B. Zhou, Sharp bounds on the spectral radius of a nonnegative matrix, Linear Algebra its Appl., 439 (2013), 2961-2970. doi: 10.1016/j.laa.2013.08.026
    [5] G. Frobenius, Über Matrizen aus nicht negativen Elementen, Sitzungsber, Kön. Preuss. Akad. Wiss. Berlin, (1912), 465-477.
    [6] J. He, Y.M. Liu, J.K. Tian, et al., Some new sharp bounds for the spectral radius of a nonnegative matrix and its application, J. Inequalities Appl., 260 (2017), 1-6.
    [7] W. Hong, L. You, Further results on the spectral radius of matrices and graphs, Appl. Math. Comput., 239 (2014), 326-332.
    [8] R.A. Horn, C.R. Johnson, Matrix Analysis, Cambridge University Press, Cambridge, second edition, 2013.
    [9] H. Lin, B. Zhou, On sharp bounds for spectral radius of nonnegative matrices, Linear Multilinear Algebra, 65 (2017), 1554-1565. doi: 10.1080/03081087.2016.1246514
    [10] A. Melman, Upper and lower bounds for the Perron root of a nonnegative matrix, Linear Multilinear Algebra, 61 (2013), 171-181. doi: 10.1080/03081087.2012.667096
    [11] R. Xing, B. Zhou, Sharp bounds for the spectral radius of nonnegative matrices, Linear Algebra its Appl., 449 (2014), 194-209. doi: 10.1016/j.laa.2014.02.031
    [12] C. Wen, T. Z. Huang, A modified algorithm for the Perron root of a nonnegative matrix, Appl. Math. Comput., 217 (2011), 4453-4458.
    [13] P. Walters, An introduction to ergodic theory, New York (NY), Springer-Verlag, 1982.
  • This article has been cited by:

    1. Fotis Babouklis, Maria Adam, Nicholas Assimakis, 2020, Bounds on the Spectral Radius of Nonnegative Matrices, 978-1-7281-6695-7, 51, 10.1109/MACISE49704.2020.00016
    2. Maria Adam, Aikaterini Aretaki, Bounds for the spectral radius of nonnegative matrices and generalized Fibonacci matrices, 2022, 10, 2300-7451, 308, 10.1515/spma-2022-0165
    3. Gang Wang, Jinfa Liu, An Improved Diagonal Transformation Algorithm for the Maximum Eigenvalue of Zero Symmetric Nonnegative Matrices, 2022, 14, 2073-8994, 1707, 10.3390/sym14081707
    4. Maria Adam, Iro Oikonomou, Aikaterini Aretaki, Sequences of lower and upper bounds for the spectral radius of a nonnegative matrix, 2023, 00243795, 10.1016/j.laa.2023.03.006
    5. Zulfiya R. Gabidullina, 2024, Chapter 3, 978-3-031-62791-0, 39, 10.1007/978-3-031-62792-7_3
    6. Chuang Lü, Lihua You, Yufei Huang, Sharp Bounds for the Spectral Radii of Nonnegative Tensors, 2023, 18, 2731-8648, 883, 10.1007/s11464-020-0006-2
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4614) PDF downloads(539) Cited by(6)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog