Systems & Control Letters

Size: px
Start display at page:

Download "Systems & Control Letters"

Transcription

1 Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: A converse to the eterministic separation principle Jochen Trumpf a,, Harry L. Trentelman b a esearch School of Engineering, Australian National University, anberra, Australia b Johann Bernoulli Institute for Mathematics an omputer Science, University of Groningen, The Netherlans article info abstract Article history: Available online xxxx Keywors: Separation principle Output feeback Linear systems Behavioral approach In the classical theory of finite-imensional linear time-invariant systems in state space form the term eterministic separation principle refers to the observation that a stabilizing output feeback controller can be constructe by first constructing an asymptotic state observer that is then couple to a stabilizing state feeback controller. In this paper we iscuss the following converse problem: an every stabilizing output feeback controller be realize as interconnection of an asymptotic state observer an a stabilizing state feeback controller? We will provie an affirmative answer to this question (moulo a number of technicalities) in a behavioral setting an with the help of rational representations. 216 Elsevier B.V. All rights reserve. 1. Introuction The classical eterministic separation principle says that, given the plant ẋ Ax + Bu, y x, an asymptotic full state observer ˆx (A G)ˆx + Bu + Gy (1) with A G Hurwitz, an a stabilizing static full state feeback controller u F ˆx (2) with A + BF Hurwitz, the close-loop ynamics is given by x A + BF BF x, (3) e A G e where e ˆx x is the observer error [1]. It follows from the form of the system matrix in (3) that lim t1 x(t), i.e. the observerbase output feeback controller (1) an (2) is stabilizing. In fact, it is even internally or totally stabilizing since x implies both y an ˆx (since also e ), an hence also u. orresponing author. aresses: Jochen.Trumpfanu.eu.au (J. Trumpf), h.l.trentelmanrug.nl (H.L. Trentelman) / 216 Elsevier B.V. All rights reserve. The above principle is calle a separation principle because it allows to complete the task of constructing an output feeback controller with esirable properties (namely stability) by separately constructing a state observer an a full state feeback controller with that property. Another classical but unrelate separation principle is that of optimal stochastic control, see e.g. [2]. An obvious converse question is whether there are any other constructions of (totally) stabilizing output feeback controllers, or whether any such controller permits an interpretation as a series connection of a full state observer followe by a (possibly ynamic) full state feeback controller. A partial answer to this question was given by Schumacher using the geometric notion of compensator couples at the beginning of the 198s [3], but a full answer remaine elusive to this ate. In this paper we aress the converse question in a behavioral framework using both polynomial an rational representations of linear ifferential systems. We show that, uner mil assumptions on the to be controlle system with variables (x, u, y), any controllable, regular, totally stabilizing controller through the variables (u, y) can be separate into an asymptotic i/o-observer for x from (u, y) with variables (ˆx, u, y) an a regular, totally stabilizing controller with variables (ˆx, u) in the sense that the controllable part of the observer/(ˆx, u)-controller interconnection coincies with the given (u, y)-controller. The paper is organize as follows. Section 2 introuces our notation an collects relevant results from the theory of behaviors incluing basics on rational representations. In Section 3 we review the require material on stabilization in a behavioral framework. In Section 4, we evelop a convenient system representation that

2 2 J. Trumpf, H.L. Trentelman / Systems & ontrol Letters ( ) is aapte to the problem treate in this paper. Section 5 contains the main result, Section 6 iscusses the special case of state space systems an Section 7 conclues the paper. 2. Behaviors of linear ifferential systems In this paper we will make heavy use of the mathematical machinery of the behavioral approach to linear ifferential systems. A linear ifferential system is efine as a triple (, w, B) whose behavior B 1 (, w ) is the solution space of a finite set of higher orer constant coefficient linear ifferential equations Polynomial an rational kernel representations Behaviors of linear ifferential systems can be represente in terms of a real polynomial matrix (s) with w columns as ( )w, so that B w 2 1 (, w ) w. (4) The representation (4) is calle a polynomial kernel representation of B, an we often write B ker ( ). If 1 an 2 are two full row rank polynomial matrices, then they represent the same behavior B, i.e. B ker 1 ( ) ker 2( ), if an only if there exists a unimoular polynomial matrix U such that 2 U 1. For an extensive treatment of polynomial representations of behaviors we refer to [4]. Behaviors also amit representations in terms of real rational matrices. A etaile exposition on rational representations can be foun in [5]. Here we will give a brief review. ecall that any given real rational matrix amits a left coprime factorization into polynomial matrices. A factorization of a real rational matrix as P 1 Q with P, Q real polynomial matrices is calle a left coprime factorization if P Q is left prime (meaning that it has a polynomial right inverse) an et(p) 6. Following [5], if P 1 Q is such a left coprime factorization then we efine w to be a solution of ( )w if it is a solution of the ifferential equation Q ( )w. In other wors, we efine ker : ker Q, (5) which is well-efine since any two left coprime factorizations of iffer by a unimoular polynomial factor. For a given rational matrix, we call a representation of B as ( )w arational kernel representation of B an write B ker ( ). For aitional material on rational representations we refer to [6,7]. In this paper we will often assume that the rational matrices (s) use in kernel representations have full row rank over the fiel of real rational functions. This is equivalent to saying that the kernel representation is minimal, see [5,6]. As note before, two minimal polynomial kernel representations iffer by a unimoular polynomial factor. A similar statement oes not hol for rational representations. We will come back to this in the next subsection ontrollability an controllable part In the behavioral approach an important role is playe by the property of controllability. The efinition of controllability of a behavior B is well known, an can be foun in [4]. ontrollability of a behavior can be teste in terms of its full row rank rational kernel representations as follows: If B ker ( ) where (s) is a rational matrix, then B is controllable if an only if has no zeros. A behavior B is calle autonomous if it is a finite imensional subspace of 1 (, w ). In terms of its rational kernel representations B ker ( ) this property requires that has full column rank. Any behavior B amits a irect sum ecomposition as B B cont B aut, where B cont, calle the controllable part of B, is the largest controllable subbehavior of B, an B aut, calle an autonomous part, is an autonomous subbehavior of B. The controllable part is uniquely etermine by B. In terms of its rational kernel representations ( )w, the controllable part of B can be foun by factorizing Q with Q nonsingular rational an a left prime polynomial matrix. For any such factorization we have B cont ker ( ), see [5]. It was shown in [6] that if 1 an 2 are full row rank rational matrices, then there exists a nonsingular rational matrix Q such that 2 Q 1 if an only if 1 an 2 represent behaviors with the same controllable part, i.e. (ker 1 ( )) cont (ker 2 ( )) cont Elimination of variables We will now review the basics of elimination of variables. Suppose we have a behavior B in which the manifest variable is partitione into two parts as w (v, c). Let v ( )v + c( )c be a polynomial kernel representation of B. The space of trajectories that satisfies this equation is calle the full behavior. The space of tractories c that are compatible with the equation of the full behavior is calle the behavior with v eliminate an is given by B c : c there exists v such that v v + c c. (6) The elimination problem is to obtain a kernel representation of (6). Such a kernel representation can be obtaine as follows: first fin a unimoular polynomial matrix U such that v,1 U v where v,1 has full row rank. Next, apply the same unimoular matrix to c to obtain c,1 U c. c,2 Since ker v c ( ) ker U v c ( ), a new, more structure, polynomial kernel representation of B is then given by v,1 c,1 v. c A kernel representation of the behavior (6) is now given by c,2 ( )c (see [4]). The above construction to obtain the eliminate behavior (6) is only vali for polynomial kernel representations an uses unimoular premultiplication. Its counterpart for the case that we eal with rational kernel representations an, instea of unimoular premultiplication, we use premultiplication with a nonsingular rational matrix is more subtle an is ealt with in the following lemma. Lemma 2.1. Let v ( )v + c( )c be a rational kernel representation of B. Let Q be a nonsingular rational matrix such that v,1 Q v

