1. Introduction. The problem of transforming the nonlinear control singleinput. ξ = f(ξ)+g(ξ)u. x = φ(ξ), u = α(ξ)+β(ξ)v

Size: px
Start display at page:

Download "1. Introduction. The problem of transforming the nonlinear control singleinput. ξ = f(ξ)+g(ξ)u. x = φ(ξ), u = α(ξ)+β(ξ)v"

Transcription

1 SIAM J CONTROL OPTIM Vol 0, No 0, pp c 2003 Society for Industrial and Applied Mathematics FEEDBACK CLASSIFICATION OF NONLINEAR SINGLE-INPUT CONTROL SYSTEMS WITH CONTROLLABLE LINEARIZATION: NORMAL FORMS, CANONICAL FORMS, AND INVARIANTS ISSA AMADOU TALL AND WITOLD RESPONDEK Abstract We study the feedback group action on single-input nonlinear control systems We follow an approach of Kang and Krener based on analyzing, step by step, the action of homogeneous transformations on the homogeneous part of the same degree of the system We construct a dual normal form and dual invariants with respect to those obtained by Kang We also propose a canonical form and a dual canonical form and show that two systems are equivalent via a formal feedback if and only if their canonical forms resp, their dual canonical forms coincide We give an explicit construction of transformations bringing the system to its normal, dual normal, canonical, and dual canonical forms We illustrate our results by simple examples on R 3 and R 4 Key words feedback equivalence, normal forms, canonical forms, invariants AMS subect classifications 93B11, 93B17, 93B27 PII S Introduction The problem of transforming the nonlinear control singleinput system Σ : ξ = fξ+gξu by a feedback transformation of the form Γ : x = φξ, u = αξ+βξv to a simpler form has been extensively studied during the last twenty years transformation Γ brings Σ to the system The Σ : ẋ = fx+ gxv, whose dynamics are given by f = φ f + gα, g = φ gβ, where for any vector field f and any diffeomorphism φ we denote φ fx =dφφ 1 x fφ 1 x A natural question to ask is whether we can find a transformation Γ such that the Received by the editors November 28, 2000; accepted for publication in revised form April 21, 2002; published electronically DATE Institut National des Sciences Appliquées de Rouen, Laboratoire de Mathématiques de l INSA, Pl Emile Blondel, Mont Saint Aignan Cedex, France tall@insa-rouenfr, wresp@insarouenfr The first author worked on this paper while on leave from the Department of Mathematics, Université Cheikh Anta Diop, Dakar, Senegal 1

2 2 ISSA AMADOU TALL AND WITOLD RESPONDEK transformed system Σ is linear, that is, whether we can linearize the system Σ via feedback Necessary and sufficient geometric conditions for this to be the case have been given in [13] and [18] Those conditions, except for the planar case, turn out to be restrictive, and a natural problem that arises is to find normal forms for nonlinearizable systems Although natural, this problem is very involved and has been extensively studied during the last twenty years Four basic methods have been proposed for studying feedback equivalence problems The first method is based on the theory of singularities of vector fields and distributions, and their invariants, and using this method a large variety of feedback classification problems have been solved; see, eg, [4], [7], [14], [15], [18], [19], [27], [29], [32], [38] The second approach, proposed by Gardner [9], uses Cartan s method of equivalence [6] and describes the geometry of feedback equivalence [10], [11], [12], [28] The third method, inspired by the Hamiltonian formalism for optimal control problems, was developed by Bonnard [3], [4] and Jakubczyk [16], [17] and has led to a very nice description of feedback invariants in terms of singular extremals Finally, a very fruitful approach was proposed by Kang and Krener [26] and then followed by Kang [21], [22] Their idea, which is closely related with Poincaré s classical technique for linearization of dynamical systems see, eg, [1], is to analyze the system Σ and the feedback transformation Γ step by step and, as a consequence, to produce a simpler equivalent system Σ also step by step Our paper is deeply inspired by those of Kang and Krener [26], [21] and can be considered as a completion of their results In [21], Kang constructed a normal form for single-input nonlinear control systems with controllable linearization using successively homogeneous feedback transformations, and he proved that the homogeneous terms of a given degree of his normal form are unique under homogeneous feedback transformations of the same degree He also showed that a nonlinear system can admit different normal forms under feedback resulting from the action of lower order terms of the feedback transformation on higher order terms of the system The main goal of our paper is to propose a canonical form for the class of single-input systems with controllable linearization and to prove that two systems are equivalent, via a formal feedback, if and only if their canonical forms coincide In [26] Kang and Krener constructed two normal forms for the quadratic part of a single-input system In the first normal form, all components of the linear part of the control vector field are annihilated and all nonremovable quadratic nonlinearities are grouped in the drift; in the second normal form, all quadratic terms of the drift are annihilated and all nonremovable nonlinearities are present in the control vector field Kang normal form is a generalization, for higher order terms, of the first normal form In this paper, we generalize the second one and produce a dual normal form for higher order terms We also construct dual invariants of homogeneous feedback transformations They contain the same information, as Kang invariants, encoded in a different way We also give a dual canonical form and prove that two systems are equivalent, via formal feedback, if and only if their dual canonical forms coincide The third aim of the paper is to construct explicit homogeneous feedback transformations which bring the homogeneous part of the system of the same degree into its normal, or dual normal, form For any fixed degree, our transformations are easily computable via differentiation and integration of polynomials A successive application of those transformations gives formal feedbacks that bring any system to its normal form, dual normal form, canonical form, and dual canonical form The theory of normal forms initialized and developed by Kang and Krener [26] and Kang [21], [22] and continued in the present paper and in [33], [34] has proved to be very useful in analyzing structural properties of nonlinear control systems It

3 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 3 has been used to study bifurcations of nonlinear systems [23], [24], [25], has led to a complete description of symmetries around equilibrium [30], [31], and has allowed us to characterize systems equivalent to feedforward forms [35], [36], [37] The paper is organized as follows In section 2 we will introduce, following [21] and [26], homogeneous feedback transformations We give a normal form obtained by Kang and discuss invariants of homogeneous transformations, also obtained by him We provide an explicit construction of transformations bringing the system to Kang normal form In section 3 we construct a canonical form and give one of our main results stating that two control systems are feedback equivalent if and only if their canonical forms coincide Proofs of results presented in sections 2 and 3 are given in section 4 Section 5 dualizes the main results of section 2: we give a dual normal form, explicitly construct transformations bringing the system to that form, and define dual invariants of homogeneous transformations Similarly to normal forms, a given system can admit different dual normal forms In section 6 we thus dualize the results of section 3 by constructing a dual canonical form and proving that two control systems are feedback equivalent if and only if their dual canonical forms coincide Section 7 contains proofs of results presented in sections 5 and 6 Throughout the paper, we illustrate our results by simple examples on R 3 and R 4 2 Normal form and m-invariants All obects, that is, functions, maps, vector fields, control systems, etc, are considered in a neighborhood of 0 R n and assumed to be C -smooth Let h be a smooth R-valued function By hx =h [0] x+h [1] x+h [2] x+ = h x we denote its Taylor series expansion at 0 R n, where h x stands for a homogeneous polynomial of degree m Similarly, for a map φ of an open subset of R n to R n resp, for a vector field f on an open subset of R n we will denote by φ resp, by f the homogeneous term of degree m of its Taylor series expansion at 0 R n, that is, each component φ of φ resp, f of f is a homogeneous polynomial of degree m in x We will denote by H x the space of homogeneous polynomials of degree m of the variables x 1,,x n and by H m x the space of formal power series of the variables x 1,,x n starting from terms of degree m Analogously, we will denote by R x the space of homogeneous vector fields whose components are in H x and by R m x the space of vector fields formal power series whose components are in H m x Consider the Taylor series expansion of the system Σ given by m=0 21 Σ : ξ = Fξ + Gu + m=2 f ξ+g [m 1] ξu, where F = f ξ 0 and G = g0 We will assume throughout the paper that f0=0 and g0 0 Consider also the Taylor series expansion Γ of the feedback transformation Γ

4 4 ISSA AMADOU TALL AND WITOLD RESPONDEK given by 22 Γ : x = Tξ + φ ξ, m=2 u = Kξ + Lv + m=2 α ξ+β [m 1] ξv, where T is an invertible matrix and L 0 Let us analyze the action of Γ on the system Σ step by step To start with, consider the linear system ξ = Fξ + Gu Throughout the paper we will assume that it is controllable It can be thus transformed by a linear feedback transformation of the form Γ 1 : x = Tξ, u = Kξ + Lv to the Brunovský canonical form A, B; see, eg, [20] Assuming that the linear part F, G, of the system Σ given by 21, has been transformed to the Brunovský canonical form A, B, we follow an idea of Kang and Krener [26], [21] and apply successively a series of transformations 23 Γ m : x = ξ + φ ξ, u = v + α ξ+β [m 1] ξv for m =2, 3, A feedback transformation defined as a series of successive compositions of Γ m, m =1, 2,, will also be denoted by Γ because, as a formal power series, it is of the form 22 We will not address the problem of convergence and will call such a series of successive compositions a formal feedback transformation Observe that each transformation Γ m for m 2 leaves invariant all homogeneous terms of degree smaller than m of the system Σ, and we will call Γ m a homogeneous feedback transformation of degree m We will study the action of Γ m on the following homogeneous system: 24 Σ : ξ = Aξ + Bu + f ξ+g [m 1] ξu Consider another homogeneous system, Σ, given by 25 Σ : ẋ = Ax + Bv + f x+ g [m 1] xv We will say that the homogeneous system Σ is feedback equivalent to the homogeneous system Σ if there exists a homogeneous feedback transformation of the form 23, which brings Σ into Σ modulo terms in R m+1 x, v

