Next Article in Journal
A Survey on Symmetry Group of Polyhedral Graphs
Previous Article in Journal
Implementing a Novel Use of Multicriteria Decision Analysis to Select IIoT Platforms for Smart Manufacturing
Previous Article in Special Issue
An Efficient Algorithm for Eigenvalue Problem of Latin Squares in a Bipartite Min-Max-Plus System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Spectral Properties of Doubly Stochastic Matrices

1
Department of Mathematics, Sukkur IBA University, Sukkur 65200, Pakistan
2
Department of Mathematics and General Sciences, Prince Sultan University, Riyadh 11586, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(3), 369; https://doi.org/10.3390/sym12030369
Submission received: 7 February 2020 / Revised: 22 February 2020 / Accepted: 27 February 2020 / Published: 2 March 2020
(This article belongs to the Special Issue Symmetry in Numerical Linear and Multilinear Algebra)

Abstract

:
The relationship among eigenvalues, singular values, and quadratic forms associated with linear transforms of doubly stochastic matrices has remained an important topic since 1949. The main objective of this article is to present some useful theorems, concerning the spectral properties of doubly stochastic matrices. The computation of the bounds of structured singular values for a family of doubly stochastic matrices is presented by using low-rank ordinary differential equations-based techniques. The numerical computations illustrating the behavior of the method and the spectrum of doubly stochastic matrices is then numerically analyzed.

1. Introduction

A real stochastic matrix M is a matrix whose row sums or column sums are equal to 1. All the entries of a real stochastic matrix are non-negative. A real symmetric matrix with non-negative entries with row sums and column sums equal to 1 is called doubly stochastic matrix.
A list of n real numbers, i.e., 1 , λ 2 , λ 3 , λ 4 , , λ n is s . d . s realizable if there exists a symmetric doubly stochastic matrix M with its spectrum denoted by σ ( M ) . Doubly stochastic matrix describes the transitions corresponding to finite state symmetric Markov chains and this transition acts as a special class of this family. Doubly stochastic matrices are the convex hull for transition matrices with element set [1].
The inverse eigenvalue problems for non-negative doubly stochastic matrices have its origin in work of [2,3,4,5]. For more details on inverse eigenvalue problems, we refer the reader to [6,7,8,9,10] and references therein.
The non-negative matrix M has a real eigenvalue λ ^ such that λ ^ | λ i ^ | for all i. The eigenvalues λ i ^ for all i are the eigenvalues of M other than λ ^ . From the Perron–Frobenius theorem an eigenvector corresponding to λ ^ is such that each of its entries are non-negative and sums to 1. For more details, we refer to [11,12,13].
[14] showed that the eigenvector has the form x = 1 n ( 1 , 1 , 1 , , 1 ) t corresponding to eigenvalue λ ^ for x R n , where R denotes the real line. The spectrum of a doubly stochastic matrix is bounded by 1, that is λ i 1 for all i. Three important eigenvalue problems for doubly stochastic matrices are considered in [14] whenever there is a possibility that eigenvalues can be placed in complex plane, denoted by C . The first problem deals with necessary and sufficient conditions for n tuples to be the spectrum of a given doubly stochastic matrix. The second problem is about the fact that which real numbers acts as the spectrum of the doubly stochastic matrix. The last problem deals with the study in which set of n real numbers act as the spectrum of the symmetric doubly stochastic inverse eigenvalue problems.
The matrix M over a field F such that each row sum and column sum is λ ^ , is called the generalization of doubly stochastic matrices by means of Lie theory and algebra Ω ( n ) of set of generalized doubly stochastic matrices studied in [15].
In the literature [6,7,10,14,15], much attention has been payed to study the eigenvalues or eigen-spectrum of doubly stochastic matrices. According to the best of our knowledge, however, the study of spectral properties, such as structured singular values for doubly stochastic matrices, is utterly missing from the literature. The current paper deals with the study of the spectral properties, such as singular values and structured singular values for a class of doubly stochastic matrices. The main contribution is to fill this gap which reflects the novelty of our results presented in paper.
The paper is arranged as follows: In Section 2 we provide definitions of symmetric stochastic matrices, singular values and structured singular values. In Section 3 we give a detailed explanation of the computation of singular values for symmetric doubly stochastic matrices. Section 4 of this article contains the geometrical interpretation of eigenvalues and singular values of doubly stochastic matrices. The computation of structured singular values for doubly stochastic matrices is addressed in Section 5. Whereas the numerical experimentation is discussed in Section 6. Section 7 summarizes the conclusions.

