Reduced-Eort Control Laws for. Panagiotis Tsiotras and Jihao Luo y. University of Virginia, Charlottesville, VA

Size: px
Start display at page:

Download "Reduced-Eort Control Laws for. Panagiotis Tsiotras and Jihao Luo y. University of Virginia, Charlottesville, VA"

Transcription

1 Journal of Guiance, Control, an Dynamics ol. XX, No. X, July-August XXXX Reuce-Eort Control Laws for Uneractuate Rigi Spacecraft Panagiotis Tsiotras an Jihao Luo y University of irginia, Charlottesville, A 9-44 Abstract Recent results show that a nonsmooth, time-invariant feeback control law can be use to rotate an axisymmetric rigi spacecraft to the zero equilibrium using only two control torques. This metho, however, may require a signicant amountofcontrol eort, especially for initial conitions close to an equilibrium manifol corresponing to rotations about the unactuate principal axis. In this paper a control law is propose which reuces the control eort require to perform rest-to-rest maneuvers for initial conitions close to this equilibrium manifol. Specically, the phase space of the system is ivie into two parts, one corresponing to initial conitions proucing large control eort (the \ba" region) an the other corresponing to initial conitions proucing small control signals (the \goo" region). The propose control law then reners this unesirable equilibrium manifol unstable riving the trajectories of the close-loop system into the \goo" region, where the original control law is subsequently use. Numerical simulations inicate reuction of the control magnitue at the orer of 8-9 % for initial conitions close to the equilibrium manifol. Introuction The problem of stabilization of a rigi boy using less than three control inputs has receive consierable attention in the recent literature. Both the problems of the stabilization of the ynamics, an the stabilization of the kinematics have been treate in the literature {6. The stabilization problem of the complete system, i.e., the ynamics an the kinematics, has been aresse in Refs. 7{. The attitue stabilization of an axially symmetric rigi boy using two inepenent control torques was stuie by Krishnan, et al. 8,9 an Tsiotras et al.. If the uncontrolle principal axis is not the axis of symmetry, the system is strongly accessible an small time locally controllable 9. When the uncontrolle axis coincies with the axis of symmetry, the complete system fails to be controllable or even accessible. However, the system equations are strongly accessible an small time locally controllable in the case of zero spin rate. A nonlinear control approach was evelope in Ref. 8, which achieves arbitrary reorientation for this restricte case. In Refs. 4-5, the authors presente a new formulation of the attitue kinematics which was use in Ref. to solve the same problem avoiing the successive switchings of Ref. 8. References 8 an treate the axisymmetric case, whereas the non-symmetric case was ealt with in Refs. { an 6. In this paper, a moication of the control law presente in Ref. for the attitue stabilization of an axi-symmetric rigi boy using two inepenent control torques is propose. Because Brockett's necessary conition for smooth stabilizability is not satise for this system, any stabilizing (time-invariant) control law is necessarily nonsmooth. (Stabilizing time-varying Assistant Professor, Department of Mechanical, Aerospace an Nuclear Engineering. Member AIAA. y Grauate Stuent, Department of Mechanical, Aerospace an Nuclear Engineering. smooth control laws may still exist, however.) This nonsmoothness is evient in Ref. in the form of the nonierentiability of the control law at the origin. Because of the singularity at the origin, this control law may saturate the actuators, especially for initial conitions close to the equilibrium manifol. It is therefore esirable to moify the control law of Ref. to reuce the require control signals. Compare to the control law in Ref., the moie control law propose here remeies this large control input problem by riving the trajectories of the close-loop system away from the singular equilibrium manifol, towars a region in the state space where the \high authority" part of the control input remains small an boune. The proceure is simple an can be easily valiate using phase portrait consierations. A numerical example illustrates the control eort improvement using the new control law. The Uneractuate Spacecraft The ynamics of a rigi spacecraft with two controls can be written as _! = a!! + u (a) _! = a!! + u (b) _! = a!! (c) where a i are the inertia parameters satisfying a + a + a + a a a =. Herewe assume a boy-xe reference frame along the principal axes of inertia. Equations () escribe an uneractuate spacecraft with no control authority about the r principal axis. Notice in this case,! canbecontrolle only inirectly through juicious choice of the time histories of! (t) an! (t). In case of an axi-symmetric boy (about the -axis), a =ana = ;a = a an Eqs. () reuce

