Direct Methods of Optimal Control for Space Trajectory Optimization

Size: px
Start display at page:

Download "Direct Methods of Optimal Control for Space Trajectory Optimization"

Transcription

1 Direc Mehods o Opimal Conrol or Space Trajecory Opimiaion Rober G. Melon Deparmen o Aerospace Engineering Compuaional Science Seminar February 6, 7

2 Trajecory Opimiaion A B (min. ime and/or propellan) subjec o consrains Venus Earh Sun

3 Solar Sail Propulsion Use phoon pressure on large relecive sail n F ( + r) PA c ( r n) n releciviy r r ideal sail spacecra

4 Realisic sail

5 A simple approach Discreie eqs o moion in ime ( j, j,n) Unknowns: conrols u j a j ( known?) Esablish perormance measure J arge, sae (pos., vel.) Inegrae eqs o moion; use shooing mehod and LP o adjus he u j : minimie J LP nonlinear programming minimie (u) subjec o consrains Too sensiive o small changes in u j

6 Opimal Conrol Problem subjec o: minimie J h(,,, ) + g( ( ), u ( ), ) d ( ) ( ( ), u( ), ) b L ( ) g ( ), ( ),, ) b L ( g( ( ), u( ), ) b U b U and simple bounds on he saes and conrols

7 Indirec Mehods Form Hamilionian H g(, u, ) + p T ecessary condiions or opimaliy o u: opimal )* ( ) * ( ) *, *, *, ( ) *( ) *, ( ) *, *, *, ( cosae ) *, *, *, ( * sae ) *, *, *, ( * + + T h H h H H H δ δ p u p p u u p p u p p u p Diicul o ormulae necessary condiions Soln o eqs requires muliple ieraions; sensiive o i.c. s

8 Direc Collocaion wih LP (DCLP) Use implici soln o eqs o moion o ormulaion o necessary condiions rqd. Bu require soln o large-scale LP wih large number o equaliy consrains Discreie eqs o moion in ime ( j, j,n n number o nodes) Unknowns: conrols u j and saes j a nodes ( known?) Use splines or j, inerpolaion or u j min J ( c) c [ subjec o : Δ( c) M deecs Δ ( c) n w w( c) w L L ( c,, L, n ), u U U, L, u, inerior consrains erminal consrains and simple bounds on he conrols n ]

9 Calculaion o he Deecs slope ds/d Δ j d i /d -ds/d i slope d i /d i (,u,) cener poin node j node j+ spline S() j c j+

10 A Pseudospecral Mehod (also a collocaion mehod) Epands and u as inie sums o Lagrange polynomials evaluaed a Chebyshe-Gauss- Lobao (CGL) poins k cos(πk/) k,, [-,], scaled ime True imeτ [(τ -τ ) + (τ +τ )]/ ( ) φ ( ) u ( ) u φ ( ) j j j j j j where φ j Lagrange inerpolaing polynomial

11 ow solve he LP (variables: k, u k, and possibly τ ) ) (, ),,, ) ( /( subjec o /,,, cos 4 4 ),, ( ), ( " / k k U k k k k L j s s k k k k k s js j w w w w w g h J d u d u L + τ τ τ π τ τ τ τ

12 Solar Sail Problems Conrol auhoriy srong uncion o sail orienaion (angles α, δ ) D case easily solved by several direc opimiaion mehods (and including realisic orce models) 3D case -- plane change creaes challenges or some mehods Eample: compare perormance o wo mehods: DCLP and pseudospecral (DIDO)

13 Eample Cases Ideal sail a T PA mc ( ) n r n A 3 m, m 5 kg a T.5 mm/s a au rom Sun Hohmann-like ransers (no deparure or capure deails) Planeary orbis assumed circular D cases (Earh and Venus in same plane) 3D cases (plane change required)

14 Sail angles: Consrains π α, δ π OK or DCLP, no or DIDO Trig. repres. o crls: cosα, sinα and cosδ, sinδ), wih addiional consrain cos ( ) + sin ( ) Simple bounds cos( ), sin( ).

15 3D Final Consrains + r V ( ) r V μ θ + + ( ) i V cos θ θ θ ( ) ( ) Ω V cos sin cos + + θ θ θ θ ( ) ( ) Ω V sin cos sin + θ θ θ θ

16 D case wih nodes -5 - DIDO DCLP -5-3 Saii angle α (deg.) Time (years)

17 D case wih nodes (DIDO) Venus' orbi

18 D case wih nodes (DCLP) Venus' orbi

19 D case wih nodes -5 - DIDO DCLP -5-3 Saii angle α (deg.) Time (years)

20 D case wih nodes (DIDO) Venus' orbi

21 D case wih nodes (DCLP) Venus' orbi

22 Sail Angles or 3D case wih nodes 3 α (in-plane) δ (ou-o-plane) Saii angles (deg.) Time (years)

23 Trajecory or 3D case wih nodes.5 Earh Venus

24 Sail angles or 3D case wih nodes 3 α (in-plane) δ (ou-o-plane) Saii angles (deg.) Time (years)

25 Trajecory or 3D case wih nodes.5 Earh Venus

26 Processing Time Eac comparison no possible (dieren processors) D cases DCLP: 4 DIDO:.67 (sparse mari algorihms) 3D cases DCLP: eremely sensiive o iniial guess -- had o use perurbed DIDO soluion as irs guess o obain convergence DIDO:.9

27 Conclusions DIDO less sensiive o iniial guess or o addiional inal consrains in 3D problem eed a horough benchmark sudy suie o problems wih variey o consrains run under idenical condiions useul guide or selecing mehod appropriae or given ype o rajecory problem

The Pennsylvania State University. The Graduate School. Department of Aerospace Engineering APPLICATION OF A HOMOTOPY METHOD

The Pennsylvania State University. The Graduate School. Department of Aerospace Engineering APPLICATION OF A HOMOTOPY METHOD The Pennsylvania Sae Universiy The Graduae School Deparmen o Aerospace Engineering APPLICATION OF A HOMOTOPY METHOD TO LOW-THRUST TRAJECTORY OPTIMIZATION A Thesis in Aerospace Engineering by Chenmeng Tu

More information

Symbolic-Numeric Sparse Interpolation of Multivariate Polynomials

Symbolic-Numeric Sparse Interpolation of Multivariate Polynomials Symbolic-Numeric Sparse Inerpolaion o Mulivariae Polynomials Mark Giesbrech George Labahn Wen-shin Lee School o Compuer Science Universiy o Waerloo Approimae inerpolaion o polynomials Suppose we can sample

More information

TRAJECTORY SHAPING OF SURFACE-TO-SURFACE MISSILE WITH TERMINAL IMPACT ANGLE CONSTRAINT

TRAJECTORY SHAPING OF SURFACE-TO-SURFACE MISSILE WITH TERMINAL IMPACT ANGLE CONSTRAINT MAKARA, TEKNOLOGI, VOL., NO., NOVEMBER 7: 65-7 TRAJECTORY SHAPING OF SURFACE-TO-SURFACE MISSILE WITH TERMINAL IMPACT ANGLE CONSTRAINT S. Subchan Deparmen o Mahemaics, Insiu Teknologi Sepuluh Nopember,

More information

1 Subdivide the optimization horizon [t 0,t f ] into n s 1 control stages,

1 Subdivide the optimization horizon [t 0,t f ] into n s 1 control stages, Opimal Conrol Formulaion Opimal Conrol Lecures 19-2: Direc Soluion Mehods Benoî Chachua Deparmen of Chemical Engineering Spring 29 We are concerned wih numerical soluion procedures for

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

CS537. Numerical Analysis

CS537. Numerical Analysis CS57 Numerical Analsis Lecure Numerical Soluion o Ordinar Dierenial Equaions Proessor Jun Zang Deparmen o Compuer Science Universi o enuck Leingon, Y 4006 0046 April 5, 00 Wa is ODE An Ordinar Dierenial

More information

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN

More information

On Solving the Perturbed Multi- Revolution Lambert Problem: Applications in Enhanced SSA

On Solving the Perturbed Multi- Revolution Lambert Problem: Applications in Enhanced SSA On Solving he Perurbed Muli- Revoluion Lamber Problem: Applicaions in Enhanced SSA John L. Junkins and Robyn M. Woollands Texas A&M Universiy Presened o Sacie Williams (AFOSR/RT) AFOSR REMOTE SENSING PORTFOLIO

More information

Low-thrust Trajectory Optimization Yuan Ren

Low-thrust Trajectory Optimization Yuan Ren Low-hrus rajecory Opimizaion Yuan Ren Deparamen de Maemica Aplicada I Universia Policnica de Caalunya Index he Acion o rajecory Opimizaion Calculus o Variaions & Ponryagin s Minimum Principle Dynamical

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

I. OBJECTIVE OF THE EXPERIMENT.

I. OBJECTIVE OF THE EXPERIMENT. I. OBJECTIVE OF THE EXPERIMENT. Swissmero raels a high speeds hrough a unnel a low pressure. I will hereore undergo ricion ha can be due o: ) Viscosiy o gas (c. "Viscosiy o gas" eperimen) ) The air in