3 J. Trumpf, H.L. Trentelman / Systems & ontrol Letters ( ) 3 where v,1 is a full row rank polynomial matrix. Partition c,1 Q c c,2 an assume that c,1 is also a polynomial matrix. Then the controllable parts of ker c,2 ( ) an B c are equal. Proof. Define the rational matrix by v,1 : c,1. (7) Since Q ( v c ), we have that B cont (ker ( )) cont. Let c,2 P 1 2 Q 2 be a left coprime factorization. Then ker c,2 ( ) ker Q 2 ( ), an also 1 I v,1 c,1 P 2 Q 2 is a left coprime factorization, so we have v,1 ker ker c,1. Q 2 Since v,1 has full row rank, we have (ker ( )) c ker Q 2 ( ) ker c,2 ( ). This yiels (B c ) cont (B cont ) c ((ker ) cont ) c ((ker ) c ) cont ker c,2 cont, where we have omitte the ifferentiation symbol. Here we have use the fact that the operation of taking the controllable part of a behavior, an eliminating a variable from a behavior commute, see e.g. Lemma in [8]. We now stuy the relate question whether the triangular structure in (7) amits a similar triangular structure in a representation of the controllable part. More specifically, assume that the behavior B is represente by the rational kernel representation associate with (7), oes there exist a representation of its controllable part with compatible triangular structure. This inee turns out to be the case as is shown in the next lemma. Lemma 2.2. onsier the behavior B represente by the full row rank rational representation 1 v,1 c,1 B v A. (8) c Assume that v,1 an c,1 are polynomial, an v,1 has full row rank. Then there exist polynomial matrices v,1, c,1, c,2, an rational matrices Q 11, Q 12 an Q 22, with Q 11 an Q 22 square nonsingular, such that v,1 c,1 Q 11 Q 12 v,1 c,1 Q 22 an such that 1 v,1 c,1 B v A c is a representation of the controllable part B cont of B. Proof. First note that since v,1 an c,1 are polynomial matrices, an v,1 has full row rank, c,2 ( )c is a kernel representation of B c. Factorize v,1 c,1 Q11 Q 12 v,1 c,1 (9) Q 21 Q 22 v,2 c,2 where the first factor is nonsingular rational an the secon left prime polynomial. The secon factor yiels a kernel representation of B cont. We now eliminate v from B cont : by premultiplying the secon factor in (9) with a suitable unimoular polynomial matrix we obtain v,1 c,1 U11 U 12 v,1 c,1, (1) U 21 U 22 v,2 c,2 with v,1 full row rank. Then c,2 ( )v is a representation of (B cont ) c. By combining (9) an (1) we obtain v,1 c,1 Q 11 Q 12 v,1 c,1 (11) Q 21 Q 22 for suitable rational matrices Q ij. From (11) we obtain Q 21 v,1. Since v,1 has full row rank this yiels Q 21. Also we rea that c,2 Q 22 c,2. Since (B c ) cont (B cont ) c, the controllable part of ker c,2 ( ) is equal to ker c,2( ). Since both c,2 an c,2 have full row rank, this implies that Q 22 must be square. Then also Q 11 must be square, an both must be nonsingular Inputs an outputs Often, the manifest variable w of a behavior B is partitione as w col(u, y), an accoringly the rational kernel representation takes the form u ( )u + y( )y. This partitioning is calle an input output partition of the behavior if u is input an y is output, see [5]. In terms of the full row rank rational matrices this is equivalent to the property that y is a nonsingular rational matrix. In general a behavior amits many input output partitions. However, the sizes of u an y (enote by m(b) an p(b), respectively) are uniquely etermine by the behavior. It is well known that if ( )w is a kernel representation of B (either polynomial or rational), then the number of outputs p(b) of B is equal to rank. 3. Stabilization in a behavioral framework 3.1. Stabilization by full interconnection Given a behavior B (calle the plant) with polynomial kernel representation ( )w an a behavior (calle a controller) with polynomial kernel representation ( )w, the interconnection of B an is simply their intersection B \. learly, the interconnection has the kernel representation 1 B A w. The interconnection is calle a regular interconnection if p(b\) p(b) + p(). If both an have full row rank, then the interconnection is regular if an only if the matrix (12)

