Characterization of Positive Invariance of Quadratic Convex Sets for Discrete-Time Systems Using Optimization Approaches

: A positively invariant set is an important concept in dynamical systems. The study of positively invariant set conditions for discrete-time systems is one interesting topic in both theoretical studies and practical applications research. Different methods for characterizing the invariance of different types of sets have been established. For example, the ellipsoidal and the Lorenz cone, which are quadratic convex sets, have different properties from a polyhedral set. This paper presents an optimization method and a dual optimization method to characterize the positive invariance of the ellipsoidal and the Lorenz cone. The proposed methods are applicable to both linear and nonlinear discrete-time systems. Using nonlinear programming and an induced norm, the positive invariance condition problems are transformed into optimization problems, and the dual optimization method is also used to give equivalent dual forms. Fewer results on the positive invariance condition of Lorenz cones can be found than for the other type of set; this paper fulﬁlls the results of this problem. In addition, the proposed methods in this paper provide more options for checking the positive invariance of quadratic convex sets from the perspective of optimization and dual optimization. The effectiveness of this method is demonstrated by numerical examples.


Introduction
A positively invariant set is an important concept in dynamical systems.As long as the initial state and the subsequent trajectory of the system are always in a particular positively invariant set, the state quantity of the system can be guaranteed to remain in the positively invariant set.Due to their good properties, positively invariant sets play an important role in the study of system stability analysis and feedback controller design [1][2][3][4][5].Lyapunov's stability theory provides theoretical support for the study of the stability of dynamical systems.The characterization of the positively invariant set is one of the main topics in the study of positively invariant sets.Bitsoris first proposed in [6] the sufficient and necessary conditions for the polyhedral set to be the positively invariant set of a linear discrete-time system.Ellipsoids and Lorenz cones are general convex sets, both of which have quadratic forms, and the Lorenz cone itself has a constraint, which makes the study of its positive invariance more difficult than those of polyhedral sets.The Riccati equation approach is proposed in [7] for the positive invariance of an ellipsoidal set of a linear discrete-time system.The maximum invariant ellipsoids for discrete-time systems are also obtained in [8,9] using linear matrix inequalities and bilinear matrix inequalities.A special class of Lorenz cones is constructed using Dickin ellipsoids and some hyperplanes in [10], and the invariant cones of a given system are studied; the structure of the constructed cones is studied by the eigenvalue structure of the matrix in the ellipsoidal expression.The ellipsoidal positively invariant set of the Lorenz system has been estimated in [11] for all positive values of the parameters of the Lorenz system, and the minimum volume value of the ellipsoid has been obtained.There are many ways to determine whether a set is a positively invariant set of the dynamical system.The Lyapunov function is a classical method for studying the stability of systems, and a new method for solving the positively invariant set of Lorenz chaotic systems was obtained in [12] by constructing the Lyapunov function.The problem of the finite-time stability of the closed invariant sets of a class of nonlinear systems is discussed in [13] using Lyapunov functions.Ref. [14] is an excellent review paper about the conditions for positive invariance, or, from an algebraic perspective, the sufficient and necessary conditions for ellipsoidal sets to be the positively invariant sets of linear systems.Ref. [15] investigated sufficient conditions for the existence of robust positively invariant sets for switching systems with an average dwell time based on a novel sequence-based technique.In [16], a sufficient and necessary condition for an ordered class of sets to be a positively invariant set for nonlinear discrete-time systems is proposed from the point of view of the existence of monotone mappings.The condition for ellipsoidal invariance is proposed in [17,18] by the method of linear matrix inequality (LMI), and the optimal solution to the LMI problem is combined to determine whether the maximum ellipsoid is obtained.The S-procedure is an effective method for studying the positive invariance conditions for quadratic convex sets, and the sufficient and necessary conditions for ellipsoidal sets and Lorenz cones to be the positively invariant sets of linear systems are given in [19] based on the Lyapunov stability and S-procedure, but the invariance conditions for Lorenz cones are more complicated.The nonlinear programming techniques in [20,21] provide an alternative way to study positive invariance.The problem of estimating the maximum robust invariant set for discrete-time nonlinear regenerative systems in an optimal control framework is considered in [22].The study of polyhedral sets' positive invariances via optimization is investigated in [23,24].In [25], Nagumo's theorem is used to study the sufficient and necessary positive invariance conditions of the ellipsoid and Lorenz cone for continuous-time systems using optimization techniques.The existence of the solution to the optimization problem is discussed by the Karush-Kuhn-Tucker (KKT) condition.
There are many results on the positive invariance set of a dynamical system, though most of the studies are about polyhedral sets, and there are fewer discussions on the positive invariance of Lorenz cones, which contain a constraint of their own, thereby making the problem more difficult.The main contribution of this paper is the presentation of the positive invariance condition of ellipsoids and Lorenz cones by virtue of the optimization approach and dual form.For convex sets of quadratic forms, such as Lorenz cones, three positive invariance conditions are proposed in this paper using optimization methods.The proposed Lagrange and Wolfe dual optimization methods in [26][27][28] can sometimes simplify the primal problem.The method proposed in this paper establishes a connection among positive invariance conditions, optimization, dual optimization, and the induced norm, which provides an optional method for the positive invariance characterization of quadratic convex sets for nonlinear and linear discrete-time dynamical systems.
The optimization method and dual form have been successfully used in our previous work [24].Additionally, in [24], three sufficient and necessary positive invariance conditions are presented for polyhedral sets and polyhedral cones of discrete-time systems.In this paper, the optimization approach and dual form are also used, however, they are used for ellipsoidal sets combined with the induced norm parametrization.The positive invariance condition for the Lorenz cone is also studied.To check the positive invariance, one only needs to check the sign of the optimization problem's maximum value, which makes it simpler.Moreover, in numerical examples, this paper avoids the complexity of checking the positive invariance of the max-min, as in [24].
The rest of this paper is organized as follows: Section 2 provides some preliminary knowledge and definitions.In Section 3, the sufficient and necessary conditions for the positive invariances of ellipsoids in nonlinear and linear discrete-time systems are studied, respectively.In Section 4, the conditions of the positive invariance of the Lorenz cone are studied.An illustrative example is given in Section 5.The conclusions of this paper are summarized in Section 6.
Notations: x k , x k+1 ∈ R n denotes the state vector, where k ∈ N ∪ {0}.The set of real numbers is given by R. R n denotes a column vector in dimension n × 1. R n×n denotes a real square matrix of an n × n dimension.Q 0(≺ 0) means Q is a positive(negative) definite matrix.