5 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 5 Notation Because of various normal forms and various transformations that are used throughout the paper, we will keep the following notation We will denote, respectively, by Σ and Σ the following systems: Σ : ξ = Aξ + Bu + f ξ+g [m 1] ξu, Σ : ξ = Aξ + Bu + f [k] ξ+g [k 1] ξu k=2 The systems Σ and Σ will stand for the systems under consideration Their state vector will be denoted by ξ and their control by u The system Σ resp, the system Σ transformed via feedback will be denoted by Σ resp, by Σ Its state vector will be denoted by x, its control by v, and the vector fields, defining its dynamics, by f [k] and g [k 1] Feedback equivalence of homogeneous systems Σ and Σ will be established via a smooth feedback, that is, precisely, via a homogeneous feedback Γ m On the other hand, feedback equivalence of systems Σ and Σ will be established via a formal feedback Γ We will introduce two kinds of normal forms, Kang normal forms and dual normal forms, as well as canonical forms and dual canonical forms The bar symbol will correspond to the vector field f defining the Kang normal forms Σ NF and Σ NF and the canonical form Σ CF as well as to the vector field ḡ[m 1] defining the dual normal forms Σ DNF and Σ DNF and the dual canonical form Σ DCF Analogously, the m-invariants resp, dual m-invariants of the system Σ will be denoted by a,i+2 resp, by b [m 1] and the m-invariants resp, dual m-invariants of the normal form Σ NF resp, dual normal form Σ DNF byā,i+2 [m 1] resp, by b The starting point is the following result, proved by Kang [21] Proposition 1 The homogeneous feedback transformation Γ m, defined by 23, brings the system Σ, given by 24, into Σ, given by 25, if and only if the relations L Aξ φ φ +1 ξ = f ξ f ξ, 26 L B φ n L Aξ φ n L B φ + α ξ = ξ = g [m 1] f n ξ+β [m 1] ξ = g [m 1] n ξ g [m 1] ξ, ξ f n ξ, ξ g n [m 1] ξ hold for any 1 n 1, where φ are the components of φ This proposition represents the essence of the method developed by Kang and Krener and used in our paper The problem of studying the feedback equivalence of two systems Σ and Σ requires, in general, solving a system of first order partial differential equations On the other hand, if we perform the analysis step by step, then the problem of establishing the feedback equivalence of two systems Σ and Σ reduces to solving the algebraic system 26 Using the above proposition, Kang [21] proved the following result Theorem 1 The homogeneous system Σ can be transformed, via a homogeneous feedback transformation Γ m, into the normal form

6 6 ISSA AMADOU TALL AND WITOLD RESPONDEK 27 Σ NF : ẋ 1 = x 2 + n x 2 i P [m 2] 1,i x 1,,x i, i=3 n ẋ = x +1 + x 2 i P [m 2],i x 1,,x i, i=+2 ẋ n 2 = x n 1 + x 2 np [m 2] n 2,n x 1,,x n, ẋ n 1 = x n, ẋ n = v, where P [m 2],i x 1,,x i are homogeneous polynomials of degree m 2 depending on the indicated variables In order to construct invariants of homogeneous feedback transformations, let us define X m 1 i ξ = 1 i ad i Aξ+f ξ B + g[m 1] ξ and let X [m 1] i be its homogeneous part of degree m 1 By π i we will denote the proection on the subspace 28 W i = {ξ R n : ξ i+1 = = ξ n =0}, that is, π i ξ =ξ 1,,ξ i, 0,,0 Following Kang [21], we denote by a,i+2 ξ the homogeneous part of degree m 2 of the polynomials CA 1 [ X m 1 i,x m 1 ] i+1 πn i ξ, where C =1, 0,,0 T R n and, i N N, defined by ={, i N N : 1 n 2, 0 i n 2} The homogeneous polynomials a,i+2 for, i will be called m-invariants of Σ The following result of Kang [21] asserts that m-invariants a,i+2 for, i are complete invariants of homogeneous feedback and, moreover, illustrates their meaning for the homogeneous normal form Σ NF Consider two homogeneous systems Σ and Σ and let and { a,i+2 :, i } { ã,i+2 :, i } denote, respectively, their m-invariants The following theorem was proved by Kang [21] Theorem 2 The m-invariants have the following properties:

7 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 7 i Two homogeneous systems Σ and Σ are equivalent via a homogeneous feedback transformation Γ m if and only if a,i+2 =ã,i+2 for any, i ii The m-invariants ā,i+2 of the homogeneous normal form Σ NF, defined by 27, are given by 29 ā,i+2 x = 2 x 2 n i x 2 n ip [m 2],n i x 1,,x n i for any, i Our first aim is to find explicitly feedback transformations bringing the homogeneous system Σ to its normal form Σ NF Define the homogeneous polynomials ψ [m 1],i ξ by setting ψ [m 1],0 ξ =ψ [m 1] 1,1 ξ =0, n i 210 ψ [m 1],i ξ = CA 1 ad n i Aξ g[m 1] + 1 t ad t 1 Aξ ad A n i t Bf if 1 <i n and t=1 211 ψ [m 1],i ξ =L A n i Bf 1 π iξ + L Aξ ψ [m 1] 1,i π iξ + ψ [m 1] 1,i 1 π i 1ξ + ξi 0 L A n i+1 Bψ [m 1] 1,i π iξdξ i if 1 i, where ψ [m 1],i π i ξ is the restriction of ψ [m 1],i ξ to the submanifold W i Define the components φ of φ for 1 n and the feedback α,β [m 1] by 212 φ ξ = φ n i=1 ξi 0 ξ =f n 1 ξ+l Aξφ α ξ = β [m 1] ξ = f n g [m 1] n ψ [m 1],i π i ξdξ i, 1 n 1, n 1 ξ, ξ+l Aξ φ n ξ ξ+l B φ n ξ, We have the following result Theorem 3 The homogeneous feedback transformation Γ m : x = ξ + φ ξ, u = v + α ξ+β [m 1] ξv, where α, β [m 1], and the components φ of φ are defined by 212, brings the homogeneous system Σ into its normal form Σ NF given by 27 Proof of Theorem 3 Denote by Σ : ẋ = Ax + Bv + f x+ g [m 1] xv

8 8 ISSA AMADOU TALL AND WITOLD RESPONDEK the system Σ transformed via the feedback transformation Γ m defined by 212 From the expressions of α ξ and β [m 1] ξ given by 212 and the last two equations of 26, we get f n x = 0 and g [m 1] n x =0 Plugging φ, defined by 212, into the second equation of 26 gives ψ [m 1],n x = g [m 1] x g [m 1] x, which, by 210, implies g [m 1] x = 0 for 1 n 1 Now we consider the first equation of 26 From the expression of φ n we get f n 1 x = 0, and for any 1 i n, we obtain by differentiating 213 f x i = f x i φ + L Ax x i + φ +1 φ x i 1 x i In the above formula, the term φ x i 1 is not present in the case i =1 If i +1,weget f x i π i 1 x = f x i = f x i φ + L Ax x i Hence, by an induction argument, we obtain f x i + φ +1 φ x i 1 x i π i 1 x + L Ax ψ [m 1],i π i 1 x + ψ [m 1],i 1 π i 1x ψ [m 1] +1,i π i 1x π i 1 x π i 1 x + L Ax ψ [m 1],i π i 1 x + ψ [m 1],i 1 π i 1x ψ [m 1] +1,i π i 1x=0 and, finally, we get 214 f π i 1 x=0 x i If 1 i, then, using 212 and 213, we obtain f x i π i x = f x i = f + x i xi 0 φ + L Ax x i + φ +1 φ x i 1 x i π i x π i x + L Ax ψ [m 1],i π i x + ψ [m 1],i 1 π i 1x ψ [m 1],i x i 1 π i x dx i ψ [m 1] +1,i π ix

9 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 9 Using the expression 211, it follows that 215 f π i x=0 x i From the relations 214 and 215, we conclude that f x = x 2 i P [m 2],i x 1,,x i, i=+2 which proves that Σ is a normal form satisfying 27 Thus the system Σ given by 24 is feedback equivalent to the normal form Σ NF given by 27 Example 1 To illustrate results of this section, we consider the system Σ, given by 24 on R 3 Theorem 1 implies that the system Σ is equivalent, via a homogeneous feedback transformation Γ m defined by 212, to its normal form Σ NF see 27 ẋ 1 = x 2 + x 2 3P [m 2] x 1,x 2,x 3, ẋ 2 = x 3, ẋ 3 = v, where P [m 2] x 1,x 2,x 3 is a homogeneous polynomial of degree m 2 of the variables x 1,x 2,x 3 We would like now to discuss the interest of Theorem 3 As we have already mentioned, Poincaré s method allows us to replace a partial differential equation by solving successively linear algebraic equations defined by the homological equation 26; see [26] and [21], and Proposition 1 The solvability of this equation was proved in [26] and [21], while Theorem 3 provides an explicit solution in the form of the transformations 212, which are easily computable via differentiation and integration of homogeneous polynomials to the homological equation As a consequence, for any given control system, Theorem 3 gives transformations bringing the homogeneous part of the system to its normal form For example, if the system is feedback linearizable, up to order m 0 1 see [27], then a diffeomorphism and a feedback compensating all nonlinearities of degree lower than m 0 can be calculated explicitly without solving partial differential equations More generally, by a successive application of transformations given by 212 we can bring the system, without solving partial differential equations, to its normal form given in Theorem 4 below Consider the system Σ of the form 21 and recall that we assume the linear part F, G to be controllable Apply successively to Σ a series of transformations Γ m, m =1, 2,, such that each Γ m brings Σ to its normal form Σ NF ; for instance we can take a series of transformations defined by 212 Successive repeating of Theorem 1 gives the following result of Kang [21] Theorem 4 There exists a formal feedback transformation Γ which brings the

10 10 ISSA AMADOU TALL AND WITOLD RESPONDEK system Σ to a normal form Σ NF given by ẋ 1 = x 2 + x 2 i P 1,i x 1,,x i, 216 Σ NF : i=3 ẋ = x +1 + x 2 i P,i x 1,,x i, i=+2 ẋ n 2 = x n 1 + x 2 np n 2,n x 1,,x n, ẋ n 1 = x n, ẋ n = v, where P,i x 1,,x i are formal power series depending on the indicated variables Example 2 Consider a system Σ defined on R 3 whose linear part is controllable Theorem 4 implies that the system Σ is equivalent, via a formal feedback transformation Γ, to its normal form Σ NF ẋ 1 = x 2 + x 2 3P x 1,x 2,x 3, ẋ 2 = x 3, ẋ 3 = v, where P x 1,x 2,x 3 is a formal power series of the variables x 1,x 2,x 3 3 Canonical form As proved by Kang and recalled in Theorem 2, the normal form Σ NF is unique under homogeneous feedback transformation Γm The normal form Σ NF is constructed by a successive application of homogeneous transformations Γ m for m 1 which bring the corresponding homogeneous systems Σ into their normal forms Σ NF Therefore a natural and fundamental question which arises is whether the system Σ can admit two different normal forms, that is, whether the normal forms given by Theorem 4 are in fact canonical forms under a general formal feedback transformations of the form Γ It turns out that a given system can admit different normal forms, as the following example of Kang [21] shows The main reason for the nonuniqueness of the normal form Σ NF is that, although the normal form Σ NF is unique, homogeneous feedback transformation Γ m bringing Σ into Σ NF is not It is this small group of homogeneous feedback transformations of order m that preserve Σ NF described by Proposition 2 below, which causes the nonuniqueness of Σ NF The aim of this section is thus to construct a canonical form for Σ under feedback transformation Γ Example 3 Consider the system 31 ξ 1 = ξ 2 + ξ3 2 2ξ 1 ξ3, 2 ξ 2 = ξ 3, ξ 3 = u on R 3 Clearly, this system is in Kang normal form compare with Theorem 4 The