4 4 J. Trumpf, H.L. Trentelman / Systems & ontrol Letters ( ) has full row rank. The problem of stabilization by full interconnection is to fin, for B, a controller such that their interconnection is regular, an lim t1 w(t) for all w 2 B \. The controller is then calle a stabilizing controller for B. If both an have full row rank, then is a stabilizing controller for B if an only if the polynomial matrix (12) is Hurwitz, see [9,1]. Next we will prove that for each stabilizing controller for a given plant, also its controllable part is a stabilizing controller. Lemma 3.1. Let be a stabilizing controller for B. Then cont is a stabilizing controller for B. Proof. Let ( )w be a minimal polynomial kernel representation of B an let ( )w be a minimal polynomial kernel representation of. Then the polynomial matrix (12) is Hurwitz. amits a factorization Q, with Q nonsingular polynomial an where ( )w is a minimal representation of the controllable part cont of. Then from the fact that I Q it is easily seen that the matrix must be Hurwitz as well. The above can be extene to stable rational kernel representations for the controller. We o however ahere in this paper to polynomial kernel representations of the plant B. It is easily seen that if the plant B has full row rank polynomial kernel representation ( )w an the controller has full row rank rational kernel representation ( )w then is a stabilizing controller if an only if the rational matrix (12) is nonsingular an has all its zeros in. A rational matrix is calle a stable rational matrix if all its poles lie in. A square nonsingular stable rational matrix is calle miniphase if its inverse is again stable, equivalently if all its zeros lie in. If ( )w is a rational kernel representation of the controller with a full row rank stable rational matrix, then it is a stabilizing controller if an only if the rational matrix (12) is miniphase, see [5]. For a given plant B, a stabilizing controller exists if an only if it is stabilizable. A efinition of stabilizability of behaviors can be foun in [4]. If B is given by the full row rank rational kernel representation ( )w, then B is stabilizable if an only if ( ) has full row rank for all 2 +, see [5]. In the following lemma we will give a parametrization of all stabilizing controllers for a given plant B. Lemma 3.2. Let B be stabilizable, an let ( )w be a minimal polynomial kernel representation. Let be a stable rational matrix such that (13) is miniphase. Then a controller ker ( ) with full row rank rational is a stabilizing controller for B if an only if there exists a rational matrix F an a square nonsingular rational matrix G with all its zeros in such that F + G. Proof. Let F + G with F rational an G square nonsingular an all its zeros in. We have I, F G which is clearly square nonsingular an has all its zeros in. Hence ker ( ) is a stabilizing controller. onversely, let ker ( ) be a stabilizing controller, an full row rank. Then is square nonsingular an has all its zeros in. Define 1 Q11 Q 12 :. Q 22 Q 21 This prouct is clearly square nonsingular an has all its zeros in. Also, it yiels Q 11 + Q 12, which implies that (I Q 11 Q 12 ). This clearly implies Q 11 I an Q 12. Now efine F : Q 21, an G : Q 22. Then G is square nonsingular an has all its zeros in an, finally, F + G Stabilization by partial interconnection In general, the plant B has two types of variables, the variable w to be controlle, an the variable c, calle the interconnection variable, through which the plant can be interconnecte with a controller. More specific, if the plant B has the polynomial kernel representation w ( )w + c( )c, an if a controller has the polynomial kernel representation ( )c, then the interconnection of B an through c is efine as the behavior B ^c : {(w, c) (w, c) 2 B an c 2 }. As before, the interconnection is calle regular if p(b) + p() p(b ^c ). In this paper we will consier the problem of total stabilization for the plant B. A controller is calle a totally stabilizing controller for B if the interconnection is regular an if lim t1 (w(t), c(t)) for all (w, c) 2 B ^c. If both plant an controller are represente by a minimal polynomial kernel representation, then is a totally stabilizing controller for B if an only if w c (14) is Hurwitz. Again, the above can be extene to stable rational kernel representations for the controller. We o however again ahere to polynomial kernel representations of the plant B. If ( )w is a rational kernel representation of the controller with a full row rank stable rational matrix, then it is a totally stabilizing controller if an only if the rational matrix (14) is miniphase. If is a full row rank rational matrix (not necessarily stable) then it is a totally stabilizing controller if an only if (14) is a square nonsingular rational matrix with all its zeros in. We will now eal with the question uner what conitions there exists a totally stabilizing controller for a given plant B. These conitions involve stabilizability an etectability. Given the plant behavior B with variable (w, c), we say that w is etectable from c if (w, ) 2 B implies lim t1 w(t). If B is represente by the polynomial kernel representation w ( )w + c( )c, then w is etectable from c if an only if w ( ) has full column rank for all 2 +, see [4]. The following lemma gives necessary an sufficient conitions on the plant B for the existence of a totally stabilizing controller. Lemma 3.3. Given a plant B with variable (w, c), there exists a totally stabilizing controller if an only if B is stabilizable an w is etectable from c in B.

