A New Predictor-Corrector Method for Optimal Power Flow

Size: px
Start display at page:

Download "A New Predictor-Corrector Method for Optimal Power Flow"

Transcription

1 1 st Brazilian Workshop on Interior Point Methods Federal Univesity of Technology Paraná (UTFPR) A New Predictor-Corrector Method for Optimal Power Flow Roy Wilhelm Probst April, Campinas, Brazil Roy Wilhelm Probst A New PC Method for Optimal Power Flow 1 / 20

2 Linear Programming Nonlinear Programming Complete PC (CPC) Affine-Scaling and Central Path Methods min s.t. c t x Ax = b x 0 A A t I A A t I x y z x y z or = = b Ax c A t z XZu b Ax c A t z µu XZu Roy Wilhelm Probst A New PC Method for Optimal Power Flow 2 / 20

3 Linear Programming Nonlinear Programming Complete PC (CPC) Affine-Scaling and Central Path Methods min s.t. c t x Ax = b x 0 A A t I A A t I x y z x y z or = = b Ax c A t z XZu b Ax c A t z µu XZu Roy Wilhelm Probst A New PC Method for Optimal Power Flow 2 / 20

4 Predictor-Corrector Method Linear Programming Nonlinear Programming Complete PC (CPC) min s.t. c t x Ax = b x 0 A A t I x ỹ z = r 1 r 2 r 3 A A t I ˆx ŷ ẑ = 0 0 µu X Zu Roy Wilhelm Probst A New PC Method for Optimal Power Flow 3 / 20

5 Linear Programming Nonlinear Programming Complete PC (CPC) PC Method for Nonlinear Programming min f (x) s.t. g(x) = 0 x 0 x g(x) xxl(x, y, z) x g(x) t I x ỹ z = r 1 r 2 r 3 x g(x) xxl(x, y, z) x g(x) t I x ỹ z = 0 0 µu X Zu Roy Wilhelm Probst A New PC Method for Optimal Power Flow 4 / 20

6 Linear Programming Nonlinear Programming Complete PC (CPC) PC Method for Nonlinear Programming min f (x) s.t. g(x) = 0 x 0 x g(x) xxl(x, y, z) x g(x) t I x ỹ z = r 1 r 2 r 3 x g(x) xxl(x, y, z) x g(x) t I x ỹ z = 0 0 µu X Zu Roy Wilhelm Probst A New PC Method for Optimal Power Flow 4 / 20

7 Complete PC Corrections Linear Programming Nonlinear Programming Complete PC (CPC) min s.t. x t Hx XAx = b x 0 XZu = 0 Z x + X z = XZu (X + X )(Z + Z)u = X Zu Roy Wilhelm Probst A New PC Method for Optimal Power Flow 5 / 20

8 Complete PC Corrections Linear Programming Nonlinear Programming Complete PC (CPC) min s.t. x t Hx XAx = b x 0 XZu = 0 Z x + X z = XZu (X + X )(Z + Z)u = X Zu Roy Wilhelm Probst A New PC Method for Optimal Power Flow 5 / 20

9 Primal Corrections Interior-Point Methods Linear Programming Nonlinear Programming Complete PC (CPC) min s.t. x t Hx XAx = b x 0 XAx b = 0 (XA + diag(ax)) x = b XAx (X + X )A(x + x) b = XA x Roy Wilhelm Probst A New PC Method for Optimal Power Flow 6 / 20

10 Dual Corrections Interior-Point Methods Linear Programming Nonlinear Programming Complete PC (CPC) min s.t. x t Hx XAx = b x 0 2Hx + A t Yx + YAx = 0 (2H + A t Y + YA) x + (A t X + diag(ax)) y = 2Hx A t Yx YAx = A t Y x + YA x Roy Wilhelm Probst A New PC Method for Optimal Power Flow 7 / 20

11 Linear Programming Nonlinear Programming Complete PC (CPC) Complete Predictor-Corrector Method min s.t. x t Hx XAx = b x 0 XA + diag(ax) 0 0 2H + YA + A t Y A t X + diag(ax) I x ỹ z = r 1 r 2 r I ˆx ŷ ẑ = X A X u Ỹ A X u A t Ỹ X u µu X Zu Roy Wilhelm Probst A New PC Method for Optimal Power Flow 8 / 20

12 IEEE 30 bus system Interior-Point Methods Example Model T k = e k + jf k e k : Voltage Real part f k : Voltage Complex part 3 6 P k : Active Power liquid injection Barras de Geração Barras de Carga Q k : Reactive Power liquid injection Roy Wilhelm Probst A New PC Method for Optimal Power Flow 9 / 20

13 Model Interior-Point Methods Example Model min e t Ge + f t Gf s.t. P k (e, f ) + P Ck P Gk = 0 k C Q k (e, f ) + Q Ck Q Gk = 0 k C (vk min ) 2 V k (e, f ) (vk max ) 2 k N Pk min P k (e, f ) Pk max k G Qk min Q k (e, f ) Qk max k R where: P = EGe + FGf + FBe EBf Q = FGe EGf EBe FBf V = Ee + Ff Roy Wilhelm Probst A New PC Method for Optimal Power Flow 10 / 20

14 CPC Corrections Interior-Point Methods Example Model [ x P( x) ˆr 1 = t ỹ p + x P( x) t z 1p x P( x) t z 2p x Q( x) t ỹ q + x Q( x) t z 1q x Q( x) t z 2q ] ˆr 2 = ˆr 4 = S 1v z 1v µu S 1p z 1p µu S 1q z 1q µu V ( x) P( x) Q( x) ˆr 5 = ˆr 3 = V ( x) P( x) Q( x) S 2v z 2v µu S 2p z 2p µu S 2q z 2q µu ˆr 6 = V ( x) P( x) Q( x) Roy Wilhelm Probst A New PC Method for Optimal Power Flow 11 / 20