More information

Space-time Galerkin POD for optimal control of Burgers equation. April 27, 2017 Absolventen Seminar Numerische Mathematik, TU Berlin

Space-time Galerkin POD for optimal control of Burgers equation. April 27, 2017 Absolventen Seminar Numerische Mathematik, TU Berlin Space-ime Galerkin POD for opimal conrol of Burgers equaion Manuel Baumann Peer Benner Jan Heiland April 27, 207 Absolvenen Seminar Numerische Mahemaik, TU Berlin Ouline. Inroducion 2. Opimal Space Time

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

DYNAMICS AND REAL-TIME OPTIMAL CONTROL OF SATELLITE ATTITUDE AND SATELLITE FORMATION SYSTEMS

DYNAMICS AND REAL-TIME OPTIMAL CONTROL OF SATELLITE ATTITUDE AND SATELLITE FORMATION SYSTEMS DYAMICS AD REAL-TIME OPTIMAL COTROL OF SATELLITE ATTITUDE AD SATELLITE FORMATIO SYSTEMS A Disseraion by HUI YA Submied o he Oice o Graduae Sudies o Texas A&M Universiy in parial ulillmen o he requiremens

More information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

Trajectory planning in Cartesian space

Trajectory planning in Cartesian space Roboics 1 Trajecory planning in Caresian space Prof. Alessandro De Luca Roboics 1 1 Trajecories in Caresian space in general, he rajecory planning mehods proposed in he join space can be applied also in

