VISA: Versatile Impulse Structure Approximation for Time-Domain Linear Macromodeling
|
|
- Horace Weaver
- 5 years ago
- Views:
Transcription
1 VISA: Versatile Impulse Structure Approximation for Time-Domain Linear Macromodeling Chi-Un Lei and Ngai Wong Dept. of Electrical and Electronic Engineering The University of Hong Kong Presented by C.U. Lei ASP-DAC 21
2 Outline Digital IIR filter approximation technique Eigenvalue-free, initial-guess guess-free, linear macromodeling technique Introduction Macromodeling Process Time Domain Macromodeling Algorithm (VISA) Numerical Results Conclusion
3 Introduction Accurate interconnect/via/package models esp. in high frequency simulation to capture EM effects Compact & efficient macromodels for simulation
4 Introduction Import data Data preprocessing Tabulated data of input-output response Electromagnetic data Time / frequency Calculation / measurement Math. model Export netlist Model postprocessing Macromodel requirement: Accurate, causal, stable, passive
5 Existing algorithms Frequency-domain Vector Fitting (VF) [99] Iterative continuous frequency domain approximation (p input ports, q output ports) ( ) ( ) ( ) ( ) P s p m s ( ) M m uv, uv, uv, uv, ( ) = m Q s = m= q m s f s f s Time-domain Macromodeling Techniques Approximation using (truncated) time-sampled system response Easier to capture e.g. 4SID [99], GPOF [99],, TD-VF(z VF(z) ) [3,8]
6 Problems Issues and Alternatives Numerical-sensitive calculation Affected by the initial poles, noisy signal Require eigenvalues calculation Unstable pole flipping Not an accurate approx. High (multiport) computation complexity Other techniques has much higher computation complexity (Expensive SVD computation) Digital IIR filter approximation Macromodeling SISO LS [ISCAS 8] VISA [ASPDAC 1] Implication: finite-length discrete response sequences... like FIR filter responses!!!
7 Outline Introduction Algorithm (VISA) Walsh s s theorem & Complementary signal Formulation Convergence: SM iteration, model order reduction, P-norm approximation Numerical Results Conclusion
8 Least-Squares Approximation Macromodeling of a time-sampled response = IIR Approx. of FIR filter Problem: Approximation Interpolation problem Input / Output description of an allpass filtering operation Designing an allpass operator 1 L F( z) = f + f1z + + flz 1 N Pz ( ) p + p1z + + pn z 1 N Qz ( ) 1+ qz qn z H ( z) = =, N < L Δ ( z) = F( z) H( z) 1 l Δ z = Δ z Δ z * 2norm: ( ) ( ) ( ) 2 2π j z = 1 dz z
9 Find P(z) Suppose Q(z) Known Walsh s Theorem: z =, z = 1/ α, z = 1/ α,, z = 1/ α * * * N N Pz ( ) p+ pz + + p z = = F( zi ), i =,1,, N Qz ( ) 1 qz q z 1 N i 1 i N i 1 N i + 1 i + + N i Pz ( i ) Δ ( zi) = F( zi) = Δ( z) minimized! 2 Qz ( ) ( N + 1) 1 z Q( z ) Δ ( z) = R( z) where R( z) = r + rz Qz ( ) + + r z ( N + 1) 1 z Q( z ) P( z) Rz ( ) = Fz ( ) Qz ( ) Qz ( ) i Finding z i s can be numerically ill. Smarter way making use of the FIR nature of F(z), we can rewrite Δ(z) as 1 ( L 1) 1 L 1
10 Equivalent Filtering Problem Complementary Signal: All-pass filter! f L f L-1.. f 1 f 1 L-1 L p N p N-1.. p 1 p 1 N-1 N z -L F(z -1 ) z -N Q(z -1 ) / Q(z) z -N P(z -1 ) 1/ Q(z) z -L + z -(L-1) R(z -1 ) - r L-1 r L-2.. r 1 r 1 L-2 L-1 For the first L time instances rl-1 rl-2.. r1 r= z -N Q(z -1 ) / Q(z) f L f L-1.. f 1 f
11 Finding Numerator P(z) ( N + 1) 1 z Qz ( ) Pz ( ) Δ ( z) = R( z) = F( z) Qz ( ) Qz ( ) Easily prove that 2 L 1 r 2 i i= = Δ () z P(z) easily found!! For the first time L instances: f L f L-1.. f 1 f z -N Q(z -1 ) / Q(z) 1 L-1 L r L-1 r L-2.. r 1 r - 1 L-2 L-1
12 Complementary Signal Adopted from [YH62] Original Signal Approximating Signal Allpass Filtering Error Signal (Need to minimize!!!)
13 Finding Denominator Q(z) For the first L time instances f L f L-1.. f 1 f 1 L-1 L z -N Q(z -1 ) / Q(z) r L-1 r L-2.. r 1 r 1 L-2 L-1 *Find Q(z) such that energy output of the all-pass filter A(z)=z -N Q(z -1 )/Q(z) is concentrated in time L, Recall Define 1 N Q( z) 1 q1z = qn z () Q ( z) = 1 Q ( z) = 1+ q z + q z + + q z ( k) ( k) 1 ( k) 2 ( k) N 1 2 N 1 ( ) ( ) 1 ( ) ( 1) ( k k k N ) 1 2 N = 1+ z q + q z + + q z = 1 + z Q ( z) 1 ( k ) 1
14 Finding Q(z) by Iterations () Q ( z) = 1 Q ( z) = 1+ q z + q z + + q z ( k) ( k) 1 ( k) 2 ( k) N 1 2 N 1 ( ) ( ) 1 ( ) ( 1) k k k N ( 1 2 N ) = 1+ z q + q z + + q z = 1 + z Q ( z) 1 ( k ) 1 f L f L-1.. f 1 f 1 L-1 L z -N Q(z -1 ) / Q(z) r L-1 r L-2.. r 1 r 1 L-2 L-1 Define A ( k ) N ( k) 1 z Q ( z ) ( z) = ( k 1) Q ( z) N ( k) 1 ( k) z Q ( z ) N ( k) 1 ( k) U ( z) = X( z) = z Q ( z ) X ( z) ( k 1) Q ( z) ( k) ( N 1) ( k) 1 ( k) N ( k) X ( z) z Q1 ( z ) = U ( z) z X ( z)
15 Finding Q(z) by Iterations ( k) ( N 1) ( k) 1 ( k) N ( k) X ( z) z Q1 ( z ) = U ( z) z X ( z) ( k ) ( k ) x () u () ( k) ( k) ( k ) ( k) x (1) x () q N u (1) ( k ) q N 1 ( k) ( k) = ( k ) x ( N 1) x () x () ( k ) q 1 ( k) ( k) ( k) ( k ) x ( L 1) x ( L N) u ( L 1) x ( L N 1) A q = u + b ( k) ( k) ( k) ( k) ( k) ( k) ( k) Automatically uk ( ): = A q b has minimum norm By iterations, Q (k) (z) converges to Q(z) What s amazing: Q (k) (z) is always stable!
16 Finding Q(z) by Iterations ( k) ( N 1) ( k) 1 ( k) N ( k) X ( z) z Q1 ( z ) = U ( z) z X ( z) ( k ) ( k ) x () u () ( k) ( k) ( k ) ( k) x (1) x () q N u (1) ( k ) q N 1 ( k) ( k) = ( k ) x ( N 1) x () x () ( k ) q 1 ( k) ( k) ( k) ( k ) x ( L 1) x ( L N) u ( L 1) x ( L N 1) A q = u + b ( k) ( k) ( k) ( k) ( k) ( k) ( k) Automatically uk ( ): = A q b has minimum norm By iterations, Q (k) (z) converges to Q(z) What s amazing: Q (k) (z) is always stable!
17 Convergence of VISA Calculation of Q(z) ) is a simplification of Steiglitz- McBride (SM) iteration without initial guess P(z) ) is not required for the calculation of Q(z) min = N s k = 1 min N s k = () i 1 ( ) z Q z N k k L Q Q ( i 1 ) ( z ) 1 k ( i 1 ) ( z ) k 2 z 1 ( k ) F z 2 () i 1 L 1 ( ) ( ) z Q z z F z N k k k 2 Ns 1 () i = min Q 2 ( z ) F( z ) P z k = Q ( i 1 ) ( z ) k () i ( ) k k k 2 SM iteration
18 Model order selection Appropriate model order (N) for efficient and accurate simulation Similar pattern between approximant error and Hankel singular value (HSV) of the impulse response HSV distribution and the ratio of first and last HSV H = eig f f f f f 1 2 L 2 f3 L
19 P-norm approximation Modify the minimization criteria to suit different macromodeling requirements P-norm approximation User-defined norm approximation Norm-constrained approximation min L 1 n= ( ) ( z ) ( k 1 ) ( z ) ( ) ( ) ( ) Δ = min B q d k k k k n p L 1 1 ( k = min Q ) ( z ) H( z ) P z n= Q n p ( k ) ( ) n n n Applying convex programming for solving p p
20 MIMO VISA Extension Walsh s s theorem can be applied to MIMO situation Stacking the responses into a column for LS solving ( k) Less computation complexity TD-VF: O((pq+1) (pq+1)n 2 Lpq) Much higher! VISA: O(N 2 Lpq) ( k) B 1,1 d 1,1 ( k) ( k) B1,2 ( k ) d1,2 q = ( k) ( k) B d p, q p, q Single-port response
21 VISA Features No eigenvalue computations No initial guess required Less initial-pole pole-sensitive calculation Guaranteed stable pole calculation Converge to the near-l 2 optimal solution Low (multi-port) computation complexity P-norm approximation for different modeling requirement Quasi-error bound for model order selection
22 Outline Introduction Algorithm (VISA) Numerical Results Three-port counterwise RF circulator Benchmark examples Conclusion
23 Numerical Examples Three-port counterwise RF circulator More accurate : 18% less avgerage R.M.S. error after convergence Faster : >15X faster for convergence and >17X faster to achieve a -4dB accuracy L 2 error 1 (a) L 2 error 1 2 (b) 1 TD-VF MIMO VISA CPU Time CPU Time
24 Mag. (db) Mag. (db) Mag. (db) Mag. (db) Mag. (db) S(1,1) Mag. (db) Normalized Frequency S(1,3) Normalized Frequency S(2,2) Normalized Frequency S(3,1) Normalized Frequency S(3,3) Normalized Frequency Mag. (db) Mag. (db) Mag. (db) -1-2 S(1,2).5 1 Normalized Frequency S(2,1) Normalized Frequency S(2,3).5 1 Normalized Frequency S(3,2) Normalized Frequency MIMO VISA Original Magnitude Magnitude Magnitude Magnitude Magnitude S(1,1) Time [ns] S(1,3) Time [ns] S(2,2) Time [ns] S(3,1) Time [ns] S(3,3) Time [ns] Magnitude Magnitude Magnitude Magnitude.2 S(1,2) Time [ns] S(2,1) Time [ns] S(2,3) Time [ns] S(3,2) Time [ns] MIMO VISA Original
25 Numerical Examples HSV (error bound) computation: 58 sec. CPU time P-norm approximation 3.5% L 2 error, 29.5% L inf error and 24.5% CPU time reduction Four benchmark examples [IG8] 51% less L 2 error, >21 X faster, comparing to TD-VF Magnitude of HSV 1 (a) Value Magnitude (b) HSV MIMO VISA Model Order
26 Conclusion VISA: Linear Macromodeling using time- sampled response Simplified Steiglitz-McBride algorithm No initial-pole assignment No eigenvalues calculation Robust and efficient computation P-norm approximation
27 Thank you! Thank you! Questions are welcome Authors: C. U. Lei * culei@eee.hku.hk Dr. N. Wong nwong@eee.hku.hk Department of Electrical and Electronic Engineering, The University of Hong Kong
28 Back-Up Slides
29 Walsh Theorem Among the set of rational functions H(z)= )=P(z)/Q(z), with prescribed poles α 1, α 2,..., α n that are fixed and located in z < 1, the best approximation in the LS sense to F(z) ) (analytic in z > 1 and continuous in z 1) is the unique function that interpolates to F (z)) in all the points. z =, 1/α 1, 1/α 2,,..., 1/α n,, where * denotes complex conjugate Approximation Interpolation problem Input / Output description of an allpass filtering operation Design of an allpass operator
30 Complementary Signal Adopted from [YH62] Original Signal Approximating Signal Allpass Filtering Error Signal (Need to minimize!!!)
Abstract of thesis entitled. Submitted by. Lei Chi-Un. for the degree of Doctor of Philosophy at the University of Hong Kong in January 2011
Abstract of thesis entitled VLSI Macromodeling and Signal Integrity Analysis via Digital Signal Processing Techniques Submitted by Lei Chi-Un for the degree of Doctor of Philosophy at the University of
More informationEfficient Linear Macromodeling via Discrete-Time Time-Domain Vector Fitting
st International Conference on VLSI Design Efficient Linear Macromodeling via Discrete-Time Time-Domain Vector Fitting Chi-Un Lei and Ngai Wong Department of Electrical and Electronic Engineering The University
More informationEfficient 2D Linear-Phase IIR Filter Design and Application in Image Filtering
Efficient D Linear-Phase IIR Filter Design and Application in Image Filtering Chi-Un Lei, Chung-Man Cheung, and Ngai Wong Abstract We present an efficient and novel procedure to design two-dimensional
More informationEnhanced Steiglitz-McBride Procedure for. Minimax IIR Digital Filters
Enhanced Steiglitz-McBride Procedure for Minimax IIR Digital Filters Wu-Sheng Lu Takao Hinamoto University of Victoria Hiroshima University Victoria, Canada Higashi-Hiroshima, Japan May 30, 2018 1 Outline
More informationUse: Analysis of systems, simple convolution, shorthand for e jw, stability. Motivation easier to write. Or X(z) = Z {x(n)}
1 VI. Z Transform Ch 24 Use: Analysis of systems, simple convolution, shorthand for e jw, stability. A. Definition: X(z) = x(n) z z - transforms Motivation easier to write Or Note if X(z) = Z {x(n)} z
More informationDavid Weenink. First semester 2007
Institute of Phonetic Sciences University of Amsterdam First semester 2007 Digital s What is a digital filter? An algorithm that calculates with sample values Formant /machine H 1 (z) that: Given input
More informationResponses of Digital Filters Chapter Intended Learning Outcomes:
Responses of Digital Filters Chapter Intended Learning Outcomes: (i) Understanding the relationships between impulse response, frequency response, difference equation and transfer function in characterizing
More informationSubspace Identification
Chapter 10 Subspace Identification Given observations of m 1 input signals, and p 1 signals resulting from those when fed into a dynamical system under study, can we estimate the internal dynamics regulating
More informationMagnitude Vector Fitting to interval data
Available online at wwwsciencedirectcom Mathematics and Computers in Simulation 80 (2009) 572 580 Magnitude Vector Fitting to interval data Wouter Hendrickx a,, Dirk Deschrijver b, Luc Knockaert b, Tom
More informationModel Uncertainty and Robust Stability for Multivariable Systems
Model Uncertainty and Robust Stability for Multivariable Systems ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Devron Profile Control Solutions Outline Representing model uncertainty.
More informationECE503: Digital Signal Processing Lecture 6
ECE503: Digital Signal Processing Lecture 6 D. Richard Brown III WPI 20-February-2012 WPI D. Richard Brown III 20-February-2012 1 / 28 Lecture 6 Topics 1. Filter structures overview 2. FIR filter structures
More informationArchitecture-level Thermal Behavioral Models For Quad-Core Microprocessors
Architecture-level Thermal Behavioral Models For Quad-Core Microprocessors Duo Li Dept. of Electrical Engineering University of California Riverside, CA 951 dli@ee.ucr.edu Sheldon X.-D. Tan Dept. of Electrical
More informationLecture 19 IIR Filters
Lecture 19 IIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/10 1 General IIR Difference Equation IIR system: infinite-impulse response system The most general class
More informationCSE 554 Lecture 7: Alignment
CSE 554 Lecture 7: Alignment Fall 2012 CSE554 Alignment Slide 1 Review Fairing (smoothing) Relocating vertices to achieve a smoother appearance Method: centroid averaging Simplification Reducing vertex
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #21 Friday, October 24, 2003 Types of causal FIR (generalized) linear-phase filters: Type I: Symmetric impulse response: with order M an even
More informationTunable IIR Digital Filters
Tunable IIR Digital Filters We have described earlier two st-order and two nd-order IIR digital transfer functions with tunable frequency response characteristics We shall show now that these transfer
More informationWe have shown earlier that the 1st-order lowpass transfer function
Tunable IIR Digital Filters We have described earlier two st-order and two nd-order IIR digital transfer functions with tunable frequency response characteristics We shall show now that these transfer
More informationNEW STEIGLITZ-McBRIDE ADAPTIVE LATTICE NOTCH FILTERS
NEW STEIGLITZ-McBRIDE ADAPTIVE LATTICE NOTCH FILTERS J.E. COUSSEAU, J.P. SCOPPA and P.D. DOÑATE CONICET- Departamento de Ingeniería Eléctrica y Computadoras Universidad Nacional del Sur Av. Alem 253, 8000
More informationTopic 7: Filter types and structures
ELEN E481: Digital Signal Processing Topic 7: Filter tpes and structures 1. Some filter tpes 2. Minimum and maximum phase 3. Filter implementation structures 1 1. Some Filter Tpes We have seen the basics
More informationFREQUENCY DOMAIN IDENTIFICATION TOOLBOX
NASA Technical Memorandum 109039 FREQUENCY DOMAIN IDENTIFICATION TOOLBOX Lucas G. Horta and Jer-Nan Juang NASA Langley Research Center Hampton, Virginia 23681-0001 Chung-Wen Chen North Carolina State University
More informationGeneralized Pencil Of Function Method
Generalized Pencil Of Function Method In recent times, the methodology of approximating a function by a sum of complex exponentials has found applications in many areas of electromagnetics For example
More informationECE 451 Black Box Modeling
Black Box Modeling Jose E. Schutt-Aine Electrical & Computer Engineering University of Illinois jose@emlab.uiuc.edu Simulation for Digital Design Nonlinear Network Nonlinear Network 3 Linear N-Port with
More informationAutomatic Control Systems theory overview (discrete time systems)
Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations
More informationPart 4: IIR Filters Optimization Approach. Tutorial ISCAS 2007
Part 4: IIR Filters Optimization Approach Tutorial ISCAS 2007 Copyright 2007 Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 24, 2007 Frame # 1 Slide # 1 A. Antoniou Part4: IIR Filters
More informationStability preserving post-processing methods applied in the Loewner framework
Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 1 / 20 Stability preserving post-processing methods applied in the Loewner framework Ion Victor
More informationCOMS 4721: Machine Learning for Data Science Lecture 19, 4/6/2017
COMS 4721: Machine Learning for Data Science Lecture 19, 4/6/2017 Prof. John Paisley Department of Electrical Engineering & Data Science Institute Columbia University PRINCIPAL COMPONENT ANALYSIS DIMENSIONALITY
More informationQualification of tabulated scattering parameters
Qualification of tabulated scattering parameters Stefano Grivet Talocia Politecnico di Torino, Italy IdemWorks s.r.l. stefano.grivet@polito.it 4 th IEEE Workshop on Signal Propagation on Interconnects
More informationAnalog LTI system Digital LTI system
Sampling Decimation Seismometer Amplifier AAA filter DAA filter Analog LTI system Digital LTI system Filtering (Digital Systems) input output filter xn [ ] X ~ [ k] Convolution of Sequences hn [ ] yn [
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #24 Tuesday, November 4, 2003 6.8 IIR Filter Design Properties of IIR Filters: IIR filters may be unstable Causal IIR filters with rational system
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationState will have dimension 5. One possible choice is given by y and its derivatives up to y (4)
A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To
More informationIEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 49, NO. 2, FEBRUARY H(z) = N (z) D(z) =
IEEE TRASACTIOS O CIRCUITS AD SYSTEMS II: AALOG AD DIGITAL SIGAL PROCESSIG, VOL. 49, O. 2, FEBRUARY 2002 145 Transactions Briefs Design of a Class of IIR Eigenfilters With Time- Frequency-Domain Constraints
More informationLecture 7 - IIR Filters
Lecture 7 - IIR Filters James Barnes (James.Barnes@colostate.edu) Spring 204 Colorado State University Dept of Electrical and Computer Engineering ECE423 / 2 Outline. IIR Filter Representations Difference
More informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture 4: Inverse z Transforms & z Domain Analysis Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 008 K. E. Barner
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 informationA Constraint Generation Approach to Learning Stable Linear Dynamical Systems
A Constraint Generation Approach to Learning Stable Linear Dynamical Systems Sajid M. Siddiqi Byron Boots Geoffrey J. Gordon Carnegie Mellon University NIPS 2007 poster W22 steam Application: Dynamic Textures
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION CONTENTS VOLUME VII
CONTENTS VOLUME VII Control of Linear Multivariable Systems 1 Katsuhisa Furuta,Tokyo Denki University, School of Science and Engineering, Ishizaka, Hatoyama, Saitama, Japan 1. Linear Multivariable Systems
More informationA Fast Augmented Lagrangian Algorithm for Learning Low-Rank Matrices
A Fast Augmented Lagrangian Algorithm for Learning Low-Rank Matrices Ryota Tomioka 1, Taiji Suzuki 1, Masashi Sugiyama 2, Hisashi Kashima 1 1 The University of Tokyo 2 Tokyo Institute of Technology 2010-06-22
More informationDigital Signal Processing:
Digital Signal Processing: Mathematical and algorithmic manipulation of discretized and quantized or naturally digital signals in order to extract the most relevant and pertinent information that is carried
More informationQuadrature-Mirror Filter Bank
Quadrature-Mirror Filter Bank In many applications, a discrete-time signal x[n] is split into a number of subband signals { v k [ n]} by means of an analysis filter bank The subband signals are then processed
More informationAnalysis of Finite Wordlength Effects
Analysis of Finite Wordlength Effects Ideally, the system parameters along with the signal variables have infinite precision taing any value between and In practice, they can tae only discrete values within
More informationLet H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 )
Review: Poles and Zeros of Fractional Form Let H() = P()/Q() be the system function of a rational form. Let us represent both P() and Q() as polynomials of (not - ) Then Poles: the roots of Q()=0 Zeros:
More informationEEL3135: Homework #4
EEL335: Homework #4 Problem : For each of the systems below, determine whether or not the system is () linear, () time-invariant, and (3) causal: (a) (b) (c) xn [ ] cos( 04πn) (d) xn [ ] xn [ ] xn [ 5]
More informationS. Gugercin and A.C. Antoulas Department of Electrical and Computer Engineering Rice University
Proceedings of the 39" IEEE Conference on Decision and Control Sydney, Australia December, 2000 A Comparative Study of 7 Algorithms for Model Reduct ion' S. Gugercin and A.C. Antoulas Department of Electrical
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 informationChapter 5 THE APPLICATION OF THE Z TRANSFORM. 5.3 Stability
Chapter 5 THE APPLICATION OF THE Z TRANSFORM 5.3 Stability Copyright c 2005- Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org February 13, 2008 Frame # 1 Slide # 1 A. Antoniou Digital Signal
More informationLagrange Duality. Daniel P. Palomar. Hong Kong University of Science and Technology (HKUST)
Lagrange Duality Daniel P. Palomar Hong Kong University of Science and Technology (HKUST) ELEC5470 - Convex Optimization Fall 2017-18, HKUST, Hong Kong Outline of Lecture Lagrangian Dual function Dual
More informationIterative Learning Control Analysis and Design I
Iterative Learning Control Analysis and Design I Electronics and Computer Science University of Southampton Southampton, SO17 1BJ, UK etar@ecs.soton.ac.uk http://www.ecs.soton.ac.uk/ Contents Basics Representations
More information1396 IEEE TRANSACTIONS ON POWER DELIVERY, VOL. 24, NO. 3, JULY /$ IEEE
1396 IEEE TRANSACTIONS ON POWER DELIVERY, VOL. 24, NO. 3, JULY 2009 Fast Realization of the Modal Vector Fitting Method for Rational Modeling With Accurate Representation of Small Eigenvalues Bjørn Gustavsen,
More informationOptimum Ordering and Pole-Zero Pairing. Optimum Ordering and Pole-Zero Pairing Consider the scaled cascade structure shown below
Pole-Zero Pairing of the Cascade Form IIR Digital Filter There are many possible cascade realiations of a higher order IIR transfer function obtained by different pole-ero pairings and ordering Each one
More informationStability Condition in Terms of the Pole Locations
Stability Condition in Terms of the Pole Locations A causal LTI digital filter is BIBO stable if and only if its impulse response h[n] is absolutely summable, i.e., 1 = S h [ n] < n= We now develop a stability
More informationOptimum Ordering and Pole-Zero Pairing of the Cascade Form IIR. Digital Filter
Optimum Ordering and Pole-Zero Pairing of the Cascade Form IIR Digital Filter There are many possible cascade realiations of a higher order IIR transfer function obtained by different pole-ero pairings
More informationThe HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
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 informationModel reduction of large-scale systems by least squares
Model reduction of large-scale systems by least squares Serkan Gugercin Department of Mathematics, Virginia Tech, Blacksburg, VA, USA gugercin@mathvtedu Athanasios C Antoulas Department of Electrical and
More informationDesign of Stable IIR filters with prescribed flatness and approximately linear phase
Design of Stable IIR filters with prescribed flatness and approximately linear phase YASUNORI SUGITA Nagaoka University of Technology Dept. of Electrical Engineering Nagaoka city, Niigata-pref., JAPAN
More information2. Typical Discrete-Time Systems All-Pass Systems (5.5) 2.2. Minimum-Phase Systems (5.6) 2.3. Generalized Linear-Phase Systems (5.
. Typical Discrete-Time Systems.1. All-Pass Systems (5.5).. Minimum-Phase Systems (5.6).3. Generalized Linear-Phase Systems (5.7) .1. All-Pass Systems An all-pass system is defined as a system which has
More informationClosed-loop system 2/1/2016. Generally MIMO case. Two-degrees-of-freedom (2 DOF) control structure. (2 DOF structure) The closed loop equations become
Closed-loop system enerally MIMO case Two-degrees-of-freedom (2 DOF) control structure (2 DOF structure) 2 The closed loop equations become solving for z gives where is the closed loop transfer function
More informationTHEORY OF MIMO BIORTHOGONAL PARTNERS AND THEIR APPLICATION IN CHANNEL EQUALIZATION. Bojan Vrcelj and P. P. Vaidyanathan
THEORY OF MIMO BIORTHOGONAL PARTNERS AND THEIR APPLICATION IN CHANNEL EQUALIZATION Bojan Vrcelj and P P Vaidyanathan Dept of Electrical Engr 136-93, Caltech, Pasadena, CA 91125, USA E-mail: bojan@systemscaltechedu,
More information( ) John A. Quinn Lecture. ESE 531: Digital Signal Processing. Lecture Outline. Frequency Response of LTI System. Example: Zero on Real Axis
John A. Quinn Lecture ESE 531: Digital Signal Processing Lec 15: March 21, 2017 Review, Generalized Linear Phase Systems Penn ESE 531 Spring 2017 Khanna Lecture Outline!!! 2 Frequency Response of LTI System
More informationOn Input Design for System Identification
On Input Design for System Identification Input Design Using Markov Chains CHIARA BRIGHENTI Masters Degree Project Stockholm, Sweden March 2009 XR-EE-RT 2009:002 Abstract When system identification methods
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science. Fall Solutions for Problem Set 2
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Issued: Tuesday, September 5. 6.: Discrete-Time Signal Processing Fall 5 Solutions for Problem Set Problem.
More informationMinimax Design of Complex-Coefficient FIR Filters with Low Group Delay
Minimax Design of Complex-Coefficient FIR Filters with Low Group Delay Wu-Sheng Lu Takao Hinamoto Dept. of Elec. and Comp. Engineering Graduate School of Engineering University of Victoria Hiroshima University
More informationECE 451 Macromodeling
ECE 451 Macromodeling Jose E. Schutt-Aine Electrical & Computer Engineering University of Illinois jesa@illinois.edu ECE 451 Jose Schutt Aine 1 Blackbox Macromodeling Nonlinear Network 1 Nonlinear Network
More informationQuality metrics for S-parameter models
Quality metrics for S-parameter models DesignCon IBIS Summit, Santa Clara, February 4, 2010 Yuriy Shlepnev shlepnev@simberian.com Copyright 2010 by Simberian Inc. Reuse by written permission only. All
More informationUNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 9th, 011 Examination hours: 14.30 18.30 This problem set
More informationz-transform Chapter 6
z-transform Chapter 6 Dr. Iyad djafar Outline 2 Definition Relation Between z-transform and DTFT Region of Convergence Common z-transform Pairs The Rational z-transform The Inverse z-transform z-transform
More informationControl Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:
More informationEL2520 Control Theory and Practice
So far EL2520 Control Theory and Practice r Fr wu u G w z n Lecture 5: Multivariable systems -Fy Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden SISO control revisited: Signal
More informationDigital Filter Structures
Chapter 8 Digital Filter Structures 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 8-1 Block Diagram Representation The convolution sum description of an LTI discrete-time system can, in principle, be used to
More informationFilter Analysis and Design
Filter Analysis and Design Butterworth Filters Butterworth filters have a transfer function whose squared magnitude has the form H a ( jω ) 2 = 1 ( ) 2n. 1+ ω / ω c * M. J. Roberts - All Rights Reserved
More informationNeed for transformation?
Z-TRANSFORM In today s class Z-transform Unilateral Z-transform Bilateral Z-transform Region of Convergence Inverse Z-transform Power Series method Partial Fraction method Solution of difference equations
More informationTwo-Sided Arnoldi in Order Reduction of Large Scale MIMO Systems
Two-Sided Arnoldi in Order Reduction of Large Scale MIMO Systems B. Salimbahrami and B. Lohmann Institute of Automation, University of Bremen Otto-Hahn-Allee, NW1 28359 Bremen, Germany e-mail: salimbahrami@iat.uni-bremen.de
More informationOptimal Common Sub-Expression Elimination Algorithm of Multiple Constant Multiplications with a Logic Depth Constraint
3568 IEICE TRANS. FUNDAMENTALS, VOL.E91 A, NO.12 DECEMBER 2008 PAPER Special Section on VLSI Design and CAD Algorithms Optimal Common Sub-Expression Elimination Algorithm of Multiple Constant Multiplications
More informationZ-Transform. 清大電機系林嘉文 Original PowerPoint slides prepared by S. K. Mitra 4-1-1
Chapter 6 Z-Transform 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 4-1-1 z-transform The DTFT provides a frequency-domain representation of discrete-time
More informationParameterized Model Order Reduction via Quasi-Convex Optimization. Kin Cheong Sou with Luca Daniel and Alexandre Megretski
Parameterized Model Order Reduction via Quasi-Convex Optimization Kin Cheong Sou with Luca Daniel and Alexandre Megretski Systems on Chip or Package RF Inductors Interconnect & Substrate MEM resonators
More informationDiscrete-time signals and systems
Discrete-time signals and systems 1 DISCRETE-TIME DYNAMICAL SYSTEMS x(t) G y(t) Linear system: Output y(n) is a linear function of the inputs sequence: y(n) = k= h(k)x(n k) h(k): impulse response of the
More information(i) Represent discrete-time signals using transform. (ii) Understand the relationship between transform and discrete-time Fourier transform
z Transform Chapter Intended Learning Outcomes: (i) Represent discrete-time signals using transform (ii) Understand the relationship between transform and discrete-time Fourier transform (iii) Understand
More informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Transform The z-transform is a very important tool in describing and analyzing digital systems. It offers the techniques for digital filter design and frequency analysis of digital signals. Definition
More informationAnalysis Of Ill-Conditioning Of Multi-Channel Deconvolution Problems
Analysis Of Ill-Conditioning Of Multi-Channel Deconvolution Problems Ole Kirkeby, Per Rubak, and Angelo Farina * Department of Communication Techology, Fredrik Bajers Vej 7, Aalborg University, DK-9220
More informationEVALUATION OF (UNSTABLE) NON-CAUSAL SYSTEMS APPLIED TO ITERATIVE LEARNING CONTROL
EVALUATION OF (UNSTABLE) NON-CAUSAL SYSTEMS APPLIED TO ITERATIVE LEARNING CONTROL M.G.E. Schneiders, M.J.G. van de Molengraft and M. Steinbuch Eindhoven University of Technology, Control Systems Technology,
More informationRobust Control of Time-delay Systems
Robust Control of Time-delay Systems Qing-Chang Zhong Distinguished Lecturer, IEEE Power Electronics Society Max McGraw Endowed Chair Professor in Energy and Power Engineering Dept. of Electrical and Computer
More informationChapter 10: Rational Functions and the Riemann Sphere. By a rational function we mean a function f which can be expressed in the form
Chapter 10: Rational Functions and the Riemann Sphere By a rational function we mean a function f which can be expressed in the form f(z) = p(z) q(z) = a nz n +a n 1 z n 1 + +a 1 z +a 0 b m z m +b m 1
More information# FIR. [ ] = b k. # [ ]x[ n " k] [ ] = h k. x[ n] = Ae j" e j# ˆ n Complex exponential input. [ ]Ae j" e j ˆ. ˆ )Ae j# e j ˆ. y n. y n.
[ ] = h k M [ ] = b k x[ n " k] FIR k= M [ ]x[ n " k] convolution k= x[ n] = Ae j" e j ˆ n Complex exponential input [ ] = h k M % k= [ ]Ae j" e j ˆ % M = ' h[ k]e " j ˆ & k= k = H (" ˆ )Ae j e j ˆ ( )
More informationLTI Systems, Additive Noise, and Order Estimation
LTI Systems, Additive oise, and Order Estimation Soosan Beheshti, Munther A. Dahleh Laboratory for Information and Decision Systems Department of Electrical Engineering and Computer Science Massachusetts
More informationConvex Optimization Problems. Prof. Daniel P. Palomar
Conve Optimization Problems Prof. Daniel P. Palomar The Hong Kong University of Science and Technology (HKUST) MAFS6010R- Portfolio Optimization with R MSc in Financial Mathematics Fall 2018-19, HKUST,
More informationMulti-rate Signal Processing 7. M-channel Maximally Decmiated Filter Banks
Multi-rate Signal Processing 7. M-channel Maximally Decmiated Filter Banks Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes
More informationNew Design of Orthogonal Filter Banks Using the Cayley Transform
New Design of Orthogonal Filter Banks Using the Cayley Transform Jianping Zhou, Minh N. Do and Jelena Kovačević Department of Electrical and Computer Engineering University of Illinois at Urbana-Champaign,
More informationInitial Temperature Reconstruction for a Nonlinear Heat Equation: Application to Radiative Heat Transfer.
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 Initial Temperature Reconstruction for a Nonlinear Heat Equation:
More informationEnforcement Passivity. Frequency Data. Wenliang Tseng, Sogo Hsu, Frank Y.C. Pai and Scott C.S. Li. Asian IBIS Summit, Taipei, Taiwan November 12, 2010
Enforcement Passivity of S-parameter S Sampled Frequency Data Wenliang Tseng, Sogo Hsu, Frank Y.C. Pai and Scott C.S. Li Asian IBIS Summit, Taipei, Taiwan November 12, 2010 Agenda Causality and passivity
More informationAnalysis of TMT Primary Mirror Control-Structure Interaction (SPIE )
Analysis of TMT Primary Mirror Control-Structure Interaction (SPIE 7017-41) Douglas MacMynowski (Caltech) Peter Thompson (Systems Tech.) Mark Sirota (TMT Observatory) Control Problems TMT.SEN.PRE.08.046.REL01
More informationECE503: Digital Signal Processing Lecture 5
ECE53: Digital Signal Processing Lecture 5 D. Richard Brown III WPI 3-February-22 WPI D. Richard Brown III 3-February-22 / 32 Lecture 5 Topics. Magnitude and phase characterization of transfer functions
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 informationA Comparative Study on Automatic Flight Control for small UAV
Proceedings of the 5 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'18) Niagara Falls, Canada June 7 9, 18 Paper No. 13 DOI: 1.11159/cdsr18.13 A Comparative Study on Automatic
More informationChapter 7: IIR Filter Design Techniques
IUST-EE Chapter 7: IIR Filter Design Techniques Contents Performance Specifications Pole-Zero Placement Method Impulse Invariant Method Bilinear Transformation Classical Analog Filters DSP-Shokouhi Advantages
More informationAnother Walkthrough of Variational Bayes. Bevan Jones Machine Learning Reading Group Macquarie University
Another Walkthrough of Variational Bayes Bevan Jones Machine Learning Reading Group Macquarie University 2 Variational Bayes? Bayes Bayes Theorem But the integral is intractable! Sampling Gibbs, Metropolis
More informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture : Design of Digital IIR Filters (Part I) Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 008 K. E. Barner (Univ.
More informationPower Transformer Transient Modeling Using Frequency Response Analysis
Power Transformer Transient Modeling Using Frequency Response Analysis by Hosam Salem Alharbi A thesis submitted to the Faculty of Graduate Studies of The University of Manitoba in partial fulfilment of
More informationLecture 9 Infinite Impulse Response Filters
Lecture 9 Infinite Impulse Response Filters Outline 9 Infinite Impulse Response Filters 9 First-Order Low-Pass Filter 93 IIR Filter Design 5 93 CT Butterworth filter design 5 93 Bilinear transform 7 9
More informationDiscrete-time Symmetric/Antisymmetric FIR Filter Design
Discrete-time Symmetric/Antisymmetric FIR Filter Design Presenter: Dr. Bingo Wing-Kuen Ling Center for Digital Signal Processing Research, Department of Electronic Engineering, King s College London. Collaborators
More information