Stochastic Dual Dynamic integer Programming
|
|
- Peregrine Philip Horton
- 5 years ago
- Views:
Transcription
1 Stochastic Dual Dynamic integer Programming Shabbir Ahmed Georgia Tech Jikai Zou Andy Sun
2 Multistage IP Canonical deterministic formulation ( X T ) f t (x t,y t ):(x t 1,x t,y t ) 2 X t 8 t x t min x,y t=1 = state variables; = local/stage variables y t Linear objective and constraints A t x t 1 + B t x t + C t y t b t Bounded (mixed integer) variables
3 Multistage Stochastic IP Stochastic Data: { t =(f t,x t )} T t=1 Dynamics: t #... (x t 1,y t 1 )! (x t,y t )... Formulation: n min f 1 (x 1,y 1 )+E (x 1,y 1 )2X 1 h h + E n min f 2 (x 2,y 2 )+ (x 2,y 2 )2X 2 (x 1 ) min (x T,y T )2X T (x T 1 ) n oioio f T (x T,y T )
4 SIP Applications Energy Natural resources Healthcare Logis4cs Telecommunica4ons Finance 7 of 21 applica4ons are SIP
5 Trifecta of Complexities Integrality Dynamics Mul5stage IP Integer Programming Mul5stage LP MSIP Sta5c Stochas5c IP Mul5stage Stochas5c LP Sta5c Stochas5c LP Uncertainty
6 Scenario Tree Dynamic uncertainty modeling: Scenario tree Explicit construction/discretization Monte Carlo Sampling Markov Chain (implicit)
7 Extensive Form min x n,y n ( X p n f n (x n,y n ): (x a(n),x n,y n ) 2 X n 8 n 2 T n2t Explicit decision variables and constraints for each node Very large scale but deterministic Decomposition methods (MSLP): Nested Benders (Birge`85) Progressive Hedging (Rockafellar and Wets `91) )
8 Dynamic Programming Formulation: min x 1,y 1 f 1 (x 1,y 1 )+Q 1 (x 1 ): (x a(1),x 1,y 1 ) 2 X 1 Expected Cost-to-go (ECTG) function: Q n (x n ):= X q nm Q m (x n ) m2c(n) Cost-to-go function: Q m (x n )= min x m,y m {f m (x m,y m )+Q m (x m ): (x n,x m,y m ) 2 X m }
9 Dynamic Programming Formulation: min x 1,y 1 f 1 (x 1,y 1 )+Q 1 (x 1 ): (x a(1),x 1,y 1 ) 2 X 1 Expected Cost-to-go (ECTG) function: Q n (x n ):= X q nm Q m (x n ) m2c(n) Cost-to-go function: Polyhedral in LP setting Q m (x n )= min x m,y m {f m (x m,y m )+Q m (x m ): (x n,x m,y m ) 2 X m }
10 Stage-wise Independence N t Stage t has independent realizations Recombining scenario tree One expected cost-to-go function per stage Q n ( ) Q t ( ) 8 n 2 S t Q t ( ) implicitly encodes solution (policy): action under j-th realization in stage t given previous action x t 1 is (x t,y t ) = argmin x,y {f j,t (x, y)+q t (x) :(x t 1,x,y) 2 X j,t }
11 Stochastic Dual Dynamic Programming
12 Illustra5on of SDDP Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
13 Iter 1: Forward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
14 Iter 1: Forward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
15 Iter 1: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) Benders Cut x 3 x 2 x 1
16 Iter 1: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
17 Iter 1: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
18 Iter 2: Forward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
19 Iter 2: Forward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
20 Iter 2: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
21 Iter 2: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
22 Iter 2: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
23 Iter 3: Forward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
24 Iter 3: Forward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
25 Iter 3: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
26 Iter 3: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
27 Iter 3: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
28 Iter 4: Forward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
29 Iter 4: Forward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
30 Iter 4: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
31 Iter 4: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
32 Iter 4: Backward pass Q 2 (x 2 ) Q 3 (x 3 ) Q 1 (x 1 ) x 3 x 2 x 1
33 SDDP Bounds Stochastic algorithm Forward pass generates candidate solutions along sample paths, average of these objective values provides a statistical upper bound Optimal value of stage 1 problem (with lower approximation of ECTGF) provides a deterministic lower bound
34 SDDP Convergence Theorem: If sampling is done with replacement and basic dual optimal solutions are used to construct Benders cuts, then w.p. 1 SDDP converges in a finite number of iterations to an optimal policy for a multistage stochastic LP. Structure of LP value function structure is crucial Chen & Powell (1999) PhilpoL & Guan (2008) Shapiro (2011) Girardeau et al. (2014)
35 MSIP with Binary State Variables Integer variables destroy convexity of ECTG function Stage-wise decomposition relies on approximating ECTG functions by cuts i.e. convex polyhedral functions A function of binary variables can always be represented as a convex polyhedral function
36 Binary State Variables min x n,y n s.t. X p n f n (x n,y n ) n2t (x a(n),x n,y n ) 2 X n 8 n 2 T x n 2 {0, 1} d 8 n 2 T
37 Binary State Variables min x n,y n s.t. X p n f n (x n,y n ) n2t (x a(n),x n,y n ) 2 X n 8 n 2 T x n 2 {0, 1} d 8 n 2 T Theorem: [Under complete continuous recourse] A MSIP with d general state variables can be approximated to " accuracy by a MSIP with k binary state variables and p! d k = O d log " Proof: Sensitivity theorems for mixed integer linear programming (Cook et al. (1986))
38 Local Convexification (P ): min x n,y n s.t. X p n f n (x n,y n ) n2t (x a(n),x n,y n ) 2 X n 8 n 2 T x n 2 {0, 1} d 8 n 2 T
39 Local Convexification (P ): min x n,y n s.t. X p n f n (x n,y n ) n2t (x a(n),x n,y n ) 2 X n 8 n 2 T (Q) : x n 2 {0, 1} d 8 n 2 T X min p n f n (x n,y n ) x n,y n,z n s.t. n2t (z n,x n,y n ) 2 X n 8 n 2 T z n 2 [0, 1] d 8 n 2 T z n = x a(n) 8 n 2 T x n 2 {0, 1} d 8 n 2 T
40 Local Convexification (P ): min x n,y n s.t. X p n f n (x n,y n ) n2t (x a(n),x n,y n ) 2 X n 8 n 2 T x n 2 {0, 1} d 8 n 2 T (Q) : min x n,y n,z n s.t. X p n f n (x n,y n ) n2t (z n,x n,y n ) 2 X 0 n 8 n 2 T z n = x a(n) 8 n 2 T x n 2 {0, 1} d 8 n 2 T
41 Local Convexification (P ): min x n,y n s.t. X p n f n (x n,y n ) n2t (x a(n),x n,y n ) 2 X n 8 n 2 T (Q 0 ): min x n,y n,z n s.t. x n 2 {0, 1} d 8 n 2 T X p n f n (x n,y n ) n2t (z n,x n,y n ) 2 conv(x 0 n) 8 n 2 T z n = x a(n) 8 n 2 T x n 2 {0, 1} d 8 n 2 T Theorem: (P) and (Q ) are equivalent.
42 Binary State variables and a reformulation Node (forward) problem: Q j,t (x t 1 ) := min x,y,z s.t. f j,t (x, y)+q t (x) (z,x,y) 2 X j,t z = x t 1 z 2 [0, 1] d x 2 {0, 1} d j =1 Expected Cost-to-go function: Q t (x) = 1 NX Q j,t+1 (x) N j=1 t 1 t t +1 j =. j = N
43 ECTG Function Approximation Node (forward) problem: Q i (x j,t t 1) := min x,y,z s.t. f j,t (x, y)+ i t (x) (z,x,y) 2 X 0 j,t z = x t 1 Approx. Expected Cost-to-go function: i t(x) = min s.t. 1 N NX j=1 v`j,t+1 +( `j,t+1) > x 8` =1,...,i 1
44 ECTG Function Approximation Node (forward) problem: Q i (x j,t t 1) := min x,y,z s.t. f j,t (x, y)+ i t (x) (z,x,y) 2 X 0 j,t z = x t 1 Approx. Expected Cost-to-go function: i t(x) = min s.t. 1 N NX j=1 v`j,t+1 +( `j,t+1) > x 8` =1,...,i 1 Cut coefficients obtained from (relaxations of) node problems in stage t+1 in previous iterations
45 Cut Conditions and Convergence Valid: Q j,t (x) v i j,t +( i j,t) > x 8 x 2 {0, 1} d Tight: Q i+1 j,t (xi t 1) =v i j,t +( i j,t) > x i t 1 Finite: Possible cut coefficients finite Theorem: If sampling is done with replacement and cuts are valid, tight and finite, then w.p. 1 SDDiP converges in a finite number of iterations to an optimal policy for MSIP with binary state variables.
46 Convergence Proof Consider extensive tree structure since solutions are path dependent A solution/policy {x i n} n2t generated at iteration i is i specified by the approximate ECTG functions { t( )} T t=1 x i n 2 argminx,y s.t. i f n (x, y)+ tn (x) (x i a(n),x,y) 2 X n, x 2 {0, 1} d Lemma 1: If is optimal. i t n (x i n)=q tn (x i n) 8 n 2 T Proof: Backward induction of DP equations then {x i n} n2t
47 Convergence Proof (contd.) Let K =sup{i : {x i n} n2t is not optimal} Two types of iterations: (a) At least one approx ECTG function changes (b) Approx ECTG functions do not change so policy does not change K = K a + K b K = K a + P x Kx b Finite cuts mean finitely many approximations of ECTG K a < +1 Finite solution set means the summation is finite
48 Convergence Proof (contd.) Lemma 2: Proof: Pr[K x b < +1] =1 If solution is not optimal then by Lemma 1 there is a (last) node where approximation is loose By Borel-Cantelli lemma, w.p.1 forward pass will hit this node within finite iterations, then a tight cut is generated, and we have new approximation of ECTGF Thus 1 = Pr[Kb x < +1] apple Pr[K = Ka + P x Kx b < +1] Q.E.D
49 Relaxations and Cuts Benders Cut 1 N NX j=1 h i Q LP j,t+1 (xi t)+( j,t+1 ) > (x x i t) Solve LP relaxation Use LP dual solutions Valid and finite (using basic dual solutions) Not tight
50 Relaxations and Cuts Integer Cut (L vt+1 i+1 ) x xi t + vt+1 i+1 where v i+1 t+1 = 1 N NX j=1 Q i+1 j,t+1 (xi t) Laporte and Louveaux (1993) Solve IP Valid, tight and finite
51 Improved LL Cuts Angulo, A. and Dey. Improving the integer L-shaped Method, IJOC, Gustavo Angulo Santanu Dey
52 Relaxations and Cuts Lagrangian Cut NX Lj,t+1 ( j,t+1)+( j,t+1) > x t 1 N j=1 where L j,t+1 ( ) = Solve Lagrangian Dual min x,y,z s.t. j,t+1 2 argmax f j,t+1 (x, y)+ Valid and finite (using basic dual solutions) i+1 t+1 (x) > z (x, y, z) 2 X, z 2 [0, 1] d L j,t+1( )+ > x i t
53 Lagrangian Cut Theorem Lagrangian cut is tight. Primal characterization (convexification) due Lagrangian dual Local copies of state variables Binary nature of state variables (facial property) Recursive application
54 Lagrangian Cut Q(x) Q(x) = min y 1 + y 2 s.t. 2y 1 + y 2 3x 0 apple y 1 apple 2 0 apple y 2 apple 3 y 1 2 Z 1 x
55 Lagrangian Cut Cerisola et al. (2009) Thome et. al. (2013) Q(x) Q(x) = min y 1 + y 2 s.t. 2y 1 + y 2 3x 0 apple y 1 apple 2 0 apple y 2 apple 3 y 1 2 Z 1.5x 1 x
56 Lagrangian Cut Q(x) Q(x) = min y 1 + y 2 s.t. 2y 1 + y 2 3z z = x 0 apple y 1 apple 2 0 apple y 2 apple 3 0 apple z apple 1 y 1 2 Z 1 x
57 Lagrangian Cut Q(x) Q(x) = min y 1 + y 2 s.t. 2y 1 + y 2 3z z = x 0 apple y 1 apple 2 0 apple y 2 apple 3 0 apple z apple 1 y 1 2 Z 1 1+3x x
58 Lagrangian Cut Q(x) Q(x) = min y 1 + y 2 s.t. 2y 1 + y 2 3z z = x 0apple y 1 apple 2 0apple y 2 apple 3 z0 2 apple{0, z apple1} 1 y 1 2 Z 2x 1 x
59 Relaxations and Cuts Strengthened Bender s Cut 1 N NX j=1 Lj,t+1 ( LP j,t+1)+( LP j,t+1) > x t Solve LP relaxation Solve Lagrangian relaxation with LP dual values Valid and finite Not tight Dominates Benders cut Incomparable to Lagrangian cut
60 Computations: Test Case Generation Expansion Planning min s.t. P T t=1 (a> t x t + b > t y t ) A t y t apple P t =1 x 8 t e > y t = d t 8 t x t 2 Z n +, y t 2 R n + 8 t Uncertain demand and price Sampled scenarios Data from Jin et al. (2011) Six types of generators
61 Computations Small case: 10 Stages 3 nodes per stage (3^9 scenarios) Extensive form: >600K integer variables, >200K constraints CPLEX: after 2 hrs, gap = 7% Comparison of Cuts in SDDiP: Benders (B) Integer (I) Lagrangian (L) Strengthened Benders (SB)
62 Cut Comparison (Time) time FW I L B+I B+L S+I S+L S+I+L method
63 Cut Comparison (Gap) 0.6 gap 0.4 FW I L B+I B+L S+I S+L S+I+L method
64 Cut Comparison (Convergence)
65 Scalability Extensive form has > 11 trillion variables!
66 Summary SDDiP: Extension of SDDP to MSIP with binary state vars Exploit binary state variables and a reformulation with local copies of state variables Computations indicate Strengthened Benders + Integer Cuts work best Recent Applications: Hydroscheduling (M. Hjelmeland and A. Helseth) Unit commitment (J. Zou and X. Sun)
67 References G. Angulo, S. Ahmed, S.S. Dey, V. Kaibel. Forbidden vertices, Mathematics of Operations Research, vol.40, pp , G. Angulo, S. Ahmed, S.S. Dey. Improving the integer L-shaped method, INFORMS Journal on Computing, vol. 28, pp , M.N. Hjelmeland, J. Zou, A. Helseth, S. Ahmed. Nonconvex medium-term hydropower scheduling by Stochastic Dual Dynamic integer Programming, optimization-online.org, J. Zou, S. Ahmed and X. Sun. Multistage stochastic unit commitment using Stochastic Dual Dynamic Integer Programming, optimization-online.org, J. Zou, S. Ahmed, X. Sun. Stochastic Dual Dynamic Integer Programming, optimization-online.org, 2016.
Robust Dual Dynamic Programming
1 / 18 Robust Dual Dynamic Programming Angelos Georghiou, Angelos Tsoukalas, Wolfram Wiesemann American University of Beirut Olayan School of Business 31 May 217 2 / 18 Inspired by SDDP Stochastic optimization
More informationDASC: A DECOMPOSITION ALGORITHM FOR MULTISTAGE STOCHASTIC PROGRAMS WITH STRONGLY CONVEX COST FUNCTIONS
DASC: A DECOMPOSITION ALGORITHM FOR MULTISTAGE STOCHASTIC PROGRAMS WITH STRONGLY CONVEX COST FUNCTIONS Vincent Guigues School of Applied Mathematics, FGV Praia de Botafogo, Rio de Janeiro, Brazil vguigues@fgv.br
More informationWorst-case-expectation approach to optimization under uncertainty
Worst-case-expectation approach to optimization under uncertainty Wajdi Tekaya Joint research with Alexander Shapiro, Murilo Pereira Soares and Joari Paulo da Costa : Cambridge Systems Associates; : Georgia
More informationStochastic Dual Dynamic Programming Algorithm for Multistage Stochastic Programming
Stochastic Dual Dynamic Programg Algorithm for Multistage Stochastic Programg Final presentation ISyE 8813 Fall 2011 Guido Lagos Wajdi Tekaya Georgia Institute of Technology November 30, 2011 Multistage
More informationMultistage Stochastic Programming
IE 495 Lecture 21 Multistage Stochastic Programming Prof. Jeff Linderoth April 16, 2003 April 16, 2002 Stochastic Programming Lecture 21 Slide 1 Outline HW Fixes Multistage Stochastic Programming Modeling
More informationMultistage risk-averse asset allocation with transaction costs
Multistage risk-averse asset allocation with transaction costs 1 Introduction Václav Kozmík 1 Abstract. This paper deals with asset allocation problems formulated as multistage stochastic programming models.
More informationMultistage Stochastic Programming
Multistage Stochastic Programming John R. Birge University of Michigan Models - Long and short term - Risk inclusion Approximations - stages and scenarios Computation Slide Number 1 OUTLINE Motivation
More informationAssessing Policy Quality in Multi-stage Stochastic Programming
Assessing Policy Quality in Multi-stage Stochastic Programming Anukal Chiralaksanakul and David P. Morton Graduate Program in Operations Research The University of Texas at Austin Austin, TX 78712 January
More informationMultistage Stochastic Demand-side Management for Price-Making Major Consumers of Electricity in a Co-optimized Energy and Reserve Market
Multistage Stochastic Demand-side Management for Price-Making Major Consumers of Electricity in a Co-optimized Energy and Reserve Market Mahbubeh Habibian Anthony Downward Golbon Zakeri Abstract In this
More informationStochastic Dual Dynamic Programming
1 / 43 Stochastic Dual Dynamic Programming Operations Research Anthony Papavasiliou 2 / 43 Contents [ 10.4 of BL], [Pereira, 1991] 1 Recalling the Nested L-Shaped Decomposition 2 Drawbacks of Nested Decomposition
More informationThe value of multi-stage stochastic programming in capacity planning under uncertainty
The value of multi-stage stochastic programming in capacity planning under uncertainty Kai Huang and Shabbir Ahmed School of Industrial & Systems Engineering Georgia Institute of Technology April 26, 2005
More informationOptimal energy management and stochastic decomposition
Optimal energy management and stochastic decomposition F. Pacaud P. Carpentier J.P. Chancelier M. De Lara JuMP-dev workshop, 2018 ENPC ParisTech ENSTA ParisTech Efficacity 1/23 Motivation We consider a
More informationOn solving multistage stochastic programs with coherent risk measures
On solving multistage stochastic programs with coherent risk measures Andy Philpott Vitor de Matos y Erlon Finardi z August 13, 2012 Abstract We consider a class of multistage stochastic linear programs
More informationContents Critique 26. portfolio optimization 32
Contents Preface vii 1 Financial problems and numerical methods 3 1.1 MATLAB environment 4 1.1.1 Why MATLAB? 5 1.2 Fixed-income securities: analysis and portfolio immunization 6 1.2.1 Basic valuation of
More informationScenario reduction and scenario tree construction for power management problems
Scenario reduction and scenario tree construction for power management problems N. Gröwe-Kuska, H. Heitsch and W. Römisch Humboldt-University Berlin Institute of Mathematics Page 1 of 20 IEEE Bologna POWER
More information6.231 DYNAMIC PROGRAMMING LECTURE 10 LECTURE OUTLINE
6.231 DYNAMIC PROGRAMMING LECTURE 10 LECTURE OUTLINE Rollout algorithms Cost improvement property Discrete deterministic problems Approximations of rollout algorithms Discretization of continuous time
More informationProgressive Hedging for Multi-stage Stochastic Optimization Problems
Progressive Hedging for Multi-stage Stochastic Optimization Problems David L. Woodruff Jean-Paul Watson Graduate School of Management University of California, Davis Davis, CA 95616, USA dlwoodruff@ucdavis.edu
More informationReport for technical cooperation between Georgia Institute of Technology and ONS - Operador Nacional do Sistema Elétrico Risk Averse Approach
Report for technical cooperation between Georgia Institute of Technology and ONS - Operador Nacional do Sistema Elétrico Risk Averse Approach Alexander Shapiro and Wajdi Tekaya School of Industrial and
More informationDM559/DM545 Linear and integer programming
Department of Mathematics and Computer Science University of Southern Denmark, Odense May 22, 2018 Marco Chiarandini DM559/DM55 Linear and integer programming Sheet, Spring 2018 [pdf format] Contains Solutions!
More informationContract Theory in Continuous- Time Models
Jaksa Cvitanic Jianfeng Zhang Contract Theory in Continuous- Time Models fyj Springer Table of Contents Part I Introduction 1 Principal-Agent Problem 3 1.1 Problem Formulation 3 1.2 Further Reading 6 References
More informationSolving real-life portfolio problem using stochastic programming and Monte-Carlo techniques
Solving real-life portfolio problem using stochastic programming and Monte-Carlo techniques 1 Introduction Martin Branda 1 Abstract. We deal with real-life portfolio problem with Value at Risk, transaction
More informationWe formulate and solve two new stochastic linear programming formulations of appointment scheduling
Published online ahead of print December 7, 2011 INFORMS Journal on Computing Articles in Advance, pp. 1 17 issn 1091-9856 eissn 1526-5528 http://dx.doi.org/10.1287/ijoc.1110.0482 2011 INFORMS Dynamic
More informationOutline. 1 Introduction. 2 Algorithms. 3 Examples. Algorithm 1 General coordinate minimization framework. 1: Choose x 0 R n and set k 0.
Outline Coordinate Minimization Daniel P. Robinson Department of Applied Mathematics and Statistics Johns Hopkins University November 27, 208 Introduction 2 Algorithms Cyclic order with exact minimization
More information6.231 DYNAMIC PROGRAMMING LECTURE 3 LECTURE OUTLINE
6.21 DYNAMIC PROGRAMMING LECTURE LECTURE OUTLINE Deterministic finite-state DP problems Backward shortest path algorithm Forward shortest path algorithm Shortest path examples Alternative shortest path
More informationEssays on Some Combinatorial Optimization Problems with Interval Data
Essays on Some Combinatorial Optimization Problems with Interval Data a thesis submitted to the department of industrial engineering and the institute of engineering and sciences of bilkent university
More informationOnline Appendix: Extensions
B Online Appendix: Extensions In this online appendix we demonstrate that many important variations of the exact cost-basis LUL framework remain tractable. In particular, dual problem instances corresponding
More informationSOLVING ROBUST SUPPLY CHAIN PROBLEMS
SOLVING ROBUST SUPPLY CHAIN PROBLEMS Daniel Bienstock Nuri Sercan Özbay Columbia University, New York November 13, 2005 Project with Lucent Technologies Optimize the inventory buffer levels in a complicated
More informationApproximations of Stochastic Programs. Scenario Tree Reduction and Construction
Approximations of Stochastic Programs. Scenario Tree Reduction and Construction W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany www.mathematik.hu-berlin.de/~romisch
More informationBirkbeck MSc/Phd Economics. Advanced Macroeconomics, Spring Lecture 2: The Consumption CAPM and the Equity Premium Puzzle
Birkbeck MSc/Phd Economics Advanced Macroeconomics, Spring 2006 Lecture 2: The Consumption CAPM and the Equity Premium Puzzle 1 Overview This lecture derives the consumption-based capital asset pricing
More informationOptimal prepayment of Dutch mortgages*
137 Statistica Neerlandica (2007) Vol. 61, nr. 1, pp. 137 155 Optimal prepayment of Dutch mortgages* Bart H. M. Kuijpers ABP Investments, P.O. Box 75753, NL-1118 ZX Schiphol, The Netherlands Peter C. Schotman
More informationStochastic Programming: introduction and examples
Stochastic Programming: introduction and examples Amina Lamghari COSMO Stochastic Mine Planning Laboratory Department of Mining and Materials Engineering Outline What is Stochastic Programming? Why should
More informationThe Optimization Process: An example of portfolio optimization
ISyE 6669: Deterministic Optimization The Optimization Process: An example of portfolio optimization Shabbir Ahmed Fall 2002 1 Introduction Optimization can be roughly defined as a quantitative approach
More informationRisk aversion in multi-stage stochastic programming: a modeling and algorithmic perspective
Risk aversion in multi-stage stochastic programming: a modeling and algorithmic perspective Tito Homem-de-Mello School of Business Universidad Adolfo Ibañez, Santiago, Chile Joint work with Bernardo Pagnoncelli
More informationApproximation of Continuous-State Scenario Processes in Multi-Stage Stochastic Optimization and its Applications
Approximation of Continuous-State Scenario Processes in Multi-Stage Stochastic Optimization and its Applications Anna Timonina University of Vienna, Abraham Wald PhD Program in Statistics and Operations
More informationMartingale Pricing Theory in Discrete-Time and Discrete-Space Models
IEOR E4707: Foundations of Financial Engineering c 206 by Martin Haugh Martingale Pricing Theory in Discrete-Time and Discrete-Space Models These notes develop the theory of martingale pricing in a discrete-time,
More informationDynamic Programming (DP) Massimo Paolucci University of Genova
Dynamic Programming (DP) Massimo Paolucci University of Genova DP cannot be applied to each kind of problem In particular, it is a solution method for problems defined over stages For each stage a subproblem
More informationOn Complexity of Multistage Stochastic Programs
On Complexity of Multistage Stochastic Programs Alexander Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0205, USA e-mail: ashapiro@isye.gatech.edu
More informationCHAPTER 5: DYNAMIC PROGRAMMING
CHAPTER 5: DYNAMIC PROGRAMMING Overview This chapter discusses dynamic programming, a method to solve optimization problems that involve a dynamical process. This is in contrast to our previous discussions
More informationRisk Management for Chemical Supply Chain Planning under Uncertainty
for Chemical Supply Chain Planning under Uncertainty Fengqi You and Ignacio E. Grossmann Dept. of Chemical Engineering, Carnegie Mellon University John M. Wassick The Dow Chemical Company Introduction
More informationInvestigation of the and minimum storage energy target levels approach. Final Report
Investigation of the AV@R and minimum storage energy target levels approach Final Report First activity of the technical cooperation between Georgia Institute of Technology and ONS - Operador Nacional
More informationAsset-Liability Management
Asset-Liability Management John Birge University of Chicago Booth School of Business JRBirge INFORMS San Francisco, Nov. 2014 1 Overview Portfolio optimization involves: Modeling Optimization Estimation
More informationLecture 17: More on Markov Decision Processes. Reinforcement learning
Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture
More informationDynamic Risk Management in Electricity Portfolio Optimization via Polyhedral Risk Functionals
Dynamic Risk Management in Electricity Portfolio Optimization via Polyhedral Risk Functionals A. Eichhorn and W. Römisch Humboldt-University Berlin, Department of Mathematics, Germany http://www.math.hu-berlin.de/~romisch
More informationApproximate Revenue Maximization with Multiple Items
Approximate Revenue Maximization with Multiple Items Nir Shabbat - 05305311 December 5, 2012 Introduction The paper I read is called Approximate Revenue Maximization with Multiple Items by Sergiu Hart
More informationDynamic sampling algorithms for multi-stage stochastic programs with risk aversion
Dynamic sampling algorithms for multi-stage stochastic programs with risk aversion A.B. Philpott y and V.L. de Matos z March 28, 2011 Abstract We consider the incorporation of a time-consistent coherent
More informationShape-Preserving Dynamic Programming
Shape-Preserving Dynamic Programming Kenneth Judd and Yongyang Cai July 20, 2011 1 Introduction The multi-stage decision-making problems are numerically challenging. When the problems are time-separable,
More informationMultistage Stochastic Mixed-Integer Programs for Optimizing Gas Contract and Scheduling Maintenance
Multistage Stochastic Mixed-Integer Programs for Optimizing Gas Contract and Scheduling Maintenance Zhe Liu Siqian Shen September 2, 2012 Abstract In this paper, we present multistage stochastic mixed-integer
More informationInteger Programming Models
Integer Programming Models Fabio Furini December 10, 2014 Integer Programming Models 1 Outline 1 Combinatorial Auctions 2 The Lockbox Problem 3 Constructing an Index Fund Integer Programming Models 2 Integer
More informationStochastic Programming in Gas Storage and Gas Portfolio Management. ÖGOR-Workshop, September 23rd, 2010 Dr. Georg Ostermaier
Stochastic Programming in Gas Storage and Gas Portfolio Management ÖGOR-Workshop, September 23rd, 2010 Dr. Georg Ostermaier Agenda Optimization tasks in gas storage and gas portfolio management Scenario
More informationDynamic Programming: An overview. 1 Preliminaries: The basic principle underlying dynamic programming
Dynamic Programming: An overview These notes summarize some key properties of the Dynamic Programming principle to optimize a function or cost that depends on an interval or stages. This plays a key role
More informationFinancial Optimization ISE 347/447. Lecture 15. Dr. Ted Ralphs
Financial Optimization ISE 347/447 Lecture 15 Dr. Ted Ralphs ISE 347/447 Lecture 15 1 Reading for This Lecture C&T Chapter 12 ISE 347/447 Lecture 15 2 Stock Market Indices A stock market index is a statistic
More informationDynamic sampling algorithms for multi-stage stochastic programs with risk aversion
Dynamic sampling algorithms for multi-stage stochastic programs with risk aversion A.B. Philpott y and V.L. de Matos z October 7, 2011 Abstract We consider the incorporation of a time-consistent coherent
More informationTHE TRAVELING SALESMAN PROBLEM FOR MOVING POINTS ON A LINE
THE TRAVELING SALESMAN PROBLEM FOR MOVING POINTS ON A LINE GÜNTER ROTE Abstract. A salesperson wants to visit each of n objects that move on a line at given constant speeds in the shortest possible time,
More informationEE365: Risk Averse Control
EE365: Risk Averse Control Risk averse optimization Exponential risk aversion Risk averse control 1 Outline Risk averse optimization Exponential risk aversion Risk averse control Risk averse optimization
More informationInformation Relaxations and Duality in Stochastic Dynamic Programs
Information Relaxations and Duality in Stochastic Dynamic Programs David Brown, Jim Smith, and Peng Sun Fuqua School of Business Duke University February 28 1/39 Dynamic programming is widely applicable
More informationMulti-stage Stochastic Linear Programming: An Approach by Events
Multi-stage Stochastic Linear Programming: An Approach by Events C. Beltran-Royo L. F. Escudero R. E. Rodriguez-Ravines 28/07/2008 Abstract To solve the multi-stage linear programming problem, one may
More informationApproximate Dynamic Programming for the Merchant Operations of Commodity and Energy Conversion Assets
Approximate Dynamic Programming for the Merchant Operations of Commodity and Energy Conversion Assets Selvaprabu (Selva) Nadarajah, (Joint work with François Margot and Nicola Secomandi) Tepper School
More informationOptimal Security Liquidation Algorithms
Optimal Security Liquidation Algorithms Sergiy Butenko Department of Industrial Engineering, Texas A&M University, College Station, TX 77843-3131, USA Alexander Golodnikov Glushkov Institute of Cybernetics,
More informationStochastic Optimal Control
Stochastic Optimal Control Lecturer: Eilyan Bitar, Cornell ECE Scribe: Kevin Kircher, Cornell MAE These notes summarize some of the material from ECE 5555 (Stochastic Systems) at Cornell in the fall of
More informationSTOCHASTIC PROGRAMMING FOR ASSET ALLOCATION IN PENSION FUNDS
STOCHASTIC PROGRAMMING FOR ASSET ALLOCATION IN PENSION FUNDS IEGOR RUDNYTSKYI JOINT WORK WITH JOËL WAGNER > city date
More informationOPTIMIZATION METHODS IN FINANCE
OPTIMIZATION METHODS IN FINANCE GERARD CORNUEJOLS Carnegie Mellon University REHA TUTUNCU Goldman Sachs Asset Management CAMBRIDGE UNIVERSITY PRESS Foreword page xi Introduction 1 1.1 Optimization problems
More informationDynamic Asset and Liability Management Models for Pension Systems
Dynamic Asset and Liability Management Models for Pension Systems The Comparison between Multi-period Stochastic Programming Model and Stochastic Control Model Muneki Kawaguchi and Norio Hibiki June 1,
More informationWhere Has All the Value Gone? Portfolio risk optimization using CVaR
Where Has All the Value Gone? Portfolio risk optimization using CVaR Jonathan Sterbanz April 27, 2005 1 Introduction Corporate securities are widely used as a means to boost the value of asset portfolios;
More information1 Precautionary Savings: Prudence and Borrowing Constraints
1 Precautionary Savings: Prudence and Borrowing Constraints In this section we study conditions under which savings react to changes in income uncertainty. Recall that in the PIH, when you abstract from
More informationForecast Horizons for Production Planning with Stochastic Demand
Forecast Horizons for Production Planning with Stochastic Demand Alfredo Garcia and Robert L. Smith Department of Industrial and Operations Engineering Universityof Michigan, Ann Arbor MI 48109 December
More informationOPTIMIZATION MODELING FOR TRADEOFF ANALYSIS OF HIGHWAY INVESTMENT ALTERNATIVES
IIT Networks and Optimization Seminar OPTIMIZATION MODEING FOR TRADEOFF ANAYSIS OF HIGHWAY INVESTMENT ATERNATIVES Dr. Zongzhi i, Assistant Professor Dept. of Civil, Architectural and Environmental Engineering
More informationOptimal Production-Inventory Policy under Energy Buy-Back Program
The inth International Symposium on Operations Research and Its Applications (ISORA 10) Chengdu-Jiuzhaigou, China, August 19 23, 2010 Copyright 2010 ORSC & APORC, pp. 526 532 Optimal Production-Inventory
More informationPricing Problems under the Markov Chain Choice Model
Pricing Problems under the Markov Chain Choice Model James Dong School of Operations Research and Information Engineering, Cornell University, Ithaca, New York 14853, USA jd748@cornell.edu A. Serdar Simsek
More informationStochastic Approximation Algorithms and Applications
Harold J. Kushner G. George Yin Stochastic Approximation Algorithms and Applications With 24 Figures Springer Contents Preface and Introduction xiii 1 Introduction: Applications and Issues 1 1.0 Outline
More informationEARLY EXERCISE OPTIONS: UPPER BOUNDS
EARLY EXERCISE OPTIONS: UPPER BOUNDS LEIF B.G. ANDERSEN AND MARK BROADIE Abstract. In this article, we discuss how to generate upper bounds for American or Bermudan securities by Monte Carlo methods. These
More informationCasino gambling problem under probability weighting
Casino gambling problem under probability weighting Sang Hu National University of Singapore Mathematical Finance Colloquium University of Southern California Jan 25, 2016 Based on joint work with Xue
More informationOptimization in Finance
Research Reports on Mathematical and Computing Sciences Series B : Operations Research Department of Mathematical and Computing Sciences Tokyo Institute of Technology 2-12-1 Oh-Okayama, Meguro-ku, Tokyo
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationStochastic Programming and Financial Analysis IE447. Midterm Review. Dr. Ted Ralphs
Stochastic Programming and Financial Analysis IE447 Midterm Review Dr. Ted Ralphs IE447 Midterm Review 1 Forming a Mathematical Programming Model The general form of a mathematical programming model is:
More informationModeling Time-dependent Randomness in Stochastic Dual Dynamic Programming
Modeling Time-dependent Randomness in Stochastic Dual Dynamic Programming Nils Löhndorf Department of Information Systems and Operations Vienna University of Economics and Business Vienna, Austria nils.loehndorf@wu.ac.at
More informationThe Uncertain Volatility Model
The Uncertain Volatility Model Claude Martini, Antoine Jacquier July 14, 008 1 Black-Scholes and realised volatility What happens when a trader uses the Black-Scholes (BS in the sequel) formula to sell
More informationScenario tree generation for stochastic programming models using GAMS/SCENRED
Scenario tree generation for stochastic programming models using GAMS/SCENRED Holger Heitsch 1 and Steven Dirkse 2 1 Humboldt-University Berlin, Department of Mathematics, Germany 2 GAMS Development Corp.,
More informationMultistage Stochastic Programs
Multistage Stochastic Programs Basic Formulations Multistage Stochastic Linear Program with Recourse: all functions are linear in decision variables Problem of Private Investor Revisited Horizon and Stages
More informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
More information1 Shapley-Shubik Model
1 Shapley-Shubik Model There is a set of buyers B and a set of sellers S each selling one unit of a good (could be divisible or not). Let v ij 0 be the monetary value that buyer j B assigns to seller i
More informationarxiv: v1 [math.oc] 30 Mar 2018
Convex Hull, IP and European Electricity Pricing in a arxiv:1804.00048v1 [math.oc] 30 Mar 2018 European Power Exchanges setting with efficient computation of Convex Hull Prices Mehdi Madani 1, Carlos Ruiz
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More information1 Overview. 2 The Gradient Descent Algorithm. AM 221: Advanced Optimization Spring 2016
AM 22: Advanced Optimization Spring 206 Prof. Yaron Singer Lecture 9 February 24th Overview In the previous lecture we reviewed results from multivariate calculus in preparation for our journey into convex
More informationNotes. Cases on Static Optimization. Chapter 6 Algorithms Comparison: The Swing Case
Notes Chapter 2 Optimization Methods 1. Stationary points are those points where the partial derivatives of are zero. Chapter 3 Cases on Static Optimization 1. For the interested reader, we used a multivariate
More informationCS 188: Artificial Intelligence
CS 188: Artificial Intelligence Markov Decision Processes Dan Klein, Pieter Abbeel University of California, Berkeley Non-Deterministic Search 1 Example: Grid World A maze-like problem The agent lives
More informationI R TECHNICAL RESEARCH REPORT. A Framework for Mixed Estimation of Hidden Markov Models. by S. Dey, S. Marcus T.R
TECHNICAL RESEARCH REPORT A Framework for Mixed Estimation of Hidden Markov Models by S. Dey, S. Marcus T.R. 98-31 I R INSTITUTE FOR SYSTEMS RESEARCH ISR develops, applies and teaches advanced methodologies
More informationHandout 8: Introduction to Stochastic Dynamic Programming. 2 Examples of Stochastic Dynamic Programming Problems
SEEM 3470: Dynamic Optimization and Applications 2013 14 Second Term Handout 8: Introduction to Stochastic Dynamic Programming Instructor: Shiqian Ma March 10, 2014 Suggested Reading: Chapter 1 of Bertsekas,
More informationThe Irrevocable Multi-Armed Bandit Problem
The Irrevocable Multi-Armed Bandit Problem Ritesh Madan Qualcomm-Flarion Technologies May 27, 2009 Joint work with Vivek Farias (MIT) 2 Multi-Armed Bandit Problem n arms, where each arm i is a Markov Decision
More informationCOMP331/557. Chapter 6: Optimisation in Finance: Cash-Flow. (Cornuejols & Tütüncü, Chapter 3)
COMP331/557 Chapter 6: Optimisation in Finance: Cash-Flow (Cornuejols & Tütüncü, Chapter 3) 159 Cash-Flow Management Problem A company has the following net cash flow requirements (in 1000 s of ): Month
More informationRegret Minimization against Strategic Buyers
Regret Minimization against Strategic Buyers Mehryar Mohri Courant Institute & Google Research Andrés Muñoz Medina Google Research Motivation Online advertisement: revenue of modern search engine and
More informationDynamic - Cash Flow Based - Inventory Management
INFORMS Applied Probability Society Conference 2013 -Costa Rica Meeting Dynamic - Cash Flow Based - Inventory Management Michael N. Katehakis Rutgers University July 15, 2013 Talk based on joint work with
More informationLecture 14: Basic Fixpoint Theorems (cont.)
Lecture 14: Basic Fixpoint Theorems (cont) Predicate Transformers Monotonicity and Continuity Existence of Fixpoints Computing Fixpoints Fixpoint Characterization of CTL Operators 1 2 E M Clarke and E
More informationLecture 10: The knapsack problem
Optimization Methods in Finance (EPFL, Fall 2010) Lecture 10: The knapsack problem 24.11.2010 Lecturer: Prof. Friedrich Eisenbrand Scribe: Anu Harjula The knapsack problem The Knapsack problem is a problem
More informationROBUST OPTIMIZATION OF MULTI-PERIOD PRODUCTION PLANNING UNDER DEMAND UNCERTAINTY. A. Ben-Tal, B. Golany and M. Rozenblit
ROBUST OPTIMIZATION OF MULTI-PERIOD PRODUCTION PLANNING UNDER DEMAND UNCERTAINTY A. Ben-Tal, B. Golany and M. Rozenblit Faculty of Industrial Engineering and Management, Technion, Haifa 32000, Israel ABSTRACT
More informationComparative Study between Linear and Graphical Methods in Solving Optimization Problems
Comparative Study between Linear and Graphical Methods in Solving Optimization Problems Mona M Abd El-Kareem Abstract The main target of this paper is to establish a comparative study between the performance
More informationAll Investors are Risk-averse Expected Utility Maximizers
All Investors are Risk-averse Expected Utility Maximizers Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel) AFFI, Lyon, May 2013. Carole Bernard All Investors are
More informationInterpolation. 1 What is interpolation? 2 Why are we interested in this?
Interpolation 1 What is interpolation? For a certain function f (x we know only the values y 1 = f (x 1,,y n = f (x n For a point x different from x 1,,x n we would then like to approximate f ( x using
More informationAction Selection for MDPs: Anytime AO* vs. UCT
Action Selection for MDPs: Anytime AO* vs. UCT Blai Bonet 1 and Hector Geffner 2 1 Universidad Simón Boĺıvar 2 ICREA & Universitat Pompeu Fabra AAAI, Toronto, Canada, July 2012 Online MDP Planning and
More information6.231 DYNAMIC PROGRAMMING LECTURE 3 LECTURE OUTLINE
6.21 DYNAMIC PROGRAMMING LECTURE LECTURE OUTLINE Deterministic finite-state DP problems Backward shortest path algorithm Forward shortest path algorithm Shortest path examples Alternative shortest path
More informationCONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES
CONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES D. S. SILVESTROV, H. JÖNSSON, AND F. STENBERG Abstract. A general price process represented by a two-component
More information