REALIZATION THEORY OF NASH SYSTEMS∗ ˇ ´ † , MIHALY ´ JANA NEMCOV A PETRECZKY‡ , AND JAN H. VAN SCHUPPEN§ Abstract. This paper deals with realization theory of so-called Nash systems, i.e. nonlinear systems the right-hand sides of which are defined by Nash functions. A Nash function is a semialgebraic analytic function. The class of Nash systems is an extension of the class of polynomial and rational systems and it is a subclass of analytic nonlinear systems. Nash systems occur in many applications, including systems biology. Formulation of the realization problem for Nash systems and a partial solution to it are presented. More precisely, necessary and sufficient conditions for realizability of a response map by a Nash system are provided. The concepts of semi-algebraic observability and semi-algebraic reachability are formulated and their relationship with minimality is explained. In addition to their importance for systems theory, the obtained results are expected to contribute to system identification and model reduction of Nash systems. Key words. Nash systems, realization theory, reachability, observability AMS subject classifications. 93B15, 93C10, 93B25, 14P20

Realization theory is one of the central topics of system theory. It serves as a theoretical foundation for model reduction, system identification and filtering/observer design. Its aim is to answer the following questions: (1) Under which conditions is it possible to construct a (preferably minimal) system of a certain class generating the specified input/output behaviour? (2) How to characterize minimal systems of a certain class which generate the specified input/output behaviour? In this paper we investigate realization theory of the class of Nash systems. Nash systems are continuous-time dynamical systems with Nash submanifolds as statespaces and such that the right-hand sides of the differential equations which determine their dynamics, and the output functions are Nash functions. By a Nash function we mean an analytic function satisfying an algebraic equation. A Nash submanifold of Rn is a smooth manifold which is defined by polynomial equalities and inequalities. Nash functions and Nash manifolds are named after John Nash, who used them in his seminal paper on embedding of manifolds in eucledian space [24]. Our approach to realization theory for Nash systems is a continuation of the approach followed by [17] for nonlinear systems, [36, 4] for polynomial systems and the one in [41, 3, 28, 29] for rational systems. It is closely related to the more recent work [6, 5], although we consider a different system class. We introduce the framework of Nash systems and the concepts of semi-algebraic reachability and semi-algebraic observability. We link semi-algebraically reachable, semi-algebraically observable, and minimal Nash realizations. We expect that the results on realization theory for Nash systems will be useful for system identification, model reduction, filtering and control design of Nash systems. ∗ This work was partially supported by the ITEA project Twins 05004, NWO project 613.000.442 ˇ project 13-16764P. and GACR † Institute of Chemical Technology, Prague, Department of Mathematics, Technick´ a 5, 166 28 Prague 6, Czech republic ([email protected]). ‡ Ecole des Mines de Douai, Department of Computer Science and Automatic Control, F-59500 Douai, France, [email protected]. § Centrum Wiskunde & Informatica (CWI), P.O. Box 94079, 1090 GB Amsterdam, The Netherlands ([email protected]).

1

2

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

To the best of our knowledge, both the class of Nash systems and the presented results on realization theory are new. This paper is based on the Ph.D. thesis of one of the authors, see [26]. All the major proofs and results of the paper have appeared in [26]. The results of the paper were also announced in [27]. The main difference between [27] and the current paper is that the latter provides all the detailed proofs. In addition, with respect to [27] the presentation of the current paper is more detailed and certain technical concepts, such as semi-algebraic observability of Nash systems, have been further refined. The idea of studying semi-algebraic systems, i.e. systems described by polynomial equalities and inequalities, is quite well-established in the hybrid systems community, see [23, 22, 31]. In [30] realization theory of discrete-time semi-algebraic hybrid systems is investigated. Since the class of systems in [30] is different from the class of Nash systems considered in this paper, the results of this paper neither imply nor are implied by those of [30]. For semi-algebraic hybrid systems, problems related to realization theory are investigated in [10, 11, 12, 13, 21]. For polynomial and rational systems, realization theory, reachability and observability are treated in [1, 2, 3, 4, 41, 28, 29]. The study of Nash systems is motivated in Section 1. The necessary background on real algebraic geometry is presented in Section 2. In Section 3 the class of Nash systems is introduced. The realization problem for the class of Nash systems is formulated in Section 4. In addition, Section 4 presents the notions of minimality, semi-algebraic reachability and semi-algebraic observability for Nash systems. The results on the existence and the minimality of Nash realizations are presented in Section 5. The proofs are presented in Section 6 and Section 7. The appendix contains the proofs of a number of auxiliary results. 1. Motivation. The motivation for studying Nash systems can be summarized as follows: • Nash systems are both computable and general. • Nash systems are universal approximators. • Nash systems are relevant for systems biology. Let us elaborate on each of these aspects. Generality and computability. The class of Nash systems lies between polynomial/rational systems and analytic systems. While more general than the former, it still allows a constructive description by means of finitely many polynomial equalities and inequalities. Hence, it might still be possible to derive computational methods for control and analysis of Nash systems. Universal approximation property. It is well-known that (piecewise-)polynomial systems can be used to approximate the behavior of nonlinear systems, see [19]. Since polynomial systems are Nash, it follows that Nash systems might be used to approximate general nonlinear systems. Relevance for systems biology. Polynomial and rational systems, and thus Nash systems, are commonly used in systems biology to model metabolic, signaling, and genetic networks. In [32, 33] Savageau proposes an alternative (so-called power-law) framework, also known as Biochemical Systems Theory, for modeling metabolic and gene-regulatory networks by dynamical systems which still belong to the class of Nash systems, but are more general than polynomial and rational systems. In this framework, see also [34, 39], all processes are modeled by products of power-law functions, i.e. sums of products xq11 xq22 · · · xqnn of state variables x1 , . . . , xn taken to rational powers (qi , i = 1, . . . , n are rational numbers). Such functions are a special

NASH REALIZATIONS

3

case of Nash functions. Note that the values of rational exponents (kinetic orders) in power-law systems can be related to parameters of different rate laws such as for example Michaelis-Menten kinetics, see [39]. The role of power-law systems in modeling of signal transduction pathways is discussed in [37]. For further examples and details on application of power-law, and thus Nash, systems in biology see [40, 18]. There are several frameworks for modeling metabolic networks. Some of them (including power-law framework) are compared in [14]. The tendency modeling framework which extends the power-law framework combines mass-action and power-law kinetics into tendency kinetics, [38]. It also leads to Nash systems. Further, considering Nash submanifolds as the state-spaces allows a natural implementation of possible conservation laws and restrictions on the state variables. For example, one can capture the assumption on positivity of state variables which represent the concentrations of reactants of a chemical system by defining the state-space as the positive orthant of Rn which is a Nash submanifold. 2. Real algebraic geometry. In this paper we follow the notation and terminology of [42, 20, 8] on commutative algebra and real algebraic geometry. However, we do not state all results in their full generality (over arbitrary real closed fields), we work only over the field R of real numbers. First, let us recall the notion of transcendence degree. Let A be an integral domain and let R be a sub-ring of A. Recall from [42] the notion of quotient field (field of fractions) QA of A. Then the transcendence degree of A over R, denoted in this paper by trdeg A, is defined as the transcendence degree over R of the field of fractions QA of A. The transcendence degree of QA (and hence of A) equals the greatest number of algebraically independent elements of F over R. In accordance with the terminology of [42], if S, R are integral domains, S is a sub-ring of R, then we say that R is algebraic over S, if the field of fractions QR of R is an algebraic extension of the field of fractions QS of S. This is equivalent to saying that for any element r ∈ R there exists a polynomial 0 6= Q ∈ S[X] such that Q(r) = 0. If trdeg S = d and x1 , . . . , xd is a finite transcendence basis of S, then r ∈ R is algebraic over S if and only if there exists a polynomial in d + 1 variables P ∈ R[X0 , . . . , Xd ] such that P 6= 0, P (X0 , x1 , . . . , xd ) 6= 0 and P (r, x1 , . . . , xd ) = 0. A subset S ⊆ Rn is called semi-algebraic if it is a set of points of Rn which satisfy finitely many polynomial equalities and inequalities, or if it is a finite union of such sets. We say that a semi-algebraic subset S ⊆ Rn is semi-algebraically connected if it cannot be written as a union of two disjoint closed semi-algebraic sets in S. In [8] it is shown that semi-algebraically connected sets are also connected. Let S ⊆ Rn and S 0 ⊆ Rm be semi-algebraic sets. A function f : S → S 0 is a semi-algebraic function if its graph is a semi-algebraic set in Rn+m . In the sequel we will need the following technical lemma on the existence of a partition of a semi-algebraic set. Lemma 2.1 ([8]). Let S ⊆ Rn be a semi-algebraic set and let f : S → R be a semialgebraic function. There exist an integer m > 0 and semi-algebraic subsets S1 S , . . . , Sm m of S and polynomials gi ∈ R[X1 , . . . , Xn , Y ], i = 1, . . . , m such that S = i=1 Si , Si ∩ Sj = ∅ for all i 6= j ∈ {1, . . . , m} and such that for all x ∈ Si , gi (x, Y ) is not identically zero and gi (x, f (x)) = 0. Next, we define the notion of Nash function and Nash manifold. Definition 2.2. A Nash function on an open semi-algebraic set S ⊆ Rn is an analytic and semi-algebraic function from S to R. We denote the ring of Nash

4

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

functions on S by N (S). We say that a function f : S → Rk is a Nash function, if the coordinate functions of f are Nash functions, i.e. if for some f1 , . . . , fk ∈ N (S) it holds that ∀x ∈ S : f (x) = (f1 (x), . . . , fk (x)). If U and V are open semi-algebraic subsets of Rn and f : U → V is a Nash function (i.e. f is a Nash function if viewed as a map f : U → Rn ) and the inverse f −1 exists and it is a Nash function, then f is called a Nash diffeomorphism. Definition 2.3. A Nash submanifold X of Rn is a semi-algebraic set which is also an analytic manifold. More precisely, X is a Nash submanifold of Rn of dimension d, if the following holds. For every point x ∈ X there exists 1. an open semi-algebraic set Ωx ⊆ Rn such that 0 ∈ Ωx , 0 2. an open semi-algebraic set Ω0x ⊆ Rn such that x ∈ Ωx , and 0 3. a Nash diffeomorphism ϕ : Ωx → Ωx such that ϕ(0) = x and ϕ((Rd × {0}) ∩ Ωx ) = X ∩ Ω0x , where Rd × {0} = {(x1 , . . . , xd , 0, . . . , 0) ∈ Rn | (x1 , . . . , xd ) ∈ Rd }. The sets Ωx will be called coordinate charts and the maps ϕ will be called coordinate maps. We will denote the dimension of X by dim(X). A map f : X → R is said to be Nash function on X, if for every x ∈ X, the map g : Ωx 3 (x1 , . . . , xn ) 7→ f (ϕ(x1 , . . . , xd , 0, . . . , 0)) is a Nash function. We denote the set of Nash functions of the form f : X → R by N (X). A map f : X → Rk is said to be a Nash function, if every coordinate function of f is a Nash function. If M, N are two Nash manifolds, N ⊆ Rk , then a map f : M → N is said to be Nash function, if f viewed as a map f : M → Rk is a Nash function. The Nash function f : M → N is said to be a Nash diffeomorphism, if its inverse f −1 : N → M exists and it is a Nash map. We will view each Nash manifold M ⊆ Rn as a topological space, with topology inherited from Rn : we will say that a subset U of M is open if there exists an open subset V of Rn such that V ∩ M = U . In fact, the definition of Nash manifolds implies that Nash manifolds are embedded analytic sub-manifolds of some Euclidian space (see [9] for the background material on smooth and analytic manifolds). The topology described above is just the corresponding topology for embedded submanifolds. If Nash manifolds are viewed as analytic embedded submanifolds of some Euclidian space, then Nash maps defined on them are analytic functions on manifolds. Further, the Zariski closure of a Nash submanifold has, by [8, Proposition 8.4.1], the following property: Proposition 2.4. Let X ⊆ Rn be a semi-algebraically connected Nash submanifold of dimension d. Then Z-cl(X) is an irreducible real affine variety of dimension d. Recall that the dimension of an affine variety equals the Krull dimension of the ring of polynomial functions of this variety. If X is as in Proposition 2.4, then the ring of polynomial functions on Z-cl(X) equals R[X1 , . . . , Xn ]/I(X), where I(X) = {P ∈ R[X1 , . . . , Xn ] | P |X = 0} is the ideal of polynomials vanishing on X. We will call R[X1 , . . . , Xn ]/I(X) the ring of polynomials on X and denote by R[X]. If X is connected, then, by Proposition 2.4, Z-cl(X) is an irreducible affine variety and thus R[X] is an integral domain. Hence, according to [20], the Krull dimension of R[X] equals its transcendence degree. Corollary 2.5. If X ⊆ Rn is a semi-algebraically connected Nash submanifold of dimension d, then trdeg R[X] = d. Further, we will need the following important property of Nash functions on Nash

NASH REALIZATIONS

5

