Next Article in Journal
Existence Results for a New Class of Boundary Value Problems of Nonlinear Fractional Differential Equations
Previous Article in Journal
Solution of Excited Non-Linear Oscillators under Damping Effects Using the Modified Differential Transform Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Inverse Eigenvalue Problems for Two Special Acyclic Matrices

by
Debashish Sharma
1,*,† and
Mausumi Sen
2,†
1
Department of Mathematics, Gurucharan College, College Road, Silchar 788004, India
2
Department of Mathematics, National Institute of Technology Silchar, Silchar 788010, India
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2016, 4(1), 12; https://doi.org/10.3390/math4010012
Submission received: 17 December 2015 / Revised: 5 February 2016 / Accepted: 26 February 2016 / Published: 3 March 2016

Abstract

:
In this paper, we study two inverse eigenvalue problems (IEPs) of constructing two special acyclic matrices. The first problem involves the reconstruction of matrices whose graph is a path, from given information on one eigenvector of the required matrix and one eigenvalue of each of its leading principal submatrices. The second problem involves reconstruction of matrices whose graph is a broom, the eigen data being the maximum and minimum eigenvalues of each of the leading principal submatrices of the required matrix. In order to solve the problems, we use the recurrence relations among leading principal minors and the property of simplicity of the extremal eigenvalues of acyclic matrices.

1. Introduction

The problems of reconstruction of specially structured matrices from a prescribed set of eigen data are collectively known as inverse eigenvalue problems (IEPs). The level of difficulty of an IEP depends on the structure of the matrices which are to be reconstructed and on the type of eigen information available. M.T. Chu in [1] gave a detailed characterization of inverse eigenvalue problems. A few special types of inverse eigenvalue problems have been studied in [2,3,4,5,6,7,8]. Inverse problems for matrices with prescribed graphs have been studied in [9,10,11,12,13,14]. Inverse eigenvalue problems arise in a number of applications such as control theory, pole assignment problems, system identification, structural analysis, mass spring vibrations, circuit theory, mechanical system simulation and graph theory [1,12,15,16].
In this paper, we study two IEPs, namely IEPP (inverse eigenvalue problem for matrices whose graph is a path) and IEPB (inverse eigenvalue problem for matrices whose graph is a broom). Similar problems were studied in [5], for arrow matrices. The usual process of solving such problems involves the use of recurrence relations among the leading principal minors of λ I - A where A is the required matrix. However, we have included graphs in our analysis by bringing in the requirement of constructing matrices which are described by graphs. In particular, we have considered paths and brooms. Thus, in addition to recurrence relations among leading principal minors, we have used spectral properties of acyclic matrices to solve the problems. Particularly, the strict interlacing of the eigenvalues in IEPB could be proved because of the fact that the minimal and maximal eigenvalues of an acyclic matrix are simple.
The paper is organized as follows : In Section 2, we discuss some preliminary concepts and clarify the notations used in the paper. In Section 3, we define the inverse problems to be studied, namely IEPP and IEPB. Section 4 deals with the analysis of IEPP, the main result being presented as Theorem 4. Section 5 deals with the analysis of IEPB, the main result being presented as Theorem 9. In Section 6, we present some numerical examples to illustrate the solutions of IEPP and IEPB.

2. Preliminary Concepts

Let V be a finite set and let P be the set of all subsets of V which have two distinct elements. Let E P . Then G = ( V , E ) is said to be a graph with vertex set V and edge set E. To avoid confusion, the vertex set of a graph G is denoted by V ( G ) and the edge set is denoted by E ( G ) . Our choice of P implies that the graphs under consideration are free of multiple edges or loops and are undirected. If u , v V and { u , v } E , then we say that u v is an edge and u and v are then called adjacent vertices. The degree of a vertex u is the number of edges which are incident on u. A vertex of degree one is called a pendant vertex. A path P of G is a sequence of distinct vertices v 1 , v 2 , , v n such that consecutive vertices are adjacent. The path on n vertices is denoted by P n . A graph is said to be connected if there exists a path between every pair of its vertices. A cycle is a connected graph in which each vertex is adjacent to exactly two other vertices. A connected graph without any cycles is called a tree.
Given an n × n symmetric matrix A, the graph of A, denoted by G ( A ) , has vertex set V ( G ) = { 1 , 2 , 3 , , n } and edge set { i j : i j , a i j 0 } . For a graph G with n vertices, S ( G ) denotes the set of all n × n symmetric matrices which have G as their graph. A matrix whose graph is a tree is called an acyclic matrix. Some simple examples of acyclic matrices are the matrices whose graphs are paths or brooms (Figure 1).
Throughout this paper, we shall use the following notation:
  • Matrix of a path P n will be a tridiagonal matrix with non zero off-diagonal entries :
    A n = a 1 b 1 0 0 0 b 1 a 2 b 2 0 0 0 b 2 a 3 0 0 0 0 0 a n - 1 b n - 1 0 0 0 b n - 1 a n n × n ,
    where the b i s are non-zero.
  • Matrix of a broom B n + m will be of the following form:
    A n + m = a 1 b 1 0 0 0 0 0 0 b 1 a 2 b 2 0 0 0 0 0 0 b 2 a 3 0 0 0 0 0 0 0 0 a n - 1 b n - 1 0 0 0 0 0 0 b n - 1 a n b n b n + 1 b n + m - 1 0 0 0 0 b n a n + 1 0 0 0 0 0 0 b n + 1 0 a n + 2 0 0 0 0 0 b n + m - 1 0 0 a n + m ( n + m ) × ( n + m ) ,
    where the b i s are non zero.
  • A i will denote the ith leading principal submatrix of the required matrix ( A n or A n + m ).
  • P i ( λ ) = d e t ( λ I i - A i ) , i = 1 , 2 , , n (respectively i = 1 , 2 , n + m ) i.e., the ith leading principal minor of λ I n - A n (respectively λ I n + m - A n + m ), I i being the identity matrix of order i. For the sake of writing the recurrence relations with ease, we define P 0 ( λ ) = 1 , b 0 = 0 .

