Abstract

We study a class of multiobjective bilevel programs with the weights of objectives being uncertain and assumed to belong to convex and compact set. To the best of our knowledge, there is no study about this class of problems. We use a worst-case weighted approach to solve this class of problems. Our “worst-case weighted multiobjective bilevel programs” model supposes that each player (leader or follower) has a set of weights to their objectives and wishes to minimize their maximum weighted sum objective where the maximization is with respect to the set of weights. This new model gives rise to a new Pareto optimum concept, which we call “robust-weighted Pareto optimum”; for the worst-case weighted multiobjective optimization with the weight set of each player given as a polytope, we show that a robust-weighted Pareto optimum can be obtained by solving mathematical programing with equilibrium constraints (MPEC). For an application, we illustrate the usefulness of the worst-case weighted multiobjective optimization to a supply chain risk management under demand uncertainty. By the comparison with the existing weighted approach, we show that our method is more robust and can be more efficiently applied to real-world problems.

1. Introduction

Multiobjective bilevel programs are a generalization of the scalar criterion bilevel programs [1] and used to model situations where decision-makers at both upper and lower levels, called players, take actions by considering their individual multiple objectives [24]. Many papers have been published in the last two decades about (multiobjective) bilevel programs because of their possible applications [1, 4, 5].

In this paper, we consider a special class of multiobjective bilevel programs with the weights of objectives being uncertain and assumed to belong to convex and compact set. To the best our knowledge, there is no study about this class of problems and the existing approaches in multiobjective bilevel programs can not solve this class of problems. Therefore, new methods are necessary to be proposed. In this paper, we present a robust-weighted approach to solve this class of multiobjective bilevel programs. In our model, we suppose that players at both upper and lower levels are risk-averse and a player uses the robust optimization approach to manage the weight uncertainty, assuming that the other player is a robust optimizer as well. Note that the robust optimization approach is not concerning problem data or the other players’ strategies but concerning the weights to the objectives. We note that this model is closely related to robust optimization [6, 7]. It can be seen as a natural extension of the robust modeling technique by replacing a decision variable with a general deterministic function [6, 7]. However, our model is also different from the studies in [6, 7] that we study the bilevel problem, while the model in [6, 7] is single level.

Multiobjective optimization has been extensively studied and a number of different approaches have also been proposed [8]. In multiobjective problems, since there are several competing objectives to be considered and it is not possible to simultaneously optimize all objectives, a commonly accepted approach for coping with this setting is the weighted approach by assigning a nonnegative weight by considering the importance of the corresponding objective function. And then the decision-maker can make a decision by optimizing a weighted sum objective [912]. A weighted Pareto optimal point can be obtained. In the setting of multiobjective models including input from multiple experts, a decision can often be obtained only after considering different opinions from different experts [11].

However current approaches to model-based multiexpert multiobjective decision-making have several shortcomings. How to choose a solution among (possibly infinite) generated Pareto optimal solutions? There is no guidance. In multiexpert optimization, the relative weights given by experts can differ significantly as experts with differing opinions often assign different relative weights to objectives. As shown in applications, the weights are not known in advance and the player has to choose them. Ambiguity often exists in the choice of the weights to objectives, as it is not easy to decide relative weights for each objective. In addition, as shown in multiobjective optimization in the literature, relative weights given by the same decision-maker may rely on the elicitation methods [12, 13]. In the bilevel setting considered in this paper, it is more difficult for the decision-maker at the upper level to exactly assume the weights of the decision-maker’s objectives at the lower level. Therefore, in this paper, we focus on a special class of multiobjective bilevel programs with the weights of objectives being uncertain and assumed to belong to convex and compact set. There is no study about this class of problems and there no methods in multiobjective bilevel programs can solve this class of problems. It is necessary to provide a new approach to cope with these issues.