2. Preliminaries

Definition 1.
The n dimensional matrix M = ( m i j ) is said to be a doubly stochastic matrix if
( i ) m i j 0 , i , j = 1 : n
( i i ) i = 1 n m i j = 1 , j = 1 : n ; j = 1 n m i j = 1 , i = 1 : n .
Definition 2.
The n dimensional matrix M = ( m i j ) is said to be symmetric doubly stochastic matrix if its transpose is doubly stochastic matrix.
Definition 3.
The singular values of a matrix M are the non-negative real numbers σ i = λ i and σ 1 σ 2 σ n 0 , i = 1 : n .
Definition 4.
The set of block diagonal matrices is denoted by B and is defined as:
B : = { d i a g ( δ 1 I 1 , , δ S I S ; Δ 1 , Δ F ) : δ i C ( R ) , Δ j C m j , m j ( R m j , m j ) i = 1 : S & j = 1 : F } .
In the above definition S and F represent the number of repeated real or complex scalar blocks and the number of full real or complex blocks respectively.
Definition 5.
The structured singular value of M C n , n with respect to set B is denoted by μ B ( M ) and is defined as:
μ B ( M ) = 0 , i f d e t ( I M Δ ) 0 , Δ B 1 m i n Δ { Δ 2 : d e t ( I M Δ ) = 0 , Δ B } o t h e r w i s e .

3. Computing Singular Values of Doubly Stochastic Matrices

