303 research outputs found

    Optimality-based bound contraction with multiparametric disaggregation for the global optimization of mixed-integer bilinear problems

    Get PDF
    We address nonconvex mixed-integer bilinear problems where the main challenge is the computation of a tight upper bound for the objective function to be maximized. This can be obtained by using the recently developed concept of multiparametric disaggregation following the solution of a mixed-integer linear relaxation of the bilinear problem. Besides showing that it can provide tighter bounds than a commercial global optimization solver within a given computational time, we propose to also take advantage of the relaxed formulation for contracting the variables domain and further reduce the optimality gap. Through the solution of a real-life case study from a hydroelectric power system, we show that this can be an efficient approach depending on the problem size. The relaxed formulation from multiparametric formulation is provided for a generic numeric representation system featuring a base between 2 (binary) and 10 (decimal)

    An improved multi-parametric programming algorithm for flux balance analysis of metabolic networks

    Full text link
    Flux balance analysis has proven an effective tool for analyzing metabolic networks. In flux balance analysis, reaction rates and optimal pathways are ascertained by solving a linear program, in which the growth rate is maximized subject to mass-balance constraints. A variety of cell functions in response to environmental stimuli can be quantified using flux balance analysis by parameterizing the linear program with respect to extracellular conditions. However, for most large, genome-scale metabolic networks of practical interest, the resulting parametric problem has multiple and highly degenerate optimal solutions, which are computationally challenging to handle. An improved multi-parametric programming algorithm based on active-set methods is introduced in this paper to overcome these computational difficulties. Degeneracy and multiplicity are handled, respectively, by introducing generalized inverses and auxiliary objective functions into the formulation of the optimality conditions. These improvements are especially effective for metabolic networks because their stoichiometry matrices are generally sparse; thus, fast and efficient algorithms from sparse linear algebra can be leveraged to compute generalized inverses and null-space bases. We illustrate the application of our algorithm to flux balance analysis of metabolic networks by studying a reduced metabolic model of Corynebacterium glutamicum and a genome-scale model of Escherichia coli. We then demonstrate how the critical regions resulting from these studies can be associated with optimal metabolic modes and discuss the physical relevance of optimal pathways arising from various auxiliary objective functions. Achieving more than five-fold improvement in computational speed over existing multi-parametric programming tools, the proposed algorithm proves promising in handling genome-scale metabolic models.Comment: Accepted in J. Optim. Theory Appl. First draft was submitted on August 4th, 201

    Efficiently Solving Repeated Integer Linear Programming Problems by Learning Solutions of Similar Linear Programming Problems using Boosting Trees

    Get PDF
    It is challenging to obtain online solutions of large-scale integer linear programming (ILP) problems that occur frequently in slightly different forms during planning for autonomous systems. We refer to such ILP problems as repeated ILP problems. The branch-and-bound (BAB) algorithm is commonly used to solve ILP problems, and a significant amount of computation time is expended in solving numerous relaxed linear programming (LP) problems at the nodes of the BAB trees. We observe that the relaxed LP problems, both within a particular BAB tree and across multiple trees for repeated ILP problems, are similar to each other in the sense that they contain almost the same number of constraints, similar objective function and constraint coefficients, and an identical number of decision variables. We present a boosting tree-based regression technique for learning a set of functions that map the objective function and the constraints to the decision variables of such a system of similar LP problems; this enables us to efficiently infer approximately optimal solutions of the repeated ILP problems. We provide theoretical performance guarantees on the predicted values and demonstrate the effectiveness of the algorithm in four representative domains involving a library of benchmark ILP problems, aircraft carrier deck scheduling, vehicle routing, and vehicle control

    Advanced multiparametric optimization and control studies for anaesthesia

    Get PDF
    Anaesthesia is a reversible pharmacological state of the patient where hypnosis, analgesia and muscle relaxation are guaranteed and maintained throughout the surgery. Analgesics block the sensation of pain; hypnotics produce unconsciousness, while muscle relaxants prevent unwanted movement of muscle tone. Controlling the depth of anaesthesia is a very challenging task, as one has to deal with nonlinearity, inter- and intra-patient variability, multivariable characteristics, variable time delays, dynamics dependent on the hypnotic agent, model analysis variability, agent and stability issues. The modelling and automatic control of anaesthesia is believed to (i) benefit the safety of the patient undergoing surgery as side-effects may be reduced by optimizing the drug infusion rates, and (ii) support anaesthetists during critical situations by automating the drug delivery systems. In this work we have developed several advanced explicit/multi-parametric model predictive (mp-MPC) control strategies for the control of depth of anaesthesia. State estimation techniques are developed and used simultaneously with mp-MPC strategies to estimate the state of each individual patient, in an attempt to overcome the challenges of inter- and intra- patient variability, and deal with possible unmeasurable noisy outputs. Strategies to deal with the nonlinearity have been also developed including local linearization, exact linearization as well as a piece-wise linearization of the Hill curve leading to a hybrid formulation of the patient model and thereby the development of multiparametric hybrid model predictive control methodology. To deal with the inter- and intra- patient variability, as well as the noise on the process output, several robust techniques and a multiparametric moving horizon estimation technique have been design and implemented. All the studies described in the thesis are performed on clinical data for a set of 12 patients who underwent general anaesthesia.Open Acces

    An Algorithm for Biobjective Mixed Integer Quadratic Programs

    Get PDF
    Multiobjective quadratic programs (MOQPs) are appealing since convex quadratic programs have elegant mathematical properties and model important applications. Adding mixed-integer variables extends their applicability while the resulting programs become global optimization problems. Thus, in this work, we develop a branch and bound (BB) algorithm for solving biobjective mixed-integer quadratic programs (BOMIQPs). An algorithm of this type does not exist in the literature. The algorithm relies on five fundamental components of the BB scheme: calculating an initial set of efficient solutions with associated Pareto points, solving node problems, fathoming, branching, and set dominance. Considering the properties of the Pareto set of BOMIQPs, two new fathoming rules are proposed. An extended branching module is suggested to cooperate with the node problem solver. A procedure to make the dominance decision between two Pareto sets with limited information is proposed. This set dominance procedure can eliminate the dominated points and eventually produce the Pareto set of the BOMIQP. Numerical examples are provided. Solving multiobjective quadratic programs (MOQPs) is fundamental to our research. Therefore, we examine the algorithms for this class of problems with different perspectives. The scalarization techniques for (strictly) convex MOPs are reviewed and the available algorithms for computing efficient solutions for MOQPs are discussed. These algorithms are compared with respect to four properties of MOQPs. In addition, methods for solving parametric multiobjective quadratic programs are studied. Computational studies are provided with synthetic instances, and examples in statistics and portfolio optimization. The real-life context reveals the interplay between the scalarizations and provides an additional insight into the obtained parametric solution sets

    Optimal Scheduling of Combined Heat and Power Generation Considering Heating Grid Dynamics

    Get PDF
    As the share of renewable generation increases in electric grids, the traditionally heat driven operation of combined heat and power plants (CHPs) reaches its limits. Thermal storage is required for a flexible operation of CHPs. This work proposes three novel methods to use a heating grid as thermal storage by exploiting its thermal dynamics. These include the first approach proving global optimality, a novel linear formulation of grid dynamics and an easily real world applicable approach
    corecore