Elsevier

Linear Algebra and its Applications

Volume 630, 1 December 2021, Pages 204-224
Linear Algebra and its Applications

Generalization of the concept of diagonal dominance with applications to matrix D-stability

https://doi.org/10.1016/j.laa.2021.08.004Get rights and content

Abstract

In this paper, we introduce the class of diagonally dominant with respect to a given LMI region DC matrices. They are shown to possess the analogues of well-known properties of (classical) diagonally dominant matrices, e.g. their spectra are localized inside the region D. Moreover, we show that in some cases, diagonal D-dominance implies (D,D)-stability (i.e. the preservation of matrix spectra localization under multiplication by a positive diagonal matrix).

Introduction

Here, we introduce the following notations:

  • Mn×n for the set of real n×n matrices;

  • Mn×n(C) for the set of complex n×n matrices;

  • σ(A) for the spectrum of a matrix A Mn×n (i.e. the set of all eigenvalues of A defined as zeroes of its characteristic polynomial fA(λ):=det(λIA));

  • C for the open left-hand side of the complex plane C, i.e.C={zC:Re(z)<0};

  • C+ for the closed right-hand side of the complex plane C, i.e.C+={zC:Re(z)0};

  • D+Mn×n for the set of all positive diagonal matrices (i.e. matrices with positive entries on the principal diagonal while the entries outside the principal diagonal are zeroes);

  • D(0,1]+Mn×n for the subclass of D+, defined as follows:D(0,1]+={D=diag{d11,,dnn}:0<dii1,i=1,,n};

  • D1+Mn×n for the subclass of D+, defined as follows:D1+={D=diag{d11,,dnn}:1dii<+,i=1,,n}.

In this paper, we generalize the concept of diagonal dominance. First, let us recall the classical definition of diagonally dominant matrices (see, for example, [7], also [15], [16]).

Definition 1

A matrix AMn×n(C) is called strictly row diagonally dominant if the following inequalities hold:|aii|>ij|aij|i=1,,n. A matrix A is called strictly column diagonally dominant if AT is strictly row diagonally dominant.

Definition 1'. A matrix AMn×n(C) is called generalized diagonally dominant if there exist positive scalars (weights) mi, i=1,,n such thatmi|aii|>ijmj|aij|,i=1,,n (i.e. if there is a positive diagonal matrix M=diag{m1,,mn}, such that AM is strictly row diagonally dominant). If, in addition, aii<0, i=1,,n, then A is called negative diagonally dominant (NDD).

It is well-known (see, for example, [18], p. 382, Theorem 4.C.2) that if a complex matrix A is NDD, then A is Hurwitz stable, i.e. every λσ(A) satisfies Re(λ)<0. The stability of A implies the Lyapunov asymptotic stability of a system of ODE with the system matrix A. Being a sufficient for stability condition, negative diagonal dominance is of particular interest by itself due to its connection to the properties of nonlinear systems (see [17]). Less known is the fact that negative diagonal dominance implies some special concepts of matrix and system stability, that are stronger than just stability and shows stability preservation under specific matrix (system) perturbations. Such concepts include multiplicative D-stability (see, for example, [8], [10]).

Definition 2

A matrix A Mn×n is called (multiplicative) D-stable if Re(λ)<0 for all λσ(DA), where D is any matrix from D+.

Matrix D-stability plays an important role in the theory of stability (see, for example, [14], [9], [12] and references therein) and has numerous applications in the economics, mathematical ecology, mechanics and other branches of science. The proof of the fact that NDD matrices are D-stable is based on the simple observance that the negative diagonal dominance implies stability and is preserved under multiplication by a positive diagonal matrix (see [8], p. 54, Observation (i)). Note, that the property of diagonal dominance can be proved (disproved) by a finite number of steps while the conditions of D-stability of A involve checking all the products of the form DA, where D runs along the infinite set of positive diagonal matrices.

In this paper, we are concerned with the problem of a matrix spectra localization inside a prescribed convex region DC (so called D-stability problem) and its robust aspects. Many problems of the system dynamics lead to establishing D-stability of the system matrix with respect to a specified region DC (see, for example, [5], [6]). The regions of particular interest are LMI (Linear Matrix Inequality) regions, introduced in [3], that include the shifted left half-plane, the unit disk and the conic sector around the negative direction of the real axis. Due to the rapid advance in control theory and its application, we face the problem of finding some easy-to-verify conditions, which allow us to establish:

  • -

    D-stability of a given matrix A;

  • -

    the preservation of D-stability under some specific perturbations of A, e.g. multiplication by a diagonal matrix.

