Bonnel 2006
Bonnel 2006
DOI: 10.1007/s10957-006-9150-4
Communicated by H. P. Benson
Published Online: 5 December 2006
1. Introduction
365
0022-3239/06/1200-0365/0
C 2006 Springer Science+Business Media, Inc.
366 JOTA: VOL. 131, NO. 3, DECEMBER 2006
with f1 : Y → R, F1 : Y → Z, S1 ⊆ Y .
This problem is an important tool in the decision-making theory. The main
difficulty to solve Problem (OES) is due to its nonconvexity. Even when all the
functions are linear and S1 is a polyhedron, Problem (OES) is not convex. Another
difficulty is given by the fact that 1 , the feasible set of Problem (OES), is not
known explicitly. For these reasons, Problem (OES) was studied intensively during
the last two decades by many authors. Thus, the all linear case (i.e., when 1 is
a linear functional and the vector optimization problem has linear objectives and
linear constraints) has been studied theoretically and some algorithms have been
proposed in Refs. 1–5. Some nonlinear cases dealing also with theoretical aspects
and algorithms can be found in Refs. 6–13. For a survey paper and an extensive
bibliography, see Ref. 14.
Problem (OES) is obtained as a particular case of Problem (BLw ) considering
X = {0}, U = X × Y , and for every y ∈ Y ,
f (0, y) = f1 (y), F (0, y) = F1 (y), S(0) = S1 .
with
x ∈ Rnx , y ∈ Rny , f, F : Rnx +ny → R, G : Rnx +ny → Rnu ,
H : Rnx +ny → Rnl .
We obtain Problem (SBP) from Problem (BLw ) considering X = Rnx , Y =
R , Z = R, C = R+ ,
ny
2. Preliminaries
Let
·
denote the norm of Z and let ·, · denote the duality scalar product
between Z ∗ (the topological dual of Z) and Z. The extended space Z̄ = Z ∪
{−∞C , +∞C } has been introduced in Ref. 21. Recall that a neighborhood of +∞C
is a set N ⊆ Z̄ containing a + C ∪ {+∞C }, for some a ∈ Z, and its opposite −N
is a neighborhood of −∞C . The partial-order relation, compatible with the linear
structure of Z, given by
∀z, z ∈ Z, z ≤C z ⇐⇒ z − z ∈ C,
and the transitive relation
∀z, z ∈ Z, z <C z ⇐⇒ z − z ∈ int C,
are extended respectively to Z̄ by
∀z ∈ Z, −∞C <C z <C +∞C , −∞C ≤C z ≤C +∞C .
Note that the embedding Z ⊆ Z̄ is continuous and dense.
Let G be an extended-valued map from Y to Z ∪ {+∞C }, which is proper
(i.e., not identically equal to +∞C ) and let S1 ⊆ Y be a nonempty set.
The effective domain of G is
dom (G) = {y ∈ Y | G(y) = +∞C },
368 JOTA: VOL. 131, NO. 3, DECEMBER 2006
a point a ∈ Y is called:
(i) weakly-efficient (or weakly-Pareto) solution if a ∈ S1 and there is no
y ∈ S1 verifying
G(y) <C G(a) [in other words, G(a) ∈ w-INFC (G(S1 ))];
(ii) efficient (or Pareto) solution if a ∈ S1 and there is no y ∈ S1 verifying
G(y) ≤C G(a) and G(y) = G(a).
The set of weakly-efficient solutions is denoted by w-ARGMINC {G(y) | y ∈
S1 }. To solve the vector optimization problem
(VOP) MINIMIZE C G(y), s.t. y ∈ Y,
y
From now on, we suppose that the ordering cone C is also closed. For any
→ B, where A and B are some sets, we denote by
set-valued map : A →
Gr() := {(a, b) ∈ A × B | b ∈ (a)}
the graph of and by
dom() := {a ∈ A | (a) = ∅}
the domain of .
Note first that Problem (BLW ) is equivalent to the problem
(GMW ) min f (x, y), s.t. (x, y) ∈ Gr (w ) ∩ U,
x,y
In this section, we use a penalty function (inspired from Ref. 9) which allows
us to describe the graph of the set-valued map x → w (x) as an equality constraint.
Note that Reference 9 deals with the problem of optimizing a scalar function over
the weakly-efficient set [Problem OES, see Introduction] in finite-dimensional
spaces, while the present paper deals with the much more general problem of
semivectorial bilevel optimization, where the function F depends also on the
variable x of the upper level. Moreover, in this paper, we work in non-necessarily
finite-dimensional spaces.
In the sequel, e is a fixed vector belonging to int C. We consider the function
h : X × Y → R ∪ {+∞} given, for each (x, y) ∈ X × Y , by
h(x, y) = sup{t ∈ R | ∃y ∈ S(x) : F (x, y ) ≤C F (x, y) − te}. (1)
370 JOTA: VOL. 131, NO. 3, DECEMBER 2006
Proof. Let (x, y) ∈ Gr(S). It follows that (0, y) is feasible for the maxi-
mization problem defining the function h; hence, (i) holds.
Let (x, y) ∈ Gr(W ). We have that y ∈ W (x); i.e., y is a weakly-efficient
solution to Problem (LL)(x). It follows that y ∈ S(x); i.e., (x, y) ∈ Gr(S). From
(i), we have that h(x, y) ≥ 0. If h(x, y) = 0, then we can find y ∈ S(x) and a real
number t > 0 such that F (x, y ) ≤C F (x, y) − te, hence F (x, y ) <C F (x, y).
This is in contradiction to the fact that y is a weakly-efficient solution for Problem
(LL)(x); therefore, h(x, y) = 0.
Conversely, let (x, y) ∈ Gr(S) such that h(x, y) = 0. If (x, y) ∈/ Gr(W ), then
y is not a weakly-efficient solution for Problem (LL)(x); hence, there exists some
y ∈ S(x) such that F (x, y ) <C F (x, y). Thus, using the continuity of the map
t → t · e, there exists some t > 0 such that F (x, y ) ≤C F (x, y) − te. Therefore,
h(x, y) > 0 and we get a contradiction.
Using Proposition 2.1 and the definition of h, we obtain easily the following
remark.
Theorem 3.1. Let Z be endowed with the norm (strong) topology. Let X
and Y be Hausdorff topological spaces. Assume the following:
Proof. We will use the following classical result (see Ref. 23, Proposition
3.2; see also Ref. 24).
3 This
means that:
(a) S is upper semicontinuous at each x0 ∈ X; i.e., for any neighborhood N of S(x0 ), there exists a
neighborhood N of x0 such that
x ∈ N =⇒ S(x) ⊆ N.
(b) S is lower semicontinuons at each x0 ∈ X; i.e., for any y0 ∈ S(x0 ) and any neighborhood N of
y0 , there exists a neighborhood N of x0 such that
x ∈ N =⇒ S(x) ∩ N = ∅.
4 According to (A1), we have in fact that S(x) is compact.
5 Of course, α and K may depend on u0 .
372 JOTA: VOL. 131, NO. 3, DECEMBER 2006
Let
u = (x, y) ∈ Vu0 .
Since
lev0 ϕ(·, u) = {(t, y ) ∈ R × Y | y ∈ S(x), t ≥ 0, F (x, y ) ≤C F (x, y) − te},
and since (x, y) ∈ Gr(S), we have that (0, y) ∈ lev0 ϕ(·, u); hence, lev0 ϕ(·, u) is not
empty. Let (t, y ) ∈ lev0ϕ (·, u). For any λ ∈ C + ,
λ
∗ = 1(
·
∗ denotes the norm
of the dual space Z ∗ of the Banach space Z), we have λ, e ≥ d(e, Z\C) > 0
[see Ref. 21, Lemma 2.2]; hence,
t ≤ λ, F (x, y) − F (x, y )/λ, e ≤
F (x, y) − F (x, y )
/λ, e ≤ β. (4)
Since y ∈ S(x) ⊆ KYx0 , we obtain that (t, y ) ∈ K; hence, (3) holds.
To verify (ii), remark first that M(u0 ) is not empty because ϕ(·, u0 ) is contin-
uous and lev0 ϕ(·, u0 ) is compact.
Let NX̃ be a neighborhood of M(u0 ) and let (t0 , y0 ) ∈ M(u0 ). There exist
δ > 0 and a neighborhood Vy0 of y0 in Y such that
[t0 − δ, t0 + δ] × Vy0 ⊆ NX̄ . (5)
Since e ∈ int C, there exists r > 0 such that e − (1/r)BZ ⊂ C, where BZ = −BZ
is the open unit ball of Z. We can write BZ ⊂ re − C; i.e.,
∀h0 ∈ Z,
h0
< 1 =⇒ h0 ≤C re. (6)
Let ε = δ/2r. The continuity of F [restricted to Gr(S)] at (x0 , y0 ) and (x0 , y0 ),
the
lower semicontinuity of S at x0 , and the fact that y0 ∈ S(x0 ) ∩ Vy0 = ∅ ensure the
existence of a neighborhood Vx0 of x0 and of a neighborhood Vy0 of y0 such that:
(a) For any x ∈ Vx0 , we have S(x) ∩ Vy0 = ∅.
(b) For any (x, y) ∈ (Vx0 × Vy0 ) ∩ Gr(S), we have
F (x, y ) − F (x0 , y0 )
<
ε.
(c) For any (x, y ) ∈ (Vx0 × Vy0 ) ∩ Gr(S), we have
F (x, y ) − F (x0 , y0 )
<
ε.
Take
NV = (Vx0 × Vy0 ) ∩ Gr(S).
Let (x, y) ∈ NV and consider according to (a) an element ŷ0 ∈ S(x) ∩ Vy0 . From
(c), there exists h1 ∈ BZ such that
F (x, ŷ0 ) = F (x0 , y0 ) + εh1 .
From (b), there exists h2 ∈ BZ such that
F (x, y) = F (x0 , y0 ) + εh2 .
374 JOTA: VOL. 131, NO. 3, DECEMBER 2006
Theorem 3.2. Let Z be endowed with the norm (strong) topology. Let X
and Y be Hausdorff topological spaces. Assume the following:
(B1) The set Gr(S) ∩ U is compact.
(B2) The restriction of the map F to Gr(Ŝ) is continuous.
JOTA: VOL. 131, NO. 3, DECEMBER 2006 375
Proof. It is obvious from (B1) that the topological spaces X̂ and Ŷ are
compact spaces. Now, we apply Theorem 3.1 replacing the spaces X and Y by X̂
and Ŷ respectively. Consequently, F is restricted to X̂ × Ŷ and S is replaced by
Ŝ : X̂ →→ Ŷ . Note that (A1) of Theorem 3.1 is satisfied taking Nx = Kx = X̂ and
NYx = KYx = Ŷ .
Therefore, the restriction of h to Gr(Ŝ) is continuous. On the other hand,
Gr(Ŝ) ⊇ Gr(Ŝ) ∩ U and Gr(Ŝ) ∩ U = Gr(S) ∩ U according to the definition of
X̂, Ŷ , Ŝ.
Next, for each positive real number r, we consider the following penalized
problem:
(Pr ) min (f (x, y) + rh(x, y)) , s.t. (x, y) ∈ U ∩ Gr(S),
x,y
Theorem 3.3. Let X and Y be two reflexive Banach spaces. Assume that
X ⊆ X and Y ⊆ Y are weakly closed subsets. We consider X (resp. Y ) as a
Hausdorff topological space endowed with the trace of the weak topology of X
(resp. Y). We make all the assumptions of Theorem 3.1. Moreover, we suppose
that:
(A) The set Gr(S) ∩ U is closed in X × Y (i.e., with respect to the product
topology).
(B) The restriction of f to Gr(S) ∩ U is continuous (i.e., it is continuous
on the topological subspace Gr(S) ∩ U of the topological space X × Y
endowed with the product topology).
376 JOTA: VOL. 131, NO. 3, DECEMBER 2006
Proof. We adapt the classical proof (see for example Ref. 26, pp. 366–
369) to our special problem in not-necessarily finite-dimensional spaces. Problem
(BLw ) is equivalent to the problem of minimizing f over the set
A := Gr(W ) ∩ U,
which has been assumed to be nonempty in the Introduction. Let us denote
B = Gr(S) ∩ U.
Obviously, A ⊆ B; hence, B is nonempty. From Theorem 3.1, we have that the
restriction hB of h to B is continous. Moreover, according to Lemma 3.1, the
function h is positive over B. Also, we assumed that the restriction fB of f to B
is continuous. Thus, Problems (Pr ) and (BLW ) can be written equivalently as
min(fB (·, ·) + rhB (·, ·)), over B, (7)
respectively
min fB (·, ·), over A. (8)
Note that A = h−1 B ({0}); hence, A is closed in B; therefore, A is closed in X × Y
because B is closed in X × Y .
Let (x̂, ŷ) ∈ A and consider the sets
 = {(x, y) ∈ A | fB (x, y) ≤ fB (x̂, ŷ)},