15 CPC Corrections (cont.) Example Model where: P( x) = ẼG ẽ + F G f + F B ẽ ẼB f Q( x) = F G ẽ ẼG f ẼB ẽ F B f V ( x) = Ẽ ẽ + F f [ ] ẼG + F B + diag(g ẽ B f ) x P( x) = F G ẼB + diag(g f + B ẽ) [ ] F G ẼB diag(g f + B ẽ) x Q( x) = B + diag(g ẽ B f ) Roy Wilhelm Probst A New PC Method for Optimal Power Flow 12 / 20

16 Example Model OPF using Cartesian coordinates Quadratic objective function Quadratic constraints Constant Hessian Nonlinear corrections in all optimality conditions Roy Wilhelm Probst A New PC Method for Optimal Power Flow 13 / 20

17 Test cases Results MATLAB 7.8 (R2009a) Windows XP Intel Core 2 Quad Q9550 2,83 GHz 3,23 GB of RAM Roy Wilhelm Probst A New PC Method for Optimal Power Flow 14 / 20

18 Power systems test cases Test cases Results Buses and lines System N A R C M IEEE IEEE IEEE BRAZIL Roy Wilhelm Probst A New PC Method for Optimal Power Flow 15 / 20

19 Iteration count and CPU time Test cases Results v k [0.90, 1.10] Iterations CPU Time System PD PC CPC PD PC CPC IEEE IEEE IEEE Roy Wilhelm Probst A New PC Method for Optimal Power Flow 16 / 20

20 Iteration count and CPU time Test cases Results v k [0.95, 1.05] Iterations CPU Time System PD PC CPC PD PC CPC IEEE IEEE IEEE Roy Wilhelm Probst A New PC Method for Optimal Power Flow 17 / 20

21 Test cases Results Active constraints v k at optimal solution v k [0.90, 1.10] v k [0.95, 1.05] System v max v min v max v min IEEE IEEE IEEE Roy Wilhelm Probst A New PC Method for Optimal Power Flow 18 / 20

22 Iteration count and CPU time Test cases Results BRAZIL system P max / P C Iterations CPU Time PD PC CPC PD PC CPC Roy Wilhelm Probst A New PC Method for Optimal Power Flow 19 / 20

23 Interior-Point Methods Quadratic objective function and constraints OPF Problem Complete Predictor-Corrector Nonlinear Corrections Robustness and Speed Heuristic Roy Wilhelm Probst A New PC Method for Optimal Power Flow 20 / 20

Part 5: Penalty and augmented Lagrangian methods for equality constrained optimization. Nick Gould (RAL)

Part 5: Penalty and augmented Lagrangian methods for equality constrained optimization. Nick Gould (RAL) Part 5: Penalty and augmented Lagrangian methods for equality constrained optimization Nick Gould (RAL) x IR n f(x) subject to c(x) = Part C course on continuoue optimization CONSTRAINED MINIMIZATION x

More information

A Regularized Interior-Point Method for Constrained Nonlinear Least Squares

A Regularized Interior-Point Method for Constrained Nonlinear Least Squares A Regularized Interior-Point Method for Constrained Nonlinear Least Squares XII Brazilian Workshop on Continuous Optimization Abel Soares Siqueira Federal University of Paraná - Curitiba/PR - Brazil Dominique

More information

Interior Point Methods in Mathematical Programming

Interior Point Methods in Mathematical Programming Interior Point Methods in Mathematical Programming Clóvis C. Gonzaga Federal University of Santa Catarina, Brazil Journées en l honneur de Pierre Huard Paris, novembre 2008 01 00 11 00 000 000 000 000

More information

POWER SYSTEMS in general are currently operating

POWER SYSTEMS in general are currently operating TO APPEAR IN IEEE TRANSACTIONS ON POWER SYSTEMS 1 Robust Optimal Power Flow Solution Using Trust Region and Interior-Point Methods Andréa A. Sousa, Geraldo L. Torres, Member IEEE, Claudio A. Cañizares,

More information

Levy Differential Evolutionary Particle Swarm Optimization (LEVY DEEPSO)

Levy Differential Evolutionary Particle Swarm Optimization (LEVY DEEPSO) 1 Levy Differential Evolutionary Particle Swarm Optimization (LEVY DEEPSO) Developers: Kartik S. Pandya, CHARUSAT-INDIA S.K. Joshi, MSU-INDIA S.N. Singh, IIT-K-INDIA 2 Particle Swarm Optimization[1] Particles:

More information

SF2822 Applied Nonlinear Optimization. Preparatory question. Lecture 10: Interior methods. Anders Forsgren. 1. Try to solve theory question 7.

SF2822 Applied Nonlinear Optimization. Preparatory question. Lecture 10: Interior methods. Anders Forsgren. 1. Try to solve theory question 7. SF2822 Applied Nonlinear Optimization Lecture 10: Interior methods Anders Forsgren SF2822 Applied Nonlinear Optimization, KTH 1 / 24 Lecture 10, 2017/2018 Preparatory question 1. Try to solve theory question

More information

Minimization of Energy Loss using Integrated Evolutionary Approaches

Minimization of Energy Loss using Integrated Evolutionary Approaches Minimization of Energy Loss using Integrated Evolutionary Approaches Attia A. El-Fergany, Member, IEEE, Mahdi El-Arini, Senior Member, IEEE Paper Number: 1569614661 Presentation's Outline Aim of this work,