11 feedback transformation Γ 3 : FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 11 x 1 = ξ 1 ξ ξ3 2, x 2 = ξ 2 2ξ 1 ξ 2, x 3 = ξ 3 2ξ2 2 + ξ 1 ξ 3 2ξ 2 ξ3, 2 u = v +6ξ 2 ξ 3 +12ξ 1 ξ 2 ξ 3 4ξ3 3 +2ξ 1 +2ξ1 2 +2ξ 2 ξ 3 v brings the system 31 into the form ẋ 1 = x 2 + x 2 3, ẋ 2 = x 3, ẋ 3 = v modulo terms in R 4 x, v Applying successively homogeneous feedback transformations Γ m given, for any m 4, by 212, we transform the above system into the normal form 32 ẋ 1 = x 2 + x x 2 3P x, ẋ 2 = x 3, ẋ 3 = v, where P is a formal power series whose 1-et at 0 R 3 vanishes The systems 31 and 32 are in their normal forms and, moreover, feedback equivalent, but the latter system does not contain any term of degree 3 As a consequence, the normal form Σ NF is not unique under formal feedback transformations Consider the system Σ of the form 21 Since its linear part F, G is assumed to be controllable, we bring it, via a linear transformation and linear feedback, to the Brunovský canonical form A, B Let the first homogeneous term of Σ which cannot be annihilated by a feedback transformation be of degree m 0 As proved by Krener [27], the degree m 0 is given by the largest integer p such that all distributions D k = span {g,,ad k 1 f g} for 1 k n 1 are involutive modulo terms of order p 2 We can thus, due to Theorems 1 and 2, assume that, after applying a suitable feedback, Σ takes the form ξ = Aξ + Bu + f [m0] ξ+ m=m 0+1 f ξ+g [m 1] ξu, where A, B is in Brunovský canonical form and the first nonvanishing homogeneous vector field f [m0] is of the form ξ f [m0] i 2 P [m0 2],i ξ 1,,ξ i, 1 n 2, ξ = i=+2 0, n 1 n Let i 1,,i n s, where i i n s = m 0 and i n s 2, be the largest, in the lexicographic ordering, n s-tuple of nonnegative integers such that for some 1 n 2, we have m0 [m0] f 0 ξ i1 1 ξin s

12 12 ISSA AMADOU TALL AND WITOLD RESPONDEK Define { = sup =1,,n 2 : m0 [m0] f ξ i1 1 ξin s } 0 We have the following result Theorem 5 The system Σ given by 21 is equivalent by a formal feedback Γ to a system of the form 33 Σ CF : ẋ = Ax + Bv + m=m 0 f x, where, for any m m 0, 34 f x = additionally, we have 35 i=+2 and, moreover, for any m m 0 +1, x 2 i P [m 2],i x 1,,x i, 1 n 2, 0, n 1 n; m0 [m0] f x i1 1 xin s = ±1 36 m0 f x 1, 0,,0=0 x i1 1 xin s The form Σ CF satisfying 34, 35, and 36 will be called the canonical form of Σ The name is ustified by the following Theorem 6 Two systems Σ 1 and Σ 2 are formally feedback equivalent if and only if their canonical forms Σ 1,CF and Σ 2,CF coincide Proofs of Theorems 5 and 6 are given in section 4 Kang [21], generalizing [26], proved that any system Σ can be brought by a formal feedback into the normal form 33 for which 34 is satisfied He also observed that his normal forms are not unique; see Example 3 Our results, Theorems 5 and 6, complete his study We show that for each degree m of homogeneity we can use a onedimensional subgroup of feedback transformations which preserves the triangular structure of 34 and at the same time allows us to normalize one higher order term The form of 35 and 36 is a result of this normalization These one-dimensional subgroups of feedback transformations are given by the following proposition Proposition 2 The transformation Γ m given by 23 leaves invariant the system Σ defined by 24 if and only if 37 φ = a m L 1 Aξ ξm 1, 1 n, α = a m L n Aξξ m 1, β [m 1] = a m L B L n 1 Aξ ξm 1,

13 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 13 where a m is an arbitrary real parameter Proof of Proposition 2 Observe that, following Proposition 1, the transformation Γ m leaves invariant the system Σ if and only if it satisfies the following system of equations: L Aξ φ φ +1 ξ =0, 1 n 1, L B φ =0, 1 n 1, L Aξ φ n + α ξ =0, L B φ n + β [m 1] ξ =0 In order to solve the above system, let us remark, using the second equation of the system, that for any such that 1 n 1, the component φ does not depend to the variable ξ n Putting = n 2 into the first equation, we get φ n 2 ξ 1 ξ φ n 2 ξ n 1 ξ n = φ n 1 Since φ n 1 and φ n 2 do not depend on the variable ξ n, we conclude that φ n 2 does not depend on the variable ξ n 1 An inductive argument shows that φ 1 depends only on the variable ξ 1, that is, φ 1 ξ =a m ξ1 m Now, all equations of 37 follow easily Theorem 5 establishes an equivalence of the system Σ with its canonical form Σ CF via a formal feedback Its direct corollary yields the following result for equivalence under a smooth feedback of the form Γ : x = φξ, u = αξ+βξv, up to an arbitrary order Corollary 1 Consider a smooth control system Σ : ξ = fξ+gξu For any positive integer k we have the following: i There exists a smooth feedback Γ transforming Σ, locally around 0 R n, into its canonical form Σ k CF given by Σ k CF : ẋ = Ax + Bv + k m=m 0 f x, modulo Ox, v k+1, where f x, for any m 0 m k, satisfies 34, 35, 36 ii Feedback equivalence of Σ and Σ k CF, modulo Ox, vk+1, can be established via a polynomial feedback transformation Γ k of degree k iii Two smooth systems Σ 1 and Σ 2 are feedback equivalent modulo terms of order Ox, v k+1 if and only if their canonical forms Σ k 1,CF and Σ k 2,CF coincide This corollary follows directly from Theorem 5 and its proof, given in section 4, which provides explicit polynomial transformations bringing, step by step, the system into its canonical form

14 14 ISSA AMADOU TALL AND WITOLD RESPONDEK We will illustrate results of this section by two examples Example 4 Let us reconsider the system Σ given by Example 2 It is feedback equivalent to the normal form ẋ 1 = x 2 + x 2 3P x 1,x 2,x 3, ẋ 2 = x 3, ẋ 3 = v, where P x 1,x 2,x 3 is a formal power series Assume, for simplicity, that m 0 =2, which is equivalent to the following generic condition: g, ad f g, and [g, ad f g] are linearly independent at 0 R 3 This implies that we can express P = P x 1,x 2,x 3 as P = c + P 1 x 1 +x 2 P 2 x 1,x 2 +x 3 P 3 x 1,x 2,x 3, where c 0 and P 1 0 = 0 Observe that any P x 1,x 2,x 3, of the above form, gives a normal form Σ NF In order to get the canonical form Σ CF, we use Theorem 5, which ensures the existence of a feedback transformation Γ of the form x = φx, v = αx+βxṽ, which normalizes the constant c and annihilates the formal power series P 1 x 1 More precisely, Γ transforms Σ into its canonical form Σ CF, x 1 = x 2 + x 2 P 3 x 1, x 2, x 3, x 2 = x 3, x 3 =ṽ, where the formal power series P x 1, x 2, x 3 is of the form P x 1, x 2, x 3 =1+ x 2 P2 x 1, x 2 + x 3 P3 x 1, x 2, x 3 Now, we give an example of constructing the canonical form for a physical model of variable length pendulum Example 5 Consider the variable length pendulum of Bressan and Rampazzo [5] see also [2] and [8] We denote by ξ 1 the length of the pendulum, by ξ 2 its velocity, by ξ 3 the angle with respect to the horizontal, and by ξ 4 the angular velocity The control u = ξ 4 = ξ 3 is the angular acceleration The mass is normalized to 1 The equations are compare [5] and [8] ξ 1 = ξ 2, ξ 2 = g sin ξ 3 + ξ 1 ξ4, 2 ξ 3 = ξ 4, ξ 4 = u, where g denotes the gravity Notice that if we suppose to control the angular velocity ξ 4 = ξ 3, which is the case of [5] and [8], then the system is three-dimensional but the control enters nonlinearly At any equilibrium point ξ 0 =ξ 10,ξ 20,ξ 30,ξ 40 T =ξ 10, 0, 0, 0 T, the linear part of the system is controllable Our goal is to produce, for the variable length pendulum,

15 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 15 a normal form and the canonical form as well as to answer the question whether the systems corresponding to various values of the gravity constant g are feedback equivalent In order to get a normal form, put The system becomes x 1 = ξ 1, x 2 = ξ 2, x 3 = g sin ξ 3, x 4 = gξ 4 cos ξ 3, v = gξ4 2 sin ξ 3 ug cos ξ 3 ẋ 1 = x 2, x 1 ẋ 2 = x 3 + x 2 4 g 2 x 2, 3 ẋ 3 = x 4, ẋ 4 = v, which gives a normal form Indeed, we rediscover Σ NF, given by 216, with P 1,3 =0, P 1,4 = 0, and P 2,4 = x 1 g 2 x 2 3 In order to bring the system to its canonical form Σ CF, first observe that m 0 =3 Indeed, the function x 2 x 1 4 starts with third order terms, which corresponds to the g 2 x 2 3 fact that the invariants a [2],i+2 vanish for any 1 2 and any 0 i 2 The only nonzero component of f [3] is f [3] 2 = x 2 4P [1] 2,4 Hence = 2 and the only, and thus largest, quadruplet i 1,i 2,i 3,i 4 of nonnegative integers, satisfying i 1 +i 2 +i 3 +i 4 =3 and such that 3 f [3] 2 x i1 1 xi4 4 0, is i 1,i 2,i 3,i 4 =1, 0, 0, 2 In order to normalize f [3] 2, put x i = a 1 x i, 1 i 4, ṽ = a 1 v, where a 1 =1/g We get the following canonical form for the variable length pendulum: x 1 = x 2, x 2 = x 3 + x 2 x x 2, 3 x 3 = x 4, x 4 =ṽ Independently of the value of the gravity constant g, all systems are feedback equivalent to each other

