The MOSEK Python optimizer API manual Version 7.0 (Revision 141)

Optimizer API for Python - Documentation - Mosek Optimizer API for Python - Documentation - Mosek

25.11.2015 Views

188 CHAPTER 15. SENSITIVITY ANALYSIS In summary, the basis type sensitivity analysis is computationally cheap but does not provide complete information. Hence, the results of the basis type sensitivity analysis should be used with care. 15.4.3 The optimal partition type sensitivity analysis Another method for computing the complete linearity interval is called the optimal partition type sensitivity analysis. The main drawback of the optimal partition type sensitivity analysis is that it is computationally expensive compared to the basis type analysts. This type of sensitivity analysis is currently provided as an experimental feature in MOSEK. Given the optimal primal and dual solutions to (15.1), i.e. x ∗ and ((s c l )∗ , (s c u) ∗ , (s x l )∗ , (s x u) ∗ ) the optimal objective value is given by The left and right shadow prices σ 1 and σ 2 for l c i z ∗ := c T x ∗ . are given by this pair of optimization problems: and σ 1 = minimize e T i s c l subject to A T (s c l − s c u) + s x l − s x u = c, (l c ) T (s c l ) − (u c ) T (s c u) + (l x ) T (s x l ) − (u x ) T (s x u) = z ∗ , s c l , s c u, s c l , s x u ≥ 0 σ 2 = maximize e T i s c l subject to A T (s c l − s c u) + s x l − s x u = c, (l c ) T (s c l ) − (u c ) T (s c u) + (l x ) T (s x l ) − (u x ) T (s x u) = z ∗ , s c l , s c u, s c l , s x u ≥ 0. These two optimization problems make it easy to interpret the shadow price. Indeed, if ((s c l )∗ , (s c u) ∗ , (s x l )∗ , (s x u) ∗ ) is an arbitrary optimal solution then Next, the linearity interval [β 1 , β 2 ] for l c i (s c l ) ∗ i ∈ [σ 1 , σ 2 ]. is computed by solving the two optimization problems and β 1 = minimize β subject to l c + βe i ≤ Ax ≤ u c , c T x − σ 1 β = z ∗ , l x ≤ x ≤ u x , β 2 = maximize β subject to l c + βe i ≤ Ax ≤ u c , c T x − σ 2 β = z ∗ , l x ≤ x ≤ u x .

15.4. SENSITIVITY ANALYSIS FOR LINEAR PROBLEMS 189 The linearity intervals and shadow prices for u c i , lx j , and ux j are computed similarly to lc i . The left and right shadow prices for c j denoted σ 1 and σ 2 respectively are computed as follows: and σ 1 = minimize e T j x subject to l c + βe i ≤ Ax ≤ u c , c T x = z ∗ , l x ≤ x ≤ u x σ 2 = maximize e T j x subject to l c + βe i ≤ Ax ≤ u c , c T x = z ∗ , l x ≤ x ≤ u x . Once again the above two optimization problems make it easy to interpret the shadow prices. Indeed, if x ∗ is an arbitrary primal optimal solution, then x ∗ j ∈ [σ 1 , σ 2 ]. The linearity interval [β 1 , β 2 ] for a c j is computed as follows: and β 1 = minimize β subject to A T (s c l − s c u) + s x l − s x u = c + βe j , (l c ) T (s c l ) − (u c ) T (s c u) + (l x ) T (s x l ) − (u x ) T (s x u) − σ 1 β ≤ z ∗ , s c l , s c u, s c l , s x u ≥ 0 β 2 = maximize β subject to A T (s c l − s c u) + s x l − s x u = c + βe j , (l c ) T (s c l ) − (u c ) T (s c u) + (l x ) T (s x l ) − (u x ) T (s x u) − σ 2 β ≤ z ∗ , s c l , s c u, s c l , s x u ≥ 0. 15.4.4 Example: Sensitivity analysis As an example we will use the following transportation problem. Consider the problem of minimizing the transportation cost between a number of production plants and stores. Each plant supplies a number of goods and each store has a given demand that must be met. Supply, demand and cost of transportation per unit are shown in Figure 15.2. If we denote the number of transported goods from location i to location j by x ij , problem can be formulated as the linear optimization problem minimize subject to 1x 11 + 2x 12 + 5x 23 + 2x 24 + 1x 31 + 2x 33 + 1x 34

15.4. SENSITIVITY ANALYSIS FOR LINEAR PROBLEMS 189<br />

<strong>The</strong> linearity intervals and shadow prices for u c i , lx j , and ux j are computed similarly to lc i .<br />

<strong>The</strong> left and right shadow prices for c j denoted σ 1 and σ 2 respectively are computed as follows:<br />

and<br />

σ 1 = minimize e T j x<br />

subject to l c + βe i ≤ Ax ≤ u c ,<br />

c T x = z ∗ ,<br />

l x ≤ x ≤ u x<br />

σ 2 = maximize e T j x<br />

subject to l c + βe i ≤ Ax ≤ u c ,<br />

c T x = z ∗ ,<br />

l x ≤ x ≤ u x .<br />

Once again the above two optimization problems make it easy to interpret the shadow prices. Indeed,<br />

if x ∗ is an arbitrary primal optimal solution, then<br />

x ∗ j ∈ [σ 1 , σ 2 ].<br />

<strong>The</strong> linearity interval [β 1 , β 2 ] for a c j is computed as follows:<br />

and<br />

β 1 = minimize β<br />

subject to A T (s c l − s c u) + s x l − s x u = c + βe j ,<br />

(l c ) T (s c l ) − (u c ) T (s c u) + (l x ) T (s x l ) − (u x ) T (s x u) − σ 1 β ≤ z ∗ ,<br />

s c l , s c u, s c l , s x u ≥ 0<br />

β 2 = maximize β<br />

subject to A T (s c l − s c u) + s x l − s x u = c + βe j ,<br />

(l c ) T (s c l ) − (u c ) T (s c u) + (l x ) T (s x l ) − (u x ) T (s x u) − σ 2 β ≤ z ∗ ,<br />

s c l , s c u, s c l , s x u ≥ 0.<br />

15.4.4 Example: Sensitivity analysis<br />

As an example we will use the following transportation problem. Consider the problem of minimizing<br />

the transportation cost between a number of production plants and stores. Each plant supplies a<br />

number of goods and each store has a given demand that must be met. Supply, demand and cost of<br />

transportation per unit are shown in Figure 15.2. If we denote the number of transported goods from<br />

location i to location j by x ij , problem can be formulated as the linear optimization problem minimize<br />

subject to<br />

1x 11 + 2x 12 + 5x 23 + 2x 24 + 1x 31 + 2x 33 + 1x 34

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

Saved successfully!

Ooh no, something went wrong!