Stability and stabilization of linear positive systems on time scales

TL;DR

This paper establishes a criterion for the exponential stability of positive linear systems on time scales based on polynomial coefficients and derives conditions for their positive stabilizability, unifying continuous and discrete cases.

Contribution

It provides a new stability criterion for positive systems on time scales and characterizes positive stabilizability with necessary and sufficient conditions.

Findings

Positive linear systems are exponentially stable iff polynomial coefficients are positive.

Derived necessary and sufficient conditions for positive stabilizability.

Unified stability analysis applicable to continuous and discrete systems.

Abstract

It is shown that a positive linear system on a time scale with a bounded graininess is uniformly exponentially stable if and only if the characteristic polynomial of the matrix defining the system has all its coefficients positive. Then this fact is used to find necessary and sufficient conditions of positive stabilizability of a positive control system on a time scale.

Equations

$$x^{\Delta }(t)=q(t)x(t),x(t_0)=1,$$

$$x^{\Delta }(t)=Ax(t),$$

$$X^{\Delta }(t)=AX(t)$$

$$x^{\Delta }=Ax,x(t_0)=x_0$$

$$x(t)=e_A(t,t_0)x_0.$$

$$x^{\Delta }(t)=Ax(t)+Bu(t),$$

$$\mathrm{spec}(M+\alpha I)=\mathrm{spec}(M)+\alpha ,$$

$$\Vert e_A(t,t_0)\Vert \le K{e}^{-\alpha (t-t_0)}$$

$$x^{\Delta }=\lambda x,$$

$${\chi }_A(\lambda )=(\lambda -{\lambda }_1)\dots (\lambda -{\lambda }_n).$$

$$r(A+\gamma (\mathbb{T})I)=\eta (A)+\gamma (\mathbb{T})<\gamma (\mathbb{T}),$$

$$x(k+1)=\tilde{A}x(k),$$

$$x^{\Delta }(k)=x(k+1)-x(k)=Ax(k),$$

$$x^{\Delta }(kh):=\frac{x((k+1)h)-x(kh)}{h}=Ax(kh).$$

$$\Sigma :x^{\Delta }=Ax+Bu,$$

$$\forall \lambda \in \mathrm{spec}(A):\lambda \notin S(\mathbb{T})\Rightarrow \mathrm{rank}[\lambda I-A,B]=n.$$

$$\tilde{A}:=T^{-1}AT=\left(\begin{array}{ccccccc}0& 1& 0& \cdots & 0& \ast & \\ 0& 0& 1& \cdots & 0& \ast & \\ ⋮& ⋮& ⋮& \ddots & ⋮& ⋮& \\ 0& 0& 0& \dots & 1& \ast & \\ -a_1& -a_2& -a_3& \dots & -a_k& \ast & \\ 0& 0& 0& \dots & 0& A_{22}& \end{array}\right),\tilde{b}:=T^{-1}b=\left(\begin{array}{c}0\\ 0\\ ⋮\\ 0\\ 1\\ 0\end{array}\right),$$

$$k_j\ge {\alpha }_j,j=1,\dots ,n,K{v}_i<a_i,i=1,\dots ,s.$$

$$|(f({\sigma }_{\mathbb{T}}(t))-f(s))-f^{\Delta }(t)({\sigma }_{\mathbb{T}}(t)-s)|\le \epsilon |{\sigma }_{\mathbb{T}}(t)-s|$$

$${\int }_a^{b}f(\tau )\Delta \tau :={\int }_{[a,b{)}_{\mathbb{T}}}f(\tau )\Delta \tau :=F(b)-F(a).$$

Full text

1 Introduction

In positive systems the state variables take nonnegative values. Such systems appear in biology, medicine and economics [1, 2]. We study here linear positive systems on time scales. A time scale is a model of time. Time may be continuous, discrete or mixed – partly continuous and partly discrete. Delta derivative, which is used in delta differential equations that model the positive systems on time scales, may be equal to ordinary derivative or may be equal to a difference quotient, depending on the time scale and a particular point (see Appendix for the precise definitions).

It is known that in the continuous-time case the positive system $\dot{x}=Ax$ is exponentially stable if and only if the characteristic polynomial of the matrix $A$ has all its coefficients positive (see [1, 2]). We show that this is true for a system on an arbitrary time scale as long as the graininess function of the time scale (which measures the distance between a particular time instant and the next instant) is bounded and uniform exponential stability is considered. Thus any discretization (in particular nonuniform) of a positive continuous-time uniformly exponentially stable system that preserves positivity gives an uniformly exponentially stable system. We rely here on [3] and [4], where uniform exponential stability of linear systems was studied. Stability of linear systems on time scales was also considered in [5, 6, 7, 8, 9, 10, 11, 12, 13, 14]. Various authors considered different concepts of stability and different variants of the same concept. For example, in the definition of exponential stability, to estimate the solutions, either the standard exponential function was used or the exponential function on the time scale.

If a control system is not uniformly exponentially stable, we can try to use feedback to stabilize the system. If our system is positive, it is natural to require that the feedback preserves positivity. This procedure is called positive stabilization. We give necessary and sufficient conditions for positive stabilizability of positive systems on time scales with bounded graininess. One of the conditions is equivalent to standard stabilizability of a system on a time scale. The other consists of inequalities that must be satisfied by the matrix that defines the feedback.

Stabilizability of positive linear systems were investigated in e.g. [1, 2], separately for continuous-time and discrete-time systems. Other problems for positive linear systems on time scales like reachability, observability and realizations were studied in [15, 16, 17, 18, 19, 20].

2 Preliminaries

We shall study control systems on time scales. A short introduction to the calculus on time scales is provided in Appendix. More information can be found in e.g. [21].

Let $\mathbb{T}$ be a time scale with the forward graininess $\mu_{\mathbb{T}}$. We shall assume that $\sup\mathbb{T}=+\infty$ and denote $\mu_{\mathbb{T}}:=\sup\{\mu_{\mathbb{T}}(t):t\in\mathbb{T}\}$.

We shall write $e_{q}(t,t_{0})$ to denote the generalized exponential function of $q$, initialized at $t_{0}$ and evaluated at $t$. It is the unique forward solution ($t\geq t_{0}$) of the initial value problem

[TABLE]

where $x^{\Delta}(t)$ means the delta derivative of the function $x$ at time $t\in\mathbb{T}$.

For a square matrix $A$ with complex or real elements, $\operatorname{spec}(A)$ will mean its spectrum, i.e. the set of its eigenvalues and $\|A\|$ will mean the standard spectral norm of $A$.

Let $A$ be a real ${n\times n}$ matrix. Consider the time-invariant linear system on the time scale $\mathbb{T}$

[TABLE]

where $x(t)\in\mathbb{R}^{n}$ and $t\in\mathbb{T}$.

Theorem 2.1** ([21]).**

Let $t_{0}\in\mathbb{T}$ and $x_{0}\in\mathbb{R}^{n}$. Then system (1) with the initial condition $x(t_{0})=x_{0}$ has a unique solution $x:[t_{0},+\infty)\cap\mathbb{T}\rightarrow\mathbb{R}^{n}$.

This result can be extended to matrix-valued solutions of (1), which leads to the following definition.

Let $t_{0}\in\mathbb{T}$. A function $X:[t_{0},+\infty)\cap\mathbb{T}\rightarrow\mathbb{R}^{n\times n}$ that satisfies the matrix delta differential equation

[TABLE]

and the initial condition $X(t_{0})=I$, where $I$ is the $n\times n$ identity matrix, is called the matrix exponential function (corresponding to $A$) initialized at $t_{0}$. Its value at $t\in\mathbb{T}$, $t\geq t_{0}$, is denoted by $e_{A}(t,t_{0})$. Then the solution of the initial value problem

[TABLE]

can be written as

[TABLE]

All the definitions and statements of this section can be naturally extended to complex-valued functions and matrices.

3 Positive control systems on time scales

Let $\mathbb{R}^{n}_{+}$ mean the set of all vectors from $\mathbb{R}^{n}$ with nonnegative components.

Consider now a linear control system on a time scale $\mathbb{T}$:

[TABLE]

where $t\in\mathbb{T}$, $x(t)\in\mathbb{R}^{n}$ and $u(t)\in\mathbb{R}^{m}$. We assume that the control $u$ is a piecewise continuous function of time.

Definition 3.1**.**

The system (3) is positive if for all $t_{0}\in\mathbb{T}$ and any control $u:[t_{0},t_{1}]_{\mathbb{T}}\rightarrow\mathbb{R}^{m}_{+}$ the trajectory starting from any $x(t_{0})=x_{0}\in\mathbb{R}^{n}_{+}$ and corresponding to $u$ stays in $\mathbb{R}^{n}_{+}$ for all $t\in[t_{0},t_{1}]_{\mathbb{T}}$.

We extend $\mathbb{R}$ adding $+\infty$: $\bar{\mathbb{R}}:=\mathbb{R}\cup\{+\infty\}$. For $a\in\mathbb{R}$ we set $a+\infty:=+\infty$, $1/0=+\infty$, $1/+\infty=0$ and $+\infty>a$.

Let $\gamma:=\gamma(\mathbb{T}):={1}/{\bar{\mu}_{\mathbb{T}}}$. Then $\gamma(\mathbb{R})=+\infty$, $\gamma(h\mathbb{Z})=1/h$ and $\gamma(2^{\mathbb{N}})=0$.

Let $A_{\mathbb{T}}:=A+\gamma(\mathbb{T})I$, where $\gamma(\mathbb{T})I$ is the diagonal matrix with $\gamma(\mathbb{T})$ on the diagonal.

Proposition 3.2** ([16, 4]).**

The system $x^{\Delta}=Ax+Bu$ is positive if and only if all elements of $A_{\mathbb{T}}$ and $B$ are nonnegative.

The condition that $A_{\mathbb{T}}$ is nonnegative implies that $A$ is a Metzler matrix, i.e. its off-diagonal elements are nonnegative. Let $M$ be an $n\times n$ Metzler matrix. Let us set $c(M):=\min\{a\geq 0:A+aI\geq 0\}$. From the definition of $c(M)$ we get:

Proposition 3.3**.**

$A_{\mathbb{T}}$* is nonnegative if and only if $A$ is Metzler and $c(A)\leq\gamma(\mathbb{T})$.*

For a real $n\times n$ matrix $M$ let $r(M):=\max\{|\lambda|:\lambda\in\operatorname{spec}(M)\}$ be the spectral radius of $M$ and $\eta(M):=\max\{\operatorname{Re}\lambda:\lambda\in\operatorname{spec}(M)\}$ be the spectral abscissa.

Theorem 3.4** (Perron-Frobenius).**

Let $P\in\mathbb{R}^{n\times n}_{+}$. Then $r(P)$ is an eigenvalue of $P$.

The proof of Theorem 3.4 can be found e.g. in [22].

Proposition 3.5** ([4]).**

Let $M$ be a real $n\times n$ Metzler matrix. Then $\eta(M)=r(M+\alpha I)-\alpha$ for all $\alpha\geq c(M)$, and $\eta(M)$ is an eigenvalue of $M$.

Proof.

Since $M+\alpha I\in\mathbb{R}^{n\times n}_{+}$ for $\alpha\geq c(M)$, then Theorem 3.4 says that $r(M+\alpha I)$ is an eigenvalue of $M+\alpha I$. Thus $\eta(M+\alpha I)=r(M+\alpha I)$. But for any $\alpha\in\mathbb{R}$

[TABLE]

so $\eta(M+\alpha I)=\eta(M)+\alpha$. Thus $r(M+\alpha I)=\eta(M)+\alpha$ and this gives the required equality. In particular, taking $\alpha=c(M)$ we get $\eta(M)=r(M+c(M)I)-c(M)$. Moreover, since $r(M+c(M)I)$ is an eigenvalue of $M+c(M)I$, then, by (4), $r(M+c(M)I)-c(M)$ is an eigenvalue of $M$. ∎

4 Stability

We consider here uniform exponential stability.

Definition 4.1**.**

([3]) System $x^{\Delta}=Ax$ is uniformly exponentially stable if there are constants $K\geq 1$ and $\alpha>0$, and an open neighborhood $V$ of [math] in $\mathbb{R}^{n}$ such that for every $t_{0},t\in\mathbb{T}$ with $t\geq t_{0}$ and every $x_{0}\in\mathbb{R}^{n}\cap V$, the forward trajectory $x$ of the system, corresponding to the initial condition $x(t_{0})=x_{0}$, satisfies $\|x(t)\|\leq K\exp(-\alpha(t-t_{0}))\|x_{0}\|$.

It is known that the condition $\sup\{\mu_{\mathbb{T}}(t):t\in\mathbb{T}\}<+\infty$ is necessary for uniform exponential stability of the system $x^{\Delta}=f(x)$ (see e.g. [23]). Therefore, from now on we shall assume that $\bar{\mu}_{\mathbb{T}}:=\sup\{\mu_{\mathbb{T}}(t):t\in\mathbb{T}\}<+\infty$.

For linear systems on time scales uniform exponential stability has been well investigated. The following proposition follows directly from the definition.

Proposition 4.2** ([3]).**

*The following conditions are equivalent:

(i) system $x^{\Delta}=Ax$ is uniformly exponentially stable,

(ii) there exist constants $\alpha>0$ and $K\geq 1$ such that for all $t_{0}\in\mathbb{T}$*

[TABLE]

for all $t\in\mathbb{T}$ such that $t\geq t_{0}.$

We will say that $\lambda\in\mathbb{C}$ is uniformly exponentially stable if the scalar equation

[TABLE]

where $x\in\mathbb{C}^{n}$, is uniformly exponentially stable.

The set of all uniformly exponentially stable $\lambda\in\mathbb{C}$ will be denoted by $\mathcal{S}_{\mathbb{T}}$. It depends on the time scale $\mathbb{T}$. For $\mathbb{T}=\mathbb{R}$ it is equal to $\mathbb{C}_{-}$, while for $\mathbb{T}=h\mathbb{Z}$ it is the open disc of the radius $1/h$ with the center $-1/h$. But for other time scales such sets are often unknown. More details can be found in [13, 3].

In [4] the following has been shown:

Theorem 4.3**.**

Let $\mu(\mathbb{T})<+\infty$. Then $\mathcal{S}_{\mathbb{T}}\subset\mathbb{C}_{-}$ and for $\gamma=\gamma(\mathbb{T})<+\infty$ the set $\mathcal{S}_{\mathbb{T}}$ contains the open disc of the radius $\gamma$ and the center at $-\gamma$. For $\gamma(\mathbb{T})=+\infty$, $\mathcal{S}_{\mathbb{T}}=\mathbb{C}_{-}$.

For linear systems uniform exponential stability can be characterized by the eigenvalues of the matrix $A$.

Theorem 4.4** ([3]).**

*The following conditions are equivalent:

(i) system $x^{\Delta}=Ax$ is uniformly exponentially stable,

(ii) every $\lambda\in\operatorname{spec}(A)$ is uniformly exponentially stable.*

Corollary 4.5**.**

System $x^{\Delta}=Ax$ is uniformly exponentially stable if and only if $\operatorname{spec}(A)\subset\mathcal{S}_{\mathbb{T}}$.

For positive systems we are interested in what happens with the trajectories that start from initial states belonging to $\mathbb{R}^{n}_{+}$. This leads to the following definition.

Definition 4.6**.**

Assume that system $x^{\Delta}=Ax$ is positive. We say that $x^{\Delta}=Ax$ is positively uniformly exponentially stable if there are constants $K\geq 1$ and $\alpha>0$, and an open neighborhood $V$ of [math] in $\mathbb{R}^{n}$ such that for every $t_{0},t\in\mathbb{T}$ with $t\geq t_{0}$ and every $x_{0}\in\mathbb{R}^{n}_{+}\cap V$, the forward trajectory $x$ of the system, corresponding to the initial condition $x(t_{0})=x_{0}$, satisfies $\|x(t)\|\leq K\exp(-\alpha(t-t_{0}))\|x_{0}\|$ for all $t\in\mathbb{T}$ such that $t\geq t_{0}$.

It appears that for linear positive systems both properties coincide.

Proposition 4.7** ([24]).**

A positive system $x^{\Delta}=Ax$ is positively uniformly exponentially stable if and only if it is uniformly exponentially stable.

Proposition 4.7 implies that for positive linear systems positive uniform exponential stability can be characterized with the aid of the spectrum of the matrix $A$. However this spectrum for a positive system has a specific structure, so other tools for checking (positive) uniform exponential stability can be employed.

Theorem 4.8**.**

A linear positive system $x^{\Delta}=Ax$ is positively uniformly exponentially stable if and only if all the coefficients of the characteristic polynomial $\chi_{A}$ of the matrix $A$ are positive.

Proof.

Necessity. Assume that the system $x^{\Delta}=Ax$ is positively uniformly exponentially stable. From Proposition 4.7, Theorem 4.4 and Theorem 4.3 it follows that all eigenvalues $\lambda_{1},\ldots,\lambda_{n}$ of $A$ have negative real parts. Then

[TABLE]

If $\lambda_{i}$ is complex then $\lambda_{i}=-a+bi$ for $a>0$. Since $A$ is real, $\lambda_{j}=-a-bi$ for some $j\neq i$. Then the polynomial $(\lambda-\lambda_{i})(\lambda-\lambda_{j})=\lambda^{2}+2a\lambda+b^{2}$ has all its coefficients positive. If $\lambda_{k}$ is real, then $\lambda_{k}=-c$ for $c>0$, so the polynomial $\lambda-\lambda_{k}=\lambda+c$ has all its coefficients positive as well. Since $\chi_{A}(\lambda)$ is a product of polynomials of these two types, it also has all its coefficients positive.

Sufficiency. Let us assume that $\chi_{A}(\lambda)=\lambda^{n}+a_{n-1}\lambda^{n-1}+\ldots+a_{0}$ has all its coefficients positive. Since $A$ is a Metzler matrix, from Proposition 3.5 it follows that $\lambda_{1}:=\eta(A)\in\mathbb{R}$ is an eigenvalue of $A$. If it were nonnegative, then $\chi_{A}(\lambda_{1})$ would be greater than [math], and this would contradict the fact that $\lambda_{1}$ is an eigenvalue of $A$. Thus $\eta(A)<0$, which means that all eigenvalues of $A$ have negative real parts. Assume now that $\gamma(\mathbb{T})<+\infty$. Since $c(A)\leq\gamma(\mathbb{T})$ (from positivity of the system and Proposition 3.3), using again Proposition 3.5 we get

[TABLE]

which implies that $\operatorname{spec}(A+\gamma(\mathbb{T})I)$ is contained in the disc of the radius $\gamma(\mathbb{T})$ and the center at [math]. By (4), $\operatorname{spec}(A)=\operatorname{spec}(A+\gamma(\mathbb{T})I)-\gamma(\mathbb{T})$, so $\operatorname{spec}(A)$ is contained in in the disc of the radius $\gamma(\mathbb{T})$ and the center at $-\gamma(\mathbb{T})$. From Theorem 4.3 we get that $\operatorname{spec}(A)\subset\mathcal{S}(\mathbb{T})$, which, by Corollary 4.5, implies that the system is uniformly exponentially stable. For $\gamma(\mathbb{T})=+\infty$, $\mathcal{S}(\mathbb{T})=\mathbb{C}_{-}$, so $\operatorname{spec}(A)\subset\mathcal{S}(\mathbb{T})$ as well. ∎

Remark 4.9*.*

This fact has long been known for continuous-time systems [1, 2], when $\mathbb{T}=\mathbb{R}$. For a discrete-time positive system of the form

[TABLE]

(uniform) exponential stability has been characterized by the condition that the characteristic polynomial of $A:=\tilde{A}-I$ has positive coefficients [1, 2]. But (7) is equivalent to the system on the time scale $\mathbb{T}=\mathbb{Z}$:

[TABLE]

so this characterization of uniform exponential stability agrees with the characterization presented in Theorem 4.8 for systems on arbitrary time scales with the bounded graininess. In particular, this characterization is valid for $\mathbb{T}=h\mathbb{Z}$, where $h>0$. Consider an Euler discretization of the positive continuous-time system $\dot{x}=Ax$ with the step $h$:

[TABLE]

Since the matrix $A$ is the same for both systems, from Theorem 4.8 we conclude that uniform exponential stability of the continuous-time system implies uniform exponential stability of the discretized system, provided the latter is positive. This holds if and only if $h\leq 1/c(A)$.

5 Stabilization

The main goal of this section is to find conditions for feedback stabilization of a positive linear system on a time scale $\mathbb{T}$:

[TABLE]

with $x(t)\in\mathbb{R}^{n}$, $u(t)\in\mathbb{R}^{m}$. However we want to preserve positivity of the system.

Definition 5.1**.**

System 8 is positively stabilizable if there is a feedback $u=Kx$, such that the closed-loop system $x^{\Delta}=(A+BK)x$ is positive and positively uniformly exponentially stable.

Remark 5.2*.*

We do not assume that $K$ is nonnegative, so $u=Kx$ may have negative components, even for $x\in\mathbb{R}^{n}_{+}$. This may be interpreted as enlarging the set of admissible control values and making it depend on the state $x$. But we still have nonnegativity of trajectories of the system that start from the points of $\mathbb{R}^{n}_{+}$ and correspond to controls defined by the feedback. On the other hand, if the open loop system is not positively uniformly exponentially stable, then by applying the feedback $u=Kx$ with a nonnegative $K$ we cannot stabilize the system (see [2]).

Since positive stabilizability implies usual stabilizability (where we do not require that the closed-loop system is positive), we have the following:

Proposition 5.3**.**

If the positive system (8) is positively stabilizable, then the following condition holds:

[TABLE]

Condition (9) is equivalent to standard stabilizability [10], i.e. it does not guarantee existence of $K$ such that the system $x^{\Delta}=(A+BK)x$ is uniformly exponentially stable and positive at the same time. To find necessary and sufficient condition for positive stabilizability we need to add another condition to (9). For simplicity we assume now that $m=1$, so $B=b=(b_{1},\ldots,b_{n})^{T}$ and $K=(k_{1},\ldots,k_{n})$

Let us define for $j=1,\ldots,n$:

$\alpha_{j}:=-\infty$ if $\gamma(\mathbb{T})=+\infty$ and for every $i\neq j$, $b_{i}=0$,

$\alpha_{j}:=\max_{i\neq j,b_{i}\neq 0}\{\frac{-a_{ij}}{b_{i}}\}$, if $b_{j}=0$ or $\gamma(\mathbb{T})=+\infty$ and there is $i\neq j$ such that $b_{i}\neq 0$,

$\alpha_{j}:=\max\{\max_{i\neq j,b_{i}\neq 0}\{\frac{-a_{ij}}{b_{i}}\},\frac{-a_{jj}-\gamma(\mathbb{T})}{b_{j}}\}$ , otherwise.

Observe that $\alpha_{j}\leq 0$ for any $j=1,\ldots,n$.

Now assume that $\operatorname{rank}[b,Ab,\ldots,A^{n-1}b]=k$. Then $A^{k}b=-a_{1}b-\ldots-a_{k}A^{k-1}b$ for some $a_{1},\ldots,a_{k}\in\mathbb{R}$. Let us define a basis of $\mathbb{R}^{n}$ as follows:

$v_{k-i}=A^{i}b+a_{k}A^{i-1}b+\ldots+a_{k-i+1}b$ for $i=0,\ldots,k-1$, and, if $k<n$, $v_{k+1},\ldots,v_{n}$ chosen arbitrarily so that $v_{1},\ldots,v_{n}$ are linearly independent. In particular, $v_{k}=b$.

Then for $i=0,\ldots,k-2$, $Av_{k-i}=v_{k-i-1}-a_{k-i}b$ and $Av_{1}=a_{1}b$. Thus, letting $T=(v_{1},\ldots,v_{n})$ we get

[TABLE]

where $A_{22}$ is a $(n-k)\times(n-k)$ matrix and [math]’s in the last rows of $\tilde{A}$ and $\tilde{b}$ mean zero $(n-k)\times 1$ matrices.

Now we can state a characterization of positive stabilizability of a positive system.

Theorem 5.4**.**

The positive system (8) is positively stabilizable if and only if condition (9) is satisfied and the following set of linear inequalities for $K=(k_{1},\ldots,k_{n})$ is consistent:

[TABLE]

Proof.

Necessity. From Proposition 5.3 it follows that condition (9) is necessary for positive stabilizability of system (8). Moreover, positive stabilizability implies existence of $K=(k_{1},\ldots,k_{n})$ such that the closed-loop system $x^{\Delta}=(A+bK)x$ is positive. Thus such $K$ must satisfy $a_{ij}+b_{i}k_{j}\geq 0$ for $i\neq j$ and $a_{jj}+b_{j}k_{j}+\gamma(\mathbb{T})\geq 0$ for $j=1,\ldots,n$. If $\gamma(\mathbb{T})=+\infty$, the second inequality is always satisfied. Easy calculation shows that these inequalities are equivalent to the inequalities $k_{j}\geq\alpha_{j}$, $j=1,\ldots,n$. To show that the inequalities $Kv_{i}\leq a_{i}$ are also necessary let us transform system (8) to $\tilde{x}^{\Delta}=\tilde{A}\tilde{x}+\tilde{b}u$, where $\tilde{A}$ and $\tilde{b}$ are given by (10). Positive stabilizability of system (8) implies existence of $K$ such that $\chi_{A+bK}$ has positive coefficients. This is equivalent to existence of $\tilde{K}=KT$ such that $\chi_{\tilde{A}+\tilde{b}\tilde{K}}$ has positive coefficients (since $\tilde{A}+\tilde{b}\tilde{K}=T^{-1}(A+bK)T$ and thus the characteristic polynomials of $\tilde{A}+\tilde{b}\tilde{K}$ and $A+bK$ coincide). Observe that $\chi_{\tilde{A}+\tilde{b}\tilde{K}}(\lambda)=(\lambda^{s}+(a_{s}-\tilde{k}_{s})\lambda^{s-1}+\ldots+(a_{2}-\tilde{k}_{2})\lambda+(a_{1}-\tilde{k}_{1}))\chi_{A_{22}}(\lambda)$. Positivity of its coefficients implies that $\tilde{k}_{i}<a_{i}$ for $i=1,\ldots,s$. Since $\tilde{k}_{i}=Kv_{i}$, the last inequality is equivalent to $Kv_{i}<a_{i}$. Thus there must exist $K$ that satisfies (11), so the set of these inequalities is consistent.

Sufficiency. Let $K$ satisfy (11). As in the proof of Necessity, the condition $k_{j}\geq\alpha_{j}$, $j=1,\ldots,n$, means positivity of the closed-loop system $x^{\Delta}=(A+bK)x$. To show that this system is positively uniformly exponentially stable it is enough to verify that the coefficients of $\chi_{A+bK}$ are positive. As in the proof of Necessity, this is equivalent to positivity of coefficients of $\chi_{\tilde{A}+\tilde{b}\tilde{K}}$, with $\tilde{A}$ and $\tilde{b}$ given by  (10) and $\tilde{K}=KT$. As before, $\chi_{\tilde{A}+\tilde{b}\tilde{K}}(\lambda)=(\lambda^{s}+(a_{s}-\tilde{k}_{s})\lambda^{s-1}+\ldots+(a_{2}-\tilde{k}_{2})\lambda+(a_{1}-\tilde{k}_{1}))\chi_{A_{22}}(\lambda)$ and the condition $Kv_{i}<a_{i}$, $i=1,\ldots,k$, implies $a_{i}-\tilde{k}_{i}>0$ for $i=1,\ldots,s$. Observe that the matrix $A_{22}$ defines the uncontrollable subsystem of the system (8), so from  (9) the eigenvalues of $A_{22}$ must belong $S(\mathbb{T})$. Thus their real parts are negative, so $\chi_{A_{22}}$ has positive coefficients (see the proof of Theorem 4.8). Hence, $\chi_{\tilde{A}+\tilde{b}\tilde{K}}$ has positive coefficients as well. ∎

6 Conclusion

We provided a characterization of uniform exponential stability of a positive linear system on a time scale. Surprisingly, regardless of the time scale, the condition is the same, which confirms usefulness of the theory and of the language of time scales. This characterization was used then to develop criteria for positive feedback stabilizability. They guarantee existence of a linear feedback that stabilizes the system and preserves its positivity. This was done for one-dimensional controls only, so a natural extension of this part would be to characterize positive stabilizability for the multi-control case. The next step could be studying stabilization of nonlinear systems via linearization.

Acknowledgment

This work has been supported by the Bialystok University of Technology grant No. S/WI/1/2016.

Appendix

**Calculus on time scales

**A time scale $\mathbb{T}$ is an arbitrary nonempty closed subset of the set $\mathbb{R}$ of real numbers. In particular $\mathbb{R}$, $h\mathbb{Z}$ for $h>0$ and $q^{\mathbb{N}}:=\{q^{k},k\in\mathbb{N}\}$ for $q>1$ are time scales. We assume that $\mathbb{T}$ is a topological space with the relative topology induced from $\mathbb{R}$. If $t_{0},t_{1}\in\mathbb{T}$, then $[t_{0},t_{1}]_{\mathbb{T}}$ denotes the intersection of the ordinary closed interval with $\mathbb{T}$. Similar notation is used for open, half-open or infinite intervals.

For $t\in\mathbb{T}$ we define the forward jump operator $\sigma_{\mathbb{T}}:\mathbb{T}\rightarrow\mathbb{T}$ by $\sigma_{\mathbb{T}}(t):=\inf\{s\in\mathbb{T}:s>t\}$ if $t\neq\sup\mathbb{T}$ and $\sigma_{\mathbb{T}}(\sup\mathbb{T})=\sup\mathbb{T}$ when $\sup\mathbb{T}$ is finite; the backward jump operator $\rho_{\mathbb{T}}:\mathbb{T}\rightarrow\mathbb{T}$ by $\rho_{\mathbb{T}}(t):=\sup\{s\in\mathbb{T}:s<t\}$ if $t\neq\inf\mathbb{T}$ and $\rho_{\mathbb{T}}(\inf\mathbb{T})=\inf\mathbb{T}$ when $\inf\mathbb{T}$ is finite; the forward graininess function $\mu_{\mathbb{T}}:\mathbb{T}\rightarrow[0,\infty)$ by $\mu_{\mathbb{T}}(t):=\sigma_{\mathbb{T}}(t)-t$; the backward graininess function $\nu_{\mathbb{T}}:\mathbb{T}\rightarrow[0,\infty)$ by $\nu_{\mathbb{T}}(t):=t-\rho_{\mathbb{T}}(t)$.

If $\sigma_{\mathbb{T}}(t)>t$, then $t$ is called right-scattered, while if $\rho_{\mathbb{T}}(t)<t$, it is called left-scattered. If $t<\sup\mathbb{T}$ and $\sigma_{\mathbb{T}}(t)=t$ then $t$ is called right-dense. If $t>\inf\mathbb{T}$ and $\rho_{\mathbb{T}}(t)=t$, then $t$ is left-dense.

The time scale $\mathbb{T}$ is homogeneous, if $\mu_{\mathbb{T}}$ and $\nu_{\mathbb{T}}$ are constant functions. When $\mu_{\mathbb{T}}\equiv 0$ and $\nu_{\mathbb{T}}\equiv 0$, then $\mathbb{T}=\mathbb{R}$ or $\mathbb{T}$ is a closed interval (in particular a half-line). When $\mu_{\mathbb{T}}$ is constant and greater than [math], then $\mathbb{T}=\mu_{\mathbb{T}}\mathbb{Z}$.

Let $\mathbb{T}^{\kappa}:=\{t\in\mathbb{T}:t\;\text{ is nonmaximal or left-dense}\}$. Thus $\mathbb{T}^{\kappa}$ is got from $\mathbb{T}$ by removing its maximal point if this point exists and is left-scattered.

Let $f:\mathbb{T}\rightarrow\mathbb{R}$ and $t\in\mathbb{T}^{\kappa}$. The delta derivative of $f$ at $t$, denoted by $f^{\Delta}(t)$, is the real number with the property that given any $\varepsilon$ there is a neighborhood $U=(t-\delta,t+\delta)_{\mathbb{T}}$ such that

[TABLE]

for all $s\in U$. If $f^{\Delta}(t)$ exists, then we say that $f$ is delta differentiable at $t$. Moreover, we say that $f$ is delta differentiable on $\mathbb{T}^{k}$ provided $f^{\Delta}(t)$ exists for all $t\in\mathbb{T}^{k}$.

Example 6.1**.**

If $\mathbb{T}=\mathbb{R}$, then $f^{\Delta}(t)=f^{\prime}(t)$. If $\mathbb{T}=h\mathbb{Z}$, then $f^{\Delta}(t)=\frac{f(t+h)-f(t)}{h}$. If $\mathbb{T}=q^{\mathbb{N}}$, then $f^{\Delta}(t)=\frac{f(qt)-f(t)}{(q-1)t}$.

A function $f:\mathbb{T}\rightarrow\mathbb{R}$ is called rd-continuous provided it is continuous at right-dense points in $\mathbb{T}$ and its left-sided limits exist (finite) at left-dense points in $\mathbb{T}$. If $f$ is continuous, then it is rd-continuous.

A function $f:\mathbb{T}\rightarrow\mathbb{R}$ is called regressive, if $1+\mu(t)f(t)\neq 0$ for all $t\in\mathbb{T}$.

A function $F:\mathbb{T}\rightarrow\mathbb{R}$ is called an antiderivative of $f:\mathbb{T}\rightarrow\mathbb{R}$ provided $F^{\Delta}(t)=f(t)$ holds for all $t\in\mathbb{T}^{\kappa}$. Let $a,b\in\mathbb{T}$. Then the delta integral of $f$ on the interval $[a,b)_{\mathbb{T}}$ is defined by

[TABLE]

Riemann and Lebesgue delta integrals on time scales have been also defined (see e.g. [25]). It can be shown that every rd-continuous function has an antiderivative and its Riemann and Lebesgue integrals agree with the delta integral defined above.

Example 6.2**.**

If $\mathbb{T}=\mathbb{R}$, then $\int\limits_{a}^{b}f(\tau)\Delta\tau=\int\limits_{a}^{b}f(\tau)d\tau$, where the integral on the right is the usual Riemann integral. If $\mathbb{T}=h\mathbb{Z}$, $h>0$, then $\int\limits_{a}^{b}f(\tau)\Delta\tau=\sum\limits_{t=\frac{a}{h}}^{\frac{b}{h}-1}f(th)h$ for $a<b$.