A New Predictor-Corrector Method for Optimal Power Flow
|
|
- Virgil Stevens
- 5 years ago
- Views:
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) x IR n f(x) subject to c(x) = Part C course on continuoue optimization CONSTRAINED MINIMIZATION x
More informationA 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 informationInterior 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 informationPOWER 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 informationLevy 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 informationSF2822 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 informationMinimization 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 informationConstraint 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 informationInterior 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 informationOperation 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 informationInterior 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 informationEntropy 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 informationInterior-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 informationDevelopment 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 informationAn 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 informationWeighted 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 informationCS711008Z 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 informationWhat 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 informationSMO 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 informationLecture 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 informationMohd 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 informationA 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 informationAC 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 informationOptimal 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 informationCHAPTER 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 informationP 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 informationSF2822 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 informationTHE 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 informationMarcus 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 informationEfficient 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 informationFixed 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 informationAlgorithms 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 informationOptimization 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 informationPower 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 informationPrimal-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 informationRelaxed 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 informationA 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 informationLecture 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 informationA 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 informationAC 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 informationNonlinear 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 information18. 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 informationBranch 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 informationAlgorithms 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 informationExperience 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 informationNumerical 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 informationNonlinear 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 informationFast 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 informationA 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 informationAgenda. 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 informationTransient 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 informationAUGMENTED 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 informationA 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 informationLecture 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 information2.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 informationNewton 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 informationInfeasible 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 informationInterior 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 informationSemidefinite 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 informationDELFT 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 informationECE580 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 informationOptimal 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 informationSolving 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 informationPreliminaries 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 informationConvex 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 informationConvex 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 informationOpen 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 informationMATH 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 informationIncremental 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 informationEE613 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 informationA 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 informationComputational 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 informationKrylov 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 informationA 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 informationITERATIVE 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 informationNonsymmetric 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 informationTHE 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 informationA 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 informationFast 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 informationResearch 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 informationSF2822 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 informationNonlinear 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 informationarxiv: 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 informationAlgorithms 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 informationNumerical 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 informationModule 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 information2-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 informationCopositive 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 informationContinuous 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 informationInfeasible 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 informationNewton 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 informationAlternating 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 information14. 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 informationReal 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 informationSparse 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 informationLecture 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 informationParallel 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 informationAM 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 informationLocal 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 informationInterior-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