724
Views
0
CrossRef citations to date
0
Altmetric
Research Articles

Two inverse eigenvalue problems for matrices whose graphs are trees

, &
Pages 120-124 | Received 01 Jun 2023, Accepted 02 Jul 2023, Published online: 19 Jul 2023

Abstract

Inverse eigenvalue problem refers to the problem of reconstructing a matrix of a desired structure from a prescribed eigendata. In this paper, we discuss two additive inverse eigenvalue problems for matrices whose graph is tree. In order to analyze the problems, the vertices of the given tree are labeled in a suitable way so that concrete recurrence relations between the characteristic polynomials of the leading principal submatrices of the matrix can be obtained in a simple form. The method of obtaining the entries of the required matrix is constructive and provides an algorithm for computing the same. We provide some numerical examples to illustrate the results. The computations are done using SCILAB by feeding the eigendata and the adjacency pattern of the tree as inputs.

AMS CLASSIFICATION:

1 Introduction

Inverse eigenvalue problem (IEP) is the problem of reconstructing matrices of pre-assigned structure using the prescribed eigendata. It appears in various applications of mathematics and engineering [Citation3, Citation11, Citation12, Citation21]. One of the ways of describing the structure of the matrix is to represent it by its corresponding graph, where the adjacencies of vertices of the graph correspond to the nonzero off-diagonal entries of the matrix. Inverse problems of constructing matrices described by graphs have been studied by several authors. The off-diagonal entries of such matrices are governed by the adjacency pattern of the vertices in the graph. Graphs are helpful in representing various physical systems in control theory, vibration of string, geophysics, circuit theory, graph theory etc. as can be seen in the papers [Citation3–5] and the eigenvalues of the matrices associated with the graphs govern the behavior of the concerned systems. A detailed characterization of inverse eigenvalue problems has been made by [Citation4]. Inverse eigenvalue problems evolved from the discretization of inverse Strum-Liouville problems and this lead to the study of IEPs of Jacobi matrices [Citation7, Citation8], Pseudo-Jacobi matrices [Citation22], Periodic Jacobi matrices [Citation23]. IEPs corresponding to different trees have been studied extensively in recent times by several authors. IEPs of constructing matrices using the maximal and minimal eigenvalues of each leading principal submatrices have been solved for bordered diagonal matrices [Citation13, Citation16], doubly arrow matrices [Citation15]. The problem has been solved for various matrices whose graphs are special trees [Citation1, Citation2, Citation6, Citation17, Citation19]. A general solution of the extremal IEP for the case of distinct eigenvalues has been given by [Citation20]. [Citation14] studied the reconstruction of a special type of Leslie and doubly Leslie matrices, companion and doubly companion matrices, using maximal eigenvalue of each of the leading principal submatrix and its corresponding eigenvector.

In this paper, we have discussed two additive inverse eigenvalue problems for matrices whose graphs are trees. For a given tree T and a given diagonal matrix D, the first problem is to construct a hollow matrix A corresponding to T, when the largest eigenvalues of each of the leading principal submatrices of the matrix A + D are given. The second problem is to construct a matrix A corresponding to T, when the largest eigenvalues of each of the leading principal submatrices of the matrix A + D are given.

In Section 2, the basic definitions and primary results necessary for solving the above inverse eigenvalue problems have been briefly discussed. Section 4 defines the problems to be discussed in this paper. To solve the problems, the vertices of the trees have been labeled in a suitable manner. Solution of the Problems 1, 2 have been discussed in Sections 3.1 and 3.2, respectively. In Section 4, numerical solutions of the problems are obtained in SCILAB by feeding the eigendata, entries of the diagonal matrices and the adjacency pattern as inputs.

2 Preliminaries

Throughout this paper, the term graph will mean simple undirected graph without any multiple edges and self loops. A graph G consists of a finite set V of elements called vertex set and a set E of unordered pairs of distinct vertices called edge set. The order of the graph G is the number of elements in the vertex set V. Vertices v1 and v2 are adjacent or neighbors, if e={v1,v2}E. A walk in G, joining vertices u and v, is a sequence γ of vertices u=v0,v1,,vk1,vk=v such that {vi,vi+1} is an edge for each i=0,1,,k1. If the vertices u=v0,v1,,vk1,vk=v are distinct, then γ is called a path joining u and v. If u=v and the vertices u=v0,v1,,vk1 are distinct, then γ is a cycle. A graph G is called connected if for each pair of distinct vertices u and v, there is a path joining u and v. A tree is a connected graph without any cycle.