3. IEPs to be Studied

In this paper we shall study the following two inverse eigenvalue problems :
IEPP Given n real numbers λ j , 1 j n and a real vector X n = ( x 1 , x 2 , , x n ) T find an n × n matrix A n S ( P n ) such that λ j is an eigenvalue of A j , j = 1 , 2 , , n and ( λ n , X n ) is an eigenpair of A n .
IEPB Given 2 n + 2 m - 1 real numbers λ 1 ( j ) , 1 j n + m and λ j ( j ) , 2 j n + m , find an ( n + m ) × ( n + m ) matrix A n + m S ( B n + m ) such that λ 1 ( j ) and λ j ( j ) are respectively the minimal and maximal eigenvalues of A j , j = 1 , 2 , , n + m .

4. Solution of IEPP

The following Lemma gives the relation between successive leading principal minors of λ I n - A n :
Lemma 1. 
The sequence { P j ( λ ) = d e t ( λ I j - A j ) } j = 1 n of characteristic polynomials of A j satisfies the following recurrence relations:
  • P 1 ( λ ) = ( λ - a 1 )
  • P j ( λ ) = ( λ - a j ) P j - 1 ( λ ) - b j - 1 2 P j - 2 ( λ ) , 2 j n .
Here A j denotes the jth leading principal submatrix of A n , the matrix corresponding to the path on n vertices.
Lemma 2. 
For any λ R and 1 j n , P j ( λ ) and P j + 1 ( λ ) cannot be simultaneously zero.
Proof. 
If P 1 ( λ ) = 0 = P 2 ( λ ) , then ( λ - a 2 ) P 1 ( λ ) - b 1 2 = 0 , which implies b 1 = 0 , but this contradicts the restriction on A n that b 1 0 . Once again, for 2 < j n , if P j - 1 ( λ ) = 0 = P j ( λ ) , then the recurrence relation (ii) from Lemma 1, ( λ - a j + 1 ) P j ( λ ) - b j 2 P j - 1 ( λ ) = 0 , which gives P j - 1 ( λ ) = 0 . This will in turn imply that P j - 2 ( λ ) = 0 . Thus, we will end up with P 2 ( λ ) = 0 , implying that b 1 = 0 which is a contradiction. □
Lemma 3. 
If X = ( x 1 , x 2 , , x n ) T is an eigenvector of A n corresponding to an eigenvalue λ, then x 1 0 and the components of this eigenvector are given by
x j = P j - 1 ( λ ) i = 1 j - 1 b i x 1 , j = 2 , 3 , , n .
Proof. 
Since ( λ , X ) is an eigenpair of A n , we have A n X = λ X . Comparing the first n - 1 rows of this matrix equation on both sides, we have
( a 1 - λ ) x 1 + b 1 x 2 = 0 ,
b j - 1 x j - 1 + ( a j - λ ) x j + b j x j + 1 = 0 , j = 2 , , n - 1 .
By the second recurrence relation from Lemma 1,
P j ( λ ) = ( λ - a j ) P j - 1 ( λ ) - b j - 1 2 P j - 2 ( λ ) , j = 2 , 3 , , n .
We define the quantities v 1 , v 2 , , v n as
v 1 = x 1 , v j = x j i = 1 j - 1 b i , 2 j n .
Multiplying Equation (2) by i = 1 j - 1 b i , we get
b j - 1 x j - 1 i = 1 j - 1 b i + ( a j - λ ) x j i = 1 j - 1 b i + b j x j + 1 i = 1 j - 1 b i = 0 b j - 1 2 v j - 1 + ( a j - λ ) v j + v j + 1 = 0 ,
which gives
v j + 1 = ( λ - a j ) v j - b j - 1 2 v j - 1 , j = 2 , 3 , , n - 1 .
Now, from Equation (1), we have v 2 = ( λ - a 1 ) x 1 = x 1 P 1 ( λ ) . Again from Equation (4),
v 3 = ( λ - a 2 ) v 2 - b 1 2 v 1 = x 1 { ( λ - a 2 ) P 1 ( λ ) - b 1 2 } = x 1 P 2 ( λ ) .
Proceeding this way, we see that v j + 1 = x 1 P j ( λ ) , j = 1 , 2 , , n - 1 which can also be written as v j = x 1 P j - 1 ( λ ) , j = 2 , 3 , , n . This further implies that
x j = P j - 1 ( λ ) i = 1 j - 1 b i x 1 , j = 2 , 3 , , n .
Since X is an eigenvector, X 0 . If x 1 = 0 , then from Equation (5), we see that all the other components of X become zero. Thus, x 1 0 . □
Theorem 4. 
The IEPP has a unique solution if and only if x j 0 for all j = 1 , 2 , , n . The unique solution is given by
a 1 = λ 1 , a j = λ j - b j - 1 2 P j - 2 ( λ j ) P j - 1 ( λ j ) , j = 2 , , n
b 1 = x 1 x 2 P 1 ( λ n ) a n d b j - 1 = x 1 P j - 1 ( λ n ) x j i = 1 j - 2 b i , j = 3 , 4 , , n .
Proof. 
Let x j 0 for all j = 1 , 2 , , n . As per the conditions of IEPP, λ j is an eigenvalue of A j for each j = 1 , 2 , , n . Thus, P 1 ( λ 1 ) = 0 a 1 = λ 1 .
P j ( λ j ) = 0 ( λ j - a j ) P j - 1 ( λ j ) - b j - 1 2 P j - 2 ( λ j ) = 0 a j = λ j - b j - 1 2 P j - 2 ( λ j ) P j - 1 ( λ j ) ,
which gives the expression for calculating a j . The expression is valid as P j - 1 ( λ j ) 0 , because by Lemma 2, P j - 1 ( λ j ) and P j ( λ j ) cannot be simultaneously zero.
Now, since ( λ n , X ) is an eigenpair of A n , so by Equation (5), x j = P j - 1 ( λ n ) i = 1 j - 1 b i x 1 , which implies that
b 1 = x 1 x 2 P 1 ( λ n ) and b j - 1 = x 1 P j - 1 ( λ n ) x j i = 1 j - 2 b i .
Since x j 0 hence it follows that , P j - 1 ( λ n ) 0 . Hence the above expression for b j - 1 is valid and b j - 1 0 for all j = 2 , 3 , , n . Successive use of Equations (6) and (7) will give us the values of a j and b j - 1 for j = 1 , 2 , , n .
Conversely, if there exists a unique solution for IEPP, then since X is an eigenvector of A n , so by Lemma 3, x 1 0 . The existence of a solution implies that b j - 1 0 for j = 2 , 3 , , n . It then follows from the expressions in Equation (7) that x j 0 for j = 2 , 3 , , n . □