More information

OPTIMAL FEEDBACK GUIDANCE ALGORITHMS FOR PLANETARY LANDING AND ASTEROID INTERCEPT

OPTIMAL FEEDBACK GUIDANCE ALGORITHMS FOR PLANETARY LANDING AND ASTEROID INTERCEPT AAS 11-588 OPIMAL FEEDBACK GUIDANCE ALGORIHMS FOR PLANEARY LANDING AND ASEROID INERCEP Yanning Guo, * Ma Hawkins, and Bong Wie INRODUCION Opimal eedback guidance alrihms are invesigaed or planeary landing

More information

Optimal trajectory generation for wheeled mobile robot

Optimal trajectory generation for wheeled mobile robot Opimal rajecory generaion or wheeled mobile robo Bin Qin Yeng Chai Soh Ming Xie Danwei Wang Nanyang echnological Universiy, Nangyang Avenue, Singapore mbqin@nu.edu.sg, eycsoh@nu.edu.sg Absrac. he problem

More information

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,

More information

Linear Optimal Missile Guidance using LQR

Linear Optimal Missile Guidance using LQR Opimal Conrol, Guidance and Esimaion Lecure 18 Linear Opimal issile Guidance using LQR Pro. Radhakan Padhi Dep. o Aerospace Engineering Indian Insiue o Science - Bangalore Topics Basic Fundamenals Philosophy

More information

CS376 Computer Vision Lecture 6: Optical Flow

CS376 Computer Vision Lecture 6: Optical Flow CS376 Compuer Vision Lecure 6: Opical Flow Qiing Huang Feb. 11 h 2019 Slides Credi: Krisen Grauman and Sebasian Thrun, Michael Black, Marc Pollefeys Opical Flow mage racking 3D compuaion mage sequence

More information

Stochastic Model Predictive Control for Gust Alleviation during Aircraft Carrier Landing

Stochastic Model Predictive Control for Gust Alleviation during Aircraft Carrier Landing Ouline Moivaion Sochasic formulaion Sochasic Model Predicive Conrol for Gus Alleviaion during Aircraf Carrier Landing Aircraf and gus modeling Resuls Deparmen of Mechanical and Aerospace Engineering Rugers,

More information

We here collect a few numerical tests, in order to put into evidence the potentialities of HBVMs [4, 6, 7].

We here collect a few numerical tests, in order to put into evidence the potentialities of HBVMs [4, 6, 7]. Chaper Numerical Tess We here collec a few numerical ess, in order o pu ino evidence he poenialiies of HBVMs [4, 6, 7]. Tes problem Le us consider he problem characerized by he polynomial Hamilonian (4.)

More information

Book Corrections for Optimal Estimation of Dynamic Systems, 2 nd Edition

Book Corrections for Optimal Estimation of Dynamic Systems, 2 nd Edition Boo Correcions for Opimal Esimaion of Dynamic Sysems, nd Ediion John L. Crassidis and John L. Junins November 17, 017 Chaper 1 This documen provides correcions for he boo: Crassidis, J.L., and Junins,

More information

Human-Robot Cooperative Manipulation with Motion Estimation