16 16 ISSA AMADOU TALL AND WITOLD RESPONDEK 4 Proofs of Theorems 5 and 6 Proof of Theorem 5 The proof of this theorem will be done in two steps In the first step we will deal with terms of degree m 0 Then we will prove the general step by an induction argument First step Let us consider the system Σ given by 21 and let m 0 be the degree of the first nonlinearizable homogeneous part We can assume that see Theorems 1 and 2, after applying a suitable feedback transformation, the system Σ given by 21 takes the form 41 ξ = Aξ + Bu + f [m0] ξ+ m=m 0+1 f ξ+g [m 1] ξu, where A, B is in Brunovský canonical form and the first nonvanishing vector field f [m0] is of the form f [m0] ξ = i=+2 Notice that the linear feedback transformation ξ 2 i P [m0 2],i ξ 1,,ξ i, 1 n 2, 0, n 1 n Γ 1 : x = a 1 ξ, u = 1 a 1 v, where a 1 R and a 1 0, brings the system 41 into the following one: ẋ = Ax + Bv + 1 a m0 1 1 f [m0] x+ m=m 0+1 By the definitions of i 1,,i n s and,wehave m0 [m0] f 0, x i1 1 xin s and thus we can suitably choose the parameter a 1 such that f x+ g [m 1] xv m0 [m0] f = ±1 x i1 1 xin s General step Now, we assume that, for some l 1, the system Σ given by 21, takes the form 42 m 0+l 1 Σ : ξ = Aξ + Bu + m=m 0 f ξ+f [m0+l] ξ+g [m0+l 1] ξu + rξ,u, where rξ,u R m0+l+1 ξ,u and the vector fields f ξ for any m such that m 0 m m 0 + l 1 satisfy the conditions 34, 35, and 36 We will construct a transformation Γ which preserves all terms of degree smaller than m 0 + l while taking those of degree m 0 + l into the canonical form defined by 34 and 36

17 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 17 Consider the following feedback transformation 43 Γ : x = ξ + u = v + m=l+1 m=l+1 φ ξ, α ξ+β [m 1] ξv, where, for any m such that m 0 m m 0 + l 1, the triplet φ,α,β [m 1] is given by 37 and φ =0,α = 0, and β [m 1] = 0 for m m 0 + l +1 The transformation Γ is actually a polynomial transformation Γ m0+l and can be viewed as a composition of a transformation Γ m0+l 1 and a homogeneous transformation Γ m0+l defined, respectively, by 44 and 45 Γ m0+l 1 : Γ m0+l : m 0+l 1 y = ξ + φ ξ, u = w + m=l+1 m 0+l 1 m=l+1 x = y + φ [m0+l] y, α ξ+β [m 1] ξw w = v + α [m0+l] y+β [m0+l 1] yv Let us denote by Σ the system Σ, given by 42, transformed via Γ m0+l 1 Since f [m0] ξ = f [m0] y φ [l+1] y = f [m0] [m0] f y φ [l+1] y+r 1 y, y where r 1 y R m0+l+1 y and for any m m 0 +1, f ξ = f y φ [l+1] y = f y+r 2 y, where r 2 y R m0+l+1 y, we get 46 Σ : m 0+l 1 ẏ = Ay + Bw + f y+ f [m0+l] y+ g [m0+l 1] yw + r 3 y, w, m=m 0 where r 3 y, w R m0+l+1 y, w and f [m0+l] = f [m0+l] +[ f [m0],φ [l+1] ], g [m0+l 1] = g [m0+l 1] Let { a [m0+l],i+2 } :, i

18 18 ISSA AMADOU TALL AND WITOLD RESPONDEK and { ã [m0+l],i+2 } :, i denote, respectively, the sets of m 0 + l-invariants associated with the homogeneous parts of degree m 0 + l of the systems 42 and 46 We have 47 where ã [m0+l],i+2 = a [m0+l],i+2 +â [m0+l],i+2, [ i [ â [m0+l],i+2 = CA 1 1 i+k ad A i Bad i k Aξ ad A k B f [m 0],φ [l+1]] π n i ξ k=0 i 1 ] [ + 1 i+k ad A i+1 Bad i k 1 ad A k B f [m 0],φ [l+1]] π n i ξ k=0 Since the identity ad k A n 1 B adi Aξh = ad i Aξad k A n 1 B h holds for any vector field h and any k, i 0, we get by differentiating 48 L i1+l An 1Bâ[m0+l],i+2 [ i [ = CA 1 1 i+k ad A i Bad i k Aξ ad A Bad i1+l k A n 1 B f [m 0],φ [l+1]] π n i ξ k=0 i 1 ] [ + 1 i+k ad A i+1 Bad i k 1 Aξ ad A Bad i1+l k A n 1 B f [m 0],φ [l+1]] π n i ξ k=0 Due to the definition of the n s-tuple i 1,,i n s, we obtain [ ad i1+l A n 1 B f [m 0],φ [l+1]] = c 1 [ad i1a f n 1B [m0],ad l A n 1 B φ[l+1]] 49 [ + c 2 ad i1 1 f [m0] A n 1 B,ad l+1 A n 1 B φ[l+1]], where c 1 and c 2 are strictly positive integers From the relations ad l A n 1 B φ[l+1] = a l+1 l + 1!ξ 1,ξ 2,,ξ n T, ad l+1 A n 1 B φ[l+1] = a l+1 l + 1!1, 0,,0 T, one can easily deduce that identity 410 can be rewritten as [ ad i1+l A n 1 B f [m 0],φ [l+1]] = γ l ad i1 f [m0] A n 1 B, where we set γ l = a l+1 l + 1! c 1 m 0 i 1 +1+c 2 Plugging the above identity into the formula 48, we obtain L i1+l An 1Bâ[m0+l],i+2 [ i = γ l CA 1 1 i+k ad i1 A n 1 B ad A i Bad i k Aξ ad f [m0] A k B π n i ξ k=0 i 1 ] + 1 i+k ad i1 A n 1 B ad A i+1 Bad i k 1 ad f [m0] A k B π n i ξ k=0 Aξ Aξ

19 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 19 Since f [m0] is of the form 34, we get for any k such that 0 k i 1, ad Ak B f [m0] π n i ξ=0 and for any t 0, ad t Aξad A k B f [m0] π n i ξ=0 Thus, we can deduce the relation i1+l â [m0+l],i+2 = γ l CA ξ i1+l 1 1 i1+2 f [m0] ξ i1 1 ξ2 n i π n i ξ, which leads, after differentiating and setting = and i = s, to the following one: m0+l 2 â [m0+l],s+2 ξ i1+l 1 ξ i2 2 ξin s 2 m0 [m0] f = γ l ξ i1 1 ξi2 2 ξin s Differentiating 47 and using the above identity, we get m0+l 2 ã [m0+l],s+2 ξ i1+l 1 ξ i2 2 ξin s 2 m0 [m0] f + γ l 2 ξin s 2 ξ i1 1 ξi2 2 ξin s = m0+l 2 a [m0+l],s+2 ξ i1+l 1 ξ i2 n s We can choose suitably the parameter a l+1 recall the definition of γ l such that we obtain m0+l 2 ã [m0+l],s+2 =0 ξ i1+l 1 ξ i2 2 ξin i 2 Now, transforming the homogeneous part of degree m 0 + l of the system 46 to its normal form via a homogeneous transformation Γ m0+l and taking into account Theorem 2, we bring the system 46 into the form Σ : ẋ = Ax + Bv + m 0+l m=m 0 f x+rx, v, where rx, v R m0+l+1 x, v and the vector fields f, for any m such that m 0 m m 0 + l, satisfy the conditions 34, 35, and 36 This completes the proof of Theorem 5 In our proof of Theorem 6, we will use the following result Lemma 1 A transformation Γ leaves invariant all terms of degree smaller than m 0 + l of the system 42 if and only if Γ is of the form 410 Γ : x = ξ + u = v + m=l+1 m=l+1 φ ξ, α ξ+β [m 1] ξv, where, for any m such that m 0 m m 0 + l 1, the triplet φ,α,β [m 1] is given by 37

20 20 ISSA AMADOU TALL AND WITOLD RESPONDEK Proof of Lemma 1 We have shown, when proving Theorem 5, that the transformation Γ, defined by 410 and 37, leaves invariant all terms of degree smaller than m 0 + l of the system 42 Conversely, assume that the transformation Γ leaves invariant all terms of degree smaller than m 0 + l of the system 42 Without loss of generality, we can take Γ : x = ξ + φ ξ, u = v + m=k+1 m=k+1 α ξ+β [m 1] ξv, where k+1 denotes the smallest degree among degrees of all nonvanishing components φ of the transformation Γ There is nothing to prove if k +1 m 0 + l Wethus focus our attention on the case k +2 m 0 + l Since Γ leaves invariant all terms of degree smaller than m 0 + l of the system 42, in particular it leaves invariant terms of degree k +1, which implies that φ [k+1],α [k+1],β [k] satisfies the condition 37 By induction, we show that φ,α,β [m 1] also satisfies the condition 37 for any m such that k +1 m m 0 + k 1 Thus it remains only to prove that k l Assume this is false; that is, suppose k l 1 We can see that the transformation Γ brings the system 42 into the following one: 411 m 0+k 1 ẋ = Ax + Bv + f x+ f [m0+k] x+rx, v, m=m 0 where rx, v R m0+k+1 x, v and the vector field f x, for any m such that m 0 m m 0 + k 1, is of the form 34 and 36 and f [m0+k] = f [m0+k] +[ f [m0],φ [k+1] ] Since the transformation Γ leaves invariant all terms of degree smaller than m 0 + l of the system 42, in particular it leaves invariant all terms of degree m 0 + k, which is equivalent to [ f [m0],φ [k+1] ]=0 Repeating the calculations done in the proof of Theorem 5 we deduce, by differentiating, the identity m0+k CA 1 [ f [m0],φ [k+1] ] x i1+k 1 x i2 2 xin s Thus, due to the fact that m0 [m0] f x i1 1 xi2 2 xin s m0 [m0] f = γ k x i1 1 xi2 2 xin s we obtain γ k = 0 and hence φ [k+1],α [k+1],β [k] = 0, which contradicts the definition of k + 1 As a conclusion, it follows that the transformation Γ is of the form 410 and 37 0, =0

21 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 21 Proof of Theorem 6 Let us consider two systems Σ 1 and Σ 2 and let Σ 1,CF : ẋ = Ax + Bv + m=m 0 f x and Σ 2,CF : ż = Az + Bw + m=m 1 f z denote, respectively, their canonical forms, where m 0 and m 1 denote the degrees of the first nonlinearizable homogeneous parts It is obvious that Σ 1 and Σ 2 are feedback equivalent if their canonical forms Σ 1,CF and Σ 2,CF coincide To prove the converse, we assume that the systems Σ 1 and Σ 2 are feedback equivalent while their canonical forms fail to coincide Since Σ 1 and Σ 2 are feedback equivalent, so are their canonical forms Σ 1,CF and Σ 2,CF It means that there exists a transformation Γ which brings Σ 1,CF into Σ 2,CF First remark that, from the definition of the integers m 0 and m 1, we necessarily have m 0 = m 1 Then, Theorem 2 and the fact that the components f [m0] [m0] and f are normalized see 35 ensure that f [m0] = f [m1] Let l be the largest integer such that for any i l, wehave f [m0+i 1] = f [m0+i 1] This means that the transformation Γ leaves invariant all terms of degree smaller than m 0 + l of the system Σ 1,CF Then Lemma 1 shows that the transformation Γ is of the form 410 Since the transformation Γ brings Σ 1,CF into Σ 2,CF, we deduce that 412 f [m0+l] = f [m0+l] +[ f [m0],φ [l+1] ] Following arguments in the proof of Theorem 5, we obtain where m0+l 2 ã [m0+l],s+2 x i1+l 1 x i2 2 xin s 2 m0 [m0] f + γ l, 2 xin s 2 x i1 1 xi2 2 xin s = m0+l 2 ā [m0+l],s+2 x i1+l 1 x i2 n s { ā [m0+l],i+2 } :, i and { ã [m0+l],i+2 } :, i denote, respectively, the set of m 0 + l-invariants associated with the homogeneous parts of degree m 0 + l of the systems Σ 1,CF and Σ 2,CF Using Theorem 2, the last identity can be rewritten as 413 Since m0+l [m0+l] f x i1+l 1 x i2 2 xin s = m0 [m0+l] f x 1, 0,,0 = x i1 1 xin s m0+l [m0+l] f x i1+l 1 x i2 2 xin s m0 f [m0+l] x i1 1 xin s m0 [m0] f + γ l x i1 1 xi2 2 xin s x 1, 0,,0=0,