5. Solution of IEPB

Lemma 5. 
Let P ( λ ) be a monic polynomial of degree n with all real zeros and λ m i n and λ m a x be the minimal and maximal zero of P respectively.
  • If μ < λ m i n , then ( - 1 ) n P ( μ ) > 0 .
  • If μ > λ m a x , then P ( μ ) > 0 .
The proof immediately follows after expressing the polynomial as a product of its linear factors.
Lemma 6. 
If T is a tree, then the minimal and maximal eigenvalues of any matrix A S ( T ) are simple i.e., of multiplicity one. [Corollary 6 of Theorem 2 in [17]]
In other words, this Lemma says that the minimal and maximal eigenvalues of an acyclic matrix are simple. Again, since for each j, the leading principal submatrix A j corresponds to a tree so by Lemma 6 the minimal and maximal eigenvalues of A j must be simple i.e., in particular λ 1 ( j ) λ j ( j ) .
Lemma 7. 
The sequence { P j ( λ ) = d e t ( λ I j - A j ) } j = 1 n + m of characteristic polynomials of A j satisfies the following recurrence relations:
  • P 1 ( λ ) = ( λ - a 1 ) .
  • P j ( λ ) = ( λ - a j ) P j - 1 ( λ ) - b j - 1 2 P j - 2 ( λ ) , 2 j n + 1 .
  • P n + j ( λ ) = ( λ - a n + j ) P n + j - 1 ( λ ) - b n + j - 1 2 P n - 1 ( λ ) i = 1 j - 1 ( λ - a n + i ) , 2 j m .
