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
[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 |
Let Mn(R) be the algebra of n×n real matrices. We refer to A=(aij)ni,j=1∈Mn(R) as a nonnegative or a positive matrix, when each aij≥0 or aij>0, respectively, denoted by writing A≥0 or A>0. The matrix A∈Mn(R) is called irreducible if and only if (I+A)n−1>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 1≤i≤n, the quantities
ri(A)=n∑j=1aijandMi(A)=n∑j=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=max1≤i≤n{aii}andN=max1≤i,j≤ni≠j{aij}, |
respectively, as well as, its smallest diagonal and off-diagonal elements,
S=min1≤i≤n{aii}andT=min1≤i,j≤ni≠j{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 A∈Mn(R), A≥0 with i-th row sum ri(A), i=1,…,n and largest diagonal element M. Then
ρ(A)≥M, |
and
min1≤i≤n{ri(A)}≤ρ(A)≤max1≤i≤n{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 A∈Mn(R), A≥0. Then for any positive vector x∈Rn we have
min1≤i≤n{1xin∑j=1aijxj}≤ρ(A)≤max1≤i≤n{1xin∑j=1aijxj}. |
If A is also irreducible, then
ρ(A)=maxx>0min1≤i≤n{1xin∑j=1aijxj}=minx>0max1≤i≤n{1xin∑j=1aijxj}. |
Combining the definition of wi(A) in (1.1) with Lemma 2 by taking a vector x with components mi(A)>0, 1≤i≤n, 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 A∈Mn(R), A≥0. Then
min1≤i≤n{wi(A)}≤ρ(A)≤max1≤i≤n{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 A∈Mn(R), A≥0 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):1≤i,j≤n,i≠j}, consider
Ψℓ=12(wℓ(A)+M−Nbm+√Δℓ),ℓ=1,…,n, | (2.1) |
where
Δℓ=(wℓ(A)−M+Nbm)2+4Nbmℓ−1∑j=1(wj(A)−wℓ(A)). | (2.2) |
Then
ρ(A)≤min{Ψℓ:1≤ℓ≤n}. | (2.3) |
Proof. To simplify the exposition of our calculations, we let mi=mi(A), and wi=wi(A) for 1≤i≤n.
Consider ℓ=1. Due to Lemma 1, Proposition 3 and our assumptions,
M≤ρ(A)≤max1≤i≤n{wi}≡w1, |
which yield that
w1−M≥0⇒w1−M+Nbm≥Nbm>0. |
Then the quantities in (2.1), (2.2) give
Ψ1=12(w1+M−Nbm+√(w1−M+Nbm)2)=12(w1+M−Nbm+w1−M+Nbm)=w1. |
Then ρ(A)≤Ψ1.
Consider 2≤ℓ≤n. Let U=diag(m1x1,…,mℓ−1xℓ−1,mℓ,…,mn) be an n×n diagonal matrix, where xj≥1 is a variable to be determined later for 1≤j≤ℓ−1 and let B=U−1AU. Due to similarity, A and B have the same eigenvalues, hence ρ(A)=ρ(B).
For 1≤i≤ℓ−1, we derive
ri(B)=ri(U−1AU)=1xi(ℓ−1∑j=1aijmjmixj+n∑j=ℓaijmjmi)=1xi(ℓ−1∑j=1j≠iaijmjmixj+aiixi+n∑j=1aijmjmi−ℓ−1∑j=1aijmjmi)=1xi(n∑j=1aijmjmi+ℓ−1∑j=1j≠iaijmjmixj−ℓ−1∑j=1j≠iaijmjmi+aiixi−aii)=1xi(n∑j=1aijmjmi+ℓ−1∑j=1j≠iaijmjmi(xj−1)+aii(xi−1)). | (2.4) |
Obviously, aii≤M, aij≤N, and mjmi≤bm for 1≤j≤ℓ−1 and j≠i. Combining these inequalities with the definition of wi in (1.1), the equality (2.4) is formulated as
ri(B)≤1xi(wi+Nbmℓ−1∑j=1j≠i(xj−1)+M(xi−1)),i=1,…,ℓ−1. | (2.5) |
For ℓ≤i≤n, we derive
ri(B)=ri(U−1AU)=ℓ−1∑j=1aijmjmixj+n∑j=ℓaijmjmi=ℓ−1∑j=1aijmjmixj+n∑j=1aijmjmi−ℓ−1∑j=1aijmjmi=n∑j=1aijmjmi+ℓ−1∑j=1aijmjmi(xj−1) | (2.6) |
≤wℓ+Nbmℓ−1∑j=1(xj−1). | (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ℓ+M−Nbm)Ψℓ+wℓ(M−Nbm)−Nbmℓ−1∑j=1(wj−wℓ)=0. | (2.8) |
In particular, the trinomials in (2.8) have discriminant
Δℓ≡(wℓ+M−Nbm)2−4(wℓ(M−Nbm)−Nbmℓ−1∑j=1(wj−wℓ))=(wℓ−M+Nbm)2+4Nbmℓ−1∑j=1(wj−wℓ). |
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ℓ+M−Nbm+√Δℓ),ℓ=2,…,n. | (2.9) |
Now, for 1≤j≤ℓ−1, we consider
xj=1+wj−wℓΨℓ−M+Nbm⇔xj−1=wj−wℓΨℓ−M+Nbm, | (2.10) |
where Ψℓ is given by (2.9). If ∑ℓ−1j=1(wj−wℓ)>0, it is clear by relation (2.9) that
Ψℓ>12(wℓ+M−Nbm+|wℓ−M+Nbm|)≥12(wℓ+M−Nbm−(wℓ−M+Nbm))=M−Nbm, |
otherwise, w1=⋯=wℓ≥M⇒wℓ−M+Nbm>0 and (2.9) yields
Ψℓ=12(wℓ+M−Nbm+|wℓ−M+Nbm|)>12(wℓ+M−Nbm−(wℓ−M+Nbm))=M−Nbm. |
Both cases ensure xj−1≥0 in (2.10).
Additionally, by the expression (2.8), we may write
Nbmℓ−1∑j=1(wj−wℓ)=Ψ2ℓ−(wℓ+M−Nbm)Ψℓ+wℓ(M−Nbm)=Ψℓ(Ψℓ−wℓ−M+Nbm)+wℓ(M−Nbm)=Ψℓ(Ψℓ−M+Nbm)−Ψℓwℓ+wℓ(M−Nbm)=Ψℓ(Ψℓ−M+Nbm)−wℓ(Ψℓ−M+Nbm)=(Ψℓ−M+Nbm)(Ψℓ−wℓ). | (2.11) |
Overall, we take xj−1≥0, j=1,…,ℓ−1 from (2.10) and Nbm∑ℓ−1j=1(wj−wℓ) from (2.11) and substitute them into the inequality (2.5). Hence, for 1≤i≤ℓ−1, we obtain
ri(B)≤1xi(wi+Nbmℓ−1∑j=1j≠i(xj−1)+M(xi−1))=1xi(wi+Nbmℓ−1∑j=1(xj−1)+M(xi−1)−Nbm(xi−1))=1xi(wi+Nbmℓ−1∑j=1(xj−1)+(M−Nbm)(xi−1))=1xi(wi+Nbmℓ−1∑j=1wj−wℓΨℓ−M+Nbm+(M−Nbm)wi−wℓΨℓ−M+Nbm)=1xi(wi+(Ψℓ−M+Nbm)(Ψℓ−wℓ)Ψℓ−M+Nbm+(M−Nbm)wi−wℓΨℓ−M+Nbm)=wi(Ψℓ−M+Nbm)+(Ψℓ−M+Nbm)(Ψℓ−wℓ)+(M−Nbm)(wi−wℓ)xi(Ψℓ−M+Nbm)=(Ψℓ−M+Nbm)Ψℓ+(wi−wℓ)(Ψℓ−M+Nbm+M−Nbm)xi(Ψℓ−M+Nbm)=(Ψℓ−M+Nbm)Ψℓ+(wi−wℓ)ΨℓΨℓ−M+Nbm+wi−wℓΨℓ−M+Nbm(Ψℓ−M+Nbm)=(Ψℓ−M+Nbm+wi−wℓ)ΨℓΨℓ−M+Nbm+wi−wℓ=Ψℓ, | (2.12) |
where ℓ=2,…,n.
Similarly, for ℓ≤i≤n and 1≤j≤ℓ−1 we substitute the relations (2.10) and (2.11) into the inequality (2.7), which can be written as
ri(B)≤wi+Nbmℓ−1∑j=1(xj−1)=wi+NbmΨℓ−M+Nbmℓ−1∑j=1(wj−wℓ)=wi+(Ψℓ−M+Nbm)(Ψℓ−wℓ)Ψℓ−M+Nbm=wi+Ψℓ−wℓ≤wℓ+Ψℓ−wℓ=Ψℓ. | (2.13) |
Thus, for all 2≤ℓ≤n and 1≤i≤n the inequalities (2.12) and (2.13) verify
ri(B)≤Ψℓ, |
and by Lemma 1,
ρ(A)=ρ(B)≤max1≤i≤n{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 1≤i≤ℓ−1, (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 1≤j≤ℓ−1 and j≠i.
(ⅱ) For ℓ≤i≤n, (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 1≤j≤ℓ−1,
(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 A∈Mn(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 1≤i≤t−1,
(ii) aij=N and mj(A)mi(A)=bm, for 1≤i≤n, 1≤j≤t−1 and j≠i,
(iii) wt(A)=⋯=wn(A).
Proof. Suppose that A is nonnegative and irreducible with ρ(A)=Ψℓ for some 2≤ℓ≤n. Consider B=U−1AU as constructed in the proof of Theorem 4, then B is also nonnegative and irreducible with ρ(B)=max1≤i≤nri(B)=ρ(A)=Ψℓ. By Lemma 1, r1(B)=⋯=rn(B)=Ψℓ, and thus, for 1≤i≤ℓ−1 cases (a) and (b) and for ℓ≤i≤n 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 2≤t≤ℓ such that wt(A)=wℓ(A). Clearly, wi(A)>wℓ(A) for integers 1≤i≤t−1, which implies that xi>1. Therefore, conditions (ⅰ) and (ⅱ) follow from the corresponding cases (a) and (b) of Remark 1 for 1≤i≤ℓ−1 and the case (c) for ℓ≤i≤n. 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 1≤ℓ≤n we obtain Ψℓ=w1(A) and by Proposition 3, ρ(A)=w1(A)=Ψℓ. If conditions (ⅰ)–(ⅲ) hold for some 2≤t≤ℓ≤n, then (a) and (b) hold for 1≤i≤ℓ−1, and (c) and (d) hold for ℓ≤i≤n, implying that ri(B)=Ψℓ for 1≤i≤n, 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 A∈Mn(R), A≥0 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 1≤ℓ≤n, let
ˆΦℓ=rℓ(A)+M−N2+√(rℓ(A)−M+N2)2+Nℓ−1∑i=1(ri(A)−rℓ(A)). | (2.15) |
Then, ρ(A)≤ˆΦℓ for 1≤ℓ≤n.
Theorem 7. ([11,Theorem 2.1]) Let A∈Mn(R), A≥0 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):1≤i,j≤n} and
˜Φℓ=mℓ(A)+M−Nb2+√(mℓ(A)−M+Nb2)2+Nbℓ−1∑i=1(mi(A)−mℓ(A)), |
for 1≤ℓ≤n. Then, ρ(A)≤˜Φℓ for 1≤ℓ≤n.
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 Ψℓ, 1≤ℓ≤4, 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 i≠j and vivj∈E(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 vi∈V(G), mi=d−1i∑j:vivj∈E(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 Ψℓ, 1≤ℓ≤4, 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 A∈Mn(R), A≥0 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):1≤i,j≤n,i≠j}, consider
ψn=12(wn(A)+S−Tcm+√Δn), | (3.1) |
where
Δn=(wn(A)−S+Tcm)2+4Tcmn−1∑j=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 1≤i≤n.
If T=0, equality (3.1) degenerates to ψn=wn due to the fact that wn≥ann≥S and Theorem 10 is apparent from Proposition 3. Consequently, in what follows we assume T>0.
Let U=diag(m1x1,m2x2,…,mn−1xn−1,mn) be an n×n diagonal matrix, where xj≥1 for 1≤j≤n−1 is a variable to be determined later and let B=U−1AU. Due to similarity, A and B have the same eigenvalues, hence, ρ(A)=ρ(B).
For 1≤i≤n−1 we refer to the equality (2.4) with ℓ=n, that is,
ri(B)=ri(U−1AU)=1xi(n∑j=1aijmjmi+n−1∑j=1j≠iaijmjmi(xj−1)+aii(xi−1)). |
Moreover, aii≥S, aij≥T, and mjmi≥cm, for 1≤j≤n−1 and j≠i. Combining these inequalities with the definition of wi in (1.1), the latter equality is formulated as
ri(B)≥1xi(wi+Tcmn−1∑j=1j≠i(xj−1)+S(xi−1)). | (3.4) |
Furthermore, for the special case i=n the equality (2.6) can be written as
rn(B)=rn(U−1AU)=n−1∑j=1anjmjmn(xj−1)+n∑j=1anjmjmn=n−1∑j=1anjmjmn(xj−1)+wn≥Tcmn−1∑j=1(xj−1)+wn. | (3.5) |
Now, the variable xj will be formed by the roots of the equation
ψ2n−(wn+S−Tcm)ψn+wn(S−Tcm)−Tcmn−1∑j=1(wj−wn)=0. | (3.6) |
The quadratic equation in (3.6) has real roots, since its discriminant
Δn≡(wn+S−Tcm)2−4(wn(S−Tcm)−Tcmn−1∑j=1(wj−wn))=(wn−S+Tcm)2+4Tcmn−1∑j=1(wj−wn) |
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+S−Tcm+√Δn), | (3.7) |
which is used in the construction of
xj=1+wj−wnψn−S+Tcm⇔xj−1=wj−wnψn−S+Tcm, | (3.8) |
for 1≤j≤n−1. If ∑n−1j=1(wj−wn)>0, it is clear by relation (3.7) that
ψn>12(wn+S−Tcm+|wn−S+Tcm|)≥12(wn+S−Tcm−(wn−S+Tcm))=S−Tcm, |
otherwise, w1=⋯=wn≥S⇒wn−S+Tcm>0 and (3.7) yields
ψn=12(wn+S−Tcm+|wn−S+Tcm|)>12(wn+S−Tcm−(wn−S+Tcm))=S−Tcm. |
Both cases ensure xj−1≥0 in (3.8).
Moreover, from (3.6) we derive
Tcmn−1∑j=1(wj−wn)=ψ2n−(wn+S−Tcm)ψn+wn(S−Tcm)=ψn(ψn−wn−S+Tcm)+wn(S−Tcm)=ψn(ψn−S+Tcm)−ψnwn+wn(S−Tcm)=ψn(ψn−S+Tcm)−wn(ψn−S+Tcm)=(ψn−S+Tcm)(ψn−wn). | (3.9) |
Overall, we substitute xj−1≥0, j=1,…,n−1, from (3.8) and Tcm∑n−1j=1(wj−wn) from (3.9) into the inequality (3.4). Hence, for 1≤i≤n−1, we have
ri(B)≥1xi(wi+Tcmn−1∑j=1j≠i(xj−1)+S(xi−1))=1xi(wi+Tcmn−1∑j=1(xj−1)+S(xi−1)−Tcm(xi−1))=1xi(wi+Tcmn−1∑j=1(xj−1)+(S−Tcm)(xi−1))=1xi(wi+Tcmn−1∑j=1wj−wnψn−S+Tcm+(S−Tcm)wi−wnψn−S+Tcm)=1xi(wi+(ψn−S+Tcm)(ψn−wn)ψn−S+Tcm+(S−Tcm)(wi−wn)ψn−S+Tcm)=wi(ψn−S+Tcm)+(ψn−S+Tcm)(ψn−wn)+(S−Tcm)(wi−wn)xi(ψn−S+Tcm)=(ψn−S+Tcm)(ψn+wi−wn)+(S−Tcm)(wi−wn)xi(ψn−S+Tcm)=ψn(ψn−S+Tcm)+ψn(wi−wn)xi(ψn−S+Tcm)=ψn(ψn−S+Tcm+wi−wn)xi(ψn−S+Tcm)=ψn(ψn−S+Tcm+wi−wn)ψn−S+Tcm+wi−wnψn−S+Tcm(ψn−S+Tcm)=ψn. | (3.10) |
Also, by substituting xj−1≥0 and Tcm∑n−1j=1(wj−wn) from (3.8) and (3.9), respectively, into the inequality (3.5), we may write
rn(B)≥wn+Tcmn−1∑j=1wj−wnψn−S+Tcm=wn+Tcmψn−S+Tcmn−1∑j=1(wj−wn)=wn+(ψn−S+Tcm)(ψn−wn)ψn−S+Tcm=ψn. | (3.11) |
Hence, for all 1≤i≤n the inequalities (3.10) and (3.11) confirm
ri(B)≥ψn. |
By Lemma 1,
ρ(A)=ρ(B)≥min1≤i≤n{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 1≤i≤n−1, (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 1≤j≤n−1 and j≠i.
(ⅱ) 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 1≤j≤n−1.
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 A∈Mn(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 1≤i≤t−1,
(ii) aij=T and mj(A)mi(A)=cm, for 1≤i≤n, 1≤j≤t−1 and j≠i,
(iii) wt(A)=⋯=wn(A),
Proof. Suppose that A is nonnegative and irreducible with ρ(A)=ψn. Consider B=U−1AU as constructed in the proof of Theorem 10, then B is also nonnegative and irreducible with ρ(B)=max1≤i≤nri(B)=ρ(A)=Ψℓ. By Lemma 1, r1(B)=⋯=rn(B)=ψn, and thus, for 1≤i≤n−1 cases (a) and (b) and for i=n case (c) in Remark 2 hold. If w1(A)>wn(A), consider the smallest integer 2≤t≤n such that wt(A)=wn(A). Clearly, wi(A)>wn(A) for 1≤i≤t−1, which implies that xi>1. Therefore, conditions (ⅰ)-(ⅲ) of Theorem 10 follow from cases (a) and (b) for 1≤i≤n−1 and (c) of Remark 2.
Conversely, if w1(A)=⋯=wn(A), then by Proposition 3, ρ(A)=wn(A)=ψn. If conditions (ⅰ)–(ⅲ) hold for some 2≤t≤n, then (a) and (b) for 1≤i≤n−1, and (c) hold, implying that ri(B)=ψn for 1≤i≤n, 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 A∈Mn(R), A≥0 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)+S−T2+√(rn(A)−S+T2)2+Tn−1∑i=1(ri(A)−rn(A)). | (3.12) |
Then, ρ(A)≥ˆϕn.
Theorem 13. ([11,Theorem 2.3]) Let A∈Mn(R), A≥0 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):1≤i,j≤n} and by
˜ϕn=mn(A)+S−Tc2+√(mn(A)−S+Tc2)2+Tcn−1∑i=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. |
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 |