5 J. Trumpf, H.L. Trentelman / Systems & ontrol Letters ( ) 5 Proof. Let w ( )w + c( )c be a minimal polynomial kernel representation of B. By unimoular premultiplication we can obtain an alternative kernel representation in the triangular form 1 w,1 c,1 B w A, (15) c with w,1 full row rank. Using the assumption that w is etectable from c we fin that, in fact, w,1 ( ) has full column rank for all 2 +. As a consequence, w,1 is square an Hurwitz. By stabilizability of B, c,2 ( ) has full row rank for all 2 +. Thus, the behavior represente by c,2 ( )c itself is stabilizable, so there exists a full row rank polynomial matrix such that c,2 is Hurwitz. learly then w,1 c,1 is Hurwitz, so the controller ker ( ) totally stabilizes B. onversely, assume that there exists such that w c (16) is Hurwitz. Then obviously ( w ( ) c ( )) has full row rank for all 2 +, so B is stabilizable. Now, assume that (w, ) 2 B. Then, clearly, also (w, ) 2 B ^c. Since the controller is totally stabilizing we must have w(t) as t 1, so in B, w is etectable from c. 4. A convenient system representation As announce in Section 1, in this paper we will consier linear ifferential systems P with system variable (x, u, y), where x is interprete as the variable to be controlle, an (u, y) as the interconnection variable through which P can be interconnecte with a controller. We will make the following assumptions on P. (A1) P is stabilizable, (A2) x is etectable from (u, y), (A3) u is input with (x, y) output, (A4) im(y) apple im(x). We will now briefly iscuss the above assumptions. First note that by Lemma 3.3, assumptions (A1) an (A2) are necessary an sufficient for the existence of a totally stabilizing controller for P. In Section 2.3 it was explaine that starting from a minimal polynomial kernel representation x x + u u + y y of P, we can always obtain a new, more structure, minimal polynomial kernel representation for P of the form x x,1 u,1 y,1 u, (17) y in which x,1 has full row rank. We will stuy how our assumptions are reflecte in the polynomial matrices in the representation (17). First note that (A2) hols if an only if x,1 has full column rank for all 2 +. Since x,1 has also full row rank we therefore have that assumption (A2) is equivalent with x,1 being Hurwitz. Next, assumption (A3) is equivalent to the conition that the submatrix x,1 y,1 y,2 is square an nonsingular. In particular this implies that y,2 is square an nonsingular. Finally, as a consequence of assumption (A4) we have that the matrix y,1 is tall, meaning that its number of columns oes not excee its number of rows. By applying Lemma 2.2, there exist polynomial matrices x,1, u,1, y,1, u,2, y,2 an rational matrices Q 11, Q 12 an Q 22, with Q 11 an Q 22 square nonsingular, such that x,1 u,1 y,1 Q 11 Q 12 x,1 u,1 y,1 Q 22 an such that v x,1 u,1 y,1 u (18) y is a representation of the controllable part P cont of P. Obviously, uner the assumptions mae above x,1 an y,2 are both nonsingular an y,1 is tall. emark 4.1. eferring to the iscussion at the beginning of Section 3.2, every regular, totally stabilizing controller for P through (u, y) has a polynomial (an hence stable rational) representation ker u y such that x,1 u,1 y,1 u y is Hurwitz (an hence miniphase). Since 1 1 x,1 u,1 y,1 Q 11 Q 12 Q 22 A x,1 u,1 y,1 A u y I u y this implies that u,2 u,2 u y has full row rank. We will use this fact in the proof of Proposition A converse to the separation principle The following theorem is the main result of this paper. It provies a converse to the classical eterministic separation principle. Theorem 5.1. onsier a system with variables (x, u, y) such that Assumptions (A1) (A4) from Section 4 hol. Then for every controllable, regular, totally stabilizing controller through (u, y) there exists an asymptotic i/o-observer O for x from (u, y) with variables (ˆx, u, y) an a regular, totally stabilizing controller K for the plant system with variables (ˆx, u) such that O ^(ˆx,u) K (u,y),cont. (19) A few remarks are in orer before we prove this theorem. Fig. 1 illustrates the separation of the controller into the two blocks O an K. Assumptions (A1) an (A2) are clearly necessary for totally stabilizing controllers through (u, y) to exist (Lemma 3.3). Assumptions (A3) an (A4) will turn out to be sufficient conitions