More information

Constraint Reduction for Linear Programs with Many Constraints

Constraint Reduction for Linear Programs with Many Constraints Constraint Reduction for Linear Programs with Many Constraints André L. Tits Institute for Systems Research and Department of Electrical and Computer Engineering University of Maryland, College Park Pierre-Antoine

More information

Interior Point Methods for Mathematical Programming

Interior Point Methods for Mathematical Programming Interior Point Methods for Mathematical Programming Clóvis C. Gonzaga Federal University of Santa Catarina, Florianópolis, Brazil EURO - 2013 Roma Our heroes Cauchy Newton Lagrange Early results Unconstrained

More information

Operation of the Power Distribution System via Sensitivity Analysis

Operation of the Power Distribution System via Sensitivity Analysis Memorias del XVI Congreso Latinoamericano de Control Automático, CLCA 2014 Operation of the Power Distribution System via Sensitivity Analysis William M. da Rosa, Priscila Rossoni, Julio C. Teixeira e

More information

Interior Point Methods for LP

Interior Point Methods for LP 11.1 Interior Point Methods for LP Katta G. Murty, IOE 510, LP, U. Of Michigan, Ann Arbor, Winter 1997. Simplex Method - A Boundary Method: Starting at an extreme point of the feasible set, the simplex

More information

Entropy Enhanced Covariance Matrix Adaptation Evolution Strategy (EE_CMAES)

Entropy Enhanced Covariance Matrix Adaptation Evolution Strategy (EE_CMAES) 1 Entropy Enhanced Covariance Matrix Adaptation Evolution Strategy (EE_CMAES) Developers: Main Author: Kartik Pandya, Dept. of Electrical Engg., CSPIT, CHARUSAT, Changa, India Co-Author: Jigar Sarda, Dept.

More information

Interior-point based algorithms for the solution of optimal power flow problems

Interior-point based algorithms for the solution of optimal power flow problems Electric Power Systems Research 77 (2007) 508 517 Interior-point based algorithms for the solution of optimal power flow problems Florin Capitanescu, Mevludin Glavic, Damien Ernst, Louis Wehenkel Department

More information

Development of the new MINLP Solver Decogo using SCIP - Status Report

Development of the new MINLP Solver Decogo using SCIP - Status Report Development of the new MINLP Solver Decogo using SCIP - Status Report Pavlo Muts with Norman Breitfeld, Vitali Gintner, Ivo Nowak SCIP Workshop 2018, Aachen Table of contents 1. Introduction 2. Automatic

More information

An Active Set Strategy for Solving Optimization Problems with up to 200,000,000 Nonlinear Constraints

An Active Set Strategy for Solving Optimization Problems with up to 200,000,000 Nonlinear Constraints An Active Set Strategy for Solving Optimization Problems with up to 200,000,000 Nonlinear Constraints Klaus Schittkowski Department of Computer Science, University of Bayreuth 95440 Bayreuth, Germany e-mail:

More information

Weighted Least Squares Topology Error Detection And Identification

Weighted Least Squares Topology Error Detection And Identification Weighted Least Squares Topology Error Detection And Identification A THESIS SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY Jason Glen Lindquist IN PARTIAL FULFILLMENT

More information

CS711008Z Algorithm Design and Analysis

CS711008Z Algorithm Design and Analysis CS711008Z Algorithm Design and Analysis Lecture 8 Linear programming: interior point method Dongbo Bu Institute of Computing Technology Chinese Academy of Sciences, Beijing, China 1 / 31 Outline Brief

More information

What s New in Active-Set Methods for Nonlinear Optimization?

What s New in Active-Set Methods for Nonlinear Optimization? What s New in Active-Set Methods for Nonlinear Optimization? Philip E. Gill Advances in Numerical Computation, Manchester University, July 5, 2011 A Workshop in Honor of Sven Hammarling UCSD Center for

More information

SMO vs PDCO for SVM: Sequential Minimal Optimization vs Primal-Dual interior method for Convex Objectives for Support Vector Machines

SMO vs PDCO for SVM: Sequential Minimal Optimization vs Primal-Dual interior method for Convex Objectives for Support Vector Machines vs for SVM: Sequential Minimal Optimization vs Primal-Dual interior method for Convex Objectives for Support Vector Machines Ding Ma Michael Saunders Working paper, January 5 Introduction In machine learning,

More information

Lecture 10. Primal-Dual Interior Point Method for LP

Lecture 10. Primal-Dual Interior Point Method for LP IE 8534 1 Lecture 10. Primal-Dual Interior Point Method for LP IE 8534 2 Consider a linear program (P ) minimize c T x subject to Ax = b x 0 and its dual (D) maximize b T y subject to A T y + s = c s 0.

More information

Mohd Jamil Khan and 2 Yogesh Kumar. Churu, Raj., India.

Mohd Jamil Khan and 2 Yogesh Kumar. Churu, Raj., India. International Journal of Mathematics Research. ISSN 0976-5840 Volume 8, Number 3 (2016), pp. 251-263 International Research Publication House http://www.irphouse.com Optimal Power Flow (OPF) formulation

More information

A preventive control model for static voltage stability and thermal stability based on power transfer capabilities of weak branches

A preventive control model for static voltage stability and thermal stability based on power transfer capabilities of weak branches A model for static voltage stability and thermal stability based on power transfer capabilities of weak branches Wei Yan 1, Xiuqiong Hu 1, Juan Yu 1 and Wenyuan Li 1,2 (1)The State Key Laboratory of Power

More information