Given an n×n symmetric matrix A=(aij)n×n, the graph of A, denoted by G(A), is a graph with vertex set V={1,2,,n} and edge set E={{i,j}:aij0,ij}. Since we are considering graphs without any self loops, so the diagonal entries do not have any role to play in describing the adjacency pattern of the graph corresponding to the matrix. For a graph G=(V,E), where |V|=n, S(G) denotes the set of all n×n symmetric matrices that have G as their graph. Hollow matrix is a matrix whose all diagonal entries are zero. The j×j leading principal submatrix means the j×j submatrix obtained by retaining the first j rows and j columns of the matrix. shows a special tree named broom with vertex set V={1,2,,6} and its corresponding matrix.

Figure 1: A broom and its matrix.

Figure 1: A broom and its matrix.

The expression for the determinant ([Citation9]) of a matrix A is given by, (1) detA=σ(sgnσi=1naiσ(i))(1) where the summation runs over all the n! permutations of {1,2,,n} and sgnσ, the sign of permutation σ, is +1 if the permutation σ is even or 1 if σ is odd.

Lemma 2.1.

For a tree, the determinant of the corresponding matrix A is given by (2) detA=σ(sgnσi=1naiσ(i))(2) where the summation runs over all those permutations of {1,2,,n} which are either a transposition or a product of disjoint transpositions in the cycle decomposition of σ.

Proof.

For any square matrix A, its determinant is given by Equationequation (1). We claim that the product i=1naiσ(i)=0 whenever σ has a cycle of length more than 2. Let some permutation σ in the summation 1 contains a cycle (i1,i2,,ir) of length r>2. Now, if the product term corresponding to this permutation is non-zero, then aikσ(ik)0 for all k. Thus the vertices vi1,vi2,,vir form a cycle in the tree, a contradiction. Hence, the summation runs over all those permutations of {1,2,,n} which are either a transposition or a product of disjoint transpositions only in the cycle decomposition of σ. ▪

Further, given any iV={1,2,,n}, there are two possibilities for any permutation σ of V, either σ(i)=i or σ(i)i. Thus, the summation in Equationequation (2) can be split as (3) detA=σ(i)=i(sgnσi=1naiσ(i))+σ(i)i(sgnσi=1naiσ(i))(3)

For permutations σ with σ(i)=i, the sign of σ restricted to V{i} is same as that of σ. For σ(i)i, aiσ(i)0 if and only if σ(i)i, that is σ(i) is adjacent to i. Thus, we can extract aii from all the terms in the first summation of Equationequation (3). We can also extract aijaji for each ji in the second summation. Let A(i) denote the principal submatrix obtained from A by deleting its ith row and column and A(i,j) is the principal submatrix obtained from A by deleting its i-th and j-th rows and columns. Thus, the expression for the determinant of a matrix whose graph is a tree can be reduced to (4) detA=aiidetA(i)jiaij2detA(i,j)(4)

In this paper, we will solve the following two problems:

Problem 1.

Given a tree on n vertices and an n×n diagonal matrix D and n real numbers λ1,λ2,,λn, find a hollow matrix AS(T) such that each λj is the largest eigenvalue of the j×j leading principal submatrix of A+D.

Problem 2.

Given a tree on n vertices and an n×n diagonal matrix D and n real numbers λ1,λ2,,λn, find a matrix AS(T) such that each λj is the largest eigenvalue of the j×j leading principal submatrix of A+D.

3 Solution to the IEPs

We adopt a standard way of labeling the vertices of the tree in a suitable manner which provides a simplified recurrence relation between the characteristic polynomials of each leading principal submatrix. This scheme of labeling was referred to as successively connected labeling by [Citation20]. For a given tree T, it is always possible to label the vertices of the tree in such a way that each vertex j remains adjacent to exactly one of the vertices from the set {1,2,,j1}. The scheme of labeling is given below :

  1. We choose any vertex of T and label it as 1.

  2. Then we choose any unlabeled vertex which is adjacent to 1 and label it as 2.

  3. Proceeding this way, at the jth step, we label as j any unlabeled vertex which is adjacent to one of the already labeled vertices.

  4. We continue this process till all the vertices of T are labeled from 1 to n.

In the above scheme of labeling, the induced subgraph Tj corresponding to the vertices {1,2,,j} is a tree for each j=1,2,,n. Also, each vertex j is adjacent to exactly one of the vertices of {1,2,,j1}. Let v(j) denote the unique vertex among {1,2,,j1} that is adjacent to j. An example of successively connected labeling of a tree is illustrated in the following :

Figure 2: Successively connected labelling.