2 to ^ ^ i = i _! = a!! + u (a) _! = ;a!! + u (b) _! = (c) c z θ where! () =! is constant. Introucing the complex variables! =! +i! an u = u +iu (with i = p ;) the previous equations can be written as ^ b _! = ;ia!! + u () u^ Kinematics of the Attitue Motion a b ^ i The orientation of a rigi spacecraft can be specie using various parameterizations, for example, Eulerian Angles, Euler Parameters, Cayley-Rorigues Parameters, Cayley-Klein parameters, etc 7. Recently, a new parameterization using a pair of complex an real coorinates was introuce 4,5. Accoring to these results, the relative orientation between twogiven reference frames can be represente by two rotations, one corresponing to the real coorinate (z) an the other corresponing to the complex coorinate (w). Specifically, one can align an (inertial) reference frame to a boy-xe frame by rst performing an initial rotation of magnitue z about, say, the inertial -axis an then performing a secon rotation to move theintermei- ate -axis to the boy -axis. The secon rotation can be completely characterize by the complex coorinate w = w + i w. It is a rotation of magnitue = arccos about the unit vector ^u = w + w jwj ;jwj +jwj (4) i( w ; w) ^i + ^i jwj (5) This situation is epicte in Fig., where (^i ^i ^i ) is the intermeiate reference frame resulting from the rotation z about the inertial ^i axis, an where (a b c) enote the coorinates of the unit vector along the ^i axis in the boy frame, that is, ^i = a ^b + b ^b + c ^b (6) It can be shown 5 that the coorinates of the ^b axis in the ^i frame are also relate to (a b c) ^b = ;a^i ; b^i + c^i (7) With this notation, w represents the stereographic coorinates corresponing to the unit vector (a b c) ene by 8,5 w = b ; ia (8) +c Alternatively, the equations a = i (w ; w) +jwj b = w + w +jwj c = ;jwj +jwj (9) ^ i Fig. Attitue representation using (w z) coorinates. an can be use to n a b an c once w is known. Here jj enotes the absolute value of a complex number, i.e., w w = jwj, w C. The kinematic equations, which provie the geometric constraints of the motion an relate the rates of the kinematic parameters w an z to the angular velocity vector, can be written as follows,5 _w = ;i! w +! +! w (a) _z =! + Im(! w) (b) where! =! + i! an w = w + iw. Notice these equations take the convenient form t jwj = ( + jwj )Re(! w) (a) _z =! + Im(! w) (b) where bar enotes complex conjugate, Re() an Im() enote the real an imaginary parts of a complex number respectively. In Eq. (b) only the imaginary part of the prouct! w appears, while in Eq. (a) only the real part appears. This uality (or anti-symmetry) of Eqs. (a) an (b) is esirable an can be use to erive stabilizing control laws for the kinematics escribe by Eqs. (). Clearly, w = if an only if jwj =, an stabilization of the system in Eqs. () is equivalent to stabilization of the system in Eqs. (). References,9,8 inicate that the coorinates (w z) oer some signicant avantages for attitue analysis an control problems. 4 Problem Statement Consier an axi-symmetric boy with the applie torque vector in the plane which is perpenicular to the symmetry axis. In such a case, the system is escribe by Eqs. () an thus! remains constant. If initially! () 6=,nocontrol input can bring the system to the equilibrium. The system is not controllable