Theorem 1.
Let D 1 , D 2 be n dimensional symmetric stochastic matrices with row and column sum equals to 1. Let { σ i } and { σ i ^ } , i = 1 : n are singular values with { u i } and { v i } , i = 1 : n as left and right singular vectors respectively with u i 2 = 1 = v i 2 for { σ i } and { u i ^ } , { v i ^ } are the left and right singular vectors respectively with u i ^ 2 = 1 = v i ^ 2 for { σ i ^ } . The leading singular vectors u 1 and u 2 ^ are orthogonal to { u i } and { v i } for all i = 2 : n , respectively. Let e n = 1 n ( 1 , 1 , 1 , , 1 ) t be the singular vector corresponding to σ 1 and σ 1 ^ then any vector η = { η 1 , η 2 , η 3 , , η n } which is orthogonal to e n is not a singular vector to σ 1 and σ ^ 1 .
Proof. 
The result is proved by expanding and taking sum of all components of the singular value problem of the form:
D η = σ η .
In Equation (1), η is a singular vector corresponding to a singular value σ . The singular vector η is not a singular vector corresponding to singular values σ 1 and σ 1 ^ .
We write
D η = d 11 d 12 d 1 n d n 1 d n 2 d n n η 1 η 2 η n = d 11 η 1 + d 12 η 2 + + d 1 n η n d n 1 η 1 + d n 2 η 2 + + d n n η n .
Taking sum S u m of all components of D η , we have
S u m ( D η ) = S u m d 11 η 1 + d 12 η 2 + + d 1 n η n d n 1 η 1 + d n 2 η 2 + + d n n η n = η 1 S u m ( d i 1 ) + + η n S u m ( d i n ) = η 1 ( σ 1 ) + + η 1 ( σ 1 ) ,
because D is a symmetric stochastic matrix with row and column sums equal to 1. Equation (2) implies that,
S u m ( D η ) = σ 1 ( η 1 + η 2 + + η n ) , = σ 1 S u m ( η ) .
Now, by taking the sum S u m of righthand side of Equation (1) gives,
S u m ( σ η ) = σ S u m ( η ) .
From Equation (3) and (4), we have ( σ σ 1 ) S u m ( η ) = 0 . In turn, this implies that σ σ 1 and S u m ( η ) = 0 , which proves that any vector η = { η 1 , η 2 , η 3 , , η n } orthogonal to vector e n = 1 n ( 1 , 1 , 1 , , 1 ) t which is not a singular vector corresponding to singular values σ 1 and σ 1 ^ . ☐
Theorem 2.
Let D 1 , D 2 be n dimensional symmetric stochastic matrices with row and column sum equals to 1. Let σ i and σ i ^ be singular values corresponding to singular vectors as defined in Theorem 1. The matrix D ˜ is the matrix with D 1 and D 2 along the main diagonal. The singular values of D ˜ do not contain original singular values σ 1 and σ 1 ^ appearing along the main diagonal of a matrix D 1 ˜ with
D 1 ˜ = σ 1 1 1 σ 1 ^ .
The off diagonal of D 1 ˜ contains the rank-1 matrices u v t and v u t consisting of leading left and right singular vectors corresponding to σ 1 and σ 1 ^ , respectively.
Proof. 
We prove the result by computing the singular vectors of the singular value problem of the form:
D ˜ u i = σ i u i , i = 2 : n .
In Equation (6), the set of vectors { u i } for all i = 2 : n . is the singular vectors corresponding to singular values { σ i } i = 2 : n . The vector ( u i , 0 ) t for all i = 2 : n . acts as a singular vector corresponding to σ i i = 2 : n . The vector ( 0 , u i ^ ) t acts as a singular vector corresponding to singular value σ i ^ i = 2 : n . From above discussion it is clear that the vectors of the form ( α i u i , β i u i ^ ) t for all i = 2 : n . act as a singular vectors for matrix D ˜ while the vectors ( α i , β i ) t acts as singular vectors corresponding to singular values σ 1 and σ i ^ of D 1 ˜ .  ☐
Theorem 3.
Let D 1 , D 2 be two symmetric doubly stochastic matrices. Let σ i and σ i ^ be leading singular values corresponding to singular vectors u 1 and u 1 ^ , respectively.
Let
D ˜ = D 1 + ρ I 2 ρ u 1 v 1 t 2 ρ v 1 u 1 t D 2 + ρ I .
with I is an identity matrix with the same dimension as of D 1 and D 2 and ρ is any constant. The singular values of D ˜ does not contains the leading singular values of D 1 and D 2 . The leading singular values are contained in D 1 ˜ with
D 1 ˜ = 3 ρ + σ 1 0 0 ρ σ 1 ^ .
Proof. 
We prove the result by computing singular values of singular value problem,
D ˜ u i = ( σ i + ρ ) u i , i = 2 : n .
The singular values of D ˜ are α , σ 2 + ρ , σ 3 + ρ , , σ n + ρ , β , σ 2 ^ + ρ , σ 3 ^ + ρ , , σ n ^ + ρ where α , β { 3 ρ + σ 1 , ρ σ 1 ^ } which are the singular values of
3 ρ + σ 1 0 0 ρ σ 1 ^ .
Because { u i } i = 2 : n represent the singular vectors corresponding to singular values { σ i } i = 2 : n . can be treated as u i 0 i = 2 : m . Similarly for { σ i ^ } i = 2 : n . the singular vector 0 u i ^ holds true. This show that the singular vector corresponding to D ˜ can be expressed as α i u i β i u i ^ and finally the vectors α i β i i = 2 : n . act as singular vectors to 3 ρ + σ 1 and ρ σ ^ .  ☐
Theorem 4.
A two dimensional symmetric doubly stochastic matrix has singular values 1 and σ 1 if and only if 0 σ 1 1 .
Proof. 
To complete the prove its sufficient to see the fact that for 0 σ 1 1 the matrix
1 + σ 1 2 1 σ 1 2 1 σ 1 2 1 + σ 1 2 ,
is symmetric and doubly stochastic matrix. ☐
Theorem 5.
A three dimensional symmetric doubly matrix has singular values 1 , σ 1 , μ 1 if and only if 0 σ 1 1 , 0 μ 1 1 , σ 1 + 3 μ 1 + 2 0 and 3 σ 1 + μ 1 + 2 0 .
Proof. 
The proof is similar to the one in [15] and hence is omitted. ☐