Discrete-Time Dynamical Systems
Discrete-time systems are the systems mainly considered in this paper, and the forms of linear and nonlinear discrete-time dynamical systems are given in the following form: (1) where, A is an n by n dimensional matrix.x k , x k+1 ∈ R n denotes the state vector, where and the ellipsoid and the Lorenz cone are nonempty.
Definition 1 (Positively invariant set).The set D is a positively invariant set of discrete-time systems if and only if x k ∈ D implies x k+1 ∈ D for all k ∈ N ∪ {0}.

Convex Sets
Convex sets in a quadratic form are the main study object of this paper, namely, ellipsoidal sets and Lorenz cones.The ellipsoidal set is defined as where Q ∈ R n×n , and Q 0. The ellipsoidal set is also often written as a unit ball in the form of a quadratic norm, namely, Let Any ellipsoid whose center is not at the origin can be transformed into an ellipsoid whose center is at the origin.Therefore, this paper discusses the invariance condition for the ellipsoidal set whose center is at the origin.The Lorenz cone is defined as where P ∈ R n×n is a symmetric nonsingular matrix with only one negative eigenvalue, λ n .U n is the eigenvector corresponding to the negative eigenvalue λ n .

Lagrange Function
The optimization problem is written in the following form: where the functions f : R n → R, g i (x) : R n → R, i = 1, 2, . . ., m, h j (x) : R n → R, j = 1, 2, . . ., n are continuously differentiable.Additionally, g i (x) denotes the inequality constraints, and h j (x) denotes all of the equality constraints for all i, j.Define the Lagrange function as follows: where λ i , µ j are called the Lagrange operator, and it is required that λ i ≥ 0, i = 1, 2, . . ., m.

Wolfe Dual Theory
Let f (x) be a convex differentiable function with respect to x ∈ R n , and g i (x) is a differentiable convex function.Denote the gradient of the function g(x) by ∇g(x), i.e., ∇g( ∂x n ) T .The primal problem is in the following form: The Wolfe dual form is:

Slater Condition
For the convex programming problem min{ f (x) | g(x) ≤ 0, h(x) = 0, x ∈ R n }, if there exists a feasible point x such that g( x) < 0, the programming problem is said to satisfy the Slater constraint qualification, which is also known as the Slater condition.Hypothesis 2. Suppose the Slater condition can be satisfied for the optimization problem in this paper, which means the inner point of the constraints is nonempty.