Hence, our motivation to utilize a worst-case weighted approach is that it provides an alternative way to deal with the weights ambiguity. Furthermore, if each player in multiobjective bilevel problems chooses the worst-case weighted approach, then we show that the computation for robust-weighted Pareto optimum, with the choice of polytope weight set for every player, is reformulated as a solution to mathematical programing with equilibrium constraints (MPEC) which can be solved by the existing methods (e.g., the sequential quadratic programing (SQP) methods).

Though multiobjective bilevel programs have not attracted much attention in the literature, there are some interesting potential applications. One example is the multicriteria Stackelberg competition of a supply chain containing a manufacturer who supplies a set of products to a risk-averse retailer satisfying uncertain consumer demand. In our model, the manufacturer decides on the quantity for each product so as to maximize the profit and minimize the cost simultaneously by forecasting the order quantity from the retailer and the wholesale prices resulting from market clearing conditions. The retailer also decides their wholesale market order quantity for each product in order to simultaneously maximize the mean profit and minimize the standard deviation of the profit.

Since the possible applications, multiobjective bilevel programs have attracted some attention. For example, Yin [14] considers a multiobjective bilevel model for transportation planning and management problems, where genetic algorithms are proposed to solve the resulting model. Deb and Sinha [15] present evolutionary algorithms for solving multiobjective optimization problems. Pieume et al. [16] develop two methods for solving bilevel linear multiobjective optimization problems. Eichfelder [17] gives a solution method for solving nonlinear multiobjective bilevel problems based on a scalarization approach and the sensitivity analysis of adaptive parameters. To reduce traffic congestion as well as to improve workforce productivity, a bilevel multiobjective model is proposed for an urban logistics metropolis [18]. A discrete approach is proposed to solve the resulting model. The convergence result and numerical tests are also provided. Some special multiobjective bilevel programs with linear objectives are also studied extensively. For example, a linear bilevel optimization problem with multiple objectives at the upper level is studied in [19], where the original problem reduces to solving a series of linear bilevel problems with a single objective function at each level; a linear bilevel optimization problem with multiple objectives at the lower level is studied in [20], where the original problem is reformulated as an optimization problem over a nonconvex region given by a union of faces of the polyhedron defined by all constraints.

Our study is different from the above papers in that we focus on a special class of multiobjective bilevel nonlinear programs—the weights of objectives are uncertain and assumed to belong to convex and compact set. To deal with this uncertainty, a robust optimization approach is used. Then we reformulated the original problem as a robust bilevel nonlinear optimization problem. To the best of our knowledge, this is the first paper to consider these special multiobjective bilevel problems.

We feel that the primary contributions of this paper are as follows. We deal with a special class of multiobjective bilevel programs, where the weights of objectives are uncertain and assumed to belong to convex and compact set. To the best of our knowledge, this is the first paper to consider such multiobjective bilevel programs. In our model, we suppose that players at both upper and lower levels are risk-averse and a player uses the robust optimization approach to manage the weight uncertainty, assuming that the other player is a robust optimizer as well. As we know that, there is no method to solve such problem. We propose a worst-case weighted approach for solving such multiobjective bilevel programs, extending the notion of robust-weighted multiobjective optimization models [21] to multiobjective bilevel problems. We show that a robust-weighted Pareto optimal point can be calculated by solving MPEC when the weight sets chosen by the players are polytopes. We note that the MPEC has been extensively studied and can be solved by the sequential quadratic programing (SQP) method [22]. We demonstrate the usefulness of the worst-case weighted method in a bilevel multiobjective competition problem within supply chain. We note that, compared with the existing weighted approach [9], our method is more “robust" in that there are different solutions for choosing different weights by using weighted approach [9], but there is a unique solution by using our approach.