Figure 2: Successively connected labelling.

Here, v(2)=v(3)=v(4)=1, v(5)=4, v(6)=5, , v(11)=9.

Following Lemmas will be useful in solving the IEPs:

Lemma 3.1.

[Citation20] Let P(λ) be a monic polynomial of degree n with all real zeros and λmin and λmax be the minimal and maximal zeros of P, then

  1. If μ<λmin, then (1)nP(μ)>0.

  2. If μ>λmax, then P(μ)>0.

Lemma 3.2.

[Citation10] If T is a tree, the largest and smallest eigenvalues of each matrix A in S(T), have multiplicity 1. Moreover, the largest or smallest eigenvalue of a matrix A in S(T) cannot occur as an eigenvalue of a submatrix A(v), for any vertex v of T .

Theorem 3.3.

(Cauchy’s Interlacing Theorem) [Citation9] Let λ1λ2λn be the eigenvalues of any n×n real symmetric matrix A and μ1μ2μn1 be the eigenvalues of any of its submatrix of order (n1)×(n1) then λ1μ1λ2μ2λn1μn1λn

By Cauchy’s Interlacing Theorem 3.3, the largest eigenvalue of any matrix is greater than or equal to the eigenvalues of its submatrices. As an immediate consequence, we have the following necessary condition for the solution of the problems of section 4

Lemma 3.4.

A necessary condition for the existence of solutions to the problems 1 and 2 is λ1λ2λn.

Lemma 3.5.

Let A be a matrix in S(T) corresponding to a given successively connected labeling of T. Then the largest eigenvalues of the leading principal submatrices of A are distinct from each other.

Proof.

Since the tree has successively connected labeling, the induced submatrix corresponding to the vertices {1,2,,j} is a tree for each j=1,2,,n. Thus, Tj is a tree for each j. Then, by Lemma 3.2, the largest eigenvalue of Tj cannot occur as the largest eigenvalue of Tj1. Hence, the largest eigenvalues of each leading principal submatrix of the corresponding matrix of T are distinct. ▪

3.1 Solution to the problem 1

Let T denote a tree with successively connected labeling. A be any hollow matrix whose corresponding graph is T and D=diag(d1,d2,,dn) be the given diagonal matrix. The entries of A are to be computed.

Lemma 3.6.

Let X=A+D and Xj be the j×j leading principal submatrix of X. The sequence Pj(λ)=det(λIjXj) satisfies the following recurrence relation:

  1. P1(λ)=λd1

  2. Pj(λ)=(λdj)Pj1(λ)ajv(j)2Qj(λ)forj=2,3,,n

where Qj(λ) is the characteristic polynomial of the principal submatrix of Xj obtained by deleting the rows and columns indexed by j and v(j), D=diag(d1,d2,,dn) and for convenience of notation, we take Q2(λ)=1.

Proof.

The expression for Pj(λ) can be obtained by expanding the determinant of λIjXj using Equationequation (4). ▪

Theorem 3.7.

For distinct λi‘s, Problem 1 has a solution iff

  1. λ1=d1.

  2. λj>dj, for j=2,3,,n.

Proof.

Since λj is a zero of Pj for each j, so P1(λ1)=0λ1=d1. Also, Pj(λj)=0ajv(j)2=(λjdj)Pj1(λj)Qj(λj)

Here, Qj(λj)0 follows from Cauchy’s interlacing theorem and the fact that each λj is the largest eigenvalue of Xj. By Lemma 3.1, Pj1(λj)>0 and Qj(λj)>0. Thus, ajv(j)2>0 if and only if λjdj>0 i.e. λj>dj. ▪

The above results leads to Algorithm 3.1 which can be compiled into a SCILAB code to get the solutions:

Algorithm 3.1

  1. Input an n×n diagonal matrix D, eigendata λ1<λ2<<λn.

  2. Label the vertices of T following successively connected labeling and input the adjacency pattern of the corresponding matrix A.

  3. If λ1=d1 and λj>dj, for j=2,3,,n Then ajv(j)2=(λjdj)Pj1(λj)Qj(λj)

Else No solution

EndIf

3.2 Solution to the problem 2

Let T denote a tree with successively connected labeling. Let A be any matrix whose corresponding graph is T and D=diag(d1,d2,,dn) be the given diagonal matrix.

Lemma 3.8.

Let X=A+D and Xj be the j×j leading principal submatrix of X. The sequence Pj(λ)=det(λIjXj) satisfies the following recurrence relation:

  1. P1(λ)=λd1a11

  2. Pj(λ)=(λdjajj)Pj1(λ)ajv(j)2Qj(λ)forr=2,3,,n

