Observers for systems with invariant outputs

Size: px
Start display at page:

Download "Observers for systems with invariant outputs"

Transcription

1 Observers for systems with invariant outputs C. Lageman, J. Trumpf an R. Mahony Abstract In this paper we introuce a general esign approach for observers for left-invariant systems on a Lie group with right-invariant outputs on a manifol. This situation is encountere for example in attitue estimation problems where only partial information of the state can be measure. We propose a graient-like observer esign on the output space which can be lifte to a graient-like observer on the state space. The observer on the output space converges for generic initial conitions an has local exponential convergence properties. The lifte observer retains the same convergence properties with respect to the fibers of the output map. Our observer constructions are illustrate on the attitue estimation problem for a rigi boy. I. INTRODUCTION Nonlinear observer esign for invariant systems on Lie groups has receive a growing amount of interest in the last years, riven mainly by applications in state estimation for autonomouos robotic vehicles. The interest in nonlinear observers is motivate by their potential lower computational complexity an their increase robustness compare to stanar filtering approaches. Invariant systems on Lie groups have been stuie in a series of works starting in the seventies, focusing mainly on general, theoretical observability an controllablity questions [1], [2], [3], [4], [5]. Nonlinear observer esign approaches for invariant systems have been generally more closely connecte to specific applications. A large boy of works consiers nonlinear attitue observers for rigi boy ynamics [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17] or the full pose estimation problem [18], [19], [20], [21]. However, theoretical aspects of observer esign for invariant systems have receive little attention in the past. Only recent works have begun to tackle such questions an consier the general structure of observers for invariant systems an propose general approaches to observer esign [22], [23], [14], [24], [25], [26], [27]. In this paper, we exten the graient-like observer construction for full state observations propose in [25], [27] to systems with partial state observations given by a rightinvariant output map. This system set up has been consiere in [26], where the authors iscuss a general observer structure which has autonomous error ynamics an yiels a local convergence result for a suitable choice of gain coefficients. C. Lageman is with the Department of Electrical Engineering an Computer Science, Université e Liège, B-4000 Liege Sart-Tilman, Belgium christian.lageman@montefiore.ulg.ac.be J. Trumpf is with the Department of Information Engineering, The Australian National University, Canberra ACT 0200, Australia Jochen.Trumpf@anu.eu.au R. Mahony is with the Department of Engineering, The Australian National University, Canberra ACT 0200, Australia Robert.Mahony@anu.eu.au In this paper, we approach the observer construction problem by first esigning a graient-like observer on the output space. This is justifie by the fact that the system projects to a system with full state observations on the output space. For this projecte system a moification of the graient-like error construction for full state observations on a Lie group [25], [27] can be applie, which yiels an output observer for the system on the group. To obtain an observer for the group state space we use a lift of this output observer. This lift can be interprete as a graientlike observer on the Lie group. This links our approach both to [26] an the graient-like observers for full state observation [25], [27] as the lifte observer fits the observer structures iscusse in these papers. We show that uner mil conitions the errors for the output observer an the group state observer converge to a reference output or the stabilizer of this reference output, respectively. Furthermore, this convergence is locally exponential. To illustrate our concepts, we iscuss an attitue estimation problem for a rigi boy with partial state measurements. In Section II we introuce our notation an the system set up. Section III iscusses the observer construction on the output space an the convergence properties. The lift to the group an the relate convergence results are consiere in Section IV. Section V illustrates the methos on the attitue estimation problem. The paper ens with conclusions in Section VI. II. INVARIANT SYSTEMS WITH OBSERVATIONS ON A A. Notation HOMOGENEOUS SPACE We start with the notation use in the paper. We enote by G a connecte, finite-imensional Lie group with Lie algebra g an ientity e. M will enote a manifol on which G acts by a right action h: G M M, i.e. h(e, y) = y for all y M, e G the ientity an h(x, h( ˆX, y)) = h( ˆXX, y) for all X, ˆX G an y M. The action h is transitive if for all y, ŷ M there exists an X G with h(y, X) = ŷ. The group acts effectively on M if for all X G, X e exists a y M (possibly epening on X) with h X (y) y. To simplify our notation we use h y for the map h y : G M, h y (X) = h(x, y) an h X for the map h X : M M, hx (y) = h(x, y). We enote by stab(y) for y M the stabilizer group stab(y) = {X G h(x, y) = y}. We use the notation T X G an T y M for the tangent spaces of G an M at X an y, respectively. The tangent map for a function φ: G M at X G is enote by T X φ. A Riemannian metric, on G is calle right-invariant if for any X, Y G, v, w g an the

