On Nonlinear Controllability of Homogeneous Systems Linear in Control

Size: px
Start display at page:

Download "On Nonlinear Controllability of Homogeneous Systems Linear in Control"

Transcription

1 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 1, JANUARY On Nonlinear Controllability of Hoogeneous Systes Linear in Control Jaes Melody, Taer Başar, and Francesco Bullo Abstract This work considers sall-tie local controllability (STLC) of single- and ultiple-input systes, _ = ( )+ where ( ) contains hoogeneous polynoials and... are constant vector fields. For single-input systes, it is shown that even-degree hoogeneity precludes STLC if the state diension is larger than one. This, along with the obvious result that for odd-degree hoogeneous systes STLC is equivalent to accessibility, provides a coplete characterization of STLC for this class of systes. In the ultiple-input case, transforations on the input space are applied to hoogeneous systes of degree two, an exaple of this type of syste being otion of a rigid-body in a plane. Such input transforations are related via consideration of a tensor on the tangent space to congruence transforation of a atrix to one with zeros on the diagonal. Conditions are given for successful neutralization of bad type (1,2) brackets via congruence transforations. Index Ters Controllability, Lie algebras, nonlinear systes. I. INTRODUCTION Various concepts of controllability for nonlinear systes were initially explored in [1] [3]. In particular, [2] is priarily concerned with the property of accessibility of the analytic control syste _x = F (x; u), naely that the set of points attainable fro a given initial point via application of feasible input is full in the sense of having a nonepty interior. Sussann and Jurdjević deonstrated in [2] that a necessary and sufficient condition for accessibility of these systes is that the Lie algebra generated by the syste have full rank, the so-called Lie algebra rank condition (LARC). In [4], Sussann explored the property of sall-tie local controllability (STLC) for affine analytic single-input systes _x = f0(x) + f1(x)u with juj 1. A syste is said to be STLC at a point x0 if that initial point is in the interior of the set of points attainable fro it in tie T for all T > 0. In this case, the Lie algebra generated by the syste, denoted by L(ff0;f1g), is the sallest involutive distribution containing ff0;f1g or, equivalently, the distribution spanned by iterated Lie brackets of f0 and f1. Sussann gave various necessary and sufficient conditions for STLC in [4]. For exaple, a necessary condition for STLC is that the tangent vector [f1; [f1;f0]] x be in the subspace spanned by all tangent vectors at x0 generated by brackets with only one occurrence of f1, which is denoted by L 1 (ff0;f1g) x. More iportantly, the conditions conjectured by Heres were proved to be sufficient conditions for STLC. These Heres local controllability conditions (HLCC) consist of (1) x0 is a (regular) equilibriu point, (2) the LARC is satisfied, and (3) S k (ff0;f1g) x S k01 (ff0;f1g) x for all even k>1where S k (ff0;f1g) x denotes the span of all tangent vectors at x0 generated by brackets with k or less occurrences of f1. Stefani [5] provided an extension of Sussann s necessary condition by deonstrating that STLC iplies (ad 2 f f0) x 2 S 201 x for all 2f1; 2;...g. For an excellent suary and tutorial of these as well as other results in the single-input case, the inquisitive reader is directed to Kawski [6]. Manuscript received January 30, 2001; revised October 17, 2002 and Septeber 13, Recoended by Associate Editor H. Wang. This work was supported in part by the U.S. Departent of Energy and in part by the U.S. National Science Foundation under Grant CMS The authors are with the Coordinated Science Laboratory, University of Illinois, Urbana, IL USA (e-ail: jelody@uiuc.edu; tbasar@control.csl.uiuc.edu; bullo@uiuc.edu) Digital Object Identifier /TAC STLC of ultiple-input affine analytic control systes was addressed by Sussann in [7], where a general sufficiency theore was proven for analytic systes of the for _x = f0(x) + f i(x)u i with the constraints ju ij 1 for all i 2 f1;...;g. In order to understand this result, it is necessary to distinguish between foral brackets and evaluated brackets. On the one hand, a foral bracket is a pairwise parenthesized word (i.e., an eleent of the free aga) with well-defined left and right factors, nuber of factors, and degree. On the other hand, an evaluated bracket is the vector field that results fro an iterated Lie bracketing of particular vector fields. When we speak of the vector field generated by a foral bracket, we ean specifically the vector field that results fro evaluating the foral bracket with respect to particular vector fields, an operation that is theoretically captured by the evaluation ap in [7]. This distinction is captured by the following notation, which we eploy both in the stateent of Sussann s general sufficiency theore and throughout the sequel: if B =(i1;...;(i k01;i k )...)represents a foral bracket of indeterinates f0;...;g, then f B = [f i ;...[f i ;f i ]...] denotes the corresponding vector field generated by the evaluation ap fro the foral bracket B with respect to a particular set of vector fields ff0;...;f g.1 Hence, for exaple, we have the foral bracket B =(1; (1; 0)) that generates the vector field resulting fro evaluating the corresponding iterated Lie bracket [f1; [f1; f0]] for a particular pair of vector fields f0 and f1. Several results were presented in [7], but in the context of this note the ost appropriate result is based on the degree of a foral bracket. For a given foral bracket B of indeterinates fig i=0, jbj i is used to denote the nuber of occurrences of i as a factor in B.For 2 [0; 1], is defined by (B) :=jbj0 + jbj i. The general theore states that systes that satisfy the LARC at x0 and have f0(x0) =0are STLC if there exists a 2 [0; 1] such that the tangent vector at x0 generated by each foral bracket B with jbj0 odd and jbj1;...; jbj all even can be expressed as a linear cobination of tangent vectors at x0 generated by soe set of brackets fb k g N k=1 with (B k ) < (B) for all k 2f1;...;Ng. It has becoe conventional to refer to the foral brackets with jbj0 odd and jbj1;...; jbj all even as bad brackets, and if these bad brackets have corresponding vector fields that are not contained in the span of vector fields generated by good brackets of lower degree at x0 then they are referred to as potential obstructions to STLC. The obstructions are only potential because they only obstruct the known sufficient conditions. In [6], Kawski presents several exaples of systes with potential obstructions that are known to be STLC. Both Sussann in [7] and Kawski in [6] apply a generalized definition of hoogeneity to STLC. This concept of hoogeneity begins with definition of a dilation as a paraeterized ap of IR n to IR n of the for (x) =( r x1; r x2;... r x n ) where r i are nonnegative integers. A polynoial p :IR n! IR is then said to be hoogeneous of degree k with respect to the dilation, sybolically p 2 H k if p( (x)) = k p(x). Traditional hoogeneity is recovered via the dilation with r1 = 111 = r n = 1. The definition of hoogeneity is then extended to vector fields in the following anner: a vector field f is said to be hoogeneous of degree j if fp 2 H k0j whenever p 2 H k for all k 0. A related area of research that capitalizes on this generalized concept of hoogeneity is that of nilpotent and high-order approxiation of control systes presented by Heres, for exaple, in [9]. One pertinent outcoe of Heres s research is that a syste is STLC if its Taylor approxiation is STLC. However, the converse question of whether STLC can be deterined fro a finite nuber of differentiations is still open [10]. 1 Readers interested in the rich detail of foral Lie algebras ay refer to [8] /03$ IEEE