6 6 J. Trumpf, H.L. Trentelman / Systems & ontrol Letters ( ) Fig. 1. The separate controller of Theorem 5.1. for solvability of the matrix equation that is equivalent to (19), cf. Proposition 5.2. Ieally, we woul like to rop the assumption of controllability for the controller an the corresponing restriction of the theorem statement to the controllable part of the separate controller. This assumption an restriction are, however, necessitate by our use of rational representations. The general theory will likely require more avance algebraic methos an our inability to generalize the current proof using stanar polynomial or rational methos might help to explain why such a result was never obtaine in the classical state space literature (cf. [3]). Note, however, that if a controller is totally stabilizing then its controllable part is also totally stabilizing (Lemma 3.1). In orer to prove Theorem 5.1, we first erive a matrix equation formulation of conition (19). Proposition 5.2. onsier a system P with variables (x, u, y) an such that assumptions (A1) an (A2) from Section 4 hol. epresent P as in (17). onsier a controllable, regular, totally stabilizing controller ker u y through (u, y), where the polynomial representation of has been chosen as in emark 4.1. Then there exists an asymptotic i/o-observer O for x from (u, y) with variables (ˆx, u, y) an a regular, totally stabilizing controller K for P with variables (ˆx, u) such that (19) hols if an only if there exists a rational matrix Y an a stable rational matrix T such that y,1 Y YT y. (2) y,2 Proof. The proof of Proposition 5.2 procees in four major steps. In a first step we parametrize all observers O(S 1, S 2 ) for x from (u, y) using the authors recent internal moel principle for observers [11] an the algebraic generalization thereof [12]. Here, S 1 an S 2 are matrix parameters. In a secon step, we use Lemma 3.2 to provie a parametrization K(Y 1, Y 2, X) of all fully interconnecte regular stabilizing controllers for P where Y 1, Y 2 an X are matrix parameters. In a thir step we use Lemma 2.1 to show that the separate controller resulting from the interconnection of O(S 1, S 2 ) an K(Y 1, Y 2, X) has the same controllable part as if an only if we choose Y 2 an X appropriately. In a fourth an last step we characterize when the parameter Y 1 can be chosen such that K(Y 1, Y 2, X) is a controller with variables (ˆx, u) as require. Step 1. Since the system P is stabilizable, its anti-stabilizable part is equal to its controllable part, cf. Theorem 2.3 in [11]. But then O ker ˆˆx ˆu ˆy is a full row rank polynomial representation of an asymptotic i/o-observer for x from (u, y) if an only if ˆˆx is Hurwitz an P cont O, cf. Theorem 5.6 in [11] an orollary 16 in [12]; equivalently, ˆˆx ˆu ˆy S 1 S 2 x,1 u,1 y,1 S 1 u,1 + S 2 u,2 S 1 y,1 + S 2 y,2 with S 1 polynomial an Hurwitz, an S 2 polynomial. We write O(S 1, S 2 ) to inicate the epenence of O on the matrix parameters S 1 an S 2. Step 2. By Lemma 3.2, K ker K x K u K y is a fully interconnecte, regular, stabilizing controller for P if an only if Q K x K u K y Z 1 Z 11 Q 12 2 Q 22 x,1 u,1 y,1 + X u y, where Z 1 an Z 2 are rational an X is rational an non-singular square with only stable zeros. Define Y 1 Y 2 : Z 1 Z 2 Q 11 Q 12 Q 22 x,1 I an observe that the reparametrization Z 1 Z 2 $ Y 1 Y 2 is bijective. Hence K K(Y 1, Y 2, X) ker K x K u K y is a fully interconnecte, regular, stabilizing controller for P if an only if K x K u K y Y 1 Y 1 1 x,1 u,1 + Y 2 u,2 + X u Y 1 1 x,1 y,1 + Y 2 y,2 + X y (21) with Y 1 an Y 2 rational an X rational an non-singular square with only stable zeros. Step 3. ompute I Y I S 1 u,1 + S 2 u,2 S 1 y,1 + S 2 y,2 Y 1 Y 1 1 x,1 u,1 + Y 2 u,2 + X u Y 1 1 x,1 y,1 + Y 2 y,2 + X y S 1 u,1 + S 2 u,2 S 1 y,1 + S 2 y,2 (Y 2 Y S 2) u,2 + X u (Y 2 Y S 2) y,2 + X y By Lemma 2.1, (O(S 1, S 2 ) \ K(Y 1, Y 2, X)) (u,y),cont ker (Y 2 Y S 2) u,2 + X u (Y 2 Y S 2) y,2 + X y an since u,2 y,2 u y has full row rank, (O(S 1, S 2 ) \ K(Y 1, Y 2, X)) (u,y),cont if an only if Y 2 Y S 2 an X is rational an non-singular square with only stable zeros. Step 4. It only remains to characterize when K(Y 1, Y 1 1 x,1 S 1 1 S 2, X) has variables (ˆx, u) only, i.e. when K y Y 1 1 x,1 y,1 + Y 2 y,2 + X y Y 1 1 x,1 y,1 + Y S 2 y,2 + X y. Setting Y X 1 Y 1 1 x,1 an T S 1 1 S 2 shows that this is equivalent to Eq. (2). We can now complete the proof of Theorem 5.1 by showing that, uner the conitions of the theorem, Eq. (2) has a solution. Proof of Theorem 5.1. Assumption (A3) implies that y,2 is invertible. By Assumption (A4), the matrix y,1 1 y,2 is tall an hence there exists a left-invertible rational matrix B such that T : B y,1 1 y,2

7 J. Trumpf, H.L. Trentelman / Systems & ontrol Letters ( ) 7 is stable rational: let p 1 1 q U 1 y,1 1 y,2 V p r B q r A be in Smith McMillan form then B can be chosen as p 1 1 q B : U p r V 1, B q r A I making T polynomial an hence stable rational. Let B L be a leftinverse of B, i.e. B L B I. Let Y : y 1 Y YT as require. y,1 y,2 y,2 BL then y 1 y,2 BL y,1 y 1 y,2 BL (B y,1 1 y,2 ) y,2 y In the above theorem, the only constraints place on the controller K are that it is itself regular an totally stabilizing, an that it connects to the plant only through (x, u); K is written here with variables (ˆx, u) since it is intene to be combine with the observer O. In particular, it is not assume a priori that the controller K acts by feeback or that it is even static in the classical sense. The latter woul amount to the aitional conitions Kˆx I an K u F where F is a constant matrix, yieling the controller representation u F ˆx. It can be shown that Eq. (2) oes not always have a solution uner these aitional conitions. 6. Example: the state space case eturning to the state space case that we alreay briefly iscusse in Section 1, assume that our plant P with system variable (x, u, y) is represente by ẋ Ax + Bu, y x, where x(t) 2 n, u(t) 2 m an y(t) 2 p. For simplicity, assume that the pair (A, B) is controllable an that the pair (, A) is etectable. Let be a regular totally stabilizing controller through (u, y) an assume that is controllable. Let have polynomial kernel representation u u + y y. In this section we will ientify an asymptotic i/o-observer O for x through (u, y) an a regular totally stabilizing controller K through (x, u) such that the controllable part of their interconnection equals. learly, a minimal polynomial kernel representation of P is given by x I A B u (22) I y an this reveals that, uner the above assumptions on the matrix pairs (A, B) an (, A), the behavior P is controllable (so also stabilizable) an that x is etectable from (u, y) in P. Obviously, in P, u is input an (x, y) is output, so the assumptions (A1), (A2) an (A3) of Section 4 hol. Assume that in aition we have p apple n so that also (A4) hols. Since P is controllable it is equal to its controllable part P cont. We will now first bring the kernel representation (22) in the require upper iagonal form. Let (si A) 1 L 1 2 L 1 be a polynomial left coprime factorization, an let N 1 an N 2 be polynomial matrices such that N1 N 2 L 1 L 2 is unimoular. Pre-multiplying (22) by this unimoular matrix, we see that an upper iagonal polynomial kernel representation of P P cont is given by N 1 ( x I A) + N 2 N 1 B N 2 u. L 1 B L 2 y Denote : N 1 (si A)+N 2 an note that this polynomial matrix is Hurwitz. In orer to ientify a suitable observer/controller pair we first solve the nonlinear equation (2), that in this case takes the form N2 Y YT y. L 2 Since assumptions (A1) (A4) hol, this equation inee has a solution pair (Y, T) with Y rational an T stable rational (cf. the Proof of Theorem 5.1). Next, factorize T S 1 2 S 1, with S 1 polynomial an Hurwitz, an S 2 polynomial. This yiels an asymptotic i/o-observer O for x from (u, y) with polynomial kernel representation ˆx S 1 (S 1 N 1 + S 2 L 1 )B S 1 N 2 + S 2 L 2 u. y It is inee easily verifie that if (x, u, y) 2 P an (ˆx, u, y) 2 O, then the error e : ˆx x satisfies S 1 e, so e(t) as t 1. Next we ientify a regular totally stabilizing controller K through (x, u). In general, all totally stabilizing controllers for P are parametrize by (21), with Y 1, Y 2 rational an X nonsingular rational with only stable zeros. Here we take X : I m, Y 1 : Y, Y 2 : YT. This yiels the regular totally stabilizing controller K with rational kernel representation x Y Y(N 1 + TL 1 )B + u. u Now interconnect the observer O S 1 (S 1 N 1 + S 2 L 1 )B S 1 N 2 + S 2 L 2 ˆx u y with K ˆx Y Y(N 1 + TL 1 )B + u u