22 22 ISSA AMADOU TALL AND WITOLD RESPONDEK the identity 413 gives γ l m0 [m0] f x i1 1 xi2 2 xin s which implies γ l =0, that is recall the definition of γ l, we have a l+1 = 0, and consequently φ [l+1],α [l+1],β [l] =0 Then the identity 412 reduces to f [m0+l] = f [m0+l], =0, which contradicts the definition of l We conclude that the canonical forms Σ 1,CF and Σ 2,CF coincide 5 Dual normal form and dual m-invariants In the normal form Σ NF given by 27, all the components of the control vector field g [m 1] are annihilated and all nonremovable nonlinearities are grouped in f Kang and Krener in their pioneering paper [26] showed that it is possible to transform, via a transformation Γ 2 of degree 2, the homogeneous system Σ [2] : ξ = Aξ + Bu + f [2] ξ+g [1] ξu into a dual normal form In that form the components of the drift f [2] are annihilated, while this time all nonremovable nonlinearities are present in g [1] The aim of this section is to propose, for an arbitrary m, a dual normal form for the system Σ and a dual normal form for the system Σ Our dual normal form on the one hand generalizes, for higher order terms, that given in [26] for second order terms, and on the other hand dualizes the normal form Σ NF The structure of this section will follow that of section 2: we will give the dual normal form, then define and study dual m-invariants; finally, we give an explicit construction of transformations bringing the system into its dual normal form Our first result asserts that we can always bring Σ to a dual normal form Theorem 7 feedback transformation Γ m, to the dual normal form Σ DNF 51 Σ DNF : The homogeneous system Σ is equivalent, via a homogeneous given by ẋ 1 = x 2, ẋ 2 = x 3 + vx n Q [m 2] 2,n x 1,,x n, ẋ = x +1 + v i=n +2 x i Q [m 2],i x 1,,x i, ẋ n 1 = x n + v n x i Q [m 2],i x 1,,x i, ẋ n = v, i=3 where Q [m 2],i x 1,,x i are homogeneous polynomials of degree m 2 depending on the indicated variables Theorem 7 follows from Theorem 9, which gives explicit transformation bringing Σ to its dual normal form Σ DNF, and thus we omit its proof

23 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 23 Now we will define dual m-invariants To start with, recall that the homogeneous vector field X [m 1] i is defined by taking the homogeneous part of degree m 1ofthe vector field X m 1 i = 1 i ad i Aξ+f B + g [m 1] By X [m 1] i π i ξ we will denote the vector field X [m 1] i evaluated at the point π i ξ = ξ 1,,ξ i, 0,,0 of the submanifold W i = {ξ R n : ξ i+1 = = ξ n =0} Consider the system Σ and, for any such that 2 n 1, define the polynomial b [m 1] by setting b [m 1] 1 + = g [m 1] k=1 L B L k 1 Aξ f k i=1 ξi L B L 1 Aξ CX [m 1] n i π i ξdξ i 0 The homogeneous polynomials b [m 1] for 2 n 1 will be called the dual m- invariants of the homogeneous system Σ Consider two systems Σ and Σ of the form 24 and 25 Let and { b [m 1] : 2 n 1 } [m 1] { b : 2 n 1 } denote, respectively, their dual m-invariants The following result gives a dualization of Theorem 2 Theorem 8 The dual m-invariants have the following properties: i Two systems Σ and Σ are equivalent via a homogeneous feedback transformation Γ m if and only if b [m 1] [m 1] = b for any 2 n 1 ii The dual m-invariants 51, are given by b [m 1] of the dual normal form Σ DNF, defined by b[m 1] x = i=n +2 x i Q [m 2],i x 1,,x i for any 2 n 1 The above result asserts that the dual m-invariants, as do the m-invariants, form a set of complete invariants of the homogeneous feedback transformation Notice, however, that the same information is encoded in both sets of invariants in different ways We will give a proof of Theorem 8 in section 7

24 24 ISSA AMADOU TALL AND WITOLD RESPONDEK 52 Now, we define the following homogeneous polynomials: ξi φ 1 = CX [m 1] n i π i ξdξ i, i=1 φ +1 = f α = β [m 1] = f 0 + L Aξ φ, 1 n 1, n + L Aξ φ n, + L B φ n g [m 1] n The next result gives an explicit construction of feedback transformations bringing the system Σ to its dual normal form Σ DNF Theorem 9 The feedback transformation Γ m : x = ξ + φ ξ, u = v + α ξ+β [m 1] ξv, where α, β [m 1], and the components φ of φ are defined by 52, brings the system Σ into its dual normal form Σ DNF given by 51 6 Dual canonical form Consider the system Σ of the form 21 and assume that its linear part F, G is controllable Apply successively to it a series of transformations Γ m, m =1, 2,, such that each Γ m brings Σ to its dual normal form Σ DNF ; for instance we can take a series of transformations defined by 52 Successive repeating of Theorem 9 gives the following dual normal form Theorem 10 The system Σ can be transformed via a formal feedback transformation Γ into the dual normal form Σ DNF given by 61 Σ DNF : ẋ 1 = x 2, ẋ 2 = x 3 + vx n Q 2,n x 1,,x n, ẋ = x +1 + v x i Q,i x 1,,x i, ẋ n 1 = x n + v ẋ n = v, i=n +2 x i Q,i x 1,,x i, i=3 where Q,i x 1,,x i are formal power series depending on the indicated variables Naturally, as with normal forms, a given system can admit different dual normal forms We are thus interested in constructing a dual canonical form Assuming that the linear part F, G of the system Σ, of the form 21, is controllable, we denote by m 0 the degree of the first homogeneous term of the system Σ which cannot be annihilated by a feedback transformation Thus, using Theorems 8 and 9, we can assume, after applying a suitable feedback, that Σ takes the form Σ : ξ = Aξ + Bu +ḡ [m 0 1] ξu + f ξ+g [m 1] ξu, m=m 0+1

25 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 25 where A, B is in Brunovský canonical form and the first nonvanishing homogeneous vector field ḡ [m0 1] is of the form ξ i Q [m0 2] ḡ [m0 1],i ξ 1,,ξ i, 2 n 1, ξ = i=n +2 0, = 1 and = n Define { = inf =2,,n 1 : ḡ [m0 1] } ξ 0 and let i 1,,i n such that i 1 + +i n = m 0 1 be the largest, in the lexicographic ordering, n-tuple of nonnegative integers such that m0 1 ḡ [m0 1] 0 ξ i1 1 ξin n We get the following result Theorem 11 There exists a formal feedback transformation Γ which brings the system Σ into the following one: where for any m m 0, 62 Moreover, ḡ [m 1] = Σ DCF :ẋ = Ax + Bv + i=n +2 m=m 0 ḡ [m 1] xv, x i Q [m 2],i x 1,,x i, 2 n 1, 0, =1and = n 63 and for any m m 0 +1 m0 1 ḡ [m0 1] x i1 1 xin n = ±1, 64 m0 1 ḡ [m 1] x 1, 0,,0=0 x i1 1 xin n The form Σ DCF, which satisfies 62, 63, and 64, will be called the dual canonical form of Σ The name is ustified by the following Theorem 12 The two systems Σ 1 and Σ 2 are formally feedback equivalent if and only if their dual canonical forms Σ 1,DCF and Σ 2,DCF coincide Example 6 Let us consider the system Σ : ξ = fξ+gξu, ξ R 3,u R,

26 26 ISSA AMADOU TALL AND WITOLD RESPONDEK whose linear part is assumed to be controllable Theorem 10 ensures that the system Σ is formally feedback equivalent to the dual normal form Σ DNF given by ẋ 1 = x 2, ẋ 2 = x 3 + vx 3 Qx 1,x 2,x 3, ẋ 3 = v, where Qx 1,x 2,x 3 is a formal power series of variables x 1,x 2,x 3 Assume for simplicity that m 0 = 2, which is equivalent to the condition that g, ad f g, and [g, ad f g] are linearly independent at 0 R 3 This implies that we can represent Q = Qx 1,x 2,x 3 by Q = c + x 1 Q 1 x 1 +x 2 Q 2 x 1,x 2 +x 3 Q 3 x 1,x 2,x 3, where c R, c 0 Observe that any Q of the above form gives a dual normal form Σ DNF In order to get the dual canonical form we use Theorem 11, which ensures that the system Σ is formally feedback equivalent to its dual canonical form Σ DCF defined by x 1 = x 2, x 2 = x 3 +ṽ x 3 Q x1, x 2, x 3, x 3 =ṽ, where Q x 1, x 2, x 3 is a formal power series such that Q x 1, x 2, x 3 =1+ x 2 Q2 x 1, x 2 + x 3 Q3 x 1, x 2, x 3 7 Proofs of dual results In this section, we prove our dual results The proof of Theorem 7 will be omitted because in the proof of Theorem 9 we give an explicit homogeneous feedback transformation bringing a given homogeneous system into its dual normal form Theorem 10 follows from a successive application of Theorem 7 We will thus prove Theorems 8, 9, 11, and 12 Proof of Theorem 8 i We will prove that if the system Σ is equivalent to b [m 1] Σ via a transformation Γ m, then their dual m-invariants b [m 1] and The action of Γ m can be decomposed into that of a pure feedback of the form followed by that of a diffeomorphism u = v + α ξ+β [m 1] ξv x = ξ + φ ξ coincide of the state space Since the first n 1 components of the vector fields f and g [m 1], as well as those of X [m 1] n i, are invariant under pure feedback, we can conclude that the functions b [m 1] for 2 n 1 are invariant under pure feedback It remains to prove that they are also invariant under any diffeomorphism x =Φξ of the form Φξ =ξ + φ ξ The diffeomorphism Φ brings the system Σ into the form Σ : ẋ = Ax + Bu + f x+ g [m 1] xu,