2 140 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 1, JANUARY 2003 Fig. 2. Motion of a rigid body in a plane expressed in body-fixed coordinates. Fig. 1. Graphical depiction of brackets that generate nonzero vector fields of polynoial syste (1) with hoogeneity degree. In the reainder of this note, we restrict our attention to hoogeneous nonlinear systes that are linear in control, using the traditional definition of hoogeneity. We begin by addressing single-input systes. Building on Stefani s necessary condition [5] and the concepts of good and bad brackets of Sussann s general sufficiency theore [7], we deonstrate that such single-input hoogeneous systes of even degree are STLC if and only if they have a scalar state. This result, cobined with the obvious fact that for odd-degree hoogeneous single-input systes STLC is equivalent to the LARC, copletely characterizes STLC for these systes. Next, we address ultiple-input systes with the additional restriction that they be hoogeneous of degree two. In this case, we extend the applicability of Sussann s general sufficiency theore by incorporating a linear transforation on the ultidiensional input in order to neutralize potential obstructions that arise fro type (1,2) bad brackets (i.e., brackets with jbj 0 =1and jbj!0 =2). In particular, we present a foral ethod for neutralizing these type (1,2) potential obstructions wherein the proble of finding the desired linear transforation on the input space is reduced to finding a particular atrix congruence transforation. II. PROBLEM EXPOSITION In this note, we address STLC of systes of the for _x = f 0 (x) + f i u i (1) where ju i j1and x 2 IR n. f i for i 2f1;...;g are assued to be constant vector fields, i.e., f i (x) f i, and the coponents of f 0 (x) are hoogeneous polynoials of degree k 1. We use the traditional definition of hoogeneous polynoial p, naely that p(x) = k p(x). The set of such hoogeneous vector fields is denoted by H k. Our definition of degree-k hoogeneous vector fields is equivalent to that used in [6], [7], and [9] in the following anner: take : x 7! x and then H k is in the general fraework the set of vector fields hoogeneous of degree 1 0 k. With this traditional definition, we have the following eleentary facts: i) f (0) = 0 for f 2H k with k 1, and ii) [f; g] 2H j+k01 for all f 2H j and g 2H k, where H01 is interpreted as the singleton containing the zero vector field. Systes of this for are theoretically interesting because their Lie algebra at x 0 = 0has a diagonal structure, as depicted in Fig. 1. In particular, the only brackets B that generate vector fields with nonzero value at x 0 are those with jbj!0 =(k01)jbj 0 +1, where jbj!0 := jbji. This follows directly fro the eleentary facts previously given. Keeping in ind that f 0 2H k and f i 2H 0 for all i 6= 0, fro fact ii), it is clear that brackets above the diagonal have hoogeneity degree greater than zero and, hence, by fact i) have zero value at x 0. Siilarly, fro fact ii), we have that brackets below this line have hoogeneity degree less than zero and, hence, by definition are identically zero. Furtherore, systes of the for (1) coonly arise in echanics. An exaple of such a syste is the otion of a rigid body in a plane expressed in body-fixed coordinates, as depicted in Fig. 2. The equations of otion for this syste are _! =u 2 + hu 1 _v x = 0!v y _v y =!vx + u1: (2) The state consists of rotational velocity! and the two body-fixed translation velocities v x and v y. The input consists of the torque u 2 and the force u 1 applied at a oent ar of h. Hence, f 0 (x) = (0; 0x 1 x 3 ;x 1 x 2 ), f 1 = (h; 0; 1) for soe constant h, f 2 = (1; 0; 0), and the syste is of the for (1) with hoogeneity degree two. This provides a siple exaple of a syste for which Sussann s sufficient condition in [7] is not invariant with respect to input transforations. In particular, if a pure force (h = 0) and a torque are used as inputs, then the syste satisfies Sussann s sufficient condition for the ultiple-input case. 2 However, if an offset force (h 6= 0) is used, then potential obstructions appear as vector fields generated fro the type (1,2) brackets, i.e., brackets with jbj 0 =1and jbj!0 =2. In general, we eploy the phrase type (k; `) brackets to refer to all foral brackets B of indeterinates f0;...;g with jbj 0 = k and jbj!0 = `, and denote the distribution spanned by such brackets as L (k;`) (F). 3 Using this syste as a otivating exaple, we explore the neutralization via congruence transforation of potential obstructions generated by bad brackets of type (1,2) with vector fields generated by other brackets (perhaps also bad) of type (1,2). III. SINGLE-INPUT SYSTEMS In this section, we consider the single-input syste _x = f 0 (x) +f 1 u (3) where x 2 IR n, u 2 [01; 1], f 1 2 IR n, and f 0 2H k. In light of HLCC and Sussann s general sufficiency result [7], it is clear that for a syste as in (3) with odd hoogeneity degree, accessibility is equivalent to STLC, since there are no nonzero brackets with j1j 0 odd and j1j 1 even. In other words, the question of STLC reduces to the LARC. The 2 A generalized force on a rigid body consists of a pure force coponent which induces only a translational otion and a torque coponent which induces only a rotational otion. 3 The notation is used instead of to ephasize that ( ) is not necessarily a Lie subalgebra.

3 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 1, JANUARY following lea asserts that if the vector fields generated by brackets of type (1;k) do not add to the Lie algebra rank, then neither do the vector fields generated by higher-degree brackets. Lea 1: For (3) with hoogeneity degree k>0,if(ad k f f 0) x 2 spanff 1 g, then L(ff 0 ;f 1 g) x = spanff 1 g. Proof: Since the Lie algebra structure is invariant to (analytic) coordinate transforations, without loss of generality we can take f 1 to be the basis vector e 1. Let Hk 11 be the set of vector fields of hoogeneity degree k with the for (Cx k (x); 2 (x);...; n (x)) where C is any scalar constant and the power of x 1 in i is less than k for all i 2f1;...;g. Then, (adf k f 0 ) x 2 spanfe 1 g iplies that f 0 2 Hk 11, since adf k corresponds to the partial derivative k =@x k 1. Furtherore, if g 2 H 11 for soe 1, then ad f g 2 H01 11 and ad f g 2H+k Since the Lie subalgebra of (3) is spanned by vector fields of the for [f i ; [f i ;...[f i ;f i ]...]] (apply, for exaple, [11, Prop. 3.8]), all vector fields in L(ff 0;f 1g) x are a ultiple of e 1. Turning our attention to systes with even hoogeneity degree k, we see that in general L(ff 0 ;f 1 g) x does include potential obstructions generated fro bad brackets, i.e., there are bad brackets along the diagonal of Fig. 1. In particular, type (; (k 0 1) +1)brackets are (odd, even) when is odd. We ake use of the necessary condition of Stefani restated here for convenience. Theore 2 Stefani [5]: If the syste _x = f 0(x) +f 1(x)u 1 is STLC, then (adf 2k f 0 ) x 2S 2k01 (ff 0 ;f 1 g) x. When applied to even systes, Theore 2 states that (adf 2k f 0) x 2 spanff 1g is necessary for STLC. However, if (adf 2k f 0 ) x 2 spanff 1 g, then Lea 1 asserts dil(ff 0 ;f 1 g)=1, and systes with n 2 cannot be STLC (for otherwise LARC is violated). This reasoning and the fact that the syste _x = x 2k + u with x 2 IR is STLC provides the following result. Proposition 3: If the syste in (3) has odd hoogeneity degree, then it is STLC if and only if it satisfies the LARC. On the other hand, if the syste has even hoogeneity degree 2k >0, then it is STLC if and only if the state x is scalar. the unchanged equations of otion, we have the potential obstruction [f 1 ; [f 1 ;f 0 ]](0) = (0; 02h; 0) 62 spanff 1 (0);f 2 (0)g. This potential obstruction prevents the use of Sussann s sufficient condition for any h 6= 0. On the other hand, if h =0(corresponding to a force input u 1 through the center of ass) we can then apply the sufficient condition. This is so because we have the vector field [f 1 ; [f 2 ;f 0 ]] generated by the good bracket (1; (2; 0)) being (0; 01; 0), bringing the Lie algebra to full rank, and STLC is deonstrated. Furtherore, the syste for h 6= 0can be transfored into the pure-force syste (h =0)via the input transforation u 1 u 2 = 1 0 0h 1 u 1 u 2 : Since STLC is clearly invariant to full-rank input transforations, we see that a particular choice of input transforation ay provide a eans of reoving a potential obstruction, thus extending the applicability of Sussann s general theore for this class of systes. Reark 4: It is worth noting that Kawski [6] has considered techniques for neutralizing and balancing bad brackets. However, this technique is not related to ours. Kawski s technique applies in the singleinput setting, and neutralizes brackets possibly with brackets of different degree via paraeterized failies of controls carefully tailored to the syste and the brackets in question. In soe cases, these failies of controls involve switching between control liits, with the paraeter affecting the switching ties. Our technique utilizes the freedo of ultiple independent inputs to enforce a linear relation between the inputs in order to neutralize brackets of the sae type. A. Neutralization via Congruence Transfor Moving to the generic ultiple-input hoogeneous syste of degree two, suppose that there is soe bad bracket of the for (i; (i; 0)) that generates a potential obstruction, i.e., [f i; [f i;f 0]] x =2 spanff k (x 0)gk=1. We would like to find a full-rank, linear transforation T on the input space such that the transfored syste f IV. MULTIPLE-INPUT SYSTEMS We now return to consideration of the syste in (1) where f 0 2 H 2 (x). An extension of the previous results to this ultiple-input case is probleatic. In particular, the necessary condition of [5] runs into the proble of a possibility of balancing between potential obstructions generated fro bad brackets of the sae degree. This consideration, along with the otivating exaple of planar rigid-body otion lead us to investigate neutralization of potential obstructions by vector fields generated fro brackets of the sae type. Of course, for the syste in (1) the general sufficient condition of [7] can be applied to deterine STLC. Since the Lie algebra has a diagonal structure, the choice of 2 [0; 1] in the theore is iaterial. Using Sussann s concepts of good and bad brackets, the sufficient condition allows us to neutralize potential obstructions fro bad brackets with vector fields generated by good brackets of lower degree. Our goal with this section is to address the case where there are potential obstructions that cannot be neutralized in this anner, and to neutralize these potential obstructions with vector fields generated by other brackets of the sae degree via appropriate choice of linear transforation on the input space. In this endeavor, the diagonal structure of the Lie algebra will be particularly useful. Returning to the otivating exaple of planar otion in (2), it is clear that this syste is STLC. (For exaple, use the feedback transforation u 1 =!v x +u 1 and u 2 = u 2 0 hu 1 to obtain the syste described by _! = u 2, _v x = 0!v y, and _v y = u 1). However, in attepting to apply Sussann s general sufficiency result directly on _x = f 0(x)+ j=1 f it ij has the corresponding potential obstruction reoved. Of course, the full-rank input transforation will not affect spanff k gk=1. We ake the restriction ju j j1= where is the spectral radius of T in order to have the resulting u i satisfy the bounds ju ij1. 4 Suppose that there is at least one type (1,2) potential obstruction, i.e., there is soe 2 f1;...; g such that [[f ; [f ;f 0]] x =2 L (0;1) (F) x. Consider the codistribution KerL (0;1) (F) that annihilates the distribution L (0;1) (F). Then there ust exist soe differential one-for 2 KerL (0;1) (F) such that ([f ; [f ;f 0]]) x 6= 0. Let ) be the codistribution containing all covectors 2 KerL (0;1) (F) with the property ([f j ; [f j ;f 0 ]]) x 6=0 for soe j. Let us suppose that x has exactly diension one. Choosing any nonzero 2 ), we define the ap fro T IR n 2 T IR n to IR by : (f; g) 7! ([f; [g; f 0 ]]) x. This ap inherits bilinearity fro the Lie bracket and, hence, is a tensor of covariant order two at x 0. Next, we derive a atrix 9 2 IR 2 fro via (9 ) ij := (f i;f j) (4) 4 While odification of the control bound can result in difficulties in the balancing of brackets in [6], this is not a concern in our case, since the neutralization that we achieve is independent of the relative agnitudes of. u j

4 142 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 1, JANUARY 2003 for i; j 2f1;...;gand where f i, f j are the input vector fields of (1). By eploying the Jacobi identity and the fact that the input vector fields coute, it is clear that 9 is also syetric. Denoting by ^9 the corresponding atrix for the transfored syste, it is easy to see that ^9 = T T 9 T. In this anner, the question of whether the obstructing brackets can be neutralized is reduced to the following linear algebra question: given a syetric atrix 9 6=0, is there a full-rank, square atrix T such that the congruence transforation of 9, ^9 = T T 9 T, has all zeros along the diagonal? Supposing for a oent that such a congruence transfor exists. Since it is full rank, it ust be true that there is soe particular ~{ and ~ with ~{ 6= ~ such that ( ^9 ) ~{~ 6= 0. In siplified ters, such an input transforation not only neutralizes the potential obstructions along x, but also replaces the with vector fields generated by good brackets along x. Furtherore, the input transforation will not create type (1,2) potential obstructions that are annihilated by x. (This would be tantaount to T T 0T 6= 0). Of course, the input transforation will also affect the vector fields generated by higher degree brackets, possibly creating potential obstructions. Recalling that a syetric atrix is called indefinite if it has at least one positive eigenvalue and at least one negative eigenvalue, we have the following answer to the posed question. Lea 5: Given a atrix 9 =9 T 6= 0, there exists a full-rank atrix T such that ^9 = T T 9 T has all zeros on the diagonal if and only if 9 is indefinite. Proof: First, recall that by virtue of the syetry of the atrix 9, there exists a choice of orthonoral eigenvectors V := (v 1 ;...;v ) such that V T 9 V =diag( 1 ;...; ) where i are the real eigenvalues of 9. Expressing the coluns t i of T in ters of the orthonoral eigenvectors, we have ( ^9 ) ii = j=1 j (t T i v j ) 2. If 9 6=0is seidefinite, then without loss of generality, we can take 9 0, and hence ( ^9 ) ii = j=1 j (tt i v j ) 2 0. For necessity, we ust show that ( ^9 ) ii > 0. For any full-rank T there is soe colun t i and soe eigenvalue j such that j (ti T v j ) 2 > 0, and since j 0 for all j, wehave( ^9 ) ii > 0. Suppose 9 is indefinite, and group the eigenvalues into those which are positive f + i g, those which are negative f0 j g, and j=1 those which are zero f 0 kg k=1. The eigenvectors are siilarly grouped into fv i + g, fv0 j g, and j=1 fv0 kg k=1. We proceed by constructing the atrix T. The first 0 coluns t i of T are chosen so that t i = v 0 i, achieving ti T 9 t i = 0for i 2f1;...; 0 g. The next + coluns t j are chosen according to t j+ = v j + =(+ j )1=2 0 v 0 =(0 1 1 )1=2 for all j 2f1;...; +g. For this choice, t i?t j for all i 2f1;...; 0g and all j 2 f1;...; + g, and tj+ T 9 t j+ = 0 for all j 2 f1;...; +g. The final 0 coluns t k are chosen to be t k+ + = v + 1 =(+ 1 )1=2 + v 0 k =(0 k )1=2 for all k 2f1;...;0g. Siilarly, this final group of coluns is orthogonal to ft i g and has the property tk+ T + 9 t k+ + =0for all k 2f1;...;0g. Furtherore, ft`g `= +1 is linearly independent. This copletes the construction of T. B. Planar Vehicle Exaple Revisited Applying this line of reasoning to the planar vehicle exaple ) previously presented, the codistribution is spanned by = (0; 02h; 0). In this exaple, if we use the canonical isoorphiss T IR n T 3 IR n IR n and furtherore if we endow IR n with the natural inner product based on a particular choice of basis, we can interpret ) as the projection of the vector fields generated by the type (1,2) brackets onto the orthogonal copleent of L (0;1) (F). It represents the direction in which the vector fields generated by the type (1,2) brackets cannot be neutralized by vector fields generated by the type (0,1) brackets, i.e., the direction of potential obstruction. The tensor in coordinates is (0; 0; 2h;0; 0; 0; 2h; 0; 0). The associated atrix 9 = (4h 2 ; 2h;2h; 0) has eigenvalues 2h 2 6 2p h 4 + h 2. Of course h 6= 0 is assued, for otherwise there is no potential obstruction. It is easy to see that 9 is sign indefinite and, hence, the construction in the proof of Lea 5 provides the transforation T = 1(h)h 0 p p 2(h) 1+h 2 2(h)h 0 1(h) 1+h 2 1 (h) 2 (h) where 1 and 2 are continuous functions of h>0with i > 0 for all finite h > 0. This transforation yields ^9 = (0; 02; 02; 0), and the resulting type (1,2) brackets of the transfored syste generate [ f 1 ; [ f 1 ;f 0 ]] = (0; 0; 0), [ f 1 ; [ f 2 ;f 0 ]] = (0; 1=h; 0), and [ f 2 ; [ f 2 ;f 0 ]] = (0; 0; 0). Thus the potential obstruction to STLC is reoved, and the Lie subalgebra generated by the syste at x 0 is spanned by vector fields corresponding to good brackets, naely ff 1 ; f 2 ; [ f 1 ; [ f 2 ;f 0 ]]g. Hence, application of Sussann s general result [7] deonstrates STLC. It is interesting to note that while T is not equal to T as previously deterined, it does transfor the syste into one with a pure force and a torque input for any h>0. C. Exaple of Balancing Two Bad Brackets Not only can this technique neutralize a potential obstruction fro a bad type (1,2) bracket with a vector field generated by a good type (1,2) bracket, but it ay also balance two potential obstructions generated by two type (1,2) brackets. Consider the two input exaple with f 0(x) = (x 2x 3;x 1x 3;x 2 1 0x 2 2), f 1 =(1; 0; 0), and f 2 =(0; 1; 0). This syste has two potential obstructions of type (1,2), naely [f 1 ; [f 1 ;f 0 ]] = (0; 0; 2) and [f 2; [f 2;f 0]] = (0; 0; 02). Furtherore, the good bracket (1; (2; 0)) evaluates as [f 1 ; [f 2 ;f 0 ]] = (0; 0; 0). For this exaple, 9 = (4; 0; 0; 04) is clearly indefinite and, hence, we have the desired transforation T = (0:5; 00:5; 0:5; 0:5). The resulting type (1,2) brackets for the transfored syste evaluate as [ f 1 ; [ f 1 ;f 0 ]] = (0; 0; 0), [ f 1 ; [ f 2 ;f 0 ]] = (0; 0; 1), and [ f 2 ; [ f 2 ;f 0 ]] = (0; 0; 0), and again STLC is achieved. An interesting variation on this exaple is obtained if we replace f 0 (x) with (x 2 x 3 ;x 1 x 3 ;x 2 1 +x 2 2 +x 1 x 2 ), where 2 IR. This syste has 9 =(4; 2;2; 4) indefinite for jj > 2. This condition has the interpretation that even when the vector fields generated by the two bad brackets have values along with the sae sign, they can still be neutralized with a vector field generated by a good bracket provided that its value along is large enough. D. Effect of Neutralization on Other Directions Next, we consider the effect of neutralization of potential obstructions along one direction within ) on the value of the brackets along another direction within ). We consider the syste that evolves on x 2 IR 4 described by f 0 (x) =(x 2 x 4 ; 0;x 2 1 +x 1 x 2 ;x 1 x 2 ), f 1 = (1; 0; 0; 0), and f 2(0; 1; 0; 0). The type (1,2) brackets for this syste are [f 1 ; [f 1 ;f 0 ]] = (0; 0; 2; 0), [f 2 ; [f 2 ;f 0 ]] = (0; 0; 0; 0), and [f 1 ; [f 2 ;f 0 ]] = (0; 0; 1; 1) and, hence, ) is spanned by the covectors (0; 0; 1; 0) and (0; 0; 0; 1). If we concentrate on neutralizing the bad bracket in the direction =(0; p 0; 2; 0), then we have 9 = (4; 2; 2; 0) with eigenvalues = The constructed transforation atrix T =(02 01=4 ; 2 01=4 ;0; 2 1=4 ) neutralizes the bad bracket along. However, the transforation produces a potential obstruction along the covector (0; 0; 0; 1), as evidenced by the resulting vector fields [ f 1 ; [ f 1 ;f 0 ]] = (0; 0; 0; 0 p 2), [ f 1 ; [ f 2 ;f 0 ]] = (0; 0; 01; 01), and [ f 2 ; [ f 2 ;f 0 ]] = (0; 0; 0; 0). In fact, this syste is known to be not STLC, as can be seen by applying the coordinate transforation y 3 := x 3 0 x 4 and y i = x i for i 6= 3.

5 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 1, JANUARY TABLE I APPLICABILITY OF NEUTRALIZATION VIA CONGRUENCE TRANSFORMATION E. Interpretation and Ipact We have developed a ethodology for neutralizing potential obstructions generated by bad brackets of type (1,2) for hoogeneous degree-two systes that requires indefiniteness of the atrix 9 defined in (4). To interpret this requireent, we recall that a atrix is positive definite if and only if all of its principal inors are positive definite. On the other hand, a atrix is negative definite if and only if all of its principal inors are negative definite when of odd diension and positive definite when of even diension. Hence, 9 will be indefinite if soe principal inor is itself indefinite. Recalling that the ijth entry of 9 is the value of the evaluated bracket [f i ; [f j ;f0]] along the covector, the iplications of the indefiniteness test becoe intuitively clear. For the oent, let us restrict our attention to cases where x has diension one, say x =spanf x g for soe. If there is a potential obstruction fro a single, type (1,2) bad bracket and a type (1,2) good bracket is not annihilated at x0 by, then the obstruction can always be reoved since the principal inor corresponding to these two brackets is always indefinite (i.e., the atrix (2a; b; b; 0) has eigenvalues a 6 p a 2 +4b 2 ). If two or ore evaluated bad brackets are not annihilated at x0 by, then they can all be siultaneously neutralized so long as a pair of the evaluated bad brackets provide opposite signs when operated on by. On the other hand, if one or ore evaluated bad brackets all have the sae sign at x0 along the covector and all good brackets are annihilated by at x0, then the technique fails. Notice that while the exaples all had just two inputs, the technique applies without odification to hoogeneous degree-two systes with ore than two inputs ( >2). When other directions are involved, the neutralization ay encounter difficulties. Supposing that x is spanned by f ig k with k 2, the question of neutralization of potential obstructions becoes one of siultaneously transforing the atrices 9 so that they all have zeros on their diagonals. If the ranges of the atrices 9 are orthogonal, then the proble can be solved with a block diagonal transforation T, where each block appropriately transfors each 9. This procedure requires a straightforward odification of the construction of T. These interpretations are suarized in Table I. Finally, notice that the hoogeneity of f0 is not essential to the developent of neutralization via congruence transfor, the construction of the atrix 9 being sufficiently general that it applies to any nonlinear syste that is linear in control. For exaple, neutralization of potential obstructions fro type (1,2) brackets for the syste f0(x) = (x2x3;x1x3; sin 2 x1 0 sin 2 x2), f1 = (1; 0; 0), and f2 = (0; 1; 0) proceeds identically to that of the previous exaple with f0(x) = (x2x3;x1x3;x1 0 x2). Clearly the proposed technique provides for neutralization of potential obstructions fro type (1,2) brackets for these ore general systes. A generalization of neutralization via congruence transforation to inhoogeneous nonlinear systes would involve incorporation of the rich differential geoetry of nilpotent and higher order approxiations and foliations described, for exaple, in [9]. V. CONCLUSION We have presented a coplete characterization of STLC for the class of single-input, hoogeneous polynoial systes linear in control, where hoogeneous is used in the traditional sense. Specifically for odd-degree systes, STLC is equivalent to the LARC, while even-degree systes are never STLC except for the degenerate case of a scalar state. For ultiple-input hoogeneous systes linear in control, we have investigated neutralization of bad brackets with brackets of the sae type. The ethodology presented in this note provides a eans of neutralizing bad brackets of type (1,2). By consideration of the tensor generated fro the bracket structure [1; [1; f0]] applied to the direction containing a potential obstruction, we have reduced the question of neutralizing an obstruction to that of finding a congruence transfor that results in a atrix with all zeros along its diagonal. It is shown that such a transforation exists if and only if the atrix in question is indefinite. When this test is translated back to type (1,2) brackets, it has intuitive iplications, which are illustrated with several siple exaples. The ethodology presented is liited in its effectiveness by the fact that it reoves a potential obstruction only along a particular direction in the tangent space, although an extension to ultiple directions appears attainable. Although the neutralization via congruence transforation result has been presented in the context of hoogeneous systes, its developent does not rely on the hoogeneity of the drift vector field and, hence, applies to neutralization of type (1,2) brackets for any nonlinear syste that is linear in control. REFERENCES [1] H. G. Heres, On the structure of attainable sets for generalized differential equations and control systes, J. Diff. Equations, vol. 9, pp , Jan [2] H. J. Sussann and V. Jurdjevic, Controllability of nonlinear systes, J. Diff. Equations, vol. 12, pp , July [3] R. Herann and A. J. Krener, Nonlinear controllability and observability, IEEE Trans. Autoat. Contr., vol. AC-22, pp , Oct [4] H. J. Sussann, Lie brackets and local controllability: A sufficient condition for scalar-input systes, SIAM J. Control Opti., vol. 21, pp , Sept [5] G. Stefani, On the local controllability of a scalar input control syste, in Theory and Applications of Nonlinear Control Systes, C. I. Byrnes and A. Lindquist, Eds. New York: Elsevier, [6] M. Kawski, Higher-order sall-tie local controllability, in Nonlinear Controllability and Optial Control, H. J. Sussann, Ed. New York: Marcel Dekker, 1990, vol. 133, Pure and Applied Matheatics, ch. 14, pp [7] H. J. Sussann, General theore on local controllability, SIAM J. Control Opti., vol. 25, pp , Jan [8] J.-P. Serre, Lie Algebras and Lie Groups: Lectures Given at Harvard University, W. A. Benjain, Ed. Cabridge, MA, [9] H. G. Heres, Nilpotent and high-order approxiations of vector field systes, SIAM Rev., vol. 33, pp , June [10] A. A. Agračhev, Is it possible to recognize controllability in a finite nuber of differentiations?, in Open Probles in Matheatical Systes and Control Theory, V. D. Blondel, Ed. New York: Springer-Verlag, 1999, Counications and Control Engineering, ch. 4, pp [11] H. Nijeijer and A. J. van der Schaft, Nonlinear Dynaical Control Systes. New York: Springer-Verlag, 1990.

Multi-Dimensional Hegselmann-Krause Dynamics

Multi-Dimensional Hegselmann-Krause Dynamics Multi-Diensional Hegselann-Krause Dynaics A. Nedić Industrial and Enterprise Systes Engineering Dept. University of Illinois Urbana, IL 680 angelia@illinois.edu B. Touri Coordinated Science Laboratory

More information

Chapter 6 1-D Continuous Groups

Chapter 6 1-D Continuous Groups Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:

More information

Block designs and statistics

Block designs and statistics Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent

More information

Polygonal Designs: Existence and Construction

Polygonal Designs: Existence and Construction Polygonal Designs: Existence and Construction John Hegean Departent of Matheatics, Stanford University, Stanford, CA 9405 Jeff Langford Departent of Matheatics, Drake University, Des Moines, IA 5011 G

More information

Physics 139B Solutions to Homework Set 3 Fall 2009

Physics 139B Solutions to Homework Set 3 Fall 2009 Physics 139B Solutions to Hoework Set 3 Fall 009 1. Consider a particle of ass attached to a rigid assless rod of fixed length R whose other end is fixed at the origin. The rod is free to rotate about

More information

THE POLYNOMIAL REPRESENTATION OF THE TYPE A n 1 RATIONAL CHEREDNIK ALGEBRA IN CHARACTERISTIC p n

THE POLYNOMIAL REPRESENTATION OF THE TYPE A n 1 RATIONAL CHEREDNIK ALGEBRA IN CHARACTERISTIC p n THE POLYNOMIAL REPRESENTATION OF THE TYPE A n RATIONAL CHEREDNIK ALGEBRA IN CHARACTERISTIC p n SHEELA DEVADAS AND YI SUN Abstract. We study the polynoial representation of the rational Cherednik algebra

More information

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay

More information

. The univariate situation. It is well-known for a long tie that denoinators of Pade approxiants can be considered as orthogonal polynoials with respe

. The univariate situation. It is well-known for a long tie that denoinators of Pade approxiants can be considered as orthogonal polynoials with respe PROPERTIES OF MULTIVARIATE HOMOGENEOUS ORTHOGONAL POLYNOMIALS Brahi Benouahane y Annie Cuyt? Keywords Abstract It is well-known that the denoinators of Pade approxiants can be considered as orthogonal

More information

Physics 215 Winter The Density Matrix

Physics 215 Winter The Density Matrix Physics 215 Winter 2018 The Density Matrix The quantu space of states is a Hilbert space H. Any state vector ψ H is a pure state. Since any linear cobination of eleents of H are also an eleent of H, it

More information

Feature Extraction Techniques

Feature Extraction Techniques Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that

More information

Mechanics Physics 151

Mechanics Physics 151 Mechanics Physics 5 Lecture Oscillations (Chapter 6) What We Did Last Tie Analyzed the otion of a heavy top Reduced into -diensional proble of θ Qualitative behavior Precession + nutation Initial condition

More information

A new type of lower bound for the largest eigenvalue of a symmetric matrix

A new type of lower bound for the largest eigenvalue of a symmetric matrix Linear Algebra and its Applications 47 7 9 9 www.elsevier.co/locate/laa A new type of lower bound for the largest eigenvalue of a syetric atrix Piet Van Mieghe Delft University of Technology, P.O. Box

More information

Least Squares Fitting of Data

Least Squares Fitting of Data Least Squares Fitting of Data David Eberly, Geoetric Tools, Redond WA 98052 https://www.geoetrictools.co/ This work is licensed under the Creative Coons Attribution 4.0 International License. To view a

More information

Four-vector, Dirac spinor representation and Lorentz Transformations

Four-vector, Dirac spinor representation and Lorentz Transformations Available online at www.pelagiaresearchlibrary.co Advances in Applied Science Research, 2012, 3 (2):749-756 Four-vector, Dirac spinor representation and Lorentz Transforations S. B. Khasare 1, J. N. Rateke

More information

A1. Find all ordered pairs (a, b) of positive integers for which 1 a + 1 b = 3

A1. Find all ordered pairs (a, b) of positive integers for which 1 a + 1 b = 3 A. Find all ordered pairs a, b) of positive integers for which a + b = 3 08. Answer. The six ordered pairs are 009, 08), 08, 009), 009 337, 674) = 35043, 674), 009 346, 673) = 3584, 673), 674, 009 337)

More information

3.3 Variational Characterization of Singular Values

3.3 Variational Characterization of Singular Values 3.3. Variational Characterization of Singular Values 61 3.3 Variational Characterization of Singular Values Since the singular values are square roots of the eigenvalues of the Heritian atrices A A and

More information

i ij j ( ) sin cos x y z x x x interchangeably.)

i ij j ( ) sin cos x y z x x x interchangeably.) Tensor Operators Michael Fowler,2/3/12 Introduction: Cartesian Vectors and Tensors Physics is full of vectors: x, L, S and so on Classically, a (three-diensional) vector is defined by its properties under

More information

New Classes of Positive Semi-Definite Hankel Tensors

New Classes of Positive Semi-Definite Hankel Tensors Miniax Theory and its Applications Volue 017, No., 1 xxx New Classes of Positive Sei-Definite Hankel Tensors Qun Wang Dept. of Applied Matheatics, The Hong Kong Polytechnic University, Hung Ho, Kowloon,

More information

Chaotic Coupled Map Lattices

Chaotic Coupled Map Lattices Chaotic Coupled Map Lattices Author: Dustin Keys Advisors: Dr. Robert Indik, Dr. Kevin Lin 1 Introduction When a syste of chaotic aps is coupled in a way that allows the to share inforation about each

More information

The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters

The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters journal of ultivariate analysis 58, 96106 (1996) article no. 0041 The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Paraeters H. S. Steyn

More information

Reed-Muller Codes. m r inductive definition. Later, we shall explain how to construct Reed-Muller codes using the Kronecker product.

Reed-Muller Codes. m r inductive definition. Later, we shall explain how to construct Reed-Muller codes using the Kronecker product. Coding Theory Massoud Malek Reed-Muller Codes An iportant class of linear block codes rich in algebraic and geoetric structure is the class of Reed-Muller codes, which includes the Extended Haing code.

More information

RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE

RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE Proceedings of ICIPE rd International Conference on Inverse Probles in Engineering: Theory and Practice June -8, 999, Port Ludlow, Washington, USA : RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS

More information

Generalized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili,

Generalized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili, Australian Journal of Basic and Applied Sciences, 5(3): 35-358, 20 ISSN 99-878 Generalized AOR Method for Solving Syste of Linear Equations Davod Khojasteh Salkuyeh Departent of Matheatics, University

More information

Curious Bounds for Floor Function Sums

Curious Bounds for Floor Function Sums 1 47 6 11 Journal of Integer Sequences, Vol. 1 (018), Article 18.1.8 Curious Bounds for Floor Function Sus Thotsaporn Thanatipanonda and Elaine Wong 1 Science Division Mahidol University International

More information

The Fundamental Basis Theorem of Geometry from an algebraic point of view

The Fundamental Basis Theorem of Geometry from an algebraic point of view Journal of Physics: Conference Series PAPER OPEN ACCESS The Fundaental Basis Theore of Geoetry fro an algebraic point of view To cite this article: U Bekbaev 2017 J Phys: Conf Ser 819 012013 View the article

More information

lecture 36: Linear Multistep Mehods: Zero Stability

lecture 36: Linear Multistep Mehods: Zero Stability 95 lecture 36: Linear Multistep Mehods: Zero Stability 5.6 Linear ultistep ethods: zero stability Does consistency iply convergence for linear ultistep ethods? This is always the case for one-step ethods,

More information

Hybrid System Identification: An SDP Approach

Hybrid System Identification: An SDP Approach 49th IEEE Conference on Decision and Control Deceber 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Hybrid Syste Identification: An SDP Approach C Feng, C M Lagoa, N Ozay and M Sznaier Abstract The

More information

MA304 Differential Geometry

MA304 Differential Geometry MA304 Differential Geoetry Hoework 4 solutions Spring 018 6% of the final ark 1. The paraeterised curve αt = t cosh t for t R is called the catenary. Find the curvature of αt. Solution. Fro hoework question

More information

ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD

ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN

More information

Lecture #8-3 Oscillations, Simple Harmonic Motion

Lecture #8-3 Oscillations, Simple Harmonic Motion Lecture #8-3 Oscillations Siple Haronic Motion So far we have considered two basic types of otion: translation and rotation. But these are not the only two types of otion we can observe in every day life.

More information

G G G G G. Spec k G. G Spec k G G. G G m G. G Spec k. Spec k

G G G G G. Spec k G. G Spec k G G. G G m G. G Spec k. Spec k 12 VICTORIA HOSKINS 3. Algebraic group actions and quotients In this section we consider group actions on algebraic varieties and also describe what type of quotients we would like to have for such group

More information

NORMAL MATRIX POLYNOMIALS WITH NONSINGULAR LEADING COEFFICIENTS

NORMAL MATRIX POLYNOMIALS WITH NONSINGULAR LEADING COEFFICIENTS NORMAL MATRIX POLYNOMIALS WITH NONSINGULAR LEADING COEFFICIENTS NIKOLAOS PAPATHANASIOU AND PANAYIOTIS PSARRAKOS Abstract. In this paper, we introduce the notions of weakly noral and noral atrix polynoials,

More information

Distributed Subgradient Methods for Multi-agent Optimization

Distributed Subgradient Methods for Multi-agent Optimization 1 Distributed Subgradient Methods for Multi-agent Optiization Angelia Nedić and Asuan Ozdaglar October 29, 2007 Abstract We study a distributed coputation odel for optiizing a su of convex objective functions

More information

A Bernstein-Markov Theorem for Normed Spaces

A Bernstein-Markov Theorem for Normed Spaces A Bernstein-Markov Theore for Nored Spaces Lawrence A. Harris Departent of Matheatics, University of Kentucky Lexington, Kentucky 40506-0027 Abstract Let X and Y be real nored linear spaces and let φ :

More information

TEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES

TEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES TEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES S. E. Ahed, R. J. Tokins and A. I. Volodin Departent of Matheatics and Statistics University of Regina Regina,

More information

A Quantum Observable for the Graph Isomorphism Problem

A Quantum Observable for the Graph Isomorphism Problem A Quantu Observable for the Graph Isoorphis Proble Mark Ettinger Los Alaos National Laboratory Peter Høyer BRICS Abstract Suppose we are given two graphs on n vertices. We define an observable in the Hilbert

More information

Characterization of the Line Complexity of Cellular Automata Generated by Polynomial Transition Rules. Bertrand Stone

Characterization of the Line Complexity of Cellular Automata Generated by Polynomial Transition Rules. Bertrand Stone Characterization of the Line Coplexity of Cellular Autoata Generated by Polynoial Transition Rules Bertrand Stone Abstract Cellular autoata are discrete dynaical systes which consist of changing patterns

More information

Principal Components Analysis

Principal Components Analysis Principal Coponents Analysis Cheng Li, Bingyu Wang Noveber 3, 204 What s PCA Principal coponent analysis (PCA) is a statistical procedure that uses an orthogonal transforation to convert a set of observations

More information

3.8 Three Types of Convergence

3.8 Three Types of Convergence 3.8 Three Types of Convergence 3.8 Three Types of Convergence 93 Suppose that we are given a sequence functions {f k } k N on a set X and another function f on X. What does it ean for f k to converge to

More information

Quantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search

Quantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search Quantu algoriths (CO 781, Winter 2008) Prof Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search ow we begin to discuss applications of quantu walks to search algoriths

More information

A Note on Scheduling Tall/Small Multiprocessor Tasks with Unit Processing Time to Minimize Maximum Tardiness

A Note on Scheduling Tall/Small Multiprocessor Tasks with Unit Processing Time to Minimize Maximum Tardiness A Note on Scheduling Tall/Sall Multiprocessor Tasks with Unit Processing Tie to Miniize Maxiu Tardiness Philippe Baptiste and Baruch Schieber IBM T.J. Watson Research Center P.O. Box 218, Yorktown Heights,

More information

NON-COMMUTATIVE GRÖBNER BASES FOR COMMUTATIVE ALGEBRAS

NON-COMMUTATIVE GRÖBNER BASES FOR COMMUTATIVE ALGEBRAS PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volue 126, Nuber 3, March 1998, Pages 687 691 S 0002-9939(98)04229-4 NON-COMMUTATIVE GRÖBNER BASES FOR COMMUTATIVE ALGEBRAS DAVID EISENBUD, IRENA PEEVA,

More information

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either

More information

A := A i : {A i } S. is an algebra. The same object is obtained when the union in required to be disjoint.

A := A i : {A i } S. is an algebra. The same object is obtained when the union in required to be disjoint. 59 6. ABSTRACT MEASURE THEORY Having developed the Lebesgue integral with respect to the general easures, we now have a general concept with few specific exaples to actually test it on. Indeed, so far

More information

Decentralized Adaptive Control of Nonlinear Systems Using Radial Basis Neural Networks

Decentralized Adaptive Control of Nonlinear Systems Using Radial Basis Neural Networks 050 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 44, NO., NOVEMBER 999 Decentralized Adaptive Control of Nonlinear Systes Using Radial Basis Neural Networks Jeffrey T. Spooner and Kevin M. Passino Abstract

More information

A note on the realignment criterion

A note on the realignment criterion A note on the realignent criterion Chi-Kwong Li 1, Yiu-Tung Poon and Nung-Sing Sze 3 1 Departent of Matheatics, College of Willia & Mary, Williasburg, VA 3185, USA Departent of Matheatics, Iowa State University,

More information

Linear Algebra (I) Yijia Chen. linear transformations and their algebraic properties. 1. A Starting Point. y := 3x.

Linear Algebra (I) Yijia Chen. linear transformations and their algebraic properties. 1. A Starting Point. y := 3x. Linear Algebra I) Yijia Chen Linear algebra studies Exaple.. Consider the function This is a linear function f : R R. linear transforations and their algebraic properties.. A Starting Point y := 3x. Geoetrically

More information

MULTIPLAYER ROCK-PAPER-SCISSORS

MULTIPLAYER ROCK-PAPER-SCISSORS MULTIPLAYER ROCK-PAPER-SCISSORS CHARLOTTE ATEN Contents 1. Introduction 1 2. RPS Magas 3 3. Ites as a Function of Players and Vice Versa 5 4. Algebraic Properties of RPS Magas 6 References 6 1. Introduction

More information

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators Suppleentary Inforation for Design of Bending Multi-Layer Electroactive Polyer Actuators Bavani Balakrisnan, Alek Nacev, and Elisabeth Sela University of Maryland, College Park, Maryland 074 1 Analytical

More information

A note on the multiplication of sparse matrices

A note on the multiplication of sparse matrices Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani

More information

On Certain C-Test Words for Free Groups

On Certain C-Test Words for Free Groups Journal of Algebra 247, 509 540 2002 doi:10.1006 jabr.2001.9001, available online at http: www.idealibrary.co on On Certain C-Test Words for Free Groups Donghi Lee Departent of Matheatics, Uni ersity of

More information

LEADING COEFFICIENTS AND THE MULTIPLICITY OF KNOWN ROOTS

LEADING COEFFICIENTS AND THE MULTIPLICITY OF KNOWN ROOTS LEADING COEFFICIENTS AND THE MULTIPLICITY OF KNOWN ROOTS GREGORY J. CLARK AND JOSHUA N. COOPER Abstract. We show that a onic univariate polynoial over a field of characteristic zero, with k distinct non-zero

More information

Numerically repeated support splitting and merging phenomena in a porous media equation with strong absorption. Kenji Tomoeda

Numerically repeated support splitting and merging phenomena in a porous media equation with strong absorption. Kenji Tomoeda Journal of Math-for-Industry, Vol. 3 (C-), pp. Nuerically repeated support splitting and erging phenoena in a porous edia equation with strong absorption To the eory of y friend Professor Nakaki. Kenji

More information

RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS

RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS BIT Nuerical Matheatics 43: 459 466, 2003. 2003 Kluwer Acadeic Publishers. Printed in The Netherlands 459 RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS V. SIMONCINI Dipartiento di

More information

Topic 5a Introduction to Curve Fitting & Linear Regression

Topic 5a Introduction to Curve Fitting & Linear Regression /7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline

More information

Kinetic Theory of Gases: Elementary Ideas

Kinetic Theory of Gases: Elementary Ideas Kinetic Theory of Gases: Eleentary Ideas 17th February 2010 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion

More information

Fast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials

Fast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials Fast Montgoery-like Square Root Coputation over GF( ) for All Trinoials Yin Li a, Yu Zhang a, a Departent of Coputer Science and Technology, Xinyang Noral University, Henan, P.R.China Abstract This letter

More information

13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization

13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization 3 Haronic oscillator revisited: Dirac s approach and introduction to Second Quantization. Dirac cae up with a ore elegant way to solve the haronic oscillator proble. We will now study this approach. The

More information

PHYSICS 110A : CLASSICAL MECHANICS MIDTERM EXAM #2

PHYSICS 110A : CLASSICAL MECHANICS MIDTERM EXAM #2 PHYSICS 110A : CLASSICAL MECHANICS MIDTERM EXAM #2 [1] Two blocks connected by a spring of spring constant k are free to slide frictionlessly along a horizontal surface, as shown in Fig. 1. The unstretched

More information

ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE

ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE CHRISTOPHER J. HILLAR Abstract. A long-standing conjecture asserts that the polynoial p(t = Tr(A + tb ] has nonnegative coefficients whenever is

More information

Interactive Markov Models of Evolutionary Algorithms

Interactive Markov Models of Evolutionary Algorithms Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary

More information

Lecture 21. Interior Point Methods Setup and Algorithm

Lecture 21. Interior Point Methods Setup and Algorithm Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and

More information

On Lotka-Volterra Evolution Law

On Lotka-Volterra Evolution Law Advanced Studies in Biology, Vol. 3, 0, no. 4, 6 67 On Lota-Volterra Evolution Law Farruh Muhaedov Faculty of Science, International Islaic University Malaysia P.O. Box, 4, 570, Kuantan, Pahang, Malaysia

More information

The Simplex Method is Strongly Polynomial for the Markov Decision Problem with a Fixed Discount Rate

The Simplex Method is Strongly Polynomial for the Markov Decision Problem with a Fixed Discount Rate The Siplex Method is Strongly Polynoial for the Markov Decision Proble with a Fixed Discount Rate Yinyu Ye April 20, 2010 Abstract In this note we prove that the classic siplex ethod with the ost-negativereduced-cost

More information

Work, Energy and Momentum

Work, Energy and Momentum Work, Energy and Moentu Work: When a body oves a distance d along straight line, while acted on by a constant force of agnitude F in the sae direction as the otion, the work done by the force is tered

More information

Algebraic Montgomery-Yang problem: the log del Pezzo surface case

Algebraic Montgomery-Yang problem: the log del Pezzo surface case c 2014 The Matheatical Society of Japan J. Math. Soc. Japan Vol. 66, No. 4 (2014) pp. 1073 1089 doi: 10.2969/jsj/06641073 Algebraic Montgoery-Yang proble: the log del Pezzo surface case By DongSeon Hwang

More information

Chapter 1: Basics of Vibrations for Simple Mechanical Systems

Chapter 1: Basics of Vibrations for Simple Mechanical Systems Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,

More information

Kinetic Theory of Gases: Elementary Ideas

Kinetic Theory of Gases: Elementary Ideas Kinetic Theory of Gases: Eleentary Ideas 9th February 011 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion of

More information

#A52 INTEGERS 10 (2010), COMBINATORIAL INTERPRETATIONS OF BINOMIAL COEFFICIENT ANALOGUES RELATED TO LUCAS SEQUENCES

#A52 INTEGERS 10 (2010), COMBINATORIAL INTERPRETATIONS OF BINOMIAL COEFFICIENT ANALOGUES RELATED TO LUCAS SEQUENCES #A5 INTEGERS 10 (010), 697-703 COMBINATORIAL INTERPRETATIONS OF BINOMIAL COEFFICIENT ANALOGUES RELATED TO LUCAS SEQUENCES Bruce E Sagan 1 Departent of Matheatics, Michigan State University, East Lansing,

More information

The Methods of Solution for Constrained Nonlinear Programming

The Methods of Solution for Constrained Nonlinear Programming Research Inventy: International Journal Of Engineering And Science Vol.4, Issue 3(March 2014), PP 01-06 Issn (e): 2278-4721, Issn (p):2319-6483, www.researchinventy.co The Methods of Solution for Constrained

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I Contents 1. Preliinaries 2. The ain result 3. The Rieann integral 4. The integral of a nonnegative

More information

Uniform Approximation and Bernstein Polynomials with Coefficients in the Unit Interval

Uniform Approximation and Bernstein Polynomials with Coefficients in the Unit Interval Unifor Approxiation and Bernstein Polynoials with Coefficients in the Unit Interval Weiang Qian and Marc D. Riedel Electrical and Coputer Engineering, University of Minnesota 200 Union St. S.E. Minneapolis,

More information

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co

More information

Divisibility of Polynomials over Finite Fields and Combinatorial Applications

Divisibility of Polynomials over Finite Fields and Combinatorial Applications Designs, Codes and Cryptography anuscript No. (will be inserted by the editor) Divisibility of Polynoials over Finite Fields and Cobinatorial Applications Daniel Panario Olga Sosnovski Brett Stevens Qiang

More information

Using EM To Estimate A Probablity Density With A Mixture Of Gaussians

Using EM To Estimate A Probablity Density With A Mixture Of Gaussians Using EM To Estiate A Probablity Density With A Mixture Of Gaussians Aaron A. D Souza adsouza@usc.edu Introduction The proble we are trying to address in this note is siple. Given a set of data points

More information

Linear recurrences and asymptotic behavior of exponential sums of symmetric boolean functions

Linear recurrences and asymptotic behavior of exponential sums of symmetric boolean functions Linear recurrences and asyptotic behavior of exponential sus of syetric boolean functions Francis N. Castro Departent of Matheatics University of Puerto Rico, San Juan, PR 00931 francis.castro@upr.edu

More information

A Study on B-Spline Wavelets and Wavelet Packets

A Study on B-Spline Wavelets and Wavelet Packets Applied Matheatics 4 5 3-3 Published Online Noveber 4 in SciRes. http://www.scirp.org/ournal/a http://dx.doi.org/.436/a.4.5987 A Study on B-Spline Wavelets and Wavelet Pacets Sana Khan Mohaad Kaliuddin

More information

Exact tensor completion with sum-of-squares

Exact tensor completion with sum-of-squares Proceedings of Machine Learning Research vol 65:1 54, 2017 30th Annual Conference on Learning Theory Exact tensor copletion with su-of-squares Aaron Potechin Institute for Advanced Study, Princeton David

More information

Optimal Jamming Over Additive Noise: Vector Source-Channel Case

Optimal Jamming Over Additive Noise: Vector Source-Channel Case Fifty-first Annual Allerton Conference Allerton House, UIUC, Illinois, USA October 2-3, 2013 Optial Jaing Over Additive Noise: Vector Source-Channel Case Erah Akyol and Kenneth Rose Abstract This paper

More information

Lecture 13 Eigenvalue Problems

Lecture 13 Eigenvalue Problems Lecture 13 Eigenvalue Probles MIT 18.335J / 6.337J Introduction to Nuerical Methods Per-Olof Persson October 24, 2006 1 The Eigenvalue Decoposition Eigenvalue proble for atrix A: Ax = λx with eigenvalues

More information

ON REGULARITY, TRANSITIVITY, AND ERGODIC PRINCIPLE FOR QUADRATIC STOCHASTIC VOLTERRA OPERATORS MANSOOR SABUROV

ON REGULARITY, TRANSITIVITY, AND ERGODIC PRINCIPLE FOR QUADRATIC STOCHASTIC VOLTERRA OPERATORS MANSOOR SABUROV ON REGULARITY TRANSITIVITY AND ERGODIC PRINCIPLE FOR QUADRATIC STOCHASTIC VOLTERRA OPERATORS MANSOOR SABUROV Departent of Coputational & Theoretical Sciences Faculty of Science International Islaic University

More information

Module #1: Units and Vectors Revisited. Introduction. Units Revisited EXAMPLE 1.1. A sample of iron has a mass of mg. How many kg is that?

Module #1: Units and Vectors Revisited. Introduction. Units Revisited EXAMPLE 1.1. A sample of iron has a mass of mg. How many kg is that? Module #1: Units and Vectors Revisited Introduction There are probably no concepts ore iportant in physics than the two listed in the title of this odule. In your first-year physics course, I a sure that

More information

HESSIAN MATRICES OF PENALTY FUNCTIONS FOR SOLVING CONSTRAINED-OPTIMIZATION PROBLEMS

HESSIAN MATRICES OF PENALTY FUNCTIONS FOR SOLVING CONSTRAINED-OPTIMIZATION PROBLEMS R 702 Philips Res. Repts 24, 322-330, 1969 HESSIAN MATRICES OF PENALTY FUNCTIONS FOR SOLVING CONSTRAINED-OPTIMIZATION PROBLEMS by F. A. LOOTSMA Abstract This paper deals with the Hessian atrices of penalty

More information

arxiv:math/ v1 [math.nt] 15 Jul 2003

arxiv:math/ v1 [math.nt] 15 Jul 2003 arxiv:ath/0307203v [ath.nt] 5 Jul 2003 A quantitative version of the Roth-Ridout theore Toohiro Yaada, 606-8502, Faculty of Science, Kyoto University, Kitashirakawaoiwakecho, Sakyoku, Kyoto-City, Kyoto,

More information

Vulnerability of MRD-Code-Based Universal Secure Error-Correcting Network Codes under Time-Varying Jamming Links

Vulnerability of MRD-Code-Based Universal Secure Error-Correcting Network Codes under Time-Varying Jamming Links Vulnerability of MRD-Code-Based Universal Secure Error-Correcting Network Codes under Tie-Varying Jaing Links Jun Kurihara KDDI R&D Laboratories, Inc 2 5 Ohara, Fujiino, Saitaa, 356 8502 Japan Eail: kurihara@kddilabsjp

More information

Tutorial Exercises: Incorporating constraints

Tutorial Exercises: Incorporating constraints Tutorial Exercises: Incorporating constraints 1. A siple pendulu of length l ass is suspended fro a pivot of ass M that is free to slide on a frictionless wire frae in the shape of a parabola y = ax. The

More information

2.003 Engineering Dynamics Problem Set 2 Solutions

2.003 Engineering Dynamics Problem Set 2 Solutions .003 Engineering Dynaics Proble Set Solutions This proble set is priarily eant to give the student practice in describing otion. This is the subject of kineatics. It is strongly recoended that you study

More information

CONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. IX Uncertainty Models For Robustness Analysis - A. Garulli, A. Tesi and A. Vicino

CONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. IX Uncertainty Models For Robustness Analysis - A. Garulli, A. Tesi and A. Vicino UNCERTAINTY MODELS FOR ROBUSTNESS ANALYSIS A. Garulli Dipartiento di Ingegneria dell Inforazione, Università di Siena, Italy A. Tesi Dipartiento di Sistei e Inforatica, Università di Firenze, Italy A.

More information

Generalized eigenfunctions and a Borel Theorem on the Sierpinski Gasket.

Generalized eigenfunctions and a Borel Theorem on the Sierpinski Gasket. Generalized eigenfunctions and a Borel Theore on the Sierpinski Gasket. Kasso A. Okoudjou, Luke G. Rogers, and Robert S. Strichartz May 26, 2006 1 Introduction There is a well developed theory (see [5,

More information

Fixed-to-Variable Length Distribution Matching

Fixed-to-Variable Length Distribution Matching Fixed-to-Variable Length Distribution Matching Rana Ali Ajad and Georg Böcherer Institute for Counications Engineering Technische Universität München, Gerany Eail: raa2463@gail.co,georg.boecherer@tu.de

More information

Physics 2210 Fall smartphysics 20 Conservation of Angular Momentum 21 Simple Harmonic Motion 11/23/2015

Physics 2210 Fall smartphysics 20 Conservation of Angular Momentum 21 Simple Harmonic Motion 11/23/2015 Physics 2210 Fall 2015 sartphysics 20 Conservation of Angular Moentu 21 Siple Haronic Motion 11/23/2015 Exa 4: sartphysics units 14-20 Midter Exa 2: Day: Fri Dec. 04, 2015 Tie: regular class tie Section

More information

Lecture 9 November 23, 2015

Lecture 9 November 23, 2015 CSC244: Discrepancy Theory in Coputer Science Fall 25 Aleksandar Nikolov Lecture 9 Noveber 23, 25 Scribe: Nick Spooner Properties of γ 2 Recall that γ 2 (A) is defined for A R n as follows: γ 2 (A) = in{r(u)

More information

Page 1 Lab 1 Elementary Matrix and Linear Algebra Spring 2011

Page 1 Lab 1 Elementary Matrix and Linear Algebra Spring 2011 Page Lab Eleentary Matri and Linear Algebra Spring 0 Nae Due /03/0 Score /5 Probles through 4 are each worth 4 points.. Go to the Linear Algebra oolkit site ransforing a atri to reduced row echelon for

More information

Finite fields. and we ve used it in various examples and homework problems. In these notes I will introduce more finite fields

Finite fields. and we ve used it in various examples and homework problems. In these notes I will introduce more finite fields Finite fields I talked in class about the field with two eleents F 2 = {, } and we ve used it in various eaples and hoework probles. In these notes I will introduce ore finite fields F p = {,,...,p } for

More information

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control An Extension to the Tactical Planning Model for a Job Shop: Continuous-Tie Control Chee Chong. Teo, Rohit Bhatnagar, and Stephen C. Graves Singapore-MIT Alliance, Nanyang Technological Univ., and Massachusetts

More information

12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015

12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015 18.354J Nonlinear Dynaics II: Continuu Systes Lecture 12 Spring 2015 12 Towards hydrodynaic equations The previous classes focussed on the continuu description of static (tie-independent) elastic systes.

More information

Lower Bounds for Quantized Matrix Completion

Lower Bounds for Quantized Matrix Completion Lower Bounds for Quantized Matrix Copletion Mary Wootters and Yaniv Plan Departent of Matheatics University of Michigan Ann Arbor, MI Eail: wootters, yplan}@uich.edu Mark A. Davenport School of Elec. &

More information

A Simple Regression Problem

A Simple Regression Problem A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where

More information

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving

More information

c 2010 Society for Industrial and Applied Mathematics

c 2010 Society for Industrial and Applied Mathematics SIAM J MATRIX ANAL APPL Vol 31, No 4, pp 1848 1872 c 21 Society for Industrial and Applied Matheatics STOCHASTIC GALERKIN MATRICES OLIVER G ERNST AND ELISABETH ULLMANN Abstract We investigate the structural,

More information