8 8 J. Trumpf, H.L. Trentelman / Systems & ontrol Letters ( ) through the variables (ˆx, u). Then the controllable part of the (u, y) behavior of this interconnection is equal to the given controller as esire. 7. onclusions We have shown that uner the mil assumptions (A1) (A4) from Section 4 the converse of the eterministic separation principle hols true, namely that for a linear ifferential system with variables (x, u, y) any controllable, regular, totally stabilizing controller through the variables (u, y) can be separate into an asymptotic i/o-observer for x from (u, y) with variables (ˆx, u, y) an a regular, totally stabilizing controller with variables (ˆx, u) in the sense that the controllable part of the observer/(ˆx, u)-controller interconnection coincies with the given (u, y)-controller. Together with the fact that the controllable part of a totally stabilizing controller is also totally stabilizing an combine with the internal moel principle for observers an the structure of the interconnection, this result implies that any such controller must contain an internal moel of the controllable part of the to be controlle system. This observation shoul have interesting consequences for robustly stabilizing output feeback controllers. A characterization of those (u, y)-controllers that can be separate into an asymptotic observer an a feeback (ˆx, u)- controller is the topic of future work as is a generalization to not necessarily controllable (u, y)-controllers an the case im(y) > im(x). Acknowlegment This work was supporte in part by the Australian esearch ouncil through the A Discovery Project DP Geometric observer theory for mechanical control systems. eferences [1] H.L. Trentelman, P. Antsaklis, Observer-base control, in: Encyclopeia of Systems an ontrol, Springer, 214, p. 6. (no pagination). [2] T. Georgiou, A. Linquist, The separation principle in stochastic control, eux, IEEE Trans. Automat. ontrol 58 (1) (213) [3] J.M. Schumacher, Dynamic feeback in finite- an infinite-imensional linear systems (Ph.D. thesis), Vrije Universiteit te Amsteram, [4] J.. Willems, J.W. Polerman, Introuction to Mathematical Systems Theory: A Behavioral Approach, Springer, [5] J.. Willems, Y. Yamamoto, Behaviors efine by rational functions, Linear Algebra Appl. 425 (2) (27) [6] S.V. Gottimukkala, S. Fiaz, H.L. Trentelman, Equivalence of rational representations of behaviors, Systems ontrol Lett. 6 (2) (211) [7] S.V. Gottimukkala, H.L. Trentelman, S. Fiaz, ealization an elimination in rational representations of behaviors, Systems ontrol Lett. 62 (8) (213) [8] M.N. Belur, ontrol in a behavioral context (Ph.D. thesis), ijksuniversiteit Groningen, 23. [9] J.. Willems, On interconnections, control, an feeback, IEEE Trans. Automat. ontrol 42 (3) (1997) [1] M.N. Belur, H.L. Trentelman, Stabilization, pole placement, an regular implementability, IEEE Trans. Automat. ontrol 47 (5) (22) [11] J. Trumpf, H.L. Trentelman, J.. Willems, Internal moel principles for observers, IEEE Trans. Automat. ontrol 59 (214) [12] I. Blumthaler, J. Trumpf, A new parametrization of linear observers, IEEE Trans. Automat. ontrol 59 (214)

The canonical controllers and regular interconnection

