25.11.2015 Views

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

Optimizer API for Python - Documentation - Mosek

Optimizer API for Python - Documentation - Mosek

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.

188 CHAPTER 15. SENSITIVITY ANALYSIS<br />

In summary, the basis type sensitivity analysis is computationally cheap but does not provide complete<br />

information. Hence, the results of the basis type sensitivity analysis should be used with care.<br />

15.4.3 <strong>The</strong> optimal partition type sensitivity analysis<br />

Another method for computing the complete linearity interval is called the optimal partition type<br />

sensitivity analysis. <strong>The</strong> main drawback of the optimal partition type sensitivity analysis is that it<br />

is computationally expensive compared to the basis type analysts. This type of sensitivity analysis is<br />

currently provided as an experimental feature in <strong>MOSEK</strong>.<br />

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<br />

objective value is given by<br />

<strong>The</strong> left and right shadow prices σ 1 and σ 2 for l c i<br />

z ∗ := c T x ∗ .<br />

are given by this pair of optimization problems:<br />

and<br />

σ 1 = minimize e T i s c l<br />

subject to A T (s c l − s c u) + s x l − s x u = c,<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) = z ∗ ,<br />

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

σ 2 = maximize e T i s c l<br />

subject to A T (s c l − s c u) + s x l − s x u = c,<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) = z ∗ ,<br />

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

<strong>The</strong>se 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) ∗ )<br />

is an arbitrary optimal solution then<br />

Next, the linearity interval [β 1 , β 2 ] for l c i<br />

(s c l ) ∗ i ∈ [σ 1 , σ 2 ].<br />

is computed by solving the two optimization problems<br />

and<br />

β 1 = minimize β<br />

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

c T x − σ 1 β = z ∗ ,<br />

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

β 2 = maximize β<br />

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

c T x − σ 2 β = z ∗ ,<br />

l x ≤ x ≤ u x .

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

Saved successfully!

Ooh no, something went wrong!