Optimization of a point-mass walking model using direct collocation and sequential quadratic programming
|
|
- Opal Pitts
- 5 years ago
- Views:
Transcription
1 Optimization of a point-mass walking moel using irect collocation an sequential quaratic programming Chris Dembia June 5, 5
2 Telescoping actuator y Stance leg Point-mass boy m (x,y) Swing leg x Leg uring stance F(t) F(t) l(t) F(t) mg Boy uring stance mg Boy uring flight Figure : The point mass moel of walking from [SR6]. a) Some possible gaits b) Inverte penulum walk c) Impulsive run ) Hybri intermeiate gait: penular run Stance Flight Inverte penulum Bounce Inverte penulum Flight Flight l max Bounce lmax Figure : (a) Possible perioic gaits for the point mass moel. (b-) Optimal gaits for the point mass moel, epening on the prescribe spee an step step length of the gait [SR6]. Introuction It is possible to uncover some funamental principles of human locomotion with a moel that consists of only a point mass an a massless telescoping (compressible) leg. Srinivasan an Ruina foun that such a moel chooses to walk if aske to move at slow spees, an to run if aske to move at high spees [SR6]. The moel they use is shown in Figure. The moel interacts with its environment via a telescoping linear actuator (its leg). The motion consists of two phases: a stance phase, uring which the stance leg is in contact with the groun; an a swing phase uring which the point mass is in free flight. Specifically, their moel chose between three istinct gaits, epening on the spee an step length of the gait: an inverte penulum walk, an impulsive run, an a penular run; these are shown in Figure. Srinivasan an Ruina foun these gaits by solving an optimal control problem in which the objective was to maximize the energy efficiency of the gait. The parameters in the optimization consiste of the initial state of the moel an the force applie by the actuator throughout the stance phase. To evaluate the objective function, they performe a forwar simulation (integration) of the moel s equations of motion, starting from the initial state. Their metho is calle single shooting; the optimizer chooses an initial state an shoots to the final state by integrating. This metho is typically challenging for optimizers, since the motion of a nonlinear ynamical system can be very sensitive to the initial state. Another
3 approach is to iscretize the state of the moel through time, an to enforce the ynamics at these time points through explicit constraints in the optimization problem. This metho is calle irect collocation [Bet]. When formulate via irect collocation, optimal control problems are far less sensitive to variations in the parameters. Direct collocation problems have many more parameters than single shooting problems, but are extremely sparse. In this work, we replicate the work of Srinivasan an Ruina using irect collocation instea of single shooting. We also solve the problem with a suite of ifferent cost functions, an explore the effect of these ifferent costs on the resulting motion. The problem Here, we escribe the optimal control problem that Srinivasan an Ruina [SR6] solve. The variables consist of the moel s state s(t) = (x(t), y(t), ẋ(t), ẏ(t)), where (x(t), y(t)) is the position of the point mass; the force applie to the groun by the leg F (t); an the uration of the stance phase t s. The step length, spee v, an step uration t step = /v are specifie. The leg length l(t) = x (t) + y (t) has a maximum value l max uring the stance phase, an the leg also has a finite strength. The motion must also be perioic. Formally, the optimal control problem is minimize C = tstep E(s(t), F (t)) t mg cost subject to mẍ = F x/l t t s stance ynamics mÿ = F y/l mg t t s stance ynamics mẍ = t s < t t step swing ynamics mÿ = mg t s < t t step swing ynamics l(t) l max t t s leg length F (t) t t s leg strength < t s t step stance uration x(t step ) x() = istance travele y() y(t step ) ẋ() = ẋ(t step ). perioicity ẏ() ẏ(t step ) () The point mass has mass m, an g is the acceleration ue to Earth s gravity on Earth s surface. For the integran cost E, Srinivasan an Ruina use the positive work performe by the leg, E = [F (t) l(t)] + ; with this E, then C is the cost of transport, which is a imensionless measure of the energetic cost of traveling a certain istance. In their work, the ynamics constraints were enforce via forwar integration. In this work, we use m = 5kg, g = 9.8m/s, l max =.5m, = N. Different gaits are foun by varying the nonimensional step length D = /l max an spee V = v/ gl max. 3
4 . Methos Direct collocation We iscretize Problem into T step mesh points t, t,..., t Tstep = t step. The transition from stance to swing occurs at time inex i = T s ; that is, t s = t Ts T step. In orer to express the secon orer ynamics of the system using first orer ifferential equations, we introuce the variables w = ẋ an z = ẏ for the horizontal an vertical spee of the point mass. The iscretization yiels state variables s, s,..., s Tstep (with s i = (x i, y i, w i, z i )) an leg forces F, F,..., F Ts. We enforce the ifferential constraints using a numerical metho for solving orinary ifferential equations, namely the implicit Euler metho: for a ifferential equation ẏ = f(y, t), we obtain y via y i+ = y i +h i f(y i+, t i+ ). We use an implicit integrator because they are unconitionally stable (the solution won t blow up ). The resulting finite-imensional optimization problem is minimize C = T h i [E(s i+, F i+ ) + E(s i, F i )] i= subject to x i+ x i = h i w i+ i =,,..., T step y i+ yi = h i z i+ i =,,..., T step x i+ F i+ w i+ w i = h i i =,,..., T s ( m l i+ ) Fi+ y i+ z i+ z i = h i g i =,,..., T s m l i+ w i+ w i = i = T s,..., T step z i+ z i = h i g i = T s,..., T step li lmax i =,,..., T s F i i =,,..., T s < t s t step x Tstep x = y y Tstep w = w Tstep. z z Tstep () The scalars h i = t i+ t i are the integration step sizes. The performance cost is now approximate by the iscretize C, which is compute using the trapezoial rule. This problem has a total n = 4T step + T s + variables (4 states, the leg force uring stance, an t s ) an 4(T step ) ynamics constraints. For all results shown in this report, we use T s =, T step =. The value of T step is chosen to be sufficiently large so that the error in the ODE solution is below some tolerance, typically 3. For all of the constraints, the parameters at only one or two of the mesh points is present; therefore, this problem is very sparse. Since E may be non-convex, an since the thir an fourth constraints in Problem are non-affine, this problem is not convex. 4
5 Sequential quaratic programming We solve Problem using sequential quaratic programming [GMS94]. At each iteration k of the metho, we approximate the problem as a quaratic program about our current iterate q (k). The solution of this QP, q (k), is use to compute a search irection q (k) = q (k) q (k) along which to perform a line search. The result of this line search is the next iterate q (k+). Within the QP, the quaratic objective is q C (k) (q q(k) ) + (q q(k) ) qc (k) (q q(k) ), where q C (k) an qc (k) are the graient an Hessian evaluate at q (k), respectively. The only constraints that must be approximate for the QP are the thir an fourth from Problem. We efine the following: x i+ F i+ a i = w i+ w i h i, ( m l i+ Fi+ y i+ = z i+ z i h i b i m u i = (w i+, w i, x i+, y i+, F i+ ), v i = (z i+, z i, x i+, y i+, F i+ ). ) g, l i+ The quantities a i an b i are the resiuals in thir an fourth constraints, an u i an v i are the parameters containe in the thir an fourth constraints. The affine approximations to these constraints are a (k) i + ( ui a (k) i ) T (u i u (k) i ) = b (k) i + ( vi b (k) i ) T (v i v (k) i ) =. The graients in these equations are [ ( ( ui a i ) T F i+ x = h i+ i ) F i+ x i+ y i+ h m li+ 3 i [ l i+ ( m li+ 3 ( vi b i ) T F i+ x i+ y i+ F i+ y = h i h i+ m li+ 3 i ) m li+ 3 l i+ h i m h i m x i+ l i+ To solve the SQP, we use the SNOPT package [GMS5]. SNOPT is a package for solving large sparse non-convex problems, an approximates the Hessian using a BFGS upate. For the problems we explore, SNOPT requires between 5 to 35 iterations to converge; on a moern personal computer, this takes about minute. During the first few outer iterations, the QP solver takes potentially hunres of steps. After about 5 outer iterations, when the iterate is nearly feasible, the QP usually only requires two steps. Transcribing the problem The proceure of converting Problem into Problem is calle transcription. We use the PSOPT optimal control package to perform this transcription for us [Bec]. PSOPT also iteratively solves the problem on larger meshes until an ODE error tolerance is satisfie. 5 y i+ l i+ ] ].
6 b) Penular walk c) Impulsive run ) Penular run Leg-force F Leg-length L Inverte penulum e e e e τ step τ τ s τ step τ τ step τ Bounce Inverte Flight penulum Flight F' max F' max F' max τ τ τ τ s 3 Results Figure 3: Leg length an leg force for the three gaits prouce by [SR6]. Reproucing the results of Srinivasan an Ruina, 6 By varying the values of the non-imensional step length D an non-imensional spee V, Srinivasan an Ruina iscovere three istinct gaits (Figure ). Figure 3 shows the leg force F (t) an leg length l(t) for a specific instance of each of these three gaits. The penular walk emerges from small D an small V, the impulsive run emerges for all D (> ) with V above, an the hybri penular run emerges from D above for V <. We reprouce these results qualitatively using the methos presente in section.; the three gaits are shown in Figure 4. The force F is not smooth as a result of the non-smooth cost E. Notably, the penular walk oes not contain a swing phase, an the force for both penular gaits is U-shape. Exploring ifferent cost functions We explore three ifferent cost functions an analyze their effect on the resulting walking solution, similar to work one by [Sri]. For all three, we use D =.5 an V =.5. Minimize force over the gait cycle. With E = F (t), the solution minimizes the amount of force use throughout the whole motion. The solution shows that this causes the introuction of a swing phase that lasts as long as the stance phase. Minimize square leg force. The cost E = F (t) penalizes large forces. This objective might capture the esire to reuce joint loaing for those with joint pain. With this gait, the moel oes not change its height at all, an the leg force remains close to only what is necessary to counteract gravity; F (t) mg. Minimize gross work. The cost E = F (t) l(t) essentially allows the moel to propel itself using the energy absorbe from breaking. Inee, with this cost function, the net work over the gait is. The resulting motion is spring-like. 6
7 E =[F L] +, D =.5, V =.5 E =[F L] +, D =.5, V =.5 E =[F L] +, D =.5, V =.5.8 t step.8 t s t step.8 t step t step t s t step t step Figure 4: Three istinct gaits (left to right): penular walk, impulsive run, penular run. E =F, D =.5, V =.5 E =F, D =.5, V =.5 E =F L, D =.5, V =.5.8 t step.8 t step.6 t step t step t step t step Figure 5: Walking with ifferent costs (left to right): leg force, square leg force, gross work. 7
8 References [Bec] [Bet] V. M. Becerra. Solving complex optimal control problems at no cost with PSOPT. Proc. IEEE Multi-conference on Systems an Control, Yokohama, Japan, pages ,. John Betts. Practical Methos for Optimal Control an Estimation using Nonlinear Programming. Society for Inustrial an Applie Mathematics, secon eition,. [GMS94] Philip E Gill, Walter Murray, an Michael A Sauners. Large-scale sqp methos an their application in trajectory optimization. In Computational optimal control, pages 9 4. Springer, 994. [GMS5] Philip E. Gill, Walter Murray, an Michael A. Sauners. Snopt: An SQP Algorithm For Large-Scale Constraine Optimization. SIAM Review, 47():99 3, 5. [SR6] [Sri] Manoj Srinivasan an Any Ruina. Computer optimization of a minimal bipe moel iscovers walking an running. Nature, 439(77):7 75, 6. Manoj Srinivasan. Fifteen observations on the structure of energy-minimizing gaits in many simple bipe moels. Journal of The Royal Society Interface, 8(54):74 98,. 8
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 informationChapter 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 informationNumerical Integrator. Graphics
1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,
More informationu t v t v t c a u t b a v t u t v t b a
Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying
More informationLecture 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 informationTable 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 informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationMath 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 informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More information6 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 informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationMechatronics by Analogy and Application to Legged Locomotion
Mechatronics by Analogy an Application to Legge Locomotion by Victor Ragusila A thesis submitte in conformity with the requirements for the egree of Doctor of Philosophy Grauate Department of Aerospace
More informationOrdinary Differential Equations: Homework 1
Orinary Differential Equations: Homework 1 M. Gameiro, J.-P. Lessar, J.D. Mireles James, K. Mischaikow January 12, 2017 2 Chapter 1 Motivation 1.1 Exercises Exercise 1.1.1. (Frictionless spring) Consier
More informationText S1: Simulation models and detailed method for early warning signal calculation
1 Text S1: Simulation moels an etaile metho for early warning signal calculation Steven J. Lae, Thilo Gross Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, 01187 Dresen, Germany
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationMath 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v
Math Fall 06 Section Monay, September 9, 06 First, some important points from the last class: Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v passing through
More informationCalculus 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 informationThe 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 informationAssignment 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 informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationECE 422 Power System Operations & Planning 7 Transient Stability
ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of
More informationAdjoint Transient Sensitivity Analysis in Circuit Simulation
Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:
More informationLagrangian and Hamiltonian Dynamics
Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion
More informationTopic 2.3: The Geometry of Derivatives of Vector Functions
BSU Math 275 Notes Topic 2.3: The Geometry of Derivatives of Vector Functions Textbook Sections: 13.2 From the Toolbox (what you nee from previous classes): Be able to compute erivatives scalar-value functions
More informationSemi passive walking of a 7 DOF biped robot ABSTRACT 1 INTRODUCTION
Semi passive walking of a 7 DOF bipe robot N.Khraief & N.K.M Siri khraief@lrv.uvsq.fr, nkms@free.fr Laboratoire e Robotique e Versailles Université e Versailles Saint Quentin en Yvelines 0, Avenue e l
More informationMinimum-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 informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
More informationBEYOND 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 informationwhere: q = [q 1 ] T : vector of generalized coordinates and q and q are vectors of joint angular velocity and acceleration.
J Mech Behav Mater 15; 4(1-): 53 57 Maime Sare* Control of a legge robot Abstract: This paper eals with the control of hopping an running systems that interact intermittently with the environment. The
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More information19 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 information2.25 m. (a) Using Newton s laws of motion, explain why the student can gain an initial speed to leave the ground vertically.
NAME : F.5 ( ) MARS: /70 FORM FIVE PHYSICS TEST on MECHANICS Time Allowe: 70 minutes This test consists of two sections: Section A (structure type questions, 50 marks); Section B (multiple choice, 20 marks)
More informationChapter 9 Method of Weighted Residuals
Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe
More informationStatics. There are four fundamental quantities which occur in mechanics:
Statics Mechanics isabranchofphysicsinwhichwestuythestate of rest or motion of boies subject to the action of forces. It can be ivie into two logical parts: statics, where we investigate the equilibrium
More informationORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa
Manuscript submitte to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS Gianluca Crippa Departement Mathematik
More informationthere 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 informationIntroduction 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 informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationThe dynamics of the simple pendulum
.,, 9 G. Voyatzis, ept. of Physics, University of hessaloniki he ynamics of the simple penulum Analytic methos of Mechanics + Computations with Mathematica Outline. he mathematical escription of the moel.
More informationEVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES
MATHEMATICS OF COMPUTATION Volume 69, Number 231, Pages 1117 1130 S 0025-5718(00)01120-0 Article electronically publishe on February 17, 2000 EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION
More informationMechanical energy transfer by internal force during the swing phase of running
Available online at www.sciencedirect.com Procedia Engineering 34 (2012 ) 772 777 9 th Conference of the International Sports Engineering Association (ISEA) Mechanical energy transfer by internal force
More informationExperimental 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'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 informationarxiv: v1 [math.oc] 5 Oct 2016
CONTROLLING A SELF-ORGANIZING SYSTEM OF INDIVIDUALS GUIDED BY A FEW EXTERNAL AGENTS PARTICLE DESCRIPTION AND MEAN-FIELD LIMIT MARTIN BURGER, RENÉ PINNAU, ANDREAS ROTH, CLAUDIA TOTZECK, AND OLIVER TSE arxiv:161.1325v1
More informationDynamical Systems and a Brief Introduction to Ergodic Theory
Dynamical Systems an a Brief Introuction to Ergoic Theory Leo Baran Spring 2014 Abstract This paper explores ynamical systems of ifferent types an orers, culminating in an examination of the properties
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationFrom 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 informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationA Course in Machine Learning
A Course in Machine Learning Hal Daumé III 12 EFFICIENT LEARNING So far, our focus has been on moels of learning an basic algorithms for those moels. We have not place much emphasis on how to learn quickly.
More informationGoldstein Chapter 1 Exercises
Golstein Chapter 1 Exercises Michael Goo July 17, 2004 1 Exercises 11. Consier a uniform thin isk that rolls without slipping on a horizontal plane. A horizontal force is applie to the center of the isk
More informationApproximate 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 informationA 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 informationmodel 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 informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationReal-time economic optimization for a fermentation process using Model Predictive Control
Downloae from orbit.tu.k on: Nov 5, 218 Real-time economic optimization for a fermentation process using Moel Preictive Control Petersen, Lars Norbert; Jørgensen, John Bagterp Publishe in: Proceeings of
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationRole of parameters in the stochastic dynamics of a stick-slip oscillator
Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society
More informationNonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain
Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationWUCHEN 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 information2-7. Fitting a Model to Data I. A Model of Direct Variation. Lesson. Mental Math
Lesson 2-7 Fitting a Moel to Data I BIG IDEA If you etermine from a particular set of ata that y varies irectly or inversely as, you can graph the ata to see what relationship is reasonable. Using that
More informationAPPPHYS 217 Thursday 8 April 2010
APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q
More informationACCELERATION, FORCE, MOMENTUM, ENERGY : solutions to higher level questions
ACCELERATION, FORCE, MOMENTUM, ENERGY : solutions to higher level questions 015 Question 1 (a) (i) State Newton s secon law of motion. Force is proportional to rate of change of momentum (ii) What is the
More informationThe Hamiltonian particle-mesh method for the spherical shallow water equations
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 5: 89 95 (004) Publishe online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.100/asl.70 The Hamiltonian particle-mesh metho for the spherical
More informationNew Simple Controller Tuning Rules for Integrating and Stable or Unstable First Order plus Dead-Time Processes
Proceeings of the 3th WSEAS nternational Conference on SYSTEMS New Simple Controller Tuning Rules for ntegrating an Stable or Unstable First Orer plus Dea-Time Processes.G.ARVANTS Department of Natural
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationMore from Lesson 6 The Limit Definition of the Derivative and Rules for Finding Derivatives.
Math 1314 ONLINE More from Lesson 6 The Limit Definition of the Derivative an Rules for Fining Derivatives Eample 4: Use the Four-Step Process for fining the erivative of the function Then fin f (1) f(
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationMARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ
GENERALIZED UNIFORMLY CONTINUOUS SEMIGROUPS AND SEMILINEAR HYPERBOLIC SYSTEMS WITH REGULARIZED DERIVATIVES MARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ Abstract. We aopt the theory of uniformly continuous operator
More informationTRAJECTORY 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 informationReactive Power Compensation in Mechanical Systems
Reactive Power Compensation in Mechanical Systems Carlos Rengifo, Bassel Kaar, Yannick Aoustin, Christine Chevallereau To cite this version: Carlos Rengifo, Bassel Kaar, Yannick Aoustin, Christine Chevallereau.
More informationExponential 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 informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationProof of SPNs as Mixture of Trees
A Proof of SPNs as Mixture of Trees Theorem 1. If T is an inuce SPN from a complete an ecomposable SPN S, then T is a tree that is complete an ecomposable. Proof. Argue by contraiction that T is not a
More informationensembles 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 informationFurther Differentiation and Applications
Avance Higher Notes (Unit ) Prerequisites: Inverse function property; prouct, quotient an chain rules; inflexion points. Maths Applications: Concavity; ifferentiability. Real-Worl Applications: Particle
More informationCalculus I Practice Test Problems for Chapter 3 Page 1 of 9
Calculus I Practice Test Problems for Chapter 3 Page of 9 This is a set of practice test problems for Chapter 3. This is in no wa an inclusive set of problems there can be other tpes of problems on the
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More information12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes
Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow
More informationOptimal 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 informationThe Principle of Least Action and Designing Fiber Optics
University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic
More informationMonopoly Part III: Multiple Local Equilibria and Adaptive Search
FH-Kiel University of Applie Sciences Prof Dr Anreas Thiemer, 00 e-mail: anreasthiemer@fh-kiele Monopoly Part III: Multiple Local Equilibria an Aaptive Search Summary: The information about market eman
More informationStabilization of a Class of Underactuated Mechanical Systems via Interconnection and Damping Assignment
Stabilization of a Class of Uneractuate Mechanical Systems via Interconnection an Damping Assignment Romeo Ortega Lab. es Signaux et Systèmes CNRS-SUPELEC Gif sur Yvette 99 FRANCE tel. no:(33)--69-85-7-65
More informationTime Optimal Flight of a Hexacopter
Time Optimal Flight of a Hexacopter Johan Vertens an Fabian Fischer Abstract This project presents a way to calculate a time optimal control path of a hexacopter using irect multiple shooting in Matlab
More informationLQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR
5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,
More information221A Lecture Notes Notes on Classica Mechanics I
1A Lecture Notes Notes on Classica Mechanics I 1 Precursor: Fermat s Principle in Geometric Optics In geometric optics, you talk about how light rays go. In homogeneous meiums, the light rays go straight.
More informationx f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.
Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationDot trajectories in the superposition of random screens: analysis and synthesis
1472 J. Opt. Soc. Am. A/ Vol. 21, No. 8/ August 2004 Isaac Amiror Dot trajectories in the superposition of ranom screens: analysis an synthesis Isaac Amiror Laboratoire e Systèmes Périphériques, Ecole
More informationPredictive Control of a Laboratory Time Delay Process Experiment
Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results
More informationAn Introduction to Event-triggered and Self-triggered Control
An Introuction to Event-triggere an Self-triggere Control W.P.M.H. Heemels K.H. Johansson P. Tabuaa Abstract Recent evelopments in computer an communication technologies have le to a new type of large-scale
More informationExplicit Formulas for General Integrals of Motion for a Class of Mechanical Systems Subject To Virtual Constraints
Explicit Formulas for General Interals of Motion for a Class of Mechanical Systems Subject To Virtual Constraints Shiriaev, Anton; Perram, John; Robertsson, Aners; Sanber, Aners Publishe in: Decision an
More information