27 FEEDBACK CLASSIFICATION OF NONLINEAR SYSTEMS 27 where b [m 1] f = f +[Ax, φ ], g [m 1] = g [m 1] + L B φ Denoting by b [m 1] and respectively, with the homogeneous systems Σ and Σ,weget where and ˆb[m 1] for 2 n 1 the dual m-invariants associated, b[m 1] x =L B φ i=1 = b [m 1] 2 x+ L B L 1 Ax k=0 xi 0 [m 1] + ˆb, L B L k 2 Ax CA k ad Ax φ C [m 1] ˆX n i π i xdx i ˆX [m 1] n i x = 1 n i ad n i Ax+[Ax,φ ] B + L Bφ = 1 n i Φ ad n i Aξ Bx = A n i B + L A n i Bφ x We can deduce that ˆb[m 1] which gives x =L B φ i=1 = L B φ 1 x+ L B L 1 Ax k=1 xi 0 L B L k Ax φ k 1 x+ L B L k k=1 k=2 L A n i Bφ 1 π i xdx i Ax φ k b[m 1] L B L k Ax φ k L B L k k=2 = b [m 1] Ax φ k L B L 1 Ax φ 1 x =0, Thus the functions b [m 1] are invariant under any diffeomorphism of the form x = Φξ =ξ + φ ξ Therefore they remain invariant under the transformation Γ m The fact that two homogeneous systems, whose dual m-invariants coincide, are feedback equivalent follows clearly from item ii of the theorem, which will be proved below Indeed, by item ii, both systems coincide when transformed to their canonical forms [m 1] ii Denote by b for 2 n 1 the dual m-invariants associated with the dual normal form Σ DNF They are given by b[m 1] =ḡ [m 1] i=1 xi L B L 1 [m 1] Ax C X n i π i xdx i, 0

28 28 ISSA AMADOU TALL AND WITOLD RESPONDEK where the components ḡ [m 1] C are given by 51 and X [m 1] n i = 1 n i Cad n i Ax ḡ[m 1] X [m 1] n i It suffices to observe see Lemma 2 below that, on the one hand, C is a linear combination of functions L s Axḡ[m 1] for 0 s n i and 1 n i + 1 and, on the other hand, ḡ [m 1] π i x = 0 for all such that 1 n i +1 We thus conclude that C π i x=0, which implies X [m 1] n i b[m 1] =ḡ [m 1] for any such that 2 n 1 Proof of Theorem 9 Denote by Σ : ẋ = Ax + Bv + f x+ g [m 1] xv the system Σ transformed via a homogeneous feedback transformation Γ m defined by 52 From Proposition 1, it follows that for Σ we have 71 f = 0 for 1 n, g [m 1] = 0 for = 1 and = n, g [m 1] = g [m 1] + L B φ for 2 n 1 It thus suffices to show that the components g [m 1] for 2 n 1 are in the dual normal form 51 We prove easily by an induction argument that which allows us to show that g [m 1] +1 = g [m 1] +1 + Now, from the identity φ +1 = L B L Aξ φ 1 = k=1 L k Aξ f k k=1 k L k k=0 L B L k 1 Aξ f L A Bφ k 1 = CX [m 1] π n k ξ we can deduce that k g [m 1] +1 = g [m 1] +1 + k=1 Aξ + L Aξ φ 1, 1, Aξ L A k Bφ k+1 + i=n k+1 L B L k 1 f k+1 k L k Aξ k=1 i=n k+1 ξi 0 k=0 ξi 0 k L k 1 Aξ L A k Bφ CX [m 1] n i π i ξ dξ i, ξ n k k L k Aξ CX k k=0 CX [m 1] n i ξ n k π n k ξ π i ξ dξ i

On the convergence of normal forms for analytic control systems

On the convergence of normal forms for analytic control systems Chapter 1 On the convergence of normal forms for analytic control systems Wei Kang Department of Mathematics, Naval Postgraduate School Monterey, CA 93943 wkang@nps.navy.mil Arthur J. Krener Department

More information

10. Smooth Varieties. 82 Andreas Gathmann

10. Smooth Varieties. 82 Andreas Gathmann 82 Andreas Gathmann 10. Smooth Varieties Let a be a point on a variety X. In the last chapter we have introduced the tangent cone C a X as a way to study X locally around a (see Construction 9.20). It

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

Chapter 2: Linear Independence and Bases