AC Power Flows, Generalized OPF Costs and their Derivatives using Complex Matrix Notation

AC Power Flows, Generalized OPF Costs and their Derivatives using Complex Matrix Notation AC Power Flows, Generalized OPF Costs and their Derivatives using Complex Matrix Notation Ray D. Zimmerman February 24, 2010 Matpower Technical Note 2 c 2008, 2010, 2011, 2017, 2018 Power Systems Engineering

More information

Optimal Tap Settings for Voltage Regulation Transformers in Distribution Networks

Optimal Tap Settings for Voltage Regulation Transformers in Distribution Networks Optimal Tap Settings for Voltage Regulation Transformers in Distribution Networks Brett A. Robbins Department of Electrical and Computer Engineering University of Illinois at Urbana-Champaign May 9, 2014

More information

CHAPTER 6 : LITERATURE REVIEW

CHAPTER 6 : LITERATURE REVIEW CHAPTER 6 : LITERATURE REVIEW Chapter : LITERATURE REVIEW 77 M E A S U R I N G T H E E F F I C I E N C Y O F D E C I S I O N M A K I N G U N I T S A B S T R A C T A n o n l i n e a r ( n o n c o n v e

More information

P E R E N C O - C H R I S T M A S P A R T Y

P E R E N C O - C H R I S T M A S P A R T Y L E T T I C E L E T T I C E I S A F A M I L Y R U N C O M P A N Y S P A N N I N G T W O G E N E R A T I O N S A N D T H R E E D E C A D E S. B A S E D I N L O N D O N, W E H A V E T H E P E R F E C T R

More information

SF2822 Applied nonlinear optimization, final exam Saturday December

SF2822 Applied nonlinear optimization, final exam Saturday December SF2822 Applied nonlinear optimization, final exam Saturday December 5 27 8. 3. Examiner: Anders Forsgren, tel. 79 7 27. Allowed tools: Pen/pencil, ruler and rubber; plus a calculator provided by the department.

More information

THE restructuring of the power industry has lead to

THE restructuring of the power industry has lead to GLOBALLY CONVERGENT OPTIMAL POWER FLOW USING COMPLEMENTARITY FUNCTIONS AND TRUST REGION METHODS Geraldo L. Torres Universidade Federal de Pernambuco Recife, Brazil gltorres@ieee.org Abstract - As power

More information

Marcus Pantoja da Silva 1 and Celso Pascoli Bottura 2. Abstract: Nonlinear systems with time-varying uncertainties

Marcus Pantoja da Silva 1 and Celso Pascoli Bottura 2. Abstract: Nonlinear systems with time-varying uncertainties A NEW PROPOSAL FOR H NORM CHARACTERIZATION AND THE OPTIMAL H CONTROL OF NONLINEAR SSTEMS WITH TIME-VARING UNCERTAINTIES WITH KNOWN NORM BOUND AND EXOGENOUS DISTURBANCES Marcus Pantoja da Silva 1 and Celso

More information

Efficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph

Efficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph Efficient robust optimization for robust control with constraints p. 1 Efficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph Efficient robust optimization

More information

Fixed interval scheduling problems under uncertainty

Fixed interval scheduling problems under uncertainty Fixed interval scheduling problems under uncertainty Martin Branda et al. Charles University Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics Stochastic Programming

More information

Algorithms for nonlinear programming problems II

Algorithms for nonlinear programming problems II Algorithms for nonlinear programming problems II Martin Branda Charles University in Prague Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics Computational Aspects

More information

Optimization for Communications and Networks. Poompat Saengudomlert. Session 4 Duality and Lagrange Multipliers

Optimization for Communications and Networks. Poompat Saengudomlert. Session 4 Duality and Lagrange Multipliers Optimization for Communications and Networks Poompat Saengudomlert Session 4 Duality and Lagrange Multipliers P Saengudomlert (2015) Optimization Session 4 1 / 14 24 Dual Problems Consider a primal convex

More information

Power System Reactive Power Optimization Based on Fuzzy Formulation and Interior Point Filter Algorithm

Power System Reactive Power Optimization Based on Fuzzy Formulation and Interior Point Filter Algorithm Energy and Power Engineering, 203, 5, 693-697 doi:0.4236/ee.203.54b34 Published Online July 203 (htt://www.scir.org/journal/ee Power System Reactive Power Otimization Based on Fuzzy Formulation and Interior

More information

Primal-dual IPM with Asymmetric Barrier

Primal-dual IPM with Asymmetric Barrier Primal-dual IPM with Asymmetric Barrier Yurii Nesterov, CORE/INMA (UCL) September 29, 2008 (IFOR, ETHZ) Yu. Nesterov Primal-dual IPM with Asymmetric Barrier 1/28 Outline 1 Symmetric and asymmetric barriers

More information

Relaxed linearized algorithms for faster X-ray CT image reconstruction

Relaxed linearized algorithms for faster X-ray CT image reconstruction Relaxed linearized algorithms for faster X-ray CT image reconstruction Hung Nien and Jeffrey A. Fessler University of Michigan, Ann Arbor The 13th Fully 3D Meeting June 2, 2015 1/20 Statistical image reconstruction

More information

A Load-Shedding Model Based on Sensitivity Analysis in on-line Power System Operation Risk Assessment

A Load-Shedding Model Based on Sensitivity Analysis in on-line Power System Operation Risk Assessment energies Article A Load-Shedding Model Based on Sensitivity Analysis in on-line Power System Operation Risk Assessment Zhe Zhang 1, Hang Yang 1, *, Xianggen Yin 1, Jiexiang Han 1, Yong Wang 2 and Guoyan

More information

Lecture 13: Constrained optimization

Lecture 13: Constrained optimization 2010-12-03 Basic ideas A nonlinearly constrained problem must somehow be converted relaxed into a problem which we can solve (a linear/quadratic or unconstrained problem) We solve a sequence of such problems

More information

A preventive control model for static voltage stability considering the active power transfer capabilities of weak branches

A preventive control model for static voltage stability considering the active power transfer capabilities of weak branches A preventive control model for static voltage stability considering the active power transfer capabilities of weak branches Xiuqiong Hu 1,2, Zhouyang Ren 2, Yiming Li 2 (1) School of Information and Electric

More information

AC Power Flows and their Derivatives using Complex Matrix Notation and Cartesian Coordinate Voltages

AC Power Flows and their Derivatives using Complex Matrix Notation and Cartesian Coordinate Voltages AC Power Flows and their Derivatives using Complex Matrix Notation and Cartesian Coordinate Voltages Baljinnyam Sereeter Ray D. Zimmerman April 2, 2018 Matpower echnical Note 4 c 2008, 2010, 2011, 2017,

More information

Nonlinear Programming

Nonlinear Programming Nonlinear Programming Kees Roos e-mail: C.Roos@ewi.tudelft.nl URL: http://www.isa.ewi.tudelft.nl/ roos LNMB Course De Uithof, Utrecht February 6 - May 8, A.D. 2006 Optimization Group 1 Outline for week

More information

18. Primal-dual interior-point methods

18. Primal-dual interior-point methods L. Vandenberghe EE236C (Spring 213-14) 18. Primal-dual interior-point methods primal-dual central path equations infeasible primal-dual method primal-dual method for self-dual embedding 18-1 Symmetric

More information

Branch Flow Model. Computing + Math Sciences Electrical Engineering

Branch Flow Model. Computing + Math Sciences Electrical Engineering Branch Flow Model Masoud Farivar Steven Low Computing + Math Sciences Electrical Engineering Arpil 014 TPS paper Farivar and Low Branch flow model: relaxations and convexification IEEE Trans. Power Systems,

More information

Algorithms for constrained local optimization

Algorithms for constrained local optimization Algorithms for constrained local optimization Fabio Schoen 2008 http://gol.dsi.unifi.it/users/schoen Algorithms for constrained local optimization p. Feasible direction methods Algorithms for constrained

More information

Experience in Factoring Large Integers Using Quadratic Sieve

Experience in Factoring Large Integers Using Quadratic Sieve Experience in Factoring Large Integers Using Quadratic Sieve D. J. Guan Department of Computer Science, National Sun Yat-Sen University, Kaohsiung, Taiwan 80424 guan@cse.nsysu.edu.tw April 19, 2005 Abstract

More information

Numerical Methods for Model Predictive Control. Jing Yang

Numerical Methods for Model Predictive Control. Jing Yang Numerical Methods for Model Predictive Control Jing Yang Kongens Lyngby February 26, 2008 Technical University of Denmark Informatics and Mathematical Modelling Building 321, DK-2800 Kongens Lyngby, Denmark

More information

Nonlinear optimization

Nonlinear optimization Nonlinear optimization Anders Forsgren Optimization and Systems Theory Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden evita Winter School 2009 Geilo, Norway January 11

More information

Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma

Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent

More information

A solution approach for linear optimization with completely positive matrices

A solution approach for linear optimization with completely positive matrices A solution approach for linear optimization with completely positive matrices Franz Rendl http://www.math.uni-klu.ac.at Alpen-Adria-Universität Klagenfurt Austria joint work with M. Bomze (Wien) and F.

More information

Agenda. Interior Point Methods. 1 Barrier functions. 2 Analytic center. 3 Central path. 4 Barrier method. 5 Primal-dual path following algorithms

Agenda. Interior Point Methods. 1 Barrier functions. 2 Analytic center. 3 Central path. 4 Barrier method. 5 Primal-dual path following algorithms Agenda Interior Point Methods 1 Barrier functions 2 Analytic center 3 Central path 4 Barrier method 5 Primal-dual path following algorithms 6 Nesterov Todd scaling 7 Complexity analysis Interior point

More information

Transient Finite Element Analysis of a Spice-Coupled Transformer with COMSOL-Multiphysics

Transient Finite Element Analysis of a Spice-Coupled Transformer with COMSOL-Multiphysics Presented at the COMSOL Conference 2010 Paris Thomas Bödrich 1, Holger Neubert 1, Rolf Disselnkötter 2 Transient Finite Element Analysis of a Spice-Coupled Transformer with COMSOL-Multiphysics 2010-11-17

More information

AUGMENTED LAGRANGIAN METHODS FOR CONVEX OPTIMIZATION

AUGMENTED LAGRANGIAN METHODS FOR CONVEX OPTIMIZATION New Horizon of Mathematical Sciences RIKEN ithes - Tohoku AIMR - Tokyo IIS Joint Symposium April 28 2016 AUGMENTED LAGRANGIAN METHODS FOR CONVEX OPTIMIZATION T. Takeuchi Aihara Lab. Laboratories for Mathematics,

More information

A reduced basis method and ROM-based optimization for batch chromatography

A reduced basis method and ROM-based optimization for batch chromatography MAX PLANCK INSTITUTE MedRed, 2013 Magdeburg, Dec 11-13 A reduced basis method and ROM-based optimization for batch chromatography Yongjin Zhang joint work with Peter Benner, Lihong Feng and Suzhou Li Max

More information

Lecture 3. Optimization Problems and Iterative Algorithms

Lecture 3. Optimization Problems and Iterative Algorithms Lecture 3 Optimization Problems and Iterative Algorithms January 13, 2016 This material was jointly developed with Angelia Nedić at UIUC for IE 598ns Outline Special Functions: Linear, Quadratic, Convex

More information

2.098/6.255/ Optimization Methods Practice True/False Questions

2.098/6.255/ Optimization Methods Practice True/False Questions 2.098/6.255/15.093 Optimization Methods Practice True/False Questions December 11, 2009 Part I For each one of the statements below, state whether it is true or false. Include a 1-3 line supporting sentence

More information

Newton s Method. Javier Peña Convex Optimization /36-725

Newton s Method. Javier Peña Convex Optimization /36-725 Newton s Method Javier Peña Convex Optimization 10-725/36-725 1 Last time: dual correspondences Given a function f : R n R, we define its conjugate f : R n R, f ( (y) = max y T x f(x) ) x Properties and

More information

Infeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming*

Infeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming* Infeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming* Notes for CAAM 564, Spring 2012 Dept of CAAM, Rice University Outline (1) Introduction (2) Formulation & a complexity

More information

Interior Point Methods. We ll discuss linear programming first, followed by three nonlinear problems. Algorithms for Linear Programming Problems

Interior Point Methods. We ll discuss linear programming first, followed by three nonlinear problems. Algorithms for Linear Programming Problems AMSC 607 / CMSC 764 Advanced Numerical Optimization Fall 2008 UNIT 3: Constrained Optimization PART 4: Introduction to Interior Point Methods Dianne P. O Leary c 2008 Interior Point Methods We ll discuss

More information

Semidefinite Relaxation of Optimal Power Flow for ac-dc Grids

Semidefinite Relaxation of Optimal Power Flow for ac-dc Grids Semidefinite Relaxation of Optimal Power Flow for ac-dc Grids Shahab Bahrami, Student Member, IEEE, Francis Therrien, Member, IEEE, Vincent W.S. Wong, Fellow, IEEE, and Juri Jatsevich, Senior Member, IEEE

More information

DELFT UNIVERSITY OF TECHNOLOGY

DELFT UNIVERSITY OF TECHNOLOGY DELFT UNIVERSITY OF TECHNOLOGY REPORT -09 Computational and Sensitivity Aspects of Eigenvalue-Based Methods for the Large-Scale Trust-Region Subproblem Marielba Rojas, Bjørn H. Fotland, and Trond Steihaug

More information

ECE580 Solution to Problem Set 6

ECE580 Solution to Problem Set 6 ECE580 Fall 2015 Solution to Problem Set 6 December 23 2015 1 ECE580 Solution to Problem Set 6 These problems are from the textbook by Chong and Zak 4th edition which is the textbook for the ECE580 Fall

More information

Optimal Control. Lecture 20. Control Lyapunov Function, Optimal Estimation. John T. Wen. April 5. Ref: Papers by R. Freeman (on-line), B&H Ch.

Optimal Control. Lecture 20. Control Lyapunov Function, Optimal Estimation. John T. Wen. April 5. Ref: Papers by R. Freeman (on-line), B&H Ch. Optimal Control Lecture 20 Control Lyapunov Function, Optimal Estimation John T. Wen April 5 Ref: Papers by R. Freeman (on-line), B&H Ch. 12 Outline Summary of Control Lyapunov Function and example Introduction

More information

Solving large Semidefinite Programs - Part 1 and 2

Solving large Semidefinite Programs - Part 1 and 2 Solving large Semidefinite Programs - Part 1 and 2 Franz Rendl http://www.math.uni-klu.ac.at Alpen-Adria-Universität Klagenfurt Austria F. Rendl, Singapore workshop 2006 p.1/34 Overview Limits of Interior

More information

Preliminaries Overview OPF and Extensions. Convex Optimization. Lecture 8 - Applications in Smart Grids. Instructor: Yuanzhang Xiao

Preliminaries Overview OPF and Extensions. Convex Optimization. Lecture 8 - Applications in Smart Grids. Instructor: Yuanzhang Xiao Convex Optimization Lecture 8 - Applications in Smart Grids Instructor: Yuanzhang Xiao University of Hawaii at Manoa Fall 2017 1 / 32 Today s Lecture 1 Generalized Inequalities and Semidefinite Programming

More information

Convex Relaxations for Optimization of AC and HVDC Grids under Uncertainty

Convex Relaxations for Optimization of AC and HVDC Grids under Uncertainty P L Power Systems Laboratory Center for Electric Power and Energy Department of Electrical Engineering Andreas Horst Venzke Convex Relaxations for Optimization of AC and HVDC Grids under Uncertainty Master

More information

Convex Optimization on Large-Scale Domains Given by Linear Minimization Oracles

Convex Optimization on Large-Scale Domains Given by Linear Minimization Oracles Convex Optimization on Large-Scale Domains Given by Linear Minimization Oracles Arkadi Nemirovski H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology Joint research

More information

Open Access Target Tracking Algorithm Based on Improved Unscented Kalman Filter

Open Access Target Tracking Algorithm Based on Improved Unscented Kalman Filter Send Orders for Reprints to reprints@benthamscience.ae he Open Automation and Control Systems Journal, 2015, 7, 991-995 991 Open Access arget racking Algorithm Based on Improved Unscented Kalman Filter

More information

MATH 104: INTRODUCTORY ANALYSIS SPRING 2008/09 PROBLEM SET 8 SOLUTIONS

MATH 104: INTRODUCTORY ANALYSIS SPRING 2008/09 PROBLEM SET 8 SOLUTIONS MATH 04: INTRODUCTORY ANALYSIS SPRING 008/09 PROBLEM SET 8 SOLUTIONS. Let f : R R be continuous periodic with period, i.e. f(x + ) = f(x) for all x R. Prove the following: (a) f is bounded above below

More information

Incremental Latin Hypercube Sampling

Incremental Latin Hypercube Sampling Incremental Latin Hypercube Sampling for Lifetime Stochastic Behavioral Modeling of Analog Circuits Yen-Lung Chen +, Wei Wu *, Chien-Nan Jimmy Liu + and Lei He * EE Dept., National Central University,

More information

EE613 Machine Learning for Engineers. Kernel methods Support Vector Machines. jean-marc odobez 2015

EE613 Machine Learning for Engineers. Kernel methods Support Vector Machines. jean-marc odobez 2015 EE613 Machine Learning for Engineers Kernel methods Support Vector Machines jean-marc odobez 2015 overview Kernel methods introductions and main elements defining kernels Kernelization of k-nn, K-Means,

More information

A smoothing Newton-type method for second-order cone programming problems based on a new smoothing Fischer-Burmeister function

A smoothing Newton-type method for second-order cone programming problems based on a new smoothing Fischer-Burmeister function Volume 30, N. 3, pp. 569 588, 2011 Copyright 2011 SBMAC ISSN 0101-8205 www.scielo.br/cam A smoothing Newton-type method for second-order cone programming problems based on a new smoothing Fischer-Burmeister

More information

Computational Analysis of Sparsity-Exploiting Moment Relaxations of the OPF Problem

Computational Analysis of Sparsity-Exploiting Moment Relaxations of the OPF Problem Computational Analysis of Sparsity-Exploiting Moment Relaxations of the OPF Problem Daniel K. Molzahn, Ian A. Hiskens EECS Department, University of Michigan Ann Arbor, MI USA molzahn,hiskens@umich.edu

More information

Krylov Subspace Methods for Nonlinear Model Reduction

Krylov Subspace Methods for Nonlinear Model Reduction MAX PLANCK INSTITUT Conference in honour of Nancy Nichols 70th birthday Reading, 2 3 July 2012 Krylov Subspace Methods for Nonlinear Model Reduction Peter Benner and Tobias Breiten Max Planck Institute

More information

A DC Power Flow Extension

A DC Power Flow Extension 2013 4th IEEE PES Innovative Smart Grid Technologies Europe (ISGT Europe), October 6-9, Copenhagen 1 A DC Power Flow Extension Theodoros Kyriakidis, Rachid Cherkaoui, Maher Kayal Electronics Laboratory

More information

ITERATIVE techniques based on the Newton-Raphson algorithm

ITERATIVE techniques based on the Newton-Raphson algorithm Jacobian Singularities in Optimal Power Flow Problems Caused by Intertemporal Constraints Kyri Baker, Student Member, IEEE, Dinghuan Zhu, Student Member, IEEE, Gabriela Hug, Member, IEEE, and Xin Li, Senior

More information

Nonsymmetric potential-reduction methods for general cones

Nonsymmetric potential-reduction methods for general cones CORE DISCUSSION PAPER 2006/34 Nonsymmetric potential-reduction methods for general cones Yu. Nesterov March 28, 2006 Abstract In this paper we propose two new nonsymmetric primal-dual potential-reduction

More information

THE solution of the absolute value equation (AVE) of

THE solution of the absolute value equation (AVE) of The nonlinear HSS-like iterative method for absolute value equations Mu-Zheng Zhu Member, IAENG, and Ya-E Qi arxiv:1403.7013v4 [math.na] 2 Jan 2018 Abstract Salkuyeh proposed the Picard-HSS iteration method

More information

A Constraint-Reduced MPC Algorithm for Convex Quadratic Programming, with a Modified Active-Set Identification Scheme

A Constraint-Reduced MPC Algorithm for Convex Quadratic Programming, with a Modified Active-Set Identification Scheme A Constraint-Reduced MPC Algorithm for Convex Quadratic Programming, with a Modified Active-Set Identification Scheme M. Paul Laiu 1 and (presenter) André L. Tits 2 1 Oak Ridge National Laboratory laiump@ornl.gov

More information

Fast ADMM for Sum of Squares Programs Using Partial Orthogonality

Fast ADMM for Sum of Squares Programs Using Partial Orthogonality Fast ADMM for Sum of Squares Programs Using Partial Orthogonality Antonis Papachristodoulou Department of Engineering Science University of Oxford www.eng.ox.ac.uk/control/sysos antonis@eng.ox.ac.uk with

More information

Research overview. Seminar September 4, Lehigh University Department of Industrial & Systems Engineering. Research overview.

Research overview. Seminar September 4, Lehigh University Department of Industrial & Systems Engineering. Research overview. Research overview Lehigh University Department of Industrial & Systems Engineering COR@L Seminar September 4, 2008 1 Duality without regularity condition Duality in non-exact arithmetic 2 interior point

More information

SF2822 Applied nonlinear optimization, final exam Wednesday June

SF2822 Applied nonlinear optimization, final exam Wednesday June SF2822 Applied nonlinear optimization, final exam Wednesday June 3 205 4.00 9.00 Examiner: Anders Forsgren, tel. 08-790 7 27. Allowed tools: Pen/pencil, ruler and eraser. Note! Calculator is not allowed.

More information

Nonlinear Optimization: What s important?

Nonlinear Optimization: What s important? Nonlinear Optimization: What s important? Julian Hall 10th May 2012 Convexity: convex problems A local minimizer is a global minimizer A solution of f (x) = 0 (stationary point) is a minimizer A global

More information

arxiv: v1 [math.na] 21 Oct 2014

arxiv: v1 [math.na] 21 Oct 2014 Computing Symmetric Positive Definite Solutions of Three Types of Nonlinear Matrix Equations arxiv:1410.5559v1 [math.na] 21 Oct 2014 Negin Bagherpour a, Nezam Mahdavi-Amiri a, a Department of Mathematical

More information

Algorithms for nonlinear programming problems II

Algorithms for nonlinear programming problems II Algorithms for nonlinear programming problems II Martin Branda Charles University Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics Computational Aspects of Optimization

More information

Numerical Optimization

Numerical Optimization Linear Programming - Interior Point Methods Computer Science and Automation Indian Institute of Science Bangalore 560 012, India. NPTEL Course on Example 1 Computational Complexity of Simplex Algorithm

More information

Module 04 Optimization Problems KKT Conditions & Solvers

Module 04 Optimization Problems KKT Conditions & Solvers Module 04 Optimization Problems KKT Conditions & Solvers Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September

More information

2-4 Zeros of Polynomial Functions

2-4 Zeros of Polynomial Functions Write a polynomial function of least degree with real coefficients in standard form that has the given zeros. 33. 2, 4, 3, 5 Using the Linear Factorization Theorem and the zeros 2, 4, 3, and 5, write f

More information

Copositive Programming and Combinatorial Optimization

Copositive Programming and Combinatorial Optimization Copositive Programming and Combinatorial Optimization Franz Rendl http://www.math.uni-klu.ac.at Alpen-Adria-Universität Klagenfurt Austria joint work with I.M. Bomze (Wien) and F. Jarre (Düsseldorf) IMA

More information

Continuous methods for numerical linear algebra problems

Continuous methods for numerical linear algebra problems Continuous methods for numerical linear algebra problems Li-Zhi Liao (http://www.math.hkbu.edu.hk/ liliao) Department of Mathematics Hong Kong Baptist University The First International Summer School on

More information

Infeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming*

Infeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming* Infeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming* Yin Zhang Dept of CAAM, Rice University Outline (1) Introduction (2) Formulation & a complexity theorem (3)

More information

Newton s Method. Ryan Tibshirani Convex Optimization /36-725

Newton s Method. Ryan Tibshirani Convex Optimization /36-725 Newton s Method Ryan Tibshirani Convex Optimization 10-725/36-725 1 Last time: dual correspondences Given a function f : R n R, we define its conjugate f : R n R, Properties and examples: f (y) = max x

More information

Alternating Direction Augmented Lagrangian Algorithms for Convex Optimization

Alternating Direction Augmented Lagrangian Algorithms for Convex Optimization Alternating Direction Augmented Lagrangian Algorithms for Convex Optimization Joint work with Bo Huang, Shiqian Ma, Tony Qin, Katya Scheinberg, Zaiwen Wen and Wotao Yin Department of IEOR, Columbia University

More information

14. Nonlinear equations

14. Nonlinear equations L. Vandenberghe ECE133A (Winter 2018) 14. Nonlinear equations Newton method for nonlinear equations damped Newton method for unconstrained minimization Newton method for nonlinear least squares 14-1 Set

More information

Real Time Voltage Control using Genetic Algorithm

Real Time Voltage Control using Genetic Algorithm Real Time Voltage Control using Genetic Algorithm P. Thirusenthil kumaran, C. Kamalakannan Department of EEE, Rajalakshmi Engineering College, Chennai, India Abstract An algorithm for control action selection

More information

Sparse Polynomial Optimization for Urban Distribution Networks

Sparse Polynomial Optimization for Urban Distribution Networks Sparse Polynomial Optimization for Urban Distribution Networks Martin Mevissen IBM Research - Ireland MINO/ COST Spring School on Optimization, Tilburg, March 11, 2015 The Team Brad Eck Bissan Ghaddar

More information

Lecture 8. Strong Duality Results. September 22, 2008

Lecture 8. Strong Duality Results. September 22, 2008 Strong Duality Results September 22, 2008 Outline Lecture 8 Slater Condition and its Variations Convex Objective with Linear Inequality Constraints Quadratic Objective over Quadratic Constraints Representation

More information

Parallel computing of multi-contingency optimal power flow with transient stability constraints

Parallel computing of multi-contingency optimal power flow with transient stability constraints Yang et al. Protection and Control of Modern Power Systems (2018)3:20 https://doi.org/10.1186/s41601-018-0095-z Protection and Control of Modern Power Systems ORIGINAL RESEARCH Parallel computing of multi-contingency

More information

AM 205: lecture 19. Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods

AM 205: lecture 19. Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods AM 205: lecture 19 Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods Quasi-Newton Methods General form of quasi-newton methods: x k+1 = x k α

More information

Local Analysis of the Feasible Primal-Dual Interior-Point Method

Local Analysis of the Feasible Primal-Dual Interior-Point Method Local Analysis of the Feasible Primal-Dual Interior-Point Method R. Silva J. Soares L. N. Vicente Abstract In this paper we analyze the rate of local convergence of the Newton primal-dual interiorpoint

More information

Interior-point methods Optimization Geoff Gordon Ryan Tibshirani

Interior-point methods Optimization Geoff Gordon Ryan Tibshirani Interior-point methods 10-725 Optimization Geoff Gordon Ryan Tibshirani SVM duality Review min v T v/2 + 1 T s s.t. Av yd + s 1 0 s 0 max 1 T α α T Kα/2 s.t. y T α = 0 0 α 1 Gram matrix K Interpretation

More information