4. Geometrical Interpretation of Spectrum

In this section, we present the geometrical interpretation of the spectrum of symmetric doubly stochastic matrices. In particular, we discuss the geometry of the eigenspace and singular values and left, right singular vectors of such a class of matrices.
Example 1.
We take doubly stochastic matrices M 1 and M 2 taken from [16] as
M 1 = 0.9 0.1 0.1 0.9 ; M 2 = 0.1 0.2 0.7 0.5 0.3 0.2 0.4 0.5 0.1 .
The spectrum of M 1 , M 2 is shown in Figure 1a,b of Example 1, respectively. The maximum eigenvalue for both M 1 , M 2 is 1 and lies exactly on the spectral circle. The Figure 1c,d show the geometrical interpretation of the singular values and singular vectors obtained for M 1 , M 2 of Example 1, respectively.
Example 2.
We take five dimensional doubly stochastic matrices M 3 , M 4 from [17] as
M 3 = 0 0 0 1 0 0 0 1 2 0 1 2 0 1 2 1 2 0 0 0 1 2 0 0 1 2 1 0 0 0 0 ; M 4 = 0 0 0 1 0 0 0 0.25 0 0.75 0 0.25 0.75 0 0 0 0.75 0 0 0.25 1 0 0 0 0 .
The spectrum of M 3 , M 3 is shown in Figure 2a,b of Example 2, respectively. The maximum eigenvalue for both M 3 , M 3 is 1 and lies exactly on the spectral circle. Figure 2c,d show the geometrical interpretation of the singular values and singular vectors obtained for M 3 , M 4 of Example 2, respectively.
Example 3.
We take eight and nine dimensional doubly stochastic matrices from [18]. The largest eigenvalue corresponding to each matrix attains the maximum value 1. The spectrum is shown in Figure 3a,b of Example 3, respectively. The Figure 3c,d show the geometrical interpretation of the singular values and singular vectors of Example 3, respectively.

5. Computing Structured Singular Values

In this section, our aim is to discuss the spectral properties of doubly stochastic matrices based on the computation of Structured Singular Values (SSV). For this purpose, we compute SSV for a class of matrices as considered in Section 4. SSV is the straight forward generalization of the singular values for the constant matrices. The computation of the exact value of SSV is NP-hard. For this reason, one needs to approximate its bounds, i.e., lower and upper bounds.
From an application point of view, the computation of lower bounds of SSV gives sufficient information about the instability of some feedback system while the upper bounds discuss the stability of feedback system under consideration. The computation of the bounds of SSV presented in this section is based on two powerful mathematical techniques: First technique is based on power method for approximating spectrum [19]. The upper bound of SSV is computed by means of the balanced/AMI technique [20] for computing the upper bound from [21]. The second technique [22] is based on the low rank ODEs-based techniques in order to approximate the lower bounds of SSV. This technique works on a two level algorithm, i.e., inner-outer algorithm. We give a brief description of inner-outer algorithms in the subsequent subsections.

5.1. Inner-Algorithm