Lemma 8. 
For any λ R and 1 j n , P j ( λ ) and P j + 1 ( λ ) cannot be simultaneously zero.
Proof. 
Same as Lemma 2. □
By Cauchy’s interlacing theorem ([14,18]), the eigenvalues of a symmetric matrix and those of any of its principal submatrix interlace each other. Thus, λ 1 ( j ) ’s and λ j ( j ) ’s must satisfy:
λ 1 ( n + m ) λ 1 ( n + m - 1 ) λ 1 ( 2 ) λ 1 ( 1 ) λ 2 ( 2 ) λ 3 ( 3 ) λ n + m - 1 ( n + m - 1 ) λ n + m ( n + m ) .
Each diagonal element a i is also a 1 × 1 principal submatrix of A. Hence λ 1 j a i λ j ( j ) , 1 i j . Since λ 1 ( j ) and λ j ( j ) are the minimal and maximal eigenvalues of A j , so P j ( λ 1 ( j ) ) = 0 and P j ( λ j ( j ) ) = 0 . We need to solve these equations successively using the recurrence relations in Lemma 1. For j = 1 , P 1 ( λ 1 ( 1 ) ) = 0 a 1 = λ 1 ( 1 ) . For j = 2 , P 2 ( λ 1 ( 2 ) ) = 0 , P 2 ( λ 2 ( 2 ) ) = 0 which imply that
a 2 = λ 2 ( 2 ) P 1 ( λ 2 ( 2 ) ) - λ 1 ( 2 ) P 1 ( λ 1 ( 2 ) ) P 1 ( λ 2 ( 2 ) ) - P 1 ( λ 1 ( 2 ) ) , b 1 2 = ( λ 2 ( 2 ) - λ 1 ( 2 ) ) P 1 ( λ 1 ( 2 ) ) P 1 ( λ 2 ( 2 ) ) P 1 ( λ 1 ( 2 ) ) - P 1 ( λ 2 ( 2 ) ) .
a 2 will always exist as the denominator in the above expression for a 2 can never be zero. We have λ 1 ( 2 ) λ 2 ( 2 ) and so if P 1 ( λ 2 ( 2 ) ) = P 1 ( λ 1 ( 2 ) ) , then by Rolle’s theorem ∃ c ( λ 1 ( 2 ) , λ 2 ( 2 ) ) such that P 1 ( c ) = 0 1 = 0 , which is not possible. Thus, P 1 ( λ 2 ( 2 ) ) - P 1 ( λ 1 ( 2 ) ) 0 . Also, by Lemma 5, ( - 1 ) 1 P 1 ( λ 1 ( 2 ) ) 0 and so the expression for b 1 2 is non-negative and so we can get real values of b 1 .
Now for 3 j n , we have
P j ( λ 1 ( j ) ) = 0 , P j ( λ j ( j ) ) = 0 ,
which gives
a j P j - 1 ( λ 1 ( j ) ) + b j - 1 2 P j - 2 ( λ 1 ( j ) ) - λ 1 ( j ) P j - 1 ( λ 1 ( j ) ) = 0 ,
a j P j - 1 ( λ j ( j ) ) + b j - 1 2 P j - 2 ( λ j ( j ) ) - λ j ( j ) P j - 1 ( λ j ( j ) ) = 0 .
Let D j denote the determinant of the coefficient matrix of the above system of linear equations in a j and b j - 1 2 . Then D j = P j - 1 ( λ 1 ( j ) ) P j - 2 ( λ j ( j ) ) - P j - 1 ( λ j ( j ) ) P j - 2 ( λ 1 ( j ) ) . If D j 0 , then the system will have a unique solution, given by
a j = λ 1 ( j ) P j - 1 ( λ 1 ( j ) ) P j - 2 ( λ j ( j ) ) - λ j ( j ) P j - 1 ( λ j ( j ) ) P j - 2 ( λ 1 ( j ) ) D j , b j - 1 2 = ( λ j ( j ) - λ 1 ( j ) ) P j - 1 ( λ 1 ( j ) ) P j - 1 ( λ j ( j ) ) D j .
We claim that the expression for b j - 1 2 in RHS is non negative. This follows from Lemma 5. Since λ 1 ( j ) λ 1 ( j - 1 ) and λ j - 1 ( j - 1 ) λ j ( j ) , so by Lemma 5,
( - 1 ) j - 1 D j = ( - 1 ) ( j - 1 ) P j - 1 ( λ 1 ( j ) ) P j - 2 ( λ j ( j ) ) + ( - 1 ) j - 2 P j - 2 ( λ 1 ( j ) ) P j - 1 ( λ j ( j ) ) 0 .
In addition, by Lemma 2, P j ( λ 1 ( j ) ) and P j - 1 ( λ 1 ( j ) ) cannot be simultaneously zero. Thus, P j - 1 ( λ 1 ( j ) ) 0 . Similarly, P j - 1 ( λ j ( j ) ) 0 . This implies that λ 1 ( j ) λ 1 ( j - 1 ) and λ j ( j ) λ j - 1 j - 1 . Thus, we can get non-zero real values of b j - 1 from Equation (8) if and only if λ 1 ( j ) < λ 1 ( j - 1 ) and λ j - 1 ( j - 1 ) < λ j ( j ) for all j = 2 , 3 , , n .
Now, if D j = 0 , then ( - 1 ) j - 1 D j = 0 i.e., ( - 1 ) ( j - 1 ) P j - 1 ( λ 1 ( j ) ) P j - 2 ( λ j ( j ) ) + ( - 1 ) j - 2 P j - 2 ( λ 1 ( j ) ) P j - 1 ( λ j ( j ) ) = 0 . Since both the terms in this sum are non negative, we must have P j - 1 ( λ 1 ( j ) ) P j - 2 ( λ j ( j ) ) = 0 and P j - 2 ( λ 1 ( j ) ) P j - 1 ( λ j ( j ) ) = 0 . However, from Lemma 2, P j - 1 ( λ ) and P j - 2 ( λ ) cannot be simultaneously zero. In addition, P j ( λ ) and P j - 1 ( λ ) cannot be simultaneously zero. Thus, the only possibility is that P j - 2 ( λ 1 ( j ) ) = P j - 2 ( λ j ( j ) ) = 0 . However, this will then imply that a j = λ 1 ( j ) = λ j ( j ) , which is not possible as by Lemma 6 λ 1 ( j ) λ j ( j ) . Thus, D j 0 for all j = 2 , 3 , , n .
Again, λ 1 ( n + j ) and λ n + j ( n + j ) are the eigenvalues of A n + j and so P n + j ( λ 1 ( n + j ) ) = 0 and P n + j ( λ n + j ( n + j ) ) = 0 . Hence,
( λ 1 ( n + j ) - a n + j ) P n + j - 1 ( λ 1 ( n + j ) ) - b n + j - 1 2 P n - 1 ( λ 1 ( n + j ) ) i = 1 j - 1 ( λ 1 ( n + j ) - a n + i ) = 0 , ( λ n + j ( n + j ) - a n + j ) P n + j - 1 ( λ n + j ( n + j ) ) - b n + j - 1 2 P n - 1 ( λ n + j ( n + j ) ) i = 1 j - 1 ( λ n + j ( n + j ) - a n + i ) = 0 .
so we get a system of equations linear in a n + j and b n + j - 1 2 .
a n + j P n + j - 1 ( λ 1 ( n + j ) ) + b n + j - 1 2 P n - 1 ( λ 1 ( n + j ) ) i = 1 j - 1 ( λ 1 ( n + j ) - a n + i ) = λ 1 ( n + j ) P n + j - 1 ( λ 1 ( n + j ) ) , a n + j P n + j - 1 ( λ n + j ( n + j ) ) + b n + j - 1 2 P n - 1 ( λ n + j ( n + j ) ) i = 1 j - 1 ( λ n + j ( n + j ) - a n + i ) = λ n + j ( n + j ) P n + j - 1 ( λ n + j ( n + j ) ) .
We first investigate the conditions under which the coefficient matrix of the above system is singular. By Cauchy’s interlacing property, we have
λ 1 ( n + j ) λ 1 ( n + j - 1 ) λ 1 ( n + 1 ) a n + i λ n + 1 ( n + 1 ) λ n + j - 1 ( n + j - 1 ) λ n + j ( n + j ) , for all i = 1 , 2 , , m .
Thus, i = 1 j - 1 ( λ n + j ( n + j ) - a n + i ) 0 and ( - 1 ) j - 1 i = 1 j - 1 ( λ 1 ( n + j ) - a n + i ) 0 . Let D n + j be the determinant of the coefficient matrix of Equation (10). Then,
( - 1 ) n + j - 1 D n + j = ( - 1 ) n + j - 1 P n + j - 1 ( λ 1 ( n + j ) ) P n - 1 ( λ n + j ( n + j ) ) i + 1 j - 1 ( λ n + j ( n + j ) - a n + i ) , + ( - 1 ) n + j - 2 P n + j - 1 ( λ n + j ( n + j ) ) P n - 1 ( λ 1 ( n + j ) ) i = 1 j - 1 ( λ 1 ( n + j ) - a n + i ) .
As a consequence of Lemma 5, both the products in the LHS are non-negative and so ( - 1 ) n + j - 1 D n + j 0 for all j = 1 , 2 , , m . Thus, D n + j will vanish if and only if ( - 1 ) n + j - 1 D n + j will vanish i.e., if and only if
P n + j - 1 ( λ 1 ( n + j ) ) P n - 1 ( λ n + j ( n + j ) ) i = 1 j - 1 ( λ n + j ( n + j ) - a n + i ) = 0
and
P n + j - 1 ( λ n + j ( n + j ) ) P n - 1 ( λ 1 n + j ) i = 1 j - 1 ( λ 1 ( n + j ) - a n + i ) = 0 .
If P n - 1 ( λ 1 ( n + j ) ) = 0 , then since λ 1 ( n + j ) λ 1 ( n + j - 1 ) λ 1 ( n ) λ 1 ( n - 1 ) and λ 1 ( n - 1 ) is the minimum possible zero of P n - 1 , we get λ 1 n + j = λ 1 ( n + j - 1 ) = = λ 1 ( n ) = λ 1 ( n - 1 ) . Consequently, P n - 1 ( λ 1 ( n + j ) ) = P n ( λ 1 ( n + j ) ) = 0 but this contradicts Lemma 8, according to which P n - 1 ( λ 1 ( n + j ) ) and P n ( λ 1 ( n + j ) ) cannot be simultaneously zero. Hence, P n - 1 ( λ 1 ( n + j ) ) 0 . Similarly, it can be shown that P n - 1 ( λ n + j ( n + j ) ) 0 . Thus there are the following possibilities :
i.
P n + j - 1 ( λ 1 ( n + j ) ) = 0 and P n + j - 1 ( λ n + j ( n + j ) ) = 0 .
ii.
P n + j - 1 ( λ 1 ( n + j ) ) = 0 and i = 1 j - 1 ( λ 1 ( n + j ) - a n + i ) = 0 .
iii.
P n + j - 1 ( λ n + j ( n + j ) ) = 0 and i = 1 j - 1 ( λ n + j ( n + j ) - a n + i ) = 0 .
iv.
i = 1 j - 1 ( λ 1 ( n + j ) - a n + i ) = 0 and i = 1 j - 1 ( λ n + j ( n + j ) - a n + i ) = 0 .
If (i) happens, then, since b n + j - 1 0 , so from the equations in Equation (10), i = 1 j - 1 ( λ 1 ( n + j ) - a n + i ) = 0 and i = 1 j - 1 ( λ n + j ( n + j ) - a n + i ) = 0 . This implies that a n + i = λ 1 ( n + j ) for some i = 1 , 2 , , m and a n + i = λ n + j ( n + j ) for some i = 1 , 2 , , m . However, as per the inequality Equation (11), it then follows that λ 1 ( n + j ) = λ 1 ( n + j - 1 ) = = λ 1 ( n + 1 ) and λ n + 1 ( n + 1 ) = λ n + 2 ( n + 2 ) = = λ n + j ( n + j ) . Since P n + 2 ( λ 1 ( n + 2 ) ) = 0 and P n + 1 ( λ 1 ( n + 1 ) ) = 0 , so the above equalities imply that P n + 2 ( λ 1 ( n + j ) ) = 0 and P n + 1 ( λ 1 ( n + j ) ) = 0 . Hence from the recurrence relation (3) of Lemma 7, we get
( λ 1 ( n + j ) - a n + 1 ) P n + 1 ( λ 1 ( n + j ) ) - b n + 1 2 P n - 1 ( λ 1 ( n + j ) ) ( λ 1 ( n + j ) - a n + 1 ) = 0
which implies that λ 1 ( n + j ) = a n + 1 . Similarly, it will follow that λ n + j ( n + j ) = a n + 1 . However, λ 1 ( n + j ) a n + 1 λ n + j ( n + j ) and so λ 1 ( n + j ) = λ n + j ( n + j ) , but this is not possible as λ 1 ( n + j ) and λ n + j ( n + j ) are the minimal and maximal eigenvalues of the acyclic matrix A n + j and by Lemma 6, the minimal and maximal eigenvalues of an acyclic matrix are simple. Hence (i) cannot hold. From the above arguments, it also follows that (iv) cannot hold.
If (ii) holds, then the augmented matrix of the system of Equation (10) will be of rank one and so the system will have infinite number of solutions. Similarly, if (iii) holds, then the system will have infinite number of solutions. However, if we put the additional constraint that λ 1 ( n + j ) < λ 1 ( n + j - 1 ) and λ n + j - 1 ( n + j - 1 ) < λ n + j ( n + j ) for all j = 2 , 3 , , m then P n + j - 1 ( λ 1 ( n + j ) ) 0 and P n + j - 1 ( λ n + j ( n + j ) ) 0 , so that (ii) and (iii) will not hold.
Thus, we see that D n + j 0 if and only if
λ 1 ( n + j ) < λ 1 ( n + j - 1 ) < < λ 1 ( n + 1 ) < λ n + 1 ( n + 1 ) < < λ n + j - 1 ( n + j - 1 ) < λ n + j ( n + j ) , for all i = 1 , 2 , , m
Under this constraint, the unique solution of the system Equation (10) is given by
a n + j = A j - B j D n + j , b n + j - 1 2 = ( λ n + j ( n + j ) - λ 1 ( n + j ) ) P n + j - 1 ( λ 1 ( n + j ) ) P + j - 1 ( λ n + j ( n + j ) ) D n + j ,
where A j = λ 1 ( n + j ) P n + j - 1 ( λ 1 ( n + j ) ) P n - 1 ( λ n + j ( n + j ) ) i = 1 j - 1 ( λ n + j ( n + j ) - a n + i ) and B j = λ n + j ( n + j ) P n + j - 1 ( λ n + j ( n + j ) ) P n - 1 ( λ 1 ( n + j ) ) i = 1 j - 1 ( λ 1 ( n + j ) - a n + i ) .
b n + j - 1 2 = ( λ n + j ( n + j ) - λ 1 ( n + j ) ) P n + j - 1 ( λ 1 ( n + j ) ) P + j - 1 ( λ n + j ( n + j ) ) D n + j .
The above analysis of the IEP can be framed as the following theorem :
Theorem 9. 
The IEPB has a solution if and only if
λ 1 ( n + m ) < λ 1 ( n + m - 1 ) < < λ 1 ( 2 ) < λ 1 ( 1 ) < λ 2 ( 2 ) < < λ n + m - 1 ( n + m - 1 ) < λ n + m ( n + m )
and the solution is given by
a 1 = λ 1 ( 1 ) , a 2 = λ 2 ( 2 ) P 1 ( λ 2 ( 2 ) ) - λ 1 ( 2 ) P 1 ( λ 1 ( 2 ) ) P 1 ( λ 2 ( 2 ) ) - P 1 ( λ 1 ( 2 ) ) , b 1 2 = ( λ 2 ( 2 ) - λ 1 ( 2 ) ) P 1 ( λ 1 ( 2 ) ) P 1 ( λ 2 ( 2 ) ) P 1 ( λ 1 ( 2 ) ) - P 1 ( λ 1 ( 2 ) ) , a j = λ 1 ( j ) P j - 1 ( λ 1 ( j ) ) P j - 2 ( λ j ( j ) ) - λ j ( j ) P j - 1 ( λ j ( j ) ) P j - 2 ( λ 1 ( j ) ) P j - 1 ( λ 1 ( j ) ) P j - 2 ( λ j ( j ) ) - P j - 1 ( λ j ( j ) ) P j - 2 ( λ 1 ( j ) ) , j = 3 , 4 , , n + 1 , b j - 1 2 = ( λ j ( j ) - λ 1 ( j ) ) P j - 1 ( λ 1 ( j ) ) P j - 1 ( λ j ( j ) ) P j - 1 ( λ 1 ( j ) ) P j - 2 ( λ j ( j ) ) - P j - 1 ( λ j ( j ) ) P j - 2 ( λ 1 ( j ) ) , j = 3 , 4 , , n + 1 , a n + j = A j - B j D n + j , j = 1 , 2 , , m , b n + j - 1 2 = ( λ n + j ( n + j ) - λ 1 ( n + j ) ) P n + j - 1 ( λ 1 ( n + j ) ) P + j - 1 ( λ n + j ( n + j ) ) D n + j , j = 2 , 3 , , m .
The solution is unique except for the signs of the non-zero off-diagonal entries.