Invariance Conditions for Ellipsoids
The necessary and sufficient conditions for checking the positive invariance of ellipsoids sets for linear and nonlinear discrete-time systems are investigated in this section.First, by virtue of the Lagrange duality and the Wolfe duality, three equivalent optimization models for the invariance conditions of nonlinear discrete-time systems are given.Then, the positive invariance condition of the ellipsoid is studied when f d (x k ) = Ax k .

Formulation of Positive Invariance Conditions
The invariance conditions for discrete-time systems, including linear and nonlinear systems, are discussed, respectively, in this section.Theorem 1.For the ellipsoidal set (3) and the nonlinear discrete-time system (2), the sufficient and necessary condition for the ellipsoidal set to be a positively invariant set of the nonlinear discrete-time system is that the optimal value of the following optimization problem is non-negative.
Proof.The set S is a positively invariant set of the system (2) if and only if x k ∈ S, and where x k , x k+1 ∈ R n denotes the state vector, and k ∈ N ∪ {0}.Equation ( 12) holds, which implies that (13) also holds.It is necessary that all function values of ( 13) are greater than or equal to zero, i.e., the minimum value of ( 13) is satisfied by being non-negative.Then, an optimization is formulated as When the system is a linear discrete-time system, i.e., with f d (x) = Ax, the invariance condition is presented by Theorem 2.
Theorem 2. The ellipsoidal set S is a positively invariant set of the linear discrete-time system (1) if and only if the optimal value of the following optimization problem is non-negative.
Since the proof procedure of Theorem 2 is similar to that of Theorem 1, it is not repeated here.
For the ellipsoids in (5), the sufficient and necessary conditions using the induced norm to check the positive invariance of the linear discrete-time systems are given below.Mx .Then, M is said to be the induced norm of the vector norm x .Theorem 3. The ellipsoidal set S(Q, 1) is a positively invariant set of the linear discrete-time system x k+1 = Ax k if and only if the optimal value of the optimization problem below is positive.
Proof.From the definition of the induced norm in Lemma 1, a sufficient and necessary condition for the ellipsoidal set S(Q, 1) to be a positively invariant set of a linear discretetime system is that (17) holds for all x.
That is, Transform ( 18) into an optimization problem with the following constraints: From Theorems 1 and 2, the positive invariance problem is formulated as an optimization problem.In practice, the dual approach sometimes simplifies the solving of the primal problem.Next, the Lagrange dual and the Wolfe dual forms of ( 11), (15), and ( 16) are discussed, respectively.

Lagrange Dual
The Lagrange dual is a convex optimization problem regardless of whether the primal problem is convex or not, and the dual gives at least a lower bound on the optimal value of the primal problem.In this section, the Lagrange dual form of the primal problem is considered.Theorem 4 gives the positive invariance conditions in the Lagrange dual form of the ellipsoidal set for the nonlinear dynamical systems.Theorem 4. Consider the nonlinear discrete-time system x k+1 = f d (x k ) and the ellipsoidal set be S be a continuous differentiable function with respect to x.The ellipsoidal set S is a positively invariant set of the nonlinear discrete-time system (2) if and only if there exists λ ≥ 0 such that the optimal value of the following optimization problem is non-negative.
Proof.Let (11) be the primal problem.Next, we derive its dual form M(x). Let the Lagrange multiplier be λ and λ ≥ 0.Then, the Lagrange function is The, min x∈R n max λ≥0 L(x, λ) is equivalent to the primal problem that satisfies the constraint, i.e., the Lagrange dual of the primal problem (11) is in te reason that the Lagrange dual satisfies As a result, the optimal value of the primal problem must be non-negative when the optimal value of the dual problem ( 20) is non-negative.
Remark 1. Lagrange dual forms are valid only when the optimal value function of the inner optimization problem can be reduced to an analytic formula.Therefore, in Theorem 4, it is necessary to satisfy the requirement that Similarly, sufficient conditions are given for the ellipsoidal sets to be positively invariant sets of linear discrete-time systems.
Theorem 5. Let the linear discrete-time system be x k+1 = Ax k and the ellipsoidal set be S = {x ∈ R n | x T Qx ≤ 1}, where Q ∈ R n×n , and Q 0. Then, the ellipsoidal set is a positively invariant set of the linear discrete-time system if and only if there exists λ ∈ [0, 1] such that the optimal value of the following optimization problem is non-negative.
Proof.The proof of the Lagrange dual form is similar to Theorem 3, and, here, the focus is on the range of values of the Lagrange multipliers.
Let the Lagrange function of the primal problem be Take x = 0, the positive invariance condition should also satisfy L(0, λ) ≥ 0, i.e., Theorem 6 ([19]).The ellipsoidal set S is a positively invariant set of the linear discrete-time system x k+1 = Ax k if and only if there exists λ ∈ [0, 1] such that A T QA − λQ 0.
Proof.From the proof of Theorem 5, it follows that if the ellipsoidal set S is a positively invariant set of the linear discrete-time system (1), then there exists λ ∈ [0, 1] such that the optimal value of the optimization problem ( 25) is non-negative.Then The proof is completed.
Remark 2. Theorem 5 was once proved in [19] using an S-procedure, and, in this paper, a novel proof is given from an optimization point of view that shows a direct connection between positively invariant sets and optimization.
Theorem 7. The linear discrete-time system is x k+1 = Ax k , and, when the ellipsoid is given by a quadratic norm of the form (4), then the ellipsoid is a positively invariant set of the linear discrete-time system (1) if and only if there exists λ ≤ 0 such that the optimal value of the following optimization problem is positive.
Remark 3. The proof of Theorem 7 is similar to the proof in Theorem 4; the proof is omitted here.
By comparing the primal optimization problem with the Lagrange dual optimization problem, it can be seen that the feasible domain of the dual optimization problem is simpler than that of the primal problem.Subsequently, the case in which the primal problem is a convex optimization problem is discussed, and the Wolfe dual method's application to the primal optimization problem is discussed as well.

