Techniques for Multiobjective Optimization with Discrete Variables: Boxed Line Method and Tchebychev Weight Set Decomposition

Many real-world applications involve multiple competing objectives, but due to conflict between the objectives, it is generally impossible to find a feasible solution that optimizes all, simultaneously. In contrast to single objective optimization, the goal in multiobjective optimization is to generate a set of solutions that induces the nondominated (ND) frontier. This thesis presents two techniques for multiobjective optimization problems with discrete decision variables. First, the Boxed Line Method is an exact, criterion space search algorithm for biobjective mixed integer programs (Chapter 2). A basic version of the algorithm is presented with a recursive variant and other enhancements. The basic and recursive variants permit complexity analysis, which yields the first complexity results for this class of algorithms. Additionally, a new instance generation method is presented, and a rigorous computational study is conducted. Second, a novel weight space decomposition method for integer programs with three (or more) objectives is presented with unique geometric properties (Chapter 3). The weighted Tchebychev scalarization used for this weight space decomposition provides the benefit of including unsupported ND images but at the cost of convexity of weight set components. This work proves convexity-related properties of the weight space components, including star-shapedness. Further, a polytopal decomposition is used to properly define dimension for these nonconvex components. The weighted Tchebychev weight set decomposition is then applied as a “dual” perspective on the class of multiobjective “primal” algorithms (Chapter 4). It is shown that existing algorithms do not yield enough information for a complete decomposition, and the necessary modifications required to yield the missing information is proven. Modifications for primal algorithms to compute inner and outer approximations of the weight space components are presented. Lastly, a primal algorithm is restricted to solving for a subset of the ND frontier, where this subset represents the compromise between multiple decision makers’ weight vectors.Ph.D

Optimising a nonlinear utility function in multi-objective integer programming

In this paper we develop an algorithm to optimise a nonlinear utility function of multiple objectives over the integer efficient set. Our approach is based on identifying and updating bounds on the individual objectives as well as the optimal utility value. This is done using already known solutions, linear programming relaxations, utility function inversion, and integer programming. We develop a general optimisation algorithm for use with k objectives, and we illustrate our approach using a tri-objective integer programming problem.Comment: 11 pages, 2 tables; v3: minor revisions, to appear in Journal of Global Optimizatio

Multi-objective integer programming: An improved recursive algorithm

This paper introduces an improved recursive algorithm to generate the set of all nondominated objective vectors for the Multi-Objective Integer Programming (MOIP) problem. We significantly improve the earlier recursive algorithm of \"Ozlen and Azizo\u{g}lu by using the set of already solved subproblems and their solutions to avoid solving a large number of IPs. A numerical example is presented to explain the workings of the algorithm, and we conduct a series of computational experiments to show the savings that can be obtained. As our experiments show, the improvement becomes more significant as the problems grow larger in terms of the number of objectives.Comment: 11 pages, 6 tables; v2: added more details and a computational stud

Efficiently Constructing Convex Approximation Sets in Multiobjective Optimization Problems

Convex approximation sets for multiobjective optimization problems are a well-studied relaxation of the common notion of approximation sets. Instead of approximating each image of a feasible solution by the image of some solution in the approximation set up to a multiplicative factor in each component, a convex approximation set only requires this multiplicative approximation to be achieved by some convex combination of finitely many images of solutions in the set. This makes convex approximation sets efficiently computable for a wide range of multiobjective problems - even for many problems for which (classic) approximations sets are hard to compute. In this article, we propose a polynomial-time algorithm to compute convex approximation sets that builds upon an exact or approximate algorithm for the weighted sum scalarization and is, therefore, applicable to a large variety of multiobjective optimization problems. The provided convex approximation quality is arbitrarily close to the approximation quality of the underlying algorithm for the weighted sum scalarization. In essence, our algorithm can be interpreted as an approximate variant of the dual variant of Benson's Outer Approximation Algorithm. Thus, in contrast to existing convex approximation algorithms from the literature, information on solutions obtained during the approximation process is utilized to significantly reduce both the practical running time and the cardinality of the returned solution sets while still guaranteeing the same worst-case approximation quality. We underpin these advantages by the first comparison of all existing convex approximation algorithms on several instances of the triobjective knapsack problem and the triobjective symmetric metric traveling salesman problem

An outer approximation algorithm for multi-objective mixed-integer linear and non-linear programming

