t dt t SCLP Bellman (1953) CLP (Dantzig, Tyndall, Grinold, Perold, Anstreicher 60's-80's) Anderson (1978) SCLP
|
|
- Suzan Wilkinson
- 6 years ago
- Views:
Transcription
1 Coninuous Linear Programming. Separaed Coninuous Linear Programming Bellman (1953) max c () u() d H () u () + Gsusds (,) () a () u (), < < CLP (Danzig, yndall, Grinold, Perold, Ansreicher 6's-8's) Anderson (1978) max c () u() d Gu() s ds + x() = a() Hu () + z () = b () xuz,,, < < (Anderson, Nash, Philpo, Pullan 8's-9's). SCLP Dualiy, Srucure, Soluion by Discreizaion. In 2, W. proposed a simplex algorihm o solve: max ( ) c u( ) d Gu() s ds + x() = a+ α Hu () + z () = b xuz,,, < < SCLP wih linear daa MS&E324, Sanford Universiy, Spring 22 6B-1 Gideon Weiss manufacuring & conrol
2 Pullan s Dual Formulaion: Primal Problem: max c ( ) u ( ) d Gu() s ds+ x() = a() Hu () + z () = b () uxz,,, < < Pullan s SCLP Dual Problem min a ( ) dπ( ) + b ( ) rd ( ) G π() + H r() q() = c π( ) =, π, rq,, < < * Pullan s SCLP heorem (Pullan) Assume Hu () b (), u () has bounded feasible region. For abc,, piecewise analyic, here is a soluion which is piecewise analyic, For ac, piecewise polynomial of degree n +1, b piecewise polynomial of degree n here is a soluion wih u () piecewise polynomial of degree n. For ac, piecewise linear, b piecewise consan, here is a soluion wih u () piecewise consan. Number of breakpoins finie bounded (no polynomial) Srong dualiy holds Pullan s Algorihm: Sequence of discreizaions, converging (in infinie number of seps) o he opimal soluion. Noe: Separaed is very differen from General: in he one x () = usds (), direc inegraion in he oher u () = usds () + K differenial equaions MS&E324, Sanford Universiy, Spring 22 6B-2 Gideon Weiss manufacuring & conrol
3 SCLP wih Symmeric Primal / Dual max c ( ) u ( ) d+ d ( ) y ( ) d Gu() s ds + Fy() + x() = a() Hu () = b () xu,, < < SCLP min a ( ) p ( ) d+ b ( ) r ( ) d Gpsds () + Hr () q () = c () F p() = d() pq,, < < SCLP * Here: up, are primal and dual conrols ( ), xq, are primal and dual saes ( ), yr, are primal and dual supplemenary (unresriced) Dimensions: G K J, H I J, F K L MS&E324, Sanford Universiy, Spring 22 6B-3 Gideon Weiss manufacuring & conrol
4 Weak Dualiy and Complemenary Slackness Dual Objecive = a ( ) p ( ) d+ b ( ) r ( ) d { u () s G ds+ y( ) F } p() + u( ) H r() d = u( ) { G p( s) ds + H r( )} d + y( ) F p( ) d u ( ) c ( ) d+ y ( ) d ( ) d= Primal Objecive Equaliy will hold if and only if x ( ) pd ( ) =, u ( ) qd ( ) =. Complemenary Slackness: almos everywhere x () > p ( ) = u () > q ( ) = Corollary: he following are equivalen (a) uxy,, pqr,, are complemenary slack feasible primal and dual soluions (b) hey are opimal and have he same objecive value (no dualiy gap) (Srong dualiy no dualiy gap - exiss in LP bu no in oher problems, unless addiional condiions are imposed). MS&E324, Sanford Universiy, Spring 22 6B-4 Gideon Weiss manufacuring & conrol
5 Boundary Values For infiniesimal values of we solve: max c ( ) u ( ) + d ( ) y ( ) Fy() + x() = a() Hu () = b () xu,, < < min a ( ) p ( ) + b ( ) r ( ) Hr () q () = c () F p() = d() pq,, < < hese separae ino 2 ses of primal and dual problems: max d ( ) y( ) Fy( ) + x( ) = a( ) boundary LP x( ), min b ( ) r( ) Hr ( ) q( ) = c( ) q( ), boundary LP * Assumpion (Feasibiliy and Boundedness): he boundary problems are feasible and bounded. SCLP is feasible and bounded for some range of MS&E324, Sanford Universiy, Spring 22 6B-5 Gideon Weiss manufacuring & conrol
6 he Associaed LP Our soluion of SCLP will evolve hrough soluion of he following LP / LP*, under varying sign resricions. max c «( ) u ( ) + d ( ) y «( ) Gu() + Fy«( ) + ξ() = a«( ) Hu () = b () u (), y «( )" U", < < LP min a «( ) p ( ) + b ()«( r ) G p( ) + H r«( ) + θ( ) = c«( ) F p( ) = d( ) p, r«" U", < < LP * Noe: he ime indices on he variables are only labels, we label he variables o indicae heir role: ξ = x«, y«, θ = q«, r«assumpion (Non-Degeneracy): he r.h.s. of LP / LP* is (a.e. ) in general posiion o G F I column spaces of H, G H I F SCLP is non-degenerae: Soluion is unique and a an exreme poin in he funcion space of he soluions. MS&E324, Sanford Universiy, Spring 22 6B-6 Gideon Weiss manufacuring & conrol
7 he Srucure heorem Consider any funcions u ( ), x ( ), y ( ), p ( ), q ( ), r ( ), < < which saisfy up,, xyqr,,, are absoluely coninuous. Hence xyqr,,, have derivaives a.e. x( ), y( ) solve he boundary-lp, q( ), r( ) solve he boundary-lp*. For xy, hese are values a he lef boundary, qr, run in reversed ime so his is a he righ boundary. A all imes, he rae funcions (derivaives) u(), ξ(),«( y ), p( ), θ( ), r«( ), < < are complemenary slack feasible soluions for LP / LP*. hus hey are opimal soluions for LP / LP* under some sign configuraion of ξ, θ. By non-degeneracy hey are basic and unique soluions. A all imes, if xk() >, ξ k() is basic, if q ( ) >, θ ( ) is basic. j x () q (), < <. j heorem: A soluion o SCLP / SCLP* is opimal if (and only if for linear daa) i has hese properies. MS&E324, Sanford Universiy, Spring 22 6B-7 Gideon Weiss manufacuring & conrol
8 Resricion o he Linear Daa Case Problem: max ( γ + ( ) c ) u() d+ d y() d Gu() s ds + Fy() + x() = α + a Hu () xu,, < < min ( α + ( ) a ) p() d+ b r() d Gpsds () + Hr () q () = γ + c F p() = d pq,, < < Boundary: max dy min N br Fy + x = α = b SCLP ( linear daa) SCLP Hr q = γ * ( linear daa) x, q, Associaed LP max cu + d y«min ap + br «Gu+ Fy«+ ξ = a Gp + Hr «θ = c Hu = b Fp = d u, y«" U", p, r«" U" N N N MS&E324, Sanford Universiy, Spring 22 6B-8 Gideon Weiss manufacuring & conrol
9 he Algorihm he algorihm consrucs a soluion for SCLP / SCLP* for he whole range of ime horizons, < <. Corollary 1: he soluion for any ime horizon wih he 1 excepion of a finie lis: = ( ) < ( ) L < ( R) <, is characerized by a sequence of adjacen bases of LP, B 1, K, B N. For a paricular here are ime breakpoins = < 1L < N = such ha for n < < n n n n n n u, ξ,«y, p, θ,«r are basic B n. 1 n, In he single pivo Bn Bn +1 le v leave he basis if v = ξ k hen we have equaion xk( n)=. if v= hen we have equaion q ( ) =. u j Corollary 2: he soluion of he equaions xk( n)=, qj( n) =, ogeher wih: τ1+ L+ τ N = ( τ n = n n 1) deermines he soluion for a whole range of ime horizons, ( r 1) ( r) < < <. j n MS&E324, Sanford Universiy, Spring 22 6B-9 Gideon Weiss manufacuring & conrol
10 Algorihmic Srucure of Soluion: For ime horizon, Break poins: = < 1 < L < N = τ n = n n 1 N 1 Equaions: Adjacen Bn Bn+1 : = vn leaves: n n k k n k k m v = ξ x ( ) = x + ξ τ =, m= 1 n j j n j N N m j v = u q ( ) = q + θ τ =. and equaion N: τ1 + L+ τ N = m m= n+ 1 Slack Inequaliies, a local minima of x m k, q : k n k n n + 1 k n k k m ξ <, ξ > or n= N x ( ) = x + ξ τ, n j n+ 1 j n= or θ >, θ < q ( ) = q + θ τ. j m= 1 m j n j N N m j m= n+ 1 m MS&E324, Sanford Universiy, Spring 22 6B-1 Gideon Weiss manufacuring & conrol
11 Algorihm: ( ) ( 1) Sep 1: Soluion for < < is given by he soluion of boundary-lp / boundary-lp*. Sep r: Exend he soluion o nex ime horizon range. ( r) 1L 1 A d + τ = σ B I g Sep R+1 Sop when nex can be, or when soluion infeasible. MS&E324, Sanford Universiy, Spring 22 6B-11 Gideon Weiss manufacuring & conrol
12 Algorihm Case iii Buffer empies a x 1 x 2 < (l) x 2 x 1 ξ 1 leaves = (l) x 1 x 2 u 2 eners > (l) MS&E324, Sanford Universiy, Spring 22 6B-12 Gideon Weiss manufacuring & conrol
13 Algorihm Case ii Inerval shrinks a n B' C ξ 1 leaves x 1 q C B'' u 3 leaves 3 B' C B'' < (l) B' B'' ξ 1 u 3 leave x 1 q 3 B' B'' = (l) x 1 q 3 D B'' ξ 1 leaves B' D u 3 leaves B' D B'' > (l) MS&E324, Sanford Universiy, Spring 22 6B-13 Gideon Weiss manufacuring & conrol
14 Algorihm Case ii, Subproblem Inerval shrinks a n : subproblem for addiional inervals B' C ξ 1 leaves x 1 q C B'' u 3 leaves 3 B' C B'' < (l) B' B'' ξ 1 u 3 leave x 1 q 3 B' B'' = (l) subproblem E 1 D E 2 B' B'' ξ 1 u 3 leave x 1 q 3 B' B'' = (l) E 1 D 3u leaves x 1 q E B'' 1ξ leaves 3 2 B' E 1 D E 2 B'' > (l) Subproblem solved by recursive call o he algorihm, wih smaller, modified problem. MS&E324, Sanford Universiy, Spring 22 6B-14 Gideon Weiss manufacuring & conrol
15 Example 1: Pullan solved he following ASSE DISPOSAL PROBLEM in his 1993 paper: G= 1 H, = [ 1 2], 1 a1() = 4+, a2() = 3+ 2, b1() = 1, c1() = c2() = 2, < < 2. 2 max ( 2 u ) ( ) + ( 2 u ) ( d ) 1 2 u () s ds + x () = u2() s ds + x2() = 3+ 2 u1() + u2() + z1() = 1 uxz,,, < < Our formulaion, for arbirary ime horizon is: max ( ) u1( ) + ( ) u2( ) d u1() s ds + x1() = 4 + u2() s ds + x2() = 3+ 2 u1() + u2() + z1() = 1 uxz,,, < < MS&E324, Sanford Universiy, Spring 22 6B-15 Gideon Weiss manufacuring & conrol
16 Example 1: he opimal soluion is: 1 2 I is solved in 3 seps: For ime horizons < < 4 9 : Empy buffer 1 only, u () = 1, x () = 4 9, x () = For ime horizons 4 < <2 : keep buffer 1 empy, and empy buffer 2: u1() = 1, x1() =, u2() =, x2() = For ime horizons 2 < < : keep buffers 1 and 2 empy: u () = 1, x () =, u () = 2, x () = NOW WACH PULLAN SOLVE I: MS&E324, Sanford Universiy, Spring 22 6B-16 Gideon Weiss manufacuring & conrol
17 MS&E324, Sanford Universiy, Spring 22 6B-17 Gideon Weiss manufacuring & conrol
18 Example 2: A 2 machine 3 buffer example, m2 > m1+ m3: m m1 2 m 2 m > m + m Variaion 2 1 (Soluion of his example is join wih Florin Foram, 1992) MS&E324, Sanford Universiy, Spring 22 6B-18 Gideon Weiss manufacuring & conrol
19 S o l u i o n () 1, Buffer 3 drains B = { ξ, ξ, ξ, u, z } () 1 ( 2), Buffer 2 drains Opimal basis B = { ξ, ξ, u, u, z } Subproblem B = { ξ, ξ, u, ξ, u } ( 2) ( 3), Flow ou of buffer 1 unil inerval 3 shrinks B = { ξ, u, u, u, z } ( 3) ( 4), Basis B 3 repalced by B = { ξ, u, u, ξ, u } MS&E324, Sanford Universiy, Spring 22 6B-19 Gideon Weiss manufacuring & conrol ( 4), empy B = { u, u, u, z, z }
20 Example 3: 5 machine 2 buffer re-enran example I II III IV V MS&E324, Sanford Universiy, Spring 22 6B-2 Gideon Weiss manufacuring & conrol
21 S o l u i o n : MS&E324, Sanford Universiy, Spring 22 6B-21 Gideon Weiss manufacuring & conrol
22 Example 4: A coninuous ime Leonief Sysem X MS&E324, Sanford Universiy, Spring B-22 3 Gideon Weiss 4 manufacuring & conrol 5
23 Evoluion of he soluion: x 4 q 6 r2 q 6 x 3 x 4 x 3 r 2 r 2 x 3 x 4 q 12 x 2 MS&E324, Sanford Universiy, Spring 22 6B-23 Gideon Weiss manufacuring & conrol
Control of computer chip semi-conductor wafer fabs
Conrol of compuer chip semi-conducor wafer fabs Cos: x9 $ Reurn: years Cycle ime: 6 weeks WIP: 6, wafers, 8x6 $ Challenge: Conrol he queues a ~5 work seps Conrol of Manufacuring Sysems:! Sochasic racking
More informationAn introduction to the theory of SDDP algorithm
An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking
More informationSolutions to Assignment 1
MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we
More informationLecture 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 informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More informationHybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems
Hybrid Conrol and Swiched Sysems Lecure #3 Wha can go wrong? Trajecories of hybrid sysems João P. Hespanha Universiy of California a Sana Barbara Summary 1. Trajecories of hybrid sysems: Soluion o a hybrid
More informationF This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back.
Lecure 8 Las ime: Semi-free configuraion design This is equivalen o: Noe ns, ener he sysem a he same place. is fixed. We design C (and perhaps B. We mus sabilize if i is given as unsable. Cs ( H( s = +
More information3.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 informationME 391 Mechanical Engineering Analysis
Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of
More informationENGI 9420 Engineering Analysis Assignment 2 Solutions
ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion
More information1 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 informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationBackward stochastic dynamics on a filtered probability space
Backward sochasic dynamics on a filered probabiliy space Gechun Liang Oxford-Man Insiue, Universiy of Oxford based on join work wih Terry Lyons and Zhongmin Qian Page 1 of 15 gliang@oxford-man.ox.ac.uk
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationdt = 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 informationAn Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.
1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard
More informationChapter 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 informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More informationu(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x
. 1 Mah 211 Homework #3 February 2, 2001 2.4.3. y + (2/x)y = (cos x)/x 2 Answer: Compare y + (2/x) y = (cos x)/x 2 wih y = a(x)x + f(x)and noe ha a(x) = 2/x. Consequenly, an inegraing facor is found wih
More informationPredator - 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 informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationAlgorithmic Trading: Optimal Control PIMS Summer School
Algorihmic Trading: Opimal Conrol PIMS Summer School Sebasian Jaimungal, U. Torono Álvaro Carea,U. Oxford many hanks o José Penalva,(U. Carlos III) Luhui Gan (U. Torono) Ryan Donnelly (Swiss Finance Insiue,
More informationOptima 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 information4 Sequences of measurable functions
4 Sequences of measurable funcions 1. Le (Ω, A, µ) be a measure space (complee, afer a possible applicaion of he compleion heorem). In his chaper we invesigae relaions beween various (nonequivalen) convergences
More informationStructural Dynamics and Earthquake Engineering
Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/
More informationMath 334 Fall 2011 Homework 11 Solutions
Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then
More informationTEACHER NOTES MATH NSPIRED
Naural Logarihm Mah Objecives Sudens will undersand he definiion of he naural logarihm funcion in erms of a definie inegral. Sudens will be able o use his definiion o relae he value of he naural logarihm
More information8. 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 informationUtility maximization in incomplete markets
Uiliy maximizaion in incomplee markes Marcel Ladkau 27.1.29 Conens 1 Inroducion and general seings 2 1.1 Marke model....................................... 2 1.2 Trading sraegy.....................................
More informationOptimality Conditions for Unconstrained Problems
62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x
More informationSZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1
SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationEchocardiography Project and Finite Fourier Series
Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every
More informationInventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions
Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.
More informationMA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions
MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by
More informationAnn. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:
Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS
More informationSUFFICIENT 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 informationBU Macro BU Macro Fall 2008, Lecture 4
Dynamic Programming BU Macro 2008 Lecure 4 1 Ouline 1. Cerainy opimizaion problem used o illusrae: a. Resricions on exogenous variables b. Value funcion c. Policy funcion d. The Bellman equaion and an
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationAn Introduction to Stochastic Programming: The Recourse Problem
An Inroducion o Sochasic Programming: he Recourse Problem George Danzig and Phil Wolfe Ellis Johnson, Roger Wes, Dick Cole, and Me John Birge Where o look in he ex pp. 6-7, Secion.2.: Inroducion o sochasic
More informationEMS SCM joint meeting. On stochastic partial differential equations of parabolic type
EMS SCM join meeing Barcelona, May 28-30, 2015 On sochasic parial differenial equaions of parabolic ype Isván Gyöngy School of Mahemaics and Maxwell Insiue Edinburgh Universiy 1 I. Filering problem II.
More informationKEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering
More informationWEEK-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 informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More informationThe 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 informationA New Perturbative Approach in Nonlinear Singularity Analysis
Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The
More information1 Solutions to selected problems
1 Soluions o seleced problems 1. Le A B R n. Show ha in A in B bu in general bd A bd B. Soluion. Le x in A. Then here is ɛ > 0 such ha B ɛ (x) A B. This shows x in B. If A = [0, 1] and B = [0, 2], hen
More informationarxiv:math/ v1 [math.nt] 3 Nov 2005
arxiv:mah/0511092v1 [mah.nt] 3 Nov 2005 A NOTE ON S AND THE ZEROS OF THE RIEMANN ZETA-FUNCTION D. A. GOLDSTON AND S. M. GONEK Absrac. Le πs denoe he argumen of he Riemann zea-funcion a he poin 1 + i. Assuming
More informationThe Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales
Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions
More informationProblem Set on Differential Equations
Problem Se on Differenial Equaions 1. Solve he following differenial equaions (a) x () = e x (), x () = 3/ 4. (b) x () = e x (), x (1) =. (c) xe () = + (1 x ()) e, x () =.. (An asse marke model). Le p()
More informationSliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game
Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp
More informationCorrespondence should be addressed to Nguyen Buong,
Hindawi Publishing Corporaion Fixed Poin Theory and Applicaions Volume 011, Aricle ID 76859, 10 pages doi:101155/011/76859 Research Aricle An Implici Ieraion Mehod for Variaional Inequaliies over he Se
More informationRobotics 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 informationPhysics 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 informationEffects of Coordinate Curvature on Integration
Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold
More informationThe fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation
Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationHamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:
M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno
More informationHamilton Jacobi equations
Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationSPLICING OF TIME OPTIMAL CONTROLS
Dynamic Sysems and Applicaions 21 (212) 169-186 SPLICING OF TIME OPTIMAL CONTROLS H. O. FATTORINI Deparmen of Mahemaics, Universiy of California Los Angeles, California 995-1555, USA ABSTRACT. Two conrols
More informationBasilio 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 informationExample on p. 157
Example 2.5.3. Le where BV [, 1] = Example 2.5.3. on p. 157 { g : [, 1] C g() =, g() = g( + ) [, 1), var (g) = sup g( j+1 ) g( j ) he supremum is aken over all he pariions of [, 1] (1) : = < 1 < < n =
More informationInstructor: Barry McQuarrie Page 1 of 5
Procedure for Solving radical equaions 1. Algebraically isolae one radical by iself on one side of equal sign. 2. Raise each side of he equaion o an appropriae power o remove he radical. 3. Simplify. 4.
More informationarxiv: v1 [math.fa] 9 Dec 2018
AN INVERSE FUNCTION THEOREM CONVERSE arxiv:1812.03561v1 [mah.fa] 9 Dec 2018 JIMMIE LAWSON Absrac. We esablish he following converse of he well-known inverse funcion heorem. Le g : U V and f : V U be inverse
More informationComments on Window-Constrained Scheduling
Commens on Window-Consrained Scheduling Richard Wes Member, IEEE and Yuing Zhang Absrac This shor repor clarifies he behavior of DWCS wih respec o Theorem 3 in our previously published paper [1], and describes
More informationExpert Advice for Amateurs
Exper Advice for Amaeurs Ernes K. Lai Online Appendix - Exisence of Equilibria The analysis in his secion is performed under more general payoff funcions. Wihou aking an explici form, he payoffs of he
More informationRepresentation of Stochastic Process by Means of Stochastic Integrals
Inernaional Journal of Mahemaics Research. ISSN 0976-5840 Volume 5, Number 4 (2013), pp. 385-397 Inernaional Research Publicaion House hp://www.irphouse.com Represenaion of Sochasic Process by Means of
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More informationln 2 1 ln y x c y C x
Lecure 14 Appendi B: Some sample problems from Boas Here are some soluions o he sample problems assigned for Chaper 8 8: 6 Soluion: We wan o find he soluion o he following firs order equaion using separaion
More informationPOSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION
Novi Sad J. Mah. Vol. 32, No. 2, 2002, 95-108 95 POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Hajnalka Péics 1, János Karsai 2 Absrac. We consider he scalar nonauonomous neural delay differenial
More informationComparison between the Discrete and Continuous Time Models
Comparison beween e Discree and Coninuous Time Models D. Sulsky June 21, 2012 1 Discree o Coninuous Recall e discree ime model Î = AIS Ŝ = S Î. Tese equaions ell us ow e populaion canges from one day o
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationMacroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3
Macroeconomic Theory Ph.D. Qualifying Examinaion Fall 2005 Comprehensive Examinaion UCLA Dep. of Economics You have 4 hours o complee he exam. There are hree pars o he exam. Answer all pars. Each par has
More informationLecture 4 Notes (Little s Theorem)
Lecure 4 Noes (Lile s Theorem) This lecure concerns one of he mos imporan (and simples) heorems in Queuing Theory, Lile s Theorem. More informaion can be found in he course book, Bersekas & Gallagher,
More informationSPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990),
SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F Trench SIAM J Marix Anal Appl 11 (1990), 601-611 Absrac Le T n = ( i j ) n i,j=1 (n 3) be a real symmeric
More informationOn Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems
MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 2, May 2013, pp. 209 227 ISSN 0364-765X (prin) ISSN 1526-5471 (online) hp://dx.doi.org/10.1287/moor.1120.0562 2013 INFORMS On Boundedness of Q-Learning Ieraes
More information5. Stochastic processes (1)
Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly
More informationSupplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence
Supplemen for Sochasic Convex Opimizaion: Faser Local Growh Implies Faser Global Convergence Yi Xu Qihang Lin ianbao Yang Proof of heorem heorem Suppose Assumpion holds and F (w) obeys he LGC (6) Given
More informationThe expectation value of the field operator.
The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More informationCHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang
CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS Professor Dae Ryook Yang Fall 200 Dep. of Chemical and Biological Engineering Korea Universiy CHE302 Process Dynamics and Conrol Korea Universiy
More informationEssential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems
Essenial Microeconomics -- 6.5: OPIMAL CONROL Consider he following class of opimizaion problems Max{ U( k, x) + U+ ( k+ ) k+ k F( k, x)}. { x, k+ } = In he language of conrol heory, he vecor k is he vecor
More informationO Q L N. Discrete-Time Stochastic Dynamic Programming. I. Notation and basic assumptions. ε t : a px1 random vector of disturbances at time t.
Econ. 5b Spring 999 C. Sims Discree-Time Sochasic Dynamic Programming 995, 996 by Chrisopher Sims. This maerial may be freely reproduced for educaional and research purposes, so long as i is no alered,
More informationA Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs
PROC. IEEE CONFERENCE ON DECISION AND CONTROL, 06 A Primal-Dual Type Algorihm wih he O(/) Convergence Rae for Large Scale Consrained Convex Programs Hao Yu and Michael J. Neely Absrac This paper considers
More informationMath 2214 Solution Test 1B Fall 2017
Mah 14 Soluion Tes 1B Fall 017 Problem 1: A ank has a capaci for 500 gallons and conains 0 gallons of waer wih lbs of sal iniiall. A soluion conaining of 8 lbsgal of sal is pumped ino he ank a 10 galsmin.
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More informationExam 1 Solutions. 1 Question 1. February 10, Part (A) 1.2 Part (B) To find equilibrium solutions, set P (t) = C = dp
Exam Soluions Februar 0, 05 Quesion. Par (A) To find equilibrium soluions, se P () = C = = 0. This implies: = P ( P ) P = P P P = P P = P ( + P ) = 0 The equilibrium soluion are hus P () = 0 and P () =..
More informationTheory 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 informationMath 4600: Homework 11 Solutions
Mah 46: Homework Soluions Gregory Handy [.] One of he well-known phenomenological (capuring he phenomena, bu no necessarily he mechanisms) models of cancer is represened by Gomperz equaion dn d = bn ln(n/k)
More informationMath 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 informationOscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,
More informationDISCRETE GRONWALL LEMMA AND APPLICATIONS
DISCRETE GRONWALL LEMMA AND APPLICATIONS JOHN M. HOLTE MAA NORTH CENTRAL SECTION MEETING AT UND 24 OCTOBER 29 Gronwall s lemma saes an inequaliy ha is useful in he heory of differenial equaions. Here is
More information