The canonical controllers and regular interconnection Systems & Control Letters ( www.elsevier.com/locate/sysconle The canonical controllers an regular interconnection A.A. Julius a,, J.C. Willems b, M.N. Belur c, H.L. Trentelman a Department of Applie Mathematics,

More information

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and This article appeare in a journal publishe by Elsevier. The attache copy is furnishe to the author for internal non-commercial research an eucation use, incluing for instruction at the authors institution

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

POLYNOMIAL EMBEDDING ALGORITHMS FOR CONTROLLERS IN A BEHAVIORAL FRAMEWORK

POLYNOMIAL EMBEDDING ALGORITHMS FOR CONTROLLERS IN A BEHAVIORAL FRAMEWORK 1 POLYNOMIAL EMBEDDING ALGORITHMS FOR CONTROLLERS IN A BEHAVIORAL FRAMEWORK H.L. Trentelman, Member, IEEE, R. Zavala Yoe, and C. Praagman Abstract In this paper we will establish polynomial algorithms

More information

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition

More information

Discrete Mathematics

Discrete Mathematics Discrete Mathematics 309 (009) 86 869 Contents lists available at ScienceDirect Discrete Mathematics journal homepage: wwwelseviercom/locate/isc Profile vectors in the lattice of subspaces Dániel Gerbner

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable

More information

Linear Algebra- Review And Beyond. Lecture 3

Linear Algebra- Review And Beyond. Lecture 3 Linear Algebra- Review An Beyon Lecture 3 This lecture gives a wie range of materials relate to matrix. Matrix is the core of linear algebra, an it s useful in many other fiels. 1 Matrix Matrix is the

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

u!i = a T u = 0. Then S satisfies

u!i = a T u = 0. Then S satisfies Deterministic Conitions for Subspace Ientifiability from Incomplete Sampling Daniel L Pimentel-Alarcón, Nigel Boston, Robert D Nowak University of Wisconsin-Maison Abstract Consier an r-imensional subspace

More information

Robustness and Perturbations of Minimal Bases

Robustness and Perturbations of Minimal Bases Robustness an Perturbations of Minimal Bases Paul Van Dooren an Froilán M Dopico December 9, 2016 Abstract Polynomial minimal bases of rational vector subspaces are a classical concept that plays an important

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Permanent vs. Determinant

Permanent vs. Determinant Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

A. Incorrect! The letter t does not appear in the expression of the given integral

A. Incorrect! The letter t does not appear in the expression of the given integral AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question

More information

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan

More information

Summary: Differentiation

Summary: Differentiation Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

Separation Principle for a Class of Nonlinear Feedback Systems Augmented with Observers

Separation Principle for a Class of Nonlinear Feedback Systems Augmented with Observers Proceeings of the 17th Worl Congress The International Feeration of Automatic Control Separation Principle for a Class of Nonlinear Feeback Systems Augmente with Observers A. Shiriaev, R. Johansson A.

More information

Stabilization, Pole Placement, and Regular Implementability

Stabilization, Pole Placement, and Regular Implementability IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 47, NO. 5, MAY 2002 735 Stabilization, Pole Placement, and Regular Implementability Madhu N. Belur and H. L. Trentelman, Senior Member, IEEE Abstract In this

More information

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial

More information

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica

More information

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang

More information

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 5. Levi-Civita connection From now on we are intereste in connections on the tangent bunle T X of a Riemanninam manifol (X, g). Out main result will be a construction

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables* 51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*

More information

Pure Further Mathematics 1. Revision Notes

Pure Further Mathematics 1. Revision Notes Pure Further Mathematics Revision Notes June 20 2 FP JUNE 20 SDB Further Pure Complex Numbers... 3 Definitions an arithmetical operations... 3 Complex conjugate... 3 Properties... 3 Complex number plane,

More information

Interconnected Systems of Fliess Operators

Interconnected Systems of Fliess Operators Interconnecte Systems of Fliess Operators W. Steven Gray Yaqin Li Department of Electrical an Computer Engineering Ol Dominion University Norfolk, Virginia 23529 USA Abstract Given two analytic nonlinear

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Integration Review. May 11, 2013

Integration Review. May 11, 2013 Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In

More information

A nonlinear inverse problem of the Korteweg-de Vries equation

A nonlinear inverse problem of the Korteweg-de Vries equation Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:

More information

Simultaneous Input and State Estimation with a Delay

Simultaneous Input and State Estimation with a Delay 15 IEEE 5th Annual Conference on Decision an Control (CDC) December 15-18, 15. Osaa, Japan Simultaneous Input an State Estimation with a Delay Sze Zheng Yong a Minghui Zhu b Emilio Frazzoli a Abstract

More information

Journal of Algebra. A class of projectively full ideals in two-dimensional Muhly local domains

Journal of Algebra. A class of projectively full ideals in two-dimensional Muhly local domains Journal of Algebra 32 2009 903 9 Contents lists available at ScienceDirect Journal of Algebra wwwelseviercom/locate/jalgebra A class of projectively full ieals in two-imensional Muhly local omains aymon

More information

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics

More information

DISTANCE BETWEEN BEHAVIORS AND RATIONAL REPRESENTATIONS

DISTANCE BETWEEN BEHAVIORS AND RATIONAL REPRESENTATIONS DISTANCE BETWEEN BEHAVIORS AND RATIONAL REPRESENTATIONS H.L. TRENTELMAN AND S.V. GOTTIMUKKALA Abstract. In this paper we study notions of distance between behaviors of linear differential systems. We introduce

More information

Periods of quadratic twists of elliptic curves

Periods of quadratic twists of elliptic curves Perios of quaratic twists of elliptic curves Vivek Pal with an appenix by Amo Agashe Abstract In this paper we prove a relation between the perio of an elliptic curve an the perio of its real an imaginary

More information

Left-invariant extended Kalman filter and attitude estimation

Left-invariant extended Kalman filter and attitude estimation Left-invariant extene Kalman filter an attitue estimation Silvere Bonnabel Abstract We consier a left-invariant ynamics on a Lie group. One way to efine riving an observation noises is to make them preserve

More information

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies Laplacian Cooperative Attitue Control of Multiple Rigi Boies Dimos V. Dimarogonas, Panagiotis Tsiotras an Kostas J. Kyriakopoulos Abstract Motivate by the fact that linear controllers can stabilize the

More information

BLOW-UP FORMULAS FOR ( 2)-SPHERES

BLOW-UP FORMULAS FOR ( 2)-SPHERES BLOW-UP FORMULAS FOR 2)-SPHERES ROGIER BRUSSEE In this note we give a universal formula for the evaluation of the Donalson polynomials on 2)-spheres, i.e. smooth spheres of selfintersection 2. Note that