2 right-multiplication map R Y : G G, R Y (X) = XY, one has T e R Y v, T e R Y w G = v, w G. A Riemannian metric on M is calle invariant (uner h) if for all X G, y M, v, w T y M one has v, w = T y hx v, T y hx. A istribution H(X) on G is calle right-invariant if for all X, Y G the equation H(XY ) = T X R Y H(X) hols. We say that a right-invariant Riemannian metric on G inuces an invariant Riemannian metric on M if there is a right-invariant horizontal istribution H(X) such that for all y M, X G, v, w T X G we have T X h y (Xv), T X h y (Xw) = v, w G. We will use the notation (v) H for the horizontal lift of v T y M, y M to T X G, X stab(y) with respect to a horizontal istribution H(X), i.e. the unique w H(X) with T X h y w = v. At last, we recall that the graient of a function f on G or M with respect to a Riemannian metric, is given by the ientity gra f(x), v = f(x)(v) for all v T x G (or all v T x M, respectively). B. System set-up an assumptions We introuce now the systems iscusse in this paper. It is assume that we are given a connecte Lie group G an a smooth manifol M on which G acts by a transitive right action. We consier a left-invariant system on G with rightinvariant outputs on M, i.e. a system of the form Ẋ = Xu, X G y = h(x, ), y M with u: R g a smooth input an a fixe reference output. We will refer to M as the output space. Since our observer construction will be graient-base, we assume that there is a right-invariant metric on G which inuces an invariant metric on M. This assumption is satisfie e.g. for all reuctive homogeneous spaces M (cf. [28] for a efinition) on which G acts effectively, see [29, Proposition 3.16] 1. This inclues in particular G = SO(n) with M = S n 1 an the action h(r, y) = R y for R SO(n) an y S n 1, cf. Section V. The right-invariant horizontal istribution associate with the pair of Riemannian metrics is enote by H(X). In this paper we consier the problem of constructing an observer for systems of the form (1) given only measurements of the input u an the output y. III. GRADIENT-LIKE OBSERVERS ON THE HOMOGENEOUS SPACE Our observer esign approach is base on the fact that the invariant system (1) can be reuce to a system on the output space with full state observations. We will propose an observer for this projecte system similar to the graient-like observers for invariant systems with full state observations on the group propose in [25], [27]. Let us begin with the reuction of the system to the output space. 1 Usually only left actions on M are consiere in the literature but the results exten analogously to right actions by interchanging left- an rightinvariance. (1) Theorem 1: System (1) projects to the system ẏ = T X h y0 (Xu) (2) with X h 1 (y). This system is well efine an a minimal realization of (1). This result is mainly base on the symmetry reuction methos as introuce in [30]. We omit a etaile proof of the theorem. As in the Lie group case with full state observations [25], [27], the observer on the homogeneous space consists of a copy of the system ynamics an a graient-like innovation which will rive a suitably chosen error term into a no error state. For that, we assume that we are given a smooth cost function f : M M R on M. The propose observer has the form: ŷ = T ˆXh y0 ( ˆXu) gra 1 f(ŷ, y). (3) To iscuss the convergence properties of the observer for the projecte system, i.e. the output, we have to introuce a suitable error function which compares the observer state with the state of the original system on the group. Definition 2: We efine the output error E as the function E : M G M, E(ŷ, X) := h( ˆXX 1, ) = h(x 1, ŷ) ˆX h 1 (ŷ). Note that this error is well-efine an a smooth function. The error value correspons to the no error states, i.e. that the state of the observer is the output corresponing to the state X. The use of this error is justifie by the fact that any other smooth error function Ê : M M N which is constant on two copies of the system (2) can be lifte to M G as Ê(ŷ, h (X)) an this lift is the composition g E of a smooth, stab( )-invariant function g : M N with the output error. Uner suitable invariance conitions on the cost, the output error has autonomous graient ynamics. Theorem 3: Assume that f is invariant uner the action of G on M M, i.e. f(h(s, ŷ), h(s, y)) = f(ŷ, y) an that the conitions on the Riemannian metrics of Section II hol. Then we have the following autonomous ynamics for the output error. t E(ŷ, X) = gra 1 f(e(ŷ, X), ). Proof: In this proof we will enote the left an right multiplication in G by Y G as L Y an R Y. Straightforwar calculations using the equivariance of h show that h y L X = h h(x,y) an h X h y = h y R X. Thus we have E(ŷ, X) t = T X 1hŷT e L X 1T X R X 1Xu +Tŷ hx 1T ˆXh y0 ˆXu Tŷ hx 1 gra 1 f(ŷ, y) = T X 1hŷT e R X 1u + T ˆXX 1h y0 T ˆXR X 1 ˆXu Tŷ hx 1 gra 1 f(ŷ, y) = T ˆXX 1h y0 T X 1L ˆXT e R X 1u +T ˆXX 1h y0 T ˆXR X 1 ˆXu Tŷ hx 1 gra 1 f(ŷ, y) = Tŷ hx 1 gra 1 f(ŷ, y).

