HAYHHBIH BECTHHK MrTY TA SCIENTIFIC HERALD of MSTUCA
UDC 517.935.4
FLATNESS CONDITIONS FOR SYSTEMS WITH TWO INPUTS
V.N. CHETVERIKOV1
Flat systems with 2 inputs are investigated. Our approach is based on invertible differential operators and deformations of structures on diffieties of control systems. Invertible differential operators of size 2x2 are described. The order of a flat output is estimated by the order of the corresponding invertible differential operator linearizing the control system. The minimal order of invertible differential operator linearizing the system is estimated by the order of its deformation.
Keywords: Nonlinear systems, infinite-order prolongations, flat control systems, dynamic feedback linearization, invertible differential operators.
1. INTRODUCTION
Flat control systems were introduced by Fliess, Levine, Martin, and Rouchon [1]. Later, it was shown that many classes of systems commonly used in nonlinear control theory are flat [2]. Besides, control design for flat systems was developed (see [2]). These facts explain the interest in looking for flatness conditions. In the case of one input, flatness conditions are well known (see [2]). In the general case, checking whether a given system is flat still remains up an open problem.
Using geometric methods, Aranda-Bricaire, Moog, and Pomet introduced an infinitesimal Brunovsky form for nonlinear systems and obtained a necessary and sufficient condition for flatness [3]. Namely, flatness of a control system means existence of an invertible linear differential operator of a certain type satisfying some conditions (see Theorem 2 below). In geometry of differential equations [4] the operators of this type are investigated and are called C -differential. A description of invertible linear differential operators was recently obtained in [5].
In this paper, we consider the case of two inputs. Our approach is based on results of [3] and [5]. Here, the description of invertible operators is generalized to the C -differential case. We also show that invertible C -differential 2x2 operators have a simple structure. The main results of the paper follow from this structure.
2. FLAT SYSTEMS
Consider a system of the form
x = f(t, x, u), x e Rn, u e Rm, (1)
where t is the independent variable, the vector x = (Xj,...,xn) is a state, the vector u = (uj,...,um) is a control, its coordinates Uj,...,um are inputs of the system, f = (fn)
is a smooth vector function, and X = dx / dt. Here and throughout the following,
smoothness is understood as infinite differentiability. System (1) is said to be
regular if rank (df /du) = m for all considered values of the variables t, x, and u .
1 This work was supported in part by RFBR grants 13-07-00736 and 14-01-00424.
Let k be some nonnegative integer. We treat
(k)
t, x ,... , Xn , ,... , um , ..., um, ••• , um
as independent variables and consider the space with these coordinates. By O(k) denote some domain in this space.
A regular system (1) is said to be flat [1] in the domain O(k) if there exist functions
y = h (t, x, u, u,..., u(k)), i = 1,., m, (2)
defined in O( k) such that
(i) the variables x and u can be expressed via t, the functions (2), and their time-derivatives up to some finite order;
(ii) any finite set of the functions (2), their time-derivatives, and the function t is functionally independent.
In this case, the set of functions (2) is called a flat (or linearizing) output of system (1).
We shall need the following notions and facts of geometric theory of differential equations [4]. To each system (1) assign the infinite-dimensional space R° with coordinates
(t, x, u(0), u(1)u(I),...), (3)
where u(I) denotes the vector variable corresponding to the l -th order derivative of u with respect to t. By E~ we denote a subset of the space R°° with coordinates (3) in which system (1) is defined and admissible states and controls lie. The coordinates (3) are called canonical coordinates on E~. The above-introduced open sets O(I) can be understood as the subsets of E°°. Namely, the first coordinates of any point of this subset are bounded by the conditions defining O(I) and the remaining coordinates are arbitrary. Such subsets of E~ are also denoted by O(I).
On the set E ~, one introduces [4] the structure of a topological space and the usual differential-geometric notions: smooth functions, vector fields, differential forms, etc. In particular, sets of the form O(I) and their arbitrary unions are called open sets in E" . A smooth function on E~ is defined as an infinitely differentiable function depending on a finite (but arbitrary) number of coordinates (3). The algebra of smooth functions on E~ is denoted by F(E). A vector field on E~ is a differentiation of the algebra F(E).
The vector field
d ^ „„ .m ^ (i+1) d
DZ fj(t 'x' u +zz u
dt % "/dx] 1=0 du.1 )
on E~ is called the total derivative with respect to t on E". The Lie derivative of a function g e F(E) along D coincides with the time-derivative of the function g according to system (1). Phase curves of the vector field D coincide with graphs of solutions of system (1) in E°°. Therefore, the pair (E~,D) is chosen as a geometric model of system (1). It is called a diffiety (or an infinite prolongation) of system (1).
Let <pe F(E). By dCp we denote the 1-form dp-D(p)dt on E~, and by C:A(E) we denote the F(E) -module spanned by the 1-forms dCp, pe F(E). By H1 denote the F(E) -submodule of the module C:A(E) spanned by the 1-forms dCx1,...,dCxn. By definition, put Hk+1 = (ae Hk :Dae Hk} for k > 0.
The dimension of the space of covectors 0)e, ae Hk, is finite for arbitrary k > 0 and 6e E The dimension of some F(E) -submodule H of the module of 1-forms on E~ at the point 6e E~ is defined as the dimension of the space of covectors 0)e, ae H. A point 6e E~ is said to be Brunovsky-regular (or simply B-regular) if in a neighborhood of 6 system (1) is regular and the dimensions of the modules Hk and Hk + D(Hk) are constant for each k > 1.
Theorem 1. [3] In a neighborhood of a B-regular point for system (1) there exist p > 0 functions X1,.,Xp of t,xx,..,xn and m forms a\,..,am e H1 such that
(dCx,...,dcxp}u(D1 (a)|k = 1,.,m, j > 0} is a basis of the moduleC1A(E).
In the case p = 0, the set (a\,...,om} is called a B-basis or a linearizing system of differential
forms [3] at the B-regular point. The proof of Theorem 1 gives a procedure for finding a B-basis (see [3]).
A differential operator of the type
g0 + gD + g2D2 +... + gkDk, g0,g1,...,gk e F(E),
is called a C -differential operator of order k . A C -differential operator acts on ClA(E). A matrix whose entries are C -differential operators determines the operator on the set of columns of differential forms. This operator is called a matrix C -differential operator. The order of a matrix C -differential operator is the maximal order of its entries.
A matrix C -differential operator A is called invertible if there exists a matrix C -differential operator V such that the compositions V ° A and A ° V are the identity operators. The matrix of an invertible C -differential operator is square. The size of an invertible C -differential operator is its matrix size.
Theorem 2. [3] Let 6 be a B-regular point for system (1), p = 0, and (Q,...,am} be a B-basis at 6 . System (1) is flat at a neighborhood of 6 if and only if in a neighborhood of 6 there exist m functions h1,.,hm on E~ and an invertible matrix C -differential operator A such that
( i i \ ! \ dCh1 a
: =A : .
,dChm J {°m ,
The functions hl,...,hm form a flat output of system (1).
Let us set dCa = da-D(a) adt for ae A*(E). Necessary and sufficient flatness conditions based on properties of the operator [ dC ,A-1] ° A are obtained in [6]. Let us recall the basic notations from [6]. A map A of A* (E) to A* (E) such that
^(Q) = a0 AQ + a1 aDQ +... + ak aDkQ, at e C 1A(E), is called a CA -operator of order k and of grading 1 and is denoted by
a0 a1 + a aD + ... + ak ADk.
A matrix whose entries are CA -operators of order < k is called a matrix CA -operator of order < k. It is easily shown [6] that R = [dC ,A-1] ° A is a matrix CA -operator of grading 1 satisfying the conditions
(dC - R)(a>) = 0, (dC - R) ° (dC - R) = 0.
We say that this CA -operator R is a deformation of the differential dC determined by the flat output (hp..., hm) if A is the invertible C -differential operator from Theorem 2.
3. MAIN RESULTS
Let h = {h,...,hm} be a flat output of system (1), co = {a1,.,am} a B-basis. Consider the F(E) -modules
Gp = Span f ( e ){D' (tf) | i = 1,., m,j = 0,1,., p}, Fk = Span F (E {dC (Dsh) | i = 1,., m,s = 0,1,..., k},
where p > 0, k > 0.
An integer I is called the order of the flat output h with respect to the B-basis 0), if F0 e G{ but F0 e GI-1. Similarly, an integer L is called the order of the B-basis co with respect to the flat output h , if G0 e FL but G0 e FL-1.
Since the B-basis co consists of 1-forms from H1, we see that the order of a flat output of the form (2) is more than the integer k in (2).
We say that a B-regular point 0e E" is d-regular for the pair (a>, h) if in a neighborhood of 6 the dimension of the module Fk n Gp is constant for each p = 0,1,.,l and k = 0,1,.,L, where I and L are the integers such that F0 e G{ and G0 e FL .
Theorem 3. Suppose co is a B-basis at a B-regular point 6 for a flat system with 2 inputs, h is a flat output such that the point 6 is d-regular for the pair (tf,h), and kR is order of the deformation of
the differential dC determined by the flat output h ; then
1) there exists a flat output h of order < kR with respect to the B-basis co;
2) the order of the B-basis co with respect to the flat output h is not more than k1
4. INVERTIBLE DIFFERENTIAL OPERATORS
To prove the main result, we need a description of invertible C -differential operators. Denote by Gp the set of all C -differential operators of order < p whose matrices are of size
1xm. Clearly, Gp is the F(E) -module under the multiplication defined by the equality
(f V)(Q) = f V(Q), f e F(E),Ve Gp,
where Q is an m -column of 1-forms.
Let A be an invertible C -differential operator. Consider the F (E) -modules
Fk = {V o A|Ve GJ, k > 0. If A is the invertible C -differential operator from Theorem 2, then Gp = {V(a)|Ve Gp}, Fk = {V(a)|Ve Fk},
where a is the m -column of 1-forms a\,.,am forming a B-basis.
Each invertible C -differential operator determines the sequence of nonnegative integers dkp = dim(Fk n Gp) for p, k > 0. Further, we will classify invertible C -differential operators by the
sequences {dkp}.
A point 9e E" is said to be d-regular for an invertible C -differential operator A, if in a neighborhood of 9 the dimension of the module Fk nGp is constant for each p = 0 ,1 ,...,l and
k = — , L, where l = ord A, L = ord A-1.
The sequence {dk p} uniquely determines the sequence of integers
Pk,p = ^k,p - "k-1,p-1 > k, p > (4)
Where
Kk, p = dk,p - dk-1,p - dk,p-1 + dk-1,p-1
(put dk = Kk = 0, if k < 0 or p < 0). Conversely, the sequence {pkp} uniquely determines the sequence {dk p} .
Note that {dk p} is an increasing sequence for large values of k or p, whereas the sequence {Pk p} has only a finite number of non-zero values (see below).
To describe the sequences of integers (4), we consider a finite set of squares in the first quarter of the (k , p) -plane. Let corners of the squares have integer coordinates, and the sides of the squares be parallel to the coordinate axes. Let us assign the integer -1 to the top right corner of each square, and the integer 1 to the lower left corner. Note that the squares can intersect and even coincide. If there are points that are the corners of several squares, then their integers are added together. The zero value is assigned to all other integer points. We get the sequence of integers pk , k > 0, p > 0. Suppose there exists a sequence of
nonnegative integers ak such that the sequence {pk p = pk p + ak p} satisfies the following conditions:
S Pup = 0, S Pk,; = 0,
(5)
i=0 j=0
S P,0 = m S P0, ; = m
(6)
k-1 p-1
S P, p - Zk, p. S Pk, ; - Zk, p.
(7)
where k,p are arbitrary positive integers, m is some positive integer, zkp is the number of the squares whose the top right corners are situated at the point (k, p).
A set of squares satisfying the above mentioned conditions is called a d-scheme of squares (or simply a d-scheme). The corresponding sequence {pk p} is called the m -table of the d-scheme of squares.
Note that a d-scheme of squares can define more than one m -table. Besides, if we add an arbitrary positive integer a to the value p0 0, then conditions (5) and (7) hold true, and the integer m
in (6) increases by a . Thus we can go from an m -table to an (m + a) -table.
The following theorem generalizes a result of [5] to C -differential operators.
Theorem 4. a) For any invertible C -differential operator in a neighborhood of a d-regular point the sequence of integers (4) coincides with an m -table of some d-scheme.
b) If a d-scheme of squares has an m -table, then there are an invertible C -differential operator such that its sequence of integers (4) coincides with the given m -table.
The proof of Theorem 4 is similar to the proof of the corresponding theorem of [5] and is based on the algebraic theory of chain complexes and their spectral sequences.
Note that an invertible C -differential operator is not uniquely determined by its d-scheme of squares. Let us show how to construct an invertible C -differential operator for a given d-scheme and what structures still should be given for this construction. Suppose we have a d-scheme of squares and a sequence of nonnegative integers ak p such that the sequence {pk p =pk + ak p} satisfies
conditions (5)-(7). By Z denote the set of all upper right corners of squares of the d-scheme, and by Zk the set of all elements of Z with coordinates(k ,p). Besides, consider the set B of elements of
two types. The elements of the first type are the lower left corners of the squares. For any k,p > 0 there are ak elements of the second type with coordinates (k ,p). There are no other elements in B. It follows from equations (5)-(6) that the sum of the integers ak is equal to m. So the set B has m elements of the second type. Denote by Bk the set of all elements of B with coordinates (k , p).
To each element b e Bk assign a C -differential operator Vb e Gp satisfying the following conditions. The operators Vb,be Bk 0,k > 0 , are generators of the module G0. Each element ze Zk determines a relation of the form
k p
m q , 0 ,
(8)
i=1 j=1 fieBi ;
where Dz p is a scalar C -differential operator of order < min{k - i,p - j} satisfying the following properties:
(1) for any p > 0 using (8) for z e Zkp, k > 0, the operators Vb, b e Bk , k > 0, can be expressed as
p-i
Vb = HI Ab,ßoVß, (9)
1>0 j=0 ßeB, j
where Ab ß are scalar C -differential operators;
(2) for any k > 0 using (8) for z e Zkp, p > 0, the operators Vb, b e Bk , p > 0, can be expressed as
Vb = II I BbßßOVß, (10)
1 > 0 j=0 ßeBj j
where Bb ß are scalar C -differential operators.
Further, both coordinates of any element of Z are positive. Therefore, the integer pk0 coincides with the number of elements in Bk0. From the first equality in (6) it follows that the set B contains m elements whose the second coordinate is zero. Let b1, ..., bm denote these elements. Similarly, the set B contains m elements whose the first coordinate is zero. These elements are denoted by b11, ., blm. From (9) we get
Vb = I\ ° Vb' i =m-
j=i
The matrix operator A = (AjJ) is an invertible C -differential operator such that its sequence of integers (4) coincides with the given m -table. Similarly, from (10) we have
m
vbj = I° V,' j = ^m.
i=1
The matrix operator A-1 = (A--1) is the inversion of A.
Finally, if P, Q is invertible C -differential operators of order 0 , i.e., invertible matrices of functions, then the sequences {pk p} for the invertible C -differential operators A and P ° A ° Q
coincide. It can be proved that any invertible C -differential operator can be constructed in such a way. Thus an invertible C -differential operator is uniquely determined by its d-scheme, scalar C -differential operators Dz p from (8), and matrices P , Q . However this operator does not uniquely
determine operators Hz p and matrices P, Q .
Theorem 5. If for a d-scheme of squares there exists a 2-table, then the d-scheme has the form specified on Fig. 1, where s is the number of squares (s > 0), dt is the side
length of l -th square (dt > 0), (IS_Jdi, Il dt) are the coordinates of the element Zj e Z, (I* + dt,11-1 di) are the coordinates of the element bt e B,l = 1,...,s, (1*= di,0)
are the coordinates of the element b0 e B, (0, ^s dt) are the coordinates of the element »
Theorem 5 follows from the definition of 2 -tables of the d-schemes.
bs+16 B.
S+l
s
9
b •-------c>--------<
\ b( Fig. 1. d-Scheme with 2 -table
5. PROOF OF THEOREM 3
Consider a flat system with 2 inputs. Let c = {cx , c2} be a B-basis of the system, h = {h,, h2} a flat output. Denote by A the invertible C -differential operator from Theorem 2. From Theorem 4 it follows that the sequence {pkp} for the operator A-1 coincides with a 2 -table of some d-scheme.
From Theorem 5 it follows that this d-scheme has the form specified on Fig. 1. Relation (8) for the element z e Z of this d-scheme has the form
□, b., °V b., + □. b °v b +Q b. °v ^ = 0,
(11)
where HZ b+ and HZ b. are scalar C -differential operator of order 0 , i.e., functions, HZ b is a scalar C -differential operator of order < di. Denote
Vi = V b ' ai =0,^1' Q =-ai □zi.v ci =D,i,bi_i' i = 1
s.
From (9) it follows that a ^ 0. From (10) it follows that ci ^ 0. Therefore, relation (8) can be rewritten as
Vi+1 =q oVi -«Vi-1, i = 1,..., s,
(12)
where « = ci / ai ^ 0.
One has
r i+i
4 1
-a ^
J И i -1,
i = 1,...,
Denote by a, the matrix C -differential operator
f q -aл 1 0
у
(13)
This operator is invertible and its inversion is
A-1 =
0 1 1 1
4
V Ц Ц ,
By construction, V0 = (01), V1 = (10), and
iV 1
VT1 =A, о „.о A1 V 5 1
'v 5 ^
is the invertible C -differential operator such that its sequence {pkp} coincides with the 2-table. Whence, there exist invertible matrices P, Q of functions such that
A-1 = P ° As °... o A1 o Q.
Recall that operators P, Q, and Q are not uniquely determined. Let us show that we can choose P,Q, and Q in such a way that ai = 1 in (13) for i = 1,.,s. Namely, change P to P°O- , Q to O1 ° Q, and a. to Ot ° A i ° O~1, where
a =
У+1 0' , 0 У,
i = 1,..., 5,
Y1,., ys+1 are some nonzero numbers. Then the operator A 1 does not change and we have
O ° Аг ° O-1 =
У
-a
Y
У 7,-1
Puttingy0 = 1,Y = 1, and yi+1 = a for i = 1,.,s, we get the required operators. Below we consider the case a = 1, i = 1, • •, s.
We see that the order of operators A and A 1 is d2 — dt. Therefore the order of the flat output h with respect to the B-basis a and the order of the B-basis a with respect to the flat output h are d2. Denote
A — PoAs A2, A — A1 °Q.
Then
A-1 =A° A, c = (A° A)(dch)
and hence
R = [dc,A 1]°A = dCA°A 1 + A°dCA°A 1 °A 1
(14)
The order of operators A and A-1 is dz - d1. Therefore, the order of CA -operator dCA° A-1 is not more than 2(ds - d1).
Let us show that the order of R is determined by the second summand in (14). Let the matrix CA -operator dCQ ° Q-1 be
'qii ^12'
?
v q21 q22 ,
qu» q12 » q^ ^22e C'A( E
Then an easy calculation shows that
dC a ° A1 =
-q qX2 + q -Vn
a
V11+V12 P J'
where a —q q„ - ^21 +q q -^22 Q + dC Q.
Let q — g0 + +... + gDd' and let us consider 3 cases. Case 1: q12 ^ 0. The symbol of the B-basis a is
' 0 ql2 g2 D2d1 ^ v0 0 ,
>A-1,
and kc = 2ds.
Case 2: q12 = 0,q0 = gq11 - q22g + dCg ^ 0. The symbol of c is
0 qD 0 0
d1 \
>A-1,
and ka — 2ds - d1.
Case 3: gq11 - q22g + dCg — 0 . It can be proved that in this case there exists a flat output h of
order < d2.
Continuing this line of reasoning, we obtain the cases 1, 2 or a flat output of less order. In the cases 1 and 2, kR > d2. This proves Theorem 3.
REFERENCES
1. Fliess M., Levine J., Martin Ph. and Rouchon P. C.R. Acad. Sci. Pari5, 1992, I-315, pp. 619-624.
2. Martin, Ph., Murray R., and Rouchon P. Proc. of the 4th European Control Conf. Plenary lecture5 and Mini-cour5e5, Bru55el5, 1997, pp. 211-264.
3. Pomet J.-B., Moog C. and Aranda E. Proc. of the 3l5t IEEE Conf. on Deci5ion and Control, 1992, pp. 2012-2017.
4. Krasil'shchik I.S., Vinogradov A.M. (Eds.) Symmetries and Conservation Laws for Differential Equations of Mathematical Physics. Amer. Math. Soc., Providence (RI), 1999, 464 p.
5. Chetverikov V.N. Nauka I Obrzovanije v MGTU, 2014, 7, pp. 105-127 (DOI: 10.7463/0714.0718107)
6. Chetverikov V.N. Forum Math., 2004, v. 16, pp. 903-923.
ABOUT THE AUTHOR
Chetverikov Vladimir Nikolayevich, born 1958, graduated Moscow State University named after Lomonosov (1980), Sci.D. in physical and mathematical sciences, professor of the Department of Mathematical Modelling MSTU named after Bauman, author of 54 published scientific papers, sphere of interests — transformations of partial differential equations, symmetries, conservation laws, functional-differential equations, nonlinear dynamical control systems, dynamical feedback systems, flat systems.
УСЛОВИЕ ПЛОСКОСТИ ДЛЯ СИСТЕМ С ДВУМЯ ВХОДАМИ
В.Н. Четвериков
Изучаются плоские системы с двумя входами. Наш подход основан на обратимых дифференциальных операторах и деформации структур на диффеотопе систем с управлением. Описаны обратимые дифференциальные операторы размерности 2x2. Введено понятие веса для B-базиса. Если вес B-базиса нулевой, то проверка плоскости тривиальна. Минимальное количество плоских выходов оценивается на основании порядка соответствующего обратимого дифференциального оператора, линеаризующего систему с управлением. Минимальный порядок обратимого дифференциального оператора, линеаризующего систему, оценивается через порядок его деформации.
Ключевые слова: Нелинейные системы, бесконечное продолжение, плоские системы, линеаризация динамической обратной связи, обратимые дифференциальные операторы.