More information

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign

More information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On the number of isolated eigenvalues of a pair of particles in a quantum wire On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in

More information

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

A new approach to explicit MPC using self-optimizing control

A new approach to explicit MPC using self-optimizing control 28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 28 WeA3.2 A new approach to explicit MPC using self-optimizing control Henrik Manum, Sriharakumar Narasimhan an Sigur

More information

Practical Computation of Flat Outputs for Nonlinear Control Systems

Practical Computation of Flat Outputs for Nonlinear Control Systems IOP onference Series: Materials Science an Engineering PPER OPEN ESS Practical omputation of Flat Outputs for Nonlinear ontrol Systems To cite this article: Joe Imae et al 25 IOP onf Ser: Mater Sci Eng

More information

Centrum voor Wiskunde en Informatica

Centrum voor Wiskunde en Informatica Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.

More information

Mathematical Review Problems

Mathematical Review Problems Fall 6 Louis Scuiero Mathematical Review Problems I. Polynomial Equations an Graphs (Barrante--Chap. ). First egree equation an graph y f() x mx b where m is the slope of the line an b is the line's intercept

More information

Nested Saturation with Guaranteed Real Poles 1

Nested Saturation with Guaranteed Real Poles 1 Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Dynamics of the synchronous machine

Dynamics of the synchronous machine ELEC0047 - Power system ynamics, control an stability Dynamics of the synchronous machine Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct These slies follow those presente in course

More information

APPPHYS 217 Thursday 8 April 2010

APPPHYS 217 Thursday 8 April 2010 APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Tractability results for weighted Banach spaces of smooth functions

Tractability results for weighted Banach spaces of smooth functions Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Dissipative numerical methods for the Hunter-Saxton equation

Dissipative numerical methods for the Hunter-Saxton equation Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a

More information

The effect of dissipation on solutions of the complex KdV equation

The effect of dissipation on solutions of the complex KdV equation Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,

More information

Polynomial Inclusion Functions

Polynomial Inclusion Functions Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl

More information

The Non-abelian Hodge Correspondence for Non-Compact Curves

The Non-abelian Hodge Correspondence for Non-Compact Curves 1 Section 1 Setup The Non-abelian Hoge Corresponence for Non-Compact Curves Chris Elliott May 8, 2011 1 Setup In this talk I will escribe the non-abelian Hoge theory of a non-compact curve. This was worke

More information

Minimum-time constrained velocity planning

Minimum-time constrained velocity planning 7th Meiterranean Conference on Control & Automation Makeonia Palace, Thessaloniki, Greece June 4-6, 9 Minimum-time constraine velocity planning Gabriele Lini, Luca Consolini, Aurelio Piazzi Università

More information

On Decentralized Optimal Control and Information Structures

On Decentralized Optimal Control and Information Structures 2008 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June 11-13, 2008 FrC053 On Decentralize Optimal Control an Information Structures Naer Motee 1, Ali Jababaie 1 an Bassam

More information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How the potentials in different gauges yield the same retarded electric and magnetic fields How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department

More information

SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS

SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS GEORGE A HAGEDORN AND CAROLINE LASSER Abstract We investigate the iterate Kronecker prouct of a square matrix with itself an prove an invariance

More information

arxiv: v1 [math-ph] 5 May 2014

arxiv: v1 [math-ph] 5 May 2014 DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic

More information

Optimal Control of Spatially Distributed Systems

Optimal Control of Spatially Distributed Systems Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

Lecture 6: Calculus. In Song Kim. September 7, 2011

Lecture 6: Calculus. In Song Kim. September 7, 2011 Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear

More information

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device Close an Open Loop Optimal Control of Buffer an Energy of a Wireless Device V. S. Borkar School of Technology an Computer Science TIFR, umbai, Inia. borkar@tifr.res.in A. A. Kherani B. J. Prabhu INRIA

More information

ON BEAUVILLE STRUCTURES FOR PSL(2, q)

ON BEAUVILLE STRUCTURES FOR PSL(2, q) ON BEAUVILLE STRUCTURES FOR PSL(, q) SHELLY GARION Abstract. We characterize Beauville surfaces of unmixe type with group either PSL(, p e ) or PGL(, p e ), thus extening previous results of Bauer, Catanese

More information

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp . ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay

More information

Approximate reduction of dynamic systems

Approximate reduction of dynamic systems Systems & Control Letters 57 2008 538 545 www.elsevier.com/locate/sysconle Approximate reuction of ynamic systems Paulo Tabuaa a,, Aaron D. Ames b, Agung Julius c, George J. Pappas c a Department of Electrical

More information

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH

TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH English NUMERICAL MATHEMATICS Vol14, No1 Series A Journal of Chinese Universities Feb 2005 TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH He Ming( Λ) Michael K Ng(Ξ ) Abstract We

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

International Journal of Pure and Applied Mathematics Volume 35 No , ON PYTHAGOREAN QUADRUPLES Edray Goins 1, Alain Togbé 2

International Journal of Pure and Applied Mathematics Volume 35 No , ON PYTHAGOREAN QUADRUPLES Edray Goins 1, Alain Togbé 2 International Journal of Pure an Applie Mathematics Volume 35 No. 3 007, 365-374 ON PYTHAGOREAN QUADRUPLES Eray Goins 1, Alain Togbé 1 Department of Mathematics Purue University 150 North University Street,

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

Commun Nonlinear Sci Numer Simulat

Commun Nonlinear Sci Numer Simulat Commun Nonlinear Sci Numer Simulat 14 (2009) 3901 3913 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Dynamics of a ring network

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite

More information