3 Since the Riemannian metric on M is invariant uner the action of h an f is invariant uner the action of G on M M, we get Tŷ hx 1 gra 1 f(ŷ, y) = gra 1 f( h X 1(ŷ), h X 1(y)) = gra 1 f(e(ŷ, X), ). This yiels the formula for the output error ynamics. The graient error ynamics yiels the following convergence result for the output error on M. Theorem 4: Assume that f is invariant uner the action of G an that the conitions on the Riemannian metrics of Section II hol. Furthermore let y f(y, ) be a Morse-Bott function with one global minimum an no further local minima. Then the output error on M converges to for generic initial conitions an the convergence is locally exponential. Proof: By Theorem 3 the output error has graient ynamics of the function y f(y, ). As the latter function is Morse-Bott, the graient system converges for generic initial conitions to the set of local minima an the convergence at a manifol of local minima is exponential [31]. Since is the only local minimum, we get generic convergence to at a locally exponential rate. Remark 5: Note that for all results in this section, it is sufficient that we are given an invariant Riemannian metric on M. The Riemannian metric on G will only play a role for the observer construction on G. IV. GRADIENT-LIKE OBSERVERS ON THE GROUP We turn now to the construction of an observer on the group. Since the projecte system on the output space is a minimal realization, a natural approach is to lift the observer equation for the projecte system to the Lie group. Here, we use a irect copy of the system ynamics Xu an lift the innovation term on M to the group by a horizontal lift with respect to the horizontal istribution H. This yiels the observer ˆX = ˆXu ( gra 1 f(h y0 ( ˆX), H y)), (4) where the term (v) H T ˆXG enotes the horizontal lift of v T h( ˆX,y0) M, ˆX G. Straightforwar calculations an the invariance properties of H(X) yiel the following projection result for the lifte observer. Proposition 6: Assume that the conitions on the Riemannian metrics of Section II hol. The observer (4) on G projects via the map h y0 to the observer (3) on M, i.e. ( ( T ˆXh y0 ˆXu gra 1 f(h y0 ( ˆX), ) ) H y) = T ˆXh y0 ( ˆXu) gra 1 f(ŷ, y). Furthermore, for any fixe initial values X(0), ˆX(0) G an any smooth input u the error E = ˆXX 1 of the observer on G projects via the map h y0 to the output error of the observer on M for initial values X(0), ŷ(0) = h y0 ( ˆX(0)) an the same input u. Proof: It follows from the efinition of the horizontal lift that (4) projects onto (3). For the statement on the observer error note that E(ŷ(t), X(t)) = h y0 ( ˆX(t)X(t) 1 ) by efinition of the output error, i.e. ˆXX 1 projects onto the output error. The observer (4) can be alternatively obtaine by the esign metho introuce in [25] for graient-like observers with full state observations. Our conitions on the Riemannian metrics on M an G imply that the horizontal lift of the graient of the cost function is the graient of a lifte cost function. Proposition 7: Assume that the conitions on the Riemannian metrics of Section II hol. Let F : G G R be the lift of the cost function f : M M R, i.e. F ( ˆX, X) = f(h( ˆX, ), h(x, )). The graient gra 1 F of the lifte cost function is the horizontal lift of gra 1 f. This property of the lift allows the interpretation of the lifte innovation term as a graient on G. Corollary 8: Uner the conitions of Proposition 7, the observer (4) can be written as ˆX = ˆXu gra 1 F ( ˆX, X). (5) In particular, the observer (5) projects to the observer (3) on the output manifol. Note that (5) is implementable using only observations of y an u since one has F ( ˆX, X) = f(h( ˆX, ), y) an thus the graient term oes only epen on y an ˆX. Since the observer an the error ˆXX 1 project to the observer an the output error on M, it inherits the convergence properties of the observer on M relative to the fibers of the map X h y0 (X). Here, convergence to a fiber h 1 (y) means that inf{ist( ˆX(t)X 1 (t), Z) Z h 1 (y)} converges to 0 with ist the Riemannian istance on G. Theorem 9: Assume that the conitions on f : M M R as in Theorem 4 hol. Then the error ˆXX 1 on G of the observer (4) converges to stab( ) for generic initial conitions an the convergence is locally exponential with respect to the Riemannian istance on G. Proof: Note that our conitions on the Riemannian metrics imply that ist(h y0 (X), ) = inf{ist(x, Z) Z stab( )}. Since the error Ê = ˆXX 1 projects to the output error on M, it converges to the fiber h 1 ( ) = stab( ) if an only if the output error converges to. By efinition h y0 : G M is the projection map of the fiber bunle G over M an hence any open, ense subset of M M whose complement has measure 2 0 lifts by the inverse of ( ˆX, X) (h y0 ( ˆX), h y0 (X)) to an open, ense set in G G whose complement has measure 0. Thus the generic convergence on M for the canonical error yiels the generic convergence on G for the error ˆXX 1. 2 with respect to the Riemannian measure