Throughout this paper we use the following notations. is the set of real numbers; denotes the set of nonnegative real numbers; and is the set of strictly positive real numbers. Let For any , denotewhere is the index set. Given any vector function , by , we indicate that is a continuously differentiable function from to and we use to denote the gradient of the function at ; for simplicity, in this paper we use superscript (resp., ) to denote decision-makers in the upper level (resp., decision-makers in the lower level) parameters, decision variables, cost functions, and weights of payoff functions. We note that the different mechanisms (e.g., bilevel, equilibrium, multiobjective optimization, weighting, and Pareto) are used in this paper. To specially understand the links between all these components, we provide a figure (see Figure 1) to illustrate this.

In Section 2, we review some basic concepts in multiobjective bilevel programs and then introduce a new concept of the robust-weighted Pareto optimum. We show that a robust-weighted Pareto optimal solution is also a weighted Pareto optimum. In Section 3, we show that when the weight sets are given as polytopes, the robust-weighted Pareto optimum can be cast as MPEC which can be efficiently solved by SQP methods. In Section 4, we show the efficiency of our approach by an example for analyzing a multiobjective bilevel competition problem within supply chain. Finally, we conclude this paper in Section 5.

2. Preliminaries

In this paper, we consider the following multiobjective bilevel optimization:where signifies the set of solutions of the lower-level problemand , , and , , . In what follows, we assume that all functions considered in the bilevel problem (3) are continuously differentiable and ,  . Since we consider the case of a multiobjective optimization () at the lower level, given , the solution set denotes the Pareto optimal set or weak Pareto optimal set; that is, given any and , there is no feasible point to (4) such that or .

Definition 1. is said to be feasible for problem (3) if , , and . is said to be locally Pareto optimal (resp., weak Pareto) for problem (3), if it is feasible for (3) and there exists no other feasible point in the neighborhood of such that (resp., ).

In the multiobjective optimization literature, the weighted approach is one of the most often used methods. However the weighted approach is also criticized for several shortcomings. For example, as shown in applications, weights are not known in advance and the modeler or decision-maker has to choose them while it is often difficult to make a reasonable decision for choosing appropriate weights for objectives. In many situations there may be uncertainty in the weights provided by the decision-maker. Therefore there is a need for providing a novel methodology to cope with these issues. The robust-weighted approach is one possible choice as it provides an alternative to solve the weight uncertainty in the multiobjective bilevel model. Specifically, we assume that players in both levels use the robust-weighted approach to address the weight uncertainty; that is, both the leader and the follower use a robust-weighted approach to minimize their weighted payoff functions.

Definition 2. Let be a nonempty closed convex set of weights, . is called a robust-weighted optimal solution of (3), if it is an optimum to the following robust optimization problem:where signifies the set of solutions of the lower-level problem:

Different from weighted Pareto concepts, the robust-weighted optimum proposed in this paper seeks to find an optimal solution that is robust and feasible for the worst-case weight within the family of weights. It directly follows from Theorem of Hu and Mehrotra [21] and Definitions 1 and 2 that the following theorem holds.

Theorem 3. Let be defined as in Definition 2, . if is a robust-weighted optimum of (5) with the combination of weight sets, then we have the following results: (1) is a weak Pareto optimum of (3).(2)if all of the weights in are positive, then is a Pareto optimum of (3).(3)if is the unique optimal solution of (5), then is a Pareto optimum of (3).

3. MPEC Reformulations for (5)

Now that a robust-weighted optimum of (5) is a (weak) Pareto optimum of (3), our next step is to discuss how to realize such a solution. To do this, we assume that the weight sets , , are given as polyhedral regions. Then the problem of computing a robust-weighted optimal point could be cast as MPEC which has been extensively studied and can be solved by the solver PATH [23].

We suppose that the uncertain weights can be described by a fixed reference point and a perturbation region around the point; that is, the weight sets and can be given as follows: let be the reference point of and be a coefficient matrix used to construct a perturbation region around , , and then define the perturbation region around as follows:where , are two uncertain sets which belong toand here , . Note that the above method for defining uncertain parameters has been widely used in the robust optimization literature [12, 13].