Human-Robot Cooperative Manipulation with Motion Estimation Human-Robo Cooperaive Manipulaion wih Moion Esimaion Yusuke MAEDA, Takayuki HARA and Tamio ARAI (The Universiy o Tokyo) 1. Inroducion 2. Virual Compliance Conrol 3. Esimaion o Human Moion 4. Eperimens

More information

Numerical Optimization

Numerical Optimization Numerical Opimizaion A Workshop A Deparmen o Mahemaics Chiang Mai Universiy Augus 4-5, 9 Insrucor: Elecrical Engineering and Compuer Science Case Wesern Reserve Universiy Phone: 6 368 454, Fax: 6 368 33

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

A METHOD TO SOLVE MISSILE-AIRCRAFT PURSUIT-EVASION DIFFERENTIAL GAMES. Fumiaki Imado* and Takeshi Kuroda**

A METHOD TO SOLVE MISSILE-AIRCRAFT PURSUIT-EVASION DIFFERENTIAL GAMES. Fumiaki Imado* and Takeshi Kuroda** A MEHOD O SOLVE MSSLE-ARCRAF PURSU-EVASON DFFERENAL GAMES Fumiaki mado and akeshi Kuroda Deparmen o Mechanical Sysems Engineering, Shinshu Universiy 4-17-1 Wakasao Nagano, Nagano 38-8553 JAPAN Advanced

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Math Final Exam Solutions

Math Final Exam Solutions Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,

More information

High-Order State and Parameter Transition Tensor Calculations

High-Order State and Parameter Transition Tensor Calculations AIAA-8-6453 Presened a: AIAA/AAS Asrodynamics Specialis Conerence Hawaii Convenion Cener and Hilon Hawaiian Village, Honolulu, Hawaii 18-1 Aug 8 High-Order Sae and Parameer Transiion Tensor Calculaions

More information

y = (y 1)*(y 3) t

y = (y 1)*(y 3) t MATH 66 SPR REVIEW DEFINITION OF SOLUTION A funcion = () is a soluion of he differenial equaion d=d = f(; ) on he inerval ff < < fi if (d=d)() =f(; ()) for each so ha ff

More information

Optimizing heat exchangers

Optimizing heat exchangers Opimizing hea echangers Jean-Luc Thiffeaul Deparmen of Mahemaics, Universiy of Wisconsin Madison, 48 Lincoln Dr., Madison, WI 5376, USA wih: Florence Marcoe, Charles R. Doering, William R. Young (Daed:

More information

Probabilistic Robotics

Probabilistic Robotics Probabilisic Roboics Bayes Filer Implemenaions Gaussian filers Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel Gaussians : ~ π e p N p - Univariae / / : ~ μ μ μ e p Ν p d π Mulivariae

More information

A Variation Evolving Method for Optimal Control

A Variation Evolving Method for Optimal Control A Variaion Evolving Mehod or Opimal Conrol Sheng ZHANG, En-Mi YONG, Wei-Qi QIAN, and Kai-Feng HE (7. Absrac: A new mehod or he opimal soluions is proposed. Originaing rom he coninuous-ime dynamics sabiliy

More information

Contact dynamics via implicit complementarity

Contact dynamics via implicit complementarity Conac ynamics ia implici complemenariy Emo ooro Uniersiy o Washingon ICRA 200 Emo ooro (Uniersiy o Washingon) Conac ynamics ia implici complemenariy ICRA 200 / 0 Muli-join ynamics wih conac Coninuous-ime

More information

New Seven-Step Numerical Method for Direct Solution of Fourth Order Ordinary Differential Equations

New Seven-Step Numerical Method for Direct Solution of Fourth Order Ordinary Differential Equations 9 J. Mah. Fund. Sci., Vol. 8, No.,, 9-5 New Seven-Sep Numerical Mehod for Direc Soluion of Fourh Order Ordinary Differenial Equaions Zurni Omar & John Olusola Kuboye Deparmen of Mahemaics, School of Quaniaive

More information

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I Open loop vs Closed Loop Advanced I Moor Command Movemen Overview Open Loop vs Closed Loop Some examples Useful Open Loop lers Dynamical sysems CPG (biologically inspired ), Force Fields Feedback conrol

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Probabilistic Robotics SLAM

Probabilistic Robotics SLAM Probabilisic Roboics SLAM The SLAM Problem SLAM is he process by which a robo builds a map of he environmen and, a he same ime, uses his map o compue is locaion Localizaion: inferring locaion given a map

More information

FLIGHT PATH OPTIMAL INTERVAL CONTROL FOR UNMANNED AERIAL VEHICLE

FLIGHT PATH OPTIMAL INTERVAL CONTROL FOR UNMANNED AERIAL VEHICLE UDC 629.7 FLIGHT PATH OPTIMAL INTERVAL CONTROL FOR UNMANNED AERIAL VEHICLE PH.D Proessor A.A.Lobay, PH.D Suden M.A.AL-Mashhadani Belarusian Naional Technical Universiy, Minsk, Belarus Davidmar3@mail.ru

More information

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx. . Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.

More information

Fractional Method of Characteristics for Fractional Partial Differential Equations

Fractional Method of Characteristics for Fractional Partial Differential Equations Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics

More information

Introduction to Mobile Robotics

Introduction to Mobile Robotics Inroducion o Mobile Roboics Bayes Filer Kalman Filer Wolfram Burgard Cyrill Sachniss Giorgio Grisei Maren Bennewiz Chrisian Plagemann Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel

More information

Then. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k)

Then. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k) Ger sgorin Circle Chaper 9 Approimaing Eigenvalues Per-Olof Persson persson@berkeley.edu Deparmen of Mahemaics Universiy of California, Berkeley Mah 128B Numerical Analysis (Ger sgorin Circle) Le A be

