Model reduction of nonlinear circuit equations
|
|
- Lee Dawson
- 5 years ago
- Views:
Transcription
1 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
2 T. Stykel. Model reduction of nonlinear circuit equations p.1 Outline Differential-algebraic equations in circuit simulation Model order reduction problem PAssivity-preserving Balanced Truncation method for Electrical Circuits (PABTEC) Decoupling of linear and nonlinear parts Model reduction of linear equations Recoupling Numerical examples Conclusion
3 Modified Nodal Analysis (MNA) T. Stykel. Model reduction of nonlinear circuit equations p.2 R 6 f 1 C 4 L d 2 e 3 A = [A R, A C, A L, A V, A I ] V 1 R 2 C 3 a b c 0 0 Kirchhoff s current law: Aj = 0, j = [ j T R, jt C, j T L, j T V, jt I Kirchhoff s voltage law: A T η = v, v = [ v T R, vt C, v T L, v T V, vt I Branch constitutive relations: resistors: j R = g(v R ), G(v R ) = g(v R ) v R capacitors: j C = dq(v C), C(v C ) = q(v C) dt v C inductors: v L = dφ(j L), L(j L ) = φ(j L) dt j L ] T ] T
4 MNA circuit equations T. Stykel. Model reduction of nonlinear circuit equations p.3 Consider a linear DAE system E(x) ẋ = A x + f(x) + B u y = B T x with A C(AT C C η)at C A L A V E(x)= 0 L(j L ) 0, A= A T A 0 I L 0 0, B= 0 0, A T V I A R g(a T R η) [ f(x) = 0, u = 0 j I v V ] η, x = j L j V, y = [ ] v I, j V η node potentials, j L, j V, j I currents through inductors, voltage and current sources, v V, v I voltages at voltage and current sources, A R,A C,A L,A V,A I incidence matrices of resistors, capacitors, inductors, voltage and current sources
5 T. Stykel. Model reduction of nonlinear circuit equations p.4 Index Assumptions A V has full column rank (= no V-loops) [A C, A L, A R, A V ] has full row rank (= no I-cutsets) C, L, G are symmetric, positive definite Index characterization [ Estévez Schwarz/Tischendorf 00 ] Index = 0 no voltage sources and every node has a capacitive path to a reference node Index = 1 no CV-loops except for C-loops and no LI-cutsets Index = 2, otherwise
6 T. Stykel. Model reduction of nonlinear circuit equations p.5 Model reduction problem Given a large-scale system E(x) ẋ = A x + f(x) + B u y = C x with x R n and u, y R m, find a reduced-order system Ẽ( x) x = Ã x + f( x) + B u ỹ = C x with x R r, u, ỹ R m, r n. preservation of passivity and stability small approximation error ỹ y tol u for all u U need for computable error bounds numerically stable and efficient methods
7 T. Stykel. Model reduction of nonlinear circuit equations p.6 Model reduction techniques Linear circuit equations Krylov subspace methods ( moment matching ) SyPVL for RC, RL, LC circuits [ Freund et al. 96, 97 ] PRIMA, SPRIM for RLC circuits [ Odabasioglu et al. 96, 97; Freund 04, 05 ] Positive real interpolation [ Antoulas 05, Sorensen 05, Ionutiu et al. 08 ] Balancing-related model reduction methods LyaPABTEC for RC, RL circuits [ Reis/S. 10 ] PABTEC for RLC circuits [ Reis/S. 09 ] Nonlinear circuit equations Proper orthogonal decomposition (POD) [ Verhoeven 08 ] Trajectory piece-wise linear approach (TPWL) [ Rewieński 03 ] (Quadratic) bilinearization + balanced truncation [ Benner/Breiten 10 ]
8 T. Stykel. Model reduction of nonlinear circuit equations p.7 PABTEC Tool [Er,Ar,Br,Cr,... ] = PABTEC(Incidence matrices, element matrices,... ) nonlinear linear Decoupling of inear subcircuits [Erl,Arl,Brl,Crl,... ] = PABTECL(Incidence matrices, element matrices,... ) (no dynamics) no L or no C no L and no C Topology no CVI loops no LVI cutsets no R otherwise (unreducible) Lyapunov Riccati Lur e Preprocessing (Projectors) Preprocessing (Projectors) Preprocessing (Projectors) Solving Lyapunov equations (ADI, Krylov methods) Model reduction Solving Riccati equations (Newton s method) Model reduction Solving Lur e equations (Newton s method) Model reduction Postprocessing Postprocessing Postprocessing linear nonlinear Recoupling of subsystems
9 Decoupling T. Stykel. Model reduction of nonlinear circuit equations p.8 + u in u out Large linear RLC circuits arise in modelling transmission lines and pin packages modelling circuits elements by Maxwell s equations via partial element equivalent circuits (PEEC) Assume that C(A T C η) = A C = [ A Cl, A Cn ], AL = [ A Ll, A Ln ], AR = [ A R l, A R n ], [ C l 0 0 C n (A T C n η) ], L(j L ) = [ Ll 0 0 L n (j Ln ) ] [ ] Gl, g(a T R η) = A T R l η. g n (A T R n η)
10 Replacement of nonlinear elements T. Stykel. Model reduction of nonlinear circuit equations p.9 :-( LI-cutsets may arise :-) :-( LI- or I-cutsets may arise
11 Replacement of nonlinear elements T. Stykel. Model reduction of nonlinear circuit equations p.10 :- additional variables :-) :-( CV- or V-loops may arise
12 Replacement of nonlinear elements T. Stykel. Model reduction of nonlinear circuit equations p.11 :- additional variables :-) :- additional variables
13 Decoupled system T. Stykel. Model reduction of nonlinear circuit equations p.12 Linear RLC equations: E ẋ l = Ax l + B u l y l = B T x l with A C CA T C 0 0 A R GA T R A L A V E = 0 L 0, A= A T L 0 0 A C = [ A Cl ] [, A L = A Ll 0 [ ] G A V A Cn l 0 0 A V =, G = 0 G 1 0, G 2 A T V 0 0 ] [ ] [ A R, A R = l A R n,1 A R n,2, A I = 0 I I A I 0, B= I A I A R n,2 A Ln 0 I 0 [ ] x T l = η T ηz T j T Ll jv T jt, Cn [ ] u T l = j T I jz T j T Ln u T V ut ; Cn Nonlinear equations: C n (v Cn ) d dt u C n = j Cn, L n (j Ln ) d dt j L n = A T L n η, j z = (G 1 + G 2 )G 1 1 g n(a T R n η) G 2 A T R n η, ],
14 T. Stykel. Model reduction of nonlinear circuit equations p.13 Balanced truncation System G = ( E, A, B, C ) is balanced if the controllability and observability Gramians X and Y satisfy X = Y = diag(σ 1,..., σ n ). Idea: balance the system, i.e., find an equivalence transformation ( Ê, Â, ˆB, Ĉ ) = (Wb ET b, W b AT b, W b B, CT b ) ( [ ] [ ] [ ] E = 11 E 12 A, 11 A 12 B, 1, [C 1, C 2 ] E 21 E 22 A 21 A 22 B 2 such that ˆX = Ŷ = diag(σ 1,..., σ n ) and truncate the states corresponding to small σ j G = (E 11, A 11, B 1, C 1 ). ) DAEs: G(s) = C(sE A) 1 B = G sp (s) + P(s) G(s) = G sp (s) + P(s)
15 T. Stykel. Model reduction of nonlinear circuit equations p.14 Projected Lur e equations If G=(E, A, B, C ) is passive, then there exist matrices X = X T 0, J c, K c and Y = Y T 0, J o, K o that satisfy the projected Lur e equations (A BC) XE T + EX(A BC) T + 2P l BB T P T l = 2K c K T c, X = P r XP T r, E XC T P l BM T 0 = K cj T c, I M 0 M T 0 = J cj T c, (A BC) T Y E + E T Y (A BC) + 2Pr T C T CP r = 2Ko T K o, B T Y E + M0 TCP r = Jo T K o, I M0 TM 0 = JT o J o, Y = P T l Y P l, where M 0 = I 2 lim s C(sE A + BC) 1 B, P r and P l are the spectral projectors onto the left and right deflating subspaces of the pencil λe A + BC corresponding to the finite eigenvalues. 0 X min X X max, 0 Y min Y Y max X min controllability Gramian, Y min observability Gramian
16 Passivity-preserving BT method T. Stykel. Model reduction of nonlinear circuit equations p.15 Given a passive system G = (E, A, B, C ). 1. Compute P r, P l, M Compute X min = RR T, Y min = LL T (= solve the Lur e equations). ] 3. Compute the SVD L T Π1 ER = [U 1, U 2 ][ [V 1, V 2 ] T. Π 2 4. Compute the reduced-order model [ ] [ I 0 2 W T AT 2 W T B C Ẽ =, Ã = B CT 2I B C [ ] W B T B = B / [, C = CT C / 2 ] 2 ], with I M 0 = C B, W = LU 1 Π 1/2 1 and T = RV 1 Π 1/2 1.
17 Properties T. Stykel. Model reduction of nonlinear circuit equations p.16 G = (Ẽ, Ã, B, C ) is passive G = (Ẽ, Ã, B, C ) is reciprocal ( G(s) = Σ GT (s)σ ) Error bounds: G H := sup G(iω) 2 ω R If 2 I + G H (π lf π nf ) < 1, then G G H 2 I + G 2 H (π lf π nf ). If 2 I + G H (π lf π nf ) < 1, then G G H 2 I + G 2 H (π lf π nf ).
18 Application to circuit equations T. Stykel. Model reduction of nonlinear circuit equations p.17 A CAT C C 0 0 A GA T R R AIAT I A L A V E = 0 L 0, A BC = A T A I 0 L 0 0, B= 0 0 =C T A T V 0 I 0 I Compute P r and P l using the canonical projectors technique [März 96] P r = H 5 (H 4 H 2 I) H 5 H 4 A L H H 6 0 A T V (H 4H 2 I) A T V H 4A L H 6 0 with H 1 = Z T CRIV A L L 1 A T L Z CRIV, H 2 =..., H 3 = Z T C H 2Z C, H 4 = Z C H 1 3 ZT C, H 5 = Z CRIV H 1 1 ZT CRIV A L L 1 A T L I, H 6 = I L 1 A T L Z CRIV H 1 1 ZT CRIV A L, Z C and Z CRIV are basis matrices for ker A T C and ker[a C, A R, A I, A V ] T. P l = S P T r S T with S = diag(i nη, I nl, I nv )
19 Application to circuit equations T. Stykel. Model reduction of nonlinear circuit equations p.18 A CAT C C 0 0 A GA T R R AIAT I A L A V E = 0 L 0, A BC = A T A I 0 L 0 0, B= 0 0 =C T A T V 0 I 0 I Compute P r and P l = S P T r S T with S = diag(i nr, I nl, I nv ). Compute M 0 = I 2 lim s C(sE A + BC) 1 B M 0 = [ I 2A T I ZH 1 0 ZT A I 2A T I ZH 1 0 ZT A V 2A T V ZH 1 0 ZT A I I + 2A T V ZH 1 0 ZT A V ] where H 0 = Z T (A R GA T R +A I AT I +A V AT V )Z, Z = Z CZ RIV C, Z RIV C is a basis matrix for ker [A R, A I, A V ] T Z C and [Z RIV C, Z RIV C ] is nonsingular.
20 Application to circuit equations T. Stykel. Model reduction of nonlinear circuit equations p.19 Compute M 0, P r, P l = S P T r S T with S = diag(i nr, I nl, I nv ). Compute X min =RR T, Y min =FF T (solve the projected Lur e equations) If D 0 = I M 0 M0 T is nonsingular, then the projected Lur e equations are equivalent to the projected Riccati equation (A BC)XE T + EX(A BC) T + 2P l BB T P T l + 2(EXC T P l BM T 0 )D 1 0 (EXCT P l BM T 0 ) = 0, X = P rxp T r compute a low-rank approximation X min R R T, R R n,k, k n, using the generalized low-rank Newton method [ Benner/St. 10 ] Y min = S X min S T S R R T S T = F F T D 0 is nonsingular, if the circuit contains neither CVI-loops except for C-loops nor LIV-cutsets except for L-cutsets
21 Application to circuit equations T. Stykel. Model reduction of nonlinear circuit equations p.20 Compute M 0, P r, P l = S P T r S T with S = diag(i nr, I nl, I nv ). Compute X min R R T, Y min = S X min S T S R R T S T = F F T. Compute the SVD of F T E R F T E R = R T SE R is symmetric compute the EVD R T SE R Λ1 = [U 1, U 2 ][ instead of the SVD Λ 2 ] [U 1, U 2 ] T W = S RU 1 Λ 1 1/2 and T = RU 1 Λ 1 1/2 sign(λ 1 ) with Λ 1 = diag( λ 1,..., λ lf ), sign(λ 1 ) = diag(sign(λ 1 ),..., sign(λ lf ))
22 T. Stykel. Model reduction of nonlinear circuit equations p.21 Application to circuit equations Compute M 0, P r, P l = S P T r S T with S = diag(i nr, I nl, I nv ). Compute X min R R T. Compute the EVD R T SE R = [U 1, U 2 ][ Λ1 Λ 2 ] [U 1, U 2 ] T and W = S RU 1 Λ 1 1/2, T = RU 1 Λ 1 1/2 sign(λ 1 ). Compute B and C such that C B = I M 0 (I M 0 )Σ with Σ = diag(i ni, I nv ) is symmetric compute the EVD (I M 0 )Σ = U 0 Λ 0 U T 0 B = sign(λ 0 ) Λ 0 1/2 U T 0 Σ and C = U 0 Λ 0 1/2
23 PABTECL algorithm T. Stykel. Model reduction of nonlinear circuit equations p.22 Compute M 0, P r, P l = S P T r S T with S = diag(i nr, I nl, I nv ). Compute X min R R T. Compute the EVD R T SE R = [U 1, U 2 ][ Λ1 Λ 2 ] [U 1, U 2 ] T and W = S RU 1 Λ 1 1/2, T = RU 1 Λ 1 1/2 sign(λ 1 ). Compute the EVD (I M 0 )Σ = U 0 Λ 0 U T 0 with Σ = diag(i ni, I nv ) and B = sign(λ 0 ) Λ 0 1/2 U T 0 Σ, C = U 0 Λ 0 1/2. Compute the reduced-order model [ ] I 0 Ẽ =, Ã = 1 [ 2 W T AT 2 W T BC B CT 2I B C [ W B = T B B / [ ], C = CT, C 2 / ] 2. ],
24 Recoupling T. Stykel. Model reduction of nonlinear circuit equations p.23 Linear reduced-order model: Ẽ x [ l = Ã x l + B1, B 2, B 3, B 4, B ] 5 u l, C 1 A T I η C 2 A T R n η + G 1 1 g n(a T R n η) ỹ l = C 3 x l y l = A T L n η C 4 j V C 5 j Cn Nonlinear equations: C n (v Cn ) d dt u C n = j Cn, L n (j Ln ) d dt j L n = A T L n η, j z = (G 1 + G 2 )G 1 1 g n(a T R n η) G 2 A T R n η
25 Reduced-order nonlinear system T. Stykel. Model reduction of nonlinear circuit equations p.24 We obtain the reduced-order system Ẽ( x) x = Ã x + f( x) + B u ỹ = C x with Ẽ Ã + B 2 (G 1 + G 2 ) C 2 B3 B5 B2 G 1 0 L Ẽ( x) = n ( j Ln ) C n (ũ Cn ) 0, Ã= C C , G 1 C2 0 0 G 1 B 1 B4 0 x l 0 0 B= 0 0, 0 f( x) = 0, x = j Ln ũ Cn. 0 0 g n (ũ Cn ) Remark: If n Cn = 0 and n Ln = 0, then passivity is preserved and under some additional topological conditions we have the error bound ỹ y 2 c(π lf π nf )( u 2 + ỹ 2 ). [ Heinkenschloss/Reis 09 ] ũ R n
26 Example: linear RLC circuit T. Stykel. Model reduction of nonlinear circuit equations p.25 n = , m = resistors inductors capacitors Bounded real characteristic values X min R R T, R R n,84 Reduced model: r = j Frequency responses Full order PABTEC Absolute error and error bound Error Error bound Magnitude 10 5 Magnitude Frequency w Frequency w
27 Example: nonlinear circuit T. Stykel. Model reduction of nonlinear circuit equations p.26 1 Input: voltage source 2000 linear capacitors 1990 linear resistors 991 linear inductors 10 nonlinear inductors 10 diodes 1 voltage source original reduced system system Dimension Simulation time u i V x 10 3 t Output: negative current of the voltage source x 10 5 t Error in the output orig. system red. system 4 Model reduction time 822 Error in the output 4.4e-05 i V t
28 T. Stykel. Model reduction of nonlinear circuit equations p.27 Conclusions and future work Model reduction of nonlinear circuit equations topology based partitioning balancing-related model reduction of linear subsystems with preservation of passivity and computable error bounds Exploiting the structure of MNA matrices E, A, B, C use graph algorithms for computing the basis matrices use modern numerical linear algebra algorithms for solving large-scale projected Riccati/Lyapunov equations MATLAB Toolbox PABTEC Preservation of passivity and error bounds for general circuits Numerical solution of large-scale Lur e equations
Model order reduction of electrical circuits with nonlinear elements
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
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationQUESTION BANK SUBJECT: NETWORK ANALYSIS (10ES34)
QUESTION BANK SUBJECT: NETWORK ANALYSIS (10ES34) NOTE: FOR NUMERICAL PROBLEMS FOR ALL UNITS EXCEPT UNIT 5 REFER THE E-BOOK ENGINEERING CIRCUIT ANALYSIS, 7 th EDITION HAYT AND KIMMERLY. PAGE NUMBERS OF
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 informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
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 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 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 informationMoving Frontiers in Model Reduction Using Numerical Linear Algebra
Using Numerical Linear Algebra Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory Group Magdeburg, Germany Technische Universität Chemnitz
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 informationNetwork Graphs and Tellegen s Theorem
Networ Graphs and Tellegen s Theorem The concepts of a graph Cut sets and Kirchhoff s current laws Loops and Kirchhoff s voltage laws Tellegen s Theorem The concepts of a graph The analysis of a complex
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 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 information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationFluid flow dynamical model approximation and control
Fluid flow dynamical model approximation and control... a case-study on an open cavity flow C. Poussot-Vassal & D. Sipp Journée conjointe GT Contrôle de Décollement & GT MOSAR Frequency response of an
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 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 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 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 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 informationNetwork Topology-2 & Dual and Duality Choice of independent branch currents and voltages: The solution of a network involves solving of all branch currents and voltages. We know that the branch current
More information(a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? Solution: dim N(A) 1, since rank(a) 3. Ax =
. (5 points) (a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? dim N(A), since rank(a) 3. (b) If we also know that Ax = has no solution, what do we know about the rank of A? C(A)
More informationA POD Projection Method for Large-Scale Algebraic Riccati Equations
A POD Projection Method for Large-Scale Algebraic Riccati Equations Boris Kramer Department of Mathematics Virginia Tech Blacksburg, VA 24061-0123 Email: bokr@vt.edu John R. Singler Department of Mathematics
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 informationInductance, Inductors, RL Circuits & RC Circuits, LC, and RLC Circuits
Inductance, Inductors, RL Circuits & RC Circuits, LC, and RLC Circuits Self-inductance A time-varying current in a circuit produces an induced emf opposing the emf that initially set up the timevarying
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 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 informationPassivity Assessment and Model Order Reduction for Linear Time-Invariant Descriptor Systems in VLSI Circuit Simulation.
Abstract of thesis entitled Passivity Assessment and Model Order Reduction for Linear Time-Invariant Descriptor Systems in VLSI Circuit Simulation Submitted by Zheng ZHANG for the degree of Master of Philosophy
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 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 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 informationModel order reduction for multi-terminals systems : with applications to circuit simulation Ionutiu, R.
Model order reduction for multi-terminals systems : with applications to circuit simulation Ionutiu, R. DOI: 10.6100/IR716352 Published: 01/01/2011 Document Version Publisher s PDF, also known as Version
More informationJacobi-Davidson methods and preconditioning with applications in pole-zero analysis
Nat.Lab. Unclassified Report 2002/817 Date of issue: 05/2002 Jacobi-Davidson methods and preconditioning with applications in pole-zero analysis Master s Thesis Joost Rommes Unclassified Report 2002/817
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 informationthe reference terminal. For electrotechnical reasons, the current entering terminal n is given by i n = ; P n; k= i k. The conductance matrix G(v :::
Topological index calculation of DAEs in circuit simulation Caren Tischendorf, Humboldt-University of Berlin Abstract. Electric circuits are present in a number of applications, e.g. in home computers,
More informationAlternating Current Circuits. Home Work Solutions
Chapter 21 Alternating Current Circuits. Home Work s 21.1 Problem 21.11 What is the time constant of the circuit in Figure (21.19). 10 Ω 10 Ω 5.0 Ω 2.0µF 2.0µF 2.0µF 3.0µF Figure 21.19: Given: The circuit
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 informationOutline. Week 5: Circuits. Course Notes: 3.5. Goals: Use linear algebra to determine voltage drops and branch currents.
Outline Week 5: Circuits Course Notes: 3.5 Goals: Use linear algebra to determine voltage drops and branch currents. Components in Resistor Networks voltage source current source resistor Components in
More informationEmpirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems
Empirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems Antoni Ras Departament de Matemàtica Aplicada 4 Universitat Politècnica de Catalunya Lecture goals To review the basic
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 informationA Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations
A Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations Peter Benner Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory
More informationBasics of Network Theory (Part-I)
Basics of Network Theory (PartI). A square waveform as shown in figure is applied across mh ideal inductor. The current through the inductor is a. wave of peak amplitude. V 0 0.5 t (m sec) [Gate 987: Marks]
More informationTaking the Laplace transform of the both sides and assuming that all initial conditions are zero,
The transfer function Let s begin with a general nth-order, linear, time-invariant differential equation, d n a n dt nc(t)... a d dt c(t) a 0c(t) d m = b m dt mr(t)... a d dt r(t) b 0r(t) () where c(t)
More informationStability and Passivity of the Super Node Algorithm for EM Modeling of IC s
Stability and Passivity of the Super Node Algorithm for EM Modeling of IC s M.V. Ugryumova and W.H.A. Schilders Abstract The super node algorithm performs model order reduction based on physical principles.
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 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 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 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 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 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 informationDownloaded 05/27/14 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 35, No. 5, pp. B1010 B1033 c 2013 Society for Industrial and Applied Mathematics MODEL REDUCTION OF DESCRIPTOR SYSTEMS BY INTERPOLATORY PROJECTION METHODS SERKAN GUGERCIN, TATJANA
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationModel Order Reduction via Matlab Parallel Computing Toolbox. Istanbul Technical University
Model Order Reduction via Matlab Parallel Computing Toolbox E. Fatih Yetkin & Hasan Dağ Istanbul Technical University Computational Science & Engineering Department September 21, 2009 E. Fatih Yetkin (Istanbul
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik On Combining Deflation and Iteration to Low-Rank Approximative Solution of Lur e Equations Federico Poloni and Timo Reis Nr. 2012-01 January 2012 ON COMBINING
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 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 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 informationSample ECE275A Midterm Exam Questions
Sample ECE275A Midterm Exam Questions The questions given below are actual problems taken from exams given in in the past few years. Solutions to these problems will NOT be provided. These problems and
More informationAN OVERVIEW OF MODEL REDUCTION TECHNIQUES APPLIED TO LARGE-SCALE STRUCTURAL DYNAMICS AND CONTROL MOTIVATING EXAMPLE INVERTED PENDULUM
Controls Lab AN OVERVIEW OF MODEL REDUCTION TECHNIQUES APPLIED TO LARGE-SCALE STRUCTURAL DYNAMICS AND CONTROL Eduardo Gildin (UT ICES and Rice Univ.) with Thanos Antoulas (Rice ECE) Danny Sorensen (Rice
More informationControl Systems Engineering (Chapter 2. Modeling in the Frequency Domain) Prof. Kwang-Chun Ho Tel: Fax:
Control Systems Engineering (Chapter 2. Modeling in the Frequency Domain) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 02-760-4253 Fax:02-760-4435 Overview Review on Laplace transform Learn about transfer
More informationRational interpolation, minimal realization and model reduction
Rational interpolation, minimal realization and model reduction Falk Ebert Tatjana Stykel Technical Report 371-2007 DFG Research Center Matheon Mathematics for key technologies http://www.matheon.de Rational
More informationGramian based model reduction for descriptor systems
Published in Mathematics of Control, Signals, and Systems 16, 297 319, 2004. Gramian based model reduction for descriptor systems Tatjana Stykel Abstract Model reduction is of fundamental importance in
More informationRecent developments for MOR in the electronics industry
Recent developments for MOR in the electronics industry Wil Schilders Reduced Order Models in Computational Science and Engineering Aachen, January 30-31, 2014 Full proposal oc-2013-1-15312 for a new COST
More informationBasic. Theory. ircuit. Charles A. Desoer. Ernest S. Kuh. and. McGraw-Hill Book Company
Basic C m ш ircuit Theory Charles A. Desoer and Ernest S. Kuh Department of Electrical Engineering and Computer Sciences University of California, Berkeley McGraw-Hill Book Company New York St. Louis San
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 informationSinusoidal Response of RLC Circuits
Sinusoidal Response of RLC Circuits Series RL circuit Series RC circuit Series RLC circuit Parallel RL circuit Parallel RC circuit R-L Series Circuit R-L Series Circuit R-L Series Circuit Instantaneous
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 informationGramians of structured systems and an error bound for structure-preserving model reduction
Gramians of structured systems and an error bound for structure-preserving model reduction D.C. Sorensen and A.C. Antoulas e-mail:{sorensen@rice.edu, aca@rice.edu} September 3, 24 Abstract In this paper
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical Solution Of Matrix Equations Arising In Linear-Quadratic Optimal Control Of Differential-Algebraic Equations Timo Reis, Olaf Rendel Nr. 216-36 December
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 informationAdjoint networks and other elements of circuit theory. E416 4.Adjoint networks
djoint networks and other elements of circuit theory One-port reciprocal networks one-port network is reciprocal if: V I I V = Where and are two different tests on the element Example: a linear impedance
More informationBALANCING-RELATED MODEL REDUCTION FOR LARGE-SCALE SYSTEMS
FOR LARGE-SCALE SYSTEMS Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Recent Advances in Model Order Reduction Workshop TU Eindhoven, November 23,
More informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report June 2009
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 09-21 June 2009 Model order reduction for nonlinear IC models with POD by A. Verhoeven, M. Striebel, E.J.W.
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 informationModel order reduction of large-scale dynamical systems with Jacobi-Davidson style eigensolvers
MAX PLANCK INSTITUTE International Conference on Communications, Computing and Control Applications March 3-5, 2011, Hammamet, Tunisia. Model order reduction of large-scale dynamical systems with Jacobi-Davidson
More informationSome of the different forms of a signal, obtained by transformations, are shown in the figure. jwt e z. jwt z e
Transform methods Some of the different forms of a signal, obtained by transformations, are shown in the figure. X(s) X(t) L - L F - F jw s s jw X(jw) X*(t) F - F X*(jw) jwt e z jwt z e X(nT) Z - Z X(z)
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 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 information