04.02.2015 Views

Stochastic Programming - Index of

Stochastic Programming - Index of

Stochastic Programming - Index of

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

NETWORK PROBLEMS 295<br />

0<br />

2<br />

1<br />

1<br />

(0,4)<br />

(0,3)<br />

2<br />

(0,3)<br />

1<br />

3<br />

2<br />

3<br />

3<br />

4 6<br />

4<br />

2<br />

(0,3)<br />

(0,6)<br />

8<br />

1<br />

(0,4)<br />

4 7 5<br />

5 -1<br />

2<br />

1<br />

(0,7)<br />

2<br />

(0,1)<br />

Slack<br />

3<br />

Figure 12 Example network with arc capacities and external flows<br />

corresponding to the Jensen lower bound.<br />

number <strong>of</strong> points. If there are k random variables, we must work with 2 k<br />

points. This means that with more than about 10 random variables we are<br />

not in business. Thus, since there are 11 random variable in the example, we<br />

must solve 2 11 problems to find the upper bound. We have not done that here.<br />

In what follows, we shall demonstrate how to obtain a piecewise linear upper<br />

bound that does not exhibit this exponential characterization. A weakness <strong>of</strong><br />

this bound is that it may be +∞ even if the problem is feasible. That may<br />

not happen to the Edmundson–Madansky upper bound. We shall continue to<br />

use the network in Figure 11 to illustrate the ideas.<br />

6.5.1 Piecewise Linear Upper Bounds<br />

Let us illustrate the method in a simplified setting. Define φ(ξ,η) by<br />

φ(ξ,η) =min{q T y | W ′ y = b + ξ, 0 ≤ y ≤ c + η},<br />

y<br />

where all elements <strong>of</strong> the random vectors ˜ξ = (˜ξ 1 T , ˜ξ 2 T ,...) T and ˜η =<br />

(˜η 1 T, ˜ηT 2 ,...)T are mutually independent. Furthermore, let the supports be<br />

given by Ξ(˜ξ) =[A, B] andΞ(˜η) =[0,C]. The matrix W ′ is the node–arc<br />

incidence matrix for a network, with one row removed. That row represents<br />

the slack node. The external flow in the slack node equals the negative sum <strong>of</strong><br />

the external flows in the other nodes. The goal is to create an upper bounding<br />

function U(ξ,η) that is piecewise linear, separable and convex in ξ, aswellas<br />

easily integrable in η:<br />

U(ξ,η) =φ(E ˜ξ,0) + H(η)+ ∑ i<br />

{<br />

d<br />

+<br />

i (ξ i − E ˜ξ i ) if ξ i ≥ E ˜ξ i ,<br />

d − i (E ˜ξ i − ξ i ) if ξ i

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!