We propose a new method to compute the joint spectral radius and the joint spectral subradius of a set of matrices. We first restrict our attention to matrices that leave a cone invariant. The accuracy of our algorithm, depending on geometric properties of the invariant cone, is estimated. We then extend our method to arbitrary sets of matrices by a lifting procedure, and we demonstrate the eff...

In this work, the convection-diffusion integro-differential equation with a weakly singular kernel is discussed. The  Legendre spectral tau method is introduced for finding the unknown function. The proposed method is based on expanding the approximate solution as the elements of a shifted Legendre polynomials. We reduce the problem to a set of algebraic equations by using operational matrices....

One of unsolved problems in quantum measurement theory is to characterize coexistence of quantum effects. In this paper, applying positive operator matrix theory, we give a mathematical characterization of the witness set of coexistence of quantum effects and obtain a series of properties of coexistence. We also devote to characterizing bijective morphisms on quantum effects leaving the witness...

A finite horizon optimal tracking problem is considered for linear dynamical systems subject to parametric uncertainties in the state-space matrices and exogenous disturbances. suboptimal solution proposed using a model predictive control (MPC) based explicit dual approach, which enables active uncertainty learning. novel algorithm design of robustly invariant online terminal sets controllers p...

In 1957, Chandler Davis proved that unitarily invariant convex functions on the space of hermitian matrices are precisely those which are convex and symmetrically invariant on the set of diagonal matrices. We give a simple perturbation theoretic proof of this result. (Davis’ argument was also very short, though based on completely different ideas). Consider an orthogonally invariant function f ...

Alternative Legendre polynomials (ALPs) are used to approximate the solution of a class of nonlinear Volterra-Hammerstein integral equations. For this purpose, the operational matrices of integration and the product for ALPs are derived. Then, using the collocation method, the considered problem is reduced into a set of nonlinear algebraic equations. The error analysis of the method is given an...

The numerical computation of Lagrangian invariant subspaces of large scale Hamiltonian matrices is discussed in the context of the solution of Lyapunov and Riccati equations. A new version of the low-rank alternating direction implicit method is introduced, which in order to avoid numerical difficulties with solutions that are of very large norm, uses an inverse-free representation of the subsp...

We consider the set of quadruples of matrices defining singular linear time-invariant dynamical systems and show that there is a one-to-one correspondence between this set and a subset of the set of polynomial matrices of degree two. This correspondence preserves the equivalence relations introduced in both sets (feedback-similarity and strict equivalence): two quadruples of matrices are feedba...