where Qj(λ) is the characteristic polynomial of the principal submatrix of Xj obtained by deleting the rows and columns indexed by j and v(j), D=diag(d1,d2,,dn).

Proof.

The expression for Pj(λ) can be obtained by expanding the determinant of λIjXj using Equationequation (4). ▪

Theorem 3.9.

Suppose all λi‘s are distinct, then the Problem 2 has infinite numbers of solutions.

Proof.

Since λr is a zero of Pr for each r, P1(λ1)=0a11=λ1d1. Pj(λj)=0ajv(j)2=(λjajjdj)Pj1(λj)Qj(λj)

Here, Qj(λj)0 follows from Cauchy’s interlacing theorem and the fact that each λj is the largest eigenvalue of Xj. By Lemma 3.1, Pj1(λj)>0 and Qj(λj)>0. Thus, ajv(j)2>0 if and only if λjajjdj>0 i.e. λj>ajj+dj. In order to obtain a solution, we choose ajj such that ajj<λjdj. For each ajj<λjdj we get a value of ajv(j)2. So, the above problem has infinitely many solutions. ▪

The above results leads to Algorithm 3.2 which can be compiled into a SCILAB code to get the solutions :

Algorithm 3.2

  1. Input an n×n diagonal matrix D, eigendata λ1<λ2<<λn.

  2. Label the vertices of T following successively connected labeling and input the adjacency pattern of the corresponding matrix A.

  3. Set a11=λ1d1.

  4. For j=2,,n

  5. Choose ajj such that ajj<λjdj

  6. Use ajj to compute ajv(j)2=(λjajjdj)Pj1(λj)Qj(λj)

EndFor

Solution of the above two problems depends on how the vertices of the trees have been labeled. We will show it with an example that two different labeling will give rise to two different solutions.

4 Numerical examples

Numerical solutions to some particular problems has been obtained here using SCILAB programming.

Example 4.1.

Given a tree T on 8 vertices as shown in the figure below, 8 real numbers λ1=1,λ2=3,λ3=4,λ4=5,λ5=6,λ6=7,λ7=8,λ8=9 and the diagonal matrix is D=diag(1,2,3,4,5,6,7,8), find a hollow matrix AS(T) such that λj is the largest eigenvalue of the j×j leading principal submatrix of A+D for each j.

Solution 1.

The vertices of the tree have been labeled in two ways following the scheme of successively connected labeling. For the first way of labeling (), v(2)=1,v(3)=v(4)=v(5)=2,v(6)=v(7)=v(8)=5. The required matrix obtained using SCILAB is (cccccccc01.41421360000001.414213601.15470051.35400641.496291700001.154700500000001.354006400000001.49629170001.18258511.36088461.498144500001.182585100000001.360884600000001.4981445000)

Figure 3: Refer Example 4.1 (First way).

Figure 3: Refer Example 4.1 (First way).

For the second way of labeling (), v(2)=v(6)=v(7)=v(8)=1=1,v(3)=v(4)=v(5)=2. The required matrix obtained using SCILAB is (01.41421360002.30625821.97457821.98387981.414213601.15470051.35400641.496291700001.154700500000001.354006400000001.49629170000002.306258200000001.974578200000001.98387980000000)

Figure 4: Refer Example 4.1 (Second way).

Figure 4: Refer Example 4.1 (Second way).

Example 4.2.

Given the broom T () on 7 vertices, 7 real numbers λ1=3,λ2=5,λ3=9,λ4=10,λ5=13,λ6=18,λ7=20 and the diagonal matrix D=diag(1,4,7,9,11,13,17), find AS(T) such that λj is the largest eigenvalue of the j×j leading principal submatrix of A+D for each j.

Figure 5: Refer Example 4.2.

Figure 5: Refer Example 4.2.

Solution 2.

There are infinitely many solutions of this problem. Here, we show two solutions obtained using SCILAB. As per the labeling, v(2)=1,v(3)=v(4)=v(7)=2,v(5)=4,v(6)=5. The solution obtained by taking a22=a33=0,a44=1,a55=1,a66=3 and a77=1 is A=(21.4142136000001.414213601.89736663.0047581005.441138801.89736660000003.0047581011.9745771000001.974577113.3413312000003.34133123005.441138800001)

The solution obtained by taking a22=1,a33=a44=a55=0,a66=4 and a77=2 is the following matrix A=(2200000213.26598631.6949973003.956709203.26598630000001.6949973002.6954447000002.695444702.4846656000002.48466564003.956709200002)