Chapter 2: Linear Independence and Bases MATH20300: Linear Algebra 2 (2016 Chapter 2: Linear Independence and Bases 1 Linear Combinations and Spans Example 11 Consider the vector v (1, 1 R 2 What is the smallest subspace of (the real vector space

More information

LECTURE 5: THE METHOD OF STATIONARY PHASE

LECTURE 5: THE METHOD OF STATIONARY PHASE LECTURE 5: THE METHOD OF STATIONARY PHASE Some notions.. A crash course on Fourier transform For j =,, n, j = x j. D j = i j. For any multi-index α = (α,, α n ) N n. α = α + + α n. α! = α! α n!. x α =

More information

Chap. 3. Controlled Systems, Controllability

Chap. 3. Controlled Systems, Controllability Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :

More information

A method for computing quadratic Brunovsky forms

A method for computing quadratic Brunovsky forms Electronic Journal of Linear Algebra Volume 13 Volume 13 (25) Article 3 25 A method for computing quadratic Brunovsky forms Wen-Long Jin wjin@uciedu Follow this and additional works at: http://repositoryuwyoedu/ela

More information

Introduction to Geometric Control

Introduction to Geometric Control Introduction to Geometric Control Relative Degree Consider the square (no of inputs = no of outputs = p) affine control system ẋ = f(x) + g(x)u = f(x) + [ g (x),, g p (x) ] u () h (x) y = h(x) = (2) h

More information

LIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS. Suba Thomas Harry G. Kwatny Bor-Chin Chang

LIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS. Suba Thomas Harry G. Kwatny Bor-Chin Chang LIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS Suba Thomas Harry G Kwatny Bor-Chin Chang Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, PA 1914, USA

More information

COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY

COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY BIN HAN AND HUI JI Abstract. In this paper, we provide a family of compactly supported orthonormal complex wavelets with dilation

More information

Stabilization and Passivity-Based Control

Stabilization and Passivity-Based Control DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive

More information

A linear algebra proof of the fundamental theorem of algebra

A linear algebra proof of the fundamental theorem of algebra A linear algebra proof of the fundamental theorem of algebra Andrés E. Caicedo May 18, 2010 Abstract We present a recent proof due to Harm Derksen, that any linear operator in a complex finite dimensional

More information

Algebra Review 2. 1 Fields. A field is an extension of the concept of a group.

Algebra Review 2. 1 Fields. A field is an extension of the concept of a group. Algebra Review 2 1 Fields A field is an extension of the concept of a group. Definition 1. A field (F, +,, 0 F, 1 F ) is a set F together with two binary operations (+, ) on F such that the following conditions

More information

A linear algebra proof of the fundamental theorem of algebra

A linear algebra proof of the fundamental theorem of algebra A linear algebra proof of the fundamental theorem of algebra Andrés E. Caicedo May 18, 2010 Abstract We present a recent proof due to Harm Derksen, that any linear operator in a complex finite dimensional

More information

Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone

Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 9, Number 2, pp. 227 237 (2014) http://campus.mst.edu/adsa Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone

More information

Math 210B. Artin Rees and completions

Math 210B. Artin Rees and completions Math 210B. Artin Rees and completions 1. Definitions and an example Let A be a ring, I an ideal, and M an A-module. In class we defined the I-adic completion of M to be M = lim M/I n M. We will soon show

More information

Centers of projective vector fields of spatial quasi-homogeneous systems with weight (m, m, n) and degree 2 on the sphere

Centers of projective vector fields of spatial quasi-homogeneous systems with weight (m, m, n) and degree 2 on the sphere Electronic Journal of Qualitative Theory of Differential Equations 2016 No. 103 1 26; doi: 10.14232/ejqtde.2016.1.103 http://www.math.u-szeged.hu/ejqtde/ Centers of projective vector fields of spatial

More information

Disturbance Decoupling Problem

Disturbance Decoupling Problem DISC Systems and Control Theory of Nonlinear Systems, 21 1 Disturbance Decoupling Problem Lecture 4 Nonlinear Dynamical Control Systems, Chapter 7 The disturbance decoupling problem is a typical example

More information

Feedback Linearizability of Strict Feedforward Systems

Feedback Linearizability of Strict Feedforward Systems Proceedings o the 47th IEEE Conerence on Decision and Control Cancun, Mexico, Dec 9-11, 28 WeB11 Feedback Linearizability o Strict Feedorward Systems Witold Respondek and Issa Amadou Tall Abstract For

More information

Chap. 1. Some Differential Geometric Tools

Chap. 1. Some Differential Geometric Tools Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U

More information

Pattern generation, topology, and non-holonomic systems

Pattern generation, topology, and non-holonomic systems Systems & Control Letters ( www.elsevier.com/locate/sysconle Pattern generation, topology, and non-holonomic systems Abdol-Reza Mansouri Division of Engineering and Applied Sciences, Harvard University,

More information

SYMPLECTIC GEOMETRY: LECTURE 5

SYMPLECTIC GEOMETRY: LECTURE 5 SYMPLECTIC GEOMETRY: LECTURE 5 LIAT KESSLER Let (M, ω) be a connected compact symplectic manifold, T a torus, T M M a Hamiltonian action of T on M, and Φ: M t the assoaciated moment map. Theorem 0.1 (The

More information

Abstract. + n) variables with a rank of ( n2 (n+1) )+n 2 ] equations in ( n2 (n+1)

Abstract. + n) variables with a rank of ( n2 (n+1) )+n 2 ] equations in ( n2 (n+1) SOLUTION OF SECOND ORDER LINEARIZATION R.Devanathan School of Electrical and Electronic Engineering (Block S) Nanyang Technological University Singapore 6998 E-mail:edevan@ntu.edu.sg Abstract For a nonlinear

More information

Convergence Rate of Nonlinear Switched Systems

Convergence Rate of Nonlinear Switched Systems Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the

More information

We denote the space of distributions on Ω by D ( Ω) 2.

We denote the space of distributions on Ω by D ( Ω) 2. Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study

More information

CSL361 Problem set 4: Basic linear algebra

CSL361 Problem set 4: Basic linear algebra CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices

More information

An extended complete Chebyshev system of 3 Abelian integrals related to a non-algebraic Hamiltonian system

An extended complete Chebyshev system of 3 Abelian integrals related to a non-algebraic Hamiltonian system Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 6, No. 4, 2018, pp. 438-447 An extended complete Chebyshev system of 3 Abelian integrals related to a non-algebraic Hamiltonian

More information

Equivalence of dynamical systems by bisimulation

Equivalence of dynamical systems by bisimulation Equivalence of dynamical systems by bisimulation Arjan van der Schaft Department of Applied Mathematics, University of Twente P.O. Box 217, 75 AE Enschede, The Netherlands Phone +31-53-4893449, Fax +31-53-48938

More information

CANONICAL FORMS FOR LINEAR TRANSFORMATIONS AND MATRICES. D. Katz

CANONICAL FORMS FOR LINEAR TRANSFORMATIONS AND MATRICES. D. Katz CANONICAL FORMS FOR LINEAR TRANSFORMATIONS AND MATRICES D. Katz The purpose of this note is to present the rational canonical form and Jordan canonical form theorems for my M790 class. Throughout, we fix

More information

Fourier Transform & Sobolev Spaces

Fourier Transform & Sobolev Spaces Fourier Transform & Sobolev Spaces Michael Reiter, Arthur Schuster Summer Term 2008 Abstract We introduce the concept of weak derivative that allows us to define new interesting Hilbert spaces the Sobolev

More information

MOSCOW STATE UNIVERSITY DEPARTMENT OF MECHANICS AND MATHEMATICS

MOSCOW STATE UNIVERSITY DEPARTMENT OF MECHANICS AND MATHEMATICS MOSCOW STATE UNIVERSITY DEPARTMENT OF MECHANICS AND MATHEMATICS Hildebrand Roland CLASSIFICATION OF PHASE PORTRAITS OF OPTIMAL SYNTHESIS specialization: differential equations T h e s i s s u m m a r y

More information

Lagrangian submanifolds and generating functions

Lagrangian submanifolds and generating functions Chapter 4 Lagrangian submanifolds and generating functions Motivated by theorem 3.9 we will now study properties of the manifold Λ φ X (R n \{0}) for a clean phase function φ. As shown in section 3.3 Λ

More information

Explicit design of Time-varying Stabilizing Control Laws for a Class of Controllable Systems without Drift

Explicit design of Time-varying Stabilizing Control Laws for a Class of Controllable Systems without Drift Systems & Control Lett., vol. 18, p. 147-158, 1992. Explicit design of Time-varying Stabilizing Control Laws for a Class of Controllable Systems without Drift Jean-Baptiste Pomet May 23, 1991, revised

More information

INVERSE FUNCTION THEOREM and SURFACES IN R n

INVERSE FUNCTION THEOREM and SURFACES IN R n INVERSE FUNCTION THEOREM and SURFACES IN R n Let f C k (U; R n ), with U R n open. Assume df(a) GL(R n ), where a U. The Inverse Function Theorem says there is an open neighborhood V U of a in R n so that

More information

THE GROUP OF UNITS OF SOME FINITE LOCAL RINGS I

THE GROUP OF UNITS OF SOME FINITE LOCAL RINGS I J Korean Math Soc 46 (009), No, pp 95 311 THE GROUP OF UNITS OF SOME FINITE LOCAL RINGS I Sung Sik Woo Abstract The purpose of this paper is to identify the group of units of finite local rings of the

More information

1 The Observability Canonical Form

1 The Observability Canonical Form NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)

More information

Zeros and zero dynamics

Zeros and zero dynamics CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)

More information

1 Relative degree and local normal forms

1 Relative degree and local normal forms THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a

More information

Observability and forward-backward observability of discrete-time nonlinear systems

Observability and forward-backward observability of discrete-time nonlinear systems Observability and forward-backward observability of discrete-time nonlinear systems Francesca Albertini and Domenico D Alessandro Dipartimento di Matematica pura a applicata Universitá di Padova, 35100

More information

2. Intersection Multiplicities

2. Intersection Multiplicities 2. Intersection Multiplicities 11 2. Intersection Multiplicities Let us start our study of curves by introducing the concept of intersection multiplicity, which will be central throughout these notes.

More information

Math 113 Winter 2013 Prof. Church Midterm Solutions

Math 113 Winter 2013 Prof. Church Midterm Solutions Math 113 Winter 2013 Prof. Church Midterm Solutions Name: Student ID: Signature: Question 1 (20 points). Let V be a finite-dimensional vector space, and let T L(V, W ). Assume that v 1,..., v n is a basis

More information

7 Planar systems of linear ODE

7 Planar systems of linear ODE 7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

More information

A Sufficient Condition for Local Controllability

A Sufficient Condition for Local Controllability A Sufficient Condition for Local Controllability of Nonlinear Systems Along Closed Orbits Kwanghee Nam and Aristotle Arapostathis Abstract We present a computable sufficient condition to determine local

More information

THEOREM OF OSELEDETS. We recall some basic facts and terminology relative to linear cocycles and the multiplicative ergodic theorem of Oseledets [1].

THEOREM OF OSELEDETS. We recall some basic facts and terminology relative to linear cocycles and the multiplicative ergodic theorem of Oseledets [1]. THEOREM OF OSELEDETS We recall some basic facts and terminology relative to linear cocycles and the multiplicative ergodic theorem of Oseledets []. 0.. Cocycles over maps. Let µ be a probability measure

More information

1 Lyapunov theory of stability

1 Lyapunov theory of stability M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

More information

Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra

Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra D. R. Wilkins Contents 3 Topics in Commutative Algebra 2 3.1 Rings and Fields......................... 2 3.2 Ideals...............................

More information

HIGHER ORDER CRITERION FOR THE NONEXISTENCE OF FORMAL FIRST INTEGRAL FOR NONLINEAR SYSTEMS. 1. Introduction

HIGHER ORDER CRITERION FOR THE NONEXISTENCE OF FORMAL FIRST INTEGRAL FOR NONLINEAR SYSTEMS. 1. Introduction Electronic Journal of Differential Equations, Vol. 217 (217), No. 274, pp. 1 11. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu HIGHER ORDER CRITERION FOR THE NONEXISTENCE

More information

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered

More information

i=1 α i. Given an m-times continuously

i=1 α i. Given an m-times continuously 1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable

More information

be any ring homomorphism and let s S be any element of S. Then there is a unique ring homomorphism

be any ring homomorphism and let s S be any element of S. Then there is a unique ring homomorphism 21. Polynomial rings Let us now turn out attention to determining the prime elements of a polynomial ring, where the coefficient ring is a field. We already know that such a polynomial ring is a UFD. Therefore

More information

Exercise Solutions to Functional Analysis

Exercise Solutions to Functional Analysis Exercise Solutions to Functional Analysis Note: References refer to M. Schechter, Principles of Functional Analysis Exersize that. Let φ,..., φ n be an orthonormal set in a Hilbert space H. Show n f n

More information

Decomposing Bent Functions

Decomposing Bent Functions 2004 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 49, NO. 8, AUGUST 2003 Decomposing Bent Functions Anne Canteaut and Pascale Charpin Abstract In a recent paper [1], it is shown that the restrictions

More information

Linear Codes, Target Function Classes, and Network Computing Capacity

Linear Codes, Target Function Classes, and Network Computing Capacity Linear Codes, Target Function Classes, and Network Computing Capacity Rathinakumar Appuswamy, Massimo Franceschetti, Nikhil Karamchandani, and Kenneth Zeger IEEE Transactions on Information Theory Submitted:

More information

THE NUMBER OF POLYNOMIAL SOLUTIONS OF POLYNOMIAL RICCATI EQUATIONS

THE NUMBER OF POLYNOMIAL SOLUTIONS OF POLYNOMIAL RICCATI EQUATIONS This is a preprint of: The number of polynomial solutions of polynomial Riccati equations, Armengol Gasull, Joan Torregrosa, Xiang Zhang, J. Differential Equations, vol. 261, 5071 5093, 2016. DOI: [10.1016/j.jde.2016.07.019]

More information

Disturbance Attenuation for a Class of Nonlinear Systems by Output Feedback

Disturbance Attenuation for a Class of Nonlinear Systems by Output Feedback Disturbance Attenuation for a Class of Nonlinear Systems by Output Feedback Wei in Chunjiang Qian and Xianqing Huang Submitted to Systems & Control etters /5/ Abstract This paper studies the problem of

More information

Feedback Linearization

Feedback Linearization Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)

More information

Math 443 Differential Geometry Spring Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook.

Math 443 Differential Geometry Spring Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook. Math 443 Differential Geometry Spring 2013 Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook. Endomorphisms of a Vector Space This handout discusses

More information

RIEMANN SURFACES. max(0, deg x f)x.

RIEMANN SURFACES. max(0, deg x f)x. RIEMANN SURFACES 10. Weeks 11 12: Riemann-Roch theorem and applications 10.1. Divisors. The notion of a divisor looks very simple. Let X be a compact Riemann surface. A divisor is an expression a x x x

More information

Algebraic Varieties. Notes by Mateusz Micha lek for the lecture on April 17, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra

Algebraic Varieties. Notes by Mateusz Micha lek for the lecture on April 17, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra Algebraic Varieties Notes by Mateusz Micha lek for the lecture on April 17, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra Algebraic varieties represent solutions of a system of polynomial

More information

FILTERED RINGS AND MODULES. GRADINGS AND COMPLETIONS.

FILTERED RINGS AND MODULES. GRADINGS AND COMPLETIONS. FILTERED RINGS AND MODULES. GRADINGS AND COMPLETIONS. Let A be a ring, for simplicity assumed commutative. A filtering, or filtration, of an A module M means a descending sequence of submodules M = M 0

More information

Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space

Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................

More information

ALGEBRA HW 10 CLAY SHONKWILER. Proof. Suppose that these three functions are not linearly independent. Then there exist a, b R such that

