Differential-algebraic equations, the dummy derivatives method, and AD
|
|
- Conrad Tyler
- 5 years ago
- Views:
Transcription
1 Differential-algebraic equations, the dummy derivatives method, and AD John D Pryce Nedialko S Nedialkov 7th International Conference on Algorithmic Differentiation Oxford, UK, September 2016 September 14, /21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
2 1 DAEs & why they aren t ODEs intro 2 Structural analysis of DAEs geeky 3 The dummy derivatives (DDs) method for DAEs geekier 4 How AD helps apply the DDs method geekiest 5 Broader view: the MANDAE project strategic 2/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
3 DAEs and how to represent them A DAE for present purposes is a set of n equations relating n variables x j = x j (t) and some of their time-derivatives. I assume we know DAEs are important (Modelica, etc). As with ODEs, one can always reduce to equivalent first order DAE F(t, x, ẋ) = 0 the form accepted by many codes for the DAE initial value problem, such as DASSL, SUNDIALS,... But we prefer a form with arbitrary higher derivatives: f i ( t, the x j and derivatives of them ) = 0, i = 1,..., n. This lets one formulate many problems more concisely, hence solve by the Nedialkov-Pryce IVP code DAETS more conveniently. 3/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
4 Why DAEs aren t ODEs This system is a DAE: x 1 u(t) = 0, x 1 ẋ 2 = 0 where u(t) is a given driving function. Solving it requires integrating the driving function: x 2 = u(t)dt + C where C is a constant. It is an ODE in disguise. But also a DAE is: x 1 u(t) = 0, ẋ 1 x 2 = 0, whose unique solution has x 2 = u(t). Solving requires differentiating the driving function. It is not like an ODE at all. 4/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
5 Initial values are not obvious for DAEs Simple pendulum: state variables x(t), y(t), λ(t); constants G, L. 0 = A = ẍ + xλ 0 = B = ÿ + yλ G 0 = C = x 2 + y 2 L 2, Initial (x, y) must satisfy overt constraint C = 0 of lying on circle. Initial (ẋ, ẏ) must satisfy hidden constraint C = 0, that is 2(xẋ + yẏ) = 0: motion must be tangential to circle. L (x',y') One can describe configuration at time t by a point in 4D (x, y, ẋ, ẏ)-space that must obey the 2 equations C = 0, C = 0 defining the manifold M of consistent points at a given t. M s dimension, here 4 2 = 2, equals dof, the number of degrees of freedom how many independent IVs one can prescribe. 5/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016 (x,y)
6 Index Integer index (various definitions exist) attempts to measure how far DAE is from being an ODE. Probably most used is differentiation index ν d = the number of times it is necessary to differentiate one or more of the f i = 0 (w.r.t. t) such that from the result one can derive (*) an ODE in all the components. For the examples above x 1 u(t) = 0, x 1 ẋ 2 = 0 has ν d = 1; ODE could be ẋ 1 = u(t), ẋ 2 = u(t). x 1 u(t) = 0, ẋ 1 x 2 = 0 has ν d = 2; ODE could be ẋ 1 = u(t), ẋ 2 = ü(t). Defect: (*) makes you differentiate too much! To solve these, one only need differentiate 0 and 1 times respectively, not 1 and 2. 6/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
7 Index reduction, structural analysis Many numerical methods for higher-index DAEs start with index reduction: augment DAE by derivatives of some of its equations to produce a DAE of smaller index, but larger size. Differentiation index ν d is of limited use for finding how to do this. To compute ν d from its formal definition entails heavy linear algebra. So index reduction is usually done by discrete structural analysis (SA) algorithms based on DAE s sparsity structure. The classical SA algorithm is in the 1988 Pantelides paper The consistent initialisation of differential-algebraic systems. 7/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
8 Offsets Essential idea: find a number c i of times to differentiate equation f i = 0, that gives structurally nonsingular (SNS) system for resulting highest, d j th, derivatives of the x j. That is, one can renumber equations (or variables) so that x (d i ) i occurs in f (c i ) i for i = 1 : n. Equivalently the n n signature matrix Σ = (σ ij ) can be permuted to put finite values on the diagonal, where { highest order of derivative of xj in f i ; or σ ij = if x j does not occur in f i. Such offset vectors c = (c 1,..., c n ), d = (d 1,..., d n ) always exist unless DAE is structurally ill-posed, usually due to modeling error. Unique elementwise smallest non-negative c, d exist, the canonical offsets, which we assume chosen henceforth. 8/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
9 Jacobian An SA-friendly DAE is one whose n n system Jacobian ( ) J = f (c i ) i / x (d j ) j i,j=1 : n is nonsingular at some point on consistent manifold M (above). At such a point, if f i suitably smooth, solutions exist locally through this consistent point and any nearby consistent points. Pantelides method succeeds in consistent initialisation, e.g. for pendulum says initialise (x, y, ẋ, ẏ) DAE is SA-friendly dummy derivative method can be used Daets can solve the DAE. } In practice most DAEs are SA-friendly, don t ask me why. some aren t (Tischendorf, Tan), 9/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
10 Dummy derivatives, DDs Some index reduction methods convert DAE to ODE with more degrees of freedom than the DAE. Then the solution paths of the DAE are a proper subset of those of the ODE. Problem: numerical drift from manifold M, often exponential. DDs by Mattsson & Söderlind 1993 (M&S) are a systematic way to form equivalent ODE with exactly as many DOF as the DAE. View DAE as a flow on M. Then DDs describes this flow in a local coordinate system on M formed of some x j and derivatives. Numerical drift can only be within M, where it is less harmful. If solution path leaves patch of M where coordinates are nonsingular, one must choose new coordinates. This need for DD pivoting or switching complicates a numerical algorithm. 10/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
11 My take on DDs This description is equivalent to M&S s, but supports an implementation that makes DD switching almost trivial. (As M&S) Form derivatives of each f i up to the c i th, giving system f (l) i = 0, l = 0 : c i, i = 1 : n. Its unknowns are derivatives of the x j up to the d j th. View these as unrelated algebraic items x jl : x jl renames x (l) j, l = 0 : d j, j = 1 : n. It is under-determined in general: (no. of items) (no. of equations) = dof where dof is the number of degrees of freedom. 11/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
12 A way to formulate DDs, cont. DDs squares up the system by adding dof equations ẋ jl = x j,l+1, (j, l) in suitable set S with dof elements stating a genuine differential (not algebraic) relation for the selected x jl, which are state items, grouped into the current state vector x S. x S comprises the current local coordinate system of M. Result is an implicit ODE, solvable to give ODE of size dof, locally equivalent to the original DAE: ẋ S = F(t, x S ), I omit how state items are chosen linear algebra on columns of J. Key idea: any spanning set a basis (M&S) or dually: any lin indep set a basis (JDP). 12/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
13 Example: pendulum For pendulum the process produces 5 equations in 7 unknowns. The notation, here, is to rename x as x 0, ẋ as x 1, etc. etc. Augmented system 0 = A = ẍ + xλ 0 = B = ÿ + yλ G 0 = C = x 2 + y 2 L 2 0 = C = 2(xẋ + yẏ) 0 = C = 2(xẍ+ẋ 2 +yẏ+ẏ 2 ) After renaming 0 = A 0 = x 2 + x 0 λ 0 0 = B 0 = y 2 + y 0 λ 0 G 0 = C 0 = x y 2 0 L2 0 = C 1 = 2(x 0 x 1 + y 0 y 1 ) 0 = C 2 = 2(x 0 x 2 +x 2 1 +y 0y 1 +y 2 1 ) Here are the 7 unknowns showing possible state vectors. You can t avoid involving all 7 in solution process; this is general. Producible by DDs; sensible: 1st one says d dt [ ẋ x ] = F([ ẋ x ]) x ẋ ẍ y ẏ ÿ λ (compatible with AD) x ẋ ẍ y ẏ ÿ λ Not producible by DDs; silly: x ẋ ẍ y ẏ ÿ λ x ẋ ẍ y ẏ ÿ λ (poorly compatible with AD) 13/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
14 Strategic Q: how should AD tools help to automate numerical solution of SA-friendly DAEs? First, tool should support d/dt as a first-class operator also highly desirable for representing Modelica. At present widely used tools such as ADOL-C and dcc/dco do not have this ability. Daets uses Ole Stauning s C++ AD package Fadbad++. Stauning included Diff(, q) meaning d q /dt q, at our request. 14/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
15 Pendulum code E.g. code for pendulum, as in Daets user guide: 1 template <typename T> 2 void f c n (T t, const T z, T f, void param ) { 3 // z[0], z[1], z[2] are x, y, λ. 4 const double G = , L = ; 5 f [ 0 ] = D i f f ( z [ 0 ], 2 ) + z [ 0 ] z [ 2 ] ; //A = ẍ + xλ 6 f [ 1 ] = D i f f ( z [ 1 ], 2 ) + z [ 1 ] z [2] G ; //B = ÿ + yλ G 7 f [ 2 ] = s q r ( z [ 0 ] ) + s q r ( z [1]) s q r ( L ) ; //C = x 2 + y 2 L 2 8 } Daets instantiates T with 2 classes for structural analysis and 2 for numerical solution. 15/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
16 AD helping DAE solution (2) Second, tool must be able to differentiate f i selectively. E.g. for pendulum, leave A, B alone, differentiate C twice. Requires tool based on source code transformation? Actually no, as Ned observed. Don t treat different derivatives of one variable in isolation, but together as a truncated power series. Many tools do this naturally. E.g. pendulum really uses 3 objects x = (x 0, x 1, x 2 ) order 2 power series, y = (y 0, y 1, y 2 ) order 2 power series, λ = (λ 0 ) order 0 power series. Then normal Taylor series AD gives the necessary derivatives. 16/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
17 How this AD works (1): evaluate residuals Example: doing C = x 2 + y 2 L 2 proceeds via these steps: Input x = (x 0, x 1, x 2 ), also Lsq constant = (L 2, 0, 0) y = (y 0, y 1, y 2 ) Compute v 1 = x 2 = (x0 2, 2x 0x 1, 2x 0 x 2 + x1 2) v 2 = y 2 = (y0 2, 2y 0y 1, 2y 0 y 2 + y1 2) v 3 = v 1 + v 2 = (x0 2 + y 0 2, 2(x 0x 1 + y 0 y 1 ), 2(x 0 x 2 + y 0 y 2 ) + x1 2 + y 1 2) C = v 3 Lsq = (x0 2 + y 0 2 L2, 2(x 0 x 1 + y 0 y 1 ), 2(x 0 x 2 + y 0 y 2 ) + x1 2 + y 1 2) Returns order 2 power series C holding (C, C, C) = (C 0, C 1, C 2 ). Similarly A, B are returned as order 0 series A = (A 0 ) = (2x 2 + x 0 λ 0 ), B = (B 0 ) = (2y 2 + y 0 λ 0 ). (Second-order terms disagree by factor 2 with earlier since I ve gone from derivatives to Taylor coefficients.) 17/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
18 How this AD works (2): root-find In the context of reducing the DAE to an explicit ODE, suppose DDs has chosen state items x, ẋ, renamed as x 0, x 1. We implement the function F such that ẋ S = F(t, x S ), thus: 1. Input: values of t and x S = (x 0, x 1 ) items x F = (x 2, y 0, y 1, y 2, λ 0 ) are trial values that produce 5 residual values r = r(x F ) = (A 0, B 0, C 0, C 1, C 2 ). 3. Root-find using suitable Jacobians (more AD!) to get x F that makes r = 0; the F means found. 4. This finds x F as a function of t and x S (Implicit Function Theorem). 5. Extract x 2 from x F. Output: (x 1, x 2 ), which = (ẋ 0, ẋ 1 ) = ẋ S = F(t, x S ). A general way to prepare (nonstiff) DAE to solve by e.g. RK method. To DD-switch, change the index set S that defines x S, plus associated index sets and Jacobians. Ned confirms this is a viable way to implement DDs, using methods already provided by the Daets classes. 18/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
19 Broader view: the MANDAE project I m part of MANDAE bid to run EU Innovative Training Network (ITN). 14 ESRs (PhD students) at 7 partner sites across Europe, + McMaster. Technical mission: unify and extend tools for solving DAEs. Agreed: Use an Intermediate Representation IR an abstract object, mapped via an API to a program data structure that one can manipulate. IR = (IR-dae) + (IR-sa) = (DAE as syntax tree or similar) + (results of structural analysis). This can support both overloading and source-transformation tools. Needn t store differentiated equations explicitly as shown above. The IR has a meta aspect that complicates things: Daets uses a C++ form of the DAE, so it uses the C++ preprocessor, class structure, compiler & run-time system to get its effects. Our other code Daesa does similar in Matlab. The IR is passive, it cannot use these things the tools we build must deploy them on the IR s behalf. 19/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
20 Implications: Ned s basic architecture diagram Hence, to maximise use of existing code we propose this structure: Modelica function 1. M-parse IR-dae 2. CPPparse 3. unparse C++ function IR-dae + IR-sa IR-sa 5. generate F 6. generate f F(t,x,x') = 0 x' = f(t,x) 4. do SA 8. build F 7. build f code data structure information flow process DAETS "generate": implement by source-code generation "build": implement by operator overloading 20/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
21 Implications for MANDAE project Start project by specifying IR and its operations abstractly in principle but as a C++ class (IR-dae is recursive). Implement a C++ API of functions to create, write, read, manipulate IR objects. Done by key work-package WP1. Its ESR (= PhD student) will visit all project sites to gather requirements e.g. ensure API talks to AD tools smoothly. When diagram step 3 unparse done, existing Daets tools do step 4 do SA, achieving IR-dae IR-sa. Meantime other WPs work on step 1, Modelica IR-dae, and steps 4 8 which do IR simulation-code with various flavours of AD, linear algebra, exploiting DAE structure, etc. In particular a basic DDs implementation drops out, from Ned s previous work with Daets. WP1 is crucial, but when it s done coupling is low between other parts of project, reducing overall project risk. 21/21 Introduction Structural analysis Dummy derivatives AD and DDs Broader view: MANDAE AD Sep 2016
SOLVING DIFFERENTIAL-ALGEBRAIC EQUATIONS BY TAYLOR SERIES (II): COMPUTING THE SYSTEM JACOBIAN
BIT 0006-3835/00/4004-0001 $15.00 Submitted August 2005, pp. xxx xxx c Swets & Zealander SOLVING DIFFERENTIAL-ALGEBRAIC EQUATIONS BY TAYLOR SERIES II): COMPUTING THE SYSTEM JACOBIAN Nedialko S. Nedialkov
More informationThe DAETS Differential-Algebraic Equation Solver
The DAETS Differential-Algebraic Equation Solver John Pryce 1 Ned Nedialkov 2 1 Dept of Information Systems, Cranfield U., UK j.d.pryce@ntlworld.com 2 Dept of Computing and Software, McMaster U., Canada
More informationA DAESA a Matlab Tool for Structural Analysis of Differential-Algebraic Equations: Theory
A DAESA a Matlab Tool for Structural Analysis of Differential-Algebraic Equations: Theory JOHN D. PRYCE, Cardiff University NEDIALKO S. NEDIALKOV, McMaster University GUANGNING TAN, McMaster University
More informationSymbolic-Numeric Methods for Improving Structural Analysis of DAEs
Symbolic-Numeric Methods or Improving Structural Analysis o DAEs Guangning Tan, Nedialko S. Nedialkov, and John D. Pryce Abstract Systems o dierential-algebraic equations (DAEs) are generated routinely
More informationConversion Methods for Improving Structural. Analysis of Differential-Algebraic Equation Systems
Conversion Methods for Improving Structural Analysis of Differential-Algebraic Equation Systems CONVERSION METHODS FOR IMPROVING STRUCTURAL ANALYSIS OF DIFFERENTIAL-ALGEBRAIC EQUATION SYSTEMS BY GUANGNING
More informationSymbolic-Numeric Methods for Improving Structural Analysis of Differential-Algebraic Equation Systems
arxiv:1505.03445v1 [cs.sc] 13 May 2015 Symbolic-Numeric Methods for Improving Structural Analysis of Differential-Algebraic Equation Systems Guangning Tan, Nedialko S. Nedialkov McMaster University John
More informationIndex Reduction and Discontinuity Handling using Substitute Equations
Mathematical and Computer Modelling of Dynamical Systems, vol. 7, nr. 2, 2001, 173-187. Index Reduction and Discontinuity Handling using Substitute Equations G. Fábián, D.A. van Beek, J.E. Rooda Abstract
More informationThe Unifying View on Ordinary Differential Equations and Automatic Differentiation, yet with a Gap to Fill Alexander Gofen.
The Unifying View on Ordinary Differential Equations and Automatic Differentiation, yet with a Gap to Fill Alexander Gofen Conventions Differentiation is understood as for holomorphic functions in the
More information10.34: Numerical Methods Applied to Chemical Engineering. Lecture 19: Differential Algebraic Equations
10.34: Numerical Methods Applied to Chemical Engineering Lecture 19: Differential Algebraic Equations 1 Recap Differential algebraic equations Semi-explicit Fully implicit Simulation via backward difference
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationJordan normal form notes (version date: 11/21/07)
Jordan normal form notes (version date: /2/7) If A has an eigenbasis {u,, u n }, ie a basis made up of eigenvectors, so that Au j = λ j u j, then A is diagonal with respect to that basis To see this, let
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second
More informationCHAPTER 5: Linear Multistep Methods
CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order
More informationOrdinary Differential Equation Introduction and Preliminaries
Ordinary Differential Equation Introduction and Preliminaries There are many branches of science and engineering where differential equations arise naturally. Now days, it finds applications in many areas
More informationLinear Multistep Methods I: Adams and BDF Methods
Linear Multistep Methods I: Adams and BDF Methods Varun Shankar January 1, 016 1 Introduction In our review of 5610 material, we have discussed polynomial interpolation and its application to generating
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
More informationChapter 0 of Calculus ++, Differential calculus with several variables
Chapter of Calculus ++, Differential calculus with several variables Background material by Eric A Carlen Professor of Mathematics Georgia Tech Spring 6 c 6 by the author, all rights reserved - Table of
More informationCALCULUS ON MANIFOLDS
CALCULUS ON MANIFOLDS 1. Manifolds Morally, manifolds are topological spaces which locally look like open balls of the Euclidean space R n. One can construct them by piecing together such balls ( cells
More informationM. Matrices and Linear Algebra
M. Matrices and Linear Algebra. Matrix algebra. In section D we calculated the determinants of square arrays of numbers. Such arrays are important in mathematics and its applications; they are called matrices.
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationOrdinary Differential Equations (ODEs)
Ordinary Differential Equations (ODEs) 1 Computer Simulations Why is computation becoming so important in physics? One reason is that most of our analytical tools such as differential calculus are best
More informationFormal Verification and Automated Generation of Invariant Sets
Formal Verification and Automated Generation of Invariant Sets Khalil Ghorbal Carnegie Mellon University Joint work with Andrew Sogokon and André Platzer Toulouse, France 11-12 June, 2015 K. Ghorbal (CMU,
More informationIntroduction to Initial Value Problems
Chapter 2 Introduction to Initial Value Problems The purpose of this chapter is to study the simplest numerical methods for approximating the solution to a first order initial value problem (IVP). Because
More informationIntroduction to Vectors
Introduction to Vectors K. Behrend January 31, 008 Abstract An introduction to vectors in R and R 3. Lines and planes in R 3. Linear dependence. 1 Contents Introduction 3 1 Vectors 4 1.1 Plane vectors...............................
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationchapter 12 MORE MATRIX ALGEBRA 12.1 Systems of Linear Equations GOALS
chapter MORE MATRIX ALGEBRA GOALS In Chapter we studied matrix operations and the algebra of sets and logic. We also made note of the strong resemblance of matrix algebra to elementary algebra. The reader
More informationOutline Introduction: Problem Description Diculties Algebraic Structure: Algebraic Varieties Rank Decient Toeplitz Matrices Constructing Lower Rank St
Structured Lower Rank Approximation by Moody T. Chu (NCSU) joint with Robert E. Funderlic (NCSU) and Robert J. Plemmons (Wake Forest) March 5, 1998 Outline Introduction: Problem Description Diculties Algebraic
More informationA NEW SET THEORY FOR ANALYSIS
Article A NEW SET THEORY FOR ANALYSIS Juan Pablo Ramírez 0000-0002-4912-2952 Abstract: We present the real number system as a generalization of the natural numbers. First, we prove the co-finite topology,
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationLecture 2: Computing functions of dense matrices
Lecture 2: Computing functions of dense matrices Paola Boito and Federico Poloni Università di Pisa Pisa - Hokkaido - Roma2 Summer School Pisa, August 27 - September 8, 2018 Introduction In this lecture
More informationPractical Algebra. A Step-by-step Approach. Brought to you by Softmath, producers of Algebrator Software
Practical Algebra A Step-by-step Approach Brought to you by Softmath, producers of Algebrator Software 2 Algebra e-book Table of Contents Chapter 1 Algebraic expressions 5 1 Collecting... like terms 5
More informationKeble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II
Tomi Johnson 1 Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II Prepare full solutions to the problems with a self assessment of your progress on a cover
More informationComputing validated solutions of implicit differential equations
Advances in Computational Mathematics 19: 231 253, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. Computing validated solutions of implicit differential equations Jens Hoefkens a, Martin
More informationLinear Algebra and Robot Modeling
Linear Algebra and Robot Modeling Nathan Ratliff Abstract Linear algebra is fundamental to robot modeling, control, and optimization. This document reviews some of the basic kinematic equations and uses
More informationHow to Use Calculus Like a Physicist
How to Use Calculus Like a Physicist Physics A300 Fall 2004 The purpose of these notes is to make contact between the abstract descriptions you may have seen in your calculus classes and the applications
More informationALGEBRA. 1. Some elementary number theory 1.1. Primes and divisibility. We denote the collection of integers
ALGEBRA CHRISTIAN REMLING 1. Some elementary number theory 1.1. Primes and divisibility. We denote the collection of integers by Z = {..., 2, 1, 0, 1,...}. Given a, b Z, we write a b if b = ac for some
More information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More informationPorteous s Formula for Maps between Coherent Sheaves
Michigan Math. J. 52 (2004) Porteous s Formula for Maps between Coherent Sheaves Steven P. Diaz 1. Introduction Recall what the Thom Porteous formula for vector bundles tells us (see [2, Sec. 14.4] for
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More informationDifferential Equation Types. Moritz Diehl
Differential Equation Types Moritz Diehl Overview Ordinary Differential Equations (ODE) Differential Algebraic Equations (DAE) Partial Differential Equations (PDE) Delay Differential Equations (DDE) Ordinary
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationOptimality, Duality, Complementarity for Constrained Optimization
Optimality, Duality, Complementarity for Constrained Optimization Stephen Wright University of Wisconsin-Madison May 2014 Wright (UW-Madison) Optimality, Duality, Complementarity May 2014 1 / 41 Linear
More informationGAUSSIAN ELIMINATION AND LU DECOMPOSITION (SUPPLEMENT FOR MA511)
GAUSSIAN ELIMINATION AND LU DECOMPOSITION (SUPPLEMENT FOR MA511) D. ARAPURA Gaussian elimination is the go to method for all basic linear classes including this one. We go summarize the main ideas. 1.
More information4 Film Extension of the Dynamics: Slowness as Stability
4 Film Extension of the Dynamics: Slowness as Stability 4.1 Equation for the Film Motion One of the difficulties in the problem of reducing the description is caused by the fact that there exists no commonly
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationONE-STEP 4-STAGE HERMITE-BIRKHOFF-TAYLOR DAE SOLVER OF ORDER 12
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 16, Number 4, Winter 2008 ONE-STEP 4-STAGE HERMITE-BIRKHOFF-TAYLOR DAE SOLVER OF ORDER 12 TRUONG NGUYEN-BA, HAN HAO, HEMZA YAGOUB AND RÉMI VAILLANCOURT ABSTRACT.
More informationA primer on matrices
A primer on matrices Stephen Boyd August 4, 2007 These notes describe the notation of matrices, the mechanics of matrix manipulation, and how to use matrices to formulate and solve sets of simultaneous
More informationInterpolation and Approximation
Interpolation and Approximation The Basic Problem: Approximate a continuous function f(x), by a polynomial p(x), over [a, b]. f(x) may only be known in tabular form. f(x) may be expensive to compute. Definition:
More informationThe Jacobian. Jesse van den Kieboom
The Jacobian Jesse van den Kieboom jesse.vandenkieboom@epfl.ch 1 Introduction 1 1 Introduction The Jacobian is an important concept in robotics. Although the general concept of the Jacobian in robotics
More informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationENGINEERING MATH 1 Fall 2009 VECTOR SPACES
ENGINEERING MATH 1 Fall 2009 VECTOR SPACES A vector space, more specifically, a real vector space (as opposed to a complex one or some even stranger ones) is any set that is closed under an operation of
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 information18.02 Multivariable Calculus Fall 2007
MIT OpenCourseWare http://ocw.mit.edu 18.02 Multivariable Calculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. M. Matrices and Linear Algebra
More informationManifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds
MA 755 Fall 05. Notes #1. I. Kogan. Manifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds Definition 1 An n-dimensional C k -differentiable manifold
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 informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationMECH 576 Geometry in Mechanics September 16, 2009 Using Line Geometry
MECH 576 Geometry in Mechanics September 16, 2009 Using Line Geometry 1 Deriving Equations in Line Coordinates Four exercises in deriving fundamental geometric equations with line coordinates will be conducted.
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationES.1803 Topic 13 Notes Jeremy Orloff
ES.1803 Topic 13 Notes Jeremy Orloff 13 Vector Spaces, matrices and linearity 13.1 Goals 1. Know the definition of a vector space and how to show that a given set is a vector space. 2. Know the meaning
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationOrdinary Differential Equations (ODEs)
c01.tex 8/10/2010 22: 55 Page 1 PART A Ordinary Differential Equations (ODEs) Chap. 1 First-Order ODEs Sec. 1.1 Basic Concepts. Modeling To get a good start into this chapter and this section, quickly
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space 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 An inner product
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More information10 Transfer Matrix Models
MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationMath 3361-Modern Algebra Lecture 08 9/26/ Cardinality
Math 336-Modern Algebra Lecture 08 9/26/4. Cardinality I started talking about cardinality last time, and you did some stuff with it in the Homework, so let s continue. I said that two sets have the same
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More informationLinear Algebra: A Constructive Approach
Chapter 2 Linear Algebra: A Constructive Approach In Section 14 we sketched a geometric interpretation of the simplex method In this chapter, we describe the basis of an algebraic interpretation that allows
More informationRunge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot. Department U2IS ENSTA ParisTech SCAN Uppsala
Runge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot Department U2IS ENSTA ParisTech SCAN 2016 - Uppsala Contents Numerical integration Runge-Kutta with interval
More information7 Differential Algebraic Equations
7 Differential Algebraic Equations Preview In this chapter, we shall analyze simulation problems that don t present themselves initially in an explicit state space form. For many physical systems, it is
More informationDefinition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.
5 Vector fields Last updated: March 12, 2012. 5.1 Definition and general properties We first need to define what a vector field is. Definition 5.1. A vector field v on a manifold M is map M T M such that
More informationEco517 Fall 2004 C. Sims MIDTERM EXAM
Eco517 Fall 2004 C. Sims MIDTERM EXAM Answer all four questions. Each is worth 23 points. Do not devote disproportionate time to any one question unless you have answered all the others. (1) We are considering
More informationAIBO experiences change of surface incline.
NW Computational Intelligence Laboratory AIBO experiences change of surface incline. Decisions based only on kinesthetic experience vector ICNC 07, China 8/27/07 # 51 NW Computational Intelligence Laboratory
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationReview for the First Midterm Exam
Review for the First Midterm Exam Thomas Morrell 5 pm, Sunday, 4 April 9 B9 Van Vleck Hall For the purpose of creating questions for this review session, I did not make an effort to make any of the numbers
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationDelay Differential Equations with Constant Lags
Delay Differential Equations with Constant Lags L.F. Shampine Mathematics Department Southern Methodist University Dallas, TX 75275 shampine@smu.edu S. Thompson Department of Mathematics & Statistics Radford
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationExtreme Values and Positive/ Negative Definite Matrix Conditions
Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1
More informationImplicit Functions, Curves and Surfaces
Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then
More informationLinear Algebra, Summer 2011, pt. 2
Linear Algebra, Summer 2, pt. 2 June 8, 2 Contents Inverses. 2 Vector Spaces. 3 2. Examples of vector spaces..................... 3 2.2 The column space......................... 6 2.3 The null space...........................
More informationLagrangian for Central Potentials
Physics 411 Lecture 2 Lagrangian for Central Potentials Lecture 2 Physics 411 Classical Mechanics II August 29th 2007 Here we will review the Lagrange formulation in preparation for the study of the central
More informationIN this paper, we consider the capacity of sticky channels, a
72 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 54, NO. 1, JANUARY 2008 Capacity Bounds for Sticky Channels Michael Mitzenmacher, Member, IEEE Abstract The capacity of sticky channels, a subclass of insertion
More informationVector Spaces. 9.1 Opening Remarks. Week Solvable or not solvable, that s the question. View at edx. Consider the picture
Week9 Vector Spaces 9. Opening Remarks 9.. Solvable or not solvable, that s the question Consider the picture (,) (,) p(χ) = γ + γ χ + γ χ (, ) depicting three points in R and a quadratic polynomial (polynomial
More informationIntroduction to Group Theory
Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)
More informationNumerical stroboscopic averaging for ODEs and DAEs
Numerical stroboscopic averaging for ODEs and DAEs M. P. Calvo Universidad de Valladolid, Spain Joint work with Ph. Chartier, A. Murua, J. M. Sanz-Serna 1 I. HIGHLY OSCILLATORY PROBLEMS 2 Consider the
More informationCHAPTER 2. Techniques for Solving. Second Order Linear. Homogeneous ODE s
A SERIES OF CLASS NOTES FOR 005-006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES A COLLECTION OF HANDOUTS ON SCALAR LINEAR ORDINARY DIFFERENTIAL
More informationLeast-squares Solutions of Linear Differential Equations
1 Least-squares Solutions of Linear Differential Equations Daniele Mortari dedicated to John Lee Junkins arxiv:1700837v1 [mathca] 5 Feb 017 Abstract This stu shows how to obtain least-squares solutions
More informationMath Camp II. Basic Linear Algebra. Yiqing Xu. Aug 26, 2014 MIT
Math Camp II Basic Linear Algebra Yiqing Xu MIT Aug 26, 2014 1 Solving Systems of Linear Equations 2 Vectors and Vector Spaces 3 Matrices 4 Least Squares Systems of Linear Equations Definition A linear
More informationA relaxation of the strangeness index
echnical report from Automatic Control at Linköpings universitet A relaxation of the strangeness index Henrik idefelt, orkel Glad Division of Automatic Control E-mail: tidefelt@isy.liu.se, torkel@isy.liu.se
More informationf(x 0 + h) f(x 0 ) h slope of secant line = m sec
Derivatives Using limits, we can define the slope of a tangent line to a function. When given a function f(x), and given a point P (x 0, f(x 0 )) on f, if we want to find the slope of the tangent line
More information1.3 Linear Dependence & span K
( ) Conversely, suppose that every vector v V can be expressed uniquely as v = u +w for u U and w W. Then, the existence of this expression for each v V is simply the statement that V = U +W. Moreover,
More informationDefinition: An n x n matrix, "A", is said to be diagonalizable if there exists a nonsingular matrix "X" and a diagonal matrix "D" such that X 1 A X
DIGONLIZTION Definition: n n x n matrix, "", is said to be diagonalizable if there exists a nonsingular matrix "X" and a diagonal matrix "D" such that X X D. Theorem: n n x n matrix, "", is diagonalizable
More informationVerification of analog and mixed-signal circuits using hybrid systems techniques
FMCAD, November 2004, Austin Verification of analog and mixed-signal circuits using hybrid systems techniques Thao Dang, Alexandre Donze, Oded Maler VERIMAG Grenoble, France Plan 1. Introduction 2. Verification
More informationLinear Regression 9/23/17. Simple linear regression. Advertising sales: Variance changes based on # of TVs. Advertising sales: Normal error?
Simple linear regression Linear Regression Nicole Beckage y " = β % + β ' x " + ε so y* " = β+ % + β+ ' x " Method to assess and evaluate the correlation between two (continuous) variables. The slope of
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 information