5 Conclusion

In this paper, we have provide the sufficient condition for the construction of matrix whose corresponding graph is a tree under a specific scheme of labeling. We have dealt with the case where the given eigenvalues are distinct. the case for non-distinct eigenvalues can be studied further. The problem can also be studied for different connected graphs like unicyclic, multicyclic graphs and also with other types of labelings.

References

  • Babaei Zarch, M., Fazeli, S. A. S., Karbassi, S. M. (2018). Inverse eigenvalue problem for matrices whose graph is a banana tree. J. Algorithms Comput. 50(2): 89–101.
  • Babaei Zarch, M., Fazeli, S. A. S. (2019). Inverse eigenvalue problem for a kind of acyclic matrices. Iran. J. Sci. Technol. Trans. A: Sci. 43: 2531–2539.
  • Barcilon, V. (1976). On the solution of the inverse problem with amplitude and natural frequency data, Part I. Phys. Earth Planet. Inter. 13(1): P1–P8.
  • Chu, M. T. (1998). Inverse eigenvalue problems. SIAM Rev. 40(1): 1–39.
  • Gladwell, G. M. L. (1998). Inverse problems in vibration. Appl. Mech. Rev. 39(7): 1013–1018.
  • Heydari, M., Fazeli, S. A. S., Karbassi, S. M. (2019). On the inverse eigenvalue problem for a special kind of acyclic matrices. Appl. Math. 64: 351–366.
  • Hochstadt, H. (1967). On some inverse problems in matrix theory. Archiv der Math. 18: 201–207.
  • Hochstadt, H. (1974). On the construction of a Jacobi matrix from spectral data. Linear Algebra Appl. 8: 435–446.
  • Horn, R. A., Johnson, C. R. (2012). Matrix Analysis, 2nd ed. Cambridge: Cambridge University Press.
  • Johnson, C. R., Duarte, A. L., Saiago, C. M. (2003). Inverse eigenvalue problems and lists of multiplicities of eigenvalues for matrices whose graph is a tree: the case of generalized stars and double generalized stars. Linear Algebra Appl. 373: 311–330.
  • Li, N. (1997). A matrix inverse eigenvalue problem and its application. Linear Algebra Appl. 266: 143–152.
  • Nylen, P., Uhlig, F. (1997). Inverse eigenvalue problems associated with spring-mass systems. Linear Algebra Appl. 254(1): 409–425. Proceeding of the Fifth Conference of the International Linear Algebra Society.
  • Peng, J., Hu, X.-Y., Zhang, L. (2006). Two inverse eigenvalue problems for a special kind of matrices. Linear Algebra Appl. 416(2): 336–347.
  • Pickmann-Soto, H., Arela-Prez, S., Nina, H., Valero, E. (2020). Inverse maximal eigenvalues problems for Leslie and doubly Leslie matrices. Linear Algebra Appl. 592: 93–112.
  • Pickmann, H., Egaña, J., Soto, R. (2009). Two inverse eigenproblems for symmetric doubly arrow matrices. Electron. J. Linear Algebra 18: 700–718.
  • Pickmann, H., Egaa, J., Soto, R. L. (2007). Extremal inverse eigenvalue problem for bordered diagonal matrices. Linear Algebra Appl. 427(2): 256–271.
  • Sharma, D., Sen, M. (2016). Inverse eigenvalue problems for two special acyclic matrices. Mathematics 4(1): 12.
  • Sharma, D., Sen, M. (2018). Inverse eigenvalue problems for acyclic matrices whose graph is a dense centipede. Spec. Matrices 6(1): 77–92.
  • Sharma, D., Sen, M. (2021). The minimax inverse eigenvalue problem for matrices whose graph is a generalized star of depth 2. Linear Algebra Appl. 621: 334–344.
  • Sharma, D., Kumar Sarma, B. (2022). Extremal inverse eigenvalue problem for irreducible acyclic matrices. Appl. Math. Sci. Eng. 30(1): 192–209.
  • Wei, Y., Dai, H. (2016). An inverse eigenvalue problem for the finite element model of a vibrating rod. J. Comput. Appl. Math. 300: 172–182.
  • Xu, W.-R., Bebiano, N., Chen, G.-L. (2019). An inverse eigenvalue problem for pseudo-Jacobi matrices. Appl. Math. Comput. 346: 423–435.
  • Xu, Y.-H., Jiang, E.-X. (2006). An inverse eigenvalue problem for periodic Jacobi matrices. Inverse Problems 23(1): 165–181.