![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
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.
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 and
are adjacent or neighbors, if
. A walk in G, joining vertices u and v, is a sequence
of vertices
such that
is an edge for each
. If the vertices
are distinct, then
is called a path joining u and v. If
and the vertices
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 symmetric matrix
, the graph of A, denoted by
, is a graph with vertex set
and edge set
. 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
, where
,
denotes the set of all
symmetric matrices that have G as their graph. Hollow matrix is a matrix whose all diagonal entries are zero. The
leading principal submatrix means the
submatrix obtained by retaining the first j rows and j columns of the matrix. shows a special tree named broom with vertex set
and its corresponding matrix.
The expression for the determinant ([Citation9]) of a matrix A is given by,
(1)
(1) where the summation runs over all the n! permutations of
and
, the sign of permutation
, is
if the permutation
is even or
if
is odd.
Lemma 2.1.
For a tree, the determinant of the corresponding matrix A is given by
(2)
(2) where the summation runs over all those permutations of
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)(1)
(1) . We claim that the product
whenever
has a cycle of length more than 2. Let some permutation
in the summation 1 contains a cycle
of length
. Now, if the product term corresponding to this permutation is non-zero, then
for all k. Thus the vertices
form a cycle in the tree, a contradiction. Hence, the summation runs over all those permutations of
which are either a transposition or a product of disjoint transpositions only in the cycle decomposition of
. ▪
Further, given any , there are two possibilities for any permutation
of V, either
or
. Thus, the summation in Equationequation (2)
(2)
(2) can be split as
(3)
(3)
For permutations with
, the sign of
restricted to
is same as that of
. For
,
if and only if
, that is
is adjacent to i. Thus, we can extract
from all the terms in the first summation of Equationequation (3)
(3)
(3) . We can also extract
for each
in the second summation. Let
denote the principal submatrix obtained from A by deleting its ith row and column and
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)
(4)
In this paper, we will solve the following two problems:
Problem 1.
Given a tree on n vertices and an diagonal matrix D and n real numbers
, find a hollow matrix
such that each
is the largest eigenvalue of the
leading principal submatrix of
.
Problem 2.
Given a tree on n vertices and an diagonal matrix D and n real numbers
, find a matrix
such that each
is the largest eigenvalue of the
leading principal submatrix of
.
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 . The scheme of labeling is given below :
We choose any vertex of T and label it as 1.
Then we choose any unlabeled vertex which is adjacent to 1 and label it as 2.
Proceeding this way, at the jth step, we label as j any unlabeled vertex which is adjacent to one of the already labeled vertices.
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 corresponding to the vertices
is a tree for each
. Also, each vertex j is adjacent to exactly one of the vertices of
. Let
denote the unique vertex among
that is adjacent to j. An example of successively connected labeling of a tree is illustrated in the following :
Here, ,
,
,
,
.
Following Lemmas will be useful in solving the IEPs:
Lemma 3.1.
[Citation20] Let be a monic polynomial of degree n with all real zeros and
and
be the minimal and maximal zeros of P, then
If
, then
.
If
, then
.
Lemma 3.2.
[Citation10] If T is a tree, the largest and smallest eigenvalues of each matrix A in , have multiplicity 1. Moreover, the largest or smallest eigenvalue of a matrix A in
cannot occur as an eigenvalue of a submatrix
, for any vertex v of T .
Theorem 3.3.
(Cauchy’s Interlacing Theorem) [Citation9] Let be the eigenvalues of any
real symmetric matrix A and
be the eigenvalues of any of its submatrix of order
then
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 .
Lemma 3.5.
Let A be a matrix in 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 is a tree for each
. Thus,
is a tree for each j. Then, by Lemma 3.2, the largest eigenvalue of
cannot occur as the largest eigenvalue of
. 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 be the given diagonal matrix. The entries of A are to be computed.
Lemma 3.6.
Let and
be the
leading principal submatrix of X. The sequence
satisfies the following recurrence relation:
Proof.
The expression for can be obtained by expanding the determinant of
using Equationequation (4)
(4)
(4) . ▪
Theorem 3.7.
For distinct ‘s, Problem 1 has a solution iff
.
, for
.
Proof.
Since is a zero of
for each j, so
. Also,
Here, follows from Cauchy’s interlacing theorem and the fact that each
is the largest eigenvalue of
. By Lemma 3.1,
and
. Thus,
if and only if
i.e.
. ▪
The above results leads to Algorithm 3.1 which can be compiled into a SCILAB code to get the solutions:
Algorithm 3.1
Input an
diagonal matrix D, eigendata
.
Label the vertices of T following successively connected labeling and input the adjacency pattern of the corresponding matrix A.
If
and
, for
Then
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 be the given diagonal matrix.
Lemma 3.8.
Let and
be the
leading principal submatrix of X. The sequence
satisfies the following recurrence relation:
Proof.
The expression for can be obtained by expanding the determinant of
using Equationequation (4)
(4)
(4) . ▪
Theorem 3.9.
Suppose all ‘s are distinct, then the Problem 2 has infinite numbers of solutions.
Proof.
Since is a zero of
for each r,
.
Here, follows from Cauchy’s interlacing theorem and the fact that each
is the largest eigenvalue of
. By Lemma 3.1,
and
. Thus,
if and only if
i.e.
. In order to obtain a solution, we choose
such that
. For each
we get a value of
. 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
Input an
diagonal matrix D, eigendata
.
Label the vertices of T following successively connected labeling and input the adjacency pattern of the corresponding matrix A.
Set
.
For
Choose
such that
Use
to compute
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 and the diagonal matrix is
, find a hollow matrix
such that
is the largest eigenvalue of the
leading principal submatrix of
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 (), . The required matrix obtained using SCILAB is
For the second way of labeling (), . The required matrix obtained using SCILAB is
Example 4.2.
Given the broom T () on 7 vertices, 7 real numbers and the diagonal matrix
, find
such that
is the largest eigenvalue of the
leading principal submatrix of
for each j.
Solution 2.
There are infinitely many solutions of this problem. Here, we show two solutions obtained using SCILAB. As per the labeling, . The solution obtained by taking
and
is
The solution obtained by taking and
is the following matrix
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.