manifolds. Proposition 2.6. [8, Proposition 8.18] Let X ⊆ Rn be a Nash submanifold which is semi-algebraically connected, let U be a non-empty open subset of X, and let f : X → R be a Nash function. If f |U = 0 then f = 0. The proposition above implies that if two Nash functions are equal on an open neighborhood, then they are equal everywhere. Proposition 2.6 also implies the following. Proposition 2.7. The ring N (X) of Nash functions defined on a semi-algebraically connected Nash manifold X is an integral domain. Indeed, assume f1 f2 = 0 for some f1 , f2 ∈ N (X) and f1 (z) 6= 0 for some z ∈ X. Because f1 is Nash and thus continuous, the set U = {x ∈ X | f1 (x) 6= 0} is open. Further, since z ∈ U , U is non-empty. Then, f1 f2 = 0 implies that f2 |U = 0, and hence by Proposition 2.6, f2 = 0. Finally, we will need the following corollary of Lemma 2.1. Proposition 2.8. Let X be a Nash submanifold of Rn which is semi-algebraically connected. Let f1 , . . . , fk ∈ N (X) and q ∈ N (S) where S is a Nash submanifold of Rk . Assume that for all x ∈ X, (f1 (x), . . . , fk (x)) ∈ S and consider the function g(x) = q(f1 (x), . . . , fk (x)), x ∈ X. Then g ∈ N (X) and g is algebraic over R[f1 , . . . , fk ]. In particular, the ring N (X) of Nash functions on X is algebraic over the ring R[X] of polynomials on X and trdeg N (X) = trdeg R[X] = dim(X). The proof of Proposition 2.8 is presented in Appendix A. 3. Nash systems. Nash systems are dynamical systems with inputs and outputs such that their state-spaces are Nash submanifolds and their dynamics and output functions are specified by Nash functions. We consider Nash systems with the following fixed input- and output-spaces. The input-space U is a subset of Rm . The output-space is Rr . As the space of input functions for a Nash system with a given input-space U we consider the set U pc of piecewise-constant functions u : [0, Tu ] → U , where Tu ∈ [0, +∞) is the maximal finite time instance for which u is defined. Any u ∈ U pc can be identified with a finite sequence u = (α1 , t1 ) · · · (αn , tn ) where Pi Pi+1 αi ∈ U , ti ∈ [0, ∞) for i = 1, . . . , n. Then, for t ∈ [ j=0 tj , j=0 tj ), u(t) = αi+1 ∈ U Pn for i = 0, 1, . . . , n − 1, t0 = 0, and u( j=0 tj ) = αn . Note that the same u can have several different representations as a sequence of tuples (α, t) where α ∈ U ,P t ∈ [0, ∞). n Every u = (α1 , t1 ) · · · (αn , tn ) ∈ U pc has a time-domain [0, Tu ] where Tu = j=1 tj . If u = (α1 , t1 ) · · · (αn , tn ), v = (β1 , s1 ) · · · (βk , sk ) ∈ U pc then by (u)(v) ∈ U pc we denote the input (α1 , t1 ) · · · (αn , tn )(β1 , s1 ) · · · (βk , sk ). The empty input e is the input with Te = 0. Notice that u(α, 0)v and uv represent the same piecewise-constant input for any u, v ∈ U pc , α ∈ U . To express that the input u was applied The empty input e is the input with Te = 0. only on the time-domain [0, t] ⊆ [0, Tu ] we write a subindex [0, t] to u like u[0,t] . Definition 3.1. A Nash system Σ with an input-space U ⊆ Rm and an outputspace Rr is a quadruple (X, f, h, x0 ) where (i) the state-space X is a Nash submanifold of Rn which is semi-algebraically connected, (ii) the dynamics of the system is given by x(t) ˙ = f (x(t), u(t)) for an input u ∈ U pc , where f : X × U → Rn is such that for every input value α ∈ U the function fα : X 3 x 7→ f (x, α) ∈ Rn is a Nash function. (iii) the output of the system is specified by a Nash function h : X → Rr , (iv) x0 = x(0) ∈ X is the initial state of Σ.

6

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

Note that Nash system Σ can be considered a triple Σ = (X, f, h) once one does not specify the initial state. Such definition is sufficient e.g. for the definition of observability properties of Nash systems, see Section 4.3.2. However, we keep Definition 3.1 because of the need to choose an initial state to solve the realization problem. Example 3.2. The following simplified model of glycolysis and lactate production in Lactococcus lactis is adopted from the supplements of [40]: x˙ 1 = 0.3592x−1.2906 x0.2168 Glc1.1287 − 0.3115x2.17 AT P 0.8152 , 4 1 1 x˙ 2 = 0.3115x2.17 AT P 0.8152 − 0.4698x1.0297 Pi0.2377 , 1 2 x˙ 3 = 0.9396x1.0297 Pi0.2377 + 1.1452x3.5453 − 2.167x2.1649 , 2 4 3 x˙ 4 = 2.167x2.1649 − 0.3592x−1.2906 x0.2168 Glc1.1287 − 1.1452x3.5453 − 0.2087x0.0002 3 4 4 4 1 −0.9375x0.8744 x0.0991 P i−0.0005 , 2 4 x˙ 5 = 0.3592x−1.2906 x0.2168 Glc1.1287 − 0.0417x0.6202 x0.9264 − 1.3258x1.5255 4 5 2 5 1 +0.9375x0.8744 x0.0991 P i−0.0005 , 2 4 x˙ 6 = 0.0417x0.6202 x0.9264 . 5 2 Here x1 , . . . , x6 are the concentrations of the respective metabolites (G6P, FBP, PGA3, PEP, pyruvate, lactate). The state-space X = (0, +∞)6 is a Nash submanifold of R6 which is semi-algebraically connected. The initial state can be chosen as x1 (0) = · · · = x6 (0) = 1. The output function is the outflow of the pathway, i.e. h(x1 , . . . , x6 ) = x6 . The inputs are the concentrations of Glc, AT P and P i of external glucose, ATP, and inorganic phosphate, respectively. Note that this system is neither polynomial nor rational, but still Nash. Indeed, h is linear and hence Nash function of the state. The right-hand sides of the differential equations are linear combinations of terms of the form xq11 · · · xq66 where qi , i = 1, . . . , 6 is a rational number. Since Nash functions are closed under linear combination, n multiplication and division, it remains to be shown that the map g(x) = x d , where x ∈ (0, +∞) and n, d ∈ N, is Nash. Since g(x) is analytic, it is sufficient to show that g is semi-algebraic. For all x ∈ (0, +∞) and for all n, d ∈ N it holds that y = g(x) if and only if P (y, x) = 0, where P is the polynomial defined as P (Y, X) = Y d − X n . Therefore, g is a semi-algebraic function. Note that in the example above, we did not allow zero concentrations. From a modeling point of view, this is not ideal, as in many biochemical applications zero concentrations may occur. The reason for this is that our definition of Nash systems requires the state-spaces to be Nash submanifolds without boundaries. If we allowed a state-space of the form X = [0, ∞)6 , then this state-space would be a Nash manifold with boundary. This is an inherent limitation of the class of Nash systems. At this stage, it is not clear if our results can be extended to systems, state-spaces of which are manifolds with boundaries. Perhaps one needs to follow a different approach to deal with such cases. Definition 3.3. The state trajectory of a Nash system Σ = (X, f, h, x0 ) corresponding to an input U pc 3 u = (α1 , t1 ) · · · (αk , tk ) : [0, Tu ] → U is a continuous piecewise-differentiable function xΣ (.; x0 , u) : [0, Tu ] → X such that xΣ (0; x0 , u) = x0 and d xΣ (t; x0 , u) = f (xΣ (t; x0 , u), u(t)) dt Pi Pi+1 for t ∈ ( j=0 tj , j=0 tj ), i = 0, . . . , k − 1, t0 = 0.

(3.1)

NASH REALIZATIONS

7

Let Σ = (X, f, h, x0 ) be a Nash system with an input-space U . For any α ∈ U , there exists a maximal T > 0 such that there exists a unique trajectory of Σ corresponding to the constant input u = (α, T ). This follows from the fact that f (x, α) is smooth with respect to x for every α ∈ U . Nevertheless, a unique solution of (3.1) need not exist for all u ∈ U pc . In order to deal with this phenomenon we introduce the notion of admissible inputs for a Nash system. Definition 3.4. A set U pc (Σ) of admissible inputs for a Nash system Σ = (X, f, h, x0 ) is a subset of U pc such that for all u ∈ U pc (Σ) there exists a unique trajectory xΣ (·; x0 , u) : [0, Tu ] → X of Σ corresponding to the input u. 4. Realization problem for Nash systems. In this section we formulate the realization problem for the class of Nash systems. First, in Section 4.1, we recall the notion of a response map adopted from [28, 29] and define when a Nash system is a realization of a response map. Next, in Section 4.2 and Section 4.3 we define minimality and semi-algebraic reachability and observability for Nash systems. These notions are necessary for formulating the realization problem for Nash systems. Finally, in Section 4.4 we state the realization problem for Nash systems. 4.1. Response maps. In order to formulate the realization problem, we need a formal definition of the input-output (external) behavior of the system. For our purposes, the response maps introduced in [28, 29] represent the right formalization of the external behavior of Nash systems. Below we recall their definition and basic properties. g Definition 4.1. A subset U pc of the set U pc is called a set of admissible inputs if: g g (i) ∀u ∈ U pc ∀t ∈ [0, Tu ] : u[0,t] ∈ U pc , g g (ii) ∀u ∈ U pc ∀α ∈ U ∃t > 0 : (u)(α, t) ∈ U pc , g (iii) ∀u = (α1 , t1 ) · · · (αk , tk ) ∈ U pc ∃δ > 0 ∀ti ∈ [0, ti + δ], i = 1, . . . , k : g u = (α1 , t1 ) · · · (αk , tk ) ∈ U pc . The first two conditions of the definition above assure the well-definedness of the derivative at switching times of differentiable functions which are defined on a set of admissible inputs, see Definition 5.1. The third condition allows us to prove Proposition 5.2 which is crucial for the proof of existence of Nash realizations. g In the sequel, unless stated otherwise, U pc denotes a set of admissible inputs. r g Definition 4.2. We say that a function ϕ : U pc → R is a response map if for every sequence of input values α1 , . . . , αk ∈ U , k > 0, the function ϕα1 ,...,αk , to be defined below, is analytic. Let us denote by Tα1 ,...,αk the set of all k-tuples (t1 , . . . , tk ) ∈ [0, +∞)k such that g the input (α1 , t1 )(α2 , t2 ) · · · (αk , tk ) belongs to U pc . We define the map ϕα1 ,...,αk : r Tα1 ,...,αk → R as follows: for all (t1 , . . . , tk ) ∈ Tα1 ,...,αk , ϕα1 ,...,αk (t1 , . . . , tk ) = ϕ((α1 , t1 ) · · · (αk , tk )). r r g g We denote the set of all response maps ϕ : U pc → R by A(U pc → R ). Let us point out the main difference between response maps and input-output maps, both describing external behavior of systems. Let u : [0, Tu ] → Rm be an input to a system Σ and let y : [0, Tu ] → Rr be the corresponding output. Then a response map of Σ maps u to y(Tu ) while an input-output map of Σ maps u to y. r g Consider a response map p : U pc → R and denote its coordinate functions g g by pi : U pc → R, i = 1, . . . , r, i.e. ∀u ∈ U pc : p(u) = (p1 (u), . . . , pr (u)). Then

8

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

g pi ∈ A(U pc → R) for all i = 1, . . . , r. We conclude this section by defining when a Nash system describes (i.e. realizes) a certain response map. r g Definition 4.3. Consider a response map p : U pc → R . A Nash system Σ = (X, f, h, x0 ) such that g g U pc ⊆ U pc (Σ) and p(u) = h(xΣ (Tu ; x0 , u)) for all u ∈ U pc is called a Nash realization of p. That is, Σ is a realization of p, if (a) all the inputs for which p is defined are admissible for Σ, i.e. there exists a unique state trajectory corresponding to those inputs, and (b) the value of p for the input u equals the output of Σ for u at time Tu . 4.2. Minimal Nash realizations. To state the realization problem for Nash systems formally we define the notion of dimension of a Nash system and the notion of minimal Nash realizations of response maps. These definitions are presented below. Definition 4.4. By the dimension of a Nash system Σ = (X, f, h, x0 ) we mean the dimension of its state-space X, where X is viewed as a semi-algebraically connected Nash submanifold of Rn . We denote the dimension of Σ by dim(Σ) and thus dim(Σ) = dim(X). The dimension of Σ is related to the number of state variables as follows. If dim(Σ) = d then there exists a local semi-algebraic coordinate transformation around any x ∈ X into Rd , i.e. locally, Σ can be described by d state variables. Notice that our definition of dimension is similar to the classical definition for nonlinear systems [16, 35]. Definition 4.5. We say that a Nash realization Σ = (X, f, h, x0 ) of a response map p is a minimal Nash realization of p if for any Nash realization Σ0 of p it holds that dim(Σ) ≤ dim(Σ0 ). Note that if a Nash system is minimal according to the definition above, then this does not automatically imply that the number of states of that system is the smallest one among all the possible Nash realizations of the same input-output map. The reason for this is that we define the dimension of a Nash system as the dimension of its state-space. That is, if Σ is a d-dimensional Nash system, then its state-space is a d-dimensional Nash submanifold of the n-dimensional space Rn . Hence, its description will use n variables and n can be larger than d. As it was pointed out before, if the dimension of the system is d, then there exists a local coordinate transformation such that locally the system dynamics can be described by d variables. However, for the description of the global dynamics, one still needs n state variables. More precisely, in 0 each coordinate chart (set Ωx , x ∈ X according to our notation), the system dynamics can be described by d coordinates, as long as the state stays in the coordinate chart. Once the state enters another coordinate chart, another description using d variables becomes valid. However, in order to describe the change of coordinate charts and thus describe the dynamics fully, one has to use all n original state variables. This means that there might exist a minimal realization whose state-space is a d-dimensional submanifold of the n-dimensional space Rn , and a non-minimal realization whose state space is a d < d1 -dimensional submanifold of the n1 -dimensional space Rn1 , such that n > n1 . In this case, the global description of the dynamics of the minimal system will use n states, while the global description of the non-minimal system will use n1 < n states. However, the local description of the minimal system will use less variables than the local description of the non-minimal system. Below we present an example to illustrate this point.

9

NASH REALIZATIONS