In this paper, we present the first outer approximation algorithm for multi-objective mixed-integer linear programming problems with any number of objectives. The algorithm also works for certain classes of non-linear programming problems. It produces the non-dominated extreme points as well as the facets of the convex hull of these points. The algorithm relies on an oracle which solves single-objective weighted-sum problems and we show that the required number of oracle calls is polynomial in the number of facets of the convex hull of the non-dominated extreme points in the case of multiobjective mixed-integer programming (MOMILP). Thus, for MOMILP problems for which the weighted-sum problem is solvable in polynomial time, the facets can be computed with incremental-polynomial delay. From a practical perspective, the algorithm starts from a valid lower bound set for the non-dominated extreme points and iteratively improves it. Therefore it can be used in multi-objective branch-and-bound algorithms and still provide a valid bound set at any stage, even if interrupted before converging. Moreover, the oracle produces Pareto optimal solutions, which makes the algorithm also attractive from the primal side in a multi-objective branch-and-bound context. Finally, the oracle can also be called with any relaxation of the primal problem, and the obtained points and facets still provide a valid lower bound set. A computational study on a set of benchmark instances from the literature and new non-linear multi-objective instances is provided.Comment: 21 page

Multi-objective combinatorial optimization problems in transportation and defense systems

Multi-objective Optimization problems arise in many applications; hence, solving them efficiently is important for decision makers. A common procedure to solve such problems is to generate the exact set of Pareto efficient solutions. However, if the problem is combinatorial, generating the exact set of Pareto efficient solutions can be challenging. This dissertation is dedicated to Multi-objective Combinatorial Optimization problems and their applications in system of systems architecting and railroad track inspection scheduling. In particular, multi-objective system of systems architecting problems with system flexibility and performance improvement funds have been investigated. Efficient solution methods are proposed and evaluated for not only the system of systems architecting problems, but also a generic multi-objective set covering problem. Additionally, a bi-objective track inspection scheduling problem is introduced for an automated ultrasonic inspection vehicle. Exact and approximation methods are discussed for this bi-objective track inspection scheduling problem --Abstract, page iii

Domination and Decomposition in Multiobjective Programming

During the last few decades, multiobjective programming has received much attention for both its numerous theoretical advances as well as its continued success in modeling and solving real-life decision problems in business and engineering. In extension of the traditionally adopted concept of Pareto optimality, this research investigates the more general notion of domination and establishes various theoretical results that lead to new optimization methods and support decision making. After a preparatory discussion of some preliminaries and a review of the relevant literature, several new findings are presented that characterize the nondominated set of a general vector optimization problem for which the underlying domination structure is defined in terms of different cones. Using concepts from linear algebra and convex analysis, a well known result relating nondominated points for polyhedral cones with Pareto solutions is generalized to nonpolyhedral cones that are induced by positively homogeneous functions, and to translated polyhedral cones that are used to describe a notion of approximate nondominance. Pareto-oriented scalarization methods are modified and several new solution approaches are proposed for these two classes of cones. In addition, necessary and sufficient conditions for nondominance with respect to a variable domination cone are developed, and some more specific results for the case of Bishop-Phelps cones are derived. Based on the above findings, a decomposition framework is proposed for the solution of multi-scenario and large-scale multiobjective programs and analyzed in terms of the efficiency relationships between the original and the decomposed subproblems. Using the concept of approximate nondominance, an interactive decision making procedure is formulated to coordinate tradeoffs between these subproblems and applied to selected problems from portfolio optimization and engineering design. Some introductory remarks and concluding comments together with ideas and research directions for possible future work complete this dissertation

Biobjective Optimization over the Efficient Set Methodology for Pareto Set Reduction in Multiobjective Decision Making: Theory and Application

A large number of available solutions to choose from poses a significant challenge for multiple criteria decision making. This research develops a methodology that reduces the set of efficient solutions under consideration. This dissertation is composed of three major parts: (i) the formalization of a theoretical framework; (ii) the development of a solution approach; and (iii) a case study application of the methodology. In the first part, the problem is posed as a multiobjective optimization over the efficient set and considers secondary robustness criteria when the exact values of decision variables are subjected to uncertainties during implementation. The contributions are centered at the modeling of uncertainty directly affecting decision variables, the use of robustness to provide additional trade-off analysis, the study of theoretical bounds on the measures of robustness, and properties to ensure that fewer solutions are identified. In the second part, the problem is reformulated as a biobjective mixed binary program and the secondary criteria are generalized to any convenient linear functions. A solution approach is devised in which an auxiliary mixed binary program searches for unsupported Pareto outcomes and a novel linear programming filtering excludes any dominated solutions in the space of the secondary criteria. Experiments show that the algorithm tends to run faster than existing approaches for mixed binary programs. The algorithm enables dealing with continuous Pareto sets, avoiding discretization procedures common to the related literature. In the last part, the methodology is applied in a case study regarding the electricity generation capacity expansion problem in Texas. While water and energy are interconnected issues, to the best of our knowledge, this is the first study to consider both water and cost objectives. Experiments illustrate how the methodology can facilitate decision making and be used to answer strategic questions pertaining to the trade-off among different generation technologies, power plant locations, and the effect of uncertainty. A simulation shows that robust solutions tend to maintain feasibility and stability of objective values when power plant design capacity values are perturbed