Definewhere , , , and , . We now present a reformulation for problem (5) with given as (7)–(9). Next, we show that the computation for a robust-weighted Nash equilibrium point of MBG can be casted as solving MPEC.

Lemma 4. Suppose the weight sets , , are given as (7)–(9), and then the robust-weighted bilevel optimization problem (5) is equivalent to the following bilevel optimization problem:where signifies the set of solutions of the lower-level problem:

Proof. This lemma can be proven as follows. Given the leader’s strategy , it is easy to show that the follower’s optimization problem (6) can be rewritten asFor any given , the left-hand side in the first inequality constraint in the above problem is equivalent toThe corresponding dual problem iswhere “0” and “” are the zero and unit vectors with appropriate dimensions, respectively.
Therefore it follows from (12)–(14) and the strong duality theorem that the assertion of this theorem is true.

The above lemma shows the equivalence between the bilevel problems (5) and (10); that is, if is a robust-weighted optimal solution to (5), then there is such that is an optimal solution to (10); otherwise if is an optimal solution to (10), then is a robust-weighted optimal solution to (5). So, in what follows, we focus on the bilevel problem (10).

Given , it is expected that a solution of the lower-level problem (11) is a stationary point of (11); that is, it solves the following Karush-Kuhn-Tucker (KKT) system:where is the Lagrangian of (11) defined as

Therefore replacing the lower level optimization problem (11) with (15) derives MPEC:

This problem can be efficiently solved by the solver PATH [23]. We note that the PATH algorithm is originally described in [24]. This algorithm computes a stationary solution of MPEC (since it solves an equation related to the necessary optimality conditions). If for any given functions and are convex in and is concave in , then under some constraint qualification an optimal solution of the bilevel optimization problem (10) is also a stationary solution of the MPEC (17) [25].

Theorem 5. Suppose the weight sets , , are given as (7)–(9). If for any given functions and are convex in and is concave in , then under some constraint qualification an optimal solution of the bilevel optimization problem (5) is also a stationary solution of the MPEC (17).

Proof. The proof for this theorem follows directly from Lemma 4, the above discussions, and Proposition [25].

4. An Application

4.1. Model Description

In this section, we model the multicriteria competition in a supply chain of manufacturers who cooperatively decide on volume to supply products to risk-averse retailers who also cooperatively decide on volume to satisfy uncertain consumer demand. Manufacturers have two objectives (i.e., to maximize their profit and to minimize their cost of products supplied simultaneously) by choosing the supply quantity for each of products in the anticipation of the order quantities of the retailers and also anticipating of the wholesale prices resulting from the market clearing conditions. Retailers have also two objectives (i.e., to maximize the expected utility from retail sales and to minimize the risk expressed as the standard deviation of the profit) by deciding the wholesale market order quantity for each of products.

The retail demand satisfies the following stochastic linear inverse demand function: where is the aggregate price vector, is the vector of the th retailer prices with , is the price of the th retailer for the th product, and with as a positive scalar random variable with mean and finite standard deviation , , and . The vector with gives the prices that the consumers would be willing to pay if the retail market supply was 0 and is a finite number. The matrix is the positive definite matrix of inverse demand sensitivities which can be expressed as with ,  , , and  , . The vector is the aggregate retailer order vector, where is the vector of the th retailer order vector .

Retailers make their order quantity decision as follows:where is the vector of wholesale prices and , for th retailer whose first objective is to maximize their mean profit and second objective is to minimize their standard deviation for the stochastic demand. We assume that the weight vector to the two objectives and in (20) is .