4 Note, that since the minimal realization of (1) is the projecte system on the output space, the subgroup stab( ) is unobservable an the error ˆXX 1 for any observer of (1) can only converge to the subgroup stab( ). Hence our observer provies the best possible convergence result for generic initial conitions. V. EXAMPLE In this section we will illustrate our observer construction on the attitue estimation problem. The attitue of a rigi boy, with boy-fixe-frame {B}, measure with respect to an inertial frame {A}, can be ientifie with an element R of SO(3). The left invariant ynamics Ṙ = Rω (6) on SO(3) correspon to the natural boy-fixe-frame kinematics of the system. Here ω = (ω 1, ω 2, ω 3 ) {B} is the angular velocity of the rigi-boy expresse in the boyfixe-frame an ω = 0 ω 3 ω 2 ω 3 0 ω 1 ω 2 ω 1 0. Consier the case where only partial information on the rigiboy attitue is measure. The most common situation is when an inertial irection, such as the magnetic fiel, or gravitational fiel, is measure by sensor systems such as magnetometers or accelerometers on boar a vehicle. The inertial irection {A} of the measure irection is known a-priori, for example, the gravitational irection lies in the z-axis of the inertial frame. The measure irection in the boy-fixe-frame, y = R {B}, is obtaine by using the attitue matrix to transform the inertial irection from the inertial frame into the boy fixe frame. The corresponing right action of SO(3) on S 2 is h : SO(3) S 2 S 2, h(r, y) = R y. This yiels the system Ṙ =Rω y =R (7a) (7b) on SO(3) with outputs on S 2. For the observer constructions propose in this paper, we have to choose a Riemannian metric on SO(3) an S 2. We equip the unit sphere with the Riemannian metric inuce by the Eucliean scalar prouct in R 3. To choose a Riemannian metric on SO(3) recall that the tangent spaces T R SO(3) are given by {RΩ Ω = Ω}. We use the Riemannian metric RΩ, RΘ = 1 2 tr(ωθ) on SO(3). The right-invariant horizontal istribution H(R) can be obtaine as the orthogonal complement of ker T e h y0, i.e. H(R) = T e R R (ker T e h y0 ). Straightforwar calculations show that these Riemannian metrics satisfy the conitions from Section II. By Theorem 1, the system (7) projects to a system on S 2. A straighforwar calculation shows that the projecte system can be written as ẏ = ω y, y(0) = h(x(0), ) (8) on S 2. First we esign an observer for the projecte system by the methos introuce in Section III. For constructing a graient-like observer on S 2 we have to choose a cost function. The most natural choice is the Eucliean istance, i.e. f(ŷ, y) = (k/2) ŷ y 2 with k > 0 a positive scaling constant. The constant k is introuce to provie a tunable gain for the observer. The graient of a function restricte to an embee submanifol of the Eucliean space with respect to the inuce metric is given by the Eucliean orthogonal projection of the Eucliean graient to the (embee) tangent space of the submanifol. For the unit sphere this projection is given by (I ŷŷ ). This yiels the graient gra 1 f(ŷ, y) = k(i ŷŷ )(y ŷ) = k ( I ŷŷ ) y. Thus we get from (3) the observer ŷ = ω ŷ + k(i ŷŷ )y (9) on S 2. It s easy to see that the cost function f satisfies the conitions of Theorem 4, i.e. that ŷ f(ŷ, y) is a Morse function with only one global minimum ( ) an no further local minima, an is invariant uner the action (R, (ŷ, y)) (R ŷ, R y) of SO(3) on S 2 S 2. Thus the canonical error on S 2 for this observer converges to for generic initial conitions. Some calculations show that k(i ŷŷ )y = k(ŷ y) ŷ. Hence, the observer (9) can be written as ŷ = ω ŷ + k(ŷ y) ŷ. (10) This observer was originally stuie by Metni et al. [10], [13]. The erivation given in earlier work is base on a Lyapunov analysis an inclues an integral term for compensation of gyro bias that oes not fit the analysis framework presente in this paper, however, the main proportional terms in [10], [13] are exactly (10). To obtain an observer on SO(3) we have to lift (9) to the group. Since the Riemmiann metrics on S 2 an SO(3) satisfy the conitions of Section II, we can apply Corollary 8 an construct the observer from a lift of the cost function to SO(3). The lifte cost function F of f is given by F ( ˆR, R) = k 2 ˆR R 2. To calculate the graient we nee the ifferential 1 of F with respect to the first variable. Some matrix calculations