6. Numerical Examples

We apply the results obtained in the previous section to solve the following :
Example 1. 
Given 7 real numbers λ 1 = 1 , λ 2 = 5 , λ 3 = - 4 , λ 4 = 3 , λ 5 = 9 , λ 6 = - 8 , λ 7 = - 3 and a real vector X = ( - 2 , 5 , - 7 , 3 , 1 , 4 , 8 ) T , find a matrix A 7 S ( P 7 ) such that λ j is an eigenvalue of A j for each j = 1 , 2 , , 7 and ( λ 7 , X ) is an eigenpair of A 7 .
Solution Using Theorem 1, we obtain the following matrix as the solution :
A 7 = 1 1 . 6 0 0 0 0 0 1 . 6 4 . 36 4 . 8 0 0 0 0 0 4 . 8 - 1 . 0642 - 3 . 4832 0 0 0 0 0 - 3 . 4832 2 . 0515 - 39 . 5368 0 0 0 0 0 - 39 . 5368 - 347 . 5644 115 . 7937 0 0 0 0 0 115 . 7937 - 32 . 6379 0 . 3448 0 0 0 0 0 0 . 3448 - 3 . 1724
The eigenvalues of the all the leading principal submatrices are :
σ ( A 1 ) = { 1 }
σ ( A 2 ) = { 0 . 3600 , 5 . 0000 }
σ ( A 3 ) = { - 4 . 0000 , 0 . 8364 , 7 . 4594 }
σ ( A 4 ) = { - 5 . 2369 , 0 . 5809 , 3 . 0000 , 8 . 0033 }
σ ( A 5 ) = { - 351 . 9801 , - 4 . 7789 , 0 . 7321 , 5 . 8098 , 9 . 0000 }
σ ( A 6 ) = { - 389 . 1678 , - 8 . 0000 , - 3 . 1028 , 0 . 8636 , 7 . 4210 , 18 . 1310 }
σ ( A 7 ) = { - 389 . 1678 , - 8 . 0084 , - 3 . 2694 , - 3 . 0000 , 0 . 8638 , 7 . 4211 , 18 . 1333 }
Example 2. 
Given 13 real numbers 0 . 5 , 1 , - 1 . 4 , 2 , - 2 . 2 , 3 , - 3 . 8 , 4 . 7 , - 4 . 4 , 5 , - 6 , 6 , 7 , rearrange and label them as λ 1 ( j ) , 1 j 7 and λ j ( j ) , 2 j 7 and find a matrix A 4 + 3 S ( B 4 + 3 ) such that λ 1 j and λ j ( j ) are the minimal and maximal eigenvalues of A j , the jth leading principal sub matrix of A 4 + 3 .
Solution Using Theorem 9, we rearrange the numbers in the following way
λ 1 ( 7 ) < λ 1 ( 6 ) < λ 1 ( 5 ) < λ 1 ( 4 ) < λ 1 ( 3 ) < λ 1 ( 2 ) < λ 1 ( 1 ) < λ 2 ( 2 ) < λ 3 ( 3 ) < λ 4 ( 4 ) < λ 5 ( 5 ) < λ 6 ( 6 ) < λ 7 ( 7 )
i.e.,
- 6 < - 4 . 4 < - 3 . 8 < - 2 < - 1 . 4 < 0 . 5 < 1 < 2 < 3 < 4 . 7 < 5 < 6 < 7 .
Then, using the expressions for a j , b j - 1 2 , a n + j and b n + j - 1 2 we get
A 4 + 3 = 1 . 0000 0 . 7071 0 0 0 0 0 0 . 7071 1 . 5000 1 . 9380 0 0 0 0 0 1 . 9380 - 0 . 0047 2 . 1580 0 0 0 0 0 2 . 1580 3 . 3615 1 . 8991 3 . 1884 3 . 7247 0 0 0 1 . 8991 - 3 . 1612 0 0 0 0 0 3 . 1884 0 - 1 . 5823 0 0 0 0 3 . 7247 0 0 - 3 . 1925 .
Here we have taken all the b i s as positive. We can take some of the b i s as negative also. In fact, we can construct 2 6 such matrices for the above problem, the only difference being in the signs of the non-zero off-diagonal entries.
We compute the spectra of all the leading principal submatrices of A 4 + 3 to verify the the conditions of the IEPB are satisfied. The minimal and maximal eigenvalues of each principal submatrix are shown in italics.
σ ( A 7 ) = { - 6 . 0000 , - 3 . 1677 , - 2 . 1477 , - 1 . 4356 , 0 . 8594 , 2 . 8124 , 7 . 0000 }
σ ( A 6 ) = { - 4 . 4000 , - 2 . 6174 , - 1 . 4386 , 0 . 8449 , 2 . 7245 , 6 . 0000 }
σ ( A 5 ) = { - 3 . 8000 , - 1 . 8025 , 0 . 7966 , 2 . 5016 , 5 . 0000 }
σ ( A 4 ) = { - 2 . 0000 , 0 . 7714 , 2 . 3854 , 4 . 7000 }
σ ( A 3 ) = { - 1 . 4000 , 0 . 8953 , 3 . 0000 }
σ ( A 2 ) = { 0 . 5000 , 2 . 0000 }
σ ( A 1 ) = { 1 }

