Realization theory for rational systems: the existence of rational realizations
|
|
- Gertrude Carroll
- 6 years ago
- Views:
Transcription
1 Realization theory for rational systems: the existence of rational realizations Jana Němcová CWI, Amsteram, The Netherlans Jan H. van Schuppen CWI, Amsteram, The Netherlans Abstract In this paper we solve the problem of the existence of rational realizations of response maps. Sufficient an necessary conitions for a response map to be realizable by a rational system are presente. We provie also the characterization of the existence of rationally observable an canonical rational realizations for a given response map. Keywors: Rational systems, realization theory, algebraic reachability, rational observability. Mathematics Subject Classification (2000): 93B15, 93B27, 93C10 Publishe in SIAM Journal on Control an Optimization, 48(4): , DOI: / Introuction The motivation to investigate realization theory of rational systems is the use of rational systems as moels of phenomena in life sciences, in particular in systems biology. For example, rational systems occur as moels of metabolic, genetic, an signaling networks. They can be foun also in engineering, physics, an economics. Moreover, as Bartosiewicz state in [3], the theory of rational systems coul be simpler an more powerful, once it is evelope, than the theory of smooth systems. The realization problem for rational systems consiers a map from input functions to output functions an asks whether there exists a finite-imensional rational system with an initial conition such that its input/output map is ientical to the consiere map. Such a system is then calle a realization of the consiere input/output map. A generalization, which is not aresse in this paper, is to regar any relation between observe variables an ask for a realization as a rational system. Another goal of realization theory is to characterize certain properties of realizations. One wants to fin the conitions uner which the systems realizing the consiere map are observable, controllable, or minimal. The relations between realizations having these properties are also of interest since they can be applie in control an observer synthesis an in system ientification. Polynomial an rational systems are a special class of nonlinear systems amitting a more refine algebraic structure. Realization theory for iscrete-time polynomial systems was formulate by Sontag in [16]. Later, in [18], Wang an Sontag publishe their results on realization theory for polynomial an rational continuous-time systems base on the approach of formal power series in non-commuting variables an on the relation of two characterizations of observation spaces. Another approach to realization theory for polynomial continuous-time systems, motivate by the results of Jakubczyk in [11] for nonlinear realizations, is introuce by Bartosiewicz in [1, 4]. This approach is base on [16]. Furthermore, in [3], Bartosiewicz introuces the concept of rational systems an eals with the problem of immersion of smooth systems into rational systems. Since this problem is similar to the problem of rational realization, there is an analogy between our approach to the realization theory for rational systems an Bartosiewicz s results presente in [3, 4]. Compare to the realization theory for rational systems evelope by Wang an Sontag in [18], our approach is ifferent. We apply the algebraic-geometric approach rather than techniques base on formal power series. We solve the same problem of existence of rational realizations 1
2 REALIZATION THEORY FOR RATIONAL SYSTEMS 2 (compare Theorem 5.2 in [18] an Theorem 5.16 in this paper). In aition we eal with the questions of rational observability an algebraic reachability of rational realizations which are not treate in [18]. Another major ifference is that the realizations within the class of rational systems which we consier o not have to be affine in the inputs as assume by Wang an Sontag. This is motivate by the planne application of realization theory to biochemical systems where the inputs may enter in a rational way. The first step, motivate by biochemical reaction networks, in eveloping a realization theory for rational positive systems is one in [17]. We leave the problem of rational realization with the positivity constraint for further research. The organization of the paper is as follows. Terminology, notation an mathematical preliminaries are provie in Section 2. Section 3 introuces the concept of rational systems an Section 4 introuces the classes of amissible inputs which are consiere omains for response maps. The problem of rational realization is formulate, an necessary an sufficient conitions for the existence of a rational realization of a response map are presente in Section 5. Section 6 eals with canonical rational realizations. Section 7 conclues the paper. 2 Algebraic preliminaries In this section we introuce algebraic-geometric framework in which we formulate the problem of rational realization. Note that the methos of algebraic geometry are alreay known for control an system theory, see for example [6, 8, 9]. For the basic efinitions an theorems of commutative algebra an algebraic geometry see [5, 7, 12, 13, 19]. By a polynomial in finitely many ineterminates X 1,..., X n with real coefficients we mean a sum k N n c(k) n i=1 Xk(i) i where only finitely many coefficients c(k) R are non-zero. We enote the ring of all polynomials in n variables with real coefficients by R[X 1,..., X n ]. Because the fiel R is an integral omain, so is the algebra R[X 1,..., X n ]. Therefore we can efine the fiel of quotients of R[X 1,..., X n ] as the set of fractions {p/q p, q R[X 1,..., X n ], q 0}. This fiel is enote by R(X 1,..., X n ) an we refer to its elements as rational functions. Generally we use the notation Q(S) for the fiel of quotients of an integral omain S. For example, Q(R[X 1,..., X n ]) = R(X 1,..., X n ). By R/I, for a ring R an an ieal I R, we enote the quotient (factor) ring of R moulo I. Polynomial an rational functions on varieties A real affine variety X is an algebraic variety in R n, i.e. there are finitely many f 1,..., f N R[X 1,..., X n ] such that X = {(x 1,..., x n ) R n f 1 (x 1,..., x n ) = = f N (x 1,..., x n ) = 0}. We say that a variety is irreucible if we cannot write it as an union of two non-empty varieties which are its strict subvarieties. Let I R[X 1,..., X n ] be the ieal of all polynomials which are zero at every point of the variety X. By a polynomial on a variety X we mean a map p : X R for which there exists q R[X 1,..., X n ] such that p = q on X. The algebra of all polynomials on X, enote by A, is then isomorphic to R[X 1,..., X n ]/I. From Hilbert Basis Theorem an from the fact that the ieals in the quotient ring R[X 1,..., X n ]/I are in one-to-one corresponence with the ieals of R[X 1,..., X n ] containing I, every ieal in R[X 1,..., X n ]/I is finitely generate. Therefore A is a finitely generate algebra of polynomials. Since X is an irreucible variety, A is an integral omain. Therefore we can efine the fiel Q of quotients of A. The elements of Q are calle rational functions on X. On R n we consier the Zariski topology which is a topology where the close sets are efine as real affine varieties. A variety X R n is enowe with the relate topology, the Zariski topology on X. We refer to an open/close/ense set in Zariski topology as to Z-open/Z-close/Z-ense set. More etails can be foun in [10, 12]. Rational vector fiels Let X be an irreucible real affine variety, let A be the algebra of polynomials on X, an let Q enote the fiel of rational functions on X. Definition 2.1 A rational vector fiel f on X is an R-linear map f : Q Q such that f(ϕ ψ) = f(ϕ) ψ + ϕ f(ψ) for ϕ, ψ Q.
3 REALIZATION THEORY FOR RATIONAL SYSTEMS 3 We say that the rational vector fiel f is efine at the point x X if f(o x ) O x where O x = {ϕ Q ϕ is efine at x}. The set X(f) of all points at which a rational vector fiel f is efine is given as X(f) = {x X f(o x ) O x }. Definition 2.2 The trajectory of a rational vector fiel f from a point x 0 X(f) is the map x : [0, T ) X(f) X such that for t [0, T ) an ϕ A, t (ϕ x)(t) = (fϕ)(x(t)) an x(0) = x 0. Note that it is sufficient to consier only the polynomials ϕ A in the efinition above as it is prove in [3]. Theorem 2.3 [3, Theorem 1] For any rational vector fiel f an any point x 0 X(f) there exists an unique trajectory of f from x 0 efine on the maximal interval [0, T ) (T may be infinite). Proof: The proof of this statement can be foun in [2, 3]. 3 Rational systems We consier rational systems as systems on irreucible real affine varieties with the ynamics efine by rational vector fiels an with output functions having rational components. This concept of rational systems is introuce in [3] an we recall it in this section. We consier an input space U to be an arbitrary set U R m. As an output space we consier R r. By a slight moification of [3, Definition 2] we efine rational systems as follows. Definition 3.1 A rational system Σ with an input space U an an output space R r is a quaruple Σ = (X, f, h, x 0 ) where (i) X is an irreucible real affine variety, (ii) f = {f α α U} is a family of rational vector fiels on X, (iii) h : X R r is an output map with rational components, i.e. h i Q for i = 1,..., r, (iv) x 0 X is an initial state such that all h i, i = 1,..., r an all f α, α U are efine at x 0. The states at which all components of the output function h are efine an at which at least one of the rational vector fiels f α, α U is efine is a Z-ense open subset of X, see [3]. If U is a finite set then the set of points at which all components of h an all vector fiels f α, α U are efine is also a Z-ense open subset of X. As the space of input functions we consier the set U pc of piecewise-constant functions u : [0, ] U. Let u U pc, then u = (α 1, t 1 )... (α nu, t nu ) means that for t ( i j=0 t j, i+1 j=0 t j] the input u(t) = α i+1 U for i = 0, 1,..., n u 1, t 0 = 0, an u(0) = α 1. Every input u U pc has a time omain [0, T u ] where T u = n u j=1 t j. If u = (α 1, t 1 )... (α n, t n ), v = (β 1, s 1 )... (β k, s k ) U pc, then (u)(v) U pc is an input which we get by concatenating v to u, i.e. (u)(v) = (α 1, t 1 )... (α n, t n )(β 1, s 1 )... (β k, s k ). To express that the input u was applie only on time omain [0, t] [0, T u ] we write a subinex [0, t] to u like u [0,t]. The empty input e is such input that T e = 0. Consier a rational system Σ = (X, f, h, x 0 ). The trajectory of the system Σ corresponing to a constant input u = (α, T u ) U pc is the trajectory of the rational vector fiel f α from x 0, i.e. it is the map x( ; x 0, u) : [0, T u ] X for which t (ϕ x)(t; x 0, u) = (f α ϕ)(x(t; x 0, u)) an x(0; x 0, u) = x 0 for t [0, T u ] an for ϕ A. Note that a trajectory of a vector fiel on a close interval is efine as in Definition 2.2, see [4]. The trajectory of the system Σ corresponing to an input u = (α 1, t 1 )... (α nu, t nu ) U pc with T u = n u j=1 t j is the map x( ; x 0, u) : [0, T u ] X such that x(0; x 0, u) = x 0, an x(t; x 0, u) = x αi (t i 1 j=0 t j) for t [ i 1 j=0 t j, i j=0 t j], t 0 = 0, i = 1,..., n u where x αi : [0, t i ] X is a trajectory of a vector fiel f αi from the initial state x( i 1 j=0 t j; x 0, u) = x αi 1 (t i 1 ) for i = 2,..., n u, an from the initial state x 0 for i = 1. Since a trajectory of a rational system Σ oes not nee to exist for every input u U pc, see [3], we efine the set of amissible inputs for Σ.
4 REALIZATION THEORY FOR RATIONAL SYSTEMS 4 Definition 3.2 Let Σ = (X, f, h, x 0 ) be a rational system. We efine the set of amissible inputs U pc (Σ) for the system Σ as a subset of the set of piecewise-constant inputs U pc for which there exist a trajectory of Σ. Note that for every u U pc (Σ) it hols that u [0,t] U pc (Σ) for every t [0, T u ]. We conclue the convention that a trajectory of a rational system Σ with an initial state x 0 X is for the empty input e equal to x 0, i.e. x(0; x 0, e) = x 0. The set of amissible inputs U pc (Σ) may contain only the empty input e. Remark 3.3 Every input of U pc has finitely many switching time-points. This property is neee for the proof of irreucibility of the smallest variety containing a reachable set of a rational system (see Proposition 5.5) which is use in the proof of Theorem 6.1. It is also necessary conition for proving that the set of functions analytic at the switching time-points of such inputs is an integral omain (see Theorem 4.4). This is the crucial property of a class of functions for which there exists a rational realization. 4 Amissible inputs Let U R m be an input space. We efine the sets of amissible inputs for rational systems with the values in U. These sets of inputs are the sets on which the response maps stuie with respect to the realization problem for rational systems are efine. Definition 4.1 A set Ũ pc U pc of input functions with the values in an input space U R m is calle a set of amissible inputs if: (i) u Ũ pc t [0, T u ] : u [0,t] Ũ pc, (ii) u Ũ pc α U t > 0 : (u)(α, t) Ũ pc, (iii) u = (α 1, t 1 )... (α k, t k ) Ũ pc δ > 0 t i [0, t i + δ], i = 1,..., k : u = (α 1, t 1 )... (α k, t k ) Ũ pc. Definition 4.2 Consier a set Ũ pc of amissible inputs with the values in U R m. Let u Ũ pc. We enote the erivation of a real function ϕ : Ũ pc R at the switching time-point T u of the input (u)(α, t) Ũ pc, where t > 0 is sufficiently small an α U, as (D α ϕ)(u) = t ϕ((u)(α, t)) t=0+. Let ϕ : Ũ pc R be a real function an let u = (α 1, t 1 )... (α k, t k ) Ũ pc for a set Ũ pc of amissible inputs. Consier the function ϕ u (t) = ϕ(u [0,t] ) efine for t [0, T u ]. If t [ n i=0 t i, n+1 i=0 t i], n = 0,... k 1 then ϕ u (t) = ϕ((α 1, t 1 )... (α n, t n i=0 t i)). The erivation (D α ϕ)(u) is well-efine if the function ϕ (u)(α,t) (ˆt) = ϕ((u)(α, ˆt)), ˆt [T u, T u + t] is ifferentiable at T u +. We say that the map ϕ : Ũ pc R is smooth if the erivations D α1... D αi ϕ are well-efine on Ũ pc for every i N an α j U, j = 1,..., i. To simplify the notation, the erivation D α1... D αi ϕ can be rewritten as D α ϕ where α = (α 1,..., α i ). Definition 4.3 Consier a set Ũ pc of amissible inputs with the values in U R m. We say that a function ϕ : Ũ pc R is analytic at the switching time-points of the inputs from Ũ pc if for every input u = (u 1, t 1 )... (u k, t k ) Ũ pc the function ϕ u1,...,u k (t 1,..., t k ) = ϕ((u 1, t 1 )... (u k, t k )) is analytic, i.e. we can write ϕ u1,...,u k in the form of convergent formal power series in k ineterminates. We enote the set of real functions ϕ : Ũ pc R which are analytic at the switching time-points of the inputs from Ũ pc by A(Ũ pc R). We refer to the elements of A(Ũ pc R) as to the analytic functions on Ũ pc.
5 REALIZATION THEORY FOR RATIONAL SYSTEMS 5 Let ϕ be an analytic function on Ũ pc. Then for every (u)(α, 0)(v) Ũ pc ϕ((u)(α, 0)(v)) = ϕ((u)(v)). (1) Theorem 4.4 The set A(Ũ pc R) of analytic functions on a set Ũ pc of amissible inputs with the values in U R m is an integral omain. Proof: To prove that A(Ũ pc R) is an integral omain we prove that for f, g A(Ũ pc R) it hols that if fg = 0 on Ũ pc then f = 0 on Ũ pc or g = 0 on Ũ pc. Consier f, g A(Ũ pc R) such that fg = 0. Then fg(u) = 0 for every u Ũ pc. Let u Ũ pc be an arbitrary input. Because f(u), g(u) R an because R is an integral omain, the equality fg(u) = f(u)g(u) = 0 implies that either f(u) = 0 or g(u) = 0. To complete the proof we have to prove that it is not possible that there exist u, v Ũ pc such that f(u) = g(v) = 0 an f(v), g(u) 0. Let us assume for a contraiction that there exist such u = (α 1, t u 1)... (α k, t u k ), v = (β 1, t v 1)... (β l, t v l ) Ũ pc. Since u, v Ũ pc, we erive from Definition 4.1 (ii) that t v i [0, t v i ], i = 1,..., l : w = (α 1, t u 1)... (α k, t u k)(β 1, t v 1 )... (β l, t v l ) Ũ pc. From Definition 4.1 (iii), there exists δ > 0 such that w = (α 1, t u 1 )... (α k, t u k )(β 1, t v 1 )... (β l, t v l ) Ũ pc for every t u i [0, t u i +δ), i = 1,..., k, an for every tv j [0, t v j +δ), j = 1,..., l. Because fg = 0 on Ũ pc, it follows that 0 = f(w )g(w ) an thus 0 = f α1,...,α k,β 1,...,β l (t u 1,..., t u k, tv 1,..., t v k )g α 1,...,α k,β 1,...,β l (t u 1,..., t u k, tv 1,..., t v k ) for t u i [0, t u i + δ), tv j [0, t v j + δ), i = 1,..., k, j = 1,..., l. Hence, f α1,...,α k,β 1,...,β l g α1,...,α k,β 1,...,β l = 0 on [0, t u i + δ) [0, t v i + δ). (2) i=1,...,k j=1,...,l Because f, g A(Ũ pc R) it hols that f α1,...,α k,β 1,...,β l an g α1,...,α k,β 1,...,β l are convergent formal power series in k + l ineterminates for all α 1,..., α k, β 1,..., β l U, an k, l N. Since i=1,...,k [0, tu i + δ) j=1,...,l [0, tv i + δ) is an open connecte set in R k+l, from (2) an from the fact that a ring of convergent formal power series over R in finitely many ineterminates is an integral omain (this follows from [20, Vol.2, Ch.7, Theorem 1]), we erive that either f α1,...,α k,β 1,...,β l = 0 on [0, t u i + δ) [0, t v j + δ), (3) or g α1,...,α k,β 1,...,β l = 0 on i=1,...,k i=1,...,k j=1,...,l [0, t u i + δ) j=1,...,l [0, t v j + δ). (4) By assuming that either (3) or (4) hols we come to a contraiction. This completes the proof. Let us assume that (3) hols. Therefore for τ i [0, t u i + δ), τ j [0, tv j + δ), i = 1,..., k, j = 1,..., l, we get that f α1,...,α k,β 1,...,β l (τ 1,..., τ k, τ 1,..., τ l ) = 0. Because f α 1,...,α k,β 1,...,β l (0,..., 0, τ 1,..., τ l ) = f β 1,...,β l (τ 1,..., τ l ), we erive that f β 1,...,β l (τ 1,..., τ l ) = 0 for τ j [0, tv j +δ), j = 1,..., l. Further, because v Ũ pc, from Definition 4.1 (iii) it follows that ɛ > 0 t j [0, t v j + ɛ), j = 1,..., l : (β 1, t 1 )... (β l, t l ) Ũ pc. Thus, f β1,...,β l can be represente as a convergent formal power series in l ineterminates with a convergence omain containing j=1,...,l [0, tv j + ɛ). Since f β 1,...,β l = 0 on j=1,...,l [0, tv j + δ), then f β1,...,β l = 0 also on j=1,...,l [0, tv j +ɛ). Therefore f β 1,...,β l (t v 1,..., t v l ) = f(v) = 0 which contraicts the assumption f(v) 0. Let us assume that (4) hols. Then g α1,...,α k,β 1,...,β l (τ 1,..., τ k, τ 1,..., τ l ) = 0 for all τ i [0, t u i + δ), τ j [0, tv j + δ), i = 1,..., k, j = 1,..., l, an thus especially for τ i = t u i, i = 1,..., k, an τ j = 0, j = 1,..., l. Therefore 0 = g α1,...,α k,β 1,...,β l (t 1,..., t k, 0,..., 0) = g α1,...,α k (t 1,..., t k ) = g(u) which contraicts the assumption g(u) 0. Corollary 4.5 Because the set A(Ũ pc R) is an integral omain, we can efine the fiel Q(Ũ pc R) of the quotients of elements of A(Ũ pc R).
6 REALIZATION THEORY FOR RATIONAL SYSTEMS 6 5 Rational realizations Response maps In this paper we work with response maps rather than with input/output (I/O) maps since it is technically more convenient. The I/O maps are consiere to be the maps between the spaces of input an output functions (functions of time) mapping an input to an output. We call a map which escribes the outputs immeiately after applying finite parts of the inputs a response map. Definition 5.1 Let Ũ pc be a set of amissible inputs. A map p : Ũ pc R r is calle a response map if its components p i : Ũ pc R, i = 1,..., r are such that p i A(Ũ pc R). We poste extra assumptions on response maps in the efinition above because to solve the problem of realization of a response map by a rational system we use the objects like observation algebra an observation fiel of a response map, see Definition 5.9. These assumptions are necessary for well-efineness of these objects. Problem formulation A rational system which for each input gives us the same output as a response map p is calle a rational realization of p (a rational system realizing p). The realization problem for rational systems can be unerstoo as the problem of fining such a rational system for a given map. Formally we state the problem of realization of a response map by a rational system as follows: Let Ũ pc be a set of amissible inputs. Consier a response map p : Ũ pc R r. The realization problem for rational systems consists of etermining a rational system Σ = (X, f, h, x 0 ) such that p(u) = h(x(t u ; x 0, u)) for all u Ũ pc an Ũ pc U pc (Σ). Properties of rational realizations Note that the solution of a rational system Σ = (X, f, h, x 0 ) corresponing to an input u Ũ pc is a piecewise-analytic map such that the intervals on which the solution is analytic are etermine by the switching time-points of the input u. Definition 5.2 Let Σ = (X, f, h, x 0 ) be a rational realization of a response map p : Ũ pc R r where Ũ pc is a set of amissible inputs. Let A enote the algebra of polynomial functions on X. We efine the input-to-state map τ : Ũ pc X as the map τ(u) = x(t u ; x 0, u) for u Ũ pc. The map τ etermine by τ is efine as τ : A A(Ũ pc R) such that τ (ϕ) = ϕ τ for all ϕ A. To simplify further reference we state some properties of the map τ in Proposition 5.3. The proof of this proposition is omitte because it irectly follows from the efinition of τ. Proposition 5.3 Let Σ = (X, f, h, x 0 ) be a rational realization of a response map p : Ũ pc R r where Ũ pc is a set of amissible inputs. Let A be the algebra of polynomials on X, let ϕ 1,..., ϕ k A, k < be such that A = R[ϕ 1,..., ϕ k ], an let Q enote the fiel of rational functions on X. Then the map τ : A A(Ũ pc R) efine in Definition 5.2 is a homomorphism an τ (R[ϕ 1,..., ϕ k ]) = R[τ ϕ 1,..., τ ϕ k ]. Moreover, the map τ : A/Ker τ R[τ ϕ 1,..., τ ϕ k ], efine as τ ([ϕ]) = τ ϕ for every ϕ A, is an isomorphism. The map τ can be extene to an isomorphism of the fiels Q(A/Ker τ ) an R(τ ϕ 1,..., τ ϕ k ). The efinitions of algebraic reachability (Definition 5.4) an rational observability (Definition 5.6) of rational realizations are base on [3, Definition 3,4]. Definition 5.4 Let Σ = (X, f, h, x 0 ) be a rational realization of a response map p : Ũ pc R r where Ũ pc is a set of amissible inputs. The rational system Σ is sai to be algebraically reachable (from the initial state) x 0 if the reachable set is Z-ense in X. R(x 0 ) = {x(t u ; x 0, u) X u Ũ pc U pc (Σ)}
7 REALIZATION THEORY FOR RATIONAL SYSTEMS 7 Proposition 5.5 Let Σ = (X, f, h, x 0 ) be a rational system as in Definition 5.4. Then the closure Z-cl(R(x 0 )) of the reachable set R(x 0 ) in Zariski topology on X is an irreucible variety. Proof: The Zariski closure of the reachable set R(x 0 ) is the smallest variety in X containing R(x 0 ) which is given as Z-cl(R(x 0 )) = {x X ϕ(x) = 0 for all ϕ A such that ϕ = 0 on R(x 0 )}. Therefore, by consiering τ : Ũ pc X an τ : A A(Ũ pc R) efine in Definition 5.2, we erive that Z-cl(R(x 0 )) = {x X ϕ(x) = 0 for all ϕ A such that ϕ τ = 0 on Ũ pc} = {x X ϕ(x) = 0 for all ϕ A such that τ ϕ = 0} = {x X ϕ(x) = 0 for all ϕ Ker τ }. Because τ is a homomorphism an thus τ (A) is a subalgebra of A(Ũ pc R), an because A(Ũ pc R) is an integral omain, τ (A) is an integral omain. From Proposition 5.3, the map τ : A/Ker τ τ (A) is an isomorphism. This implies, since τ (A) is an integral omain, that Ker τ is a prime ieal. Accoring to [7, Chapter 5.1, Proposition 4], an affine variety V is irreucible if an only if the ieal of polynomials which vanish on V is a prime ieal. Because, the variety Z-cl(R(x 0 )) is etermine by the polynomial ieal Ker τ, an because Ker τ is a prime ieal, the variety Z-cl(R(x 0 )) is irreucible. Definition 5.6 Let Σ = (X, f = {f α α U}, h, x 0 ) be a rational system an let Q enote the fiel of rational functions on X. The observation algebra A obs (Σ) of Σ is the smallest subalgebra of the fiel Q containing all components h i, i = 1,..., r of h, an close with respect to the erivations given by rational vector fiels f α, α U. The observation fiel Q obs (Σ) of the system Σ is the fiel of quotients of A obs (Σ). The rational system Σ is calle rationally observable if Q obs (Σ) = Q. The observation algebra of a rational system is an integral omain because it is a subalgebra of an integral omain. Therefore the observation fiel of a rational system is well-efine. The observation fiel Q obs (Σ) is also close with respect to the erivations given by rational vector fiels f α, α U. Proposition 5.7 ( [3], Proposition 1) For a rational system Σ, Q obs (Σ) is a finitely generate fiel extension of R, i.e. there exist ϕ 1,..., ϕ k Q obs (Σ) such that Q obs (Σ) = R(ϕ 1,..., ϕ k ). Definition 5.8 We call a rational realization of a response map canonical if it is both rationally observable an algebraically reachable. In realization theory one is given a response map an is suppose to fin a system within a certain class of systems which correspons to this map. This map is therefore representing the unknown system. For this reason it is useful to efine objects as observation algebra an observation fiel also for response maps. Definition 5.9 Let Ũ pc be a set of amissible inputs an let p : Ũ pc R r be a response map. The observation algebra A obs (p) of p is the smallest subalgebra of the algebra A(Ũ pc R) which contains the components p i, i = 1,..., r of p, an which is close with respect to the erivations D α, α U. The observation fiel Q obs (p) of p is the fiel of quotients of A obs (p). Note that the observation fiel Q obs (p) of p is well-efine only if A obs (p) is an integral omain. This is the case for response maps because the components of a response map are the elements of A(Ũ pc R) for a set Ũ pc of amissible inputs an because A(Ũ pc R) an thus also Q(Ũ pc R) is an integral omain. For well-efineness of the observation algebra of p it is sufficient to assume that the components of p are smooth (with respect to D α erivations). Existence of rational realizations We provie sufficient an necessary conitions for a response map to be realizable by a rational system. The realizability of response maps by a polynomial system is treate in [4, Theorem 2]. The proof of that theorem an the proof of Proposition 5.14 have the same structure. The following lemma can be foun in [4] state for polynomial systems.
8 REALIZATION THEORY FOR RATIONAL SYSTEMS 8 Lemma 5.10 Let Σ = (X, f, h, x 0 ) be a rational realization of a response map p : Ũ pc R r where Ũ pc is a set of amissible inputs an let τ : Ũ pc X be as in Definition 5.2. Then for any ϕ from the algebra A of polynomials on X an for any α U where U is the set of values of the inputs of Ũ pc it hols that D α (ϕ τ) = (f α ϕ) τ. Proof: Let u Ũ pc an let α U. Because Ũ pc U pc (Σ), the trajectories of the rational system Σ corresponing to the input (u)(α, s) with sufficiently small s > 0 an to all restrictions of (u)(α, s) to shorter time-omains are well-efine. From Definition 4.2 we get that for arbitrary ϕ A D α (ϕ τ)(u) = s (ϕ τ)((u)(α, s)) s=0+ = s ϕ(τ((u)(α, s))) s=0+. By the efinition of an input-to-state map τ, τ((u)(α, s)) = x(t u + s; x 0, (u)(α, s)). accoring to Definition 2.2, Further, D α (ϕ τ)(u) = (f α ϕ)(x(t u + s; x 0, (u)(α, s))) s=0+. Finally, by the continuity of rational function f α ϕ along the trajectory of Σ etermine by the input (u)(α, s) an by the properties of a trajectory, we get that D α (ϕ τ)(u) = (f α ϕ)(τ(u)). Proposition 5.11 Let p : Ũ pc R r be a response map realizable by a rational system Σ = (X, f, h, x 0 ). Let τ : Ũ pc X be as in Definition 5.2. Then the map τ ext : A obs (Σ) A obs (p) efine as τ extϕ = ϕ τ for every ϕ A obs (Σ) is a well-efine surjective homomorphism, i.e. τ ext(a obs (Σ)) = A obs (p). Proof: Note that τext is efine in the same way as τ but on a ifferent omain, see Definition 5.2. Obviously, τext is a homomorphism. We prove that τext is well-efine an that it is surjective. By Definition 5.6, 5.9, the observation algebras of a system Σ an of a map p are generate by h i, f α1... f αj h i an p i, D α1... D αj p i, respectively, such that j N, α 1,..., α j U, an i = 1,..., r. Since τext is a homomorphism, we prove that τext(a obs (Σ)) = A obs (p) by proving that the generators of A obs (Σ) an A obs (p) are mappe to each other by τext. To prove that τext is well-efine it is sufficient to prove that τext is well-efine for the generators of the algebra A obs (Σ). Since Σ is a rational realization of p, we know that p = h τ an that p is well-efine. Because h A obs (Σ), p = τexth which implies that τext is well-efine at h. Let h num, h en A be such that h = hnum h en. For a rational vector fiel f α f, it hols that (f α h) τ = (f α h num h en ) τ = (f αh num τ)(h en τ) (f α h en τ)(h num τ) (h en τ) 2, an further, by Lemma 5.10, that (f α h) τ = D α(h num τ)(h en τ) D α (h en τ)(h num τ) (h en τ) 2. Therefore, τext(f α h) = (f α h) τ = D α ( hnum h en τ) = D α (h τ) = D α (p). As p is a response map, the erivations D α of p are well-efine an consequently τext is well-efine at f α h A obs (Σ) Q. In the next proposition we state necessary conitions for a response map to be realizable by a rational system. Proposition 5.12 Let p : Ũ pc R r be a response map realizable by a rational system Σ = (X, f, h, x 0 ). Let τ ext : A obs (Σ) A obs (p) be as in Proposition Then (i) Q obs (p) = τ (Q(A obs (Σ)/Ker τ ext)), where τ : A obs (Σ)/Ker τ ext A obs (p) is an isomorphism erive from the map τ ext : A obs (Σ) A obs (p), (ii) Q obs (p) is finitely generate.
9 REALIZATION THEORY FOR RATIONAL SYSTEMS 9 Proof: (i) Accoring to Proposition 5.11, the map τext : A obs (Σ) A obs (p) is a surjective homomorphism which is not necessarily injective. Then the map τ : A obs (Σ)/Ker τ ext A obs (p) efine as τ ([ϕ]) = τext(ϕ) for every ϕ A obs (Σ), is an isomorphism. Since the algebras A obs (Σ)/Ker τext an A obs (p) are integral omains, we can construct the fiels of fractions of A obs (Σ)/Ker τext an A obs (p). We exten the isomorphism τ of the algebras A obs (Σ)/Ker τext an A obs (p) to the isomorphism τ of the fiels Q(A obs (Σ)/Ker τext) an Q obs (p). Then Q obs (p) = τ (Q(A obs (Σ)/Ker τext)). (ii) The fiel Q(A obs (Σ)/Ker τext) is a fiel isomorphic to a subfiel F of Q obs (Σ). From [5, Chapter V, Section 14.7, Corollary 3] it follows that if G is a finitely generate fiel containing R, then every subfiel F of G containing R is finitely generate. Then, since F Q obs (Σ) Q an Q is a finite fiel extension of R, it follows that F is also finitely generate. Because a fiel isomorphic to a finitely generate fiel is finitely generate, there exist ϕ 1,..., ϕ k Q(A obs (Σ)/Ker τext) such that Q(A obs (Σ)/Ker τext) = R(ϕ 1,..., ϕ k ). (5) From (i), Q obs (p) = τ (Q(A obs (Σ)/Ker τext)) where τ is an isomorphism. Then, by (5), Q obs (p) = τ (R(ϕ 1,..., ϕ k )) = R( τ ϕ 1,..., τ ϕ k ). Thus Q obs (p) is finitely generate. In the following proposition we prove that the generators of the observation fiel Q obs (p) of a response map p can be chosen from A(Ũ pc R). This allows us to reformulate the necessary conition for a response map p to be realizable by a rational system which is state as the conition (ii) in the proposition above in the following way: Q obs (p) is finitely generate by the elements from A(Ũ pc R). Proposition 5.13 Let p : Ũ pc R r be a response map. The observation fiel Q obs (p) is finitely generate if an only if it is finitely generate by the elements from A(Ũ pc R), i.e. there exist finitely many ϕ 1,..., ϕ k A(Ũ pc R) such that Q obs (p) = R(ϕ 1,..., ϕ k ). Proof: ( ) Let Q obs (p) be finitely generate by the elements from A(Ũ pc R). Then it is obviously finitely generate. ( ) Let Q obs (p) be finitely generate. There exist ϕ 1,..., ϕ k Q obs (p) such that Q obs (p) = R(ϕ 1,..., ϕ k ). As ϕ i Q obs (p), i = 1,..., k we know that ϕ i = ϕi,num ϕ i,en where ϕ i,num, ϕ i,en A obs (p) for i = 1,..., k. We efine the fiel F = R(ϕ 1,num, ϕ 1,en,..., ϕ k,num, ϕ k,en ). Because ϕ i,num, ϕ i,en A obs (p) for i = 1,..., k, an because A obs (p) is a subalgebra of A(Ũ pc R), it follows that ϕ i,num, ϕ i,en A(Ũ pc R) for i = 1,..., k, an (6) R[ϕ 1,num, ϕ 1,en,..., ϕ k,num, ϕ k,en ] A obs (p). (7) By (6) the fiel F is generate by the elements from A(Ũ pc R). From the efinition of F it is obvious that F Q obs (p), an from (7) we get by taking the quotients that F Q obs (p). Therefore Q obs (p) = F an thus the fiel Q obs (p) is finitely generate by the elements from A(Ũ pc R). In the following proposition we specify sufficient conitions for a response map to be realizable by a rational system. Proposition 5.14 Let p : Ũ pc R r be a response map. If there exists a fiel F Q(Ũ pc R) such that (i) F is finitely generate by the elements from A(Ũ pc R), (ii) F is close with respect to D α erivations, i.e. i N, α j U, j = 1,..., i : D α1... D αi F F,
10 REALIZATION THEORY FOR RATIONAL SYSTEMS 10 (iii) Q obs (p) F, then p has a rational realization. Proof: Consier a response map p : Ũ pc R r. We assume that there exists a fiel F Q(Ũ pc R) satisfying (i)-(iii). Let ϕ 1,..., ϕ k : Ũ pc R A(Ũ pc R) be such that F = R(ϕ 1,..., ϕ k ). Since F is close with respect to D α erivations, for any α = (α 1,..., α j ) such that α 1,..., α j U, j N, an for every ϕ i, i = 1,..., k, there exists vi α R(X 1,..., X k ) such that D α ϕ i = vi α (ϕ 1,..., ϕ k ). (8) Since Q obs (p) F, the components of the map p = (p 1,..., p r ) are elements of F an thus for every p j, j = 1,..., r there exists w j R(X 1,..., X k ) such that We prove that a rational system Σ = (X, f, h, x 0 ) where X = R k, k f α = i=1 v α i p j = w j (ϕ 1,..., ϕ k ). (9) x i, α U, h j (x 1,..., x k ) = w j (x 1,..., x k ), j = 1... r, x 0 = (ϕ 1 (e),..., ϕ k (e)), e is the empty input, is a rational system realizing p. Let us efine Ψ(t) = (ϕ 1,..., ϕ k )(u [0,t] ) for u Ũ pc, t [0, T u ]. It is well-efine because ϕ i A(Ũ pc R), i = 1,..., k are efine for every u Ũ pc an because u [0,t] Ũ pc if u Ũ pc an t [0, T u ]. Especially, Ψ(0) = (ϕ 1,..., ϕ k )(u [0,0] ) is well-efine because u [0,0] = e Ũ pc an the functions ϕ 1,..., ϕ k A(Ũ pc R) are efine properly at e. Consier a constant input u = (α, T u ) Ũ pc. Note that from Definition 4.1(i), (ii) there exists ɛ > 0 such that for every t, τ 0 such that t + τ [0, T u + ɛ] it hols that u = (α, t + τ) Ũ pc. If t + τ T u then u = u [0,t+τ], an if T u < t + τ then u = u [0,T u]. In both cases we refer to the corresponing inputs as to the inputs u [0,t+τ], t + τ [0, T u + ɛ]. Then, Ψ(0) = (ϕ 1,..., ϕ k )(u [0,0] ) = (ϕ 1 (e),..., ϕ k (e)) = x 0, an t Ψ(t) = τ Ψ(t + τ) τ=0+ = τ (ϕ 1(u [0,t+τ] ),..., ϕ k (u [0,t+τ] )) τ=0+. Because (u [0,t] )(α, τ) = u [0,t+τ] an because D α ϕ(u [0,t] ) = τ ϕ((u [0,t])(α, τ)) τ=0+, the erivation t Ψ(t) equals (D αϕ 1 (u [0,t] ),..., D α ϕ k (u [0,t] )). Furthermore, ue to (8), t Ψ(t) = (vα 1 (ϕ 1,..., ϕ k )(u [0,t] ),..., v α k (ϕ 1,..., ϕ k )(u [0,t] )). Finally, as Ψ(t) = (ϕ 1,..., ϕ k )(u [0,t] ), we get that t Ψ(t) = (vα 1 (Ψ(t)),..., v α k (Ψ(t))), an Ψ(0) = x 0. Hence Ψ(t) = x(t; x 0, u [0,t] ) for a constant input u Ũ pc an for t [0, T u ]. This also proves that for u = (α 1, t 1 )... (α j, t j ) Ũ pc an for t [0, t 1 ], we have Ψ(t) = x(t; x 0, u [0,t] ). If we consier u = (α 1, t 1 )... (α j, t j ) Ũ pc an if we take t [t 1, t 1 + t 2 ] instea of t [0, t 1 ], we get that Ψ(t 1 ) = x(t 1 ; x 0, u [0,t1]) an tψ(t) = (vα2 1 (Ψ(t)),..., vα2 k (Ψ(t))) for t [t 1, t 1 + t 2 ] with the same reasoning as before. For t [t 1, t 1 + t 2 ] we have that Ψ(t) = x(t; Ψ(t 1 ), u [t1,t 1+t]) = x(t; x(t 1 ; x 0, u [0,t1]), u [t1,t 1+t]). In the analogous way we stuy the cases for t [t 1 + t 2, t 1 + t 2 + t 3 ],..., t [t t i 1, t t i ],.... Finally, Ψ(t) = x(t; x 0, u [0,t] ) for an arbitrary u Ũ pc an t [0, T u ]. Thus the trajectories of Σ are escribe by Ψ. To prove that the rational system Σ is a realization of the response map p, we have to prove that p(u) = h(x(t u ; x 0, u)) for every u Ũ pc. Consier an arbitrary u Ũ pc. Due to (9), p(u) = (p 1,..., p r )(u) = (w 1 (ϕ 1,..., ϕ k ),..., w r (ϕ 1,..., ϕ k ))(u).
11 REALIZATION THEORY FOR RATIONAL SYSTEMS 11 Further, by the efinitions of h j, j = 1,..., r an Ψ it follows that p(u) = (h 1 (ϕ 1,..., ϕ k ),..., h r (ϕ 1,..., ϕ k ))(u) = (h 1 (Ψ(T u )),..., h r (Ψ(T u ))). Finally, since Ψ(T u ) = x(t u ; x 0, u) for u Ũ pc, we erive that p(u) = h(x(t u ; x 0, u)) for u Ũ pc. Remark 5.15 Proposition 5.14 can be state as an equivalence. The proof of the other implication is the same as the sufficiency proof of Theorem The main theorem of this section solving the problem of the existence of rational realizations for response maps is base on the three propositions above. Theorem 5.16 (Existence of rational realizations) A response map p : Ũ pc R r has a rational realization if an only if Q obs (p) is finitely generate. Proof: ( ) See the part (ii) of Proposition 5.12 for this statement an the proof. ( ) From Proposition 5.14, the existence of a fiel F Q(Ũ pc R) finitely generate by the elements from A(Ũ pc R), containing Q obs (p), an close with respect to D α erivations implies the rational realizability of p. Since we assume that Q obs (p) is finitely generate, it follows from Proposition 5.13 an Definition 5.9 that Q obs (p) has those properties. By following the steps of the proof of Proposition 5.14 for F = Q obs (p), we construct a rational realization of p. Example This example is motivate by an example state in [4]. The proceure to construct a rational system realizing a given response map is following the steps mae in the proof of Proposition Let Ũ pc be a set of amissible inputs with the values in R. We etermine a rational system Σ realizing a map p : Ũ pc R efine as p(u) = exp( T u 0 u(s) (1+s) 2 s). By u(s) we enote the value of an input u Ũ pc at a time s [0, T u ]. Firstly we compute D α erivations of p. Consier α 1, α 2 R an u Ũ pc. For t 1, t 2 > 0 sufficiently small (u)(α 1, t 1 )(α 2, t 2 ) Ũ pc. Then (D α1 p)(u) = τ p((u)(α 1, τ)) τ=0+ [ ( Tu = τ exp u(s) (1 + s) 2 s + = 0 + p(u)α 1 [exp 1 = α 1 p(u) (1 + T u ) 2, 0 ( Tu+τ Tu+τ T u α 1 T u α 1 (1 + s) 2 s (1 + s) 2 s )] [ τ=0+ )] τ τ=0+ Tu+τ T u 1 (1 + s) 2 s ] τ=0+ (D α2 D α1 p)(u) = 1 = D α2 (α 1 p(u) (1 + T u ) 2 ) = α 1 1 (1 + T u ) 2 D 1 α 2 p(u) + α 1 p(u)d α2 (1 + T u ) [ ] 2 α 1 α 2 = (1 + T u ) 4 p(u) + α 1 1p(u) τ (1 + T u + τ) 2 = α 1α 2 p(u) (1 + T u ) 4 + 2α 1p(u) (1 + T u ) 3. We can compute the erivations (D αi... D α1 p)(u) for any i N, α j R, j 1,..., i. If we efine ϕ 1 (u) = p(u) an ϕ 2 (u) = 1 + T u then for any i N an α j R, j 1,..., i it hols that (D αi... D α1 p)(u) R(ϕ 1 (u), ϕ 2 (u)). Therefore, by Definition 5.9, Q obs (p) R(ϕ 1, ϕ 2 ) an consequently Q obs (p) is finitely generate. Hence, accoring to Theorem 5.16, there exists a rational system realizing p. We construct a rational system Σ = (X, f, h, x 0 ) realizing p by following the proof of Proposition Let us consier the fiel F = R(ϕ 1, ϕ 2 ). It is finitely generate, contains Q obs (p), τ=0+
12 REALIZATION THEORY FOR RATIONAL SYSTEMS 12 an is close with respect to D α erivations. The number of generators of F equals 2 which implies that the state-space X can be taken as R 2. To etermine a family of rational vector fiels f = {f α α R} we compute v α 1 (ϕ 1, ϕ 2 ) = D α ϕ 1 = D α p = αϕ 1 1 ϕ 2 2 an v α 2 (ϕ 1, ϕ 2 ) = D α ϕ 2 = 1. The last equality D α ϕ 2 = 1 hols because D α ϕ 2 (u) = [ τ ϕ 2((u)(α, τ)) ] τ=0+ = [ τ (1 + T u+ τ)] τ=0+ = [1] τ=0+ = 1 for any u Ũ pc. The output map h is etermine by a map w such that w(ϕ 1, ϕ 2 ) = p = ϕ 1, an the initial point x 0 is given as x 0 = (ϕ 1 (e), ϕ 2 (e)) where e is the empty input. Finally, the rational realization Σ = (X, f, h, x 0 ) of p is given as X = R 2, f α (x 1, x 2 ) = v1 α (x 1, x 2 ) + v2 α (x 1, x 2 ) = α x 1 x 1 x 2 x 2 2 h(x 1, x 2 ) = w(x 1, x 2 ) = x 1, x 0 = (ϕ 1 (e), ϕ 2 (e)) = (1, 1). 6 Canonical rational realizations x 1 + x 2, α R, The systems which we call canonical are calle minimal by Bartosiewicz in [4]. Therefore [4, Theorem 3] escribing the conitions uner which a map has a minimal polynomial realization correspons to our Theorem 6.1 consiering rational systems. The proofs are analogous. Theorem 6.1 Let p : Ũ pc R r be a response map. The following statements are equivalent: (i) p has an rationally observable rational realization, (ii) p has a canonical rational realization, (iii) Q obs (p) is finitely generate. Proof: (i) (ii) Let p : Ũ pc R r be a response map. Consier a rational system Σ = (X, f, h, x 0 ) realizing p which is rationally observable. Then Σ is canonical if it is algebraically reachable from the initial state x 0, i.e. if Z-cl(R(x 0 )) = X. Let us enote X = Z-cl(R(x 0 )). From Proposition 5.5, the variety X is an irreucible real affine variety. If X = X, then the proof is complete. Let us assume that X X. Because X is an irreucible real affine variety, we can consier X to be a state-space of a rational system Σ. Let I be an ieal of polynomials of A which vanish on X. Then the quotient ring A/I can be ientifie with A, the algebra of polynomials on X. We enote the corresponing bijection by Ψ : A/I A. This is an one-to-one an onto mapping which preserves sums an proucts. If we consier ϕ A an ϕ A such that Ψ([ϕ]) = ϕ, it hols that ϕ X = ϕ. More etails can be foun in [7, Chapter 5.2]. There exist polynomials ϕ 1,..., ϕ k from A such that A = R[ϕ 1,..., ϕ k ]. Moreover, A is an integral omain an we can efine the fiel Q of rational functions on X by Q = R(ϕ 1,..., ϕ k ). We efine the rational system Σ = (X, f, h, x 0) by consiering the state-space X an by eriving the family f of rational vector fiels, an output function h an an initial state x 0 from the rational system Σ. As the initial state x 0 of Σ we consier x 0 = x 0 Z-cl(R(x 0 )) = X. The output function h of Σ is efine as ( ) h h = (h 1,..., h 1,num r) =,..., h r,num h 1,en h r,en ( Ψ([h1,num ]) = Ψ([h 1,en ]),..., Ψ([h ) r,num]) Ψ([h r,en ]) = (h 1 X,..., h r X ) = h X, ( h1,num X = h 1,en X,..., h ) r,num X h r,en X where h i,num, h i,en A an h i,num, h i,en A for i = 1,..., r are such that h i = h i,num h an i,en h i = hi,num h i,en. The output function h is efine on a Z-ense subset of X because X is irreucible
13 REALIZATION THEORY FOR RATIONAL SYSTEMS 13 an because h i,en (x 0 ) 0 for i = 1,..., r an thus h i,en / I for i = 1,..., r. We efine the relation between the rational vector fiels f = {f α : Q Q α U} an the rational vector fiels f = {f α : Q Q α U} as f α Ψ([q num ]) Ψ([q en ]) = Ψ([(f αq) num ]) Ψ([(f α q) en ]) for q = qnum q en Q where q num, q en A an q en / I. Note that we assume that (f α q) num, (f α q) en A are such that f α q = (fαq)num (f αq) en for a consiere q Q. Since for ϕ A such that Ψ([ϕ]) = ϕ it is true that ϕ X = ϕ, it follows that for q = qnum Q we have that f α Ψ([q num]) Ψ([q en ]) = (fαq)num X (f αq) en = (f X α q) X, an finally if q = Ψ([qnum]) Ψ([q en ]), q en / I, f αq = (f α q) X. The well-efineness of the rational vector fiels f α, q α U follows from the fact that (f num α q en ) X is inepenent on the choice of representants q num, q en A of classes [q num ], [q en ] A/I. To verify this we consier arbitrary ϕ 1, ϕ 2 I. Thus ϕ 1 an ϕ 2 are ientically zero on X an therefore (q num + ϕ 1 ) X = q num an (q en + ϕ 2 ) X = q en. Moreover as ϕ 1 = 0 on X, ϕ 1 is ientically zero on the trajectories of the system Σ an therefore ϕ 1 τ = 0 on Ũ pc. Then also D α (ϕ 1 τ) = 0 on Ũ pc an by Lemma 5.10, (f α ϕ 1 ) τ = 0 on Ũ pc. Therefore f α ϕ 1 = 0 on X ( an hence )(f α ϕ 1 ) ( X = 0. Analogically, (f α ϕ 2 ) X = 0. Finally, ) as f α q is an R-linear map, f num+ϕ 1 α q en +ϕ 2 X = (fαq num+f αϕ 1)(q en +ϕ 2) (f αq en +f αϕ 2)(q num+ϕ 1) (q en +ϕ 2) 2 X = ( ) ( ) (fαq num)q en (f αq en )q num q X = f num α q en X. q 2 en Consier a trajectory x of Σ from x 0 etermine by the input u = (α, T u ) Ũ pc, an a trajectory x of Σ from x 0 etermine by the input u = (α, T u ) Ũ pc. We assume without loss of generality that T u = T u. Then t (ϕ x )(t; x 0, u ) = (f αϕ )(x (t; x 0, u )), ϕ (x (0)) = x 0 for t [0, T u ], ϕ A, an t (ϕ x)(t; x 0, u) = (f α ϕ)(x(t; x 0, u)), ϕ(x(0)) = x 0 for t [0, T u ], ϕ A. Since X = Z-cl(R(x 0 )), both trajectories x an x stay in X. Let ϕ A an let ϕ A be such that ϕ = Ψ([ϕ]). Because f αψ([ϕ]) = (f α ϕ) X, x 0 = x 0, an because ϕ = ϕ on X, it follows that (f α ϕ) X (x (t; x 0, u )) = (f α ϕ)(x (t; x 0, u )) = t (ϕ x )(t; x 0, u ), an ϕ(x (0)) = x 0. Therefore, by Theorem 2.3, the trajectories of the systems Σ an Σ are the same. Then also the reachable sets of both systems coincie an thus X = Z-cl(R(x 0 )) = Z-cl(R (x 0)). So, Σ is algebraically reachable. From the equality of trajectories of Σ an Σ an from the efinition of the output function h of Σ, it follows that the system Σ is a rational realization of p. More precisely, let u Ũ pc, then p(u) = h(x(t u ; x 0, u)), because Σ realizes p, = (h X )(x(t u ; x 0, u)), because Z-cl(R(x 0 )) = X, = (h X )(x (T u ; x 0, u)), from the equalities of the trajectories of Σ an Σ, = h (x (T u ; x 0, u)), by the efinition of h. Thus, the system Σ is an algebraically reachable rational realization of p. We prove that Σ is also rationally observable. Let us compute the observation fiel Q obs (Σ ) of the system Σ. Firstly, the observation algebra A obs (Σ ) is the smallest algebra containing the elements h i, f αh i for i = 1,..., r an α = (α 1,..., α k ) such that k N, α j U, j = 1,... k. As h i = h i X an f αh i = f α h i X for i = 1,..., k an α = (α 1,..., α k ), α j U, k N, we get that ϕ = ϕnum ϕ en A obs (Σ), where ϕ num, ϕ en A, if an only if ϕ = Ψ([ϕnum]) Ψ([ϕ en ]) A obs (Σ ). Furthermore, since A obs (Σ) an A obs (Σ ) are integral omains, ϕ = ϕnum ϕ en Q obs (Σ), ϕ en / I if an Q obs (Σ ). Because the system Σ is rationally observable an because Q(Ψ(A/I)) = Q(A ), we erive that the system Σ is also rationally observable. (ii) (iii) We assume that p has a canonical rational realization Σ = (X, f, h, x 0 ). From Theorem 5.16, Q obs (p) is finitely generate. (iii) (i) We assume that the observation fiel Q obs (p) of p is finitely generate. From Proposition 5.13, Q obs (p) is finitely generate by the elements of A(Ũ pc R). Let Q obs (p) = R(ϕ 1,..., ϕ k ) where ϕ i A(Ũ pc R), i = 1,..., k. By the efinition of observation fiel, Q obs (p) is close with respect to D α erivations α = (α 1,..., α i ), i N, α j U, j = 1,..., i. To prove that there exists a only if ϕ = Ψ([ϕnum]) Ψ([ϕ en ]) q en
14 REALIZATION THEORY FOR RATIONAL SYSTEMS 14 rationally observable rational realization of p, we construct a rational realization Σ = (X, f, h, x 0 ) of p such that Q obs (Σ) = Q where Q enotes the fiel of rational functions on X. The fiel F = Q obs (p) = R(ϕ 1,..., ϕ k ) fulfills the conitions (i)-(iii) of Proposition By following the proof of Proposition 5.14, with F = Q obs (p), we construct a rational realization Σ = (X, f, h, x 0 ) of p as This realization is such that X = R k, f α = k vi α (X 1,..., X k ), α U, X i i=1 h j (X 1,..., X k ) = w j (X 1,..., X k ), j = 1... r, x 0 = (ϕ 1 (e),..., ϕ k (e)). p j = w j (ϕ 1,..., ϕ k ), j = 1... r, an D α ϕ i = v α i (ϕ 1,..., ϕ k ), i = 1... k. Because X = R k, we get for the fiel Q of rational functions on X that Q = R(X 1,..., X k ). To consier h j an f α h j, j = 1... r is the same as to consier p j an D α p j but in ifferent coorinates. Therefore, Q obs (Σ), as a fiel of quotients of the smallest subalgebra of Q containing all h j, f α h j where j = 1,..., r, α = (α 1,..., α k ), k N, α i U, i = 1,..., k, equals R(X 1,..., X k ) in analogy to the relation Q obs (p) = R(ϕ 1,..., ϕ k ). That means that Q obs (Σ) = R(X 1,..., X k ) = Q which proves the rational observability of Σ. Corollary 6.2 Let p be a response map. Accoring to Theorem 5.16 an Theorem 6.1 the following statements are equivalent: (i) p is realizable by a rational system, (ii) p has a rational realization which is rationally observable, (iii) p has a rational realization which is canonical, thus both rationally observable an algebraically reachable. 7 Concluing remarks We have provie sufficient an necessary conitions for a response map to be realizable by a rational system. Since the proof of Proposition 5.14 concerning sufficient conitions for the existence of a rational realization for a given response map is constructive, it provies an algorithm for constructing a rational system realizing a response map. We prove that for a given response map the problems of the existence of a rational/rationally observable rational/canonical rational realization are equivalent. The proof of Theorem 6.1 stating this equivalence is also constructive. It provies an algorithm for the construction of a canonical rational realization from a rationally observable rational realization. Moreover it shows that if we consier as a fiel F from Proposition 5.14 the observation fiel of a response map to be realize, then the algorithm for constructing a rational realization from Proposition 5.14 gives as a result a rationally observable rational realization. Aitional results which we got by stuying minimal rational realization can be foun in [15]. For some remarks on algebraic reachability of rational systems see [14]. We work with irreucible real affine varieties for two reasons. Firstly, irreucibility simplifies the technical etails of the proofs. Seconly, working with the real varieties allows us to have a better geometric unerstaning of the state-spaces of rational systems an it is also sufficient for applications for example in systems biology. It is possible to generalize the results of this paper for reucible varieties. Algebraic framework we use can be useful from computational point of view. The proceures to check the properties of rational systems such as rational observability an algebraic reachability, an minimality of rational realizations are formulate in [15]. They can be implemente by using alreay existing computer algebra packages. The proceures for constructing rational realizations of esire properties are still to be evelope. Further research is require on the realization theory for rational positive systems. There are several issues to be overcome in that research.
15 REALIZATION THEORY FOR RATIONAL SYSTEMS 15 Acknowlegements The authors are very grateful to Mihály Petreczky for long iscussions, useful comments an avice. References [1] Z. Bartosiewicz. Realizations of polynomial systems. In M. Fliess an M. Hazenwinkel, eitors, Algebraic an Geometric Methos in Nonlinear Control Theory, pages D. Reiel Publishing Company, Dorrecht, [2] Z. Bartosiewicz. Orinary ifferential equations on real affine varieties. Bulletin of the Polish Acaemy of Sciences Mathematics, 35(1-2):13 18, [3] Z. Bartosiewicz. Rational systems an observation fiels. Systems an Control Letters, 9: , [4] Z. Bartosiewicz. Minimal polynomial realizations. Mathematics of control, signals, an systems, 1: , [5] N. Bourbaki. Algebra II, Chapters 4-7. Elements of mathematics. Springer-Verlag, [6] C.I. Byrnes an P.L. Falb. Applications of algebraic geometry in system theory. American Journal oj Mathematics, 101(2): , [7] D. Cox, J. Little, an D. O Shea. Ieals, varieties, an algorithms: An introuction to computational algebraic geometry an commutative algebra. Springer, 3r eition, [8] P. Falb. Methos of algebraic geometry in control theory: Part 1, Scalar linear systems an affine algebraic geometry. Birkhäuser, Boston, [9] P. Falb. Methos of algebraic geometry in control theory: Part 2, Multivariable linear systems an projective algebraic geometry. Birkhäuser, Boston, [10] N. Jacobson. Basic Algebra II. W. H. Freeman an company, San Francisco, [11] B. Jakubczyk. Existence an uniqueness of realizations of nonlinear systems. SIAM J. Control an Optimization, 18(4): , July [12] E. Kunz. Introuction to commutative algebra an algebraic geometry. Birkhäuser, Boston, [13] S. Lang. Algebra. Aison-Wesley, Reaing, Massachusetts, [14] J. Němcová. Algebraic reachability of rational systems. In Proceeings of European Control Conference, Buapest, Hungary, [15] J. Němcová an J.H. van Schuppen. Realization theory for rational systems: Minimal rational realizations. to appear in Acta Applicanae Mathematicae. [16] E.D. Sontag. Polynomial response maps. Number 13 in Lecture Notes in Control an Information Sciences. Springer-Verlag, Berlin, Heielberg, [17] J.H. van Schuppen. System theory of rational positive systems for cell reaction networks. In Proceeings of MTNS, [18] Y. Wang an E.D. Sontag. Algebraic ifferential equations an rational control systems. SIAM J. Control Optim., 30(5): , [19] D.J. Winter. The structure of fiels. Springer-Verlag, New York, [20] O. Zariski an P. Samuel. Commutative algebra I, II. Springer, 1958.
6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationWitt#5: Around the integrality criterion 9.93 [version 1.1 (21 April 2013), not completed, not proofread]
Witt vectors. Part 1 Michiel Hazewinkel Sienotes by Darij Grinberg Witt#5: Aroun the integrality criterion 9.93 [version 1.1 21 April 2013, not complete, not proofrea In [1, section 9.93, Hazewinkel states
More information7. Localization. (d 1, m 1 ) (d 2, m 2 ) d 3 D : d 3 d 1 m 2 = d 3 d 2 m 1. (ii) If (d 1, m 1 ) (d 1, m 1 ) and (d 2, m 2 ) (d 2, m 2 ) then
7. Localization To prove Theorem 6.1 it becomes necessary to be able to a enominators to rings (an to moules), even when the rings have zero-ivisors. It is a tool use all the time in commutative algebra,
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationInternational Journal of Pure and Applied Mathematics Volume 35 No , ON PYTHAGOREAN QUADRUPLES Edray Goins 1, Alain Togbé 2
International Journal of Pure an Applie Mathematics Volume 35 No. 3 007, 365-374 ON PYTHAGOREAN QUADRUPLES Eray Goins 1, Alain Togbé 1 Department of Mathematics Purue University 150 North University Street,
More informationZachary Scherr Math 503 HW 3 Due Friday, Feb 12
Zachary Scherr Math 503 HW 3 Due Friay, Feb 1 1 Reaing 1. Rea sections 7.5, 7.6, 8.1 of Dummit an Foote Problems 1. DF 7.5. Solution: This problem is trivial knowing how to work with universal properties.
More informationTHE GENUINE OMEGA-REGULAR UNITARY DUAL OF THE METAPLECTIC GROUP
THE GENUINE OMEGA-REGULAR UNITARY DUAL OF THE METAPLECTIC GROUP ALESSANDRA PANTANO, ANNEGRET PAUL, AND SUSANA A. SALAMANCA-RIBA Abstract. We classify all genuine unitary representations of the metaplectic
More informationIntegration Review. May 11, 2013
Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In
More informationSystems & Control Letters
Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationThe canonical controllers and regular interconnection
Systems & Control Letters ( www.elsevier.com/locate/sysconle The canonical controllers an regular interconnection A.A. Julius a,, J.C. Willems b, M.N. Belur c, H.L. Trentelman a Department of Applie Mathematics,
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationDIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10
DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 5. Levi-Civita connection From now on we are intereste in connections on the tangent bunle T X of a Riemanninam manifol (X, g). Out main result will be a construction
More informationCHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold
CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationarxiv: v1 [math.dg] 1 Nov 2015
DARBOUX-WEINSTEIN THEOREM FOR LOCALLY CONFORMALLY SYMPLECTIC MANIFOLDS arxiv:1511.00227v1 [math.dg] 1 Nov 2015 ALEXANDRA OTIMAN AND MIRON STANCIU Abstract. A locally conformally symplectic (LCS) form is
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationDiscrete Mathematics
Discrete Mathematics 309 (009) 86 869 Contents lists available at ScienceDirect Discrete Mathematics journal homepage: wwwelseviercom/locate/isc Profile vectors in the lattice of subspaces Dániel Gerbner
More informationDiophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations
Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationON BEAUVILLE STRUCTURES FOR PSL(2, q)
ON BEAUVILLE STRUCTURES FOR PSL(, q) SHELLY GARION Abstract. We characterize Beauville surfaces of unmixe type with group either PSL(, p e ) or PGL(, p e ), thus extening previous results of Bauer, Catanese
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationNEERAJ KUMAR AND K. SENTHIL KUMAR. 1. Introduction. Theorem 1 motivate us to ask the following:
NOTE ON VANISHING POWER SUMS OF ROOTS OF UNITY NEERAJ KUMAR AND K. SENTHIL KUMAR Abstract. For xe positive integers m an l, we give a complete list of integers n for which their exist mth complex roots
More informationA nonlinear inverse problem of the Korteweg-de Vries equation
Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeare in a journal publishe by Elsevier. The attache copy is furnishe to the author for internal non-commercial research an eucation use, incluing for instruction at the authors institution
More informationSINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
More informationDECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS
DECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS ARNAUD BODIN Abstract. We state a kin of Eucliian ivision theorem: given a polynomial P (x) an a ivisor of the egree of P, there exist polynomials h(x),
More informationAlmost Split Morphisms, Preprojective Algebras and Multiplication Maps of Maximal Rank
Syracuse University SURFACE Mathematics Faculty Scholarship Mathematics 12-30-2005 Almost Split Morphisms, Preprojective Algebras an Multiplication Maps of Maximal Rank Steven P. Diaz Syracuse University
More informationStable Polynomials over Finite Fields
Rev. Mat. Iberoam., 1 14 c European Mathematical Society Stable Polynomials over Finite Fiels Domingo Gómez-Pérez, Alejanro P. Nicolás, Alina Ostafe an Daniel Saornil Abstract. We use the theory of resultants
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationLecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012
CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationZachary Scherr Math 503 HW 5 Due Friday, Feb 26
Zachary Scherr Math 503 HW 5 Due Friay, Feb 26 1 Reaing 1. Rea Chapter 9 of Dummit an Foote 2 Problems 1. 9.1.13 Solution: We alreay know that if R is any commutative ring, then R[x]/(x r = R for any r
More informationACTIONS IN MODIFIED CATEGORIES OF INTEREST WITH APPLICATION TO CROSSED MODULES Dedicated to Teimuraz Pirashvili on his 60th birthday
Theory an Applications of Categories, Vol. 30, No. 25, 2015, pp. 882 908. ACTIONS IN MODIFIED CATEGORIES OF INTEREST WITH APPLICATION TO CROSSED MODULES Deicate to Teimuraz Pirashvili on his 60th birthay
More informationDarboux s theorem and symplectic geometry
Darboux s theorem an symplectic geometry Liang, Feng May 9, 2014 Abstract Symplectic geometry is a very important branch of ifferential geometry, it is a special case of poisson geometry, an coul also
More informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
More informationLecture 10: October 30, 2017
Information an Coing Theory Autumn 2017 Lecturer: Mahur Tulsiani Lecture 10: October 30, 2017 1 I-Projections an applications In this lecture, we will talk more about fining the istribution in a set Π
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationPermanent vs. Determinant
Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More informationu!i = a T u = 0. Then S satisfies
Deterministic Conitions for Subspace Ientifiability from Incomplete Sampling Daniel L Pimentel-Alarcón, Nigel Boston, Robert D Nowak University of Wisconsin-Maison Abstract Consier an r-imensional subspace
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationTractability results for weighted Banach spaces of smooth functions
Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More information. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp
. ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationarxiv: v1 [physics.class-ph] 20 Dec 2017
arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.
More informationMath 1B, lecture 8: Integration by parts
Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS
ON THE RIEANN EXTENSION OF THE SCHWARZSCHILD ETRICS Valerii Dryuma arxiv:gr-qc/040415v1 30 Apr 004 Institute of athematics an Informatics, AS R, 5 Acaemiei Street, 08 Chisinau, olova, e-mail: valery@ryuma.com;
More information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationInterconnected Systems of Fliess Operators
Interconnecte Systems of Fliess Operators W. Steven Gray Yaqin Li Department of Electrical an Computer Engineering Ol Dominion University Norfolk, Virginia 23529 USA Abstract Given two analytic nonlinear
More informationTensors, Fields Pt. 1 and the Lie Bracket Pt. 1
Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 PHYS 500 - Southern Illinois University September 8, 2016 PHYS 500 - Southern Illinois University Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 September 8,
More informationA new proof of the sharpness of the phase transition for Bernoulli percolation on Z d
A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition
More informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.
More informationALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS
ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an
More informationREVERSIBILITY FOR DIFFUSIONS VIA QUASI-INVARIANCE. 1. Introduction We look at the problem of reversibility for operators of the form
REVERSIBILITY FOR DIFFUSIONS VIA QUASI-INVARIANCE OMAR RIVASPLATA, JAN RYCHTÁŘ, AND BYRON SCHMULAND Abstract. Why is the rift coefficient b associate with a reversible iffusion on R given by a graient?
More informationHilbert functions and Betti numbers of reverse lexicographic ideals in the exterior algebra
Turk J Math 36 (2012), 366 375. c TÜBİTAK oi:10.3906/mat-1102-21 Hilbert functions an Betti numbers of reverse lexicographic ieals in the exterior algebra Marilena Crupi, Carmela Ferró Abstract Let K be
More informationSYMPLECTIC GEOMETRY: LECTURE 3
SYMPLECTIC GEOMETRY: LECTURE 3 LIAT KESSLER 1. Local forms Vector fiels an the Lie erivative. A vector fiel on a manifol M is a smooth assignment of a vector tangent to M at each point. We think of M as
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationNotes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk
Notes on Lie Groups, Lie algebras, an the Exponentiation Map Mitchell Faulk 1. Preliminaries. In these notes, we concern ourselves with special objects calle matrix Lie groups an their corresponing Lie
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationGeneralized Nonhomogeneous Abstract Degenerate Cauchy Problem
Applie Mathematical Sciences, Vol. 7, 213, no. 49, 2441-2453 HIKARI Lt, www.m-hikari.com Generalize Nonhomogeneous Abstract Degenerate Cauchy Problem Susilo Hariyanto Department of Mathematics Gajah Maa
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationAcute sets in Euclidean spaces
Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationCharacterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations
Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationPerturbation Analysis and Optimization of Stochastic Flow Networks
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationFractional Geometric Calculus: Toward A Unified Mathematical Language for Physics and Engineering
Fractional Geometric Calculus: Towar A Unifie Mathematical Language for Physics an Engineering Xiong Wang Center of Chaos an Complex Network, Department of Electronic Engineering, City University of Hong
More informationOn Decentralized Optimal Control and Information Structures
2008 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June 11-13, 2008 FrC053 On Decentralize Optimal Control an Information Structures Naer Motee 1, Ali Jababaie 1 an Bassam
More informationJournal of Algebra. A class of projectively full ideals in two-dimensional Muhly local domains
Journal of Algebra 32 2009 903 9 Contents lists available at ScienceDirect Journal of Algebra wwwelseviercom/locate/jalgebra A class of projectively full ieals in two-imensional Muhly local omains aymon
More informationUNDERSTANDING INTEGRATION
UNDERSTANDING INTEGRATION Dear Reaer The concept of Integration, mathematically speaking, is the "Inverse" of the concept of result, the integration of, woul give us back the function f(). This, in a way,
More informationRealization theory for systems biology
Realization theory for systems biology Mihály Petreczky CNRS Ecole Central Lille, France February 3, 2016 Outline Control theory and its role in biology. Realization problem: an abstract formulation. Realization
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationA Constructive Inversion Framework for Twisted Convolution
A Constructive Inversion Framework for Twiste Convolution Yonina C. Elar, Ewa Matusiak, Tobias Werther June 30, 2006 Subject Classification: 44A35, 15A30, 42C15 Key Wors: Twiste convolution, Wiener s Lemma,
More informationLIAISON OF MONOMIAL IDEALS
LIAISON OF MONOMIAL IDEALS CRAIG HUNEKE AND BERND ULRICH Abstract. We give a simple algorithm to ecie whether a monomial ieal of nite colength in a polynomial ring is licci, i.e., in the linkage class
More informationON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS
ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS ADRIAN CONSTANTIN AND BORIS KOLEV Abstract. Accoring to the principle of least action, the spatially perioic motions of one-imensional
More informationDiscrete Hamilton Jacobi Theory and Discrete Optimal Control
49th IEEE Conference on Decision an Control December 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Discrete Hamilton Jacobi Theory an Discrete Optimal Control Tomoi Ohsawa, Anthony M. Bloch, an Melvin
More informationAlgebra IV. Contents. Alexei Skorobogatov. December 12, 2017
Algebra IV Alexei Skorobogatov December 12, 2017 Abstract This course is an introuction to homological algebra an group cohomology. Contents 1 Moules over a ring 2 1.1 Definitions an examples.........................
More informationApproximate reduction of dynamic systems
Systems & Control Letters 57 2008 538 545 www.elsevier.com/locate/sysconle Approximate reuction of ynamic systems Paulo Tabuaa a,, Aaron D. Ames b, Agung Julius c, George J. Pappas c a Department of Electrical
More informationTOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH
English NUMERICAL MATHEMATICS Vol14, No1 Series A Journal of Chinese Universities Feb 2005 TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH He Ming( Λ) Michael K Ng(Ξ ) Abstract We
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationSelection principles and the Minimal Tower problem
Note i Matematica 22, n. 2, 2003, 53 81. Selection principles an the Minimal Tower problem Boaz Tsaban i Department of Mathematics an Computer Science, Bar-Ilan University, Ramat-Gan 52900, Israel tsaban@macs.biu.ac.il,
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More information