5 yiel ( 1 F ( ˆR, R)( ˆRΩ) =k tr ˆRΩ( ˆR R) y0 ( k tr ( ˆR R)Ω ˆR ) y0 =k tr ( Ω(ŷ y)ŷ ) k tr (( ˆR R)Ωŷ(ŷ y) ) =2k tr ( Ω(ŷy yŷ ) ). From the ientity ) gra 1 F ( ˆR, R), RΩ = 2 tr ((gra 1 F ( ˆR, R)) ˆR Ω Ω ) = 1 F (R)(RΩ), = Ω, an the fact that ˆR gra1 F ( ˆR, R) is a skew symmetric matrix, one obtains the graient gra 1 F ( ˆR, R) = k ˆR(ŷy yŷ ) = k ˆR(ŷ y). Thus the graient-like observer on the Lie group is t ˆR = ˆRu k ˆR(ŷ y). Some further calculations show that the observer has the error ynamics t ( ˆRR ) = k ˆRR sk(r ˆR y 0 ) with sk(a) = 1 2 (A A ), which is the negative graient of f(e) = k 2 E 2. The observer is the explicit complementary filter propose in [14] excluing the integral introuce in that paper to compensate gyro bias. In particular, our observer esign approach illustrates the connection between the filters in [10], [13] an the complementary filter in [14]. VI. CONCLUSIONS In this paper we introuce a esign metho for constructing an observer for a left-invariant system on a Lie group with right invariant outputs on a manifol. We propose a graient-like observer for the projecte system on the output space which can be lifte to a graient-like observer on the group. The error for the observer for the projecte system converges to the reference output for generic initial states, with a locally exponential convergence rate. This lifts to convergence to the stabilizer of the reference output for the lifte observer on the Lie group, again with a locally exponential convergence rate. We have emonstrate this esign metho on the attitue estimation problem with partial state measurements for a rigi boy. In particular, this gave us new insights into the filters propose in [10], [13], [14]. VII. ACKNOWLEDGMENTS This paper presents research results of the Belgian Network DYSCO (Dynamical Systems, Control, an Optimization), fune by the Interuniversity Attraction Poles Programme, initiate by the Belgian State, Science Policy Office. The scientific responsibility rests with its authors. The research of the first author has also been partially supporte by the Australian Research Council (ARC) Centre of Excellence for Mathematics an Statistics of Complex Systems. The research of the secon an thir authors is partially supporte by the ARC through the Discovery project DP REFERENCES [1] R. Brockett, System theory on group manifols an coset spaces, SIAM J. Control, vol. 10, no. 2, pp , [2] V. Jurjevic an H. J. Sussmann, Control systems on Lie groups, Journal of Differential Equations, vol. 12, pp , [3] V. Jurjevic an I. Kupka, Control systems on semi-simple Lie groups an their homogeneous spaces, Annales e l institut Fourier, vol. 31, no. 4, pp , [4] D. Cheng, W. P. Dayawansa, an C. F. Martin, Observability of systems on Lie groups an coset spaces, SIAM J. Contr. Optim., vol. 28, pp , [5] N. Apostolou an D. Kazakos, On local observability for invariant systems on Lie groups an coset spaces, Proc. R. Soc. Lon., vol. 452, pp , [6] S. Salcuean, A globally convergent angular velocity observer for rigi boy motion, IEEE Transactions on Automatic Control, vol. 46, no. 12, pp , [7] H. Nijmeijer an T. I. Fossen, Es., New Directions in Nonlinear Observer Design, ser. Lecture Notes in Control an Information Sciences. Springer, [8] J. Thienel an R. M. Sanner, A couple nonlinear spacecraft attitue controller an observer with an unknow constant gyro bias an gyro noise, IEEE Transactions on Automatic Control, vol. 48, no. 11, pp , [9] D. Maithripala, J. Berg, an W. Dayawansa, Almost-global tracking of simple mechanical systems on a general class of Lie groups, IEEE Transactions on Automatic Control, vol. 51, no. 2, pp , [10] N. Metni, J.-M. Pflimlin, T. Hamel, an P. Soueres, Attitue an gyro bias estimation for a flying UAV, in Intelligent Robots an Systems, (IROS 2005) IEEE/RSJ International Conference on, 2-6 Aug. 2005, pp [11] R. Mahony, T. Hamel, an J.-M. Pflimlin, Complementary filter esign on the special orthogonal group SO(3), in Proceeings of the 44th Conference on Decision an Control, Seville, Spain, 2005, pp [12] S. Bonnabel, P. Martin, an P. Rouchon, A non-linear symmetrypreserving observer for velocity-aie inertial navigation, in Proceeings of the American Control Conference, June 2006, pp [13] N. Metni, J.-M. Pflimlin, T. Hamel, an P. Soueres, Attitue an gyro bias estimation for a VTOL UAV, Control Engineering Practice, vol. 14, no. 12, pp , December [14] R. Mahony, T. Hamel, an J.-M. Pflimlin, Non-linear complementary filters on the special orthogonal group, IEEE Transactions on Automatic Control, vol. 53, no. 5, pp , [15] A. Tayebi, S. McGilvray, A. Roberts, an M. Moallem, Attitue estimation an stabilization of a rigi boy using low-cost sensors, in Proceeings of the 46th IEEE Conference on Decision an Control, 2007, pp [16] P. Martin an E. Salaün, Invariant observers for attitue an heaing estimation from low-cost inertial an magnetic sensors, in Proceeings of the 46th IEEE Conference on Decision an Control, 2007, pp [17] J. Vasconcelos, C. Silvestre, an P. Oliveira, A nonlinear observer for rigi boy attitue estimation using vector observations, in Proceeings of the 17th Worl Congress The International Feeration of Automatic Control, Seoul, Korea, 2008, pp [18] B. Vik an T. Fossen, A nonlinear observer for GPS an INS integration, in Proceeings of the 40th IEEE Conference on Decision an Control, 2001, pp [19] H. Rehbiner an B. Ghosh, Pose estimation using line-base ynamic vision an inertial sensors, IEEE Transactions on Automatic Control, vol. 48, no. 2, pp , 2003.

6 [20] G. Balwin, R. Mahony, J. Trumpf, T. Hamel, an T. Cheviron, Complementary filter esign on the Special Eucliean group SE(3), in European Control Conference (ECC 2007), Kos, Greece, 2007, pp [21] J. Vasconcelos, R. Cunha, C. Silvestre, an P. Oliveira, Lanmark base nonlinear observer for rigi boy attitue an position estimation, in Proceeings of the 46th IEEE Conference on Decision an Control, 2007, pp [22] S. Bonnabel an P. Rouchon, Control an Observer Design for Nonlinear Finite an Infinite Dimensional Systems, ser. Lecture Notes in Control an Information Sciences. Springer-Verlag, 2005, vol. 322, ch. On Invariant Observers, pp [23] T. Hamel an R. Mahony, Attitue estimation on SO(3) base on irect inertial measurements, in Proceeings of the 2006 IEEE International Conference on Robotics an Automation. Orlano Fl., USA: Institue of Electrical an Electronic Engineers, April 2006, pp [24] S. Bonnabel, P. Martin, an P. Rouchon, Invariant asymptotic observers, IEEE Transactions on Automatic Control, 2008, accepte for publication. [Online]. Available: [25] C. Lageman, R. Mahony, an J. Trumpf, State observers for invariant ynamics on a Lie group, in Proceeings of the international conference on the Mathematical Theory of Networks an Systems, [26] S. Bonnabel, P. Martin, an P. Rouchon, Non-linear observer on Lie groups for left-invariant ynamics with right-left equivariant output, in Proceeings of the 17th Worl Congress The International Feeration of Automatic Control, Seoul, Korea, July [27] C. Lageman, J. Trumpf, an R. Mahony, Graient-like observers for invariant ynamics on a Lie group, Accepte for publication in IEEE Transactions on Automatic Control, [28] K. Nomizu, Invariant affine connections on homogeneous spaces, American Journal of Mathematics, vol. 76, no. 1, pp , [29] J. Cheeger an D. G. Ebin, Comparison theorems in Riemannian geometry. Amsteram: North-Hollan, [30] H. Nijmeijer an A. J. van er Schaft, Partial symmetries for nonlinear systems, Mathematical Systems Theory, vol. 18, pp , [31] U. Helmke an J. Moore, Optimization an ynamical systems. Springer, 1994.