More information

Probabilistic Robotics SLAM

Probabilistic Robotics SLAM Probabilisic Roboics SLAM The SLAM Problem SLAM is he process by which a robo builds a map of he environmen and, a he same ime, uses his map o compue is locaion Localizaion: inferring locaion given a map

More information

Haar Wavelet Operational Matrix Method for Solving Fractional Partial Differential Equations

Haar Wavelet Operational Matrix Method for Solving Fractional Partial Differential Equations Copyrigh 22 Tech Science Press CMES, vol.88, no.3, pp.229-243, 22 Haar Wavele Operaional Mari Mehod for Solving Fracional Parial Differenial Equaions Mingu Yi and Yiming Chen Absrac: In his paper, Haar

More information

From Particles to Rigid Bodies

From Particles to Rigid Bodies Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and

More information

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.7. PLANE LINEAR ELASTICITY Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Plane Linear Elasici Theor Plane Sress Simplifing Hpohesis Srain Field Consiuive Equaion Displacemen Field The Linear

More information

Distance Between Two Ellipses in 3D

Distance Between Two Ellipses in 3D Disance Beween Two Ellipses in 3D David Eberly Magic Sofware 6006 Meadow Run Cour Chapel Hill, NC 27516 eberly@magic-sofware.com 1 Inroducion An ellipse in 3D is represened by a cener C, uni lengh axes

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

Physics Notes - Ch. 2 Motion in One Dimension

Physics Notes - Ch. 2 Motion in One Dimension Physics Noes - Ch. Moion in One Dimension I. The naure o physical quaniies: scalars and ecors A. Scalar quaniy ha describes only magniude (how much), NOT including direcion; e. mass, emperaure, ime, olume,

More information

Intermediate Differential Equations Review and Basic Ideas

Intermediate Differential Equations Review and Basic Ideas Inermediae Differenial Equaions Review and Basic Ideas John A. Burns Cener for Opimal Design And Conrol Inerdisciplinary Cener forappliedmahemaics Virginia Polyechnic Insiue and Sae Universiy Blacksburg,

More information

The LUNAR ORBIT PHASING AND RENDEZVOUS TESTS OF CHANG E-5 REENTRY FLIGHT EXPERIMENT VEHICLE

The LUNAR ORBIT PHASING AND RENDEZVOUS TESTS OF CHANG E-5 REENTRY FLIGHT EXPERIMENT VEHICLE The LUNAR ORBIT PHASING AND RENDEZVOUS TESTS OF CHANG E-5 REENTRY FLIGHT EXPERIMENT VEHICLE Gefei Li, Jianfeng Xie, Jianliang Zhou, Jun Song Beijing Aerospace Conrol Cener, POB.510, NO.105, Beijing, China,

More information

Basilio Bona ROBOTICA 03CFIOR 1

Basilio Bona ROBOTICA 03CFIOR 1 Indusrial Robos Kinemaics 1 Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function Elecronic Journal of Qualiaive Theory of Differenial Equaions 13, No. 3, 1-11; hp://www.mah.u-szeged.hu/ejqde/ Exisence of posiive soluion for a hird-order hree-poin BVP wih sign-changing Green s funcion

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

On a State-Constrained PDE Optimal Control Problem arising from ODE-PDE Optimal Control

On a State-Constrained PDE Optimal Control Problem arising from ODE-PDE Optimal Control On a Sae-Consrained PDE Opimal Conrol Problem arising rom ODE-PDE Opimal Conrol S. Wendl, H. J. Pesch and A. Rund Chair o Mahemaics in Engineering Sciences, Universiy o Bayreuh, Germany sean.wendl@uni-bayreuh.de,

