870 likes | 1.11k Views
PARETO LINEAR PROGRAMMING. 3.4 Pareto Linear Programming. The Problem: P-opt Cx s.t Ax ≤ b x ≥ 0 where C is a kxn matrix so that Cx = (c (1) x, c (2) x, ..., c (k) x) where c (j) = jth row of C. . Example. z 1. z 2.
E N D
PARETO LINEAR PROGRAMMING
3.4 Pareto Linear Programming • The Problem: P-opt Cx s.t Ax ≤ b x ≥ 0 where C is a kxn matrix so that Cx = (c(1)x, c(2)x, ..., c(k)x) where c(j) = jth row of C.
Example z1 z2 • P- max {3x1 + 2x2, 2x1 – x2} • S.t. • x1 ≤ 2 • 3x1 + x2 ≤ 9 • x1, x2 ≥ 0 Here z = (z1, z2).
We refer to the Pareto solutions also as efficient points, or non-dominated points. • We shall focus on opt = max. • For simplicity we shall assume that the set of feasible solution, namely X := {x in Rn: Ax ≤ b, x ≥ 0} is bounded. • Let Z:= {Cx: x in X} observing that Z is a subset of Rk, and usually k <<<<<n.
GRAPHICAL INSIGHT
We can solve for the efficient points graphically if we only have 2 variables and 2 objective functions.This also gives us a good picture of what is going on, and ideas for generalizing. • P- max {3x1 + 2x2, 2x1 – x2} • S.t. • x1 ≤ 2 • 3x1 + x2 ≤ 9 • x1, x2 ≥ 0
So more generally:Geometry ZRk X Rn Cx x X Ax ≤ b x ≥ 0
Find X and Z for the problem below. • P- max {3x1 + 2x2, 2x1 – x2} • S.t. • x1 ≤ 2 • 3x1 + x2 ≤ 9 • x1, x2 ≥ 0
How do we generate all the Pareto Solutions? • Idea: • Generate the non-dominated extreme points of Z • Use them to generate the other efficient points of Z. • Rationale: • Other efficient points of Z are convex combinations of the efficient extreme points of Z
k=2 (i.e 2-objectives) z2 z1
k=2 z2 Z z1
k=2 Efficient points (Pareto solutions) z2 Z z1
k=2 Efficient points z2 Every efficient point can be expressed as a convex combination of the efficient extreme points. So we first aim to generate the efficient extreme points. Efficient extreme points Z z1
An extremely important theorem • Idea: • Theorem 4.6.1 • An extreme point of Z must be produced by an extreme point of X, ie if z* is an extreme point of Z then z*=Cx* where x* is an extreme point of X • Proof: • By contradiction. Assume that z* is an extreme point of Z but x* in X for which z*=Cx*, is not an extreme point of X.
Geometry z* ZRk X Rn x* Cx x X Ax ≤ b x ≥ 0
Since x* is not an extreme point of X we have x* = x’ + (1-)x” , 0 < < 1 for some points x’ and x” in X. Thus, z* = Cx* = C(x’ + (1-)x”) = Cx’ + (1-)Cx” = z’ + (1-)z” , (z’=Cx’, z”=Cx”) This, however, contradicts the assertion that z* is an extreme point of Z. (This is true only when z’ and z” are different and neither of z’ and z” is an extreme point)
Comment z* ZRk X Rn x* Cx x X Ax ≤ b x ≥ 0 If z* extreme then x* extreme. BUT not necessarily vice versa.
THEORY LEADING TO MORE GENERAL SOLUTIONS
One way to generate the efficient extreme points of Z is by deploying the following well known results: • Theorem 4.6.2: • If z* is an efficient point of Z then there must be a in Rk such that z* is an optimal solution to the problem: • max { z: z in Z} (****) i.e. max {1 z1+ 2 z2+…+ k zk : z in Z} • Furthermore, > 0 (all components are strictly positive)
Theorem 4.6.3: If > 0 then any optimal solution to max {1 z1+ 2 z2+…+ k zk : z in Z} (****) is an efficient point of Z.
Proof of Theorem 4.6.3: • Assume that > 0 and let z* be any optimal solution to (****). Contrary to the Theorem, assume that z* is not an efficient point of Z. This means that there exists some z’ in Z such that • z’j ≥ z*j , for all j=1,2,...,k • and • z’p > z*p , for some 1 ≤ p ≤ k. • (Remember that we mean Pareto-max here) • Since >0, this implies that z’ > z*, contradicting the assertion that z* is an optimal solution to (****).
Understanding Theorem 4.6.2 requires some fundamental results!!! • We will omit slides 17 - 40 (not examinable). • From slide 41 on is examinable.
Review(618-261) • Convex set: • If y’ and y” are elements of Y then the entire line segment connecting these points is also in Y. • Line segment: • Line segment connecting y’ and y” is the set of all the convex combinations of y’ and y”. • Convex combinations: • y =y’ + (1-)y” , 0 ≤ ≤ 1.
X2 Extreme Points Convex comb. of X1 and X2 X1 X2 X1 X2 • A point y in Y is an extreme point of Y if it cannot be expressed as a convex combination of two other points in Y. Not convex comb. of X1 and X2 X1
2 extreme points 4 extreme points (the corners) Entire boundary is extreme
Properties of Convex Sets • If S is a convex set, then for any in R+, S := {s: s in S} is a convex set. • If S and T are convex sets, then so is S + T : = {s+t: s in S, t in T} • The intersection of any collection of convex sets is convex.
S S A origin S+T S B T
facts about hyperplanes • Fact 1: • Let be a non-zero element of Rn and let b be a real number. Then the set H := {x in Rn: x = b} is a hyperplane in Rn. • Example: • In R3, the set of all points satisfying 3x1 + 3x2 + x3 =5 is a hyperplane (set b=5 and = (3,3,1) )
Fact 2: • Let H be a hyperplane in Rn. Then there is a non-zero vector in Rn and a constant b such that • H = S:= {x in Rn: x=b}.
bottom line • A hyperplane in Rn is the set of solutions to a single linear equation.
Half Spaces • Given a hyperplane, say H := {x in Rn : x = b} we shall consider the two closed half spaces it generates: H+ := {x in Rn: x ≥ b} H– := {x in Rn: x ≤ b}
H+ H– H := {x in Rn: x=b}
Main Results(for our purposes) • Theorem 4.6.4Separating Hyperplanes • Given a convex set S and a point x exterior to its closure*, there is a hyperplane containing x that contains S in one of its half spaces. • Closure of S: the smallest closed set containing S. • Closed set: A set with the property that any point that is arbitrarily close to it is a member of the set.
x H S
x S H x S
Theorem 4.6.5 Supporting Hyperplane • Let S be a convex set and let x be a boundary point of S. Then there is a hyperplane containing x and containing S in one of its closed half spaces. x S Supporting hyperplane at x
Theorem 4.6.6 • Let S be a convex set, H a supporting hyperplane of S and I the intersection of H and S. • Then every extreme point of I is an extreme point of S.
I One common extreme point S H Face two extreme points in common I S Facet H
And so ...... • For every extreme point in Z there is a supporting hyperplane • Each extreme point of Z is an optimal solution to max {z : z in Z} for some non zero in Rk.
And so ...... • For each efficient extreme point of Z there is a strictly positive in Rk such that the point is an optimal solution to max {z : z in Z}. • Each extreme point of Z is an optimal solution to max {z : z in Z} for some non zero in Rk. • And also 1 + + k = 1 (See supplementary notes for details).
Bottom line • We can generate the efficient extreme points associated with P-opt Cx Ax ≤ b x ≥ 0 by solving maxCx Ax ≤ b x≥ 0 for all > 0 .
k=2 • Thus we need two multipliers 1 and 2. • The objective function of the parametric linear programming problem will therefore be of the form: z() := 1c(1)x + 2c(2)x • But since the parameters are positive, we can divide say by 1, so obtain an equivalent objective function of the form z() := c(1)x + c(2)x , =2 / 1 The parametric problem is thus:
Efficient points z2 Efficient extreme points Z z1
z2 Z z1
z1 + z2 = Constant z2 Z z1
z1 + z2 = Constant z2 Z z1