Wolfe Dual Forms
When the objective function is convex and the constraint is also convex, the primal optimization problem can be solved through the Wolfe dual form.The Wolfe dual form of the ellipsoid's positive invariance condition for the discrete-time systems is studied in this section.

Remark 4.
Since the primal problem is a convex optimization problem, and the Slater condition is satisfied, the strong duality theorem holds.Therefore, the optimal value of the dual problem is equal to that of the primal problem.Then, whether the ellipsoidal set is a positively invariant set of the discrete-time system can be determined by the sign of the optimal value of the optimization problem.Theorem 8.The ellipsoidal set is expressed as S = {x ∈ R n | x T Qx ≤ 1}, where Q ∈ R n×n , and Q 0. The nonlinear discrete-time system is given by (2).Let 1 − f T d (x)Q f d (x) be a convex function differentiable with respect to x, and let x T Qx − 1 be a convex function.Then, the ellipsoidal set S is a positively invariant set of nonlinear discrete-time systems (2) if and only if there exists λ ≥ 0 such that the optimal value of the following problem is non-negative.
Proof.From the primal problem (11), both the objective function and the constraint are convex functions.Then, the Wolfe dual theorem can be applied, i.e., max Theorem 9. Let the linear discrete-time system be (1) and the ellipsoidal set be S = {x ∈ R n | x T Qx ≤ 1}, where Q ∈ R n×n , and Q 0. Assume that 1 − x T A T QAx is a convex function and that x T Qx − 1 is a convex function; then, the ellipsoidal set S is a positively invariant set of linear discrete-time systems (1) if and only if there exists λ ∈ [0, 1] such that the optimal value of the following problem is non-negative.
Proof.From the assumptions, it is clear that ( 15) is a convex optimization problem, so the Wolfe dual theory can be applied.In particular, when x = 0, there is 1 − λ ≥ 0, i.e., 0 ≤ λ ≤ 1.The Wolfe dual of ( 15) is Remark 5.By comparing Theorems 4 and 5 with Theorems 8 and 9, it can be seen that the Wolfe dual applies to the case in which the primal problem is a convex optimization problem, while the Lagrange dual applies to the more general case.