More information

Figure 1. Jaw RMS target-tracker difference for a.9hz sinusoidal target.

Figure 1. Jaw RMS target-tracker difference for a.9hz sinusoidal target. Approaches o Bayesian Smooh Unimodal Regression George Woodworh Dec 3, 999 (Draf - please send commens o george-woodworh@uiowa.edu). Background Speech Ariculaion Daa The daa in Figure were obained by asking

More information

Mesh Refinement Using Density Function for Solving Optimal Control Problems

Mesh Refinement Using Density Function for Solving Optimal Control Problems Mesh Refinemen Using Densiy Funcion for Solving Opimal Conrol Problems Yiming Zhao and Panagiois Tsioras School of Aerospace Engineering, Georgia Insiue of Technology, Alana, GA, 333-5 This paper presens

More information

Scientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS

Scientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS Scieniic Research o he Insiue o Mahemaics and Compuer Science DIERENT VARIANTS O THE BOUNDARY ELEMENT METHOD OR PARABOLIC EQUATIONS Ewa Majchrzak,, Ewa Ładyga Jerzy Mendakiewicz, Alicja Piasecka Belkhaya

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

THE SINE INTEGRAL. x dt t

THE SINE INTEGRAL. x dt t THE SINE INTEGRAL As one learns in elemenary calculus, he limi of sin(/ as vanishes is uniy. Furhermore he funcion is even and has an infinie number of zeros locaed a ±n for n1,,3 Is plo looks like his-

More information

A Shooting Method for A Node Generation Algorithm

A Shooting Method for A Node Generation Algorithm A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan

More information

Initial Value Problems

Initial Value Problems Iniial Value Problems ChEn 2450 d d f(, ) (0) 0 6 ODE.key - November 26, 2014 Example - Cooking a Lobser Assumpions: The lobser remains a a uniform emperaure. This implies ha he hermal conduciviy of he

More information

Math 221: Mathematical Notation

Math 221: Mathematical Notation Mah 221: Mahemaical Noaion Purpose: One goal in any course is o properly use he language o ha subjec. These noaions summarize some o he major conceps and more diicul opics o he uni. Typing hem helps you

More information

Inverse Heat Conduction in a Finite Slab with Measured. Back Surface Temperature and Heat Flux

Inverse Heat Conduction in a Finite Slab with Measured. Back Surface Temperature and Heat Flux Inverse Hea Conducion in a Finie Slab wih Measured Bac Surace Temperaure and Hea Flux Jianhua Zhou, Yuwen Zhang, J. K. Chen 3, and Z. C. Feng 4 Universiy o Missouri, Columbia, MO 65, USA In High-Energy

More information

UNIVERSITY OF TRENTO MEASUREMENTS OF TRANSIENT PHENOMENA WITH DIGITAL OSCILLOSCOPES. Antonio Moschitta, Fabrizio Stefani, Dario Petri.

UNIVERSITY OF TRENTO MEASUREMENTS OF TRANSIENT PHENOMENA WITH DIGITAL OSCILLOSCOPES. Antonio Moschitta, Fabrizio Stefani, Dario Petri. UNIVERSIY OF RENO DEPARMEN OF INFORMAION AND COMMUNICAION ECHNOLOGY 385 Povo reno Ialy Via Sommarive 4 hp://www.di.unin.i MEASUREMENS OF RANSIEN PHENOMENA WIH DIGIAL OSCILLOSCOPES Anonio Moschia Fabrizio

More information

Kinematics and kinematic functions

Kinematics and kinematic functions Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables and vice versa Direc Posiion

More information

Application of He s Variational Iteration Method for Solving Seventh Order Sawada-Kotera Equations

Application of He s Variational Iteration Method for Solving Seventh Order Sawada-Kotera Equations Applied Mahemaical Sciences, Vol. 2, 28, no. 1, 471-477 Applicaion of He s Variaional Ieraion Mehod for Solving Sevenh Order Sawada-Koera Equaions Hossein Jafari a,1, Allahbakhsh Yazdani a, Javad Vahidi

More information

Topic Astable Circuits. Recall that an astable circuit has two unstable states;

Topic Astable Circuits. Recall that an astable circuit has two unstable states; Topic 2.2. Asable Circuis. Learning Objecives: A he end o his opic you will be able o; Recall ha an asable circui has wo unsable saes; Explain he operaion o a circui based on a Schmi inverer, and esimae

More information

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)

More information

Maxime Chupin 1,2, Thomas Haberkorn 3 and Emmanuel Trélat 1