3 to the equilibrium, but it is controllable to the submanifol! = w = in the (!! w z)-space. For a more etaile iscussion on this issue, refer to Refs. 8{. Therefore, for an axi-symmetric boy, actively controlle rotation to the equilibrium for the system in Eqs. ()-() makes sense only if!. In this case, the system equations simplify to _! = u (a) _w =! +! w (b) _z = Im(! w) (c) This system can be stabilize to the origin, but any time-invariant stabilizing control law has to be necessarily nonsmooth, since Eqs. () fail Brockett's necessary conition for smooth stabilizability. Therefore we concentrate on using nonsmooth (albeit time-invariant) stabilizers for this system. Equations () represent a system in cascae form, with the kinematics (b)-(c) being the riven subsystem an the ynamics (a) being the riving subsystem. Methoology in Ref. use this structure to erive a non-smooth control law to stabilize Eqs. (). In essence, the controller esign consists of a two-step process. In the rst step, only stabilization of the kinematics is aresse, with the angular velocity treate as the control input. In the secon step, the control torque u is chosen to shape the esire velocity prole. Since the angular velocity in the rst step is (necessarily) a nonsmooth function of w an z, caution shoul be exercise when implementing this angular velocity in the secon step. The nonsmooth controller of Ref., along with its potential rawbacks, is summarize in the next section. 5 A Nonsmooth Controller for the Kinematics In Ref. a nonsmooth control law was propose for the kinematic system escribe by _w =! +! w (a) _z = Im(! w) (b) an was later implemente though the integrator in Eq. (a). The propose control law inref. was motivate by the ecoupling of these equations with respect to the prouct! w, as evient from the iscussion following Eqs. (). This control law is given by! = ;w ; i z w (4) where >= >. With this control law, the close loop system in terms of jwj an z is given by jwj t = ;( + jwj )jwj (5a) _z = ;z (5b) which is globally exponentially stable. As can be easily inferre by observing Eqs. (4) an (), the rst term in the control law (4) has an eect only on the ierential equation for w, whereas the secon term in Eq. (4) has an eect only on the ierential equation for z. Moreover, the secon term in Eq. (4) is a nonsmooth function of z an w. The main isavantage of the control law in Eq. (4) is that the last term, which involves the ratio z= w, may become unboune without careful choice of the gains. The previously impose gain conition >= ensures that the rate of ecay ofz is at least as large as the rate of ecay of w, such that their ratio remains boune. Actually, one can easily establish from Eqs. (5) that for >=, along the solutions of the system, z= w! as t!. Introucing the variable v = jwj, the system in Eqs. (5) takes the form _v = ;( + v)v (6a) _z = ;z (6b) This is a system which evolves on IR + IR. Typical trajectories an the vector el of the close-loop system in Eq. (6) for =an = are shown in Fig.. (Since z oes not change sign it suces to plot only the z> case.) Fig. Phase portrait of system in Eqs. (6). Although in Eq. (4), the ratio z= w, an hence the control eort!, remains boune by properchoice of control gains, the control input! may take large values in the region where w is small. From Eq. (5a), jw(t)j jw()j for all t an for small initial conitions w(), the control law may use a substantial amount of energy, especially in regions where jzj is large. In Fig., for example, the region which is close to the z axis is clearly unesirable as far as control expeniture is concerne. Moication of the control law in Eq. (4), such that the vector el close to the z axis points away from this axis, is highly esirable. In short, the iea is to ivie the (z v) phase space into two regions accoring to the value of the ratio = z jwj = z v (7)

4 This ratio is a irect inication of the relative magnitue between z an w. This ratio shoul be kept small in orer to avoi high control eort. Hence, if initially the states are in an unesirable region where attains large values, the feeback control strategy shoul rive the trajectories to a \safe" region in the state space where remains relatively small. Without loss of generality, choose as unesirable the region where jj >, leaving jj as the esirable region. These two regions, enote by D an D respectively, are therefore ene by D = f(z v) IR IR + : > jj > g (8a) D = f(z v) IR IR + : jj g (8b) These two regions are shown in Fig D D η = + η = D Fig. Regions D an D in (z v) phase space. 6 Main Results The propose moication to the control law in Eq. (4) is simple. Positive feeback is use for v when the trajectory is in region D, while z is ecreasing. This change will make the manifol v = (equivalently, w = ) unstable an the trajectories will move towars the region D an subsequently stay there. The control law in region D is essentially the same as in Eq. (4). Notice that, by enition, insie the region D we have jj, an since jzj=j wj = jjjwj we can ensure that! will not take excessive values as long as the trajectories remain in D. These statements will be mae more precise in the sequel. 6. Propose Control Law for Kinematics The propose control law for the system in Eqs. () is ene by! = ;()w ; i() z w where () an () are smooth functions satisfying (9) ; c () < () < c 8 (z v) D (a) () c c () c 8 (z v) D (b) an < c < c. One possible choice is, for example, () = c arctan ; ( ; ) (a) () = c arctan ; ( ; ) + c From Eqs. (), an are boune as (b) ; c () c an () c () for all IR. Moreover, notice that () < () for all (z v) D. The next theorem gives the main result of the paper. Theorem 6. Consier the system in Eqs. () an let the control law be as in Eqs. (9)-() with < c < c. Then for initial conitions (z() w()) IR (Cnfg), the following properties hol: (i) w(t) 6= 8 t. (ii) the trajectory(z() w()) is boune an lim (z(t) w(t)) = () t! (iii) the control law!() is boune an it has a boune erivative. With the control law in Eq. (9), the close-loop system takes the form _v = ;()( + v)v (4a) _z = ;()z (4b) where v = jwj an as in Eq. (7). From Eq. () we have that z ecays monotonically for all initial conitions, whereas v increases in the region D an ecreases in D. The result is that the trajectories of Eqs. (4) ten to D an then to the origin, as require. Before proving Theorem 6. we nee to establish the following two lemmas. Lemma 6. The region D is invariant for the system in Eqs. (4). Proof. The bounary of the set D is given by the two lines = (cf. Fig. ). On the bounary of D the feeback gains are () = an () = c=. The vector el on the bounary of D is therefore _v = (5a) _z = ; c z (5b) which points into the interior of D. Therefore trajectories in D cannot escape this region an thus it is invariant for the close-loop system in Eqs. (4). This lemma establishes that for initial conitions in D, the trajectories of the close-loop system remain in D for all times. Equivalently, if at some time t the trajectory enters D,itstays in D for all t t. Figure shows the vector el on the bounary of D. 4

5 Lemma 6. Consier the system in Eqs. (4). For all initial conitions (z v) D the trajectories enter the region D in nite time. As long as (z v) D, from Eq. (a) is Proof. boune as () < c=. This implies that z is boune. Actually, jz(t)jjz()j for all t. Note that z oes not change sign for all t. Without loss of generality, assume that z() (the case z() being similar). If (z() v()) D then, by enition () >. The erivative of in D is then _ = _z v ; z v _v = ;() + ()( + v) ;() (6) since () < anv> hence is boune in D. Let cl D enote the closure of D in IR, that is, cl D = D [f(z v) IR IR + : jj =g [f(z v) IR IR + : v =g (7) Then it is an easy exercise to show that _ 6= for all (z v) cl D nf( )g. Hence there exists >such that _ <; in D an consequently, monotonically ecreases. Thus, every trajectory starting in D will leave this set an enter D in nite time. Notice that the set f(z v) IR IR + : v = an z 6= g is an unstable manifol for the close-loop system. Figure shows the vector el on the bounary of D. The following corollary follows irectly from Lemmas 6. an 6.. Corollary 6. Consier the system in Eqs. (4). For all initial conitions (z() v()) IR (IR +nfg), is boune for all t. We are now reay to give the proof of Theorem 6.. Proof. [Theorem 6.] From Eqs. (4a) an () we have that _v ; c( + v)v (8) where c >. The solution of the ierential equation _x = ; c( + x)x x() = x > (9) is given by x(t) = () c e c t ; where c =(x +)=x. Clearly, x(t) 6= for all t an lim t! x(t) =. Therefore v() is boune below by the solutions of the ierential equation (9) subject to initial conition x = v(). Hence, jw(t)j 6= for all t an w() approaches the origin asymptotically. We now show that lim t!(z(t) v(t)) =. If (z() v()) D then accoring to Lemma 6. we have that (z(t) v(t)) D for all t an D is an invariant set for the close-loop system. Consier now the positive enite, raially unboune function :IR IR +! IR + given by (z v)= v + z 8 (z v) D () The erivative of along the trajectories of (4) is _ = ;()( + v)v ; ()z 8(z v) D () therefore, the trajectories are boune in D. Moreover, _ = if an only if ()( + v)v + ()z =. Using the enitions of () an() ind an recalling that v, one establishes that the last equality is not satise in D unless z = v =. By LaSalle's theorem, lim t!(z(t) v(t)) =, for all initial conitions in D.To nish the proof, recall from Lemma 6. that if (z() v()) D then jzj is boune by jz()j an there exist a time t > suchthat(z(t ) v(t )) D. This implies that for all t t the trajectories in D are boune, an are conne insie the strip jz(t)j jz()j. However, accoring to the previous iscussion, the trajectory with initial conition (z(t ) v(t )) satises lim t!(z(t) v(t)) =. Therefore, it has been shown that for all (z() v()) IR (IR +nfg) the trajectories remain boune an have the property that lim t!(z(t) v(t)) =. By the enition of v, this implies that lim (z(t) w(t)) = () t! In orer to show that! is boune, write the ratio z= w = w. From Eq. (9) one obtains j!j cjwj + cjjjwj (4) From Corollary 6., for all initial conitions (z() w()) IR (Cnfg) is boune. Since w is also boune, from Eq. (4) it follows that! is boune. From Eq. (a) it follows immeiately that _w is also boune. Moreover, since _ = ;() + ()( + v) (5) an () () v an are all boune, we have that _ is boune. The erivative of! is given by _! = ; _()w ; () _w ; i _() w ; i()_w ; i() _w (6) Using Eqs. () one has _() = ; 4c _() = ; c _ (7a) + ( ; ) _ (7b) + ( ; ) Since _ is boune, _() an _() are both boune. Finally, the bouneness of _! follows irectly from Eq. (6) an the fact all the terms in the right han sie of this equation are boune. The vector el an the corresponing trajectories of the close-loop system with the control law in Eq. (9) is shown in Fig. 4 (compare with Fig. ). Remark 6. Theorem 6. shows that for all initial conitions w() 6= the control law in Eq. (9) rives the system trajectories to the origin. This control law cannot be use if w() = (an z 6= ). Linearization of the 5

6 Fig. 4 Phase portrait of system in Eqs. (4). system represente by Eq. () about w =, however, shows that this system is controllable an choosing, for example, a constant control! =! c C, one can move away from the z-axis into the D region once in D, use of the control in Eq. (9) rives the system to the origin. Remark 6. Another choice of a feeback control for Eq. () is the sublinear control in terms of w,! = ; w +jwj ; i z w which reners the close-loop system (8) _v = ;v (9a) _z = ;z (9b) globally exponentially stable. The previous methoology can be applie mutatis mutanis to this control law, as well. Moreover, several other similar moications can be introuce to the control law in Eq. (4). It shoul be evient that the results in this section can be applie to these control laws with only minor moications. 6. Propose Control Law for Complete System The control law in Eq. (9) was shown to achieve lim t!(z(t) w(t)) =. Moreover, it is a boune controller with boune erivative. This allows one to implement this control through the ynamics in Eq. (a). To this en, ene the error e =! ;! (4) where! is the esire angular velocity prole given in Eq. (9). Consier the following feeback control u =_! ; (! + ()w + i()w) (4) where > an where _! is given in Eq. (6), along with Eqs. (7). The value of _ is now given by _ = ;()+()(+v)+im(e=w);(+v)re(e=w) (4) With the control law in Eq. (4) the close-loop system takes the form _e = ;e (4a) _v = ;()( + v)v +( + v)re(e w) (4b) _z = ;()z + Im(e w) (4c) Notice that for e = the system reuces to the one in Eqs. (4). For large enough, Eq. (4a) is essentially a bounary layer subsystem to the slow system given by Eqs. (4b)-(4c). Singular perturbation theory guarantees that as soon as the error becomes small enough, the (z v) trajectories of the system will follow the ones of Eqs. (4). Next we show that the control law in Eq.(4)is wellene, in the sense that it remains boune for all t. We show that with large enough w(t) 6= for all t, i.e., w(t) tens to zero only asymptotically for all initial conitions insie an a priori given compact set. Proposition 6. Consier the system in Eqs. (4) an the compact set N = f(! w z) W: jej[( + v)=v] g (44) where W = C (Cnfg) IR, an let c > c > an >( c + )=. Then for all initial conitions in N, jwj is boune belowbyanexponentially ecaying function. Proof. Equation (4b) can be re-written as t jwj = ;( + jwj )(()jwj ; Re(e w)) (45) Note that from Eq. (4a) je(t)j je()je ;t an using Eq. (44), je(t)j jw()j +jw()j Consier now the ierential equation e ;(c +)t= t (46) t j ^wj = ;( c + )( + j ^wj )j ^wj (47) The solution of this equation is given by j ^w(t)j = c e (c+)t= ; c; e ;(c+)t= (48) where c = (j ^w()j +)=j ^w()j. Comparison of Eqs. (46) an (48) implies that je(t)j j ^w(t)j 8 t (49) where j ^wj obeys Eq. (47) with j ^w()j = jw()j. 6

7 Notice now that since Re(e w) jejjwj, an using Eq. (49), one has from Eq. (45) that t jwj ;( + jwj )(()jwj + jejjwj) ;( + jwj )(()jwj + j ^wjjwj) (5) an since ; c () c nally, t jwj ;( + jwj )( cjwj + j ^wjjwj) (5) By comparing Eqs. (47) an (5) an since jw()j = j ^w()j, one obtains t jw()j t j ^w()j (5) Therefore there exist some t > such that jw(t)j j ^w(t)j for all t t. We claim that, actually, jw(t)j j^w(t)j for all t, an thus jwj is boune below by the exponentially ecaying function j ^wj. Assume that at some point t > we have that jw(t )j = j ^w(t )j an jw(t )j < j ^w(t )j see Fig. 5. t t Then the anger that w will move towars the z-axis before the control law in Eq. (4) becomes eective. This is one more reason which motivate the choice of the control law in Eq. (9). Namely, it is benecial for w to move away from the z-axis. This can reuce the value of the gain signicantly. In most situations it is not necessary to chose from Proposition 6.. Actually, asthenumerical simulations in the next section show, for most practical examples it suces to choose to be \suciently larger" than the gains c an c. From Eq. (4) it is also clear that shoul be at least as large as c=, in orer for e=w to remain boune. Remark 6. The rigi boy problem subject to two control inputs is only but one example of an uneractuate mechanical system. Systems of this form can be foun in the class of systems subject to nonholonomic, i.e., non-integrable constraints. Time-invariant control laws for these systems are necessarily nonsmooth an recently propose control laws {5 inclue singularities of the same form as in Eq. (4). It is therefore conceptually straightforwar to exten the results of this paper to this more general case. w ^ w 7 Numerical Example To illustrate the previous theoretical analysis, we have simulate the ierential equations () with the two control laws in Eqs. (4) an (9). The gains are chosen as c =:5 an c =. The value of the parameter =. The initial conitions were taken as w() = : ; i :5 an z() = :5. The results are shown in Figs. 6 an 7. Figure 6 shows the correspon- t Fig. 5 Time history of jwj an j ^wj. t.5.5 Ol Metho New Metho t jw(t )j = ;( + jw(t )j )( cjw(t )j + jw(t )j ) = ;( + j ^w(t )j )( c + )j ^w(t )j = t j ^w(t )j (5) which leas to a contraiction. Therefore jw(t)j j ^w(t)j an thus w(t) 6= for all t. In Ref. the control law in Eq. (4) was also implemente using the same methoology. That is, the control for the complete system was given by Eq. (4), where = c, = c an _ =_ =. The value of the gain increases with, which in turns increases as jwj ecreases. That is, when the initial conition is close to w = then a faster transient for! is require. This faster transient isachieve by taking large enough. A potential problem in the implementation of the control in Eq. (4) is now evient. If e oes not ecay \fast enough" so that!!! suciently fast, then there is Fig. 6 Close-loop trajectories for the two methos (kinematics only). ing close-loop trajectories, an Fig. 7 shows the magnitue of the angular velocity (control input for the kinematics) j!j. Soli lines correspon to the new control law in Eq. (9) an the ashe lines correspon to the previous control law given in Eq. (4). As it is evient 7

8 4 Ol Metho New Metho.5 α = α = 4 α = 8 ω u Time Fig. 7 Control eort for the two methos (kinematics only). from these gures, there is a substantial ecrease in control eort by using the control law in Eq. (9), especially uring the initial portion of the trajectory where z is large an jwj is small α = This control law was later implemente through the ynamics in Eq. (a). A rest-to-rest maneuver was consiere, thus!() =. Simulations for several values of are shown in Figs The trajectories in the (z v) space are very similar to the ones when! is the control input. In fact, for = the trajectories for the comα = 4 α =.5.5 Fig. 8 Close-loop trajectories for the complete system. plete system are essentially ientical to the ones with control law in Eq. (9). Figure 9 shows that increasing may increase the control eort, mainly because of the high-gain bounary layer part of the controller. At any rate, the corresponing control eort for the control law in Ref. is several orers of magnitue higher (not shown here). Moreover, for small such as = an = 4, the control eort for the controller in Ref. is not boune. In these cases the slow transients of Time Fig. 9 Control eort for the complete system. e allowe w to rift towars the z-axis before the control law in Eq. (4) is activate. On the other han, the controller in Eq. (9) forces the system trajectories away fromthez-axis, thus proviing enough time for the ynamic controller to \catch up." 8 Conclusions A nonsmooth control law has been constructe which stabilizes the kinematics of an uneractuate rigi spacecraft. It is shown that the propose control law is well ene an uses consierably less control effort than a previously erive control law. The main iea is to ivie the state space into two regions, one which inclues initial conitions resulting in high control expeniture an one which inclues initial conitions resulting in acceptable control input signals. The propose control law then forces all the close-loop system trajectories to leave the unesirable region of high control eort an subsequently use the original control law. Numerical examples inicate a signicant control eort reuction using the new control scheme. Because of the limite control torque on-boar a spacecraft, for practical situations this may be the ierence between feasibility an infeasibility of a particular reorientation maneuver. Acknowlegment: Financial support from the National Science Founation uner Grant CMS is gratefully acknowlege. References Aeyels, D., \Stabilization by SmoothFeeback ofthe Angular elocity of a Rigi Boy," Systems an Control Letters, ol. 6, No., 985, pp. 59{6. Aeyels, D. an Szafranski, M., \Comments on the Stabilizability of the Angular elocity of a Rigi Boy," Systems an Control Letters, ol., No., 988, 8

9 pp. 5{9. Outbib, R. an Sallet, G., \Stabilizability oftheangular elocity of a Rigi Boy Revisite," Systems an Control Letters, ol. 8, No., 99, pp. 9{98. 4 Sontag, E. an Sussmann, H., \Further Comments of the Stabilizability of the Angular elocity of a Rigi Boy," Systems an Control Letters, ol., No., 988, pp. {7. 5 Leonar, N. E. an Krishnaprasa, P. S., \Averaging for Attitue Control an Motion Planning," in Proceeings of the n Conference on Decision an Control, pp. 98{4, 99. San Antonio, TX. 6 Leonar, N. E. an Krishnaprasa, P. S., \High- Orer Averaging on Lie Groups an Control of an Autonomous Unerwater ehicle," in Proceeings of the American Control Conference, pp. 57{6, 994. Baltimore, MD. 7 Byrnes, C. I. an Isiori, A., \On the attitue stabilization of a rigi spacecraft," Automatica, ol. 7, No., 99, pp. 87{95. 8 Krishnan, H., McClamroch, H., an Reyhanoglu, M., \On the attitue stabilization of a rigi spacecraft using two control torques," in Proceeings of the American Control Conference, pp. 99{995, 99. Chicago, IL. 9 Krishnan, H., Reyhanoglu, M., an McClamroch, H., \Attitue stabilization of a rigi spacecraft using two control torques: A nonlinear control approach base on the spacecraft attitue ynamics," Automatica, ol., No. 6, 994, pp. {7. Tsiotras, P., Corless, M., an Longuski, M., \A Novel Approach for the Attitue Control of an Axisymmetric Spacecraft Subject to Two Control Torques," Automatica, ol., No. 8, 995, pp. 99{. Walsh, G. C., Montgomery, R., an Sastry, S. S., \Orientation control of a ynamic satellite," in Proceeings of the American Control Conference, pp. 8{ 4, 994. Baltimore, MD. Morin, P., Samson, C., Pomet, J., an Jiang,., \Time-arying Feeback Stabilization of the Attitue of a Rigi Spacecraft with Two Controls," in Proc. r IEEE Conf. on Decision an Control, pp. 94{ 95, Dec. 4-6, 994. Lake Buena ista, FL. Coron, J. M. an Kerai, E. L., \Explicit Feeback Stabilizing the Attitue of a Rigi Spacecraft with Two Control Torques," Automatica, ol. 6, No. 5, 996, pp. 669{ Tsiotras, P. an Longuski, J. M., \New kinematic relations for the large angle problem in rigi boy attitue ynamics," Acta Astronautica, ol., No., 994, pp. 8{9. 5 Tsiotras, P. an Longuski, J. M., \A New Parameterization of the Attitue Kinematics," Journal of the Astronautical Sciences, ol. 4, No., 995, pp. 4{ 6. 6 Morin, P. an Samson, C., \Time-arying Exponential Stabilization of the Attitue of a Rigi Spacecraft with Two Controls," in Proceeings of the 4th Conference on Decision an Control, pp. 988{99, 995. New Orleans, LA. 7 Shuster, M. D., \A Survey of Attitue Representations," Journal of the Astronautical Sciences, ol. 4, No. 4, 99, pp. 49{57. 8 Longuski, J. M. an Tsiotras, P., \Analytic solution of the large angle problem in rigi boy attitue ynamics," Journal of the Astronautical Sciences, ol. 4, No., 995, pp. 5{46. 9 Tsiotras, P. an Longuski, J. M., \Spin-Axis Stabilization of Symmetric Spacecraft with Two Control Torques," Systems an Control Letters, ol., No. 6, 994, pp. 95{4. Brockett, R. W., \Asymptotic Stability an Feeback Stabilization," in Dierential Geometric Control Theory (Brockett, R. W., Millman, R. S., an Sussman, H. J., Es.), Birkhauser, 98, pp. 8{8. Khalil, H. K., Nonlinear Systems. New Jersey, Prentice Hall, n e., 996, p. 8. Murray, R. M., Li,., an Sastry, S.S.,A Mathematical Introuction to Robotic Manipulation. Boca Raton, Floria, CRC Press, 994. Astol, A., \Discontinuous Control of Nonholonomic Systems," Systems an Control Letters, ol. 7, 996, pp. 7{45. 4 Khennouf, H. an Canuas e Wit, C., \On the Construction of Stabilizing Discontinuous Controllers for Nonholonomic Systems," in Proc. IFAC Nonlinear Control Systems Design Symposium, pp. 747{75, 995. Tahoe City, CA. 5 Luo, J. an Tsiotras, P., \Exponentially Convergent Control Laws for n-dimensional Nonholonomic Systems in Power Form," in Proceeings of the American Control Conference, 997. Albuquerque, NM. 9

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

Abstract: It is a well known fact that a symmetric spacecraft with two control torques supplied by gas

Abstract: It is a well known fact that a symmetric spacecraft with two control torques supplied by gas Systems & Control Letters 23 (1994) 395-402 North Holland 395 Spin-axis stabilization of symmetric spacecraft with two control torques Panagiotis Tsiotras and James M. Longuski School of Aeronautics and

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

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department

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

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

Accelerate Implementation of Forwaring Control Laws using Composition Methos Yves Moreau an Roolphe Sepulchre June 1997 Abstract We use a metho of int

Accelerate Implementation of Forwaring Control Laws using Composition Methos Yves Moreau an Roolphe Sepulchre June 1997 Abstract We use a metho of int Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 1997-11 Accelerate Implementation of Forwaring Control Laws using Composition Methos 1 Yves Moreau, Roolphe Sepulchre, Joos Vanewalle

More information

initial configuration initial configuration end

initial configuration initial configuration end Design of trajectory stabilizing feeback for riftless at systems M. FLIESS y J. L EVINE z P. MARTIN x P. ROUCHON { ECC95 Abstract A esign metho for robust stabilization of riftless at systems aroun trajectories

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

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

Problems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract

Problems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract Pseuo-Time Methos for Constraine Optimization Problems Governe by PDE Shlomo Ta'asan Carnegie Mellon University an Institute for Computer Applications in Science an Engineering Abstract In this paper we

More information

1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a

1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a Theory of the nerson impurity moel: The Schrieer{Wol transformation re{examine Stefan K. Kehrein 1 an nreas Mielke 2 Institut fur Theoretische Physik, uprecht{karls{universitat, D{69120 Heielberg, F..

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

PROPORTIONAL DERIVATIVE (PD) CONTROL ON THE EUCLIDEAN GROUP F. BULLO AND R. M. MURRAY. Division of Engineering and Applied Science

PROPORTIONAL DERIVATIVE (PD) CONTROL ON THE EUCLIDEAN GROUP F. BULLO AND R. M. MURRAY. Division of Engineering and Applied Science PROPORTIONAL DERIVATIVE (PD) CONTROL ON THE EUCLIDEAN GROUP F. BULLO AND R. M. MURRAY Division of Engineering an Applie Science California Institute of Technology Pasaena, CA 95 CDS Technical Report 95-

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

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

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

has been extensively investigated during the last few years with great success [7,4,5,9,4,7]. It can be shown that time-varying smooth control laws fo

has been extensively investigated during the last few years with great success [7,4,5,9,4,7]. It can be shown that time-varying smooth control laws fo Control Design for Chained-Form Systems with Bounded Inputs? Jihao Luo Department of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, VA 93-44, USA Panagiotis Tsiotras School

More information

*School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN

*School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN NEW CONTROL LAWS FOR THE ATTITUDE STABILIZATION OF RIGID BODIES PANAGIOTIS TSIOTRAS *School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN 7907. Abstract. This paper introduces

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

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

ATTITUDE STABILIZATION OF A SATELLITE BY MAGNETIC COILS

ATTITUDE STABILIZATION OF A SATELLITE BY MAGNETIC COILS Acta Astronautica Vol. 50, No. 12, pp. 721 728, 2002? 2002 Publishe by Elsevier Science Lt. Printe in Great Britain www.elsevier.com/locate/actaastro PII: S0094-5765(02)00011-5 0094-5765/02/$ - see front

More information

TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM

TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM 265 Asian Journal of Control, Vol. 4, No. 3, pp. 265-273, September 22 TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM Jurachart Jongusuk an Tsutomu Mita

More information

All s Well That Ends Well: Supplementary Proofs

All s Well That Ends Well: Supplementary Proofs All s Well That Ens Well: Guarantee Resolution of Simultaneous Rigi Boy Impact 1:1 All s Well That Ens Well: Supplementary Proofs This ocument complements the paper All s Well That Ens Well: Guarantee

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

'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

Stabilization of Angular Velocity of Asymmetrical Rigid Body. Using Two Constant Torques

Stabilization of Angular Velocity of Asymmetrical Rigid Body. Using Two Constant Torques Stabilization of Angular Velocity of Asymmetrical Rigid Body Using Two Constant Torques Hirohisa Kojima Associate Professor Department of Aerospace Engineering Tokyo Metropolitan University 6-6, Asahigaoka,

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

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

Control of underactuated spacecraft with bounded inputs

Control of underactuated spacecraft with bounded inputs Automatica 36 (2000) 1153}1169 Control of underactuated spacecraft with bounded inputs Panagiotis Tsiotras*, Jihao Luo School of Aerospace Engineering, Georgia Institute of Technology, Atlanta, GA 30332-0150,

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

Approximate Reduction of Dynamical Systems

Approximate Reduction of Dynamical Systems Proceeings of the 4th IEEE Conference on Decision & Control Manchester Gran Hyatt Hotel San Diego, CA, USA, December 3-, 6 FrIP.7 Approximate Reuction of Dynamical Systems Paulo Tabuaa, Aaron D. Ames,

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

Interpolated Rigid-Body Motions and Robotics

Interpolated Rigid-Body Motions and Robotics Interpolate Rigi-Boy Motions an Robotics J.M. Selig Faculty of Business, Computing an Info. Management. Lonon South Bank University, Lonon SE AA, U.K. seligjm@lsbu.ac.uk Yaunquing Wu Dept. Mechanical Engineering.

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

Implicit Lyapunov control of closed quantum systems

Implicit Lyapunov control of closed quantum systems Joint 48th IEEE Conference on Decision an Control an 28th Chinese Control Conference Shanghai, P.R. China, December 16-18, 29 ThAIn1.4 Implicit Lyapunov control of close quantum systems Shouwei Zhao, Hai

More information

Short Intro to Coordinate Transformation

Short Intro to Coordinate Transformation Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent

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

Experimental Robustness Study of a Second-Order Sliding Mode Controller

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

More information

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

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

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

State observers and recursive filters in classical feedback control theory

State observers and recursive filters in classical feedback control theory State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Adaptive Optimal Path Following for High Wind Flights

Adaptive Optimal Path Following for High Wind Flights Milano (Italy) August - September, 11 Aaptive Optimal Path Following for High Win Flights Ashwini Ratnoo P.B. Sujit Mangal Kothari Postoctoral Fellow, Department of Aerospace Engineering, Technion-Israel

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Preprints of the 8th IFAC Worl Congress Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Guoong Shi ACCESS Linnaeus Centre, School of Electrical Engineering, Royal

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

FORMATION INPUT-TO-STATE STABILITY. Herbert G. Tanner and George J. Pappas

FORMATION INPUT-TO-STATE STABILITY. Herbert G. Tanner and George J. Pappas Copyright 2002 IFAC 5th Triennial Worl Congress, Barcelona, Spain FORMATION INPUT-TO-STATE STABILITY Herbert G. Tanner an George J. Pappas Department of Electrical Engineering University of Pennsylvania

More information

the solution of ()-(), an ecient numerical treatment requires variable steps. An alternative approach is to apply a time transformation of the form t

the solution of ()-(), an ecient numerical treatment requires variable steps. An alternative approach is to apply a time transformation of the form t Asymptotic Error Analysis of the Aaptive Verlet Metho Stephane Cirilli, Ernst Hairer Beneict Leimkuhler y May 3, 999 Abstract The Aaptive Verlet metho [7] an variants [6] are time-reversible schemes for

More information

VISUAL SERVOING WITH ORIENTATION LIMITS OF A X4-FLYER

VISUAL SERVOING WITH ORIENTATION LIMITS OF A X4-FLYER VISUAL SERVOING WITH ORIENTATION LIMITS OF A X4-FLYER Najib Metni,Tarek Hamel,Isabelle Fantoni Laboratoire Central es Ponts et Chaussées, LCPC-Paris France, najib.metni@lcpc.fr Cemif-Sc FRE-CNRS 2494,

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

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Advanced Partial Differential Equations with Applications

Advanced Partial Differential Equations with Applications MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.

More information

2.5. x x 4. x x 2. x time(s) time (s)

2.5. x x 4. x x 2. x time(s) time (s) Global regulation and local robust stabilization of chained systems E Valtolina* and A Astolfi* Π *Dipartimento di Elettronica e Informazione Politecnico di Milano Piazza Leonardo da Vinci 3 33 Milano,

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210 IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc

More information

A Path Planning Method Using Cubic Spiral with Curvature Constraint

A Path Planning Method Using Cubic Spiral with Curvature Constraint A Path Planning Metho Using Cubic Spiral with Curvature Constraint Tzu-Chen Liang an Jing-Sin Liu Institute of Information Science 0, Acaemia Sinica, Nankang, Taipei 5, Taiwan, R.O.C., Email: hartree@iis.sinica.eu.tw

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

Some Remarks on the Boundedness and Convergence Properties of Smooth Sliding Mode Controllers

Some Remarks on the Boundedness and Convergence Properties of Smooth Sliding Mode Controllers International Journal of Automation an Computing 6(2, May 2009, 154-158 DOI: 10.1007/s11633-009-0154-z Some Remarks on the Bouneness an Convergence Properties of Smooth Sliing Moe Controllers Wallace Moreira

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

6.003 Homework #7 Solutions

6.003 Homework #7 Solutions 6.003 Homework #7 Solutions Problems. Secon-orer systems The impulse response of a secon-orer CT system has the form h(t) = e σt cos(ω t + φ)u(t) where the parameters σ, ω, an φ are relate to the parameters

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

Stabilization of a Class of Underactuated Mechanical Systems via Interconnection and Damping Assignment

Stabilization of a Class of Underactuated Mechanical Systems via Interconnection and Damping Assignment Stabilization of a Class of Uneractuate Mechanical Systems via Interconnection an Damping Assignment Romeo Ortega Lab. es Signaux et Systèmes CNRS-SUPELEC Gif sur Yvette 99 FRANCE tel. no:(33)--69-85-7-65

More information

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,

More information

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS

More information

Distributed coordination control for multi-robot networks using Lyapunov-like barrier functions

Distributed coordination control for multi-robot networks using Lyapunov-like barrier functions IEEE TRANSACTIONS ON 1 Distribute coorination control for multi-robot networks using Lyapunov-like barrier functions Dimitra Panagou, Dušan M. Stipanović an Petros G. Voulgaris Abstract This paper aresses

More information

Connections Between Duality in Control Theory and

Connections Between Duality in Control Theory and Connections Between Duality in Control heory an Convex Optimization V. Balakrishnan 1 an L. Vanenberghe 2 Abstract Several important problems in control theory can be reformulate as convex optimization

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

H. MOUNIER, J. RUDOLPH, M. FLIESS, AND P. ROUCHON. interior mass is considered. It is viewed as a linear delay system. A

H. MOUNIER, J. RUDOLPH, M. FLIESS, AND P. ROUCHON. interior mass is considered. It is viewed as a linear delay system. A ESAIM: Control, Optimisation an Calculus of Variations URL: http://www.emath.fr/cocv/ September 998, Vol. 3, 35{3 TRACKING CONTROL OF A VIBRATING STRING WITH AN INTERIOR MASS VIEWED AS DELAY SYSTEM H.

More information

Dissipative numerical methods for the Hunter-Saxton equation

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

More information

arxiv: v1 [math.ds] 21 Sep 2017

arxiv: v1 [math.ds] 21 Sep 2017 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION WITH POSITIVE FEEDBACK arxiv:709.07295v [math.ds] 2 Sep 207 ISTVÁN GYŐRI, YUKIHIKO NAKATA, AND GERGELY RÖST Abstract. We stuy boune, unboune

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

Local Input-to-State Stabilization of 1-D Linear Reaction-Diffusion Equation with Bounded Feedback

Local Input-to-State Stabilization of 1-D Linear Reaction-Diffusion Equation with Bounded Feedback Local Input-to-State Stabilization of -D Linear Reaction-Diffusion Equation with Boune Feeback Aneel Tanwani, Swann Marx, Christophe Prieur To cite this version: Aneel Tanwani, Swann Marx, Christophe Prieur.

More information

Left-invariant extended Kalman filter and attitude estimation

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

More information

Stable Adaptive Control and Recursive Identication. Using Radial Gaussian Networks. Massachusetts Institute of Technology. functions employed.

Stable Adaptive Control and Recursive Identication. Using Radial Gaussian Networks. Massachusetts Institute of Technology. functions employed. NSL-910901, Sept. 1991 To appear: IEEE CDC, Dec. 1991 Stable Aaptive Control an Recursive Ientication Using Raial Gaussian Networks Robert M. Sanner an Jean-Jacques E. Slotine Nonlinear Systems Laboratory

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

The effect of dissipation on solutions of the complex KdV equation

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

More information

The Sokhotski-Plemelj Formula

The Sokhotski-Plemelj Formula hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),

More information

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into

More information

Position in the xy plane y position x position

Position in the xy plane y position x position Robust Control of an Underactuated Surface Vessel with Thruster Dynamics K. Y. Pettersen and O. Egeland Department of Engineering Cybernetics Norwegian Uniersity of Science and Technology N- Trondheim,

More information

Linear Feedback Control Using Quasi Velocities

Linear Feedback Control Using Quasi Velocities Linear Feedback Control Using Quasi Velocities Andrew J Sinclair Auburn University, Auburn, Alabama 36849 John E Hurtado and John L Junkins Texas A&M University, College Station, Texas 77843 A novel approach

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

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement

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

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

RECENTLY, vertical take-off and landing (VTOL) unmanned

RECENTLY, vertical take-off and landing (VTOL) unmanned IEEE/CAA JOURNAL OF AUTOMATICA SINICA VOL. NO. JANUARY 5 65 Trajectory Tracking of Vertical Take-off an Laning Unmanne Aerial Vehicles Base on Disturbance Rejection Control Lu Wang an Jianbo Su Abstract

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

The canonical controllers and regular interconnection

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

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

A new approach to explicit MPC using self-optimizing control

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

More information

Neuro-Fuzzy Processor

Neuro-Fuzzy Processor An Introuction to Fuzzy State Automata L.M. Reyneri Dipartimento i Elettronica - Politecnico i Torino C.so Duca Abruzzi, 24-10129 Torino - ITALY e.mail reyneri@polito.it; phone ++39 11 568 4038; fax ++39

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

arxiv: v1 [math.oc] 25 Nov 2017

arxiv: v1 [math.oc] 25 Nov 2017 Constraine Geometric Attitue Control on SO(3) Shankar Kulumani* an Taeyoung Lee November 8, 017 arxiv:17199v1 [math.oc] 5 Nov 017 Abstract This paper presents a new geometric aaptive control system with

More information

Delocalization of boundary states in disordered topological insulators

Delocalization of boundary states in disordered topological insulators Journal of Physics A: Mathematical an Theoretical J. Phys. A: Math. Theor. 48 (05) FT0 (pp) oi:0.088/75-83/48//ft0 Fast Track Communication Delocalization of bounary states in isorere topological insulators

More information

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract

More information

Proof of SPNs as Mixture of Trees

Proof of SPNs as Mixture of Trees A Proof of SPNs as Mixture of Trees Theorem 1. If T is an inuce SPN from a complete an ecomposable SPN S, then T is a tree that is complete an ecomposable. Proof. Argue by contraiction that T is not a

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