Application of Taylor Models to the Worst-Case Analysis of Stripline Interconnects
|
|
- Clifford Charles
- 5 years ago
- Views:
Transcription
1 Application of Taylor Models to the Worst-Case Analysis of Stripline Interconnects Paolo Manfredi, Riccardo Trinchero, Igor Stievano, Flavio Canavero 1
2 Motivation q Ever-growing interest in the inclusion of parameters uncertainties in circuit simulations q Variability due to unknown components, manufacturing tolerances, uncertain design parameters q Key resource in the design phase parameters! #! $! " U(s) input CIRCUIT / SYSTEM Y(s,x) random output 2
3 Worst Case q Determining the worst case configuration is fundamental in a number of critical applications (transportation, industrial, medical, etc ) q Standard methods require a large number of deterministic simulations to guarantee an accurate estimation (the worst case often lies in regions of low probability!) Deterministic Simulations samples In principle an infinite number of samples is needed to explore all the possible combinations of the parameters 3
4 State-of-the-art (features) Monte Carlo Polynomial Chaos Interval methods Interval Arithmetic Affine Arithmetic Taylor Models Easy to implement Slow convergence (> samples) Underestimation (i.e., approx. of the bound from inside ) Blind analysis (no parametric model) Nice theoretical framework Improved efficiency No info on error: under/over estimation?? Inherently worst-case Promising alternative that worth to be carefully investigated 4
5 Problem statement q Numerical techniques that propagate bounded variables by redefining the standard operations producing guaranteed enclosures q E.g., the analysis of a circuit defined by components with bounded values requires to handle operations over bounded variables in a conservative way: MNA: 1 ' # ()) +, ()) = / + 1) )4 5# 7()) 0 linear operations addition/subtraction, product,.. nonlinear operations division, exponential, powers, etc complex operations matrix operations Addition/subtraction, product, inversion In order to produce a guaranteed enclosure of the output responses, each operation must provide a conservative bound of its result 5
6 Interval Methods [1] Ø Interval Arithmetic (IA) Random parameters = bounded interval e.g., subtraction: 9, ; <, = = 9 =, ; < 2,3 2,3 = 2 3,3 2 = [ 1,1] 0! [9, ;] 9 ;! No correlation among variables Overestimation, i.e. TOO PESSIMISTIC Ø Affine Arithmetic (AA) Random parameters = linear model via new symbol?!@ =! A +! #? # +! "? " ; 9 Accounts for correlation among parameters 1 1 ℇ Extension to nonlinear operations Non-trivial (even for the product) [1] T. Ding, R. Trinchero, P. Manfredi, I. S. Stievano and F. G. Canavero, "How Affine Arithmetic Helps Beat Uncertainties in Electrical Systems," in IEEE Circuits and Syst. Mag., vol. 15, no. 4, pp ,
7 Taylor model q Hybrid representation of a nonlinear function of an interval variable in terms of a Taylor expansion plus an interval remainder [2]. q Given! = [9, ;] J! N O! = P O! + + O [K, L] a Lagrange Remainder [2] M. Berz and K. Makino, Verified integration of ODEs and flows using differential algebraic methods on high-order Taylor models, Rel. Computing, vol. 4, no. 4, pp , ! I + 3! " b! P(x) gives parametric model of f(x) Lagrange remainder provides conservative error estimation
8 Worst-case bounds J(!) Q J! = [min J!, max J(!)] N O = P O! + + O Q(N O ) = Q(P O! ) ++ O True worst case Q(P O! ) ++ O a b Estimated bound: guaranteed conservative!! 8
9 Polynomial bound: Bernstein q Polynomial bound calculation is non-trivial for multivariate polynomials with order > 2 q Effective solution: conversion of P O into Bernstain polynomials [3] P O (!) X ; Y Q Y (!) q Bernstein polynomials have a notable property: min ; Y P O (!) max ; Y [3] R. Trinchero, P. Manfredi, T. Ding, and I. S. Stievano, Combined parametric and worst-case circuit analysis via Taylor models," accepted for publication in IEEE Trans. Circuits and Syst. (2016). 9
10 Taylor model algebra (basic rules) Sums, subtractions (straightforward)! + [1,2] ! + [ 1,1] = 3 + 3! + 0,3 Multiplication e.g., order 1! + [1,2] \ 3 + 2! + [ 1,1] sum of poly computed via IA e.g., order 1 = 3! + 2! " +! 1,1 + 1, ! + 1,1 1,2 B(2! " ) B(!) B(3 + 2!) computed via IA = 3! + [ 2,15] 10
11 Nonlinear operators q The functional Taylor expansion allows replacing nonlinear operators (e.g., ] N O = 1 N O ) with a mere series of products and sums Sum of scalar coefficients and products b ] N O = P^! + +^ = X ] _ < O `! _ca available!! _ N O < O + Q(/ b (d)) Taylor expansion of ] around Taylor model center < O = P O (0) Lagrange Remainder / b d = ] bg# d h + 1! N O < O bg# with d Q(N O ) [4] [4] M.A. Abramowitz and I.A. Stegun, Handbook of mathematical functions. New York, NY, USA: Dover Publications,
12 Complex and Matrix Operations q Complex calculations are reduced to real-valued operations on Re/Im E.g. 9 r + 19 s ; r + 1; s = 9 r ; r 9 s ; s + 1(9 r ; s + 9 s ; r ) q A Taylor model matrix is a matrix where all the entries are Taylor models j k = N k,## N k,#" N k,"# N k,"" N k,## = P k,##! + + k,## q Linear matrix operations: operate element-wise with scalar operations Matrix Sum/Subtraction j k ± j p Yq = N k,yq ± N p,yq scalar addition Matrix Product j k j p Yq = X N k,y_ N p,_q _ scalar multiplication 12
13 Matrix Inversion q How to invert a matrix TM and still obtain a conservative remainder? ØIdea: Sherman-Morrison formula (SM) t + u 5# = t 5# tr t 5# u (t5# ut 5# ) ü exact solution (NO approximation) ü only requires inversion of t q Requirements t must be (easily) invertible u must be rank-1 suitable iterative extension available 13
14 Solution [3]: Matrix Inversion (cont d) j O = 7 + Q constant matrix (center of the expansion) TM matrix (Poly + remainder) (j x ) 5# = t 5# tr t 5# u (t5# ut 5# ) standard matrix inversion available TM operations poly+remainder handled as it is NO INV NEEDED Rigorous solution for the inversion of an interval matrix!! [3] R. Trinchero, P. Manfredi, T. Ding, and I. S. Stievano, Combined parametric and worst-case circuit analysis via Taylor models," accepted for publication in IEEE Trans. Circuits and Syst. (2016). 14
15 Example (from [4]) R 1 R 6 R 11 R 16 R 21 In R 2 R 3 C 1 R 7 R 8 C 3 R 12 R 13 C 5 R 17 R 18 C 7 C 9 Out Modified Nodal Analysis y + z{ = } R 4 R 9 R 14 R 19 R 5 C 2 R 1 =5.4779kΩ R 6 =4.44kΩ R 15 = kΩ C 1 =12nF R 10 R 2 =2.0076kΩ R 7 =5.9999kΩ R 16 = kΩ C 3 =6.8nF C 4 Solve the system j = j y + zj { 5# } and obtain the WC R 15 C 6 R 3,4,8,9,13,14,18,19 =3.3kΩ R 10 =4.2573kΩ R 17 =1.0301kΩ C 5 =4.7nF R 20 R 11 =3.2201kΩ R 20 = kΩ C 7 =6.8nF C 8 R 5 =4.5898kΩ R 12 = kΩ R 21 =1.2201kΩ C 2,4,6,8,9 =10nF TM Matrices y j y { j { Defining the bounded parameters 2 " = 2 ",A ( ! # ) 2 = 2,A ! # 2 Ä = 2 Ä,A ( ! " ) 2 Å = 2 Å,A ! " 2 Ç = 2 Ç,A ! I where,! Y [ 1,1] [4] T.-A. Pham, E. Gad, M. S. Nakhla, and R. Achar, Decoupled polynomial chaos and its applications to statistical analysis of high-speed interconnects, IEEE Trans. Compon. Packag. Manuf. Techol., 4(10) (2014). 15
16 Application of Taylor Model to the Worst-Case Analysis of Stripline Interconnects SPI2016 Torino, May 8-11 Results (see [3]) q Worst-Case Analysis q Parametric Analysis 1.5 Pass band Ripple MC TM ( or de r 5) Frequency [Hz] MonteCarlo (10000 runs) Very Good accuracy Bernstein polys yield tight bounds TM offers a parameterized description! [3] R. Trinchero, P. Manfredi, T. Ding, and I. S. Stievano, Combined parametric and worst-case circuit analysis via Taylor models," accepted for publication in IEEE Trans. Circuits and Syst. (2016). 16
17 Application example: Coupled Stripline q First application of TMs to the WC analysis of a transmission line q Coupled stripline [4] in a homogeneous dielectric medium with relative permittivity ε r 0.6 mil 5 mil 5 mil 5 mil ε r =4 20 mil q Defining the following bounded parameters 50Ω 10 cm 1pF? É =? É,A ± Δ ÖÜ = 4 ± à = 2 à,a ± Δ âä = 1 ± 0.05 pf 50Ω 1pF [4] T.-A. Pham, E. Gad, M. S. Nakhla, and R. Achar, Decoupled polynomial chaos and its applications to statistical analysis of high-speed interconnects, IEEE Trans. Compon. Packag. Manuf. Techol., 4(10) (2014). 17
18 Application example: Coupled Stripline (ii) q The far- and near-end voltage and current variables are computed through the chain-parameter matrix: å ## = cos ïê ñ å #" = 1ó A sin ïê ò å "# = 1ô A sin ïê { å "" = cos(ïê) ñ ã 0 = å "# ç é + è, å #" å "" è, å ## ç 5# é \ (å "# è, å ## ) 0 = é ç é ã 0 ê = å ## é + å #" å ## ç é ã 0 ë ê = è, (ê) q The WC analysis requires to handle in the TM framework: Ø matrix inversion available Ø sin/cos operations details in SPI paper! 18
19 Application example: Coupled Stripline (iii) q Worst-Case Analysis q Parametric model for ö õ Far-end voltage [db] Far-end crosstalk [db] SPICE (MC) Taylor model Frequency [GHz] Far-end voltage [db] Far-end crosstalk [db] SPICE Taylor model Frequency [GHz] ε r ε r SPICE MonteCarlo (10000 runs): 680 s TM (order 3): 35.6 s 19
20 Conclusions q Promising solution for the direct worst-case analysis of the frequency-domain response of circuits and interconnects (key tool in critical applications) qtm-based approach: Ø represents any variable/response as a polynomial + interval remainder (for approximation and round-off errors) Ø Elegant and robust mathematical framework Ø Extended to handle matrix operations (inversion!!) Ø allows to jointly obtain a parameterized response and WC bounds ØTighter bounds w.r.t. other interval approaches Ø promising results in terms of both accuracy and efficiency 20
21 Thank you for your attention Q&A 21
Prediction of Stochastic Eye Diagrams via IC Equivalents and Lagrange Polynomials
Prediction of Stochastic Eye Diagrams via IC Equivalents and Lagrange Polynomials Paolo Manfredi, Igor S. Stievano, Flavio G. Canavero Department of Electronics and Telecommunications (DET) Politecnico
More informationSplit and Merge Strategies for Solving Uncertain Equations Using Affine Arithmetic
Split and Merge Strategies for Solving Uncertain Equations Using Affine Arithmetic Oliver Scharf, Markus Olbrich, Erich Barke SIMUTOOLS 2015 Outline Introduction Uncertain parameters Affine arithmetic
More informationInterval Methods and Taylor Model Methods for ODEs
Interval Methods and Taylor Model Methods for ODEs Markus Neher, Dept. of Mathematics KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 TM VII, Key West KIT University of the State of Baden-Wuerttemberg and Interval
More informationAn Introduction to Differential Algebra
An Introduction to Differential Algebra Alexander Wittig1, P. Di Lizia, R. Armellin, et al. 1 ESA Advanced Concepts Team (TEC-SF) SRL, Milan Dinamica Outline 1 Overview Five Views of Differential Algebra
More informationSOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC
SOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC Alfredo Vaccaro RTSI 2015 - September 16-18, 2015, Torino, Italy RESEARCH MOTIVATIONS Power Flow (PF) and
More informationStatistical approach in complex circuits by a wavelet based Thevenin s theorem
Proceedings of the 11th WSEAS International Conference on CIRCUITS, Agios Nikolaos, Crete Island, Greece, July 23-25, 2007 185 Statistical approach in complex circuits by a wavelet based Thevenin s theorem
More informationPolitecnico di Torino. Porto Institutional Repository
Politecnico di Torino Porto Institutional Repository [Doctoral thesis] Worst-Case Analysis of Electrical and Electronic Equipment via Affine Arithmetic Original Citation: Tongyu Ding (2015). Arithmetic.
More informationBroadband Electromagnetic and Stochastic Modelling and Signal Analysis of Multiconductor Interconnections
Broadband Electromagnetic and Stochastic Modelling and Signal Analysis of Multiconductor Interconnections Electrical & Computer Engineering University of Toronto Distinguished Lecture Series 2017-2018
More informationON TAYLOR MODEL BASED INTEGRATION OF ODES
ON TAYLOR MODEL BASED INTEGRATION OF ODES M. NEHER, K. R. JACKSON, AND N. S. NEDIALKOV Abstract. Interval methods for verified integration of initial value problems (IVPs for ODEs have been used for more
More informationVerified High-Order Optimal Control in Space Flight Dynamics
Verified High-Order Optimal Control in Space Flight Dynamics R. Armellin, P. Di Lizia, F. Bernelli-Zazzera K. Makino and M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16
More informationDeterministic Global Optimization for Dynamic Systems Using Interval Analysis
Deterministic Global Optimization for Dynamic Systems Using Interval Analysis Youdong Lin and Mark A. Stadtherr Department of Chemical and Biomolecular Engineering University of Notre Dame, Notre Dame,
More informationSecond Order ODEs. CSCC51H- Numerical Approx, Int and ODEs p.130/177
Second Order ODEs Often physical or biological systems are best described by second or higher-order ODEs. That is, second or higher order derivatives appear in the mathematical model of the system. For
More informationRigorous Global Optimization of Impulsive Space Trajectories
Rigorous Global Optimization of Impulsive Space Trajectories P. Di Lizia, R. Armellin, M. Lavagna K. Makino, M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16 19, 2006 Motivation
More informationECE 598 JS Lecture 06 Multiconductors
ECE 598 JS Lecture 06 Multiconductors Spring 2012 Jose E. Schutt-Aine Electrical & Computer Engineering University of Illinois jesa@illinois.edu 1 TELGRAPHER S EQUATION FOR N COUPLED TRANSMISSION LINES
More informationA Rigorous Implementation of Taylor Model-Based Integrator in Chebyshev Basis Representation 1
A Rigorous Implementation of Taylor Model-Based Integrator in Chebyshev Basis Representation 1 Tomáš Dzetkulič Institute of Computer Science Czech Academy of Sciences December 15, 2011 1 This work was
More informationCombined FDTD/Macromodel Simulation of Interconnected Digital Devices
Combined FDTD/Macromodel Simulation of Interconnected Digital Devices S. Grivet-Talocia, I. S. Stievano, I. A. Maio, F. G. Canavero Dip. Elettronica, Politecnico di Torino, Torino, Italy (E-mail: grivet@polito.it)
More informationValidated Solutions of Initial Value Problems for Parametric ODEs
Validated Solutions of Initial Value Problems for Parametric ODEs Youdong Lin and Mark A. Stadtherr Department of Chemical and Biomolecular Engineering University of Notre Dame, Notre Dame, IN 46556, USA
More informationLinear Algebra Practice Problems
Math 7, Professor Ramras Linear Algebra Practice Problems () Consider the following system of linear equations in the variables x, y, and z, in which the constants a and b are real numbers. x y + z = a
More informationTwo-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra
Two-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra Politecnico di Milano Department of Aerospace Engineering Milan, Italy Taylor Methods and Computer Assisted Proofs
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 informationAccurate Computation of Vector Potentials in Lossy Media
Accurate Computation of Vector Potentials in Lossy Media Swagato Chakraborty and Vikram Jandhyala {swagato,jandhyala}@ee.washington.edu Dept of EE, University of Washington Seattle WA, 98195-500 UWEE Technical
More informationHyperbolic Polynomial Chaos Expansion (HPCE) and its Application to Statistical Analysis of Nonlinear Circuits
Hyperbolic Polynomial Chaos Expansion HPCE and its Application to Statistical Analysis of Nonlinear Circuits Majid Ahadi, Aditi Krishna Prasad, Sourajeet Roy High Speed System Simulations Laboratory Department
More informationL 2 -induced Gains of Switched Systems and Classes of Switching Signals
L 2 -induced Gains of Switched Systems and Classes of Switching Signals Kenji Hirata and João P. Hespanha Abstract This paper addresses the L 2-induced gain analysis for switched linear systems. We exploit
More informationTuning PI controllers in non-linear uncertain closed-loop systems with interval analysis
Tuning PI controllers in non-linear uncertain closed-loop systems with interval analysis J. Alexandre dit Sandretto, A. Chapoutot and O. Mullier U2IS, ENSTA ParisTech SYNCOP April 11, 2015 Closed-loop
More informationCSE 190: Reinforcement Learning: An Introduction. Chapter 8: Generalization and Function Approximation. Pop Quiz: What Function Are We Approximating?
CSE 190: Reinforcement Learning: An Introduction Chapter 8: Generalization and Function Approximation Objectives of this chapter: Look at how experience with a limited part of the state set be used to
More informationComputation of Interval Extensions Using Berz-Taylor Polynomial Models
Computation of Interval Extensions Using Berz-Taylor Polynomial Models Hongkun Liang and Mark A. Stadtherr Λ Department of Chemical Engineering University of Notre Dame Notre Dame, IN 46556 USA AIChE Annual
More informationG H. Extended Unit Tests A L L. Higher Still Advanced Higher Mathematics. (more demanding tests covering all levels) Contents. 3 Extended Unit Tests
M A T H E M A T I C S H I G H E R Higher Still Advanced Higher Mathematics S T I L L Extended Unit Tests A (more demanding tests covering all levels) Contents Extended Unit Tests Detailed marking schemes
More informationPreviously Monte Carlo Integration
Previously Simulation, sampling Monte Carlo Simulations Inverse cdf method Rejection sampling Today: sampling cont., Bayesian inference via sampling Eigenvalues and Eigenvectors Markov processes, PageRank
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationNumerical solutions of nonlinear systems of equations
Numerical solutions of nonlinear systems of equations Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan E-mail: min@math.ntnu.edu.tw August 28, 2011 Outline 1 Fixed points
More informationF O R SOCI AL WORK RESE ARCH
7 TH EUROPE AN CONFERENCE F O R SOCI AL WORK RESE ARCH C h a l l e n g e s i n s o c i a l w o r k r e s e a r c h c o n f l i c t s, b a r r i e r s a n d p o s s i b i l i t i e s i n r e l a t i o n
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationMATH 423 Linear Algebra II Lecture 3: Subspaces of vector spaces. Review of complex numbers. Vector space over a field.
MATH 423 Linear Algebra II Lecture 3: Subspaces of vector spaces. Review of complex numbers. Vector space over a field. Vector space A vector space is a set V equipped with two operations, addition V V
More informationPulse Shaping and ISI (Proakis: chapter 10.1, 10.3) EEE3012 Spring 2018
Pulse Shaping and ISI (Proakis: chapter 10.1, 10.3) EEE3012 Spring 2018 Digital Communication System Introduction Bandlimited channels distort signals the result is smeared pulses intersymol interference
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon s method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationRigorous Polynomial Approximations and Applications
Rigorous Polynomial Approximations and Applications Mioara Joldeș under the supervision of: Nicolas Brisebarre and Jean-Michel Muller École Normale Supérieure de Lyon, Arénaire Team, Laboratoire de l Informatique
More informationPropagation of Uncertainties in Nonlinear Dynamic Models
Propagation of Uncertainties in Nonlinear Dynamic Models Youdong Lin 1, Scott Ferson 2, George F. Corliss 3 and Mark A. Stadtherr 1 1 Department of Chemical and Biomolecular Engineering, University of
More informationTaylor Model Range Bounding Schemes
Taylor Model Range Bounding Schemes I Kyoko Makino II Martin Berz Department of Physics and Astronomy Michigan State University Introduction LDB (Linear Dominated Bounder) QDB (Quadratic Dominated Bounder)
More informationTHEORETICAL ANALYSIS OF Ä È MEAN COMPARATORS. Constantine Kotropoulos and Ioannis Pitas
THEORETICAL ANALYSIS OF Ä È MEAN COMPARATORS Constantine Kotropoulos and Ioannis Pitas Dept. of Informatics, Aristotle University of Thessaloniki Box 451, Thessaloniki 54 6, GREECE costas,pitas@zeus.csd.auth.gr
More informationRandom Eigenvalue Problems in Structural Dynamics: An Experimental Investigation
Random Eigenvalue Problems in Structural Dynamics: An Experimental Investigation S. Adhikari, A. Srikantha Phani and D. A. Pape School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk
More informationHigh-Order Representation of Poincaré Maps
High-Order Representation of Poincaré Maps Johannes Grote, Martin Berz, and Kyoko Makino Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA. {grotejoh,berz,makino}@msu.edu
More informationChapter 2. Matrix Arithmetic. Chapter 2
Matrix Arithmetic Matrix Addition and Subtraction Addition and subtraction act element-wise on matrices. In order for the addition/subtraction (A B) to be possible, the two matrices A and B must have the
More informationVerified Global Optimization with Taylor Model based Range Bounders
Verified Global Optimization with Taylor Model based Range Bounders KYOKO MAKINO and MARTIN BERZ Department of Physics and Astronomy Michigan State University East Lansing, MI 88 USA makino@msu.edu, berz@msu.edu
More informationEssential Mathematics for Economists
Ludwig von Auer (Universität Trier) Marco Raphael (Universität Trier) October 2014 1 / 343 1 Introductory Topics I: Algebra and Equations Some Basic Concepts and Rules How to Solve Simple Equations Equations
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationLecture XI. Approximating the Invariant Distribution
Lecture XI Approximating the Invariant Distribution Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Invariant Distribution p. 1 /24 SS Equilibrium in the Aiyagari model G.
More informationCase Study: Parallel Coupled- Line Combline Filter
MICROWAVE AND RF DESIGN MICROWAVE AND RF DESIGN Case Study: Parallel Coupled- Line Combline Filter Presented by Michael Steer Reading: 6. 6.4 Index: CS_PCL_Filter Based on material in Microwave and RF
More informationAN INTERVAL METHOD FOR INITIAL VALUE PROBLEMS IN LINEAR ORDINARY DIFFERENTIAL EQUATIONS
AN INTERVAL METHOD FOR INITIAL VALUE PROBLEMS IN LINEAR ORDINARY DIFFERENTIAL EQUATIONS NEDIALKO S. NEDIALKOV Abstract. Interval numerical methods for initial value problems (IVPs) for ordinary differential
More informationDerivation of the Kalman Filter
Derivation of the Kalman Filter Kai Borre Danish GPS Center, Denmark Block Matrix Identities The key formulas give the inverse of a 2 by 2 block matrix, assuming T is invertible: T U 1 L M. (1) V W N P
More information1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det
What is the determinant of the following matrix? 3 4 3 4 3 4 4 3 A 0 B 8 C 55 D 0 E 60 If det a a a 3 b b b 3 c c c 3 = 4, then det a a 4a 3 a b b 4b 3 b c c c 3 c = A 8 B 6 C 4 D E 3 Let A be an n n matrix
More informationAccurate numerical orbit propagation using Polynomial Algebra Computational Engine PACE. ISSFD 2015 congress, Germany. Dated: September 14, 2015
Accurate numerical orbit propagation using Polynomial Algebra Computational Engine PACE Emmanuel Bignon (), Pierre mercier (), Vincent Azzopardi (), Romain Pinède () ISSFD 205 congress, Germany Dated:
More informationFast evaluation of mixed derivatives and calculation of optimal weights for integration. Hernan Leovey
Fast evaluation of mixed derivatives and calculation of optimal weights for integration Humboldt Universität zu Berlin 02.14.2012 MCQMC2012 Tenth International Conference on Monte Carlo and Quasi Monte
More informationEffects from the Thin Metallic Substrate Sandwiched in Planar Multilayer Microstrip Lines
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 115 Effects from the Thin Metallic Substrate Sandwiched in Planar Multilayer Microstrip Lines L. Zhang and J. M. Song Iowa
More informationFrom Particle Accelerators to Celestial Dynamics: An Introduction to Differential Algebra
AstroNet II Summer School 2014 Friday, 27 June 2014 Zielona Góra (Poland) From Particle Accelerators to Celestial Dynamics: An Introduction to Differential Algebra 1, Pierluigi di Lizia 1, Roberto Armellin
More informationInterval Arithmetic An Elementary Introduction and Successful Applications
Interval Arithmetic An Elementary Introduction and Successful Applications by R. Baker Kearfott Department of Mathematics University of Southwestern Louisiana U.S.L. Box 4-1010 Lafayette, LA 70504-1010
More informationGlobal Optimization of Impulsive Interplanetary Transfers
Global Optimization of Impulsive Interplanetary Transfers R. Armellin, Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano Taylor Methods and Computer Assisted Proofs Barcelona, June, 3 7, 2008
More informationCompute the behavior of reality even if it is impossible to observe the processes (for example a black hole in astrophysics).
1 Introduction Read sections 1.1, 1.2.1 1.2.4, 1.2.6, 1.3.8, 1.3.9, 1.4. Review questions 1.1 1.6, 1.12 1.21, 1.37. The subject of Scientific Computing is to simulate the reality. Simulation is the representation
More informationA Polynomial Chaos Approach to Robust Multiobjective Optimization
A Polynomial Chaos Approach to Robust Multiobjective Optimization Silvia Poles 1, Alberto Lovison 2 1 EnginSoft S.p.A., Optimization Consulting Via Giambellino, 7 35129 Padova, Italy s.poles@enginsoft.it
More informationCompact Distributed RLC Interconnect Models Part I: Single Line Transient, Time Delay, and Overshoot Expressions
2068 IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. 47, NO. 11, NOVEMBER 2000 Compact Distributed RLC Interconnect Models Part I: Single Line Transient, Time Delay, and Overshoot Expressions Jeffrey A. Davis
More informationEnabling Advanced Automation Tools to manage Trajectory Prediction Uncertainty
Engineering, Test & Technology Boeing Research & Technology Enabling Advanced Automation Tools to manage Trajectory Prediction Uncertainty ART 12 - Automation Enrique Casado (BR&T-E) enrique.casado@boeing.com
More informationAffine Arithmetic: Concepts and Applications
SCAN 2002 Affine Arithmetic: Concepts and Applications Luiz Henrique de Figueiredo (IMPA) Jorge Stolfi (UNICAMP) Outline What is affine arithmetic? The dependency problem in interval arithmetic Main concepts
More informationEstimating functional uncertainty using polynomial chaos and adjoint equations
0. Estimating functional uncertainty using polynomial chaos and adjoint equations February 24, 2011 1 Florida State University, Tallahassee, Florida, Usa 2 Moscow Institute of Physics and Technology, Moscow,
More informationCONTENTS. Preface Preliminaries 1
Preface xi Preliminaries 1 1 TOOLS FOR ANALYSIS 5 1.1 The Completeness Axiom and Some of Its Consequences 5 1.2 The Distribution of the Integers and the Rational Numbers 12 1.3 Inequalities and Identities
More informationOur goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always
October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.
More informationROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino
ROBOTICS 01PEEQW Basilio Bona DAUIN Politecnico di Torino Probabilistic Fundamentals in Robotics Gaussian Filters Course Outline Basic mathematical framework Probabilistic models of mobile robots Mobile
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More information15-780: LinearProgramming
15-780: LinearProgramming J. Zico Kolter February 1-3, 2016 1 Outline Introduction Some linear algebra review Linear programming Simplex algorithm Duality and dual simplex 2 Outline Introduction Some linear
More informationModel 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 informationCayley-Hamilton Theorem
Cayley-Hamilton Theorem Massoud Malek In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n Let A be an n n matrix Although det (λ I n A
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 informationAccounting for Variability and Uncertainty in Signal and Power Integrity Modeling
Accounting for Variability and Uncertainty in Signal and Power Integrity Modeling Andreas Cangellaris& Prasad Sumant Electromagnetics Lab, ECE Department University of Illinois Acknowledgments Dr. Hong
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon's method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationVerification of Invertibility of Complicated Functions over Large Domains
Reliable Computing 8: 67 82, 2002. 67 c 2002 Kluwer Academic Publishers. Printed in the Netherlands. Verification of Invertibility of Complicated Functions over Large Domains JENS HOEFKENS and MARTIN BERZ
More informationOp-Amp Circuits: Part 3
Op-Amp Circuits: Part 3 M. B. Patil mbpatil@ee.iitb.ac.in www.ee.iitb.ac.in/~sequel Department of Electrical Engineering Indian Institute of Technology Bombay Introduction to filters Consider v(t) = v
More informationMethods in Computer Vision: Introduction to Matrix Lie Groups
Methods in Computer Vision: Introduction to Matrix Lie Groups Oren Freifeld Computer Science, Ben-Gurion University June 14, 2017 June 14, 2017 1 / 46 Definition and Basic Properties Definition (Matrix
More informationInverses of Square Matrices
Inverses of Square Matrices A. Havens Department of Mathematics University of Massachusetts, Amherst February 23-26, 2018 Outline 1 Defining Inverses Inverses for Products and Functions Defining Inverse
More informationChapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of
Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple
More informationErrors. Intensive Computation. Annalisa Massini 2017/2018
Errors Intensive Computation Annalisa Massini 2017/2018 Intensive Computation - 2017/2018 2 References Scientific Computing: An Introductory Survey - Chapter 1 M.T. Heath http://heath.cs.illinois.edu/scicomp/notes/index.html
More informationAdaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties
Australian Journal of Basic and Applied Sciences, 3(1): 308-322, 2009 ISSN 1991-8178 Adaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties M.R.Soltanpour, M.M.Fateh
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
3 Numerical Solution of Nonlinear Equations and Systems 3.1 Fixed point iteration Reamrk 3.1 Problem Given a function F : lr n lr n, compute x lr n such that ( ) F(x ) = 0. In this chapter, we consider
More informationA Robustness Optimization of SRAM Dynamic Stability by Sensitivity-based Reachability Analysis
ASP-DAC 2014 A Robustness Optimization of SRAM Dynamic Stability by Sensitivity-based Reachability Analysis Yang Song, Sai Manoj P. D. and Hao Yu School of Electrical and Electronic Engineering, Nanyang
More informationChapter 1 Error Analysis
Chapter 1 Error Analysis Several sources of errors are important for numerical data processing: Experimental uncertainty: Input data from an experiment have a limited precision. Instead of the vector of
More informationRobust Tuning of Power System Stabilizers Using Coefficient Diagram Method
International Journal of Electrical Engineering. ISSN 0974-2158 Volume 7, Number 2 (2014), pp. 257-270 International Research Publication House http://www.irphouse.com Robust Tuning of Power System Stabilizers
More informationNumerical Algorithms. IE 496 Lecture 20
Numerical Algorithms IE 496 Lecture 20 Reading for This Lecture Primary Miller and Boxer, Pages 124-128 Forsythe and Mohler, Sections 1 and 2 Numerical Algorithms Numerical Analysis So far, we have looked
More informationAPPROXIMATE SOLUTION OF A SYSTEM OF LINEAR EQUATIONS WITH RANDOM PERTURBATIONS
APPROXIMATE SOLUTION OF A SYSTEM OF LINEAR EQUATIONS WITH RANDOM PERTURBATIONS P. Date paresh.date@brunel.ac.uk Center for Analysis of Risk and Optimisation Modelling Applications, Department of Mathematical
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationGENG2140, S2, 2012 Week 7: Curve fitting
GENG2140, S2, 2012 Week 7: Curve fitting Curve fitting is the process of constructing a curve, or mathematical function, f(x) that has the best fit to a series of data points Involves fitting lines and
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Computer Representation of Numbers Counting numbers (unsigned integers) are the numbers 0,
More informationChapter 8: Taylor s theorem and L Hospital s rule
Chapter 8: Taylor s theorem and L Hospital s rule Theorem: [Inverse Mapping Theorem] Suppose that a < b and f : [a, b] R. Given that f (x) > 0 for all x (a, b) then f 1 is differentiable on (f(a), f(b))
More informationTS Method Summary. T k (x,y j 1 ) f(x j 1,y j 1 )+ 2 f (x j 1,y j 1 ) + k 1
TS Method Summary Let T k (x,y j 1 ) denote the first k +1 terms of the Taylor series expanded about the discrete approximation, (x j 1,y j 1 ), and ẑ k,j (x) be the polynomial approximation (to y(x))
More informationSECTION 5: POWER FLOW. ESE 470 Energy Distribution Systems
SECTION 5: POWER FLOW ESE 470 Energy Distribution Systems 2 Introduction Nodal Analysis 3 Consider the following circuit Three voltage sources VV sss, VV sss, VV sss Generic branch impedances Could be
More information8.3 Partial Fraction Decomposition
8.3 partial fraction decomposition 575 8.3 Partial Fraction Decomposition Rational functions (polynomials divided by polynomials) and their integrals play important roles in mathematics and applications,
More informationA Fast Surface Method to Model Skin Effect in Transmission Lines with Conductors of Arbitrary Shape or Rough Profile
A Fast Surface Method to Model Skin Effect in Transmission Lines with Conductors of Arbitrary Shape or Rough Profile Utkarsh R. Patel and Piero Triverio Edward S. Rogers Sr. Department of Electrical and
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
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 informationSCORE BOOSTER JAMB PREPARATION SERIES II
BOOST YOUR JAMB SCORE WITH PAST Polynomials QUESTIONS Part II ALGEBRA by H. O. Aliu J. K. Adewole, PhD (Editor) 1) If 9x 2 + 6xy + 4y 2 is a factor of 27x 3 8y 3, find the other factor. (UTME 2014) 3x
More informationHow to Explain Usefulness of Different Results when Teaching Calculus: Example of the Mean Value Theorem
Journal of Uncertain Systems Vol.7, No.3, pp.164-175, 2013 Online at: www.jus.org.uk How to Explain Usefulness of Different Results when Teaching Calculus: Example of the Mean Value Theorem Olga Kosheleva
More informationCSE 190: Reinforcement Learning: An Introduction. Chapter 8: Generalization and Function Approximation. Pop Quiz: What Function Are We Approximating?
CSE 190: Reinforcement Learning: An Introduction Chapter 8: Generalization and Function Approximation Objectives of this chapter: Look at how experience with a limited part of the state set be used to
More informationFast Buffer Insertion Considering Process Variation
Fast Buffer Insertion Considering Process Variation Jinjun Xiong, Lei He EE Department University of California, Los Angeles Sponsors: NSF, UC MICRO, Actel, Mindspeed Agenda Introduction and motivation
More informationpercentage of problems with ( 1 lb/ub ) <= x percentage of problems with ( 1 lb/ub ) <= x n= n=8 n= n=32 n= log10( x )
Soft vs. Hard Bounds in Probabilistic Robustness Analysis Xiaoyun Zhu Yun Huang John Doyle California Institute of Technology, Pasadena, CA 925 Abstract The relationship between soft vs. hard bounds and
More information