ALGEBRA HW 10 CLAY SHONKWILER. Proof. Suppose that these three functions are not linearly independent. Then there exist a, b R such that ALGEBRA HW 10 CLAY SHONKWILER 1 (a): If A is a square matrix satisfying A 3 A 2 + A I = 0, find A 1 in terms of A. Answer: If A 3 A 2 + A I = 0, then (A 2 A + I)A = A 3 A 2 + A = I = A 3 A 2 + A = A(A

More information

Tangent spaces, normals and extrema

Tangent spaces, normals and extrema Chapter 3 Tangent spaces, normals and extrema If S is a surface in 3-space, with a point a S where S looks smooth, i.e., without any fold or cusp or self-crossing, we can intuitively define the tangent

More information

Takao Akahori. z i In this paper, if f is a homogeneous polynomial, the correspondence between the Kodaira-Spencer class and C[z 1,...

Takao Akahori. z i In this paper, if f is a homogeneous polynomial, the correspondence between the Kodaira-Spencer class and C[z 1,... J. Korean Math. Soc. 40 (2003), No. 4, pp. 667 680 HOMOGENEOUS POLYNOMIAL HYPERSURFACE ISOLATED SINGULARITIES Takao Akahori Abstract. The mirror conjecture means originally the deep relation between complex

More information

Math 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.

Math 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces. Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,

More information

5 Quiver Representations

5 Quiver Representations 5 Quiver Representations 5. Problems Problem 5.. Field embeddings. Recall that k(y,..., y m ) denotes the field of rational functions of y,..., y m over a field k. Let f : k[x,..., x n ] k(y,..., y m )

More information

SOME COMMENTS ON OBSERVABILITY NORMAL FORM AND STEP-BY-STEP SLIDING MODE OBSERVER

SOME COMMENTS ON OBSERVABILITY NORMAL FORM AND STEP-BY-STEP SLIDING MODE OBSERVER SOME COMMENTS ON OBSERVABILITY NORMAL FORM AND STEP-BY-STEP SLIDING MODE OBSERVER G ZHENG 1, L BOUTAT-BADDAS 2, JP BARBOT 1 AND D BOUTAT 3 1 EQUIPE COMMANDE DES SYSTÈMES (ECS), ENSEA, 6 AV DU PONCEAU,

More information

Biorthogonal Spline Type Wavelets

Biorthogonal Spline Type Wavelets PERGAMON Computers and Mathematics with Applications 0 (00 1 0 www.elsevier.com/locate/camwa Biorthogonal Spline Type Wavelets Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan

More information

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix

More information

VERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE

VERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE PHYSCON 2013 San Luis Potosí México August 26 August 29 2013 VERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE M Isabel García-Planas Departamento de Matemàtica Aplicada I Universitat

More information

Coordinate Systems and Canonical Forms

Coordinate Systems and Canonical Forms Appendix D Coordinate Systems and Canonical Forms D.1. Local Coordinates Let O be an open set in R n. We say that an n-tuple of smooth realvalued functions defined in O, (φ 1,...,φ n ), forms a local coordinate

More information

The extreme rays of the 5 5 copositive cone

The extreme rays of the 5 5 copositive cone The extreme rays of the copositive cone Roland Hildebrand March 8, 0 Abstract We give an explicit characterization of all extreme rays of the cone C of copositive matrices. The results are based on the

More information

Lecture 2: Controllability of nonlinear systems

Lecture 2: Controllability of nonlinear systems DISC Systems and Control Theory of Nonlinear Systems 1 Lecture 2: Controllability of nonlinear systems Nonlinear Dynamical Control Systems, Chapter 3 See www.math.rug.nl/ arjan (under teaching) for info

More information

Nonlinear Control Systems

Nonlinear Control Systems Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 7. Feedback Linearization IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs1/ 1 1 Feedback Linearization Given a nonlinear

More information

LECTURE 25-26: CARTAN S THEOREM OF MAXIMAL TORI. 1. Maximal Tori

LECTURE 25-26: CARTAN S THEOREM OF MAXIMAL TORI. 1. Maximal Tori LECTURE 25-26: CARTAN S THEOREM OF MAXIMAL TORI 1. Maximal Tori By a torus we mean a compact connected abelian Lie group, so a torus is a Lie group that is isomorphic to T n = R n /Z n. Definition 1.1.

More information

8. Prime Factorization and Primary Decompositions

8. Prime Factorization and Primary Decompositions 70 Andreas Gathmann 8. Prime Factorization and Primary Decompositions 13 When it comes to actual computations, Euclidean domains (or more generally principal ideal domains) are probably the nicest rings

More information

2 Two-Point Boundary Value Problems

2 Two-Point Boundary Value Problems 2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x

More information

On the Harish-Chandra Embedding

On the Harish-Chandra Embedding On the Harish-Chandra Embedding The purpose of this note is to link the Cartan and the root decompositions. In addition, it explains how we can view a Hermitian symmetric domain as G C /P where P is a

More information

Math 225B: Differential Geometry, Final

Math 225B: Differential Geometry, Final Math 225B: Differential Geometry, Final Ian Coley March 5, 204 Problem Spring 20,. Show that if X is a smooth vector field on a (smooth) manifold of dimension n and if X p is nonzero for some point of

More information

HASSE-MINKOWSKI THEOREM

HASSE-MINKOWSKI THEOREM HASSE-MINKOWSKI THEOREM KIM, SUNGJIN 1. Introduction In rough terms, a local-global principle is a statement that asserts that a certain property is true globally if and only if it is true everywhere locally.

More information

Part V. 17 Introduction: What are measures and why measurable sets. Lebesgue Integration Theory

Part V. 17 Introduction: What are measures and why measurable sets. Lebesgue Integration Theory Part V 7 Introduction: What are measures and why measurable sets Lebesgue Integration Theory Definition 7. (Preliminary). A measure on a set is a function :2 [ ] such that. () = 2. If { } = is a finite

More information

Constructions with ruler and compass.

Constructions with ruler and compass. Constructions with ruler and compass. Semyon Alesker. 1 Introduction. Let us assume that we have a ruler and a compass. Let us also assume that we have a segment of length one. Using these tools we can

More information

Linear Algebra. Min Yan

Linear Algebra. Min Yan Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................

More information

arxiv: v2 [cs.dm] 29 Mar 2013

arxiv: v2 [cs.dm] 29 Mar 2013 arxiv:1302.6346v2 [cs.dm] 29 Mar 2013 Fixed point theorems for Boolean networks expressed in terms of forbidden subnetworks Adrien Richard Laboratoire I3S, CNRS & Université de Nice-Sophia Antipolis, France.

More information

ON THE CONTINUED FRACTION EXPANSION OF THE UNIQUE ROOT IN F(p) OF THE EQUATION x 4 + x 2 T x 1/12 = 0 AND OTHER RELATED HYPERQUADRATIC EXPANSIONS

ON THE CONTINUED FRACTION EXPANSION OF THE UNIQUE ROOT IN F(p) OF THE EQUATION x 4 + x 2 T x 1/12 = 0 AND OTHER RELATED HYPERQUADRATIC EXPANSIONS ON THE CONTINUED FRACTION EXPANSION OF THE UNIQUE ROOT IN F(p) OF THE EQUATION x 4 + x 2 T x 1/12 = 0 AND OTHER RELATED HYPERQUADRATIC EXPANSIONS by A. Lasjaunias Abstract.In 1986, Mills and Robbins observed

More information

Proof. We indicate by α, β (finite or not) the end-points of I and call

Proof. We indicate by α, β (finite or not) the end-points of I and call C.6 Continuous functions Pag. 111 Proof of Corollary 4.25 Corollary 4.25 Let f be continuous on the interval I and suppose it admits non-zero its (finite or infinite) that are different in sign for x tending

More information

An introduction to Birkhoff normal form

An introduction to Birkhoff normal form An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an

More information

Bare-bones outline of eigenvalue theory and the Jordan canonical form

Bare-bones outline of eigenvalue theory and the Jordan canonical form Bare-bones outline of eigenvalue theory and the Jordan canonical form April 3, 2007 N.B.: You should also consult the text/class notes for worked examples. Let F be a field, let V be a finite-dimensional

More information

LINEAR ALGEBRA W W L CHEN

LINEAR ALGEBRA W W L CHEN LINEAR ALGEBRA W W L CHEN c W W L Chen, 994, 28. This chapter is available free to all individuals, on the understanding that it is not to be used for financial gain, and may be downloaded and/or photocopied,

More information

Remarks on Quadratic Hamiltonians in Spaceflight Mechanics

Remarks on Quadratic Hamiltonians in Spaceflight Mechanics Remarks on Quadratic Hamiltonians in Spaceflight Mechanics Bernard Bonnard 1, Jean-Baptiste Caillau 2, and Romain Dujol 2 1 Institut de mathématiques de Bourgogne, CNRS, Dijon, France, bernard.bonnard@u-bourgogne.fr

More information

Group Gradings on Finite Dimensional Lie Algebras

Group Gradings on Finite Dimensional Lie Algebras Algebra Colloquium 20 : 4 (2013) 573 578 Algebra Colloquium c 2013 AMSS CAS & SUZHOU UNIV Group Gradings on Finite Dimensional Lie Algebras Dušan Pagon Faculty of Natural Sciences and Mathematics, University

More information

Definition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.

Definition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M. 5 Vector fields Last updated: March 12, 2012. 5.1 Definition and general properties We first need to define what a vector field is. Definition 5.1. A vector field v on a manifold M is map M T M such that

More information

Solution to Homework 1

Solution to Homework 1 Solution to Homework Sec 2 (a) Yes It is condition (VS 3) (b) No If x, y are both zero vectors Then by condition (VS 3) x = x + y = y (c) No Let e be the zero vector We have e = 2e (d) No It will be false

More information

Approximation of complex algebraic numbers by algebraic numbers of bounded degree

Approximation of complex algebraic numbers by algebraic numbers of bounded degree Approximation of complex algebraic numbers by algebraic numbers of bounded degree YANN BUGEAUD (Strasbourg) & JAN-HENDRIK EVERTSE (Leiden) Abstract. To measure how well a given complex number ξ can be

More information

Some notes on Coxeter groups

Some notes on Coxeter groups Some notes on Coxeter groups Brooks Roberts November 28, 2017 CONTENTS 1 Contents 1 Sources 2 2 Reflections 3 3 The orthogonal group 7 4 Finite subgroups in two dimensions 9 5 Finite subgroups in three

More information

Solutions. Problems of Chapter 1

Solutions. Problems of Chapter 1 Solutions Problems of Chapter 1 1.1 A real square matrix A IR n n is invertible if and only if its determinant is non zero. The determinant of A is a polynomial in the entries a ij of A, whose set of zeros

More information

Putzer s Algorithm. Norman Lebovitz. September 8, 2016

Putzer s Algorithm. Norman Lebovitz. September 8, 2016 Putzer s Algorithm Norman Lebovitz September 8, 2016 1 Putzer s algorithm The differential equation dx = Ax, (1) dt where A is an n n matrix of constants, possesses the fundamental matrix solution exp(at),

More information

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/793: Analysis and Control of Nonlinear Systems MCE693/793: Analysis and Control of Nonlinear Systems Input-Output and Input-State Linearization Zero Dynamics of Nonlinear Systems Hanz Richter Mechanical Engineering Department Cleveland State University

More information

Geometric control and dynamical systems

Geometric control and dynamical systems Université de Nice - Sophia Antipolis & Institut Universitaire de France 9th AIMS International Conference on Dynamical Systems, Differential Equations and Applications Control of an inverted pendulum

More information

SYMMETRIC SUBGROUP ACTIONS ON ISOTROPIC GRASSMANNIANS

SYMMETRIC SUBGROUP ACTIONS ON ISOTROPIC GRASSMANNIANS 1 SYMMETRIC SUBGROUP ACTIONS ON ISOTROPIC GRASSMANNIANS HUAJUN HUANG AND HONGYU HE Abstract. Let G be the group preserving a nondegenerate sesquilinear form B on a vector space V, and H a symmetric subgroup

More information