State observers for invariant dynamics on a Lie group

State observers for invariant dynamics on a Lie group State observers for invariant dynamics on a Lie group C. Lageman, R. Mahony, J. Trumpf 1 Introduction This paper concerns the design of full state observers for state space systems where the state is evolving

More information

Bias Estimation for Invariant Systems on Lie Groups with Homogeneous Outputs

Bias Estimation for Invariant Systems on Lie Groups with Homogeneous Outputs Bias Estimation for Invariant Systems on Lie Groups with Homogeneous Outputs A. Khosravian, J. Trumpf, R. Mahony, C. Lageman Abstract In this paper, we provide a general method of state estimation for

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

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem Invariant Extene Kalman Filter: Theory an application to a velocity-aie estimation problem S. Bonnabel (Mines ParisTech) Joint work with P. Martin (Mines ParisTech) E. Salaun (Georgia Institute of Technology)

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

Measurement Observers for Pose Estimation on SE(3)

Measurement Observers for Pose Estimation on SE(3) Measurement Observers for Pose Estimation on SE(3) By Geoffrey Stacey u4308250 Supervised by Prof. Robert Mahony 24 September 2010 A thesis submitted in part fulfilment of the degree of Bachelor of Engineering

More information

IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. XX, MONTH YEAR 1

IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. XX, MONTH YEAR 1 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. XX, MONTH YEAR 1 Non-linear complementary filters on the special orthogonal group Robert Mahony, Member, IEEE, Tarek Hamel, Member, IEEE, an Jean-Michel

More information

Unit quaternion observer based attitude stabilization of a rigid spacecraft without velocity measurement

Unit quaternion observer based attitude stabilization of a rigid spacecraft without velocity measurement Proceedings of the 45th IEEE Conference on Decision & Control Manchester Grand Hyatt Hotel San Diego, CA, USA, December 3-5, 6 Unit quaternion observer based attitude stabilization of a rigid spacecraft

More information

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

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

More information

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

SYMPLECTIC GEOMETRY: LECTURE 3

SYMPLECTIC GEOMETRY: LECTURE 3 SYMPLECTIC GEOMETRY: LECTURE 3 LIAT KESSLER 1. Local forms Vector fiels an the Lie erivative. A vector fiel on a manifol M is a smooth assignment of a vector tangent to M at each point. We think of M as

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

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

A Complementary Filter for Attitude Estimation of a Fixed-Wing UAV

A Complementary Filter for Attitude Estimation of a Fixed-Wing UAV A Complementary Filter for Attitude Estimation of a Fixed-Wing UAV Mark Euston, Paul Coote, Robert Mahony, Jonghyuk Kim and Tarek Hamel Abstract This paper considers the question of using a nonlinear complementary

More information

Geometrical methods in observer design: observers and symmetries

Geometrical methods in observer design: observers and symmetries Geometrical methos in observer esign: observers an symmetries Silvère Bonnabel Coworkers: Pierre Rouchon, Philippe Martin, Erwan Salaun Mines ParisTech Centre e Robotique Mathématiques et Systèmes silvere.bonnabel@mines-paristech.fr

More information

Systems & Control Letters

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

More information

Geometric Nonlinear Observer Design for SLAM on a Matrix Lie Group

