Model order reduction of electrical circuits with nonlinear elements
|
|
- Denis Bennett
- 5 years ago
- Views:
Transcription
1 Model order reduction of electrical circuits with nonlinear elements Andreas Steinbrecher and Tatjana Stykel 1 Introduction The efficient and robust numerical simulation of electrical circuits plays a major role in computer aided design of electronic devices. While the structural size of such devices is decreasing, the complexity of the electrical circuits is increasing. This usually leads to a system of model equations in form of differential-algebraic equations (DAE) with a huge number of unknowns. Simulation of such models is unacceptably time and storage consuming. Model order reduction presents a way out of this dilemma. A general idea of model reduction is to replace a large-scale system by a much smaller model which approximates the input-output relation of the original system within a required accuracy. While a large variety of model reduction techniques exists for linear networks, e.g., [1, 3, 4, model reduction of nonlinear circuits is only in its infancy. In [2, model reduction of nonlinear circuit with only nonlinear resistors was considered. In this paper, we extend these results to more general circuits that may contain other nonlinear elements like nonlinear capacitors or inductors. 2 Circuit Equations A commonly used tool for modeling electrical circuits is the Modified Nodal Analysis (MNA). An electrical circuit can be modeled as a directed graph whose nodes correspond to the nodes of the circuit and whose branches correspond to the circuit elements. Using Kirchhoff s laws as well as the branch constitutive relations, the dynamics of an electrical circuit can be described by a DAE system of the form E (x) d dt x = A x+ f(x)+b u, y = BT x, (1a) Andreas Steinbrecher and Tatjana Stykel, Institut für Mathematik, TU Berlin, Straße des 17. Juni 136, 1623 Berlin, Germany, anst@math.tu-berlin.de, stykel@math.tu-berlin.de 1
2 2 Andreas Steinbrecher and Tatjana Stykel with A C C (A T C η)at C A L A V η E (x)= L (ι L ), A = A T L, x= ι L, (1b) A T V ι V A R g(a T R η) A I f(x)=, B= ιi ui, u=[, y= [. (1c) u I V ι V Here, x, u and y are state, input and output, respectively, η is the vector of node potentials, ι L, ι V and ι I are the vectors of currents through inductors, voltage sources and current sources, respectively, u V and u I are the vectors of voltages of voltage sources and current sources, respectively. We will distinguish between linear circuit elements (denoted by a bar) that are characterized by linear current-voltage relations and nonlinear circuit components (denoted by a tilde) that are characterized by nonlinear current-voltage relations. Without loss of generality, we assume that the circuit elements are ordered such that the incidence matrices describing the circuit topology have the form A C = [A C A C R n η,n C +n C, A L = [A L A L R n η,n L +n L, A R = [A R A R R n η,n R +n R, AV R n η,n V and A I R n η,n I, where the incidence matrices A C, A L and A R correspond to the linear circuit components, and A C, A L and A R correspond to the nonlinear circuit components. Furthermore, the conductance matrix-valued functionc : R n C R n C,n C, the inductance matrix-valued function L : R n L R n L,n L and the resistor relation g : R n R R n R given by [ C C (A T C η) = C (AT C η),l (ι L ) = [ L L (ι L ) [Ḡ A T R, g(a T R η) = η g(a T R η) characterize the physical properties of the capacitors, inductors and resistors, respectively. Here ι L R n L is the vector of currents through the nonlinear inductors. We will assume that the matrices A V and [A C A L A R A V have full rank, the matricesc (A T C η) andl(ι L ) are symmetric, positive definite and the function g(at R η) is monotonically increasing for all admissible η and ι L. This assumptions imply that the circuit elements do not generate energy, i.e., the circuit is passive. (1d) 3 Model Reduction for Nonlinear Circuits In this section, we present a model reduction technique for nonlinear circuits with a small number of nonlinear elements. This technique is based on decoupling of the linear and nonlinear subcircuits in a suitable way, reduction of the linear part using the PABTEC method [4 followed by an adequate recoupling of the unchanged nonlinear part and the reduced linear part to obtain a nonlinear reduced-order model.
3 Model order reduction of electrical circuits with nonlinear elements Decoupling of Linear and Nonlinear Subcircuits Consider the circuit equations (1). Let A R = A1 R + A2 R be decomposed such that A 1 R {,1}n η,n R and A 2 R { 1,} n η,n R and let G 1,G 2 R n R,n R be given such that G 1 and G 1 are both symmetric, positive definite and G 2 is positive semidefinite. Assume that u C R n C and ι z R n R satisfy the relations u C = A T C η and ι z = (G 1 )G 1 1 g(at R η) G 2AT R η. (2) Then system (1) together with the relations η z = (G 1 ) 1 ((G 1 (A 1 R )T G 2 (A 2 R )T )η ι z ), ι C = C (u C ) d dt u C (3a) (3b) for the additional unknowns η z R n R and ι C R n C, is equivalent to the system L (ι L ) d dt ι L = AT L η, (4) coupled with the linear system E d dt x l = Ax l + Bu l, y l = BT x l. (5a) Here, the system matrices are in the MNA form A C CAC T A R GA T R A L A V E = L, A= A T L AV T, B= A I (5b) I with incidence and element matrices A C =[ A C, A R =[ A R A1 R A2 R, A I I L =[ A L AV, A V =[ A C, (5c) AI Ḡ A A I =[ 2 R A L, C = C, L = L, G = G I 1, (5d) G 2 and x T l =[ η T ηz T ι T L ιv T ιt C, u T l = [ ιi T ιt z ι T L u T V ut C, y T l = [ y T 1 yt 2 yt 3 yt 4 yt 5. With (2), equivalency here means that [x T ηz T ι T C T solves (1) and (3) if and only if [x T l ιt L T solves (5) and (4). More details and the proof can be found in [ Model-Order Reduction of the Linear Subsystem We now apply the PABTEC method [4 to the linear system (5) with a transfer function G(s) = B T (se A) 1 B. The assumptions above on the nonlinear system (1) guarantee that the following projected Lur e equations E X(A BB T ) T +(A BB T )XE T + 2P l BB T P T l = 2K c K T c, EXB P l BM T = K cj T c, J c J T c = I M M T, X = P rxp T r (6)
4 4 Andreas Steinbrecher and Tatjana Stykel and E T Y(A BB T )+(A BB T ) T YE + 2Pr T BB T P r = 2Ko T K o, E T YB+Pr T BM = Ko T J o, Jo T J o = I M T M, Y = (7) PT l YP l are solvable for X, K c, J c and Y, K o, J o, respectively, Here P r and P l are the spectral projectors onto the right and left deflating subspaces of λe (A BB T ) corresponding to the finite eigenvalues and M = I 2lim s B T (se A+BB T ) 1 B. The minimal solutions X min and Y min of (6) and (7) that satisfy X min X and Y min Y for all symmetric solutions X and Y of (6) and (7), are called the controllability and observability Gramians of system (5). Using the block structure of the system matrices (5b), we can show that P l = S int P r S int and Y min = S int X min S int with S int = diag(i, I, I). Model reduction consists in approximating the large-scale DAE system (5) of order n l = n η + n R + n L + n V + n C by a reduced-order model Ê d dt ˆx l = Â ˆx l + ˆBu l, ŷ = Ĉ ˆx l, (8) where ˆx l R r and r n l. It is required that the approximate system (8) captures the input-output behavior of (5) to a required accuracy and preserves passivity. The PABTEC model reduction method is based on transforming system (5) into a balanced form whose controllability and observability Gramians are equal and diagonal. Then a reduced-order model is computed by truncation the states corresponding to the small diagonal elements of the balanced Gramians. We summarize the PABTEC method in Algorithm 1. Algorithm 1. Passivity-preserving balanced truncation for electrical circuits (PABTEC). 1. Compute the Cholesky factor R of the minimal solution X min = RR T of (6). 2. Compute the eigenvalue decompositions R T S int ER = [U 1 U 2 diag(λ 1,Λ 2 )[V 1 V 2 T and (I M )S ext = U Λ U T, where [U 1 U 2, [V 1 V 2 and U are orthogonal, Λ 1 = diag(λ 1,...,λ r f ), Λ 2 = diag(λ rf +1,...,λ q) and Λ = diag(ˆλ 1,..., ˆλ m ). 3. Compute the reduced-order system (8) with [ I Ê =, Â= 1 [ 2W T AT 2W T [ BC 2 W 2 B B T, ˆB= T [ B T 2I B C B /, Ĉ T T T B = 2 C / T, 2 where B = S Λ 1/2 U T S ext, C = U Λ 1/2, S ext = diag(i, I), W = S int RU 1 Λ 1 1/2, S = diag(sign(ˆλ 1 ),...,sign(ˆλ m )), Λ = diag( ˆλ 1,..., ˆλ m ), T = RU 1 S 1 Λ 1 1/2, S 1 = diag(sign(λ 1 ),...,sign(λ rf )), Λ 1 = diag( λ 1,..., λ rf ). One can show that the reduced model computed by the PABTEC method preserves passivity and we have the error bound G G H 2 I + G 2 H (λ rf λ q), provided 2 I + G H (λ rf λ q) < 1, see [4 for details. If D c = I M M T is nonsingular, then D o = I M T M is also nonsingular and the projected Lur e equation (6) is equivalent to the projected Riccati equation P l HP T l + FXF T + EXE T + EXQXE T =, X = P r XP T r (9) with F = A BB T 2P l BD 1 o M T BT P r, H = 2BD 1 o B T and Q = 2BD 1 c B T. Such an equation can be solved via Newton s method [4. Note that the matrix M and the
5 Model order reduction of electrical circuits with nonlinear elements 5 projectors P l and P r required in (9) can be constructed in explicit form exploiting the topological structure of the circuit equations (5), see [4. For large-scale problems, the numerical solution of projected Lur e equations is currently under investigation. 3.3 Recoupling of the Reduced Linear Subsystem and the Nonlinear Subsystem Let ˆB = [ ˆB 1 ˆB 2 ˆB 3 ˆB 4 ˆB 5 and Ĉ T = [ Ĉ1 T Ĉ2 T Ĉ3 T Ĉ4 T Ĉ5 T in (8) be partitioned in blocks according to u l and y l, respectively. Since the vector ŷ l = Ĉ ˆx l is an approximation to the output vector y l of system (5), we have (A 2 R )T η η z Ĉ 2 ˆx l, A T L η Ĉ 3 ˆx l, ι C Ĉ 5 ˆx l. (1) Then equations (4) and (3b) are approximated by L (ˆι L ) d dt ˆι L = Ĉ 3 ˆx l, C (û C ) d dt û C = Ĉ 5 ˆx l, (11) respectively, where ˆι L and û C are approximations to ι L and u C, respectively. Furthermore, for ι z defined in (2), η z defined in (3a) and u R = A T R η Rn R, we have ι z = (G 1 )G 1 1 g(u R ) G 2 u R, (A 2 R )T η η z = AT R η + G 1 g(at R η) = u R + G 1 1 g(u R ). (12) Then the first equation in (1) together with (12) imply the relation 1 = G 1 Ĉ 2 ˆx l G 1 û R + g(û R ), (13) where û R approximates u R. Combining (8), (11) and (13) and adding to ˆx l also the approximations ˆι L, û C, and û R as state variables, we get the reduced-order nonlinear DAE system with Ê E ˆ( ˆx) = L (ˆι L ) C (û C ), ˆf(ˆx)=, B= ˆ g(û R ) E ˆ( ˆx) dt d ˆx = Aˆ ˆx+ ˆf(ˆx)+ Bu, ˆ ˆB 1 ˆB 4 A ˆ=, CT r = ŷ = ˆ C ˆx Â+ ˆB 2 (G 1 )Ĉ 2 ˆB 3 ˆB 5 ˆB 2 G 1 Ĉ 3 Ĉ 5 G 1 Ĉ 2 G 1 Ĉ T 1 Ĉ T 4, ˆx l ˆx= ˆι Ḽ u C, û R, ŷ= [ ŷ1 ŷ4 that approximates the original nonlinear system (1). This reduced model can now be used instead of (1) in the analysis of the dynamical behavior of the circuit.
6 6 Andreas Steinbrecher and Tatjana Stykel 4 Numerical Experiments Consider an electrical circuit constructed by 1 repetitions of subcircuits containing one inductor, two capacitors, and two resistors, as shown in the figure on the right. In the 1st, 11st, u V g( ) L ( ) u V (t) = sin 1 (1πt) 21st, etc., subcircuits, a linear resistor is replaced by a diode. Furthermore, in the 1th, 2th, 3th, etc., subcircuits, the linear inductor is replaced by a nonlinear inductor. For more details, we refer to [5. original system: 1 voltage source, 199 linear resistors, 1 diodes, 991 linear inductors, 1 nonlinear inductors, 2 linear capacitors dimension of the original system 43 simulation time for the original 4557s system (t o ) tolerance for model reduction 1e-5 time for model reduction 822s dimension of the reduced system 23 simulation time for the reduced 67s system (t r ) absolute error in the output 4.4e-6 speedup (t o /t r ) 68.5 i V i V x 1 3 Output: negative current through voltage source orig. system red. system t x Error in the output t Fig. 1 (Left) numerical results; (right) simulation results for the original and the reduced systems. The numerical simulation is done for t [,.5 seconds using the BDF method of order 2 with fixed stepsize of length The numerical results for the prescribed tolerance is given in Figures 1. In the upper plot, the output y(t) = ι V (t) of the original system and the output ŷ(t) = ˆι V (t) of the reduced system are presented. In the lower plot, the error ι V (t) = ŷ(t) y(t) is shown. References 1. Freund, R.: SPRIM: structure-preserving reduced-order interconnect macromodeling. In: Technical Digest of the 24 IEEE/ACM International Conference on Computer-Aided Design, pp Los Alamos, CA (24) 2. Heinkenschloss, M., Reis, T.: Model order reduction for a class of nonlinear electrical circuits. In: Proc. of the 48th IEEE Conf. on Decision and Control, pp IEEE (29) 3. Odabasioglu, A., Celik, M., Pileggi, L.: PRIMA: Passive reduced-order interconnect macromodeling algorithm. IEEE Trans. Circuits Syst. 17(8), (1998) 4. Reis, T., Stykel, T.: PABTEC: Passivity-preserving balanced truncation for electrical circuits. IEEE Trans. Computer-Aided Design Integr. Circuits Syst. 29(9), (21) 5. Steinbrecher, A., Stykel, T.: Model order reduction of nonlinear circuits with a small number of nonliniear elements. Technical Report, TU Berlin, Berlin, Germany. In preparation
Model reduction of nonlinear circuit equations
Model reduction of nonlinear circuit equations Tatjana Stykel Technische Universität Berlin Joint work with T. Reis and A. Steinbrecher BIRS Workshop, Banff, Canada, October 25-29, 2010 T. Stykel. Model
More informationModel Order Reduction of Electrical Circuits with Nonlinear Elements
Model Order Reduction of Electrical Circuits with Nonlinear Elements Tatjana Stykel and Technische Universität Berlin July, 21 Model Order Reduction of Electrical Circuits with Nonlinear Elements Contents:
More informationPassive Interconnect Macromodeling Via Balanced Truncation of Linear Systems in Descriptor Form
Passive Interconnect Macromodeling Via Balanced Truncation of Linear Systems in Descriptor Form Boyuan Yan, Sheldon X.-D. Tan,PuLiu and Bruce McGaughy Department of Electrical Engineering, University of
More informationModel reduction of coupled systems
Model reduction of coupled systems Tatjana Stykel Technische Universität Berlin ( joint work with Timo Reis, TU Kaiserslautern ) Model order reduction, coupled problems and optimization Lorentz Center,
More informationThe Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications
MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department
More informationIdentification of Electrical Circuits for Realization of Sparsity Preserving Reduced Order Models
Identification of Electrical Circuits for Realization of Sparsity Preserving Reduced Order Models Christof Kaufmann 25th March 2010 Abstract Nowadays very-large scale integrated circuits contain a large
More informationSecond-Order Balanced Truncation for Passive Order Reduction of RLCK Circuits
IEEE RANSACIONS ON CIRCUIS AND SYSEMS II, VOL XX, NO. XX, MONH X Second-Order Balanced runcation for Passive Order Reduction of RLCK Circuits Boyuan Yan, Student Member, IEEE, Sheldon X.-D. an, Senior
More informationStructure-Preserving Model Order Reduction of RCL Circuit Equations
Structure-Preserving Model Order Reduction of RCL Circuit Equations Roland W. Freund Department of Mathematics, University of California at Davis, One Shields Avenue, Davis, CA 95616, U.S.A. freund@math.ucdavis.edu
More informationEE5900 Spring Lecture 5 IC interconnect model order reduction Zhuo Feng
EE59 Spring Parallel VLSI CAD Algorithms Lecture 5 IC interconnect model order reduction Zhuo Feng 5. Z. Feng MU EE59 In theory we can apply moment matching for any order of approximation But in practice
More informationTowards parametric model order reduction for nonlinear PDE systems in networks
Towards parametric model order reduction for nonlinear PDE systems in networks MoRePas II 2012 Michael Hinze Martin Kunkel Ulrich Matthes Morten Vierling Andreas Steinbrecher Tatjana Stykel Fachbereich
More informationComparison of Model Reduction Methods with Applications to Circuit Simulation
Comparison of Model Reduction Methods with Applications to Circuit Simulation Roxana Ionutiu, Sanda Lefteriu, and Athanasios C. Antoulas Department of Electrical and Computer Engineering, Rice University,
More informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More informationPassivity Preserving Model Reduction for Large-Scale Systems. Peter Benner.
Passivity Preserving Model Reduction for Large-Scale Systems Peter Benner Mathematik in Industrie und Technik Fakultät für Mathematik Sonderforschungsbereich 393 S N benner@mathematik.tu-chemnitz.de SIMULATION
More informationModel reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
More informationParallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation
Parallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation Peter Benner 1, Enrique S. Quintana-Ortí 2, Gregorio Quintana-Ortí 2 1 Fakultät für Mathematik Technische Universität
More informationFactorization of Indefinite Systems Associated with RLC Circuits. Patricio Rosen Wil Schilders, Joost Rommes
Factorization of Indefinite Systems Associated with RLC Circuits Patricio Rosen Wil Schilders, Joost Rommes Outline Circuit Equations Solution Methods (Schilders Factorization Incidence Matrix Decomposition
More informationMODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS
MODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS Ulrike Baur joint work with Peter Benner Mathematics in Industry and Technology Faculty of Mathematics Chemnitz University of Technology
More informationPractical Considerations For Passive Reduction of RLC Circuits
Practical Considerations For Passive Reduction of RLC Circuits Altan Odabasioglu and Mustafa Celik {altan, celik}@montereydesign.com Monterey Design Systems, Inc. Sunnyvale, CA 94089 Lawrence T. Pileggi
More information3 Gramians and Balanced Realizations
3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations
More informationStructured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control
Electronic Journal of Linear Algebra Volume 34 Volume 34 08) Article 39 08 Structured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control Christian Mehl Technische Universitaet
More informationModel Order Reduction using SPICE Simulation Traces. Technical Report
Model Order Reduction using SPICE Simulation Traces Paul Winkler, Henda Aridhi, and Sofiène Tahar Department of Electrical and Computer Engineering, Concordia University, Montreal, Canada pauwink@web.de,
More informationStructure preserving Krylov-subspace methods for Lyapunov equations
Structure preserving Krylov-subspace methods for Lyapunov equations Matthias Bollhöfer, André Eppler Institute Computational Mathematics TU Braunschweig MoRePas Workshop, Münster September 17, 2009 System
More informationStability and Inertia Theorems for Generalized Lyapunov Equations
Published in Linear Algebra and its Applications, 355(1-3, 2002, pp. 297-314. Stability and Inertia Theorems for Generalized Lyapunov Equations Tatjana Stykel Abstract We study generalized Lyapunov equations
More informationLINEAR ALGEBRA QUESTION BANK
LINEAR ALGEBRA QUESTION BANK () ( points total) Circle True or False: TRUE / FALSE: If A is any n n matrix, and I n is the n n identity matrix, then I n A = AI n = A. TRUE / FALSE: If A, B are n n matrices,
More informationADI-preconditioned FGMRES for solving large generalized Lyapunov equations - A case study
-preconditioned for large - A case study Matthias Bollhöfer, André Eppler TU Braunschweig Institute Computational Mathematics Syrene-MOR Workshop, TU Hamburg October 30, 2008 2 / 20 Outline 1 2 Overview
More informationDesign of State Observer for a Class of Non linear Singular Systems Described by Takagi-Sugeno Model
Contemporary Engineering Sciences, Vol. 6, 213, no. 3, 99-19 HIKARI Ltd, www.m-hikari.com Design of State Observer for a Class of Non linear Singular Systems Described by Takagi-Sugeno Model Mohamed Essabre
More informationMay 9, 2014 MATH 408 MIDTERM EXAM OUTLINE. Sample Questions
May 9, 24 MATH 48 MIDTERM EXAM OUTLINE This exam will consist of two parts and each part will have multipart questions. Each of the 6 questions is worth 5 points for a total of points. The two part of
More informationPassive reduced-order modeling via Krylov-subspace methods
Passive reduced-order modeling via Krylov-subspace methods Roland W. Freund Bell Laboratories Room 2C 525 700 Mountain Avenue Murray Hill, New Jersey 07974 0636, U.S.A. freund@research.bell-labs.com Keywords:
More informationModel reduction of interconnected systems
Model reduction of interconnected systems A Vandendorpe and P Van Dooren 1 Introduction Large scale linear systems are often composed of subsystems that interconnect to each other Instead of reducing the
More informationOn Linear-Quadratic Control Theory of Implicit Difference Equations
On Linear-Quadratic Control Theory of Implicit Difference Equations Daniel Bankmann Technische Universität Berlin 10. Elgersburg Workshop 2016 February 9, 2016 D. Bankmann (TU Berlin) Control of IDEs February
More informationVariational Integrators for Electrical Circuits
Variational Integrators for Electrical Circuits Sina Ober-Blöbaum California Institute of Technology Joint work with Jerrold E. Marsden, Houman Owhadi, Molei Tao, and Mulin Cheng Structured Integrators
More informationAn Efficient Graph Sparsification Approach to Scalable Harmonic Balance (HB) Analysis of Strongly Nonlinear RF Circuits
Design Automation Group An Efficient Graph Sparsification Approach to Scalable Harmonic Balance (HB) Analysis of Strongly Nonlinear RF Circuits Authors : Lengfei Han (Speaker) Xueqian Zhao Dr. Zhuo Feng
More informationStructure preserving model reduction of network systems
Structure preserving model reduction of network systems Jacquelien Scherpen Jan C. Willems Center for Systems and Control & Engineering and Technology institute Groningen Reducing dimensions in Big Data:
More informationParallel VLSI CAD Algorithms. Lecture 1 Introduction Zhuo Feng
Parallel VLSI CAD Algorithms Lecture 1 Introduction Zhuo Feng 1.1 Prof. Zhuo Feng Office: EERC 513 Phone: 487-3116 Email: zhuofeng@mtu.edu Class Website http://www.ece.mtu.edu/~zhuofeng/ee5900spring2012.html
More informationModel-Order Reduction of High-Speed Interconnects: Challenges and Opportunities
Model-Order Reduction of High-Speed Interconnects: Challenges and Opportunities Michel Nakhla Carleton University Canada Model Reduction for Complex Dynamical Systems Berlin 2010 EMI Delay Crosstalk Reflection
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 informationRobust Multivariable Control
Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms
More informationTransient Sensitivity Analysis CASA Day 13th Nov 2007 Zoran Ilievski. Zoran Ilievski Transient Sensitivity Analysis
CASA Day 13th Nov 2007 Talk Structure Talk Structure Introduction Talk Structure Introduction Recap Sensitivity Talk Structure Introduction Recap Sensitivity Examples and Results Talk Structure Introduction
More informationAPPLICATION TO TRANSIENT ANALYSIS OF ELECTRICAL CIRCUITS
EECE 552 Numerical Circuit Analysis Chapter Nine APPLICATION TO TRANSIENT ANALYSIS OF ELECTRICAL CIRCUITS I. Hajj Application to Electrical Circuits Method 1: Construct state equations = f(x, t) Method
More informationModel Reduction for Unstable Systems
Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model
More informationModel Reduction for Linear Dynamical Systems
Summer School on Numerical Linear Algebra for Dynamical and High-Dimensional Problems Trogir, October 10 15, 2011 Model Reduction for Linear Dynamical Systems Peter Benner Max Planck Institute for Dynamics
More informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More informationPort-Hamiltonian Realizations of Linear Time Invariant Systems
Port-Hamiltonian Realizations of Linear Time Invariant Systems Christopher Beattie and Volker Mehrmann and Hongguo Xu December 16 215 Abstract The question when a general linear time invariant control
More informationRecycling Krylov Subspaces for Solving Linear Systems with Successively Changing Right-Hand Sides Arising in Model Reduction
Recycling Krylov Subspaces for Solving Linear Systems with Successively Changing Right-Hand Sides Arising in Model Reduction Peter Benner and Lihong Feng Abstract We discuss the numerical solution of successive
More informationNumerical Treatment of Unstructured. Differential-Algebraic Equations. with Arbitrary Index
Numerical Treatment of Unstructured Differential-Algebraic Equations with Arbitrary Index Peter Kunkel (Leipzig) SDS2003, Bari-Monopoli, 22. 25.06.2003 Outline Numerical Treatment of Unstructured Differential-Algebraic
More informationTowards One-Step Multirate Methods in Chip Design
Bergische Universität Wuppertal Fachbereich Mathematik und Naturwissenschaften Lehrstuhl für Angewandte Mathematik und Numerische Mathematik Preprint BUW-AMNA 04/09 M. Striebel, M.Günther Towards One-Step
More informationSolving projected generalized Lyapunov equations using SLICOT
Solving projected generalized Lyapunov equations using SLICOT Tatjana Styel Abstract We discuss the numerical solution of projected generalized Lyapunov equations. Such equations arise in many control
More informationPeter C. Müller. Introduction. - and the x 2 - subsystems are called the slow and the fast subsystem, correspondingly, cf. (Dai, 1989).
Peter C. Müller mueller@srm.uni-wuppertal.de Safety Control Engineering University of Wuppertal D-497 Wupperta. Germany Modified Lyapunov Equations for LI Descriptor Systems or linear time-invariant (LI)
More informationLecture notes: Applied linear algebra Part 1. Version 2
Lecture notes: Applied linear algebra Part 1. Version 2 Michael Karow Berlin University of Technology karow@math.tu-berlin.de October 2, 2008 1 Notation, basic notions and facts 1.1 Subspaces, range and
More informationOptimal control of nonstructured nonlinear descriptor systems
Optimal control of nonstructured nonlinear descriptor systems TU Berlin DFG Research Center Institut für Mathematik MATHEON Workshop Elgersburg 19.02.07 joint work with Peter Kunkel Overview Applications
More informationSyreNe System Reduction for Nanoscale IC Design
System Reduction for Nanoscale Max Planck Institute for Dynamics of Complex Technical Systeme Computational Methods in Systems and Control Theory Group Magdeburg Technische Universität Chemnitz Fakultät
More informationAdvanced Computational Methods for VLSI Systems. Lecture 4 RF Circuit Simulation Methods. Zhuo Feng
Advanced Computational Methods for VLSI Systems Lecture 4 RF Circuit Simulation Methods Zhuo Feng 6. Z. Feng MTU EE59 Neither ac analysis nor pole / zero analysis allow nonlinearities Harmonic balance
More informationModel reduction of large-scale dynamical systems
Model reduction of large-scale dynamical systems Lecture III: Krylov approximation and rational interpolation Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/
More informationIterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations
Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Klajdi Sinani Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University
More informationEE 16B Midterm 2, March 21, Name: SID #: Discussion Section and TA: Lab Section and TA: Name of left neighbor: Name of right neighbor:
EE 16B Midterm 2, March 21, 2017 Name: SID #: Discussion Section and TA: Lab Section and TA: Name of left neighbor: Name of right neighbor: Important Instructions: Show your work. An answer without explanation
More informationLINEAR-QUADRATIC OPTIMAL CONTROL OF DIFFERENTIAL-ALGEBRAIC SYSTEMS: THE INFINITE TIME HORIZON PROBLEM WITH ZERO TERMINAL STATE
LINEAR-QUADRATIC OPTIMAL CONTROL OF DIFFERENTIAL-ALGEBRAIC SYSTEMS: THE INFINITE TIME HORIZON PROBLEM WITH ZERO TERMINAL STATE TIMO REIS AND MATTHIAS VOIGT, Abstract. In this work we revisit the linear-quadratic
More informationKrylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations
Krylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations Peter Benner and Tobias Breiten Abstract We discuss Krylov-subspace based model reduction
More informationWHILE designing VLSI chips, engineers need to take
IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS TurboMOR-RC: an Efficient Model Order Reduction Technique for RC Networks with Many Ports Denis Oyaro, Student Member, IEEE, and Piero
More informationElectrical Circuits I
Electrical Circuits I This lecture discusses the mathematical modeling of simple electrical linear circuits. When modeling a circuit, one ends up with a set of implicitly formulated algebraic and differential
More informationBalancing of Lossless and Passive Systems
Balancing of Lossless and Passive Systems Arjan van der Schaft Abstract Different balancing techniques are applied to lossless nonlinear systems, with open-loop balancing applied to their scattering representation.
More informationDefinite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More informationA geometric Birkhoffian formalism for nonlinear RLC networks
Journal of Geometry and Physics 56 (2006) 2545 2572 www.elsevier.com/locate/jgp A geometric Birkhoffian formalism for nonlinear RLC networks Delia Ionescu Institute of Mathematics, Romanian Academy of
More informationVariational Methods for Solving Warped Multirate PDAEs
Bergische Universität Wuppertal Fachbereich Mathematik und Naturwissenschaften Lehrstuhl für Angewandte Mathematik und Numerische Mathematik Preprint BUW-AMNA 5/ Roland Pulch Variational Methods for Solving
More informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More informationSTA141C: Big Data & High Performance Statistical Computing
STA141C: Big Data & High Performance Statistical Computing Numerical Linear Algebra Background Cho-Jui Hsieh UC Davis May 15, 2018 Linear Algebra Background Vectors A vector has a direction and a magnitude
More information18.06 Professor Johnson Quiz 1 October 3, 2007
18.6 Professor Johnson Quiz 1 October 3, 7 SOLUTIONS 1 3 pts.) A given circuit network directed graph) which has an m n incidence matrix A rows = edges, columns = nodes) and a conductance matrix C [diagonal
More informationAn Optimum Fitting Algorithm for Generation of Reduced-Order Models
An Optimum Fitting Algorithm for Generation of Reduced-Order Models M.M. Gourary 1, S.G. Rusakov 1, S.L. Ulyanov 1, M.M. Zharov 1, B.J. Mulvaney 2 1 IPPM, Russian Academy of Sciences, Moscow 1523, e-mail:
More informationKrylov Techniques for Model Reduction of Second-Order Systems
Krylov Techniques for Model Reduction of Second-Order Systems A Vandendorpe and P Van Dooren February 4, 2004 Abstract The purpose of this paper is to present a Krylov technique for model reduction of
More informationPerturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop
Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael
More informationSparse Implicit Projection (SIP) for Reduction of General Many-Terminal Networks
Sparse mplicit Projection (SP) for Reduction of General Many-Terminal Networks Zuochang Ye, Dmitry Vasilyev, Zhenhai Zhu, Joel R. Phillips Cadence Research Laboratories Massachusetts nstitute of Technology
More informationKrylov subspace methods for projected Lyapunov equations
Krylov subspace methods for projected Lyapunov equations Tatjana Stykel and Valeria Simoncini Technical Report 735-2010 DFG Research Center Matheon Mathematics for key technologies http://www.matheon.de
More informationApproximation of the Linearized Boussinesq Equations
Approximation of the Linearized Boussinesq Equations Alan Lattimer Advisors Jeffrey Borggaard Serkan Gugercin Department of Mathematics Virginia Tech SIAM Talk Series, Virginia Tech, April 22, 2014 Alan
More informationDecomposition of Linear Port-Hamiltonian Systems
American ontrol onference on O'Farrell Street, San Francisco, A, USA June 9 - July, Decomposition of Linear Port-Hamiltonian Systems K. Höffner and M. Guay Abstract It is well known that the power conserving
More informationLecture: Quadratic optimization
Lecture: Quadratic optimization 1. Positive definite och semidefinite matrices 2. LDL T factorization 3. Quadratic optimization without constraints 4. Quadratic optimization with constraints 5. Least-squares
More informationOptimal Sensor and Actuator Location for Descriptor Systems using Generalized Gramians and Balanced Realizations
Optimal Sensor and Actuator Location for Descriptor Systems using Generalized Gramians and Balanced Realizations B. MARX D. KOENIG D. GEORGES Laboratoire d Automatique de Grenoble (UMR CNRS-INPG-UJF B.P.
More informationIBIS Modeling Using Latency Insertion Method (LIM)
IBIS Modeling Using Latency Insertion Method (LIM) José E. Schutt Ainé University of Illinois at Urbana- Champaign Jilin Tan, Ping Liu, Feras Al Hawari, Ambrish arma Cadence Design Systems European IBIS
More informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Serkan Gugercin, Tatjana Stykel, Sarah Wyatt Model Reduction of Descriptor Systems by Interpolatory Projection Methods Preprint Nr. 3/213 29.
More informationAutomatic Generation of Geometrically Parameterized Reduced Order Models for Integrated Spiral RF-Inductors
Automatic Generation of Geometrically Parameterized Reduced Order Models for Integrated Spiral RF-Inductors L. Daniel and J. White Department of Electrical Engineering and Computer Science Research Lab
More informationAutomatic Generation of Geometrically Parameterized Reduced Order Models for Integrated Spiral RF-Inductors
Automatic Generation of Geometrically Parameterized Reduced Order Models for Integrated Spiral RF-Inductors Luca Daniel Massachusetts Institute of Technology Dept. of Elect. Engin. and Comp. Science Research
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Balanced Truncation Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu These slides are based on the recommended textbook: A.C. Antoulas, Approximation of
More informationWeighted balanced realization and model reduction for nonlinear systems
Weighted balanced realization and model reduction for nonlinear systems Daisuke Tsubakino and Kenji Fujimoto Abstract In this paper a weighted balanced realization and model reduction for nonlinear systems
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION
ISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION Harald K. Wimmer 1 The set of all negative-semidefinite solutions of the CARE A X + XA + XBB X C C = 0 is homeomorphic
More informationPoor Man s TBR: A Simple Model Reduction Scheme
Poor Man s TBR: A Simple Model Reduction Scheme Joel Phillips Cadence Berkeley Laboratories Cadence Design Systems San Jose, CA 95134, U.S.A. jrp@cadence.com L. Miguel Silveira INESC ID/IST, Cadence Laboratories
More informationPassivity Preserving Model Reduction and Selection of Spectral Zeros M A R Y A M S A A D V A N D I
Passivity Preserving Model Reduction and Selection of Spectral Zeros M A R Y A M S A A D V A N D I Master of Science Thesis Stockholm, Sweden 28 Passivity Preserving Model Reduction and Selection of Spectral
More informationSynopsis of Numerical Linear Algebra
Synopsis of Numerical Linear Algebra Eric de Sturler Department of Mathematics, Virginia Tech sturler@vt.edu http://www.math.vt.edu/people/sturler Iterative Methods for Linear Systems: Basics to Research
More informationModel Order Reduction for Systems with Non-Rational Transfer Function Arising in Computational Electromagnetics
Model Order Reduction for Systems with Non-Rational Transfer Function Arising in Computational Electromagnetics Lihong Feng and Peter Benner Abstract We consider model order reduction of a system described
More informationMath 265 Linear Algebra Sample Spring 2002., rref (A) =
Math 265 Linear Algebra Sample Spring 22. It is given that A = rref (A T )= 2 3 5 3 2 6, rref (A) = 2 3 and (a) Find the rank of A. (b) Find the nullityof A. (c) Find a basis for the column space of A.
More informationECE 275A Homework #3 Solutions
ECE 75A Homework #3 Solutions. Proof of (a). Obviously Ax = 0 y, Ax = 0 for all y. To show sufficiency, note that if y, Ax = 0 for all y, then it must certainly be true for the particular value of y =
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationMATH 5640: Functions of Diagonalizable Matrices
MATH 5640: Functions of Diagonalizable Matrices Hung Phan, UMass Lowell November 27, 208 Spectral theorem for diagonalizable matrices Definition Let V = X Y Every v V is uniquely decomposed as u = x +
More informationTHE EVER-INCREASING integration density of modern
IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS, VOL. 3, NO. 7, JULY 22 3 A Practical Regulariation Technique for Modified Nodal Analysis in Large-Scale Time-Domain Circuit
More informationModel Order Reduction of Fully Parameterized Systems by Recursive Least Square Optimization
Model Order Reduction of Fully Parameterized Systems by Recursive Least Square Optimization Zheng Zhang, Ibrahim (Abe) M. Elfadel, and Luca Daniel Research Lab of Electronics, Massachusetts Institute of
More informationKrylov-subspace methods for reduced-order modeling in circuit simulation
Krylov-subspace methods for reduced-order modeling in circuit simulation Roland Freund To cite this version: Roland Freund. Krylov-subspace methods for reduced-order modeling in circuit simulation. Journal
More informationA comparison of model reduction techniques from structural dynamics, numerical mathematics and systems and control
A comparison of model reduction techniques from structural dynamics, numerical mathematics and systems and control B. Besselink a, A. Lutowska b, U. Tabak c, N. van de Wouw a, H. Nijmeijer a, M.E. Hochstenbach
More informationarxiv:math/ v1 [math.ds] 5 Sep 2006
arxiv:math/0609153v1 math.ds 5 Sep 2006 A geometric Birkhoffian formalism for nonlinear RLC networks Delia Ionescu, Institute of Mathematics of the Romanian Academy P.O. Box 1-764, RO-014700, Bucharest,
More informationModel reduction via tangential interpolation
Model reduction via tangential interpolation K. Gallivan, A. Vandendorpe and P. Van Dooren May 14, 2002 1 Introduction Although most of the theory presented in this paper holds for both continuous-time
More informationWideband Modeling of RF/Analog Circuits via Hierarchical Multi-Point Model Order Reduction
2E-1 Wideband Modeling of RF/Analog Circuits via Hierarchical Multi-Point Model Order Reduction Zhenyu Qi, Sheldon X.-D. Tan Department of Electrical Engineering, University of California, Riverside, CA
More informationMATHEMATICAL ENGINEERING TECHNICAL REPORTS. Index Characterization of Differential-Algebraic Equations in Hybrid Analysis for Circuit Simulation
MATHEMATICAL ENGINEERING TECHNICAL REPORTS Index Characterization of Differential-Algebraic Equations in Hybrid Analysis for Circuit Simulation Mizuyo TAKAMATSU and Satoru IWATA (Communicated by Kazuo
More informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Nguyen Thanh Son, Tatjana Stykel Model Order Reduction of Parameterized Circuit Equations Based on Interpolation Preprint Nr. 01/2014 24. Januar
More information