Positive Invariance Conditions for Lorenz Cone
The Lorenz cone is a convex set in quadratic form, but the positive invariance is more complicated because the Lorenz cone itself contains the constraint x T PU n ≤ 0. The positive invariance condition for the Lorenz cone proposed in this paper makes it simpler.Similar to the ellipsoid, three equivalent sufficient conditions for the positive invariance of the Lorenz cone for nonlinear discrete-time systems are given.
Theorem 10.Lorenz cone (6) is a positively invariant set for a nonlinear discrete-time system (2) if and only if the optimal value of the following optimization problem is non-negative.
When f d (x) = Ax, and the Lorenz cone S L is a positively invariant set of linear discrete-time systems, then the optimal value of the following optimization problem needs to be non-negative.
Proof.First, the Lorenz cone is a positively invariant set of the nonlinear discrete-time system (2) if and only if x k ∈ S L , and x k+1 ∈ S L ; i.e., it needs to satisfy When formulate as an optimization problem, it is expressed, i.e., as follows: If the discrete-time system is linear, i.e., when f (x k ) = Ax, the positive invariance condition is min Next, the Lagrange dual optimization form for the optimization problem in (32) and ( 33) is studied.
Theorem 11.Consider that the nonlinear discrete-time systems are x k+1 = f d (x k ), and the Lorenz cone is given by (6).Let − f T d (x)P f d (x) be a continuous differentiable function with respect to x.Then the Lorenz cone S L is a positively invariant set of the nonlinear discrete-time system (2) if and only if there exists λ, µ, η ≥ 0 such that the optimal value of the following optimization problem is non-negative.
If f d (x) = Ax, and the Lorenz cone S L is a positively invariant set of linear discrete-time systems, then the optimal value of the following optimization problem needs to be non-negative.
Proof.Taking (32) as the primary problem and introducing the multiplier λ, µ, η ≥ 0, one can write its Lagrange function as follows: let the primary problem be P(x).There exists max λ,µ,η≥0 L(x, λ, µ, η) is equivalent to (32), and the Lagrange dual of (32) is max λ,µ,η≥0 Since the optimal value of the primal problem must be greater than or equal to the optimal value of its Lagrange dual problem, when the optimal value of the Lagrange dual problem is non-negative, the optimal value of the primal problem must be non-negative.When f (x k ) = Ax, the Lagrange dual problem of (33) is max λ,µ,η≥0 min When the optimization problem in (32) and ( 33) is convex optimization, it can be solved by the Wolfe dual form optimization.
Theorem 12.The Lorenz cone S L = {x ∈ R n | x T Px ≤ 0, x T PU n ≤ 0}, where P ∈ R n×n .Let (32) be a convex optimization problem.Then, the Lorenz cone S L is a positively invariant set of nonlinear discrete-time systems (2) if and only if there exists λ, µ, η ≥ 0 such that the optimal value of the following problem is non-negative.
If f d (x) = Ax, when (33) is convex optimization, the Lorenz cone S L is a positively invariant set of linear discrete-time systems if and only if there exists λ, µ, η ≥ 0 such that the optimal value of the following optimization problem needs to be non-negative.
The proof of Theorem 12 is similar to the proof of Theorem 8; it is omitted here.Remark 6.Our above conclusions can also be generalized to additional systems, such as Markovian jump systems [29].

Numerical Examples
In this section, some examples are given to illustrate the theorems given in Sections 3 and 4. The numerical examples in this paper are mainly solved by Matlab and Lingo.Note that the change in the initial state does not affect the final determination of whether the set is a positively invariant set or not.

Remark 7.
In all of the examples in this paper, the superscript k in x (k) represents the time index corresponding to the state of discrete-time systems (1) and (2).
In method 1, the positive invariance of this ellipsoidal set is checked by using (11) of Theorem 1.We transform the invariance problem into an optimization problem, i.e., Firstly, simplify the objective function.
Next, the constraint function is simplified. x The simplified optimization problem is 2 ) 2 ≤ 1.
Solve nonlinear programming problems in Matlab by applying the f mincon function.A 2 × 1-dimensional matrix is generated as the initial state by using rand (2,1).The optimal value of this optimization problem is obtained by MATLAB as 0.5 > 0.Then, the ellipsoidal set S is a positively invariant set of the nonlinear discrete-time system.
In method 2, apply (20) in Theorem 4 to check whether the ellipsoidal set is a positively invariant set of this nonlinear discrete-time system.The optimization problem in method 1 is transformed into its Lagrange dual optimization form, i.e., Firstly, let the inner optimization problem be g(x 2 ), and take the partial derivatives of x 2 ) to be equal to zero, respectively, i.e., ∂g ∂x .

