Title: Stochastic Optimization ESI 6912
1Stochastic OptimizationESI 6912
NOTES 2 FARMING EXAMPLE
- Instructor Prof. S. Uryasev
2Farming Example
Outline
- 1. Deterministic Setup of The Optimization
Problem - 2. Extensive Form of The Stochastic Program
- 3. Recourse Reformulation of The Stochastic
Problem. - 4. Expected Value of Perfect Information (EVPI)
- 5. Value of Stochastic Solution
3Initial Data
4Variables
- acres of land devoted to wheat - acres of land
devoted to corn - acres of land devoted to sugar
beets - tons of wheat sold - tons of corn sold -
tons of sugar beets sold at favorable price -
tons of sugar beets sold at lower price - tons
of wheat purchased - tons of corn purchased
5Deterministic Setup of Optimization Problem
6Optimal Solution for Deterministic Optimization
problem
7Scenario Analysis
Optimal solution based on expected yields
Scenario 2
8Scenario Analysis (contd)
Optimal solution based on above average yields
(20)
Scenario 1
Optimal solution based on below average yields
(-20)
Scenario 3
9Extensive Form of Stochastic Problem
Scenario approach
Yield (W, C, SB)
the second stage -
1st scenario
20
(3.0, 3.6, 24)
probability
1/3
2nd scenario
1/3
0
(2.5, 3.0, 20)
1/3
3d scenario
-20
(2.0, 2.4, 16)
Decision
- the first stage
10Extensive Form of Stochastic Problem(contd)
1st scenario
2nd scenario
3d scenario
1st scenario
2nd scenario
3d scenario
11Extensive Form of Stochastic Problem(contd)
Optimal Solution
12Recourse Reformulation of Stochastic Program
- s-th scenario
Random matrix
13Recourse Reformulation of Stochastic Program
(contd)
deterministic part (1st stage)
stochastic part (2nd stage)
deterministic constraints (1st stage)
stochastic constraints (2nd stage)
14Recourse Reformulation of Stochastic Program
(contd)
Recourse function
Expected value of the recourse function
15Recourse Reformulation of Stochastic Program
(contd)
General model formulation
16Uncertain variables with continuous distributions
1. - distributed independently
2. -
uniformly distributed
density
17Extensive Formulation of the Stochastic Program
deterministic part (1st stage)
stochastic part (2nd stage)
deterministic constraints (1st stage)
stochastic constraints (2nd stage)
18Decomposition of Stochastic Program
sugar beets
corn
wheat
depend only upon decision and random yield
(wheat)
depend only upon decision and random yield
(corn)
depend only upon decision and random yield
(sugar beets)
19Recourse Functions
Wheat
Corn
Sugar beets
20Explicit Form for Recourse Functions
Wheat
Corn
Sugar beets
21Recourse Formulation of Stochastic Program
22Calculation of Expected Values for Recourse
Functions
Wheat
yield is uniformly distributed
1. In case when
where is the expected value
of
23Calculation of Expected Values for Recourse
Functions (contd)
2. In the case when
3. In the case when
24Calculation of Expected Values for Recourse
Functions (contd)
Wheat
Corn
Sugar beets
25Expected Recourse Value for Wheat as a Function
of Acres Planted
26Global Formulation of Stochastic Program
Convex optimization problem
are continuous convex functions depending only
upon decision vector
27Derivation of Optimal Solution
Notations
- the dual variable
Karush-Kuhn-Tucker conditions
28Calculation of Derivatives
Wheat
Corn
Sugar beets
29Optimal solution
Assume
Using enumerative technique, it can be
established that optimal solution should satisfy
System
Optimal values
30Deterministic Optimization
random variable
expected value of random variable
31Scenario Analysis
Optimization problem corresponding to scenario
realization of random variable in the
scenario
vector in the scenario
32Extensive Form of Stochastic Program
probability of the scenario
expected loss
33Variables
- acres of land devoted to wheat - acres of land
devoted to corn - acres of land devoted to sugar
beets - tons of wheat sold - tons of corn sold -
tons of sugar beets sold at favorable price -
tons of sugar beets sold at lower price - tons
of wheat purchased - tons of corn purchased
34Deterministic Setup of Optimization Problem
f(x,y)
35Optimal Solution for Deterministic Optimization
problem
36Scenario Analysis
Optimal solution based on expected yields
Scenario 2
37Scenario Analysis (contd)
Optimal solution based on above average yields
(20)
Scenario 1
Optimal solution based on below average yields
(-20)
Scenario 3
38Extensive Form of Stochastic Problem
Scenario approach
Yield (W, C, SB)
the second stage -
1st scenario
20
(3.0, 3.6, 24)
probability
1/3
2nd scenario
1/3
0
(2.5, 3.0, 20)
1/3
3d scenario
-20
(2.0, 2.4, 16)
Decision
- the first stage
39Extensive Form of Stochastic Problem(contd)
f(x,y1)
1st scenario
f(x,y2)
2nd scenario
f(x,y3)
3d scenario
1st scenario
2nd scenario
3d scenario
40Optimal Solution for Stochastic
ProblemFormulated in Extensive Form
41Problem with Recourse
Recourse function
42Problem with Recourse (contd)
where
43Recourse Reformulation of Stochastic Program
- s-th scenario
Random matrix
44Expected Value of Perfect Information (EVPI)
Solution with perfect information
is an optimal solution of
expected performance with perfect
information
Stochastic programming solution
is an optimal solution of
EVPI
45Value of Stochastic Solution (VSS)
Expected value solution
is an optimal solution of
Stochastic programming solution
is an optimal solution of
VSS
46Stochastic Programming Solution
47Scenario Analysis
Optimal solution based on expected yields
Scenario 2
48Scenario Analysis (contd)
Optimal solution based on above average yields
(20)
Scenario 1
Optimal solution based on below average yields
(-20)
Scenario 3
49Farming Example, EVPI
Solution with perfect information
Stochastic programming solution
50Farming Example, VSS
Expected value solution
is an optimal solution of
Stochastic programming solution
VSS