Maxime Chupin 1,2, Thomas Haberkorn 3 and Emmanuel Trélat 1 ESAIM: MAN 51 (17) 965 996 DOI: 1.151/man/1644 ESAIM: Mahemaical Modelling and Numerical Analsis www.esaim-man.org LOW-THRUST LYAPUNOV TO LYAPUNOV AND HALO TO HALO MISSIONS WITH L -MINIMIZATION Maime Chupin

More information

Multiple Choice Solutions 1. E (2003 AB25) () xt t t t 2. A (2008 AB21/BC21) 3. B (2008 AB7) Using Fundamental Theorem of Calculus: 1

Multiple Choice Solutions 1. E (2003 AB25) () xt t t t 2. A (2008 AB21/BC21) 3. B (2008 AB7) Using Fundamental Theorem of Calculus: 1 Paricle Moion Soluions We have inenionally included more maerial han can be covered in mos Suden Sudy Sessions o accoun for groups ha are able o answer he quesions a a faser rae. Use your own judgmen,

More information

Sections 2.2 & 2.3 Limit of a Function and Limit Laws

Sections 2.2 & 2.3 Limit of a Function and Limit Laws Mah 80 www.imeodare.com Secions. &. Limi of a Funcion and Limi Laws In secion. we saw how is arise when we wan o find he angen o a curve or he velociy of an objec. Now we urn our aenion o is in general

More information

Optima and Equilibria for Traffic Flow on a Network

Optima and Equilibria for Traffic Flow on a Network Opima and Equilibria for Traffic Flow on a Nework Albero Bressan Deparmen of Mahemaics, Penn Sae Universiy bressan@mah.psu.edu Albero Bressan (Penn Sae) Opima and equilibria for raffic flow 1 / 1 A Traffic

More information

Interpretation of special relativity as applied to earth-centered locally inertial

Interpretation of special relativity as applied to earth-centered locally inertial Inerpreaion of special relaiviy as applied o earh-cenered locally inerial coordinae sysems in lobal osiioning Sysem saellie experimens Masanori Sao Honda Elecronics Co., Ld., Oyamazuka, Oiwa-cho, Toyohashi,

More information

Estimation of Poses with Particle Filters