The inner-algorithm is used to solve the minimization problem addressed in Equation (5). For this purpose, one needs to construct and then solve a gradient system of ordinary differential equations associated with the optimization problem. The construction of system of ODEs involves the approximation of the local extremizers of structured spectral values sets [22]. Following Theorem 6 helps us to approximate the local extremizer of structured spectral vale sets
Theorem 6.
[22]. For a perturbation Δ B with the block diagonal structure
Δ = { d i a g ( δ 1 I 1 , δ s I s , δ s + 1 I s + 1 , δ S I S ; Δ 1 , , Δ F } ,
with Δ 2 = 1 , acts as a local extremizer of structured spectral value set. For a simple smallest eigenvalue λ = | λ | e ι θ , θ R of matrix valued function ( I ϵ M Δ ) with the right and left eigenvetors x and y scaled as S = e ι θ y x and let z = M y . The non-degeneracy conditions
z k x k 0 , = 1 : S
R e ( z k x k ) 0 , = 1 : S + 1 : S
a n d | | z s + h | | . | | x s + h | | 0 , h = 1 : F ,
hold. Then the magnitude of each complex scalar δ i i = 1 : s appears to be exactly equal to 1 while each full block possesses a unit 2-norm.
Proof. 
For proof we refer to [22].
The system of ordinary differential equations corresponding to a perturbation Δ B is to approximate an extremizer of smallest eigenvalue in magnitude, i.e., λ = | λ | ϵ i θ which is obtained as,
δ i ˙ = ν i ( x i z i R e ( x i z i δ i ¯ ) δ i ) ; i = 1 : s
δ l ˙ = s i g n ( R e ( z l x l ) Ψ ( 1 , 1 ) ( δ l ) ; l = s + 1 : s
Δ j ˙ = ν j ( z s + j x s + j R e Δ j ; z s + j x s + j ) ; j = 1 : F ,
where δ i C , i = 1 : s , δ l R for l = s + 1 and Ψ ( 1 , 1 ) is a characteristic function. For further discussion on the construction of system of ordinary differential equations in above equations, we refer to [22]. ☐

5.2. Outer-Algorithm

The main aim of the outer-algorithm is to vary a small positive parameter ϵ > 0 known as the perturbation level. To vary perturbation level a fast Newton’s iteration was used in [22]. The quantity 1 / ϵ provides the approximation of lower bound to SSV.
The fast Newton’s iteration to solve a problem is
| λ ( ϵ ) | = 1 , ϵ > 0 .
To solve Equation (10), we compute the derivative
d d ϵ | λ ( ϵ ) | .
To compute d d ϵ ( | λ ( ϵ ) | ) , one need following Theorem 7
Theorem 7.
[22]. Let Δ B be the matrix valued function and let x and y as a function of ϵ > 0 are right and left eigenvectors of the perturbed matrix ( ϵ M Δ ) . Consider the scaling of vectors x and y accordingly to Theorem 6. Let z = M y and consider that the non-degenracy conditions as discussed in Theorem 6 holds true, then it yields that
d d ϵ ( | λ ( ϵ ) | ) = 1 | y ( ϵ ) x ( ϵ ) | i = 1 s | z i ( ϵ ) x i ( ϵ ) | + j = 1 F | | z s + j ( ϵ ) | | . | | y s + j ( ϵ ) | | > 0 .
For the proof of the above statement, we refer to [22].

6. Numerical Experimentation

The aim of this section is to present numerical eperimentations for lower bounds of SSV for a class of doubly stochastic matrices obtained in Section 4. The numerical experimentations show that the obtained lower bounds with the help of algorithm [22] are either tighter or equal to one approximated with MATLAB function mussv.
Example 4.
We consider a three dimensional real doubly stochastic valued matrix M 2 taken from [16]
M 2 = 0.1000 0.2000 0.7000 0.5000 0.3000 0.2000 0.4000 0.5000 0.1000 .
We take block uncertainties B = d i a g { δ 1 I 1 , δ 2 I 1 , δ 3 I 1 : δ 1 , δ 3 C , δ 2 R } . The admissible perturbation E is approximated as
E = 1.0000 + 0.0000 i 0 0 0 1.0000 0 0 0 1.0000 + 0.0000 i ,
with E 2 = 1 . The lower bound of SSV by using algorithm [22] is obtained as 1 which is equal to the lower bound approximated by mussv function.
Moreover, by using mussv function, the admissible perturbation is obtained as
= 1.0000 0 0 0 1.0000 0 0 0 1.0000 ,
such that 2 = 1 . The mussv function approximates the same lower and upper bounds of SSV, i.e., 1.
Example 5.
We consider a five dimensional real doubly stochastic valued matrix M 3 taken from [17]
M 3 = 0 0 0 1.0000 0 0 0 0.5000 0 0.5000 0 0.5000 0.5000 0 0 0 0.5000 0 0 0.5000 1.0000 0 0 0 0 .
We take block uncertainties B = d i a g { δ 1 I 1 , δ 2 I 1 , δ 3 I 1 , δ 4 I 1 , δ 5 I 1 : δ 1 , δ 3 , δ 5 C , δ 2 , δ 4 R } . The admissible perturbation E is approximated as
E = 1.0000 + 0.0000 i 0 0 0 0 0 1.0000 0 0 0 0 0 1.0000 + 0.0000 i 0 0 0 0 0 1.0000 0 0 0 0 0 1.0000 + 0.0000 i ,
with E 2 = 1 . The lower bound of SSV by using algorithm [22] is obtained as 1 which is same as the upper bound approximated by mussv function.
Moreover, by using mussv function, the admissible perturbation is obtained as
= 0 0 0 0 0 0 1.2361 0 0 0 0 0 1.2361 0 0 0 0 0 0 0 0 0 0 0 0 ,
such that 2 = 1.2361 Themussvfunction approximates the lower bounds as 0.8090 and upper bound 1 for SSV.
Example 6.
We consider a five dimensional real doubly stochastic valued matrix M 4 taken from [17]
M 4 = 0 0 0 1.0000 0 0 0 0.2500 0 0.7500 0 0.2500 0.7500 0 0 0 0.7500 0 0 0.2500 1.0000 0 0 0 0 .
We take block uncertainties B = d i a g { Δ 1 , δ 1 I 1 , δ 2 I 1 , δ 3 I 1 : Δ 1 C 2 , 2 , δ 1 , δ 3 R , δ 2 C } . The admissible perturbation E is approximated as
E = 0.5000 0.5000 0 0 0 0.5000 0.5000 0 0 0 0 0 1.0000 0 0 0 0 0 1.0000 + 0.0000 i 0 0 0 0 0 1.0000 ,
with E 2 = 1 . The lower bound of SSV by using algorithm [22] is obtained as 1.
Moreover, by using mussv function, the admissible perturbation is obtained as
= 0.6631 0.6028 0 0 0 0.6028 0.5480 0 0 0 0 0 1.2111 0 0 0 0 0 1.2111 0 0 0 0 0 0 ,
such that 2 = 1.2111 The mussv function approximates the lower bounds as 0.8257 and upper bound 1 for SSV.

7. Conclusions

In this article, we presented some useful theorems concerning the spectral properties such as singular values and structured singular values for a class of doubly stochastic matrices. We used low-rank ordinary differential equations-based techniques and MATLAB function mussv to approximate bounds of structured singular values corresponding to doubly stochastic matrices. The numerical experimentations show the behavior of singular values and structured singular values which agree with the fact that the largest value of each singular value and structured singular value for doubly stochastic matrix is bounded above by 1. The obtained results for singular values and structured singular values agree with the results obtained for eigenvalues of doubly stochastic matrices, that is:
  • The doubly stochastic matrix has an eigenvalue 1.
  • The absolute value of any eigenvalue corresponding to a doubly stochastic matrix is less than or equal to 1. The results achieved in this study for structured singular values of doubly stochastic matrices could lead the way to discuss the stability and instability analysis of:
  • Stochastic optimal control systems.
  • Linear feedback systems in control.

Author Contributions

M.-U.R. and A.H. both contributed in conceptualization, methodology, writing manuscript and validation of results, while J.H.B. and J.A. contributed in writing-review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

J. Alzabut would like to thank Prince Sultan University for funding this work through research group Nonlinear Analysis Methods in Applied Mathematics (NAMAM) group with fund number RG-DES-2017-01-17.

Acknowledgments

The authors acknowledge the anonymous reviewers for their potential inputs to increase the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gnacik, M.; Kania, T. Inverse problems for symmetric doubly stochastic matrices whose suleimanova spectra spectra are to be bounded below by 1/2. arXiv 2019, arXiv:1909.01291. [Google Scholar] [CrossRef]
  2. Perfect, H. On positive stochastic matrices with real characteristic roots. In Mathematical Proceedings of the Cambridge Philosophical Society; Cambridge University Press: Cambridge, UK, 1952; Volume 48, pp. 271–276. [Google Scholar]
  3. Perfect, H. Methods of constructing certain stochastic matrices. Duke Math. J. 1965, 69, 35–57. [Google Scholar] [CrossRef]
  4. Suleimanova, H.R. Stochastic matrices with real characteristic values. Dokl. Akad. Nauk SSSR 1949, 66, 343–345. [Google Scholar]
  5. Suleimanova, H.R. The question of a necessary and sufficient condition for the existence of a stochastic matrix with prescribed characteristic numbers. Trudy Vsesojuz. Zaocn. Energet. Inst. Vyp 1965, 28, 33–49. [Google Scholar]
  6. Ccapa, J.; Soto, R.L. On spectra perturbation and elementary divisors of positive matrices. Electron. J. Linear Algebra 2009, 18, 462–481. [Google Scholar] [CrossRef] [Green Version]
  7. Hwang, S.-G.; Pyo, S.-S. The inverse eigenvalue problem for symmetric doubly stochastic matrices. Linear Algebra Appl. 2004, 379, 77–83. [Google Scholar] [CrossRef] [Green Version]
  8. Johnson, C.R.; Paparella, P. Perron spectratopes and the real nonnegative inverse eigenvalue problem. Linear Algebra Appl. 2016, 493, 281–300. [Google Scholar] [CrossRef]
  9. Lei, Y.-J.; Xu, W.-R.; Lu, Y.; Niu, Y.-R.; Gu, X.-M. On the symmetric doubly stochastic inverse eigenvalue problem. Linear Algebra Appl. 2014, 445, 181–205. [Google Scholar] [CrossRef]
  10. Mourad, B.; Abbas, H.; Mourad, A.; Ghaddar, A.; Kaddoura, I. An algorithm for constructing doubly stochastic matrices for the inverse eigenvalue problem. Linear Algebra Appl. 2013, 439, 1382–1400. [Google Scholar] [CrossRef]
  11. Berman, A.; Plemmons, R.J. Non-Negative Matrices in the Mathematical Sciences; SIAM Publications: Philadelphia, PA, USA, 1994. [Google Scholar]
  12. Bhatia, R. Matrix Analysis; Springer: New York, NY, USA, 1997. [Google Scholar]
  13. Seneta, E. Non-Negative Matrices and Markov Chains, 2nd ed.; Springer: New York, NY, USA, 1981. [Google Scholar]
  14. Mourad, B. On a spectral property of doubly stochastic matrices and its application to their inverse eigenvalue problem. Linear Algebra Appl. 2012, 436, 3400–3412. [Google Scholar] [CrossRef] [Green Version]
  15. Mourad, B. Generalized Doubly-Stochastic Matrices and Inverse Eigenvalue Problems; University of New South Wales: Sydney, Australia, 1998. [Google Scholar]
  16. Mehlum, M. Doubly Stochastic Matrices and the Assignment Problem. Master’s Thesis, University of Oslo, Oslo, Norway, 2012. [Google Scholar]
  17. Mashreghi, J.; Rivard, R. On a conjecture about the eigenvalues of doubly stochastic matrices. Linear Multilinear Algebra 2007, 55, 491–498. [Google Scholar] [CrossRef]
  18. Xu, W.-R.; Chen, G.-L. Some results on the symmetric doubly stochastic inverse eigenvalue problem. Bull. Iran. Math. Soc. 2017, 43, 853–865. [Google Scholar]
  19. Young, P.M.; Doyle, J.C. Computation of mu with real and complex uncertainties. In Proceedings of the 29th IEEE Conference on Decision and Control, Honolulu, HI, USA, 7 December 1990; pp. 1230–1235. [Google Scholar]
  20. Young, P.M.; Newlin, M.P.; Doyle, J.C. Practical computation of the mixed μ problem. In Proceedings of the 1992 American Control Conference, Chicago, IL, USA, 26 June 1992; pp. 2190–2194. [Google Scholar]
  21. Fan, M.K.H.; Doyle, J.C.; Tits, A.L. Robustness in the presence of parametric uncertainty and unmodeled dynamics. Adv. Comput. Control. 1989, 130, 363–367. [Google Scholar]
  22. Guglielmi, N.R.; Mutti, U.; Kressner, D. A novel iterative method to approximate structured singular values. SIAM J. Matrix Anal. Appl. 2019, 38, 361–386. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Geometrical interpretation of spectrum and singular values/vectors of Example 1. (a) and (b): spectrum of M 1 , M 2 of Example 1, respectively; (c) and (d): geometrical interpretation of the singular values and singular vectors obtained for M 1 , M 2 of Example 1, respectively.
Figure 1. Geometrical interpretation of spectrum and singular values/vectors of Example 1. (a) and (b): spectrum of M 1 , M 2 of Example 1, respectively; (c) and (d): geometrical interpretation of the singular values and singular vectors obtained for M 1 , M 2 of Example 1, respectively.
Symmetry 12 00369 g001
Figure 2. Geometrical interpretation of spectrum and singular values/vectors of Example 2. (a) and (b): spectrum of M 3 , M 3 of Example 2, respectively; (c) and (d): geometrical interpretation of the singular values and singular vectors obtained for M 3 , M 4 of Example 2, respectively
Figure 2. Geometrical interpretation of spectrum and singular values/vectors of Example 2. (a) and (b): spectrum of M 3 , M 3 of Example 2, respectively; (c) and (d): geometrical interpretation of the singular values and singular vectors obtained for M 3 , M 4 of Example 2, respectively
Symmetry 12 00369 g002
Figure 3. Geometrical interpretation of spectrum and singular values/vectors of Example 3. (a) and (b): spectrum of Example 3, respectively; (c) and (d): geometrical interpretation of the singular values and singular vectors of Example 3, respectively
Figure 3. Geometrical interpretation of spectrum and singular values/vectors of Example 3. (a) and (b): spectrum of Example 3, respectively; (c) and (d): geometrical interpretation of the singular values and singular vectors of Example 3, respectively
Symmetry 12 00369 g003

Share and Cite

MDPI and ACS Style

Rehman, M.-U.; Alzabut, J.; Brohi, J.H.; Hyder, A. On Spectral Properties of Doubly Stochastic Matrices. Symmetry 2020, 12, 369. https://doi.org/10.3390/sym12030369

AMA Style

Rehman M-U, Alzabut J, Brohi JH, Hyder A. On Spectral Properties of Doubly Stochastic Matrices. Symmetry. 2020; 12(3):369. https://doi.org/10.3390/sym12030369

Chicago/Turabian Style

Rehman, Mutti-Ur, Jehad Alzabut, Javed Hussain Brohi, and Arfan Hyder. 2020. "On Spectral Properties of Doubly Stochastic Matrices" Symmetry 12, no. 3: 369. https://doi.org/10.3390/sym12030369

APA Style

Rehman, M. -U., Alzabut, J., Brohi, J. H., & Hyder, A. (2020). On Spectral Properties of Doubly Stochastic Matrices. Symmetry, 12(3), 369. https://doi.org/10.3390/sym12030369

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