By substituting x
(k) 2 into the optimization function, one can obtain the function that is only related to λ.When λ = 0.5, the optimal value of the function is 1 > 0. Therefore, the ellipsoidal set is the positively invariant set of the nonlinear discrete-time system.
In method 3, since the optimization problem in method 1 is a convex optimization problem, (29) in Theorem 8 can be applied to solve this problem.That is, it is necessary to find a u that satisfies u > 0 such that the optimal value of the objective function is non-negative.
In this example, let λ = 0.25, and substitute the objective function to obtain the value 0.25 > 0. Therefore, the ellipsoidal set is the positively invariant set of this nonlinear discrete-time system.
Remark 8. Due to the complexity of nonlinear systems, it is more complicated to theoretically derive the conditions in which quadratic convex sets are the positively invariant sets of nonlinear systems.
The effectiveness and advantage of the optimization method proposed in this paper for determining the quadratic convex set to be a positively invariant set of a nonlinear system is obviously simple.
Example 2. Let the linear discrete-time system x k+1 = Ax k , where A = 0 −3.2 −0.1 0.4 and the set of ellipsoids is S = {(x 15) in Theorem 2 to check it.Substituting the data in this example into (15) yields an optimization problem.
Firstly, the objective function is simplified.
We can solve nonlinear optimization problems using Matlab with the f mincon function.A 2 × 1-dimensional random matrix is generated as the initial state by using rand(2, 1) and the optimal value of the optimization function as −9.4002 < 0.Then, the ellipsoidal set is not a positively invariant set of this linear discrete-time system.Remark 9.In the positive invariance condition given in [19], for the ellipsoid to be a linear discretetime system, it can be verified that the ellipsoid is a positively invariant set of the discrete-time system, but it is more complex in this example.
Example 3. The linear discrete-time system is x k+1 = Ax k , where A = 0.5 0 0 0.5 .Additionally, the ellipsoidal set is (x 2 ) 2 ≤ 1.In method 1, using (15) in Theorem 2, the problem of checking that the ellipsoidal set is a positively invariant set of a linear discrete-time system is transformed into an optimization problem.
Then, the matrices A, Q are substituted into the objective function to obtain Then, the optimization function is Solve nonlinear optimization problems using Matlab with the f mincon function.The initial state is set as rand(2, 1) in MATLAB, and the optimal value of the objective function is 0.75 > 0. Therefore, the ellipsoidal set S is the positively invariant set of the linear discrete-time system.In method 2, the ellipsoidal set is written in the form of (4), and the invariance condition of the ellipsoidal set is transformed into an optimization problem by using (16) in Theorem 3.
Then, the matrix A and Q are substituted and simplified to obtain the following optimization problem, i.e., Solve nonlinear optimization problems using Matlab with the f mincon function.The initial state is set as rand(2, 1) in MATLAB, and the optimal value of the objective function is 0.75 > 0. Therefore, the ellipsoidal set S is the positively invariant set of the linear discrete-time system.In method 3, use (25) in Theorem 5 to determine whether the ellipsoidal set is a positively invariant set of the linear discrete-time system. 2 ) 2 + λ((x Let the inner optimization function be g(x 2 ), and take the partial derivatives with respect to x  Take x 1 = x 2 = 0, the optimal value is 1 > 0 by substituting it into the objective function.When x 1 , x 2 = 0, λ = 0.25 can be obtained, and the optimal value is 0.75 > 0 by substituting it into the objective function.Therefore, the set of ellipsoids is the positively invariant set of this linear discrete-time system.
In method 4, apply (28) in Theorem 7 to determine the positive invariance of the ellipsoidal set.That is, in this example, it is necessary to find a λ < 0 such that the optimal value of the optimization problem is positive.Substituting the data in this example into (28) and simplifying gives 2 ) 2 + λ((x The optimal value is 0.75 when λ = 0.75.Therefore, this ellipsoidal set is the positively invariant set of this linear discrete-time system. In method 5, apply (29) in Theorem 8 to determine the positive invariance of the ellipsoidal set.This method needs to check whether a λ ≥ 0 can be found such that the optimal value of the optimization problem is non-negative.The constraint condition is satisfied when λ = 0.25, and the value of 0.75 > 0 can be obtained by substituting λ into the objective function.Therefore, the ellipsoidal set is the positively invariant set of the linear discrete-time system.
Examples 4 and 5 apply Lingo to solve the optimization problems.2 ) 2 ≤ 0, x 2 ≥ 0}, and the nonlinear discrete-time system is x In method 1, (32) in Theorem 10 is applied to determine whether the Lorenz cone is a positively invariant set of this nonlinear discrete-time system.
According to (32) in Theorem 10, when the optimal value is great than/equal 0, the Lorenz cone is the positively invariant set of the nonlinear discrete-time system.

Lemma 1 .
A vector norm • on C n is known, and, for any square matrix M ∈ C n×n , let M = sup x =1

Example 4 .
The Lorenz cone is represented by S L = {(x min x∈R n − f d (x) T P f d (x) s.t f d (x) T PU n ≤ 0 x T Px ≤ 0 x T PU n ≤ 0,