Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications
|
|
- Francis Cannon
- 5 years ago
- Views:
Transcription
1 Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications Navin Khaneja lectre notes taken by Christiane Koch Jne 24, 29 1 Variation yields a classical Hamiltonian system Sppose that we have a system that is described by the following eqations of motion, ẋ = f(x,, t) x R n, Ω R m, (1) where x denotes the state vector, the controls, and t time. We can inclde t in the state vector with ṫ = 1. Eq. (1) is then redced to ẋ = f(x, ) x R n+1, Ω R m. (2) The cost fnction for the optimization problem is typically written as a sm of two terms, J = φ(x f ) + L(x, )dt, (3) where φ(x f ) denotes the teral cost, with x f = x(t ) the state vector at final time T. The second term in Eq. (3) corresponds to the rnning cost. Adding zero, J can be rewritten J = φ(x f ) + = φ(x f ) + L(x, )dt + { L(x, ) + λ T f(x, ) + ( λ T )x The second line is obtained by integrating - λt ẋ dt by parts. In order to search for the optimm, we vary J, δj = φ δx f + λ T ( f(x, ) ẋ ) dt } dt λ T (T )x(t ) + λ T ()x(). { ] } L L f f δx(t) + δ(t) + λt δx(t) + δ(t) + λ T δx(t) dt λ T (T )δx(t ), (4) 1
2 where we have assmed that the initial state is fixed, δx() =, and δx(t ) = δx f. Regroping terms, we obtain ] φ δj = λ T δx f (5) {( ) ( ) } L f L f + + λt δx(t) + + λt δ(t) + ( λ T )δx(t) dt We are looking for an optimm of J, i.e. the variations shold vanish. The variation with respect to x(t) and x f, respectively, detere the eqation of motion for the adjoint state, λ T (t), and its initial condition: The term inside the sqare brackets in Eq. (5) is zero if ( φ T λ(t ) = ) and the terms inside the integral in Eq. (5) that are mltiplied with δx(t) vanish if ( ) L f T λ = + λt. We are then left with δj = We can introdce a (Hamiltonian) fnction We then see that and the variation of J becomes x f, ( ) L f + λt δ(t) dt. H(x, λ, ) = L(x, ) + λ T f(x, ). (6) δj = A small variation of can be written as λ = δ = ε, H δ(t) dt.. If is optimal, then at least H = sch that δj =. We then obtain the following set of 1st order necessary conditions, ẋ = f(x, ) = with x() (7) H λ λ = with λ(t ) = φ (8) = (9) (x(t),λ(t),(t)) 2
3 This set of eqations looks like a classical Hamiltonian system, where the state vector x takes the role of the coordinates and the adjoint state λ that of the conjgate momenta. 2 The maximm principle Let s first state Pontryagin s maximm principle 1] before proving it: Along the optimal trajectory, we have (t) = arg H ( x (t), λ (t), ). (1) Note that in or context it is rather a imm than a maximm principle since we seek to imize the cost. The claim is that if (t) is an optimal control, then it imizes H along the optimal trajectory globally at each instant of time. This is a very strong statement since we have only sed first order variation. To prove this 1, let s first recast J sch that it represents a final-time cost only. This can be done by defining, in Eq. (3), L(x, ) in terms of a state variable. That is, if or state vector is x = (x 1, x 2,..., x n ) R n, we add x n+1 with ẋ n+1 = L(x, ). (11) We can then write a cost fnctional that depends on the final time only, J = φ(x f ) + x n+1 (T ). (12) Now sppose that we know the optimal control (t) which deteres the optimal x (t). Dring a very short time interval, τ dτ, τ], let s change the control drastically from to v: Ω v v (t) τ dτ τ t x(t) 1 1 x(τ) x f t This leads to the following change in x, δx(τ) = f ( x (τ), v ) f ( x (τ), (τ) )] dτ. (13) We ths obtain a change of initial condition at time t = τ for the evoltion of x(t) at times t > τ. We can write (x + δx) = f(x + δx, + δ), 1 Note that we will show the proof only for reglar extrema. Pontryagin s proof covers also abnormal extrema 1] bt is mch more difficlt to follow. Roghly speaking, abnormal extrema are those where the final state x f can be reached only for a single T. 3
4 which leads to δx = f δx + f (x (t), (t)) δ (x }{{} (t), (t)) }{{} A(t) B(t) δx = A(t)δx(t) + B(t)δ(t). (14) This can be expressed in terms of a propagator (or Green s fnction), Φ, Or claim is now or, eqivalently, δx f = δx(t ) = Φ(T, τ) δx(τ). φ δx f, φ Φ(T, τ) δx(τ). (15) Using the initial condition, λ(t ), cf. Eq. (8), we can write for λ at time τ, T φ λ(τ) = Φ(T, τ)] = ( Φ(T, τ) ) T λ(t ) (16) Inserting Eq. (16) and Eq. (13) into Eq. (15), we obtain λ T (τ) f ( x (τ), v ) f ( x (τ), (t) )] δτ or, eqivalently, λ T (τ) f ( x (τ), (t) ) λ T (τ) f ( x (τ), v ) Since Eq. (6) for a final-time only cost becomes we obtain for all v Ω H(x, λ, ) = λ T f(x, ), (17) H(x, λ, ) H(x, λ, v), (18) which proves that indeed (t) imizes H globally along the optimal trajectory, i.e. Eq. (1). This necessary condition is stronger than jst the first order variation being zero. In fact, it is so restrictive that often it allows to detere the optimal soltion. In conclsion, if (t) is the global optimm, then there is a canonical way to define H with the necessary condition, Eq. (18) and generate the state x(t) and co-state λ(t). Eq. (17) implies H f = λt, which, together with the definition of A(t), cf. Eq. (14), leads to λ = = A(t) λ. (19) 4
5 An example: Controlling the force on a particle 1] The eqation of motion is given by which can be rewritten ẍ = with 1 ẋ 1 = x 2, ẋ 2 =. The goal is to drive the system from any initial point in phase space, (x 1 (), x 2 ()), to the origin in imal time. We therefore have to formlate time as a cost, ṫ = ẋ 3 = 1. Evalating Eq. (17), the eqations of motion lead to the Hamiltonian fnction for this optimization problem, H = λ 1 x 2 + λ 2 + λ 3. The λ j are fond from the Hamiltonian eqations, cf. Eq. (8), λ 1 = H 1 = λ 1 = c = const λ 2 = λ 1 λ 2 = c t + d λ 3 = We can now detere the control from Eq. (9), or, more specifically from the condition to maximize H. Taking the restriction, 1, into accont, H takes its maximm vale for = sgn(λ 2 ). This corresponds to sing the maximm allowed force, i.e. bang-bang control. Since λ 2 is a linear fnction, it changes sign only once. We therefore obtain for possibilities for : t t t t 1 Evalating the Hamiltonian fnction, Eq. (7), we can also obtain the evoltion of the state variables. For = 1, ẋ 2 =, yields x 2 (t) = t + c and for = 1, we get correspondingly x 1 (t) = t2 2 + c t + d = x d 1, x 1 (t) = x d 1. 5
6 That is, we obtain parabolas in phase space: A(t) x 2 I II B(t) x 1 If the controls are not switched, or initial point mst be on either crve A(t) or B(t), cf. Eq. (14), since these are the only parabolas that take s to the origin. If the initial point is (x 1 (), x 2 ()) I, then the control is switched from -1 to 1 when the parabola that contains (x 1 (), x 2 ()) and describes the flow (x 1 (t), x 2 (t) crosses B(t). The phase space flow then contines on B(t) to the origin. If (x 1 (), x 2 ()) II, then the control is switched from + 1 to 1 when the parabola that contains (x 1 (), x 2 ()) and describes (x 1 (t), x 2 (t) crosses A(t). The phase space flow then contines on A(t) to the origin. We have solved this example by constrcting the Hamiltonian fnction. Alternatively, it cold also be solved by considering the final-time cost, φ(x f ), which in the example above corresponds to i.e. the total time to reach the target. 2 φ(x 1,f, x 2,f, x 3,f ) = x 3,f = T, 3 Hamilton-Jacobi-Bellman principle When we transform all rnning costs into final costs φ(x f ), the optimization problem consists in designing the trajectory from initial point x to final point x f that imizes the final cost, φ(x f ). We seek to solve this problem globally for all initial points x and introdce the optimal retrn fnction, V (x ) = φ(x f ), (2) that is V (x) corresponds to the best or imm cost starting from initial vale x. Sppose that we have solved the optimization problem and know V (x). This defines the optimal trajectory from x to x f : 2 Note that φ(x f ) doesn t inclde the target state becase we assme the target state to be fixed and the zero variations δx 1,f, δx 2,f wold leave λ 1 and λ 2 ndefined. 6
7 Chicago 11 Boston δx() x x f 11 Santa Barbara Consider now the optimal retrn for small deviations from the optimal trajectory, The claim is now that V ( x + δx() ). V ( x + δx() ) = V (x). (21) This is the principle of dynamic programg or Hamilton-Jacobi-Bellman (HJB) principle. It states that if we have fond the optimal way to go from Santa Barbara to Boston, and Chicago is located on this way, then also the way from Chicago to Boston is optimal. We can rewrite Eq. (21) by expressing δx in terms of the eqation of motion of x, Eq. (1), which leads to the Bellman PDE, V (x) + V ] f(x, ) δt, = V (x) ] V f(x, ) =. (22) The optimal cost mst satisfy Eq. (22). If we want to achieve the optimization in inm time, we need to consider t as a state variable sch that or eqations of motion become Eq. (22) is then rewritten ẋ = f(x,, t), ṫ = 1 (23) V t + V ] f(x, ) =, (24) which is the Hamilton-Jacobi-Bellman eqation. Note that if there is no restriction on time, then V does not depend on t, i.e. the optimm is achieved at infinite time (these problems are called infinite horizon problems). We now show the connection between the HJB eqation, Eq. (24), and the Pontryagin principle, Eq. (1), by considering what happens along the optimal trajectory. Evalating, along the optimal trajectory, the partial derivative of V (x) with respect to x which also depends on x, we find ( ) V x (t) = λ T (t). (25) Inserting this into Eq. (22), we see that λ T (t)f ( x (t), )] =. 7
8 Since the term inside the sqare brackets is nothing bt the Hamiltonian, cf. Eq. (17), along the optimal trajectory, x (t), we see that the optimal control imizes H at each time. That is, we recover Pontryagin s imm principle. We can therefore detere the optimal control as the argment that imizes the Hamiltonian, or, (x) = arg (x) = arg H ( x (t), λ(t), )]. (26) ] V f(x, ). (27) Eq. (27) implies V f(x, (x)) =. (28) Eq. (26) implies that the Hamiltonian is zero along the optimal trajectory, why? H ( x (t), λ(t), (t) ) =, and that the optimal cost V is constant along the optimal trajectory, dv dt =, (x (t)) since dv dt = H. We now show that Eq. (25) indeed satisfies the eqation of motion for λ T (t), Eq. (8). Let s denote partial derivatives by sbscripts, Taking the derivative of Eq. (25), we obtain V = V x. λ T = V xx ( x (t) ) f(x, ). On the other hand, we can differentiate Eq. (28) with respect to x and find = ] V x f(x, ) = V xx f ( x, (x) ) + V x f x + V x f. For optimal, f = and we can express V xx in terms of V x f x and obtain λ T = V f. With Eq. (25), we find λ(t) = f T λ. 8
9 From the constrction of the Hamiltonian, Eq. (17), it follows that H = λt f, and we indeed recover Eq. (8). The connection to Pontryagin s maximm principle is easily made by considering the optimal retrn at the final point, V (x f ) = φ(x f ). This is inserted into Eq. (25) for t = T and we obtain References λ T (T ) = φ ( ). 1] L. S. Pontryagin. The Mathematical Theory of Optimal Processes. Gordon and Breach Science Pblishers,
Essentials of optimal control theory in ECON 4140
Essentials of optimal control theory in ECON 4140 Things yo need to know (and a detail yo need not care abot). A few words abot dynamic optimization in general. Dynamic optimization can be thoght of as
More informationThe Linear Quadratic Regulator
10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationPartial Differential Equations with Applications
Universit of Leeds MATH 33 Partial Differential Eqations with Applications Eamples to spplement Chapter on First Order PDEs Eample (Simple linear eqation, k + = 0, (, 0) = ϕ(), k a constant.) The characteristic
More informationMean Value Formulae for Laplace and Heat Equation
Mean Vale Formlae for Laplace and Heat Eqation Abhinav Parihar December 7, 03 Abstract Here I discss a method to constrct the mean vale theorem for the heat eqation. To constrct sch a formla ab initio,
More informationLecture: Corporate Income Tax
Lectre: Corporate Income Tax Ltz Krschwitz & Andreas Löffler Disconted Cash Flow, Section 2.1, Otline 2.1 Unlevered firms Similar companies Notation 2.1.1 Valation eqation 2.1.2 Weak atoregressive cash
More informationLecture: Corporate Income Tax - Unlevered firms
Lectre: Corporate Income Tax - Unlevered firms Ltz Krschwitz & Andreas Löffler Disconted Cash Flow, Section 2.1, Otline 2.1 Unlevered firms Similar companies Notation 2.1.1 Valation eqation 2.1.2 Weak
More informationMath 116 First Midterm October 14, 2009
Math 116 First Midterm October 14, 9 Name: EXAM SOLUTIONS Instrctor: Section: 1. Do not open this exam ntil yo are told to do so.. This exam has 1 pages inclding this cover. There are 9 problems. Note
More informationDynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining
Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining P. Hgenin*, B. Srinivasan*, F. Altpeter**, R. Longchamp* * Laboratoire d Atomatiqe,
More informationMath 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:
Math 263 Assignment #3 Soltions 1. A fnction z f(x, ) is called harmonic if it satisfies Laplace s eqation: 2 + 2 z 2 0 Determine whether or not the following are harmonic. (a) z x 2 + 2. We se the one-variable
More informationSection 7.4: Integration of Rational Functions by Partial Fractions
Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to
More informationInformation Theoretic views of Path Integral Control.
Information Theoretic views of Path Integral Control. Evangelos A. Theodoro and Emannel Todorov Abstract We derive the connections of Path IntegralPI and Klback-LieblerKL control as presented in machine
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More information4.2 First-Order Logic
64 First-Order Logic and Type Theory The problem can be seen in the two qestionable rles In the existential introdction, the term a has not yet been introdced into the derivation and its se can therefore
More informationFormal Methods for Deriving Element Equations
Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact
More informationComplex Variables. For ECON 397 Macroeconometrics Steve Cunningham
Comple Variables For ECON 397 Macroeconometrics Steve Cnningham Open Disks or Neighborhoods Deinition. The set o all points which satis the ineqalit
More informationFRÉCHET KERNELS AND THE ADJOINT METHOD
PART II FRÉCHET KERNES AND THE ADJOINT METHOD 1. Setp of the tomographic problem: Why gradients? 2. The adjoint method 3. Practical 4. Special topics (sorce imaging and time reversal) Setp of the tomographic
More informationEN Applied Optimal Control Lecture 8: Dynamic Programming October 10, 2018
EN530.603 Applied Optimal Control Lecture 8: Dynamic Programming October 0, 08 Lecturer: Marin Kobilarov Dynamic Programming (DP) is conerned with the computation of an optimal policy, i.e. an optimal
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More informationPrandl established a universal velocity profile for flow parallel to the bed given by
EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient
More informationLinear System Theory (Fall 2011): Homework 1. Solutions
Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts
More informationReflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by
Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introdction The transmission line eqations are given by, I z, t V z t l z t I z, t V z, t c z t (1) (2) Where, c is the per-nit-length
More informationCDS 110b: Lecture 1-2 Introduction to Optimal Control
CDS 110b: Lectre 1-2 Introdction to Optimal Control Richard M. Mrray 4 Janary 2006 Goals: Introdce the problem of optimal control as method of trajectory generation State the maimm principle and give eamples
More informationECON3120/4120 Mathematics 2, spring 2009
University of Oslo Department of Economics Arne Strøm ECON3/4 Mathematics, spring 9 Problem soltions for Seminar 4, 6 Febrary 9 (For practical reasons some of the soltions may inclde problem parts that
More information1. INTRODUCTION. A solution for the dark matter mystery based on Euclidean relativity. Frédéric LASSIAILLE 2009 Page 1 14/05/2010. Frédéric LASSIAILLE
Frédéric LASSIAILLE 2009 Page 1 14/05/2010 Frédéric LASSIAILLE email: lmimi2003@hotmail.com http://lmi.chez-alice.fr/anglais A soltion for the dark matter mystery based on Eclidean relativity The stdy
More informationSTEP Support Programme. STEP III Hyperbolic Functions: Solutions
STEP Spport Programme STEP III Hyperbolic Fnctions: Soltions Start by sing the sbstittion t cosh x. This gives: sinh x cosh a cosh x cosh a sinh x t sinh x dt t dt t + ln t ln t + ln cosh a ln ln cosh
More informationCHARACTERIZATIONS OF EXPONENTIAL DISTRIBUTION VIA CONDITIONAL EXPECTATIONS OF RECORD VALUES. George P. Yanev
Pliska Std. Math. Blgar. 2 (211), 233 242 STUDIA MATHEMATICA BULGARICA CHARACTERIZATIONS OF EXPONENTIAL DISTRIBUTION VIA CONDITIONAL EXPECTATIONS OF RECORD VALUES George P. Yanev We prove that the exponential
More informationPendulum Equations and Low Gain Regime
WIR SCHAFFEN WISSEN HEUTE FÜR MORGEN Sven Reiche :: SwissFEL Beam Dynamics Grop :: Pal Scherrer Institte Pendlm Eqations and Low Gain Regime CERN Accelerator School FELs and ERLs Interaction with Radiation
More informationCHAPTER 2 THE MAXIMUM PRINCIPLE: CONTINUOUS TIME. Chapter2 p. 1/67
CHAPTER 2 THE MAXIMUM PRINCIPLE: CONTINUOUS TIME Chapter2 p. 1/67 THE MAXIMUM PRINCIPLE: CONTINUOUS TIME Main Purpose: Introduce the maximum principle as a necessary condition to be satisfied by any optimal
More informationOptimal Control, Statistics and Path Planning
PERGAMON Mathematical and Compter Modelling 33 (21) 237 253 www.elsevier.nl/locate/mcm Optimal Control, Statistics and Path Planning C. F. Martin and Shan Sn Department of Mathematics and Statistics Texas
More informationOptimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance
Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for
More informationA Single Species in One Spatial Dimension
Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,
More informationUNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL
8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a
More information1. Tractable and Intractable Computational Problems So far in the course we have seen many problems that have polynomial-time solutions; that is, on
. Tractable and Intractable Comptational Problems So far in the corse we have seen many problems that have polynomial-time soltions; that is, on a problem instance of size n, the rnning time T (n) = O(n
More informationUncertainties of measurement
Uncertainties of measrement Laboratory tas A temperatre sensor is connected as a voltage divider according to the schematic diagram on Fig.. The temperatre sensor is a thermistor type B5764K [] with nominal
More informationDiscussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli
1 Introdction Discssion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen Department of Economics, University of Copenhagen and CREATES,
More informationDecision Oriented Bayesian Design of Experiments
Decision Oriented Bayesian Design of Experiments Farminder S. Anand*, Jay H. Lee**, Matthew J. Realff*** *School of Chemical & Biomoleclar Engineering Georgia Institte of echnology, Atlanta, GA 3332 USA
More informationGround Rules. PC1221 Fundamentals of Physics I. Position and Displacement. Average Velocity. Lectures 7 and 8 Motion in Two Dimensions
PC11 Fndamentals of Physics I Lectres 7 and 8 Motion in Two Dimensions A/Prof Tay Sen Chan 1 Grond Rles Switch off yor handphone and paer Switch off yor laptop compter and keep it No talkin while lectre
More informationOptimal Mean-Variance Portfolio Selection
Research Report No. 14, 2013, Probab. Statist. Grop Manchester 25 pp Optimal Mean-Variance Portfolio Selection J. L. Pedersen & G. Peskir Assming that the wealth process X is generated self-financially
More informationSetting The K Value And Polarization Mode Of The Delta Undulator
LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions
More informationTheoretical and Experimental Implementation of DC Motor Nonlinear Controllers
Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,
More informationAffine Invariant Total Variation Models
Affine Invariant Total Variation Models Helen Balinsky, Alexander Balinsky Media Technologies aboratory HP aboratories Bristol HP-7-94 Jne 6, 7* Total Variation, affine restoration, Sobolev ineqality,
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationSubcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany
Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2
More informationOPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1
OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body
More informationA Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time
A Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time Tomas Björk Department of Finance, Stockholm School of Economics tomas.bjork@hhs.se Agatha Mrgoci Department of Economics Aarhs
More informationInformation Source Detection in the SIR Model: A Sample Path Based Approach
Information Sorce Detection in the SIR Model: A Sample Path Based Approach Kai Zh and Lei Ying School of Electrical, Compter and Energy Engineering Arizona State University Tempe, AZ, United States, 85287
More informationLecture 8: September 26
10-704: Information Processing and Learning Fall 2016 Lectrer: Aarti Singh Lectre 8: September 26 Note: These notes are based on scribed notes from Spring15 offering of this corse. LaTeX template cortesy
More informationThe Scalar Conservation Law
The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan
More information4 Exact laminar boundary layer solutions
4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved
More informationPhysicsAndMathsTutor.com
. Two smooth niform spheres S and T have eqal radii. The mass of S is 0. kg and the mass of T is 0.6 kg. The spheres are moving on a smooth horizontal plane and collide obliqely. Immediately before the
More informationFormules relatives aux probabilités qui dépendent de très grands nombers
Formles relatives ax probabilités qi dépendent de très grands nombers M. Poisson Comptes rends II (836) pp. 603-63 In the most important applications of the theory of probabilities, the chances of events
More informationSufficient Optimality Condition for a Risk-Sensitive Control Problem for Backward Stochastic Differential Equations and an Application
Jornal of Nmerical Mathematics and Stochastics, 9(1) : 48-6, 17 http://www.jnmas.org/jnmas9-4.pdf JNM@S Eclidean Press, LLC Online: ISSN 151-3 Sfficient Optimality Condition for a Risk-Sensitive Control
More informationIntegration of Basic Functions. Session 7 : 9/23 1
Integration o Basic Fnctions Session 7 : 9/3 Antiderivation Integration Deinition: Taking the antiderivative, or integral, o some nction F(), reslts in the nction () i ()F() Pt simply: i yo take the integral
More informationHADAMARD-PERRON THEOREM
HADAMARD-PERRON THEOREM CARLANGELO LIVERANI. Invariant manifold of a fixed point He we will discss the simplest possible case in which the existence of invariant manifolds arises: the Hadamard-Perron theorem.
More informationGravitational Instability of a Nonrotating Galaxy *
SLAC-PUB-536 October 25 Gravitational Instability of a Nonrotating Galaxy * Alexander W. Chao ;) Stanford Linear Accelerator Center Abstract Gravitational instability of the distribtion of stars in a galaxy
More informationPhysicsAndMathsTutor.com
C Integration - By sbstittion PhysicsAndMathsTtor.com. Using the sbstittion cos +, or otherwise, show that e cos + sin d e(e ) (Total marks). (a) Using the sbstittion cos, or otherwise, find the eact vale
More informationChapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method
CIVL 7/87 Chater - The Stiffness Method / Chater Introdction to the Stiffness (Dislacement) Method Learning Objectives To define the stiffness matrix To derive the stiffness matrix for a sring element
More informationChem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions
Chem 4501 Introdction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics Fall Semester 2017 Homework Problem Set Nmber 10 Soltions 1. McQarrie and Simon, 10-4. Paraphrase: Apply Eler s theorem
More informationMove Blocking Strategies in Receding Horizon Control
Move Blocking Strategies in Receding Horizon Control Raphael Cagienard, Pascal Grieder, Eric C. Kerrigan and Manfred Morari Abstract In order to deal with the comptational brden of optimal control, it
More informationFOUNTAIN codes [3], [4] provide an efficient solution
Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn
More informationMathematical Analysis of Nipah Virus Infections Using Optimal Control Theory
Jornal of Applied Mathematics and Physics, 06, 4, 099- Pblished Online Jne 06 in SciRes. http://www.scirp.org/jornal/jamp http://dx.doi.org/0.436/jamp.06.464 Mathematical Analysis of Nipah Virs nfections
More informationSimplified Identification Scheme for Structures on a Flexible Base
Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles
More information10.4 Solving Equations in Quadratic Form, Equations Reducible to Quadratics
. Solving Eqations in Qadratic Form, Eqations Redcible to Qadratics Now that we can solve all qadratic eqations we want to solve eqations that are not eactl qadratic bt can either be made to look qadratic
More informationρ u = u. (1) w z will become certain time, and at a certain point in space, the value of
THE CONDITIONS NECESSARY FOR DISCONTINUOUS MOTION IN GASES G I Taylor Proceedings of the Royal Society A vol LXXXIV (90) pp 37-377 The possibility of the propagation of a srface of discontinity in a gas
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationMath 273b: Calculus of Variations
Math 273b: Calcls of Variations Yacob Kreh Homework #3 [1] Consier the 1D length fnctional minimization problem min F 1 1 L, or min 1 + 2, for twice ifferentiable fnctions : [, 1] R with bonary conitions,
More informationModelling by Differential Equations from Properties of Phenomenon to its Investigation
Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University
More informationStep-Size Bounds Analysis of the Generalized Multidelay Adaptive Filter
WCE 007 Jly - 4 007 London UK Step-Size onds Analysis of the Generalized Mltidelay Adaptive Filter Jnghsi Lee and Hs Chang Hang Abstract In this paper we analyze the bonds of the fixed common step-size
More informationAnalytical Investigation of Hyperbolic Equations via He s Methods
American J. of Engineering and Applied Sciences (4): 399-47, 8 ISSN 94-7 8 Science Pblications Analytical Investigation of Hyperbolic Eqations via He s Methods D.D. Ganji, M. Amini and A. Kolahdooz Department
More informationKrauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.
Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore
More informationBridging the Gap Between Multigrid, Hierarchical, and Receding-Horizon Control
Bridging the Gap Between Mltigrid, Hierarchical, and Receding-Horizon Control Victor M. Zavala Mathematics and Compter Science Division Argonne National Laboratory Argonne, IL 60439 USA (e-mail: vzavala@mcs.anl.gov).
More informationConditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane
Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching
More informationPontryagin s Minimum Principle 1
ECE 680 Fall 2013 Pontryagin s Minimum Principle 1 In this handout, we provide a derivation of the minimum principle of Pontryagin, which is a generalization of the Euler-Lagrange equations that also includes
More informationClassify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.
Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports
More informationApproximate Solution of Convection- Diffusion Equation by the Homotopy Perturbation Method
Gen. Math. Notes, Vol. 1, No., December 1, pp. 18-114 ISSN 19-7184; Copyright ICSRS Pblication, 1 www.i-csrs.org Available free online at http://www.geman.in Approximate Soltion of Convection- Diffsion
More informationFeature Selection for Neuro-Dynamic Programming
Featre Selection for Nero-Dynamic Programming Day Hang Wei Chen Prashant Mehta Sean Meyn Amit Srana October 4, 211 Abstract Nero-Dynamic Programming encompasses techniqes from both reinforcement learning
More informationElements of Coordinate System Transformations
B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationLecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2
BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE
More information5. The Bernoulli Equation
5. The Bernolli Eqation [This material relates predominantly to modles ELP034, ELP035] 5. Work and Energy 5. Bernolli s Eqation 5.3 An example of the se of Bernolli s eqation 5.4 Pressre head, velocity
More informationRemarks on strongly convex stochastic processes
Aeqat. Math. 86 (01), 91 98 c The Athor(s) 01. This article is pblished with open access at Springerlink.com 0001-9054/1/010091-8 pblished online November 7, 01 DOI 10.1007/s00010-01-016-9 Aeqationes Mathematicae
More informationAMS 212B Perturbation Methods Lecture 05 Copyright by Hongyun Wang, UCSC
AMS B Pertrbation Methods Lectre 5 Copright b Hongn Wang, UCSC Recap: we discssed bondar laer of ODE Oter epansion Inner epansion Matching: ) Prandtl s matching ) Matching b an intermediate variable (Skip
More informationA Note on Irreducible Polynomials and Identity Testing
A Note on Irrecible Polynomials an Ientity Testing Chanan Saha Department of Compter Science an Engineering Inian Institte of Technology Kanpr Abstract We show that, given a finite fiel F q an an integer
More informationBertrand s Theorem. October 8, µr 2 + V (r) 0 = dv eff dr. 3 + dv. f (r 0 )
Bertrand s Theorem October 8, Circlar orbits The eective potential, V e = has a minimm or maximm at r i and only i so we mst have = dv e L µr + V r = L µ 3 + dv = L µ 3 r r = L µ 3 At this radis, there
More information10.2 Solving Quadratic Equations by Completing the Square
. Solving Qadratic Eqations b Completing the Sqare Consider the eqation ( ) We can see clearl that the soltions are However, What if the eqation was given to s in standard form, that is 6 How wold we go
More informationIntrodction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc
A New Architectre for a Parallel Finite Field Mltiplier with Low Complexity Based on Composite Fields Christof Paar y IEEE Transactions on Compters, Jly 996, vol 45, no 7, pp 856-86 Abstract In this paper
More informationGarret Sobczyk s 2x2 Matrix Derivation
Garret Sobczyk s x Matrix Derivation Krt Nalty May, 05 Abstract Using matrices to represent geometric algebras is known, bt not necessarily the best practice. While I have sed small compter programs to
More informationA Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations
Applied Mathematics, 05, 6, 04-4 Pblished Online November 05 in SciRes. http://www.scirp.org/jornal/am http://d.doi.org/0.46/am.05.685 A Comptational Stdy with Finite Element Method and Finite Difference
More informationTwo-media boundary layer on a flat plate
Two-media bondary layer on a flat plate Nikolay Ilyich Klyev, Asgat Gatyatovich Gimadiev, Yriy Alekseevich Krykov Samara State University, Samara,, Rssia Samara State Aerospace University named after academician
More informationElliptic Curves Spring 2013
18.783 Elliptic Crves Spring 2013 Lectre #25 05/14/2013 In or final lectre we give an overview of the proof of Fermat s Last Theorem. Or goal is to explain exactly what Andrew Wiles [14], with the assistance
More information), σ is a parameter, is the Euclidean norm in R d.
WEIGHTED GRAPH LAPLACIAN AND IMAGE INPAINTING ZUOQIANG SHI, STANLEY OSHER, AND WEI ZHU Abstract. Inspired by the graph Laplacian and the point integral method, we introdce a novel weighted graph Laplacian
More informationAn Investigation into Estimating Type B Degrees of Freedom
An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information
More informationApproach to a Proof of the Riemann Hypothesis by the Second Mean-Value Theorem of Calculus
Advances in Pre Mathematics, 6, 6, 97- http://www.scirp.org/jornal/apm ISSN Online: 6-384 ISSN Print: 6-368 Approach to a Proof of the Riemann Hypothesis by the Second Mean-Vale Theorem of Calcls Alfred
More informationDiscussion Papers Department of Economics University of Copenhagen
Discssion Papers Department of Economics University of Copenhagen No. 10-06 Discssion of The Forward Search: Theory and Data Analysis, by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen,
More informationSensitivity Analysis in Bayesian Networks: From Single to Multiple Parameters
Sensitivity Analysis in Bayesian Networks: From Single to Mltiple Parameters Hei Chan and Adnan Darwiche Compter Science Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwiche}@cs.cla.ed
More informationThe Cryptanalysis of a New Public-Key Cryptosystem based on Modular Knapsacks
The Cryptanalysis of a New Pblic-Key Cryptosystem based on Modlar Knapsacks Yeow Meng Chee Antoine Jox National Compter Systems DMI-GRECC Center for Information Technology 45 re d Ulm 73 Science Park Drive,
More informationHomework 5 Solutions
Q Homework Soltions We know that the colmn space is the same as span{a & a ( a * } bt we want the basis Ths we need to make a & a ( a * linearly independent So in each of the following problems we row
More informationNon-Lecture I: Linear Programming. Th extremes of glory and of shame, Like east and west, become the same.
The greatest flood has the soonest ebb; the sorest tempest the most sdden calm; the hottest love the coldest end; and from the deepest desire oftentimes enses the deadliest hate. Th extremes of glory and
More informationChapter 2 Difficulties associated with corners
Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces
More information