Manufacturers choose their supply quantity to maximize their profit in the wholesale market, anticipating the market clearing wholesale market price as well as the retailer order quantities. Mathematically, manufacturers face the following decision problem:where is the th manufacturer unit cost, is the supply quantity of the th manufacturer and is the th retailer order quantity for a given wholesale price , is the weight vector to manufacturer i’s profit and cost, and the first constraint is the wholesale market clearing conditions. We assume that all the uncertainty exists in the retail market and that the wholesale market is deterministic except for the weight choice. The reason for this is that the demand uncertainty in the retail demand function is realized when all wholesale market related decisions are made in our model. We also assume that the weight vector is uncertain and belongs to a known convex and compact set. By this assumption, we can utilize the proposed robust-weighted approach.

4.2. Numerical Results

We now detail the numerical tests for finding a robust-weighted manufacturer-retailer equilibrium by using the numerical solver PATH [23]. We present numerical results for the robust-weighted manufacturer-retailer equilibrium of models in the decentralized supply chain with finite weights and general polyhedral weight regions, respectively. Suppose , , and , that is, two manufacturers, two retailers, and two products. We assume that and ; that is, the mean and the standard deviation of the stochastic variable are 1 and 0.5, respectively. We give the following assumption for the parameters in supply chain:

Analytical Hierarchy Process (for brief AHP) [26] is used to assign weights for the two objectives in models (20) and (21). We assume that retailer is unsure about the relative importance of and for any . Consequently, we first consider that there are finite weights to quantify the relative importance for the two objective functions. Then we generalize the finite weights to the weight set given by a polytope. We first suppose that retailers allow three possible scale indices for the two objective functions as , 1, or 1.5 and , 1.28, or 1.69. Similar to the assumption about the retailers’ weights, we also assume that manufacturer is also unsure about the relative importance of and , . As such, manufacturers allow for two possible scale indices to the two objective functions as or 1 and or 1.2.

The numerical results are reported in Tables 14, where (manufacturer-retailer equilibriums) are the weighted Nash equilibrium by using the weighted approach [9], and here the vector is the aggregate manufacturer supply vector. The robust-weighted Nash equilibrium (RWNE) with finite weights is reported in the last row of Table 4. There are different weighted manufacturer-retailer equilibriums for a different choice of the weights. The robust-weighted retailer equilibrium is computed by using the weight region generated as the convex hull of the nine weight vectors. It is easy to see that the robust-weighted approach presents a unique equilibrium which mitigates the conflict with different weighted manufacturer-retailer equilibria for a different choice of the weights in Tables 14. We note that, compared with the existing weighted approach [9], our method is more “robust" in that there are different solutions for choosing different weights by using weighted approach [9], but there is a unique solution by using our approach.

The managerial impact of the proposed method is as follows. The proposed approach overcomes several shortcomings for the current approaches to the multicriteria competition in a supply chain. There is no clear guidance to a decision-maker (manufacturers or retailers) to choose a solution among (possibly a very large number of) generated Pareto optimal solutions. However, there is a unique solution by using our approach. In the multicriteria competition in a supply chain, since different decision-maker (manufacturers or retailers) may have different opinions, the relative weights given by different decision-maker can differ significantly. Furthermore, there may be ambiguity in the weights given by a decision-maker (manufacturer or retailer), because it is often difficult to give relative weights for each objective.

For the general polytope weight region, we suppose the reference point and matrices and are defined as follows: ,  , , , and

The perturbation subsets are defined aswhere is a parameter defined in . With the above data the robust-weighted manufacturer-retailer equilibrium with different can be calculated by utilizing the method developed in Section 5. The stability results are shown in Table 5 which contains small data errors due to the small changes in .

5. Conclusion

This paper introduces a robust-weighted approach to multicriteria bilevel problems. The computation of the robust-weighted Pareto optimal point, with the polyhedral weight sets, can be equivalently transformed into solving MPEC. As an application, the new approach is used to analyze supply chain multicriteria bilevel competition. The numerical tests show that our method is more robust than existing approaches and a unique robust-weighted Pareto optimum can always be obtained for general polyhedral weight sets.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work is supported by The Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning (no. TP2014043). The work is also supported by Natural Scientific Foundation of China (no. 71401107; no. 71571055).