Geometric Nonlinear Observer Design for SLAM on a Matrix Lie Group 218 IEEE Conference on Decision and Control (CDC Miami Beach, FL, USA, Dec. 17-19, 218 Geometric Nonlinear Observer Design for SLAM on a Matrix Lie Group Miaomiao Wang and Abdelhamid Tayebi Abstract This

More information

Calculus of Variations

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

More information

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

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

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

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

More information

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

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10

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

More information

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

SOME RESULTS ON THE GEOMETRY OF MINKOWSKI PLANE. Bing Ye Wu

SOME RESULTS ON THE GEOMETRY OF MINKOWSKI PLANE. Bing Ye Wu ARCHIVUM MATHEMATICUM (BRNO Tomus 46 (21, 177 184 SOME RESULTS ON THE GEOMETRY OF MINKOWSKI PLANE Bing Ye Wu Abstract. In this paper we stuy the geometry of Minkowski plane an obtain some results. We focus

More information

Feature-based Recursive Observer Design for Homography Estimation

Feature-based Recursive Observer Design for Homography Estimation Feature-base Recursive Observer Design for Homography Estimation Minh-Duc Hua, Jochen Trumpf, Tarek Hamel, Robert Mahony, an Pascal Morin arxiv:66.3v [cs.cv] 9 Jun 6 Abstract This paper presents a new

More information

Euler equations for multiple integrals

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

More information

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

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

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

More information

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

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

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

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

Contraction and Observer Design on Cones

Contraction and Observer Design on Cones Contraction an Observer Design on Cones Silvère Bonnabel, Alessanro Astolfi, an Roolphe Sepulchre Abstract We consier the problem of positive observer esign for positive systems efine on soli cones in

More information

there is no special reason why the value of y should be fixed at y = 0.3. Any y such that

there is no special reason why the value of y should be fixed at y = 0.3. Any y such that 25. More on bivariate functions: partial erivatives integrals Although we sai that the graph of photosynthesis versus temperature in Lecture 16 is like a hill, in the real worl hills are three-imensional

More information

Nöether s Theorem Under the Legendre Transform by Jonathan Herman

Nöether s Theorem Under the Legendre Transform by Jonathan Herman Nöether s Theorem Uner the Legenre Transform by Jonathan Herman A research paper presente to the University of Waterloo in fulfilment of the research paper requirement for the egree of Master of Mathematics

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

WUCHEN LI AND STANLEY OSHER

WUCHEN LI AND STANLEY OSHER CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability

More information

LECTURE 1: BASIC CONCEPTS, PROBLEMS, AND EXAMPLES

LECTURE 1: BASIC CONCEPTS, PROBLEMS, AND EXAMPLES LECTURE 1: BASIC CONCEPTS, PROBLEMS, AND EXAMPLES WEIMIN CHEN, UMASS, SPRING 07 In this lecture we give a general introuction to the basic concepts an some of the funamental problems in symplectic geometry/topology,

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

Darboux s theorem and symplectic geometry

Darboux s theorem and symplectic geometry Darboux s theorem an symplectic geometry Liang, Feng May 9, 2014 Abstract Symplectic geometry is a very important branch of ifferential geometry, it is a special case of poisson geometry, an coul also

More information

Consensus on homogeneous manifolds

Consensus on homogeneous manifolds Joint 48th IEEE Conference on Decision an Control an 28th Chinese Control Conference Shanghai, P.R. China, December 16-18, 2009 FrA05.1 Consensus on homogeneous manifols Alain Sarlette Roolphe Sepulchre

More information

Automated Tuning of the Nonlinear Complementary Filter for an Attitude Heading Reference Observer

Automated Tuning of the Nonlinear Complementary Filter for an Attitude Heading Reference Observer Automated Tuning of the Nonlinear Complementary Filter for an Attitude Heading Reference Observer Oscar De Silva, George K.I. Mann and Raymond G. Gosine Faculty of Engineering and Applied Sciences, Memorial

More information

Automatica. Global trajectory tracking control of VTOL-UAVs without linear velocity measurements

Automatica. Global trajectory tracking control of VTOL-UAVs without linear velocity measurements Automatica 46 (200) 053 059 Contents lists available at ScienceDirect Automatica journal homepage: www.elsevier.com/locate/automatica Brief paper Global trajectory tracking control of VTOL-UAVs without

More information

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

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

More information

arxiv: v1 [math.dg] 30 May 2012

arxiv: v1 [math.dg] 30 May 2012 VARIATION OF THE ODUUS OF A FOIATION. CISKA arxiv:1205.6786v1 [math.dg] 30 ay 2012 Abstract. The p moulus mo p (F) of a foliation F on a Riemannian manifol is a generalization of extremal length of plane

More information

Adaptive Estimation of Measurement Bias in Six Degree of Freedom Inertial Measurement Units: Theory and Preliminary Simulation Evaluation

Adaptive Estimation of Measurement Bias in Six Degree of Freedom Inertial Measurement Units: Theory and Preliminary Simulation Evaluation Adaptive Estimation of Measurement Bias in Six Degree of Freedom Inertial Measurement Units: Theory and Preliminary Simulation Evaluation Andrew R. Spielvogel and Louis L. Whitcomb Abstract Six-degree

More information

Notes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk

Notes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk Notes on Lie Groups, Lie algebras, an the Exponentiation Map Mitchell Faulk 1. Preliminaries. In these notes, we concern ourselves with special objects calle matrix Lie groups an their corresponing Lie

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

Energy-preserving affine connections

Energy-preserving affine connections 2 A. D. Lewis Enery-preservin affine connections Anrew D. Lewis 28/07/1997 Abstract A Riemannian affine connection on a Riemannian manifol has the property that is preserves the kinetic enery associate

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

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

LECTURE NOTES ON DVORETZKY S THEOREM

LECTURE NOTES ON DVORETZKY S THEOREM LECTURE NOTES ON DVORETZKY S THEOREM STEVEN HEILMAN Abstract. We present the first half of the paper [S]. In particular, the results below, unless otherwise state, shoul be attribute to G. Schechtman.

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

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

Lower bounds on Locality Sensitive Hashing

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

More information

IN the recent past, the use of vertical take-off and landing

IN the recent past, the use of vertical take-off and landing IEEE TRANSACTIONS ON ROBOTICS, VOL. 27, NO. 1, FEBRUARY 2011 129 Aaptive Position Tracking of VTOL UAVs Anrew Roberts, Stuent Member, IEEE, an Abelhami Tayebi, Senior Member, IEEE Abstract An aaptive position-tracking

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

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

Visual Servoing for Underactuated VTOL UAVs : a Linear, Homography-Based Framework

Visual Servoing for Underactuated VTOL UAVs : a Linear, Homography-Based Framework Visual Servoing for Uneractuate VTOL UAVs : a Linear, Homography-Base Framework Henry e Plinval, Pascal Morin, Philippe Mouyon, Tarek Hamel H. e Plinval an P. Mouyon are with ONERA-The French Aerospace

More information

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

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

More information

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1 Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 PHYS 500 - Southern Illinois University September 8, 2016 PHYS 500 - Southern Illinois University Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 September 8,

More information

The Variational Attitude Estimator in the Presence of Bias in Angular Velocity Measurements

The Variational Attitude Estimator in the Presence of Bias in Angular Velocity Measurements The Variational Attitude Estimator in the Presence of Bias in Angular Velocity Measurements Maziar Izadi 1, Sasi Prabhakaran 1, Amit Sanyal 2,, Carlos Silvestre 3, and Paulo Oliveira 4 Abstract Estimation

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

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

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

Characteristic classes of vector bundles

Characteristic classes of vector bundles Characteristic classes of vector bunles Yoshinori Hashimoto 1 Introuction Let be a smooth, connecte, compact manifol of imension n without bounary, an p : E be a real or complex vector bunle of rank k

More information

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

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

More information

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

ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS

ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS ADRIAN CONSTANTIN AND BORIS KOLEV Abstract. Accoring to the principle of least action, the spatially perioic motions of one-imensional

More information

The nonsmooth Newton method on Riemannian manifolds

The nonsmooth Newton method on Riemannian manifolds The nonsmooth Newton method on Riemannian manifolds C. Lageman, U. Helmke, J.H. Manton 1 Introduction Solving nonlinear equations in Euclidean space is a frequently occurring problem in optimization and

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

Minimum-time constrained velocity planning

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

More information

Function Spaces. 1 Hilbert Spaces

Function Spaces. 1 Hilbert Spaces Function Spaces A function space is a set of functions F that has some structure. Often a nonparametric regression function or classifier is chosen to lie in some function space, where the assume structure

More information

Witten s Proof of Morse Inequalities

Witten s Proof of Morse Inequalities Witten s Proof of Morse Inequalities by Igor Prokhorenkov Let M be a smooth, compact, oriente manifol with imension n. A Morse function is a smooth function f : M R such that all of its critical points

More information

IMU-Laser Scanner Localization: Observability Analysis

IMU-Laser Scanner Localization: Observability Analysis IMU-Laser Scanner Localization: Observability Analysis Faraz M. Mirzaei and Stergios I. Roumeliotis {faraz stergios}@cs.umn.edu Dept. of Computer Science & Engineering University of Minnesota Minneapolis,

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

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const. G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

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

More information

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

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

More information

Implicit Differentiation

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

More information

Chapter 2 Lagrangian Modeling

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

More information

Relative Position Sensing by Fusing Monocular Vision and Inertial Rate Sensors

Relative Position Sensing by Fusing Monocular Vision and Inertial Rate Sensors Proceeings of ICAR 2003 The 11th International Conference on Avance Robotics Coimbra, Portugal, June 30 - July 3, 2003 Relative Position Sensing by Fusing Monocular Vision an Inertial Rate Sensors Anreas

More information

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

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

More information

From Local to Global Control

From Local to Global Control Proceeings of the 47th IEEE Conference on Decision an Control Cancun, Mexico, Dec. 9-, 8 ThB. From Local to Global Control Stephen P. Banks, M. Tomás-Roríguez. Automatic Control Engineering Department,

More information

Accurate and Robust Ego-Motion Estimation using Expectation Maximization

Accurate and Robust Ego-Motion Estimation using Expectation Maximization Accurate an Robust Ego-Motion Estimation using Expectation Maximization Gijs Dubbelman 1,, Wannes van er Mark 1 an Frans C.A. Groen Abstract A novel robust visual-oometry technique, calle EM-SE(3) is presente

More information

Nonlinear Wind Estimator Based on Lyapunov

Nonlinear Wind Estimator Based on Lyapunov Nonlinear Based on Lyapunov Techniques Pedro Serra ISR/DSOR July 7, 2010 Pedro Serra Nonlinear 1/22 Outline 1 Motivation Problem 2 Aircraft Dynamics Guidance Control and Navigation structure Guidance Dynamics

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

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

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

arxiv: v1 [math.dg] 1 Nov 2015

arxiv: v1 [math.dg] 1 Nov 2015 DARBOUX-WEINSTEIN THEOREM FOR LOCALLY CONFORMALLY SYMPLECTIC MANIFOLDS arxiv:1511.00227v1 [math.dg] 1 Nov 2015 ALEXANDRA OTIMAN AND MIRON STANCIU Abstract. A locally conformally symplectic (LCS) form is

More information

Chromatic number for a generalization of Cartesian product graphs

Chromatic number for a generalization of Cartesian product graphs Chromatic number for a generalization of Cartesian prouct graphs Daniel Král Douglas B. West Abstract Let G be a class of graphs. The -fol gri over G, enote G, is the family of graphs obtaine from -imensional

More information

model considered before, but the prey obey logistic growth in the absence of predators. In

model considered before, but the prey obey logistic growth in the absence of predators. In 5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an

More information

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition

More information

Convergence of Random Walks

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

More information

Optimal Control of Spatially Distributed Systems

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

More information

The Principle of Least Action

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

More information

Relation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function

Relation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives

More information

Nonlinear attitude estimation with measurement decoupling and anti-windup gyro-bias compensation

Nonlinear attitude estimation with measurement decoupling and anti-windup gyro-bias compensation Nonlinear attitude estimation with measurement decoupling and anti-windup gyro-bias compensation Minh-Duc Hua, Konrad Rudin, Guillaume Ducard, Tarek Hamel, Robert Mahony I3S UNS CNRS, Sophia Antipolis,

More information

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

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

More information