The canonical controllers and regular interconnection

Size: px
Start display at page:

Download "The canonical controllers and regular interconnection"

Transcription

1 Systems & Control Letters ( 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, University of Twente, P.O. Box 217, 7500 AE Enschee, The Netherlans b K.U. Leuven, Kasteelpark Arenberg 10, B-3001 Leuven-Heverlee, Belgium c Department of Electrical Engineering, Inian Institute of Technology Bombay, Powai, Mumbai , Inia Department of Mathematics, University of Groningen, P.O. Box 800, 9700 AV Groningen, The Netherlans Receive 28 April 2004; receive in revise form 30 November 2004; accepte 1 December 2004 Abstract We stuy the control problem from the point of view of the behavioral systems theory. Two controller constructions, calle canonical controllers, are introuce. We prove that for linear time-invariant behaviors, the canonical controllers implement the esire behavior if an only if there exists a controller that implements it. We also investigate the regularity of the canonical controllers, an establish the fact that they are maximally irregular. This means a canonical controller is regular if an only if every other controller that implements the esire behavior is regular Elsevier B.V. All rights reserve. Keywors: Behaviors; Behavioral control; Regular interconnection; Regular controller; Canonical controller; Implementability 1. Introuction Control problems, seen from the behavioral systems theory point of view, amount to fining a controller, which when interconnecte with the plant in a specifie way yiels the esire behavior [14,1]. The problem may be formulate as follows. Consier a plant to be controlle which has two kins of : Corresponing author. Tel.: ; fax: aresses: a.a.julius@math.utwente.nl (A.A. Julius, jan.willems@esat.kuleuven.ac.be (J.C. Willems, belur@ee.iitb.ac.in (M.N. Belur, h.l.trentelman@math.rug.nl (H.L. Trentelman. an control. A controller is a evice that is attache to the control an restricts their behavior. This restriction is impose on the plant, such that it affects the behavior of the (see Fig. 1. The resulting behavior is calle the controlle system. In this paper we iscuss the properties of a special type of controllers, the so-calle canonical controllers [11,10,16]. We are particularly intereste in their regularity properties. The concepts of canonical an regular controllers will be formally introuce later in this paper. While the behavioral approach sees control as interconnection, the more common point of view in control theory is to view a controller as a feeback /$ - see front matter 2005 Elsevier B.V. All rights reserve. oi: /j.sysconle

2 2 A.A. Julius et al. / Systems & Control Letters ( W CONTROLLER full control CONTROLLED SYSTEM Fig. 1. Control in the behavioral approach. processor that accepts the plant sensor outputs as its inputs an prouces the actuator inputs as its outputs. In [14], this paraigm is calle intelligent control, as the controller acts as an intelligent agent capable of reasoning how to react to sensory observations. The main avantages of the behavioral over the classical feeback point of view are: (i Its practical generality. In many control systems, the controller (i.e. the evice ae to the plant system to obtain a esire behavior oes not act as a sensor/actuator evice. Dampers, heat fins, acoustic noise insulators are examples of such evices. (ii Its theoretical simplicity. Control in the behavioral setting has been introuce in [14], an subsequent evelopment inclues the work in [6,17,9]. The reaer is referre to [12] for further motivation an etails. Throughout this paper, we shall enote the control as c an the as w. These take their value at any given time from their respective signal spaces C an W. Let W an C enote the set of all trajectories of the w an c that are a priori possible, before we have even moelle the plant. In ynamical systems, W an C are typically the set of (smooth signals from the time axis to the signal spaces W an C. In iscrete event systems, the time axis is iscrete an the signal spaces are typically calle alphabets. We shall now iscuss our problem from a purely set theoretic point of view. The behavioral moel of the plant system that captures the relevant relation between w an c is calle the full plant behavior, which is enote by P full. Naturally, we assume that P full is C Fig. 2. The relation between P full, P, C, an K. containe in W C. The full plant behavior consists of all signal pairs (w, c compatible with the plant ynamics. If we project the full behavior on W, weobtain the so calle manifest behavior, which is enote as P. Thus, P := {w W c C such that (w, c P full }. A controller C is a subset of C, containing all signals c allowe by the controller. The controlle behavior is then efine as K := {w W c C such that (w, c P full an c C}. The relationship between the full plant behavior, the manifest behavior, the controller an the controlle behavior is capture in Fig. 2. In this framework, the control problem can be formulate as to fin a controller C that yiels a esire controlle behavior D. Hence the controller C shoul yiel the controlle behavior K = D. We call D the esire controlle behavior. If it is possible to fin a controller C that yiels K = D, then D is sai to be implementable 1 or implemente by C. Further, if a given esire D is implementable, we say that the control problem is solvable. 1 In the literature the term achievable is sometimes use instea of implementable. C

3 A.A. Julius et al. / Systems & Control Letters ( 3 The remainer of this paper is organize as follows. In Section 2 we present two constructions for the canonical controllers an their properties in the setting of general behaviors. In Section 3 we iscuss the concept of implementability, particularly for LTI behaviors. Section 4 is evote for iscussion on the concept of regularity. In Section 5, we stuy the behavior implemente by the canonical controller an use the result from Section 4 to establish its regularity. control CANONICAL CONTROLLER DESIRED CONTROLLED BEHAVIOR 2. The construction an properties of the canonical controllers In the previous section we explaine that we work in the generality in which the plant has two types of, the an the control, an further that a controller can put a restriction on just the control. This restriction is propagate through the plant to the. The iea of the canonical controller uses the internal moel principle in the following way [11,10]. We use a plant moel that has the same behavior as the plant. The propagation of information explaine in the previous paragraph is then reverse by interconnecting the plant moel to the esire behavior D using the 2. This is how we construct the canonical controller. Fig. 3 illustrates this construction. Notice that the wor is mirrore to highlight the fact that the interconnection is reverse. The behavior of the canonical controller obtaine using this construction is enote as C canonical. C canonical := {c C v W such that (v, c P full an v D}. (1 Fig. 4 provies a block iagram showing how C canonical is applie. In [11,10], it is proven that for a class of plants, which are homogeneous in the plant an control, the control problem is solvable (i.e. D is implementable if an only if the canonical controller C canonical implements D. We now efine the homogeneity property. 2 A construction similar to the canonical controller has been use in [7] Fig. 3. The construction of C canonical. Mirrore text reflects the iea that the interconnection is reverse. Definition 1. A full plant behavior P full is sai to have the homogeneity property if for any w 1,w 2 W an c 1,c 2 C, the following implication hols: (w 1,c 1, (w 1,c 2, (w 2,c 1 P full (w 2,c 2 P full. Homogeneity can also be unerstoo as follows. The behavior P full can be seen as a relation between W an C. A relation is calle inepenent if it can be written as a Cartesian prouct of its projections on the relate omains. The behavior P full has the homogeneity property if it can be written as a isjoint union of inepenent relations. In particular, if W an C are linear spaces an P full is a linear subspace, then P full has the homogeneity property. The following theorem captures an important property of the canonical controller C canonical. Theorem 2. If the full plant behavior P full has the homogeneity property, then the canonical controller C canonical implements the smallest implementable behavior containing P D. Proof. Denote the behavior implemente by C canonical as K. From (1, we can infer that K P D. To show that C canonical implements the smallest implementable behavior containing P D, consier any other controller C that implements K such that K P D. We shall prove that K K. Take any element w K. We are going to show that w K. If w P D, then w K, since K P D.

4 4 A.A. Julius et al. / Systems & Control Letters ( DESIRED CONTROLLED BEHAVIOR control CANONICAL CONTROLLER Fig. 4. The canonical controller C canonical in action. If w/ P D, there exists a c C an w P D such that both (w, c an (w,c are elements of P full. Now we are going to show that w/ K is a contraiction. Suppose that it is true, then {c C (w, c P full } C =. By the homogeneity property, we also have that {c C (w, c P full }={c C (w,c P full }. Thus, the following relation is also true: {c C (w,c P full } C =. This implies w / K, which is a contraiction. Remark3. Theorem 2 also tells us that for every control problem involving a plant with the homogeneity property, the smallest implementable behavior containing P D exists. Obviously, a necessary conition for implementability of D is D P. This fact, combine with Theorem 2 gives C canonical its special property. Corollary 4. Assume the full plant behavior has the homogeneity property. Then, the canonical controller C canonical implements D if an only if D is implementable (that is, if an only the control problem is solvable. As alreay note, if P full is a linear behavior, then it has the homogeneity property. Hence, although seemingly restrictive, the class of behaviors with the homogeneity property is, in fact, fairly large, an most importantly, it captures the class of linear time-invariant behaviors, the subject of Sections 3 5. We give LTI behaviors as examples of behaviors with the homogeneity property. For that of behaviors without the homogeneity property, refer to the plant behavior epicte in Fig. 2. There is, in fact, a secon canonical controller that is of interest, an has been introuce in [16]. The canonical controller, enote as C canonical, is efine as follows, C canonical := {c C v such that (v, c P full, an (v, c P full v D }. (2 In other wors, this canonical controller accepts a control-variable trajectory c if an only if every tobe-controlle- trajectory v that can be paire with c is accepte in the esire behavior. Clearly, whatever behavior is implemente by this controller, it must be containe in D. In fact, we have the following theorem. Theorem 5. The canonical controller C canonical implements the largest implementable behavior containe in D.

5 A.A. Julius et al. / Systems & Control Letters ( 5 DESIRED CONTROLLED BEHAVIOR control CANONICAL CONTROLLER Fig. 5. The canonical controller C canonical in action. Proof. Denote the behavior implemente by C canonical as K. It is quite obvious that K D. To show that C canonical implements the largest implementable behavior in D, consier any other controller C that implements K such that K D. We shall prove that K K. Take any element w K. There exists a c C such that K W p full (w, c P full, {v W (v, c P full } K D. (3 K From (2 an (3, we can infer that c C therefore w K. canonical an Remark6. Using Theorem 5, we can also infer that for every control problem (not necessarily the ones with homogeneity property, the largest implementable behavior containe in D exists. As a consequence of Theorem 5, the canonical controller C canonical possesses the following special property. Corollary 7. The canonical controller C canonical implements D if an only if D is implementable (that is, if the control problem is solvable. Fig. 5 illustrates the action of C canonical. Notice that the connectors are replace with symbols enoting implies. For comparison, the behaviors C canonical C canonical Fig. 6. Comparison between the C canonical an C canonical. implemente by the two canonical controllers are shown in Fig Implementability of linear behaviors In the remaining of the paper we shall restrict our attention to LTI ifferential behaviors. This class of behaviors has been iscusse quite extensively in the literature, see [13,6,12], for example. In the following we give a brief introuction to the subject to make this paper self-containe. C

6 6 A.A. Julius et al. / Systems & Control Letters ( We use the symbol L w to enote the class of LTI ifferential 3 systems with w. These are ynamical systems Σ = (R, R w, B, where the behavior B can be expresse as the solutions of a system of ifferential equations ( R w = 0. (4 Here the polynomial matrix R R w [ξ]. Furthermore, the behavior B is efine to be the set of all smooth solutions to the ifferential equations, { B := w C (R, R w ( } R w = 0. (5 The ifferential equation (4 is calle a kernel representation of B an sometimes we write B = ker R(/. A kernel representation is calle minimal if the rank of the polynomial matrix is equal to the number of its rows. Often, the behavior is efine through auxiliary. In this case, we use the term manifest for the of interest, an latent for the auxiliary ones. If B L w+l is a system involving the manifest w an the latent l then it can be proven (see [13,6] that the manifest behavior B w efine by B w := {w C (R, R w l C (R, R l such that (w, l B} is also an element of L w. This result is referre to as the elimination theorem. Remark8. The choice of the unerlying function space C (R, R w is mae for mathematical convenience. An alternative that is quite commonly use is to regar the behavior as the collection of weak solutions of the ifferential equation (4, which are elements of L loc 1 (R, Rw, the space of locally integrable functions. We refer to [6] for further exposition on this issue. We return to the control problem iscusse in Section 1. For linear time-invariant ifferential systems, the control problem can be formulate as follows. The plant behavior P full L w+c is expresse in terms of 3 The analysis also hols if ifference equations were use instea of ifferential equations; however, we restrict our attention to ifferential equations in orer to ease the exposition. the w an the control c. The controller behavior C is an element of L c. The controlle behavior K efine by K := {w C (R, R w c C such that (w, c P full }, is an element of L w (as a consequence of the elimination theorem. For linear ifferential systems, the implementability question becomes: Question. Given P full L w+c, which behaviors K L w can be implemente by using a suitable controller C L c? The answer to this question is summarize in the following theorem. Theorem 9. Given P full L w+c, the behavior K L w is implementable if an only if N K P, (6 where N L w is the hien behavior efine by N := {w C (R, R w (w, 0 P full }, an P is the manifest plant behavior efine by P := {w C (R, R w c C (R, R c such that (w, c P full }. This result was first publishe in [15] an subsequently use in [2,5,8 11]. It is important to notice that the controller that implements K is usually not unique. Generally speaking, the controllers that implement the same behavior may have very ifferent properties. 4. Regular interconnections In the previous section, we have been iscussing interconnection of linear ifferential behaviors without consiering any further restrictions. We now introuce a notion of compatibility in the control problem. In orer to motivate it, consier the following example.

7 A.A. Julius et al. / Systems & Control Letters ( 7 Let P an C be efine as the following: { P := w C (R, R w } 2 w 2 w = 0, { C := w C (R, R w w } + w = 0. We can easily verify that P is an autonomous behavior [6], that is, any trajectory in P is completely characterize by its past. Moreover, this behavior has unstable exponential trajectories. Thus, P is an unstable autonomous behavior. Now, consier the interconnection P C, which consists of the trajectories in the intersection of P an C, i.e. P C. (Here, enotes the interconnection of two systems, while enotes the intersection of the behaviors of the two systems. For the above example, notice that the unstable trajectories in P (the trajectories that are not boune as t o not belong to the interconnection P C. Moreover, since all trajectories in P C are stable (i.e., lim t w(t = 0 for all elements w in P C weinfer that the interconnection of P an C yiels a stable behavior. Therefore, if we o not a any further restrictions on the amissible controllers, it is perfectly possible that an autonomous unstable behavior is stabilize. Such controllers may be impossible to implement. More on this can be foun in [14,3,10]. In orer to cope with this, we introuce the concept of compatibility. With this concept, interconnections like the one in the previous paragraph are iscounte. The control problem then becomes to fin a compatible controller instea of just to fin any controller. A notion of compatibility for behavior interconnections in a general sense, not limite to just LTI systems, has been stuie in [3]. For LTI systems, this general notion is relate to the concept of regular interconnection, that has been introuce before in [14]. Consier a behavior B L w. Let R(/w=0 be a minimal kernel representation of B. Being minimal, R has at least as many columns as it has rows. Let g be the number of rows. The number of columns is obviously w. Therefore, g w. This means that we are always able to select g columns from R to form a square polynomial matrix with nonzero eterminant. Notice that the selection is generally not unique. If we group together the components of w corresponing to the g selecte rows an call them y, an o similarly to the remaining w g components an call them u, we en up with partitioning w into output y an input u. The reason u is calle input is because it is free, in the following sense. For any choice of C input trajectory u, we can always fin an output trajectory y such that (u, y B. Notice that the number of inputs an outputs are properties of the behavior, an are inepenent of the minimal kernel representation use to represent the behavior. Therefore, we can efine two maps m an p such that m : L w {0, 1,...,w} an p : L w {0, 1,...,w}, which give the number of inputs an the number of outputs of a given behavior, respectively. Obviously, for any behavior B L w,we have that m(b + p(b = w. Definition 10. The interconnection of two behaviors B 1 an B 2 is sai to be regular if p(b 1 + p(b 2 = p(b 1 B 2. In a sense, regularity implies that the set of equations governing the ynamics of both behaviors are inepenent of each other. We shall now apply regularity to the control problem. Recalling the efinitions of the full plant behavior P full an the controller behavior C, we efine the full controlle behavior K full as K full := {(w, c P full c C}. The interconnection between the plant an the controller is regular if p(p full + p(c = p(k full. If this is the case, then the controller C is calle a regular controller. It can be shown that a controller is regular if an only if it can be realize as a (possibly non-proper transfer function from the output to the input of P full. Therefore, a controller is regular if it can be viewe as an intelligent controller that processes sensor outputs into actuator inputs. See [14] for more etails. Given the formulation of regular interconnection, we recast the question of implementability of K in the previous section into the question of regular implementability of K.

8 8 A.A. Julius et al. / Systems & Control Letters ( Question. Given P full L w+c, which behaviors K L w can be implemente by using a suitable regular controller C L c? It turns out that regular implementability involves controllability of the plant. For more on the concept of controllability from the behavioral systems theory point of view, we refer the reaer to [6]. In fact, it has been proven in [4,14] that every implementable behavior K of a controllable plant P is regularly implementable. A necessary an sufficient conition for regular implementability, even when the plant may not be controllable, is given in [2]. Theorem 11. Given a full plant behavior P full L w+c. Denote the hien behavior an the manifest plant behavior as N an P respectively. Let P controllable be the controllable part of P. The behavior K L w is regularly implementable if an only if (1 K is implementable, i.e. N K P an (2 K + P controllable = P. Regular implementability of a behavior K implies the existence of at least one regular controller that implements it. In general, given a regularly implementable behavior, there exist irregular controllers that implement it. The question that we aress in the rest of this section is: Uner what conitions on the plant P full an the controlle behavior K, can we conclue that every controller that implements K is regular? It turns out that the answer to this question oes not epen on K, but just on the plant. Define the control variable plant behavior P c L c as follows: P c := {c w such that (w, c P full }. We have the following result. Theorem 12. Let P full L w+c be the full plant behavior. Given any controlle behavior K L w, every controller C L c that implements K is a regular controller if an only if P c = C (R, R c. Proof. Let ( R w + M ( c = 0 be a minimal kernel representation of P full. Note that P c = C (R, R c is equivalent to R having full row rank. (If Take any controller C L c. Let C(/c= 0 be its minimal kernel representation. Since R has full row rank, it follows that ( ( R M ( 0 C ( ( w = 0 c is a minimal kernel representation of K full. Therefore, p(k full = rank R + rank C = p(p full + p(c. Hence the controller is regular. (Only if Suppose that P c = C (R, R c. Let P(/c = 0 be a minimal kernel representation of P c. Note that P = 0. It follows that if we choose a controller that has the same minimal kernel representation as P c, then the resulting interconnection is not regular. 5. Control with the canonical controller Let us revisit the formulation of the control problem for linear time-invariant systems. We are given a full plant behavior P full. Let ( ( R w + M c = 0 (7 be a minimal kernel representation of P full. We are also given a esire controlle behavior D, whose minimal kernel representation is D(/w = 0. The behavior of the first canonical controller C canonical Lc is efine as C canonical := {c C (R, R c v W such that (v, c P full an v D}. Obviously, a kernel representation for this controller can be obtaine by eliminating v from the following kernel representation: ( ( R M ( ( D ( v = 0. (8 0 c

9 A.A. Julius et al. / Systems & Control Letters ( 9 The behavior of the secon canonical controller is efine as C canonical := {c C (R, R c v { } (v, c Pfull, an such that. (v, c P full v D (9 Recall the fact that LTI behaviors satisfy the homogeneity property. For linear time-invariant ifferential systems, these two canonical controllers are essentially equivalent, as shown in the following theorem. Theorem 13. The following three statements are equivalent: (i C canonical = C canonical. (ii The secon canonical controller C canonical is not empty. (iii The hien behavior N is containe in the esire controlle behavior D, i.e. N D. Proof. We prove (i (ii (iii (i. (i (ii: Notice that the zero trajectory is always containe in C canonical, therefore C canonical is never empty. (ii (iii: We first prove that for any c 1 an c 2 in C canonical, their linear combinations are also in. Notice that from efinition (9, it is not C canonical clear that this is the case. Suppose that c 1 an c 2 are in C canonical. There exist w 1 an w 2, both in D, such that (w 1,c 1 an (w 2,c 2 are both in P full. Now take any linear combination α 1 c 1 + α 2 c 2. For any w such that (w, α 1 c 1 +α 2 c 2 P full, the following reasoning hols: (w, α 1 c 1 + α 2 c 2 P full linearity of P full (w α 1 w 1 α 2 w 2 + w 1,c 1 P full, property of C canonical (w α 1 w 1 α 2 w 2 + w 1 D, linearity of D w D. Hence, if C canonical is nonempty, it is obvious that the zero trajectory is inclue in C canonical. Let K canonical be the controlle behavior implemente by C canonical. Since 0 C canonical, we have that N K canonical. From Theorem 5 we also know that K canonical D. Hence N D. (iii (i: Takeanyc C canonical. There exists a w D such that (w, c P full. Take any other w such that (w,c P full, then we also have (w w, 0 P full. Therefore (w w N D. By the linearity of D, we conclue that w D an therefore c C canonical. We have shown that C canonical C canonical. The converse is obvious from the efinitions of the canonical controllers. This theorem implies that if the control problem is solvable, then the two canonical controllers are equal. Motivate by this theorem, we shall consier in the subsequent iscussion only the first canonical controller C canonical. The question what controlle behavior is actually implemente by the canonical controller is answere by the following theorem. Theorem 14. Consier P full L w+c an D L w. The controlle behavior K implemente by the canonical controller C canonical Lc is K = N + D P with N the hien behavior an P the manifest plant behavior. Proof. Let the kernel representation of P full be given by (7, an let D be represente by D(/w = 0. We then know that a kernel representation of C canonical can be obtaine by eliminating w from (8. Therefore, K is the manifest behavior (with w as the manifest variable of the behavior represente by R ( M ( 0 0 M ( R ( [ ] w 0 0 D ( c = 0. v Notice that K is then also the manifest behavior (with w as the manifest variable of the behavior represente by R ( M ( R ( 0 0 D ( [ ] w v c = 0. (10 v Now efine w := w v, we can see from (10 that the ynamics of w is ecouple from that of c an v. Furthermore, the behavior of w is exactly N (see Section 3. The secon an thir rows of (10 inicate that the behavior of v, which is obtaine by eliminating c, isd P. From here, using the fact

10 10 A.A. Julius et al. / Systems & Control Letters ( that w = w + v, we obtain K = N + D P. This result is not unexpecte. In fact, we can see it as an application of Theorem 2 to the special case of LTI behaviors. Similarly, we can apply Corollary 4 to LTI systems to obtain the following corollary. Corollary 15. The canonical controller C canonical L c implements D L w if an only if D is implementable, i.e. N D P. So far we have seen that when D is implementable, the canonical controller C canonical implements it. However, as we have seen in Section 4, implementability alone may not be goo enough. In the following we shall aress the issue of regularity of the canonical controller C canonical. Theorem 16. Given a full plant behavior P full L w+c an a esire controlle behavior D L w. Assume that D is implementable. The canonical controller C canonical implements D regularly if an only if P c = C (R, R c. Proof. (If Follows irectly from Theorem 12. (Only if Without loss of generality, we can assume that P full has a minimal kernel representation of the following form. [ R1 ( M 1 ( 0 M 2 ( ][ ] w = 0, c with both R 1 an M 2 having full row rank. The kernel representations of N an P c are then given by R 1 (/w=0 an M 2 (/c=0, respectively. Since N D, we are able to fin a suitable full row rank matrix F(/ such that F(/R 1 (/w= 0isa minimal kernel representation of D. Therefore a kernel representation of the canonical controller C canonical can be obtaine by eliminating v from ( R ( 1 M 1 ( 0 M 2 F ( ( R1 0 [ v c ] = 0. Since R 1 has full row rank, we easily obtain the following kernel representation of C canonical (possibly non-minimal. [ ( M ] 2 F ( ( c = 0. M1 We see that C canonical always repeats some laws of P full, namely the rows in M 2. Thus C canonical is regular only if M 2 is the zero matrix, which implies P c = C (R, R c. This result, combine with Theorem 12, tells us that the canonical controller is maximally irregular, in the sense that if there exists any irregular controller that implements the esire behavior D, then the canonical controller is irregular too. 6. Concluing remarks The iea of the canonical controller in the behavioral framework is attractive because of its simplicity of construction an also since it formalizes the internal moel principle without unue recourse to the equations with which the plant is escribe. This approach of builing systems without using the equations explicitly unerlines the representation free nature of behavioral theory. Some specific issues are summarize here to highlight the main results of the paper. We efine two canonical controllers. For the case of linear time-invariant behaviors (or more generally, for plants with homogeneity property, when the esire behavior D is implementable, the canonical controllers are the same (see Theorem 13. However, they can iffer when D is not implementable. In this situation each of the two canonical controllers are extreme in a certain sense. The first canonical controller implements the smallest implementable behavior that contains P D. The secon canonical controller C canonical implements the largest implementable behavior containe in D. These statements were formulate an prove in Theorems 2 an 5, respectively. For the case of linear time-invariant behaviors, the implementability of D implies non-emptiness of C canonical (Theorem 13. When C canonical is nonempty, it is a linear subspace of C (R, R, an hence contains the zero trajectory. Thus, we have a necessary conition for implementability of D, namely if the zero trajectory belongs to C canonical.

11 A.A. Julius et al. / Systems & Control Letters ( 11 We then aresse the issue of regularity of a controller with respect to a given plant. It turns out that the conition of guarantee regularity of every controller that implements a given esire behavior, is a property of just the plant, an is inepenent of the given esire behavior. The issues of regularity of every controller an the canonical controllers are relate by the results in the final section. Here we showe that, given a plant, irregularity of any controller implies irregularity of the canonical controller, an hence we terme the canonical controller as maximally irregular. Acknowlegements This research is supporte by the Dutch NWO Project No , Compositional Analysis an Specification of Hybri Systems, by the Belgian Feeral Government uner the DWTC program Interuniversity Attraction Poles, Phase V, , Dynamical Systems an Control: Computation, Ientification an Moelling, by the KUL Concerte Research Action (GOA MEFISTO 666, an by several grants an projects from IWT-Flaners an the Flemish Fun for Scientific Research. References [1] M.N. Belur, Control in a behavioral context, Ph.D. Thesis, University of Groningen, June [2] M.N. Belur, H.L. Trentelman, Stabilization, pole placement an regular implementability, IEEE Trans. Automat. Control 47 ( [3] A.A. Julius, A.J. van er Schaft, Compatibility of behavior interconnections, in: Proceeings of the European Control Conference, Cambrige, IEE, Lonon, September [4] J.W. Polerman, Sequential continuous time aaptive control: a behavioral approach, in: Proceeings of the 39th IEEE Conference on Decision an Control, Syney, IEEE, New York, 2000, pp [5] J.W. Polerman, I. Mareels, A behavioral approach to aaptive control, in: J.W. Polerman, H.L. Trentelman (Es., Mathematics of System an Control, Festschrift on the occasion of the 60th birthay of J.C. Willems, Groningen, 1999, pp [6] J.W. Polerman, J.C. Willems, Introuction to Mathematical Systems Theory: A Behavioral Approach, Springer, New York, [7] P. Rocha, Regular implementability of nd behaviors, in: Proceeings of Mathematical Theory of Networks an Systems, University of Notre Dame, [8] H.L. Trentelman, A truly behavioral approch to the control problem, in: J.W. Polerman, H.L. Trentelman (Es., Mathematics of System an Control, Festschrift on the occasion of the 60th birthay of J.C. Willems, Groningen, 1999, pp [9] H.L. Trentelman, J.C. Willems, Synthesis of issipative systems using quaratic ifferential forms: part II, IEEE Trans. Automat. Control 47 ( [10] A.J. van er Schaft, Achievable behavior of general systems, Systems Control Lett. 49 ( [11] A.J. van er Schaft, A.A. Julius, Achievable behavior by composition, in: Proceeings of the 41st IEEE Conference on Decision an Control, Las Vegas, IEEE, New York, 2002, pp [12] J.C. Willems, jwillems. [13] J.C. Willems, Paraigms an puzzles in the theory of ynamical systems, IEEE Trans. Automat. Control 36 ( [14] J.C. Willems, On interconnections, control an feeback, IEEE Trans. Automat. Control 42 ( [15] J.C. Willems, Behaviors, latent, an interconnections, Systems Control Inform. (Japan 43 ( [16] J.C. Willems, M.N. Belur, A.A. Julius, H.L. Trentelman, The canonical controller an its regularity, in: Proceeings of the IEEE Conference on Decision an Control, Hawaii, December 2003, pp [17] J.C. Willems, H.L. Trentelman, Synthesis of issipative systems using quaratic ifferential forms, Part I, IEEE Trans. Automat. Control 47 (

Systems & Control Letters

Systems & Control Letters Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen

More 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

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

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

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

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

θ 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

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

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

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

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

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

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

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

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

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

On the enumeration of partitions with summands in arithmetic progression

On the enumeration of partitions with summands in arithmetic progression AUSTRALASIAN JOURNAL OF COMBINATORICS Volume 8 (003), Pages 149 159 On the enumeration of partitions with summans in arithmetic progression M. A. Nyblom C. Evans Department of Mathematics an Statistics

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

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

'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

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 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

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

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

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email

More information

Linear and quadratic approximation

Linear and quadratic approximation Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function

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

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower Bounds for the Smoothed Number of Pareto optimal Solutions Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.

More information

On at systems behaviors and observable image representations

On at systems behaviors and observable image representations Available online at www.sciencedirect.com Systems & Control Letters 51 (2004) 51 55 www.elsevier.com/locate/sysconle On at systems behaviors and observable image representations H.L. Trentelman Institute

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

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

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

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

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

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

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

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

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

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

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

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

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling Balancing Expecte an Worst-Case Utility in Contracting Moels with Asymmetric Information an Pooling R.B.O. erkkamp & W. van en Heuvel & A.P.M. Wagelmans Econometric Institute Report EI2018-01 9th January

More information

arxiv: v1 [physics.class-ph] 20 Dec 2017

arxiv: v1 [physics.class-ph] 20 Dec 2017 arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.

More information

Calculus and optimization

Calculus and optimization Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function

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

Lower bounds on Locality Sensitive Hashing

Lower bounds on Locality Sensitive Hashing Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,

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

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

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

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

REAL ANALYSIS I HOMEWORK 5

REAL ANALYSIS I HOMEWORK 5 REAL ANALYSIS I HOMEWORK 5 CİHAN BAHRAN The questions are from Stein an Shakarchi s text, Chapter 3. 1. Suppose ϕ is an integrable function on R with R ϕ(x)x = 1. Let K δ(x) = δ ϕ(x/δ), δ > 0. (a) Prove

More information

Convergence of Random Walks

Convergence of Random Walks Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of

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

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

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

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

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything

More 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

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

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 appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution

More information

Conservation laws a simple application to the telegraph equation

Conservation laws a simple application to the telegraph equation J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness

More 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

Witt#5: Around the integrality criterion 9.93 [version 1.1 (21 April 2013), not completed, not proofread]

Witt#5: Around the integrality criterion 9.93 [version 1.1 (21 April 2013), not completed, not proofread] Witt vectors. Part 1 Michiel Hazewinkel Sienotes by Darij Grinberg Witt#5: Aroun the integrality criterion 9.93 [version 1.1 21 April 2013, not complete, not proofrea In [1, section 9.93, Hazewinkel states

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

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

OPTIMAL CONTROL PROBLEM FOR PROCESSES REPRESENTED BY STOCHASTIC SEQUENTIAL MACHINE

OPTIMAL CONTROL PROBLEM FOR PROCESSES REPRESENTED BY STOCHASTIC SEQUENTIAL MACHINE OPTIMA CONTRO PROBEM FOR PROCESSES REPRESENTED BY STOCHASTIC SEQUENTIA MACHINE Yaup H. HACI an Muhammet CANDAN Department of Mathematics, Canaale Onseiz Mart University, Canaale, Turey ABSTRACT In this

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

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

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs

Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Ashish Goel Michael Kapralov Sanjeev Khanna Abstract We consier the well-stuie problem of fining a perfect matching in -regular bipartite

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

Two formulas for the Euler ϕ-function

Two formulas for the Euler ϕ-function Two formulas for the Euler ϕ-function Robert Frieman A multiplication formula for ϕ(n) The first formula we want to prove is the following: Theorem 1. If n 1 an n 2 are relatively prime positive integers,

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

The maximum sustainable yield of Allee dynamic system

The maximum sustainable yield of Allee dynamic system Ecological Moelling 154 (2002) 1 7 www.elsevier.com/locate/ecolmoel The maximum sustainable yiel of Allee ynamic system Zhen-Shan Lin a, *, Bai-Lian Li b a Department of Geography, Nanjing Normal Uni ersity,

More information

ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa

ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa Manuscript submitte to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS Gianluca Crippa Departement Mathematik

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

Mathematical Biosciences

Mathematical Biosciences Mathematical Biosciences 26 (2008) 40 49 Contents lists available at ScienceDirect Mathematical Biosciences journal homepage: www.elsevier.com/late/mbs Homotopy methos for counting reaction network equilibria

More information

Sliding mode approach to congestion control in connection-oriented communication networks

Sliding mode approach to congestion control in connection-oriented communication networks JOURNAL OF APPLIED COMPUTER SCIENCE Vol. xx. No xx (200x), pp. xx-xx Sliing moe approach to congestion control in connection-oriente communication networks Anrzej Bartoszewicz, Justyna Żuk Technical University

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

Generalized Tractability for Multivariate Problems

Generalized Tractability for Multivariate Problems Generalize Tractability for Multivariate Problems Part II: Linear Tensor Prouct Problems, Linear Information, an Unrestricte Tractability Michael Gnewuch Department of Computer Science, University of Kiel,

More information

Efficient Construction of Semilinear Representations of Languages Accepted by Unary NFA

Efficient Construction of Semilinear Representations of Languages Accepted by Unary NFA Efficient Construction of Semilinear Representations of Languages Accepte by Unary NFA Zeněk Sawa Center for Applie Cybernetics, Department of Computer Science Technical University of Ostrava 17. listopau

More information

Iterated Point-Line Configurations Grow Doubly-Exponentially

Iterated Point-Line Configurations Grow Doubly-Exponentially Iterate Point-Line Configurations Grow Doubly-Exponentially Joshua Cooper an Mark Walters July 9, 008 Abstract Begin with a set of four points in the real plane in general position. A to this collection

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

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

Acute sets in Euclidean spaces

Acute sets in Euclidean spaces Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of

More information

4. Important theorems in quantum mechanics

4. Important theorems in quantum mechanics TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this

More 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

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

Energy behaviour of the Boris method for charged-particle dynamics

Energy behaviour of the Boris method for charged-particle dynamics Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion

More information

TOTAL ENERGY SHAPING CONTROL OF MECHANICAL SYSTEMS: SIMPLIFYING THE MATCHING EQUATIONS VIA COORDINATE CHANGES

TOTAL ENERGY SHAPING CONTROL OF MECHANICAL SYSTEMS: SIMPLIFYING THE MATCHING EQUATIONS VIA COORDINATE CHANGES TOTAL ENERGY SHAPING CONTROL OF MECHANICAL SYSTEMS: SIMPLIFYING THE MATCHING EQUATIONS VIA COORDINATE CHANGES Giuseppe Viola,1 Romeo Ortega,2 Ravi Banavar Jose Angel Acosta,3 Alessanro Astolfi, Dipartimento

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

Section 2.7 Derivatives of powers of functions

Section 2.7 Derivatives of powers of functions Section 2.7 Derivatives of powers of functions (3/19/08) Overview: In this section we iscuss the Chain Rule formula for the erivatives of composite functions that are forme by taking powers of other functions.

More information

UNDERSTANDING INTEGRATION

UNDERSTANDING INTEGRATION UNDERSTANDING INTEGRATION Dear Reaer The concept of Integration, mathematically speaking, is the "Inverse" of the concept of result, the integration of, woul give us back the function f(). This, in a way,

More information

On classical orthogonal polynomials and differential operators

On classical orthogonal polynomials and differential operators INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 38 2005) 6379 6383 oi:10.1088/0305-4470/38/28/010 On classical orthogonal polynomials an ifferential

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

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

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information