B̂ = {(x, y) ∈ B | fB (x, y) ≤ fB (x̂, ŷ)}.
Obviously  and B̂ are nonempty (because (x̂, ŷ) ∈  ⊆ B̂) weakly closed sets
in X × Y (because they are closed in B) and bounded in X × Y according to
(CH). Thus, Â and B̂ are weakly compact sets in X × Y, hence compact sets in
the topological space X × Y . Let us denote
Lr = {(x, y) ∈ B | fB (x, y) + rhB (x, y) ≤ fB (x̂, ŷ)}.
Since hB ≥ 0, it follows that Lr ⊆ B̂. Hence, Lr is compact in X × Y and
nonempty [it contains (x̂, ŷ) because hB (x̂, ŷ) = 0]. Thus, for each r > 0, the
problem
min(fB (·, ·) + rhB (·, ·)), over Lr ,
admits an optimal solution (xr , yr ) which is obviously an optimal solution for
Problem (Pr ).
Also, there exists minimizers (x ∗ , y ∗ ) ∈ Â and (x̄, ȳ) ∈ B̂ of fB over the
compact sets  and B̂ respectively. It is obvious that (x ∗ , y ∗ ) solves Problem
(BLw ) and that
f∗ := f (x̄, ȳ) = min f (x, y).
(x,y)∈B
Since, for each r > 0, (xr , yr ) ∈ B̂, we obtain that the generalized sequence
(xr , yr )r↑∞ is bounded in X × Y.
For each r > 0, we have
f∗ ≤ f (xr , yr ) ≤ f (xr , yr ) + rh(xr , yr ) ≤ f (x ∗ , y ∗ );
hence,
0 ≤ h(xr , yr ) ≤ [f (x ∗ , y ∗ ) − f∗ ]/ r.
Thus, we obtain that h(xr , yr ) tends to 0 as r → +∞. Let (x̃, ỹ) be a weak cluster
point of the generalized sequence (xr , yr )r↑∞ in X × Y, hence (x̃, ỹ) ∈ X × Y .
Let (xrn , yrn )n∈N , with rn → +∞, be a sequence convergent to (x̃, ỹ) in X × Y
(of course, there exists such a sequence). By the continuity of hB , we obtain that
hB (x̃, ỹ) = 0, hence (x̃, ỹ) ∈ A. Also, for each (x, y) ∈ A, n ∈ N, we have
f (xrn , yrn ) ≤ f (xrn , yrn ) + rn h(xrn , yrn ) ≤ f (x, y);
hence, letting n → ∞, we obtain f (x̃, ỹ) ≤ f (x, y).
Theorem 3.4. Assume that the hypotheses of Theorem 3.2 hold and that
the restriction of f to Gr(S) ∩ U is continuous. Then:
(i) For every r > 0, Problem (Pr ) has at least a solution (xr , yr ).
378 JOTA: VOL. 131, NO. 3, DECEMBER 2006
(ii) Any cluster point (which exists) of the generalized sequence (xr , yr )r↑∞
solves Problem (BLw ).
Proof. Since now we assume that Gr(S) ∩ U is compact, one can adapt
easily and simplify the proof of Theorem 3.3. For the sake of clarity, we give the
proof.
Problem (BLW ) is equivalent to the problem of minimizing f over the set
A := Gr(W ) ∩ U , which has been assumed to be nonempty in the Introduction.
Let us denote
B = Gr(S) ∩ U.
Obviously A ⊆ B, hence B is nonempty. From Theorem 3.2, we have that the
restriction hB of h to B is continuous. Moreover, according to Lemma 3.1, the
function h is positive over B. Also, we assumed that the restriction fB of f to B
is continuous. Thus, Problems (Pr ) and (BLW ) can be written equivalently as
min(fB (·, ·) + rhB (·, ·)), over B, (9)
respectively
min fB (·, ·), over A. (10)
Note that A = h−1
B ({0}),
hence A is closed in B, therefore A is compact.
For each r > 0, Problem (Pr ) admits an optimal solution (xr , yr ) ∈ B. Also,
there exists minimizers (x ∗ , y ∗ ) ∈ A and (x̄, ȳ) ∈ B of fB over the compact sets
A and B respectively. It is obvious that (x ∗ , y ∗ ) solves Problem (BLW ) and we
denote f∗ := f (x̄, ȳ).
For each r > 0, we have
f∗ ≤ f (xr , yr ) ≤ f (xr , yr ) + rh(xr , yr ) ≤ f (x ∗ , y ∗ );
hence,
0 ≤ h(xr , yr ) ≤ [f (x ∗ , y ∗ ) − f∗ ]/r.
Thus, we obtain that h(xr , yr ) tends to 0 as r → +∞. Let (x̃, ỹ) be a cluster point
of the generalized sequence (xr , yr )r↑∞ ; let (xrn , yrn )n∈N , with rn → +∞, be a
sequence convergent to (x̃, ỹ) (of course there exists such a sequence). By the
continuity of hB , we obtain that hB (x̃, ỹ) = 0, hence (x̃, ỹ) ∈ A. Also, for each
(x, y) ∈ A, n ∈ N, we have that
f (xrn , yrn ) ≤ f (xrn , yrn ) + rn h(xrn , yrn ) ≤ f (x, y);
hence, letting n → ∞, we obtain
f (x̃, ỹ) ≤ f (x, y).
JOTA: VOL. 131, NO. 3, DECEMBER 2006 379
In the next remark, we show that, assuming X and Y are Euclidean spaces,
some hypothesis may be simplified.
Remark 3.2. In Theorem 3.3, let us assume that X and Y are Euclidean
spaces. In this case, the weak topology coincides with the strong one; thus, we have
to consider X and Y as closed subsets of X and Y. Also in this case, Assumption
(A1) of Theorem 3.1 may be simplified replacing it with the following:
(A1 ) for any
x ∈ X, there exists a bounded neighborhood Nx of x such that
the set x ∈Nx S(x ) is bounded.
in other words,
x − Ay ∈ X.
The only assumption made on the matrix A is that aij > 0 for all i, j .
Let us show that, for all x ∈ X, S(x) is a bounded, nonempty set. Denote
380 JOTA: VOL. 131, NO. 3, DECEMBER 2006
p
αi = min aij , 1 ≤ i ≤ n.
j =1
We have αi > 0 for all i. It is obvious that 0 ∈ S(x) and, for all y ∈ S(x), 1 ≤ i ≤
n, we have
p
p
αi
y
= αi yj ≤ aij yj ≤ xi ≤
x
; (11)
j =1 j =1
here,
·
stands for the l 1 norm in Rn or Rp ; i.e.,
x
= ni=1 |xi |. Thus, S(x) ⊆
[0, x1 /α1 ]p ; hence, S(x) is bounded.
Now, we show that S verifies Assumptions (A1) and (A3) of Theorem 3.1.
It is obvious that S(x) is closed for all x ∈ X. Let us show that S is Lipschitz;
i.e., there exists a constant γ > 0 such that, for any x, x ∈ X and for any y ∈ S(x),
there exists y ∈ S(x ) such that
y − y
≤ γ
x − x
. Indeed, consider
γ = max α1−1 , . . . , αn−1
and let x, x ∈ X and y ∈ S(x). From (11), we obtain that
y
≤ γ
x
. Thus, if
x = 0, then y = 0 [i.e., S(0) = {0}] and, for any y ∈ S(x ), we have
y − y
=
y
≤ γ
x
= γ
x − x
.
So, we will assume further that x = 0.
If there exists i0 such that xi0 = 0, take y = 0. We have
y − y
=
y
≤ xi0 /αi0 = (xi0 − xi0 )/αi0
≤ (1/αi0 )
x − x
≤ γ
x − x
.
If, for all i, xi > 0, then take y = ty, where t is defined by t = min{xi /xi |xi >
0} = xi0 /xi0 . We have y ∈ Y, Ay ≤ x , hence y ∈ S(x ) and
y − y
= |1 − t| ·
y
= |xi0 − xi0 |
y
/xi0
≤
x − x
y
/xi0
≤
x − x
αi−1
0
≤ γ
x − x
.
The fact that S is Lipschitz implies the continuity of the set-valued map S :
X→→ Y , hence Assumption (A3) is fulfilled.
To prove (A1’) stated in Remark 3.2 [which in our example is equivalent to
(A1)], note that 0 ≤ x ≤ x ⇒ S(x) ⊆ S(x ). Denote 1n the column vector in Rn
having all entries equal to 1. Given x ∈ X, we can take
Nx ={x ∈ X | x ≤ x + 1n }
JOTA: VOL. 131, NO. 3, DECEMBER 2006 381
={x ∈ X | ∀i : xi ≤ xi + 1}
n
= [0, xi + 1].
i=1
S(x ) = S(x + 1n ),
x ∈Nx
which is bounded.
4. Concluding Remarks
In this paper, we have shown in a very general setting that the semivectorial
bilevel optimization problem (BLw ) can be approached using a penalty method.
Of course, it is interesting and useful to see how this method applies in particular
cases, for example when all the functions are linear and we have polyhedral
constraints. Another important case is when all the functions and constraints are
convex. However, the study of all these cases require special developments. This
will be the object of future works.
References
1. BENSON, H. P., Optimization over the Efficient Set, Journal of Mathematical Analysis
and Applications, Vol. 98, pp. 562–580, 1984.
2. BENSON, H. P., A Finite, Nonadjacent Extreme Point Search Algorithm for Optimization
over the Efficient Set, Journal of Optimization Theory and Applications, Vol. 73, pp. 47–
64, 1992.
3. CRAVEN, B. D., Aspects of Multicriteria Optimization, Recent Prospects in Mathemat-
ical Programming, Gordon and Breach, Philadelphia, Pennsylvania, 1991.
4. FÜLÖP, J., A Cutting Plane Algorithm for Linear Optimization over the Efficient
Set, Generalized Convexity, Lecture Notes in Economics and Mathematical Systems,
Springer Verlag, Berlin, Germany, Vol. 405, 1994.
5. PHILIP, J., Algorithms for the Vector Maximization Problem, Mathematical Program-
ming, Vol. 2, pp. 207–229, 1972.
6. AN, L. T. H., TAO, P. D., and MUU, L. D., Numerical Solution for Optimization over the
Efficient Set by D,C, Optimization Algorithms, Operations Research Letters, Vol. 19,
pp. 117–128, 1996.
7. BOLINTINÉANU, S., Optimality Conditions for Minimization over the (Weakly or Prop-
erly) Efficient Set, Journal of Mathematical Analysis and Applications, Vol. 173,
pp. 523–541, 1993.
8. BOLINTINÉANU, S., Minimization of a Quasiconcave Function over an Efficient Set,
Mathematical Programming, Vol. 61, pp. 89–110, 1993.
382 JOTA: VOL. 131, NO. 3, DECEMBER 2006