Estimation of Poses with Particle Filters Esimaion of Poses wih Paricle Filers Dr.-Ing. Bernd Ludwig Chair for Arificial Inelligence Deparmen of Compuer Science Friedrich-Alexander-Universiä Erlangen-Nürnberg 12/05/2008 Dr.-Ing. Bernd Ludwig (FAU

More information

Lecture 10 Estimating Nonlinear Regression Models

Lecture 10 Estimating Nonlinear Regression Models Lecure 0 Esimaing Nonlinear Regression Models References: Greene, Economeric Analysis, Chaper 0 Consider he following regression model: y = f(x, β) + ε =,, x is kx for each, β is an rxconsan vecor, ε is

More information

ECE 2100 Circuit Analysis

ECE 2100 Circuit Analysis ECE 1 Circui Analysis Lesson 35 Chaper 8: Second Order Circuis Daniel M. Liynski, Ph.D. ECE 1 Circui Analysis Lesson 3-34 Chaper 7: Firs Order Circuis (Naural response RC & RL circuis, Singulariy funcions,

More information

CHAPTER 2: Mathematics for Microeconomics

CHAPTER 2: Mathematics for Microeconomics CHAPTER : Mahemaics for Microeconomics The problems in his chaper are primarily mahemaical. They are inended o give sudens some pracice wih he conceps inroduced in Chaper, bu he problems in hemselves offer

More information

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints IJCSI Inernaional Journal of Compuer Science Issues, Vol 9, Issue 1, No 1, January 2012 wwwijcsiorg 18 Applying Geneic Algorihms for Invenory Lo-Sizing Problem wih Supplier Selecion under Sorage Capaciy

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y Review - Quiz # 1 (1) Solving Special Tpes of Firs Order Equaions I. Separable Equaions (SE). d = f() g() Mehod of Soluion : 1 g() d = f() (The soluions ma be given implicil b he above formula. Remember,

More information

Simulating models with heterogeneous agents

Simulating models with heterogeneous agents Simulaing models wih heerogeneous agens Wouer J. Den Haan London School of Economics c by Wouer J. Den Haan Individual agen Subjec o employmen shocks (ε i, {0, 1}) Incomplee markes only way o save is hrough

More information

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4) Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion

More information

Simulation of magneto-hydrodynamic (MHD) flows in OpenFOAM

Simulation of magneto-hydrodynamic (MHD) flows in OpenFOAM Simulaion o magneo-hydrodynamic (MHD) los in OpenFOAM Sepember 1-15, 2010 NUMAP-FOAM, Universiy o Zagreb Ouline Problem: magneohydrodynamic (MHD) lo Magneo-hydrodynamic equaions Opimizaion o presen solver

More information

A Hop Constrained Min-Sum Arborescence with Outage Costs

A Hop Constrained Min-Sum Arborescence with Outage Costs A Hop Consrained Min-Sum Arborescence wih Ouage Coss Rakesh Kawara Minnesoa Sae Universiy, Mankao, MN 56001 Email: Kawara@mnsu.edu Absrac The hop consrained min-sum arborescence wih ouage coss problem

More information

A Discrete Variational Integrator for Optimal Control Problems on SO(3)

A Discrete Variational Integrator for Optimal Control Problems on SO(3) Proceedings of he 45h IEEE Conference on Decision & Conrol Mancheser Grand Hya Hoel San Diego, CA, USA, December 13-15, 26 A Discree Variaional Inegraor for Opimal Conrol Problems on SO3 Islam I. Hussein

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

SOME PROPERTIES OF GENERALIZED STRONGLY HARMONIC CONVEX FUNCTIONS MUHAMMAD ASLAM NOOR, KHALIDA INAYAT NOOR, SABAH IFTIKHAR AND FARHAT SAFDAR

SOME PROPERTIES OF GENERALIZED STRONGLY HARMONIC CONVEX FUNCTIONS MUHAMMAD ASLAM NOOR, KHALIDA INAYAT NOOR, SABAH IFTIKHAR AND FARHAT SAFDAR Inernaional Journal o Analysis and Applicaions Volume 16, Number 3 2018, 427-436 URL: hps://doi.org/10.28924/2291-8639 DOI: 10.28924/2291-8639-16-2018-427 SOME PROPERTIES OF GENERALIZED STRONGLY HARMONIC

More information

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation Commun Theor Phys Beijing, China 43 2005 pp 591 596 c Inernaional Academic Publishers Vol 43, No 4, April 15, 2005 An Invariance for 2+1-Eension of Burgers Equaion Formulae o Obain Soluions of KP Equaion

More information

( ) = Q 0. ( ) R = R dq. ( t) = I t

( ) = Q 0. ( ) R = R dq. ( t) = I t ircuis onceps The addiion of a simple capacior o a circui of resisors allows wo relaed phenomena o occur The observaion ha he ime-dependence of a complex waveform is alered by he circui is referred o as

More information

Exponential model. The Gibbs sampler is described for ecological examples by Clark et al. (2003, 2004) and

Exponential model. The Gibbs sampler is described for ecological examples by Clark et al. (2003, 2004) and APPEDIX: GIBBS SAMPLER FOR BAYESIA SAE-SPACE MODELS Eponenial model he Gis sampler is descried or ecological eamples Clark e al. (3 4) and Wikle (3). I involves alernael sampling rom each o he condiional

More information

Stability analysis of nonsmooth limit cycles with applications from power electronics

Stability analysis of nonsmooth limit cycles with applications from power electronics Inernaional Workshop on Resonance Oscillaions & Sabiliy of Nonsmooh Sysems Imperial College June 29 Sabiliy analysis of nonsmooh limi cycles wih applicaions from power elecronics Dr D. Giaouris damian.giaouris@ncl.ac.uk

More information

HERMITE POLYNOMIALS FOR ITERATIVE OUTPUT REPLANNING FOR FLAT SYSTEMS AFFECTED BY ADDITIVE NOISE

HERMITE POLYNOMIALS FOR ITERATIVE OUTPUT REPLANNING FOR FLAT SYSTEMS AFFECTED BY ADDITIVE NOISE Asian Journal o Conrol, Vol. 5, No., pp. 292 3, January 23 Published online 4 March 22 in Wiley Online ibrary (wileyonlinelibrary.com) DOI:.2/asjc.55 HERMITE POYNOMIAS FOR ITERATIVE OUTPUT REPANNING FOR

More information

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion

More information

Order Reduction of Large Scale DAE Models. J.D. Hedengren and T. F. Edgar Department of Chemical Engineering The University of Texas at Austin

Order Reduction of Large Scale DAE Models. J.D. Hedengren and T. F. Edgar Department of Chemical Engineering The University of Texas at Austin Order Reducion of Large Scale DAE Models J.D. Hedengren and T. F. Edgar Deparmen of Chemical Engineering The Universiy of Teas a Ausin 1 Ouline Moivaion Two Sep Process for DAE Model Reducion 1. Reducion

More information