Example 4.6. Consider the Nash system Σ = (X, f, h, x0 ) where p U = R, X = {(x1 , x2 , x3 ) ∈ R3 | √ x21 + x22 +√x23 = 1, x1√= x2 }, h(x1 , x2 , x3 ) = 1 + x23 (x21 + x22 ), f ((x1 , x2 , x3 ), u) = ( 0.5ux3 , 0.5ux3 , − 2ux1 )T and x0 = (0, 0, 1). It is easy to see that X is a Nash manifold and fu and h are Nash functions 3 and√ dim(X) √ = 1. Indeed, define the sets3 Ω1 = {(x1 , x2 , x3 ) ∈√R |√x3 > 0, x1 ∈ (− 0.5, 0.5)}, Ω2 = {(x1 , x2 , x3 ) ∈ R | x3 < 0, x1 ∈ (− 0.5, 0.5)}, Ω3 = {(x1 , x2 , x3 ) ∈ R3 | x1 > 0, x3 ∈ (−1, 1)}, Ω4 =√{(x1 ,√x2 , x3 ) ∈ R3 | x1 < 0, x3 ∈ (−1, 1)} and let Ω = {(x1 , x2 , x3 ) ∈ R3 | x1 ∈ (− 0.5, 0.5)}. Define ψi : Ω → Ωi , i = 1, . . . , 4 as follows: p  if i = 1 1 − 2x21 ) (x , x , 1 1  p   2 (x if i = 2 p1 ) p1 , x1 , − 1 − 2x √ ψi (x1 , x2 , x3 ) = 2 ), 2 ), 2x )  ( 0.5(1 − 2x 0.5(1 − 2x if i = 3 1  1 1 p p √  (− 0.5(1 − 2x21 ), − (0.5(1 − 2x21 ), 2x1 ) if i = 4 0

0

For each p = (p1 , p2 , p3 ) ∈ X, define the sets Ωp , Ωp ⊆ R2 and the map ϕp : Ωp → Ωp :  Ω1 if p ∈ Ω1    Ω if p ∈ Ω 0 2 2 √ √ Ωp = 0.5, 0.5, Ω if p = (  3  √ √ 0)  Ω4 if p = (− 0.5, − 0.5, 0) √ √  (−√0.5 + p1 , √0.5 + p1 ) × R2 if p ∈ Ω√ 1 ∪ Ω2 √ √ √ Ωp = (− 0.5 + p3 , 0.5 + p3 ) × R2 if p = ( 0.5, 0, 5, 0) or p = (− 0.5, − 0.5, 0)  ψi (p − (x1 , x2 , x3 )) if p ∈ Ω√ 1, 2  i , i =√ ψ3 ((p3 , p1 , p2 ) − (x1 , x2 , x3 )) if p = ( √ 0.5, 0.5, ϕp (x1 , x2 , x3 ) = √ 0)  ψ4 ((p3 , p1 , p2 ) − (x1 , x2 , x3 )) if p = (− 0.5, − 0.5, 0) It is then easy to see that ϕp and ϕ−1 p are Nash functions, and hence that ϕp is a Nash 0 diffeomorphism. Moreover, ϕp (0, 0, 0) = p and ϕp ((R × {0} × {0}) ∩ Ωp ) = X ∩ Ωp . Hence, X is indeed a Nash manifold and it is of dimension d = 1. Because the maps Ωp 3 (x1 , x2 , x3 ) 7→ fu (ϕp (x1 , x2 , x3 )) and Ωp 3 (x1 , x2 , x3 ) 7→ h(ϕp (x1 , x2 , x3 )) are obviously Nash, one concludes that fu and h are Nash functions. Notice that the response map of Σ can be described by s Z Tu Z Tu 2 2 p(u) = 1 + cos ( u(s)ds) sin ( u(s)ds), 0

0

for all u ∈ U pc , since xΣ (Tu ; x0 , u) = (

Z p 0, 5 sin( 0

Tu

u(s)ds),

Z p 0, 5 sin( 0

Tu

Z u(s)ds), cos(

Tu

u(s)ds)) 0

satisfies the differential equation of the system. Clearly, dim(Σ) = dim(X) = d = 1 and hence Σ is a minimal realization. Notice that one uses three state variables to describe the system. However, if (x1 , x2 , x3 ) ∈ (Ω1 ∪ Ω2 ) ∩ X, then one p the dynamics of the sysp variable (x1 ) is sufficient to√describe tem: h(x1 , x2 , x3 ) = p1 + 2x21 (1 − 2x21 ) and x˙ 1 = 0.5u 1 − 2x21 for (x1 , x2 , x3 ) ∈ √ Ω1 and x˙ 1 = − 0.5u 1 − 2x21 for (x1 , x2 , x3 ) ∈ Ω2 . Similar local equations in one variable can be used if (x1 , x2 , x3 ) ∈ (Ωi ∩ X), i = 3, 4.

10

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

ˆ 1 , x2 ) = ˆ ˆ ˆ ˆ ˆ0 ) where X ˆ = R2 , x ˆ0 = (0, 1), h(x p Consider now the system Σ = (X, f , h, x 2 2 1 + x2 x1 , fˆ(x1 , x2 , u) = (ux2 , −ux1 ). It is easy to see that Z xΣˆ (Tu ; x ˆ0 , u) = (sin(

Tu

Z

Tu

u(s)ds), cos( 0

u(s)ds)) 0

ˆ is a realization of p. and hence Σ ˆ = 2 and hence Σ ˆ is not a minimal Nash realization of p. HowNote that dim(Σ) ˆ can be described by 2 ever, while the global description of Σ involves 3 variables, Σ variables. Note that the phenomenon above is not entirely a matter of choosing the wrong definition of minimality. One expects the minimal system to be reachable and observˆ is not reachable, and hence it is able in some sense. We will also show that in fact Σ not a good candidate for being a minimal realization. Note that the situation described by Example 4.6 can also be encountered when considering analytic or smooth systems on manifolds, it is not specific to Nash systems. 4.3. Reachability and observability. In order to present the realization problem for Nash systems, we have to choose a suitable definition of reachability and observability for Nash systems. While there are standard definitions for both concepts, these definitions do not take into account the algebraic structure of the systems of interest. As a result, the precise definition of observability and reachability varies with the system class. For Nash systems, we will use semi-algebraic reachability and semialgebraic observability instead of reachability and observability. We present these notions in Subsections 4.3.1 and 4.3.2. We chose these two concepts because they seemed to be useful for characterizing minimality of Nash systems and because they are closely related to the standard concepts of reachability and observability. It remains a topic of future research whether semi-algebraic reachability and observability are also useful for control design or state estimation. 4.3.1. Semi-algebraic reachability. We define semi-algebraic reachability of Nash realizations by a slight modification of the usual concept of reachability. Namely, instead of the requirement that the whole state-space is reachable from the initial state, we only require that the set of reachable states is, in a sense, a sufficiently large subset of the state-space. r g Definition 4.7. Let p : U pc → R be a response map and let Σ = (X, f, h, x0 ) be a Nash realization of p. We denote by R(x0 ) the set of states of Σ reachable from g x0 by the inputs of U pc , i.e. g R(x0 ) = {xΣ (Tu ; x0 , u) | u ∈ U pc }. g Since U pc ⊆ U pc (Σ), the reachable set R(x0 ) is potentially smaller than the set of all states of Σ reachable from x0 by the inputs of U pc (Σ). Definition 4.8. We say that a Nash realization Σ = (X, f, h, x0 ) of a response r g map p : U pc → R is semi-algebraically reachable, if no non-zero element of N (X) vanishes on the set of states of Σ reachable from x0 , i.e. if ∀g ∈ N (X) : (g = 0 on R(x0 ) ⇒ g = 0) .

11

NASH REALIZATIONS

ˆ from Example 4.6 is not semi-algebraically reachExample 4.9. The system Σ able: it is clear that for any u ∈ U pc , xΣˆ (Tu ; x ˆ0 , u) belongs to {(x, y) ∈ R2 | x2 + y 2 = 2 2 1}. Hence, the polynomial P (x, y) = x − y − 1 is zero on the set of states reachable ˆ is not semi-algebraically from x ˆ0 . Since P 6= 0 and P ∈ N (R2 ), it follows that Σ reachable. If a Nash system Σ is reachable, i.e. if every state of Σ can be reached by a suitable g input, then Σ is also semi-algebraically reachable provided that U pc = U pc (Σ). Recall from [35, Section 4.3] that a nonlinear smooth system is called accessible if the set of reachable states contains an open set. Accessibility of nonlinear systems admits a characterization in terms of the rank of the Lie-algebra generated by the vector fields of the system. Since the definition of accessibility and the corresponding Lie-rank condition can directly be applied to Nash systems, the following proposition implies that Lie-rank condition yields a sufficient condition for semi-algebraic reachability. Proposition 4.10. Let Σ = (X, f, h, x0 ) be a Nash realization of a response map r g p:U pc → R . If Σ is accessible, i.e. if there exists a non-empty open subset S of X such that S ⊆ R(x0 ) then Σ is semi-algebraically reachable. Proof. Let g ∈ N (X) be arbitrary. Because S ⊆ R(x0 ), g = 0 on R(x0 ) implies g = 0 on S. Since S is a non-empty open subset of X and since the Nash submanifold X is semi-algebraically connected, from Proposition 2.6 it follows that g = 0. Thus, Σ is semi-algebraically reachable. We are now in position to show that the system Σ from Example 4.6 is semialgebraically reachable. Example 4.11. It is easy to see that the state trajectory of Σ from the initial state x0 = (0, 0, 1) equals Z Tu Z Tu Z Tu √ √ xΣ (Tu ; u, x0 ) = ( 0.5 sin( u(s)ds), 0.5 sin( u(s)ds), cos( u(s)ds)). 0

0

0

Hence, states arising from constant inputs already yield the whole state-space X, i.e. R(x0 ) = X and in particular, R(x0 ) contains a subset which is open in the topology of X. Hence, by Proposition 4.10 Σ is semi-algebraically reachable. The converse of Proposition 4.10 is not true, as demonstrated by the following example. That is, it can happen that the system is semi-algebraically reachable, even though the set of reachable states has an empty interior. Example 4.12. Consider the system Σ = (X, f, h, x0 ) where U = R, X = R2 , x0 = (0, 1), ∀u ∈ U , f ((x1 , x2 ), u) = (1, x2 ) and h(x1 , x2 ) = x1 + x22 . It is easy to see that xΣ (Tu ; x0 , u) = (Tu , eTu ). In other words, R(x0 ) = {(t, et ) | t ∈ [0, +∞)}. Since R(x0 ) is a graph of a smooth function of time, it cannot contain an open subset of R2 . We will show that Σ is semi-algebraically reachable. To this end, consider a Nash function 0 6= g ∈ N (R2 ) and assume that g|R(x0 ) = 0, that is, g(t, et ) = 0 for all t ≥ 0. As g is an analytic function, the map R 3 t 7→ g(t, et ) is analytic and hence g(t, et ) = 0 should hold for all t ∈ R. Since g is a semi-algebraic function, we can apply Lemma 2.1 to it. Then there exist finitely many polynomials gi and semi-algebraic sets Si , Sd i = 1, . . . , d, such that Si are disjoint, gi (x, Y ) 6= 0 for all x ∈ Si , i=1 Si = R2 and gi (x, g(x)) = 0 for all x ∈ Si . Define Q(X1 , X2 , Y ) = g1 (X1 , X2 , Y ) · · · gd (X1 , X2 , Y ). It is clear that Q 6= 0. It then follows that Q(x, g(x)) = 0 on R2 . If Q(X1 , X2 , 0) = 0, then for some k > 0, Q(X1 , X2 , Y ) = Y k R(X1 , X2 , Y ), where R(X1 , X2 , Y ) is a polynomial such that R(X1 , X2 , 0) 6= 0. It then follows that (g(x))k R(x, g(x)) = 0 on R2 . Since the ring of entire analytic functions is an integral domain, it follows that either g = 0 or R(x, g(x)) = 0 for all x. Since g 6= 0, it then follows that R(x, g(x)) =

12

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

0 for all x. Hence, by possibly replacing Q by R, we can assume that E(X1 , X2 ) = Q(X1 , X2 , 0) 6= 0. Then from Q(x, g(x)) = 0 it follows that E(t, et ) = Q(t, et , 0) = 0 for all t ∈ R. In other words, the function et is an algebraic function. But the exponential function is known to be transcendental, hence this is a contradiction. That is, there cannot exists a non-zero Nash function g ∈ N (R2 ) such that g vanishes on R(x0 ). 4.3.2. Semi-algebraic observability. Below we define the concept of semialgebraic observability for Nash systems. The proposed concept was inspired by the notion of algebraic observability of [4, 28, 29]. Since the to be defined concepts do not depend on the initial state, in this section, unless stated otherwise, a Nash system will mean a tuple Σ = (X, f, h) without the initial state. First, we will define the concept of observation algebra for Nash systems. Definition 4.13. The observation algebra Aobs (Σ) of a Nash system Σ = (X, f, h) is the smallest subalgebra of N (X) which contains hi , i = 1, . . . , r and which is closed under taking Lie-derivatives with respect to the vector fields fα : X 3 x 7→ f (x, α) ∈ Rn , α ∈ U . The definition of the observation algebra of a Nash system is analogous to the one of rational and polynomial systems presented in [3, 4, 41]. Definition 4.14. We say that a Nash system Σ = (X, f, h) is semi-algebraically observable if trdeg Aobs (Σ) = trdeg N (X). We will say that a Nash system Σ = (X, f, h, x0 ) equipped with an initial state is semi-algebraically observable, if (X, f, h) is semi-algebraically observable. That is, Σ is semi-algebraically observable, if the algebra of all Nash functions on X is algebraic over the observation algebra of Σ. This means that any Nash function is a solution of a polynomial equation, coefficients of which are elements of the observation algebra. Another way to gain intuition for the concept of semi-algebraic observability is to use the analogy with linear system theory. For linear systems, the observation algebra corresponds to the algebra generated by the rows of the observability matrix, and the transcendence degree corresponds to the rank of the observability matrix. Hence, for linear systems, semi-algebraic observability boils down to the Kalman rank condition for observability. In other words, the notion of semi-algebraic observability can be viewed as a generalization of the Kalman rank condition for observability. To illustrate the concept further, we will review the system of Example 4.6. Example 4.15. Consider the system Σ = (X, f, h, x0 ) from Example 4.6. It is easy to see that the Nash manifold X is connected and hence by Lemma 2.8, trdeg N (X) = dim(X) = 1. Hence, if we can show that trdeg Aobs (Σ) = 1, it then follows that trdeg Aobs (Σ) = trdeg N (X). To this end, it is enough to show that there exists no polynomial in one variable 0 6= Q ∈ R[T ], such that Q(h) = 0 on X. Indeed, Q(h(x)) = 0 for all x ∈ X implies that h(x) is a root of Q. But Q can have only finitely many roots and hence h(x), should take only finitely many values. However, this is clearly not the case. In order to connect the concept of semi-algebraic observability to the standard definition of observability, we introduce the notion of strong semi-algebraic observability. If Σ is strongly semi-algebraically observable, then it is observable in the sense that there are no indistinguishable states, and it is also semi-algebraically observable. The details are as follows. ash Definition 4.16. The Nash extension AN obs (Σ) of the observation algebra Aobs (Σ) of a Nash system Σ = (X, f, h) is defined as the sub-algebra of N (X) gener-

NASH REALIZATIONS

13

ated by Nash functions g : X → R of the form ∀x ∈ X : g(x) = q(ϕ1 (x), . . . , ϕk (x)) where • ϕ1 , . . . , ϕk ∈ Aobs (Σ), k > 0 • q ∈ N (S) for some Nash manifold S ⊆ Rk , such that ∀x ∈ X : (ϕ1 (x), . . . , ϕk (x)) ∈ S. The concept of Nash extension formulated in Definition 4.16 is a slight generalization of a similar concept from [27]. Definition 4.17. We say that a Nash system Σ = (X, f, h) is strongly semiash algebraically observable if AN obs (Σ) = N (X). We say that a Nash system (X, f, h, x0 ) equipped with an initial state is strongly semi-algebraically observable, if (X, f, h) is strongly semi-algebraically observable. Intuitively, strong semi-algebraic observability means that any Nash function can be obtained by substituting some elements of the observation algebra into a Nash function which is defined on some semi-algebraic set. This is analogous to the definition of algebraic observability in [4, 28, 29]; algebraic observability means that any polynomial (rational) function defined on the state-space can be obtained by substituting the elements of the observation algebra into a polynomial (rational) function. Strong semi-algebraic observability and semi-algebraic observability are related as follows. Proposition 4.18. If a Nash system is strongly semi-algebraically observable then it is semi-algebraically observable. Note that the converse of Proposition 4.18 is not true, as it is demonstrated by the following counter-example. Example 4.19. Consider the system Σ from Example 4.12. It then follows that Lfα (h) = 1 + 2x22 and Lfα1 · · · Lfαk h = 2k x22 for all k > 1. Hence, Aobs (Σ) = R[x1 , x22 ]. Since the polynomials x1 and x22 are trivially algebraically independent, it follows that trdeg Aobs (Σ) = trdeg N (R2 ) = dim(R2 ) = 2. Assume now that Σ is strongly semi-algebraically observable. Then there exists a Nash manifold S ⊆ Rk and a Nash function g ∈ N (S), and k polynomials R1 , . . . , Rk , such that x2 = g(R1 (x1 , x22 ), . . . , R2 (x1 , x22 )) for all x1 , x2 ∈ R. Since the composition of Nash functions with polynomials is again a Nash function, without loss of generality we can assume that x2 = g(x1 , x22 ) and k = 2. Since dxd 1 x2 = dxd 1 g(x1 , z) = 0, it then follows that we can take k = 1 and assume that x2 = g(x22 ). But then 1 = g(12 ) = g((−1)2 ) = ash −1, a contradiction. Hence, clearly x2 cannot belong to AN obs (Σ). This means that Σ cannot be strongly semi-algebraically observable. Notice that if in the counter-example above we take X = R × (0, +∞), then the thus obtained system is strongly semi-algebraically observable. In fact, we can show the following. Proposition 4.20. If Σ = (X, f, h) is semi-algebraically observable, then there exists an open semi-algebraic subset V of X such that the system ΣV = (V, f V , h|V ) is strongly semi-algebraically observable. We denote by fαV the restriction of fα to V , and by h|V the restriction of h to V . Recall from [8] that any open semi-algebraic subset of a Nash submanifold of dimension d is itself a Nash submanifold of dimension d. A system-theoretic interpretation of strong semi-algebraic observability is provided by the following proposition and its corollary.

14

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

Proposition 4.21. If a Nash system Σ is strongly semi-algebraically observable, then any two states x1 6= x2 of Σ are distinguishable by an element of Aobs (Σ), i.e. ∃g ∈ Aobs (Σ) : g(x1 ) 6= g(x2 ). Corollary 4.22. Let Σ = (X, f, h) be a strongly semi-algebraically observable Nash system. Then Σ is observable in the sense that it has no indistinguishable states, i.e. if x1 6= x2 ∈ X then there exists u ∈ U pc such that xΣ (Tu ; x1 , u) and xΣ (Tu ; x2 , u) are both defined and h(xΣ (Tu ; x1 , u)) 6= h(xΣ (Tu ; x2 , u)). Corollary 4.22 implies that differential-geometric conditions for observability of nonlinear systems, see for example [16], also yield necessary conditions for strong semialgebraic observability of Nash systems. The proofs of Proposition 4.18, Proposition 4.20, Proposition 4.21 and Corollary 4.22 can be found in Appendix A. Recall [16] that Σ = (X, f, h) is said to satisfy the observability rank condition, if the rank of the observability co-distribution ∀x ∈ X : ∆(x) = span{d(Lfα1 · · · Lfαk h(x) | α1 , . . . , αk ∈ U, k ≥ 0}, equals dim(X), i.e. if dim(∆(x)) = dim(X) for all x ∈ X. Recall from [15, Theorem 3.1] that if Σ satisfies the observability rank condition, then it is locally weakly observable. By local weak observability we mean that for any state x there is an open neighborhood V of x in X, such that for any open subset W ⊆ V and any x1 ∈ W , there exists an u ∈ U pc such that xΣ (.; u, x), xΣ (.; u, x1 ) exist and ∀t ∈ [0, Tu ] : xΣ (t; u, x), xΣ (t; u, x1 ) ∈ W and h(xΣ (Tu ; u, x1 )) 6= h(xΣ (Tu ; u, x2 )). For analytic systems, and hence for Nash systems, local weak observability is known ([7, Theorem 2.1]) to be equivalent to weak observabibility: Σ is called weakly observable, if for any x ∈ X there exists an open neighborhood V of x such that ∀x1 ∈ V : ∃u ∈ U pc : xΣ (.; u, x), xΣ (.; u, x1 ) exist, and h(xΣ (Tu ; u, x1 )) 6= h(xΣ (Tu ; u, x)). From [15, Theorem 3.11] it follows that if the system is locally weakly observable, then there exists an open subset W of X such that for all x ∈ W : dim(∆(x)) = dim(X). Proposition 4.23. Let Σ = (X, f, h) be a Nash system. Then Σ is semialgebraically observable if and only if there exists x ∈ X, dim(∆(x)) = dim(X). In particular, if Σ is locally weakly observable or if Σ satisfies the observability rank condition, then it is semi-algebraically observable. The proof of Proposition 4.23 is presented in Appendix A. Note that while local weak observability (observability rank condition) implies semi-algebraic observability, the converse is not true in general, see the counter-example below. Example 4.24. Consider the system Σ from Example 4.19. It is easy to see that h(xΣ (Tu ; (0, −x), u) = h(xΣ (Tu ; (0, x), u) = Tu + e2Tu x2 . Hence, in any neighborhood V of (0, 0), any two elements of the form (0, x), (0, −x) are indistinguishable. Hence, Σ is not locally weakly observable. It follows that the observability codistribution is ∆(m) = span{dx(m), ydy(m)} for any m ∈ R2 . But dim(∆(0)) = dim(span{dx(0)}) = 1 6= 2, hence Σ does not satisfy the observability rank condition. However, in Example 4.19 we have already shown that Σ is semi-algebraically observable. 4.4. Problem formulation. The realization problem for Nash systems consists of three subproblems which concern the existence and minimality of Nash realizations and the algorithms for computing them. r g Problem 4.25. Let p : U pc → R be a response map. The Nash realization problem for p consists of the following subproblems: Existence Determine necessary and sufficient conditions for the existence of a Nash realization of p. Minimality Determine necessary and sufficient conditions for the existence of a minimal Nash realization of p. Determine whether a minimal Nash realization of p is

NASH REALIZATIONS

15

unique in any sense (for example, up to an isomorphism). Determine the relationship between minimality, observability and reachability. Realization algorithm Formulate the algorithms for computing a (minimal) Nash realization of p from finite data directly obtainable from p. In addition, formulate algorithms for checking minimality of a Nash realization and for transforming a Nash realization of p to a minimal one. In this paper we deal with the existence and minimality problems for Nash realizations. We do not address the algorithmic part of the realization problem. 5. Main results. In this section we present the main results of the paper. We will start with conditions for existence of a Nash realization and we will continue by addressing minimality of Nash systems. 5.1. Existence of a Nash realization. In order to state necessary and sufficient conditions for the existence of a Nash realization of a response map we need to define a number of concepts for response maps. First, we recall the notion of a derivative of a response maps from [4, 28, 29]. g Definition 5.1 ([4]). Consider a function ϕ ∈ A(U pc → R). For any α ∈ U we define the derivative (Dα ϕ)(u) =

d g ϕ((u)(α, t))|t=0+ for all u ∈ U pc . dt

Next, we recall the fact that the set of all response maps forms an algebra, in fact it is an integral domain. g Proposition 5.2 ([28, 29]). The set A(U pc → R) with point-wise addition and multiplication forms a commutative algebra over R and it is an integral domain. Further, we recall the notion of observation algebra of a response map. r g Definition 5.3 ([4]). Let p : U pc → R be a response map. The observation g algebra Aobs (p) of p is the smallest subalgebra of the algebra A(U pc → R) which contains the components pi , i = 1, . . . , r of p, and which is closed with respect to the derivations Dα , α ∈ U . Thus, Aobs (p) is the smallest algebra such that p1 , . . . , pr ∈ Aobs (p), and if ϕ ∈ Aobs (p) then Dα ϕ ∈ Aobs (p) for all α ∈ U . Another key notion is the one of Nash extension of a finite subset of the set g A(U pc → R) of scalar-valued response maps. Recall that the concept of Nash extension of an infinite set of functions is introduced in Definition 4.16 for observation algebras. The following definition of Nash extension of a finite set of functions is then a modification of Definition 4.16. Definition 5.4. Let X be a Nash submanifold of Rn which is semi-algebraically g connected, let U pc be a set of admissible inputs, and let A = {ϕ1 , . . . , ϕn } be a subset g g of A(U → R). Assume that for all u ∈ U pc pc , (ϕ1 (u), . . . , ϕn (u)) ∈ X. The Nash N ash g extension A (X) of A with respect to X is the subalgebra of A(U pc → R) generated g by the maps U 3 u → 7 q(ϕ (u), . . . , ϕ (u)) ∈ R, where q ∈ N (X). pc 1 n Intuitively, the Nash extension of A is obtained by substituting the elements of A into Nash functions defined on X. Note that the set of substitutions into linear forms (polynomials) yields the linear space (algebra) generated by A. Then, the Nash extension of A can be thought of as the generalization of the notions of linear space and algebra generated by A. We are now able to state the following necessary and sufficient condition for the existence of a Nash realization.

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

16

r g Theorem 5.5. A response map p : U pc → R has a Nash realization if and only g if there exist a finite subset A = {ϕ1 , . . . , ϕn } of A(U pc → R) and a Nash submanifold n X ⊆ R which is semi-algebraically connected such that g (i) ∀u ∈ U pc : (ϕ1 (u), . . . , ϕn (u)) ∈ X, (ii) pi ∈ AN ash (X) for i = 1, . . . , r, (iii) ∀ϕ ∈ AN ash (X) ∀α ∈ U : Dα ϕ ∈ AN ash (X). The proof of Theorem 5.5 is presented in Section 6. Remark 5.6. The idea of the proof is very similar to that of [6, Theorem 4.6], where an analogous result for analytic systems on time scales is stated. Since the system class of [6] is much more general, the result of [6, Theorem 4.6] does not imply Theorem 5.5. More precisely, when applied to continuous-time systems, the results of [6, Theorem 4.6] provide necessary and sufficient conditions for the existence of an analytic realization. The corresponding conditions are the same as in Theorem 5.5, except for the usage of the analytic closure C ω (ϕ1 , . . . , ϕk ) = {f (ϕ1 , . . . , ϕk ) | f : g V → R analytic, V ⊆ Rk open, ∀u ∈ U pc : (ϕ1 (u), . . . , ϕk (u)) ∈ V } instead of the k N ash Nash closure A (R ). Since AN ash (Rk ) ⊆ C ω (ϕ1 , . . . , ϕk ), if the conditions (i)–(iii) of Theorem 5.5 hold, then the conditions of [6, Theorem 4.6] hold and an analytic realization exists. Since any Nash system is an analytic system on time scales (see [6] for the definition), this is not at all surprising. However, if the conditions of [6, Theorem 4.6] hold, then it is not necessarily true that (i)–(iii) of Theorem 5.5 are satisfied. Similar ideas occur in [25], when attempting to view smooth manifolds as spectra of the algebras of smooth functions. Example 5.7. Consider the response map s Z Tu Z Tu u(t)dt) cos2 ( u(t)dt), p(u) = 1 + sin2 ( 0

0

where u ∈ U pc are real-valued piecewise-constant inputs. Let us denote ϕ1 (u) = RT RT sin( 0 u u(t)dt) and ϕ2 (u) = cos( 0 u u(t)dt). Then, for all u ∈ U pc (ϕ1 (u), ϕ2 (u)) ∈ {(x, y) ∈ R2 : x2 + y 2 = 1} =: X. Note that X is semi-algebraically connected Nash submanifold of R2 , see for example Example 4.6. p Since p(u) = 1 + ϕ21 (u)ϕ22 (u), it follows that p ∈ AN ash (X). Let α ∈ U be arbitrary. Then Z Tu Z Tu +τ Z Tu d [sin( u(t)dt + αdt)]|τ =0+ = α cos( u(t)dt) = αϕ2 (u) Dα ϕ1 (u) = dτ 0 Tu 0 and similarly Z Dα ϕ2 (u) = −α sin(

Tu

u(t)dt) = −αϕ1 (u).

0

Let ϕ ∈ AN ash (X) be arbitrary. Then ϕ(u) = q(ϕ1 (u), ϕ2 (u)), where q ∈ N (X) and ∂q ∂q Dα ϕ(u) = ∂x (ϕ1 (u), ϕ2 (u))Dα ϕ1 (u) + ∂y (ϕ1 (u), ϕ2 (u))Dα ϕ2 (u) = ∂q ∂q α ∂x (ϕ1 (u), ϕ2 (u))ϕ2 (u) − α ∂y (ϕ1 (u), ϕ2 (u))ϕ1 (u), where N ash

Dα ϕ ∈ A

(X).

∂q ∂q ∂x , ∂y

∈ N (X) and thus

NASH REALIZATIONS

17

Therefore, the map p satisfies the conditions of Theorem 5.5 and hence its Nash realization does exist. We construct the realization in the same way as it is done in the proof of Theorem 5.5. Namely, X = {(x, y) ∈ R2 : x2 + y 2 = 1} is considered the state space of the realization Σ. The dynamics of Σ is defined according to the derivatives Dα ϕ1 (u) = αϕ2 (u) and Dα ϕ1 (u) = −αϕ1 (u) as follows: x˙ = uy, y˙ = −ux. p The output function h(x, y) = 1 + x2 y 2 is based on the formula defining the dependence of p on ϕ1 and ϕ2 . Finally, the initial states are x(0) = 0, y(0) = 1 which follows from ϕ1 (e) = 0, ϕ2 (e) = 1, respectively. The conditions of Theorem 5.5 are rather difficult to check. Therefore, we will derive a simpler necessary condition for the existence of a Nash realization. In order g to state this condition, recall that A(U pc → R) is an integral domain, and hence so is the observation algebra of p. Therefore, the transcendence degree trdeg Aobs (p) of the observation algebra of p is well-defined. Theorem 5.8. Let Σ = (X, f, h, x0 ) be a Nash realization of a response map p : g U pc → Rr . Then, trdeg Aobs (p) ≤ dim(X). Hence, if there exists a Nash realization of a response map p then trdeg Aobs (p) < +∞. The proof of Theorem 5.8 is presented in Section 6. This necessary condition is analogous to the finite Hankel-rank condition for linear systems. In particular, the transcendence degree of Aobs (p) can be seen as a generalization of the rank of the Hankel-matrixqfor linear systems. RT RT Example 5.9. Consider the map p(u) = 1 + sin2 ( 0 u u(t)dt) cos2 ( 0 u u(t)dt) which is realizable by a Nash system, see Example 5.7. We show that trdeg Aobs (p) < +∞. More precisely, we show that trdeg Aobs (p) = 1 = dim(X), where X stands for the state space of the realization computed in Example p 5.7. First, let us compute Aobs (p). Recall that p = 1 + ϕ21 ϕ22 and ϕ21 + ϕ22 = 1. Because p ∈ Aobs (p), one gets ϕ21 ϕ22 ∈ Aobs (p). Because Dα (ϕ21 ϕ22 ) = 2αϕ1 ϕ32 − 2αϕ31 ϕ2 , one obtains ϕ1 ϕ32 − ϕ31 ϕ2 ∈ Aobs (p). Taking derivatives again, it is easy to see that Dα (ϕ1 ϕ32 − ϕ31 ϕ2 ) = α(ϕ41 + ϕ42 ) − 6αϕ21 ϕ22 . Thus, ϕ41 + ϕ42 ∈ Aobs (p). Further, Dα (ϕ41 + ϕ42 ) = −4α(ϕ1 ϕ32 − ϕ31 ϕ2 ) ∈ Aobs (p). Consider now the algebra R = R[ϕ21 ϕ22 , ϕ1 ϕ32 − ϕ31 ϕ2 , ϕ41 + ϕ42 ]. It then follows that R is closed under derivations Dα , i.e. if r ∈ R, then Dα (r) ∈ R. It is easy to see that Dα1 · · · Dαk p are 2 2 algebraic combinations of derivatives of ϕp 1 ϕ2 and (positive or negative) powers of p 2 2 1 + ϕ1 ϕ2 . Hence, R ⊆ Aobs (p) ⊆ R[ 1 + ϕ21 ϕ22 , √ 1 2 2 ]. Since both R and 1+ϕ1 ϕ2 p 1 2 2 R[ 1 + ϕ1 ϕ2 , √ ] are finitely generated, it then follows that their transcen2 2 1+ϕ1 ϕ2

dence degrees arepfinite. Hence, the transcendence degree of Aobs (p) is finite too. It is easy to see that 1 + ϕ21 ϕ22 , √ 1 2 2 are algebraic over R: for p(x, y) = y − x2 − 1, 1+ϕ1 ϕ2 p 2 2 2 2 p( 1 + ϕ1 ϕ2 , ϕ1 ϕ2 ) = 0. Similarly, for p(x, y) = x2 (y+1)−1, p( √ 1 2 2 , ϕ21 ϕ22 ) = 0. 1+ϕ1 ϕ2 p Hence, trdeg Aobs (p) ≤ trdeg R[ 1 + ϕ21 ϕ22 , √ 1 2 2 ] = trdeg R. Finally, notice that 1+ϕ1 ϕ2

ϕ41 + ϕ42 = (ϕ21 + ϕ22 )2 − 2ϕ21 ϕ22 = 1 − ϕ21 ϕ22 and that s = ϕ1 ϕ32 − ϕ31 ϕ2 = ϕ1 ϕ2 (ϕ21 − ϕ22 ) satisfies p(s, ϕ21 ϕ22 ) = 0 for p(x, y) = x2 − y(1 − 4y). Hence, R is algebraic over R[ϕ21 ϕ22 ] and hence trdeg R = trdeg Aobs (p) = 1. We conjecture that the finiteness of the transcendence degree of Aobs (p) is also a sufficient condition for the existence of a Nash realization, at least provided that p

18

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

is defined for inputs whose values vary only within a finite input-space U . However, this statement is still an open problem. What we can prove is the following. We will say that a response map p is locally realized by a Nash system Σ = (X, f, h, x0 ) if for any integer k > 0, α1 , . . . , αk ∈ U , there exists t¯i > 0, i = 1, . . . , k such that for all ti ∈ [0, t¯i ], i = 1, . . . , k, the input g v = (α1 , t1 ) · · · (αk , tk ) satisfies v ∈ U pc (Σ), and p(v) = h(xΣ (Tv ; x0 , v)). That is, Σ is a local realization of p, if for any piecewise-constant input with small enough switching times, the response of p to this input equals the response of Σ to u r g g this input. Further, if u ∈ U pc , then denote by pu the map pu : U pc → R , where u g g g U pc = {v ∈ U pc | uv ∈ U pc }, and g ∀v ∈ U pc

u

: pu (v) = p(uv).

It is easy to see that pu is again a response map. Theorem 5.10. Assume trdeg Aobs (p) = d < +∞ and U is finite. For any  > 0 g one can choose an input u ∈ U pc and a Nash system Σ such that Tu <  and Σ is a local realization of pu . The proof of Theorem 5.10 is presented in Section 6. Theorem 5.10 implies that if the observation algebra has a finite transcendence degree, then we can find a local realization of the response map p, if we restrict attention to inputs preceded by some input u. The latter u can be chosen to be defined on an arbitrarily small time interval. Note that as Tu is getting smaller, the responses of pu will get closer to those of p. Hence, pu can be seen as an approximation of p and thus Σ can be viewed as an approximation of a Nash realization of p. Theorem 5.10 thus indicates that the necessary condition of finite transcendence degree of the observation algebra is in fact very close to being a sufficient one as well. This gives hope that the existence results can be made more sharp and practically more useful. Note that polynomial systems as defined in [4], see the definition below, are special cases of Nash systems. Hence, sufficient conditions derived in [4] for the existence of polynomial realizations yield sufficient conditions for the existence of Nash realizations. In order to present these conditions and their straightforward application to Nash systems, we first recall the notion of a polynomial system. Definition 5.11 ([4]). A polynomial system Σ with an input-space U ⊆ Rm and an output-space Rr is a quadruple (X, f, h, x0 ) where (i) the state-space X is an irreducible variety in Rn , (ii) the dynamics of the system is given by x(t) ˙ = f (x(t), u(t)) for an input u ∈ U pc , where f : X × U → Rn is such that for every input value α ∈ U the components fα,i : X → R, i = 1, . . . , n of f (x, α) = (fα,1 (x), . . . , fα,n (x)) are polynomial functions on X (fα,i is the ith coordinate of the vector field fα : X 3 x 7→ f (x, α) ∈ Rn ), (iii) the output of the system is specified by the map h : X → Rr , the components h1 , . . . , hr of h are polynomial functions on X, (iv) x0 = x(0) ∈ X is the initial state of Σ. The results on the existence of polynomial realizations stated in the following theorem are derived in [4, Theorem 3] and in [41, Theorem 5.2]. r g Theorem 5.12. Let p : U pc → R be a response map. If Aobs (p) is a finitely generated R-algebra then p has a polynomial realization, and hence a Nash realization.

NASH REALIZATIONS

19

Note that if Aobs (p) is a finitely generated algebra, then its transcendence degree is finite. 5.2. Minimality of Nash systems. By Definition 4.5, a minimal Nash realization of a response map is a Nash realization such that the dimension of its state-space is the smallest one within the class of all Nash realizations of the same response map. As a direct consequence of Theorem 5.8 we can state the following sufficient condition for minimality. Theorem 5.13. If the dimension of a Nash realization Σ of a response map p equals trdeg Aobs (p), then Σ is a minimal Nash realization of p. Theorem 5.13 is not very appealing, because it does not relate minimality to systems-theoretic properties. However, it has the following corollary, which relates minimality to semi-algebraic reachability and semi-algebraic observability. Before stating this corollary, we state the following result, which is interesting on its own right. Theorem 5.14. A Nash realization Σ of a response map p is semi-algebraically reachable and semi-algebraically observable if and only if dim(Σ) = trdeg Aobs (p). The proof of Theorem 5.14 is presented in Section 7. Theorem 5.14 and Theorem 5.13 then yield the following sufficient condition for minimality. Theorem 5.15. If a Nash realization Σ of a response map p is semi-algebraically reachable and semi-algebraically observable then Σ is minimal. In particular, if Σ is semi-algebraically reachable and strongly semi-algebraically observable, then it is minimal. Note that both Theorem 5.13 and Theorem 5.15 formulate sufficient conditions for minimality of a Nash realization. At this point we do not know if any Nash realization can be converted to a semi-algebraically reachable and observable one. One way to show this is to show that for any response map p, which is realizable by a Nash system, there exists a Nash realization Σ such that dim(Σ) = trdeg Aobs (p). We do not know if this is true. However, Theorem 5.10 indicates that this is ‘almost’ the case. It is natural to wonder what the relationship is between minimality of Nash systems and minimality of analytical non-linear systems; it is clear that if a Nash system is minimal as an analytic system, then it is also minimal as a Nash system. The converse does not hold. Indeed, recall from [16] that a minimal non-linear analytic (smooth) system must be accessible and locally weakly observable. But the system presented in Example 4.19 is neither locally weakly observable nor it is accessible (see Example 4.6), and hence it cannot be minimal as a smooth system. However, it is semi-algebraically reachable and semi-algebraically observable and hence it is minimal as a Nash system (it has the smallest possible dimension among all Nash realizations of its input-output maps). 6. Proofs of the results on the existence of a Nash realization. In this section we present the proofs of the results presented in Subsection 5.1. We start with formulating the concepts of dual input-to-state maps and stating their properties. We then use these properties to present the proofs of Theorem 5.5, Theorem 5.8 and Theorem 5.12. 6.1. Dual input-to-state maps. Below we present the concept of dual inputto-state maps. These maps and their properties will be central to the proofs of the theorems on the existence of Nash realizations. For the same concepts in case of polynomial and rational systems see [4, 3, 28].

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

20

Definition 6.1. Let Σ = (X, f, h, x0 ) be a Nash realization of a response map r g p : U pc → R and let xΣ (·; x0 , u) be the trajectory of Σ corresponding to an input g g g u∈U pc . The map τΣ : U pc → X defined as τΣ (u) = xΣ (Tu ; x0 , u) for all u ∈ U pc is called the input-to-state map. By the dual input-to-state map of Σ we mean the map g τΣ∗ : N (X) → A(U pc → R) such that for every Nash function g ∈ N (X) it holds that ∗ τΣ (g) = g ◦ τΣ , i.e. ∗ g ∀u ∈ U pc : τΣ (g)(u) = g(τΣ (u)) = g(xΣ (Tu ; x0 , u)).

Intuitively, the dual input-to-state map of a Nash system Σ = (X, f, h, x0 ) maps each Nash function g on X to the response map which is generated by the Nash system Σ2 = (X, f, g, x0 ) where g is used instead of h as the output map. The dual inputto-state map plays a role similar to the observability Grammian of linear systems. In particular, it allows us to relate the properties of a response map p to the properties of the ring of Nash functions on the state-space X of a Nash realization Σ = (X, f, h, x0 ) of p. Proposition 6.2. Let Σ = (X, f, h, x0 ) be a Nash realization of a response map g g p : U pc → Rr and let τΣ∗ : N (X) → A(U pc → R) be as in Definition 6.1. Then, ∗ (i) ∀g ∈ N (X) ∀α ∈ U : Dα τΣ (g) = τΣ∗ (Lfα g), where Lfα g is the Lie-derivative of the map g with respect to the vector field fα : X 3 x 7→ f (x, α), and Dα is the derivation defined in Definition 5.1, (ii) Aobs (p) = τΣ∗ (Aobs (Σ)) ⊆ τΣ∗ (N (X)). g (iii) ∀g ∈ N (Rk ) ∀ϕ1 , . . . , ϕk ∈ N (X) ∀u ∈ U pc : τΣ∗ (g(ϕ1 , . . . , ϕk ))(u) = g(τΣ∗ (ϕ1 )(u), . . . , τΣ∗ (ϕk )(u)). The proof of Proposition 6.2 is analogous to the proof of [4, Lemma 1]. However, [4, Lemma 1] was stated only for polynomial systems and thus it is not directly applicable to our case. Hence, for the sake of completeness, we state the proof of Proposition 6.2 below. Proof. [Proof of Proposition 6.2] ∗ ∗ g (i) We show that ∀g ∈ N (X) ∀α ∈ U ∀u ∈ U pc : Dα τΣ (g)(u) = τΣ (Lfα g)(u). Let g g ∈ N (X), α ∈ U , u ∈ U pc be arbitrary. Then, Dα τΣ∗ (g)(u) = d = g(xΣ (Tu + t; x0 , (u)(α, t)))|t=0+ = dt  n  X d d = = g(xΣ (Tu + t; x0 , (u)(α, t))) xΣ,i (Tu + t; x0 , (u)(α, t)) dxi dt t=0+ i=1 =

n X d g(xΣ (Tu ; x0 , u))fi (xΣ (Tu ; x0 , u), α) = (Lfα g)(xΣ (Tu ; x0 , u)) = dxi i=1

= (Lfα g)(τΣ (u)) = τΣ∗ (Lfα g)(u). (ii) From Definition 4.13, Aobs (Σ) ⊆ N (X). To prove (ii) it is sufficient to prove that Aobs (p) = τΣ∗ (Aobs (Σ)). r g Because Σ = (X, f, h, x0 ) is a Nash realization of a response map p : U pc → R , g it holds that pi (u) = hi (τΣ (u)) for all u ∈ U pc and i = 1, . . . , r. Thus, pi = τΣ∗ (hi ) for all i = 1, . . . , r.

(6.1)

NASH REALIZATIONS

21

Further, by using (6.1) and repeatedly applying (i) we obtain ∀i = 1, . . . , r ∀k ∈ N ∀α1 , . . . , αk ∈ U : Dαk · · · Dα1 pi = τΣ∗ (Lfαk · · · Lfα1 hi ). Since Aobs (p) is generated by the elements of the set {pi , Dα1 · · · Dαk pi | i = 1, . . . , r; k ∈ N; α1 , . . . , αk ∈ U }, and Aobs (Σ) is generated by the elements of the set {hi , Lfα1 · · · · · · Lfαk hi | i = 1, . . . , r; k ∈ N; α1 , . . . , αk ∈ U }, and τΣ∗ is a homomorphism, we can conclude that Aobs (p) = τΣ∗ (Aobs (Σ)). (iii) This statement follows directly from the definition of τΣ∗ . g Corollary 6.3. Let Σ = (X, f, h, x0 ) be a Nash system such that U pc ⊆ U pc (Σ) ∗ g g for a set U pc of admissible inputs. Let τΣ : N (X) → A(U pc → R) be as in Definition 6.1. Then, ∀g ∈ N (X) ∀k ∈ N ∀α1 , . . . , αk ∈ U : Dαk · · · Dα1 τΣ∗ (g) = τΣ∗ (Lfαk · · · Lfα1 g). Corollary 6.4. Consider the same assumptions as in Proposition 6.2. Then trdeg Aobs (p) ≤ trdeg τΣ∗ (N (X)) and trdeg Aobs (p) ≤ trdeg Aobs (Σ). 6.2. Proof of Theorem 5.5. Proof. r g (⇒) Let Σ = (X, f, h, x0 ) be a Nash realization of a response map p : U pc → R g and assume that X ⊆ Rn . Let ϕi ∈ A(U pc → R), i = 1, . . . , n be the maps defined as g ϕi (u) = xΣ,i (Tu ; x0 , u) for all u ∈ U , i.e. ϕi (u) is the ith component of the state of pc g Σ at time Tu under the input u ∈ U pc . Let us define the set A = {ϕ1 , . . . , ϕn }. Since g (ϕ1 (u), . . . , ϕn (u)) = xΣ (Tu ; x0 , u) for all u ∈ U pc , the condition (i) of the theorem is satisfied. Because Σ is a Nash realization of p, hi ∈ N (X) and pi (u) = hi (ϕ1 (u), . . . , ϕn (u)) N ash g for i = 1, . . . , r and for all u ∈ U (X) for all pc . Hence, by Definition 5.4, pi ∈ A i = 1, . . . , r, i.e. (ii) is valid. Let ϕ ∈ AN ash (X) and α ∈ U be arbitrary. To prove (iii) we show that Dα ϕ ∈ N ash A (X). From Definition 5.4 and from the definition of the maps ϕi , i = 1, . . . , n, there exists q ∈ N (X) such that ϕ = q(ϕ1 , . . . , ϕn ) = τΣ∗ (q). Then, Proposition 6.2(i) implies that Dα ϕ = Dα τΣ∗ (q) = τΣ∗ (Lfα q) = (Lfα q)(ϕ1 , . . . , ϕn ). As Lfα q ∈ N (X), from the last equality above it follows that Dα ϕ ∈ AN ash (X). n g (⇐) Let A = {ϕ1 , . . . , ϕn } ⊆ A(U pc → R) and X ⊆ R be as in the theorem. We will define a Nash system Σ realizing p. From (ii) and (iii) of the theorem it follows g that for all i = 1, . . . , r, j = 1, . . . , n, α ∈ U , and u ∈ U pc there exist hi , fα,j ∈ N (X) such that pi (u) = hi (ϕ1 (u), . . . , ϕn (u)) and Dα ϕj (u) = fα,j (ϕ1 (u), . . . , ϕn (u)). Consider the Nash system Σ = (X, f, h, x0 ) such that for all x ∈ X and for all α ∈ U it holds that f (x, α) = (fα,1 (x), . . . , fα,n (x)), h(x) = (h1 (x), . . . , hr (x)) and x0 = (ϕ1 (e), . . . , ϕn (e)). Recall that e stands for the empty input. We prove that this system is a Nash realization of p.

22

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

Let pri : X → R : (x1 , . . . , xn ) 7→ xi for i = 1, . . . , n be the restriction of the projection on the ith coordinate on X. First, notice that if pri : X → R is the ith projection, i.e. pri (x1 , . . . , xn ) = xi , i = 1, . . . , n, then Lfα pri = fα,i . g For every input u = (α1 , t1 ) · · · (αk , tk ) ∈ U pc we define the map ψu : [0, Tu ] → X as ψu (t) = (ϕ1 , . . . , ϕn )(u[0,t] ) for t ∈ [0, Tu ]. Then, for i = 1, . . . , k and for t ∈ Pi−1 Pi ( j=0 tj , j=0 tj ) (note that t0 = 0) it holds that d d d ψ˙ u (t) = ψu (t) = ψu (t + τ )|τ =0+ = (ϕ1 (u[0,t+τ ] ), . . . , ϕn (u[0,t+τ ] ))|τ =0+ dt dτ dτ i−1 X d = (ϕ1 ((α1 , t1 ) · · · (αi−1 , ti−1 )(αi , t + τ − tj )), . . . dτ j=0 . . . , ϕn ((α1 , t1 ) · · · (αi−1 , ti−1 )(αi , t + τ −

i−1 X

tj )))|τ =0+

j=0

= (Dαi ϕ1 (u[0,t] ), . . . , Dαi ϕn (u[0,t] )) = (Lfαi pr1 )(ϕ1 (u[0,t] , . . . , ϕn (u[0,t] )), . . . , (Lfαi prn )(ϕ1 (u[0,t] , . . . , ϕn (u[0,t] ))  = (Lfαi pr1 )(ψu (t)), . . . , (Lfαi prn )(ψu (t))



= (f1 (ψu (t), αi ), . . . , fn (ψu (t), αi )) = f (ψu (t), αi ) = f (ψu (t), u(t)) and ψu (0) = (ϕ1 (e), . . . , ϕn (e)) = x0 . Therefore, by Definition 3.3, ψu is the state trajectory of Σ corresponding to the input g u, i.e. ψu (t) = xΣ (t; x0 , u) for t ∈ [0, Tu ]. Hence, for arbitrary u ∈ U pc and for all i = 1, . . . , r, hi (xΣ (Tu ; x0 , u)) = hi (ψu (Tu )) = hi ((ϕ1 , . . . , ϕn )(u)) = hi (ϕ1 (u), . . . , ϕn (u)) = pi (u). Thus, Σ is indeed a Nash realization of p. 6.3. Proofs of Theorems 5.8 and 5.10. r g Proof. Let Σ = (X, f, h, x0 ) be a Nash realization of a response map p : U pc → R . By Proposition 2.8, N (X) is algebraic over R[X] and trdeg N (X) = trdeg R[X] = dim(X). By applying Corollary 6.4, we derive that trdeg Aobs (p) ≤ trdeg τΣ∗ (N (X)) ≤ trdeg N (X) = dim(X). Proof. [Proof of Theorem 5.10] Assume that ϕ1 , ϕ2 , . . . , ϕd ∈ Aobs (p) is a transcendence basis of Aobs (p). Then there exist polynomials Qα,i (Y, X1 , . . . , Xd ), i = 1, 2, . . . , d, α ∈ U and Pj (Y, X1 , . . . , Xd ), j = 1, 2, . . . , r in d + 1 variables such that Qα,i (Dα ϕi , ϕ1 , ϕ2 , . . . , ϕd ) = 0 for all α ∈ U, i = 1, 2 . . . , d Pj (pj , ϕ1 , ϕ2 , . . . , ϕd ) = 0 for all j = 1, 2, . . . , r

(6.2)

By choosing Qα,i and Pj such that they are of minimal degree with respect to Y , we d d can assume that dY Qα,i (Dα ϕi , ϕ1 , ϕ2 , . . . , ϕd ) 6= 0 and dY Pj (pj , ϕ1 , ϕ2 , . . . , ϕd ) 6= 0

23

NASH REALIZATIONS

g viewed as elements of A(U pc → R). Let R be the product of d and dY Pj (pj , ϕ1 , ϕ2 , . . . , ϕd )’s, i.e. R = (Πα∈U Πdi=1

d dY

Qα,i (Dα ϕi , ϕ1 , ϕ2 , . . . , ϕd )’s

d d Qα,i (Dα ϕi , ϕ1 , ϕ2 , . . . , ϕd ))(Πrj=1 Pj (pj , ϕ1 , ϕ2 , . . . , ϕd )). dY dY

g Clearly, R is a response map and as A(U pc → R) is an integral domain, R 6= 0, since g none of the factors of R is zero. Hence, there exists v ∈ U pc such that R(v) 6= 0. Assume that v = (β1 , s1 ) · · · (βl , sl ), β1 , . . . , βl ∈ U , s1 , . . . , sl ∈ [0, +∞). Then, from the properties of admissible inputs it follows that there exists a connected open subset V of [0, +∞)l such that (0, . . . , 0), (s1 , . . . , sl ) ∈ V and for any (τ1 , . . . , τl ) ∈ V , g (β1 , τ1 ) · · · (βl , τl ) ∈ U pc , and the map g : V 3 (τ1 , . . . , τl ) 7→ R((β1 , τ1 ) · · · (βk , τl )) is analytic. Since g(s1 , . . . , sl ) 6= 0, it then follows the set S of those elements τ ∈ V such that g(τ ) = 0 is closed and its interior is empty. Hence, for any (τ1 , . . . , τl ) ∈ V \ d S, u = (β1 , τ1 ) · · · (βl , τl ), R(u) 6= 0 and thus dY Qα,i (Dα ϕi (u), ϕ1 (u), ϕ2 (u), . . . , ϕd (u)) 6= d 0 and dY Pj (pj (u), ϕ1 (u), ϕ2 (u), . . . , ϕd (u)) 6= 0. Note that since 0 ∈ V and S has empty interior, for any  > 0, there exists a choice (τ1 , . . . , τk ), such that Tu = Pk i=1 τi < . We apply now the implicit function theorem [8, Corollary 2.9.8] to Qα,i and Pj at 0 points (yα,i , x0 ) = (Dα ϕi (u), ϕ1 (u), . . . , ϕd (u)) and (yj0 , x0 ) = (pj (u), ϕ1 (u), . . . , ϕd (u)), 0 where x = (ϕ1 (u), . . . , ϕd (u)). It then follows that there exist a semi-algebraic open subset X ⊆ Rd and Nash maps fα,i , hj ∈ N (X), i = 1, . . . , d, and semi-algebraic open subsets Vα,i , Vj of R α ∈ U , j = 1, . . . , r such that 0 1. x0 ∈ X, yα,i ∈ Vα,j , yj0 ∈ Vj , α ∈ U , i = 1, . . . , d, j = 1, . . . , r, and 2. for every (yα,i , x) ∈ Vα,i × X, (yj , x) ∈ Vj × X, Pj (yj , x) = 0 ⇐⇒ yj = hj (x) and Qα,i (yα,i , x) = 0 ⇐⇒ yα,i = fα,i (x), α ∈ U , i = 1, . . . , d, j = 1, . . . , r. We will show that Σ = (X, f, h, x0 ), where f (x, α) = (fα,1 (x), . . . , fα,d (x))T and h(x) = (h1 (x), . . . , hr (x)) is a local realization of pu . To this end, take an arbitrary sequence of inputs α1 , . . . , αk ∈ U , k ≥ 0. For any t = (t1 , . . . , tl ) ∈ [0, +∞)k , denote by v(t) the input (α1 , t1 ) · · · (αk , tk ). From the properties of admissible inputs it follows that there exists an open subset W ⊆ g [0, +∞)l , 0 ∈ W such that for all t ∈ W , uv(t) ∈ U pc . Define G : W 3 t 7→ (ϕ1 (uv(t), . . . , ϕd (uv(t)) and Gα,i : W 3 t 7→ Dα ϕi (uv(t)) and Hj : W 3 t 7→ 0 and Hj (0) = yj0 , i = 1, . . . , d, j = 1, . . . , r, pj (uv(t)). Notice G(0) = x0 , Gα,i (0) = yα,i 0 0 α ∈ U . Hence, there exists a subset W of W of the form W = [0, r)k for some r > 0, 0 such that for all t ∈ W , G(t) ⊆ X, Gα,i (t) ∈ Vα,i , Qα,i (Gα,i (t), G(t)) = 0 and Hj (t) ∈ Vj , Pj (Hj (t), G(t)) = 0, i = 1, . . . , d, j = 1, . . . , r, α ∈ U . Hence, Hj (t) = 0 hj (G(t)) and Gα,i (t) = fα,i (G(t)) for all t ∈ W and i = 1, . . . , d, j = 1, . . . , r, α ∈ U . Pk 0 For every t = (t1 , . . . , tk ) ∈ W , define xΣ (.; x0 , v(t)) : [0, i=1 ti ] 7→ X as follows: Pi Pi Pi+1 xΣ (t; x0 , v(t)) = G((t1 , . . . , ti , t − j=1 tj , 0, . . . , 0)) if t ∈ ( j=1 tj , j=1 tj ] for some 0 0 i = 0, . . . , k − 1. Note that if (t1 , . . . , tk ) ∈ W , then (t1 , . . . , ti , s, 0, . . . , 0) ∈ W for d s < ti+1 , i = 0, . . . , k − 1. Notice that the lth entry of the vector dt xΣ (t; x0 , v(t)) Pi Pi+1 Pi equals Gαi+1 ,l (t1 , . . . , ti , t − j=1 ti , 0, . . . , 0), if t ∈ ( j=1 tj , j=1 tj ]. Hence, x(t) = xΣ (t; x0 , v(t)) satisfies the differential equation x(t) ˙ = f (x(t), v(t)(t)) with initial state 0 x(0) = x0 . It follows from the definitions above that for any t ∈ W , p(uv(t)) = h(xΣ (Tv(t) ; x, v(t))).

24

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

Since the choice of α1 , . . . , αk was arbitrary, it follows that Σ is indeed a local realization of pu . 7. Proof of the results on minimality of Nash realizations. In this section we present the proof of the sufficient condition for the minimality presented in Subsection 5.2. To this end, we need the following technical results characterizing semi-algebraic reachability. Proposition 7.1. Let Σ = (X, f, h, x0 ) be a Nash realization of a response map g p : U pc → Rr . Consider the following statements. (i) Σ is semi-algebraically reachable, g (ii) the dual input-to-state map τΣ∗ : N (X) → A(U pc → R) is injective, ∗ (iii) the ideal Ker τΣ is the zero ideal in N (X). (iv) trdeg Aobs (p) = trdeg Aobs (Σ). The statements (i), (ii) and (iii) are equivalent. The statement (ii) implies (iv). If Σ is semi-algebraically observable, then (iv) and (ii) are equivalent. Proof. The equivalence (ii) ⇔ (iii) and the implication (ii) =⇒ (iv) are obvious. If Σ is semi-algebraically observable, then trdeg N (X) = trdeg Aobs (Σ). Hence, if trdeg Aobs (Σ) = trdeg Aobs (p), then trdeg N (X) = trdeg τΣ∗ (N (X)) < +∞. From [42, Chapter II, Theorem 28,29] it follows that the latter is possible only if τΣ∗ is injective. It is left to prove (i) ⇔ (iii). By Definition 4.7 and 4.8, Σ is semi-algebraically reachable if and only if   g ∀g ∈ N (X) : ∀u ∈ U pc : g(xΣ (Tu ; x0 , u)) = 0 ⇒ g = 0. Using the notation of τΣ∗ , see Definition 6.1, one can reformulate this characterization as ∀g ∈ N (X) : τΣ∗ (g) = 0 ⇒ g = 0. Thus, Σ is semi-algebraically reachable if and only if Ker τΣ∗ = (0) ⊆ N (X). Proof. [Proof of Theorem 5.14] Let Σ = (X, f, h, x0 ) be a semi-algebraically reachable and semi-algebraically observable Nash realization of a response map p. By Proposition 2.8, dim(Σ) = dim(X) = trdeg N (X). Then, because Σ is semialgebraically observable, see Definition 4.14, dim(Σ) = trdeg Aobs (Σ).

(7.1)

Since Σ is semi-algebraically reachable, it follows by Proposition 7.1 that trdeg Aobs (Σ) = trdeg Aobs (p).

(7.2)

Finally, by (7.1) and (7.2), dim(Σ) = trdeg Aobs (p). Conversely, assume that dim(Σ) = trdeg Aobs (p). If Σ is not semi-algebraically observable, then trdeg Aobs (Σ) < dim(Σ). Since τΣ∗ (Aobs (Σ)) = Aobs (p), it follows that trdeg Aobs (p) ≤ trdeg Aobs (Σ) < dim(Σ), which is a contradiction. Hence, Σ must be semi-algebraically observable. From Proposition 7.1 and trdeg Aobs (p) = dim(Σ) = trdeg Aobs (Σ) it then also follows that Σ is semi-algebraically reachable. 8. Conclusions. In this paper we have introduced the class of Nash systems and we have formulated the realization problem for it. Our motivation for studying Nash systems is their relevance both for theory and applications. The presented results on realization theory of Nash systems are not yet complete. Further research aims at extending these results to necessary and sufficient conditions for the existence of minimal Nash realizations, characterizing the relations

NASH REALIZATIONS

25

between Nash realizations of the same response maps, and obtaining a realization algorithm. In addition, we would like to investigate the applications of the results of the paper to system identification, model reduction, observer and control design of Nash systems. Another possible research direction is to explore the extension of our results to systems which are defined on Nash manifolds with boundaries. Acknowledgment The authors thank the anonymous referees for their numerous useful remarks and suggestions. REFERENCES [1] J. Baillieul. The geometry of homogeneous polynomial dynamical systems. Nonlinear Anal., Theory, Meth. and Appl., 4(5):879–900, 1980. [2] J. Baillieul. Controllability and observability of polynomial dynamical systems. Nonlinear Anal., Theory, Meth. and Appl., 5(5):543–552, 1981. [3] Z. Bartosiewicz. Rational systems and observation fields. Systems and Control Letters, 9:379– 386, 1987. [4] Z. Bartosiewicz. Minimal polynomial realizations. Mathematics of control, signals, and systems, 1:227–237, 1988. [5] Z. Bartosiewicz, U. Kotta, E. Pawluszewicz, M. Wyrwas. Control systems on regular time scales and their differential rings. MCSS, 22, pp 185 –201, 2011. [6] Z. Bartosiewicz, E. Pawluszewicz. Realizations of nonlinear control systems on time scales. IEEE Transactions on Automatic Control, 53, pp. 571 – 575, 2008. [7] Z. Bartosiewicz. Local observability of non-linear systems. Systems & Control Letters, 25, pp. 295 – 298, 1995. [8] J. Bochnak, M. Coste, and M.-F. Roy. Real algebraic geometry. Springer-Verlag, Berlin Heidelberg, 1998. [9] M.W. Boothby. An introduction to differentiable manifolds and riemannian geometry. Academic Press, 1975. [10] A. Casagrande, K. Casey, R. Falchi, C. Piazza, B. Ruperti, G. Vizzotto, and B. Mishra. Translating time-course gene expression profiles into semi-algebraic hybrid automata via dimensionality reduction. Algebraic Biology, pp 51–65, 2007. [11] A. Casagrande, P. Corvaja, C. Piazza, and B. Mishra. Composing semi-algebraic o-minimal automata. HSCC, pp 668–671, 2007. [12] K. Forsman. Hybrid control systems and comprehensive Gr¨ obner bases. In Tech. Report, Automatic Control, Link¨ oping, pp 4092 – 4097, 1994. [13] M. Fr¨ anzle. What will be eventually true of polynomial hybrid automata? In 4th Int. Symposium on Theoretical Aspects of Computer Software TACS, pp 340–359, 2001. [14] J.J. Heijnen. Approximative kinetic formats used in metabolic network modeling. Biotechnology and Bioengineering, 91(5):534–545, 2005. [15] R. Hermann, A.J. Krener. Nonlinear Controllability and Observability. IEEE Transactions on Automatic Control, AC-22(5):728 – 740, 1977. [16] A. Isidori. Nonlinear control systems. Springer, 3rd edition, 2001. [17] B. Jakubczyk. Existence and uniqueness of realizations of nonlinear systems. SIAM J. Control and Optimization, 18(4):455 – 471, 1980. [18] E. Klipp, R. Herg, A. Kowald, C. Wierling, and H. Lehrach. Systems biology in practice: Concepts, implementation and application. WILEY-VCH, Weinheim, 2005. [19] A. Krener. Bilinear and non-linear realization of input-output maps. SIAM Journal on Control, 13(4):827–834, 1975. [20] E. Kunz. Introduction to commutative algebra and algebraic geometry. Birkh¨ auser, Boston, 1985. [21] G. Lafferriere, S. Pappas, and G.J. Sastry. O-minimal hybrid systems. Mathematics of Control, Signals, and Systems, 13(1):1–21, 2000. [22] V. Mysore and B. Mishra. Algorithmic algebraic model checking IV: Characterization of metabolic networks. Algebraic Biology, pp 170–184, 2007. [23] V. Mysore, C. Piazza, and B. Mishra. Algorithmic algebraic model checking II: Decidability of semi-algebraic model checking and its applications to systems biology. In ATVA, pp 217–233, 2005. [24] J. Nash. The imbedding problem for Riemannian manifolds. Ann. Mathematics, 63, pp 20 – 63, 1956. [25] J. Nestruev. Smooth manifolds and observables Springer 2003, ISBN 0387955437.

26

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

[26] J. Nˇ emcov´ a. Rational Systems in Control and System Theory. Ph.D. thesis, Vrije Universiteit, Amsterdam, 2009. [27] J. Nˇ emcov´ a and M. Petreczky and J.H. van Schuppen. Realization Theory of Nash Systems. In Proceedings of CDC, Shanghai, China, pp. 5935-5940, 2009. [28] J. Nˇ emcov´ a and J.H. van Schuppen. Realization theory for rational systems: The existence of rational realizations. SIAM J. Control Optim., 48(4): 2840–2856, 2009. [29] J. Nˇ emcov´ a and J.H. van Schuppen. Realization theory for rational systems: Minimal rational realizations. Acta Applicandae Mathematicae, 110: 605–626, 2010. [30] M. Petreczky and R. Vidal. Realization theory for semi-algebraic hybrid systems. In Hybrid Systems: Computation and Control, volume 4981 of LNCS. Springer-Verlag, 2008. [31] C. Piazza, M. Antoniotti, V. Mysore, A. Policriti, F. Winkler, and B. Mishra. Algorithmic algebraic model checking I: Challenges from systems biology. In CAV, 2005. [32] M.A. Savageau. Biochemical systems analysis. I. some mathematical properties of the rate law for the component enzymatic reactions. J. Theor. Biol., 25(3):365–369, 1969. [33] M.A. Savageau. Biochemical systems analysis. II. the steady-state solutions for an n-pool system using a power-law approximation. J. Theor. Biol., 25(3):370–379, 1969. [34] M.A. Savageau. Biochemical systems analysis. A study of function and design in molecular biology. Addison-Wesley, Reading, Massachusetts, 1976. [35] E.D. Sontag. Mathematical control theory: Deterministic finite dimensional systems. Text in applied mathematics 6. Springer, 2nd edition, 1998. [36] E.D. Sontag. Polynomial Response Maps, volume 13 of Lecture Notes in Control and Information Sciences. Springer Verlag, 1979. [37] J. Vera, E. Balsa-Canto, P. Wellstead, J.R. Banga, and O. Wolkenhauer. Power-law models of signal transduction pathways. Cell. Signal., 19(7):1531–1541, 2007. [38] D. Visser, R. van der Heijden, K. Mauch, M. Reuss, and S. Heijnen. Tendency modeling: A new approach to obtain simplified kinetic models of metabolism applied to saccharomyces cerevisiae. Metabolic Engineering, 2:252–275, 2000. [39] E.O. Voit. Computational analysis of biochemical systems: A practical guide for biochemists and molecular biologists. Cambridge University Press, Cambridge, 2000. [40] E.O. Voit, J. Almeida, S. Marino, R. Lall, G. Goel, A.R. Neves, and H. Santos. Regulation of glycolysis in lactococcus lactis: an unfinished systems biological case study. IEE Proc.-Syst. Biol., 153(4):286–298, 2006. [41] Y. Wang and E.D. Sontag. Algebraic differential equations and rational control systems. SIAM J. Control Optim., 30(5):1126–1149, 1992. [42] O. Zariski and P. Samuel. Commutative algebra I, II. Springer, 1958.

Appendix A. Proofs of the technical results. Proof. [Proof of Proposition 2.8] Let g(x) = q(f1 (x), . . . , fk (x)), x ∈ X be as in Proposition 2.8. We prove that g is algebraic over R[f1 , . . . , fk ]. By Lemma 2.1, there exist semi-algebraic subsets S1S, . . . , Sd of S and polynomials d gi ∈ R[X1 , . . . , Xk , Y ], i = 1, . . . , d such that S = i=1 Si , Si ∩ Sj = ∅ for all i 6= j ∈ {1, . . . , d} and such that for all s ∈ Si , gi (s, Y ) 6= 0 and gi (s, q(s)) = 0. For every x ∈ X there exists a unique index i(x) ∈ {1, . . . , d} such that ξ(x) = (f1 (x), . . . , fk (x)) ∈ Si(x) . Indeed, ξ(x) ∈ S, and hence ξ(x) belongs to one of the disjoint sets Si(x) . Then gi(x) (ξ(x), Y ) 6= 0 and gi(x) (ξ(x), q(ξ(x))) = gi(x) (ξ(x), g(x)) = 0. Consider the subset I = {i(x)|x ∈ X} ⊆ {1, . . . , d}. Assume that I = {i1 , . . . , il }. Consider the polynomial P = gi1 · · · gil ∈ R[X1 , . . . , Xk , Y ]. Then P (ξ(x), g(x)) = gi1 (ξ(x), g(x)) · · · gil (ξ(x), g(x)) = 0 for all x ∈ X. The proof is complete once we show P (ξ(x), Y ) 6= 0 for some x ∈ X. The fact that there exists x ∈ X such that P (ξ(x), Y ) 6= 0 is equivalent to the requirement that the polynomial P (f1 , . . . , fk , Y ) ∈ N (X)[Y ] with coefficients in N (X) is not identically zero. Since N (X) is an integral domain, it holds that the algebra N (X)[Y ] is also an integral domain, see [42, Vol.1, Ch.1.16]. Hence, if P (f1 , . . . , fk , Y ) = Πlj=1 gij (f1 , . . . , fk , Y ), where gij (f1 , . . . , fk , Y ) ∈ N (X)[Y ], equals zero as an element of N (X)[Y ] then there exists at least one i ∈ I such that gi (f1 , . . . , fn , Y ) = 0. The latter means that gi (f1 (x), . . . , fk (x), Y ) = gi (ξ(x), Y ) = 0 for all x ∈ X. But there exists at least one x ˆ such that i = i(ˆ x), i.e. such that

NASH REALIZATIONS

27

ξ(ˆ x) ∈ Si and gi (ξ(ˆ x), Y ) 6= 0 which is a contradiction. Therefore, we get that P (f1 , . . . , fn , Y ) 6= 0. To prove that the ring N (X) is algebraic over the polynomial ring R[X] we consider arbitrary f ∈ N (X) and repeat the proof above with • fi = pri , i = 1, . . . , n, where pri : X → R : x 7→ xi is the restriction of the projection map on the ith coordinate to X, • S = X, • g = q = f. Next, we present the proof of Proposition 4.18. To this end, we need the following auxiliary result. ash Proposition A.1. Let Σ = (X, f, h, x0 ) be a Nash system. Then AN obs (Σ) is N ash algebraic over Aobs (Σ) and consequently trdeg Aobs (Σ) = trdeg Aobs (Σ). ash Proof. Consider an arbitrary g ∈ AN obs (Σ). According to Definition 4.16 there exist k ∈ N, ϕ1 , . . . , ϕk ∈ Aobs (Σ), and q ∈ N (S) for some Nash manifold S ⊆ Rk , such that g(x) = q(ϕ1 (x), . . . , ϕk (x)) for all x ∈ X. From Proposition 2.8 it follows that g is algebraic over R[ϕ1 , . . . , ϕk ] ⊆ Aobs (Σ). ash N ash Because g ∈ AN obs (Σ) was arbitrary, Aobs (Σ) is algebraic over Aobs (Σ). This N ash implies that trdeg Aobs (Σ) = trdeg Aobs (Σ). Proof. [Proof of Proposition 4.18] As Σ is strongly semi-algebraically observable, ash we derive from Definition 4.17 that N (X) = AN obs (Σ). Consequently, trdeg N (X) = ash trdeg AN (Σ). Then, from Proposition A.1, trdeg N (X) = trdeg Aobs (Σ) which obs proves that Σ is semi-algebraically observable. Proof. [Proof of Proposition 4.20] Let dim(X) = d and assume that X ⊆ Rn . Let ϕ1 , . . . , ϕd be a transcendence basis of Aobs (Σ). Let x1 , . . . , xn be the coordinate functions of R[X], i.e. πi (x) yields the ith coordinate of a point x ∈ Rn , then xi is the restriction of πi to X. It is clear that π1 , . . . , πn ∈ N (X). Since Aobs (Σ) ⊆ N (X) and trdeg Aobs (Σ) = trdeg N (X), it then follows that π1 , . . . , πn are algebraic over Aobs (Σ), i.e. there exist non-zero polynomials Q1 , . . . , Qn ∈ R[T, T1 , . . . , Td ] such that Qi (πi , ϕ1 , . . . , ϕd ) = 0 for all x ∈ X, i = 1, . . . , n, where Qi (πi , ϕ1 , . . . , ϕd ) is viewed as a map X 3 x 7→ Qi (πi (x), ϕ1 (x), . . . , ϕd (x)). The latter map clearly belongs to N (X). By choosing Q1 , . . . , Qn to be of minimal degree with respect to the indeterminant d T , we can assume that dT Qi (πi , ϕ1 , . . . , ϕd ) 6= 0, i = 1, . . . , n. Consider the product d d R = dT Q1 (π1 , ϕ1 , . . . , ϕd ) · · · dT Qn (πn , ϕ1 , . . . , ϕd ) and notice that it can be viewed as an element of N (X). Since N (X) is integral domain, it then follows that R 6= 0 and hence there exists x ∈ X such that R(x) 6= 0. In other words, there exists x ∈ d Qi (πi (x), ϕ1 (x), . . . , ϕd (x)) 6= 0 X such that Qi (πi (x), ϕ1 (x), . . . , ϕd (x)) = 0 and dT n+d n for all i = 1, . . . , n. Define G : R → R by G(y, z) = (G1 (y, z), . . . , Gn (y, z)), Gi (y1 , . . . , yn , z) = Qi (yi , z), i = 1, . . . , n and let y0 = (π1 (x), . . . , πn (x)) = x and z0 = (ϕ1 (x), . . . , ϕd (x)). It then follows that the matrix [ dyd i Gj (y0 , z0 )]i,j=1,...,n is invertible. Hence, from the implicit function theorem [8, Corollary 2.9.8] it follows that there exist open semi-algebraic subsets W1 , W2 of Rn and Rd respectively and a Nash map ψ : W2 → W1 such that z0 ∈ W2 , y0 ∈ W1 and for any (y, z) ∈ W1 × W2 , G(y, z) = 0 ⇐⇒ y = ψ(z). Define V = {s ∈ X | s ∈ W1 and (ϕ1 (s), . . . , ϕd (s)) ∈ W2 }. Then V is an open semi-algebraic subset of X and x ∈ V . Moreover, for any s ∈ V , G(s, ϕ1 (s), . . . , ϕd (s)) = 0 and hence s = ψ(ϕ1 (s), . . . , ϕd (s)) for all s ∈ V .

28

ˇ ´ AND M. PETRECZKY AND J.H. VAN SCHUPPEN J. NEMCOV A

If V is not semi-algebraically connected, then we replace V by the biggest semialgebraically connected component of V which contains x. Notice that since V is a semi-algebraic and open subset of X, V is also a Nash submanifold of dimension d. If g ∈ N (V ), then g(s) = g(ψ(ϕ1 (s), . . . , ϕd (s))) for all s ∈ V . Take then S = W2 ∩ ψ −1 (V ); since S is open subset of Rd , it is a ddimensional Nash submanifold. It then follows that q = g ◦ ψ ∈ N (S) and g(s) = q(ϕ1 (s), . . . , ϕd (s)), for all s ∈ V . Proof. [Proof of Proposition 4.21] Let Σ = (X, f, h, x0 ) be a strongly semialgebraically observable Nash system. Because Nash functions on X distinguish the points of X, it holds that all states of Σ are distinguishable by the elements ash of AN obs (Σ) = N (X). That is, ash ∀x1 , x2 ∈ X, x1 6= x2 ∃g ∈ AN obs (Σ) : g(x1 ) 6= g(x2 ).

We prove that the same statement holds even if we can choose g only from Aobs (Σ). ash Let us assume that for some x1 6= x2 ∈ X there exists only g ∈ AN obs (Σ) \ Aobs (Σ) which distinguishes x1 and x2 . Let g be such a map, i.e. g ∈ N (X) and g(x1 ) 6= g(x2 ). Then, by Definition 4.16, there exist k ∈ N, ϕ1 , . . . , ϕk ∈ Aobs (Σ), and q ∈ N (S) for some Nash manifold S ⊆ Rk such that g(x1 ) = q(ϕ1 (x1 ), . . . , ϕk (x1 )) 6= q(ϕ1 (x2 ), . . . , ϕk (x2 )) = g(x2 ). If ϕi (x1 ) = ϕi (x2 ) for all i = 1, . . . , k then it would imply that q(ϕ1 (x1 ), . . . , ϕk (x1 )) = q(ϕ1 (x2 ), . . . , ϕk (x2 )). Therefore, since g(x1 ) 6= g(x2 ), it follows that there is at least one i ∈ {1, . . . , k} such that ϕi (x1 ) 6= ϕi (x2 ). Thus, ϕi ∈ Aobs (Σ) distinguishes x1 and x2 which is a contradiction. Proof. [Proof of Corollary 4.22] Consider arbitrary x1 6= x2 ∈ X. We denote by ∗ ∗ the dual input-to-state maps corresponding to the Nash system and by τΣ,x τΣ,x 2 1 g Σ = (X, f, h) with the initial state x1 and x2 , respectively. Denote by U pc (x1 , x2 ) the set of all piecewise-constant inputs u such that both xΣ (.; x1 , u) and xΣ (.; x2 , u) are g defined. It is easy to see that U pc (x1 , x2 ) is a set of admissible inputs. Let us assume by g contradiction that h(xΣ (Tu ; x1 , u)) = h(xΣ (Tu ; x2 , u)) for all u ∈ U pc (x1 , x2 ). Then, ∗ ∗ Dαk · · · Dα1 τΣ,x (hi ) = Dαk · · · Dα1 τΣ,x (hi ) 1 2

(A.1)

g for i = 1, . . . , r and for all k ∈ N, α1 , . . . , αk ∈ U . Because U pc (x1 , x2 ) is a set of admissible inputs and because all components hi , i = 1, . . . , r of h are Nash functions on X, we derive by Corollary 6.3 that ∗ ∗ Dαk · · · Dα1 τΣ,x (hi ) = τΣ,x (Lfαk · · · Lfα1 hi ) j j

(A.2)

for j = 1, 2 and for all i = 1, . . . , r, k ∈ N, α1 , . . . , αk ∈ U . From (A.1) and (A.2), ∗ ∗ τΣ,x (Lfαk · · · Lfα1 hi ) = τΣ,x (Lfαk · · · Lfα1 hi ) for i = 1, . . . , r, k ∈ N, α1 , . . . , αk ∈ 1 2 U . Thus, (Lfαk · · · Lfα1 hi )(xΣ (Tu ; x1 , u)) = (Lfαk · · · Lfα1 hi )(xΣ (Tu ; x2 , u)) for all g g u∈U pc (x1 , x2 ) and hence also for the empty input e ∈ U pc (x1 , x2 ). Then, (Lfαk · · · Lfα1 hi )(x1 ) = (Lfαk · · · Lfα1 hi )(x2 )

(A.3)

for all i = 1, . . . , r, k ∈ N, α1 , . . . , αk ∈ U . Because R[{Lfαk · · · Lfα1 hi | i = 1, . . . , r; k ∈ N; α1 , . . . , αk ∈ U }] = Aobs (Σ), (A.3) implies that ∀g ∈ Aobs (Σ) : g(x1 ) = g(x2 ) which contradicts the fact that Σ is strongly semi-algebraically observable.

NASH REALIZATIONS

29

g Therefore, there exists u ∈ U pc (x1 , x2 ) such that h(xΣ (Tu ; x1 , u)) 6= h(xΣ (Tu ; x2 , u)). Proof. [Proposition 4.23] If Σ is semi-algebraically observable, then by Proposition 4.20 the restriction ΣV of Σ to some semi-algebraic subset V of X is strongly semi-algebraically observable. From Corollary 4.22 it follows that ΣV is observable, in particular, it is weakly observable. Then from [7, Theorem 2.1] it follows that ΣV , being an analytic system, is also locally weakly observable. Notice that the observability co-distribution of ΣV is just the restriction of ∆ to V . Hence, according to [15, Theorem 3.11], there exists an open subset W of V such that for all x ∈ W : dim(∆(x)) = dim(V ) = dim(X). If dim(∆(x)) = dim(X) for some x ∈ X, then there exists an open semi-algebraic subset V of X, such that x ∈ V and dim(∆(z)) = dim(X) for all z ∈ V . Hence, there exist ϕ1 , . . . , ϕd ∈ Aobs (Σ), such that for a dϕ1 (z), . . . , dϕd (z) are linearly independent for all z ∈ V . Consider the map g : V 3 z 7→ (ϕ1 (z), . . . , ϕd (z)). By the inverse function theorem, there exists an open semi-algebraic sets U1 ⊆ V and U2 ⊆ Rd such that g : U1 → U2 is a Nash diffeomorphism. Assume now that ϕ1 , . . . , ϕd are not algebraically independent. Then there exists a non-zero polynomial Q such that Q(g(x)) = 0 for all x ∈ U1 . But it then follows that Q(g(g −1 (y)) = Q(y) = 0 for all y ∈ U2 and hence Q is zero on U2 ⊆ Rd . Hence, Q must be a zero polynomial, which is a contradiction. But then it follows that ϕ1 , . . . , ϕd are algebraically independent and hence trdeg Aobs (Σ) = dim(X) = trdeg N (X), i.e. Σ is semi-algebraically observable.

REALIZATION THEORY OF NASH SYSTEMS ...

mean an analytic function satisfying an algebraic equation. A Nash ...... sium on Theoretical Aspects of Computer Software TACS, pp 340–359, 2001. [14] J.J. ...

423KB Sizes 3 Downloads 407 Views

Recommend Documents

Realization Theory of Nash Systems | Google Sites
sary and sufficient conditions for realizability of a response map by a Nash system. ...... hauer. Power-law models of signal transduction pathways. Cell. Signal.,.

Realization Theory of Hybrid Systems
means that as we advance in time, more and more data points are needed to ..... can construct such a hybrid representation from the columns of a suitably big ...

Realization Theory of Bilinear Hybrid Systems
The main tool used in the paper is the theory of for- mal power series. ... Finally, Section 5. develops realization theory for bi- linear hybrid systems. 2. Bilinear ...

Realization theory of linear hybrid systems
International Scientific & Technical Encyclopedia (ISTE), London,. 2006. 43. M. Petreczky and J. H. van Schuppen. Realization theory for linear hy- brid systems ...

Realization Theory of Bilinear Hybrid Systems
Page 1 ... Fliess-series expansion, the rank of its Hankel-matrix is finite and the discrete parts of ... The main tool used in the paper is the theory of for- mal power ...

Realization Theory for Linear Switched Systems
function f is said to be analytic if there exists an open set U ⊆ Rn and a function g : U → Rk ... equations are assumed to have solution on the whole time-axis.

Realization Theory For Linear Hybrid Systems, Part I ...
Zf (s) = ˜CMs ˜Bf = Cqk Aαk+1 qk. Mqk,γk,qk−1 ···Mq1,γ1,q0 Aα1 q0 µ(f). (26) for each f ∈ Φ and j = 1,...,m. If ( ¯A,ζ) is a realization of DΦ, we get that for each f ...

Affine LPV systems: realization theory, input-output ...
q ,Ba q ,Ca q )}D q=1) with na = n. A nonsingular matrix S ∈ Rn×n is said to be an ALPV ... q S = SAq, Ba q = SBq, Ca ..... Notes in Control and Information Sciences. Springer ... IEEE Transactions on Control Systems Technology, 2012. In print.

Realization Theory For Bilinear Switched Systems
Section V presents realization theory of bilinear switched systems. II. BILINEAR SWITCHED SYSTEMS. For sets A, B, denote by PC(A, B) the class of piecewise-.

Realization Theory for Discrete-Time Piecewise-Affine Hybrid Systems
May 21, 2006 - piecewise-affine hybrid system if and only if it has a realization by a ..... DTAPA and we will call Σl the linearised DTAPA associated with Σ. ...... Proceedings of 44th IEEE Conference on Decision and Control, 2005, 2005.

Realization Theory for Linear Switched Systems - Semantic Scholar
the definition of minimality for linear switched systems isn't that obvious. The approach taken in this paper is to treat switched systems as a subclass of ab-.

Realization Theory for Discrete-Time Semi-Algebraic Hybrid Systems
the same output as semi-algebraic systems and implicit polynomial systems. We then derive ...... In: Proc. of 44th IEEE Conference on Decision and Control. (2005) 690–695. 13. ... terdam (2006) available online: http://www.cwi.nl/˜mpetrec. 14.

Realization Theory For Linear Switched Systems
Nov 11, 2006 - In fact, the connection between realization theory of bilinear .... the continuous input u and the switching sequence w = (q1,t1)···(qk,tk). It.

Realization Theory of Discrete-Time Linear Hybrid ...
Keywords: Realisation theory, minimization, identification, hybrid, identifiability ..... (Rank conditions). ... for checking observability and span-reachability and for.

Realization Theory of Nonlinear Hybrid Sys- tems
In the special case of nonlinear analytic systems both the finite Lie-rank and the ... out that existence of a realization by an analytic hybrid system implies that the ...

Realization Theory of Stochastic Jump-Markov Linear ...
JMLSs is the formulation and solution of a stochastic realization problem for a ... In turn, the solution ...... Theoretical Computer Science, 138:101–112, 1995.

Realization theory of discrete-time linear switched ...
ement ai is called the ith letter of w, for i = 1,...,k and k is called the length w. We denote by ...... systems: A tutorial. ... Switched linear systems : control and design.

Realization Theory of Nonlinear Hybrid Sys- tems - Semantic Scholar
system such that the vector fields, reset maps and readout maps are in fact ... hybrid coalgebra realization We will prove that an input-output map cannot have a.

Communication Theory of Secrecy Systems
As in communication theory a language is considered to be ... 1 Shannon, C. E., “A Mathematical Theory of Communication,” Bell System Technical Journal, July. 1948 ... abilities change. If N is large enough (say 50 letters) there is usually a sin

Theory of Computing Systems - Springer Link
Bradley C. Kuszmaul is now Senior Research Scientist at Akamai. Technologies as well as a Research Affiliate at Yale University. Dana S. Henry and Gabriel H.

Nash guide sfv
There was a problem previewing this document. Retrying... Download. Connect more apps... Try one of the apps below to open or edit this item. Nash guide sfv.

The Nash Equilibrium
Bi(a-i) = fai 2 Ai : ui(ai;a-i) ¸ ui(bi;a-i. ) ... Bi(a-i. ) is a set which may or may not be a singleton.) So, for example, in a ..... final order of business in this chapter.

John Nash
May 18, 2006 - For more information regarding JSTOR, please contact support@ j ... solution. Their approach leaves the final situation only determined up.