Though an amount of research (see [5], [6]) is due to D-stability (or eigenvalue clustering in a region DC), the progress in describing the matrix classes which are D-stable, and, moreover, preserve D-stability under multiplication by a specific diagonal matrix is still very little and concerns classical regions such as the left half-plane (see [1], [4], [8]) and the unit disk (see [9]). Here, we make certain efforts to fill this gap.

The paper is organized as follows. In Section 2, we study unbounded LMI regions and their properties, focusing on the most well-known examples, such as the shifted left half-plane, the conic sector, hyperbola and parabola. We also recall the concept of (D,D)-stability for an unbounded LMI region D. In Section 3, we introduce the crucial concept of this paper, namely, diagonal dominance with respect to a given LMI region D (so-called diagonal D-dominance). We consider the particular cases of diagonal D-dominance, with respect to the most important LMI regions. Section 4 deals with the basic properties of diagonally D-dominant matrices. Section 5 gives the most important results of the paper, i.e. the implications between diagonal D-dominance, D-stability and (D,D)-stability.

Section snippets

Definition and examples of LMI regions

Consider the following type of regions, introduced in [2] (see also [3]).

Definition 3

A subset DC that can be defined asD={zC:L+Mz+MTz0}, where L,MMn×n, LT=L, is called an LMI region with the characteristic function fD(z)=L+zM+zMT and generating matrices M and L. It was shown in [11], that an LMI region D is open. In the sequel, we shall use the notation D for its closure and D for its boundary.

Later on, we shall be especially interested in the following types of LMI regions.

Example 1

The shifted left

General definition

Given a (bounded or unbounded) LMI region DC, let the intersection of D with the real axis R be denoted by(α,β):=DR, where both α and β may be infinite. According to this notation, β=xmax. Note, that DR is nonempty whenever D is nonempty (see [11], Lemma 21).

Let the part of the boundary of D which lies above the real axis (if D is bounded from above) be represented by a concave function y(x):(α,β)R+. Note, that the boundary function y(x) of an LMI region D, is defined implicitly.

We define

Properties of diagonally D-dominant matrices

Here, we consider the basic properties of diagonally D-dominant matrices, that we shall use later. First, let us state and prove the following technical lemma.

Lemma 1

Let the function r:(R+×R+)R be defined by the formula:r(a,b):=aba2+b2. Then r(a2,b2)r(a1,b1) whenever a2a1, b2b1.

Proof

Since a2a1 and b2b1, we obtain the following inequalities:a12a22(b22b12)+b12b22(a22a12)0;a22b22(a12+b12)a12b12(a22+b22)0;a22b22(a12+b12)a12b12(a22+b22);a22b22a22+b22a12b12a12+b12;a2b2a22+b22a1b1a1+b1. The last

Diagonal D-dominance implies D-stability and (D,D)-stability

First, we recall the following fundamental statement of matrix eigenvalue localization (see, for example, [7], p. 344).

Theorem 1 Gershgorin

Let A={ai,j}i,j=1nMn×n(C), defineRi:=j=1;ijn|aij|,1in. Let D(aii,Ri)C be a closed disk centered at aii with the radius Ri. Then all the eigenvalues of A are located in the union of n discsG(A):=i=1nD(aii,Ri).

A disk D(aii,Ri) defined in the statement of Theorem 1, is called a Gershgorin disk.

Let us recall a well-known fact that if a complex matrix A={aij}i,j=1n is

Conclusions

In this paper, we generalized the concept of diagonally dominance, which is known to be a sufficient condition for stability of a given matrix A. In particular, we addressed the problem of a matrix spectra localization inside a prescribed convex region DC (the so called D-stability problem), with the purpose of finding some easy-to-verify conditions, so as to establish:

  • (i)

    D-stability of a given matrix A;

  • (ii)

    the preservation of D-stability under some specific perturbations of A, e.g. multiplication by

Declaration of Competing Interest

No competing interest.

Acknowledgements

The research was supported by the National Science Foundation of China grant number 12050410229.

References (18)

There are more references available in the full text version of this article.

Cited by (4)

View full text