7. Conclusions

The inverse eigenvalue problems discussed in this paper require the construction of specially structured matrices from mixed eigendata. The results obtained here provide an efficient way to construct such matrices from given set of some of the eigenvalues of leading principal submatrices of the required matrix.
The problems IEPP and IEPB are significant in the sense that they are partially described inverse eigenvalue problems i.e., they require the construction of matrices from partial information of eigenvalues and eigenvectors. Such partially described problems may occur in computations involving a complicated physical system where it is often difficult to obtain the entire spectrum. Many times, only the minimal and maximal eigenvalues are known in advance. Thus, the study of inverse problems having such prescribed eigen structure are significant. It would be interesting to consider such IEPs for other acyclic matrices as well.

Acknowledgments

The authors are grateful to the anonymous reviewers for the valuable comments and suggestions.

Author Contributions

Both the authors have contributed equally to the work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chu, M.T. Inverse eigenvalue problems. SIAM Rev. 1998, 40, 1–39. [Google Scholar] [CrossRef]
  2. Ghanbari, K. m-functions and inverse generalized eigenvalue problem. Inverse Prob. 2001, 17, 211–217. [Google Scholar] [CrossRef]
  3. Zhang, Y. On the general algebraic inverse eigenvalue problems. J. Comput. Math. 2004, 22, 567–580. [Google Scholar]
  4. Ehlay, S.; Gladwell, G.M.; Golub, G.H.; Ram, Y.M. On some eigenvector-eigenvalue relations. SIAM J. Matrix Anal. Appl. 1999, 20, 563–574. [Google Scholar]
  5. Peng, J.; Hu, X.-Y.; Zhang, L. Two inverse eigenvalue problems for a special kind of matrices. Linear Algebra Appl. 2006, 416, 336–347. [Google Scholar] [CrossRef]
  6. Ghanbari, K.; Parvizpour, F. Generalized inverse eigenvalue problem with mixed eigen data. Linear Algebra Appl. 2012, 437, 2056–2063. [Google Scholar] [CrossRef]
  7. Pivovarchik, V.; Rozhenko, N.; Tretter, C. Dirichlet-Neumann inverse spectral problem for a star graph of Stieltjes strings. Linear Algebra Appl. 2013, 439, 2263–2292. [Google Scholar] [CrossRef]
  8. Pivovarchik, V.; Tretter, C. Location and multiples of eigenvalues for a star graph of Stieltjes strings. J. Differ. Equ. Appl. 2015, 21, 383–402. [Google Scholar] [CrossRef]
  9. Duarte, A.L. Construction of acyclic matrices from spectral data. Linear Algebra Appl. 1989, 113, 173–182. [Google Scholar] [CrossRef]
  10. Duarte, A.L.; Johnson, C.R. On the minimum number of distinct eigenvalues for a symmetric matrix whose graph is a given tree. Math. Inequal. Appl. 2002, 5, 175–180. [Google Scholar]
  11. Nair, R.; Shader, B.L. Acyclic matrices with a small number of distinct eigenvalues. Linear Algebra Appl. 2013, 438, 4075–4089. [Google Scholar] [CrossRef]
  12. Monfared, K.H.; Shader, B.L. Construction of matrices with a given graph and prescribed interlaced spectral data. Linear Algebra Appl. 2013, 438, 4348–4358. [Google Scholar] [CrossRef]
  13. Sen, M.; Sharma, D. Generalized inverse eigenvalue problem for matrices whose graph is a path. Linear Algebra Appl. 2014, 446, 224–236. [Google Scholar] [CrossRef]
  14. Hogben, L. Spectral graph theory and the inverse eigenvalue problem of a graph. Electron. J. Linear Algebra 2005, 14, 12–31. [Google Scholar] [CrossRef]
  15. Nylen, P.; Uhlig, F. Inverse Eigenvalue Problems Associated With Spring-Mass Systems. Linear Algebra Appl. 1997, 254, 409–425. [Google Scholar] [CrossRef]
  16. Gladwell, G.M.L. Inverse Problems in Vibration; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2004. [Google Scholar]
  17. Johnson, C.R.; Duarte, A.L.; Saiago, C.M. The Parter Wiener theorem: Refinement and generalization. SIAM J. Matrix Anal. Appl. 2003, 25, 352–361. [Google Scholar] [CrossRef]
  18. Horn, R.; Johnson, C.R. Matrix Analysis; Cambridge University Press: New York, NY, USA, 1985. [Google Scholar]
Figure 1. Path P n and Broom B n + m .
Figure 1. Path P n and Broom B n + m .
Mathematics 04 00012 g001

Share and Cite

MDPI and ACS Style

Sharma, D.; Sen, M. Inverse Eigenvalue Problems for Two Special Acyclic Matrices. Mathematics 2016, 4, 12. https://doi.org/10.3390/math4010012

AMA Style

Sharma D, Sen M. Inverse Eigenvalue Problems for Two Special Acyclic Matrices. Mathematics. 2016; 4(1):12. https://doi.org/10.3390/math4010012

Chicago/Turabian Style

Sharma, Debashish, and Mausumi Sen. 2016. "Inverse Eigenvalue Problems for Two Special Acyclic Matrices" Mathematics 4, no. 1: 12. https://doi.org/10.3390/math4010012

APA Style

Sharma, D., & Sen, M. (2016). Inverse Eigenvalue Problems for Two Special Acyclic Matrices. Mathematics, 4(1), 12. https://doi.org/10.3390/math4010012

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop