1. Introduction
In nature and engineering applications, there exist some quantities that potentially remain to be non-negative, such as the populations of organisms within ecological systems [
1], the propagation rate of a signal in network communication [
2], and the concentration of liquid in chemistry [
3]. Such a system characterized by non-negative quantities is commonly known to be positive for which the state variables and output trajectories reside in the first quadrant whenever the inputs and initial states are non-negative [
4]. The studies concerning positive systems are derived from [
5], based on which significant results such as reachability and controllability [
6], Kalman-Yakubovich-Popov [
7] have been developed. Due to the positivity, the state variables remain in the non-negative cone instead of the entire linear space [
8]. This property leads to the obstacle of applying the well-established general theories to positive systems, which promotes specialized research including state-feedback control [
9], output–feedback control [
10,
11], observer design [
12,
13], etc.
On the other hand, actual systems bear complex dynamics and the order of exact models is relatively high, which brings about great difficulties in the control and implementation of related systems. Naturally, the model reduction problem aiming at approximating a higher-order model by a lower-order one with a sufficiently small error has received considerable attention [
14,
15,
16,
17,
18], and various model-reduction techniques have been established. The traditional balanced-truncation method was extended for the inhomogeneous initial condition case [
19]. A lower-order stable transfer function was explored to approximate a given transfer function with the
norm [
20]. A complete characterization of all optimal Hankel-norm approximations to a rational transfer function was provided in [
21]. Nevertheless, two challenges arise when attempting to adapt the developed model reduction approaches to positive systems. On one hand, when approximating a specified positive system, we expect that the reduced-order model is preserved to be positive, which imposes additional constraints on system matrices. On the other hand, the approximation error is usually characterized by specific criteria of an associated error system, which can be transformed into a type of bilinear matrix inequalities (BMIs) involving the coupled terms between the Lyapunov variable and unknown system matrices [
22]. An intuitive idea to resolve this issue is to transform the bilinear constraints into linear ones. In addition, the congruent transformation method [
23] and the D-K iteration method [
24] have been proposed to find feasible solutions in some cases. Nevertheless, the above results did not take the intrinsic positivity constraints into account, and solving the non-convex stability condition and the positivity constraint simultaneously remains a major challenge.
This paper explores the positivity-preserving model reduction method for positive systems. A reduced-order positive model is sought to capture the dominant dynamics of the original high-order system within a prescribed error bound. The sufficient and necessary condition with the form of BMI for the positivity-preserving model reduction is developed in a way such that the approximation error is minimized within a guaranteed level. With the aid of the inner-approximation strategy, we develop a successive convex optimization (SCO) algorithm, under which the desired reduced-order model is achieved by sequentially solving the approximated convex problems. Overall, the main contributions of the paper are stated as follows:
This paper proposes a positivity-preserving model reduction scheme, which guarantees the performance of the resulted error system and preserves the positivity simultaneously.
To render the positivity-preserving model reduction problem numerically tractable, we propose an inner-approximation strategy, and then establish an SCO algorithm to solve a type of BMI problem without parametrization techniques.
To achieve a smaller approximation error, the zero initial condition can be adapted to iterate the reduced-order model, which simplifies the design process by abolishing the initialization step.
The paper is arranged as follows:
Section 2 discusses the problem statement and presents the necessary fundamental results.
Section 3 establishes the conditions required for the existence of a positive reduced-order model, based on which a significant algorithm for optimizing the desired reduced-order model is provided. In
Section 4, we provide an example to illustrate the effectiveness and potential benefits of the presented results. Finally,
Section 5 summarizes the conclusions and highlights future directions for improvement.
Notations: are introduced to denote the set of real numbers, n-dimensional column vectors and -dimensional matrices, respectively, and defines n-dimensional column vectors whose all elements are non-negative. and represent the identical and zero matrices with appropriate dimensions, respectively. For a square matrix X, indicates that X is positive-definite (negative-definite), and , where is the transpose of X. For an arbitrary matrix , implies that every elements in Y is non-negative (positive).
2. Preliminaries and Problem Formulation
Consider a class of stable systems, whose dynamics can be described as follows:
where
are the vectors of state, controlled input and output, respectively, and
are known constant matrices. In this paper, we assume that the system (
1) is positive, and the definition of positive systems is given as follows.
Lemma 1 ([
25])
. Given the initial condition and the input , the system (1) is said to be positive if . In the following, we present a well-established condition for checking whether the system (
1) is positive.
Lemma 2 ([
26])
. The system (1) is positive if and only if . In this paper, we aim at constructing a lower-order model to approximate the system (
1), and the dynamics of the reduced-order model can be formulated in the following form
where
are the states and approximated outputs of the reduced-order system (
2), and
are constant matrices to be designed.
Denote the augmented vector
and the approximation error
. Substituting (
2) into (
1) results in the following error system
where
For the error system (
3), the transfer function from
to
is given by
The primary objective of the conventional
model reduction problem is to search for a reduced-order model (
2) such that
where
is a specified scalar. Nevertheless, a positive model is more likely to capture the dynamics of (
1) with a sufficiently small error. That is, apart from the
performance (5), it is expected that the reduced-order model is preserved as positive.To ensure positivity, it follows from Lemma 2 that
.
For convenience, we denote . The problem formulation for the positivity-preserving model reduction can be stated as follows:
Positivity-preserving model reduction: Given a stable positive system (
1), construct a positive model with a lower order such that
To characterize the performance in (5), we present the following bounded real lemma.
Lemma 3 ([
27])
. Given the state-space realization (3) and a scalar , the following statements are equivalent:- 1.
The error system in (3) is asymptotically stable with an performance . - 2.
There exists matrix satisfying
To facilitate the construction of the reduced-order model, an efficient condition equivalent to (6) is provided.
Lemma 4. The error system in (3) is asymptotically stable with an performance (5) if and only if there exist matrices and G such that Proof. (Necessity) Letting , the condition (6) is expressed as (7) exactly.
(Sufficiency) Assume that the condition (7) holds, and we can infer from (7) that
, which means that
G is non-singular. Further,
yields
, and it is straightforward that (7) implies
Pre- and post-multiplying (8) with and its transpose leads to (6). The proof is completed. □
Remark 1. Lemma 3 establishes the necessary and sufficient condition ensuring that the error system (3) is stable under a given error bound. Note that (6) is hard to solve, for which there exist coupled terms between the Lyapunov variable P and , and more significantly, both P and bear structural constraints. To address this issue, an equivalent condition (7) is developed. It is observed that (7) is easier to implement than (6) by eliminating the coupled terms between P and , and thus facilitating parametrization on the positivity constraint. 3. Main Results
This section designs the desired reduced-order model to approximate the original high-order one. Firstly, we realize the parameterization of unknown matrices in (
2) by providing an alternative characterization for the system (
3). We then formulate the design conditions for a reduced-order model in terms of BMIs, guaranteeing the stability and
specification of the error system. Finally, an iterative procedure based on the inner-approximation strategy is proposed to achieve the desired reduced-order model.
Denote the augmented matrices:
and the error system (
3) can be described by
with
.
Remark 2. To render the problem of positivity-preserving model reduction tractable, the matrices to be designed are reorganized as a single matrix . Notice that is embedded between two known matrices, and by Lemma (4), the design condition for is recast into a BMI involving the coupled term between G and . Furthermore, implies , and the additional positivity constraint on sharpens the difficulty in searching for feasible solutions for the non-convex optimization problem.
Theorem 1. For a given scalar and known matrices , the problem of positivity-preserving model reduction is solvable if there exist matrices such that andwhere Proof of Theorem 1. From the above analysis, (7) is presented in the form of BMI, which can be further reformulated as
where
For given feasible matrices
, (10) is equivalent to
which contains the linear term
and the bilinear term
. By performing the matrix decomposition [
22], we have
where
U is introduced as an auxiliary variable satisfying
. Denote
. Replacing the bilinear term
with its upper-approximation and applying the Schur complement lemma yields
and we can reorganize (13) as
Recalling that
, then (14) can be conservatively expressed as
By the Schur complement lemma, one obtains
where
For a given , we obtain according to the concavity of . Replacing with its over-approximation and pre- and post-multiplying (16) with and its transpose, the convex condition (9) is established.
□
With the aid of the inner-approximation strategy (12), Theorem 1 establishes a convex approximation for the non-convex constraint (9). Note that given known matrices
satisfying
, the condition (9) is formulated in the form of linear matrix inequality (LMI), and it suffices to check the feasibility of the original non-convex problem by alternatively solving the convex surrogate optimization problem
In addition, the superscript
can be viewed as an iteration indicator, and once given the feasible initial condition
, one can obtain a sequence of reduced-order models as the iteration proceeds. Next, we summarize the above procedures into the following SCO algorithm (Algorithm 1).
Algorithm 1 Successive convex optimization algorithm for the computation of reduced-order models |
Require: : given tolerance level; : the maximum number of iterations.
- 1:
Given the initial condition , determine the initial condition by solving
and fix - 2:
while do - 3:
Solve (18) - 4:
if then - 5:
- 6:
break; - 7:
end if - 8:
- 9:
end while - 10:
Result : the reduced-order model; : the performance index.
|
Remark 3. Algorithm 1 provides a straightforward solution to the non-convex BP problem without any parametrization techniques. In addition, it is worth noting that Algorithm 1 is performed iteratively, which presents a potential computational burden relative to the conventional method. As mentioned above, an essential prerequisite for implementing Algorithm 1 is that the initial condition should satisfy . Nevertheless, finding a feasible initial solution can be as equally challenging as solving the original non-convex problem. Fortunately, the original high-order system is stable, and the reduced-order model with enables the resulted error system to be stable. In this setting, compared with [24], the proposed algorithm is easier to implement by avoiding the arduous initialization step. 4. Simulation
Compartmental networks consist of a finite number of homogeneous, well-mixed subsystems, called compartments, which exchange with each other and the environment [
28]. In this example, we compare the proposed model-reduction method with [
26]. Consider the compartmental model of six states with the following system matrices:
and we observe that the system is stable. The aim of this example is to construct reduced-order positive models with different orders to approximate the original system.
Case 1:
According to Remark 3, the initial condition
is adopted, and by Algorithm 1, one can iteratively obtain the reduced-order model as:
Moreover, by adopting the model reduction method in [
26], we obain
Similarly, the corresponding second-order model for Algorithm 1 and [
26] can be readily obtained as follows.
Moreover, the
performance indices
with respect to the above cases are listed in
Table 1. It can be readily observed from the attained results that the provided algorithm can achieve smaller
performance indices than [
26]. To evaluate the approximation effects of different model reduction methods,
Figure 1 and
Figure 2 provide the frequency response gains of the original and reduced-order systems, and
Figure 3 and
Figure 4 depict the frequency response gains of the associated error systems with different orders. We can easily infer that compared with
and
, the frequency gains of proposed forth-order
model and second-order one
are much closer to the original system, and the associated error system with
and
have smaller gains than that for corresponding
and
, which illustrates the benifits of the proposed algorithm.
To further verify the efficiency of the obtained reduced-order models, we will show their approximation performances under the initial conditions
and the input
Figure 5 and
Figure 6 depict the output trajectories of the original model and the reduced-order models with
and
, respectively, and
Figure 7 and
Figure 8 show the corresponding approximation errors between the original model and the reduced-order ones. It can be seen that the proposed method can approximate the original model without causing remarkable errors.