Exponential integrators
|
|
- Ada Stephens
- 5 years ago
- Views:
Transcription
1 Exponential integrators Marlis Hochbruck Heinrich-Heine University Düsseldorf, Germany Alexander Ostermann University Innsbruck, Austria Helsinki, May 25 p.1
2 Outline Time dependent Schrödinger equations application: quantum dynamics methods: Magnus integrators 2nd order differential equations application: laser/plasma interaction nonlinear wave equations methods: Gautschi-type integrators p.2
3 Quantum dynamics simulations Schrödinger equation: ψ (t) = ih(t)ψ(t), H(t) = U + V (t), U x (hence H(t) large, 1/ x 2 ) standard scheme: split-step Fourier method p.3
4 Quantum dynamics simulations Schrödinger equation: ψ (t) = ih(t)ψ(t), H(t) = U + V (t), U x (hence H(t) large, 1/ x 2 ) standard scheme: split-step Fourier method alternative: exponential midpoint rule (Magnus method) ψ n+1 = exp ( ihh(t n+1/2 ) ) ψ n (exact for constant H) p.3
5 Magnus integrators y = A(t)y(t), y() = y Magnus, 1954: determine Ω(t) such that y(t) = exp(ω(t))y Magnus expansion (valid for Ω(t) < π) Ω(t) = t A(τ)dτ t t t τ [ [ τ [ τ A(σ)dσ, A(τ)] dτ σ [ A(µ)dµ, A(σ)] dσ, A(τ)] dτ A(σ)dσ, [ τ A(µ)dµ, A(τ)]] dτ +... p.4
6 Numerical methods Magnus integrators y n+1 = exp(ω n )y n, Ω n Ω(h) (review: Iserles, Munthe-Kaas, Nørsett, Zanna, ) Approximation involves truncating the Magnus expansion (after k terms) t k = 1 : k = 2 : Ω(t) Ω(t) t A(t n + τ)dτ A(t n + τ)dτ 1 2 t τ [ A(t n + σ)dσ, A(t n + τ)] dτ approximating integrals by replacing A(t) by interpolation polynomial Â(t) for quadrature nodes t n + c j h p.5
7 Examples of Magnus integrators k = 1, exponential midpoint rule Ω n = ha(t n + h/2). k = 2, two-point Gauß quadrature rule: Ω n = h 2 (A 1 + A 2 ) + 3h 2 12 [A 2, A 1 ], A j = A(t n + c j h), c j nodes of Gauß quadrature rule k = 2, method by Blanes, Casas, Ros Ω n = h 6 ( A(tn )+4A(t n+1/2 )+A(t n+1 ) ) h2 12 [A(t n), A(t n+1 )]. p.6
8 Convergence for Schrödinger equations Theorem (H., Lubich, 23) Exponential midpoint rule y n y(t n ) Ch 2 t n 4th-order Gauß method: for h D c y n y(t n ) Ch 4 t n max Dy(t) t t n max D 3 y(t) t t n p.7
9 Convergence for Schrödinger equations Theorem (H., Lubich, 23) Exponential midpoint rule y n y(t n ) Ch 2 t n 4th-order Gauß method: for h D c y n y(t n ) Ch 4 t n max Dy(t) t t n max D 3 y(t) t t n error bound for classical implicit midpoint rule: y n y(t n ) Ch 2 t n max d3 t t n dt 3 y(t) p.7
10 Numerical experiment i ψ t = 1 2 ψ + b(x, t)ψ, x = (x 1,..., x d ) R d, t > discretized s. t. h D 3.5 (N = 32,..., 248 Fourier modes) midpoint Gauss 4 Blanes et al time steps versus error p.8
11 2nd order differential equations y = f(y), y() = y, y () = y standard scheme: leap-frog, Störmer-Verlet y n+1 2y n + y n 1 = h 2 f(y n ) nice properties: order two, symmetric, symplectic Problem: stability requires hω max 2 (1/2) derive alternative without this restriction for y + Ω 2 (t, y)y = g(y) which reduces to leap-frog for Ω = p.9
12 Application: Laser plasma interaction a A =.1, Q =.3 density profile pulses at different times z/λ a A =.12, Q =.3 density profile pulses at different times z/λ p.1
13 Physical model a vectorpotential, E parallel component of electrical field 2 E t 2 2 a z 2 2 a t 2 = Q (n + δn) γ + Q n E = n γ z a where n plasma profile and γ 2 = 1 + a 2 δn = E z relativistic factor density variation p.11
14 2nd order differential equations y + Ω 2 (t, y)y = g(y), y() = y, y () = y for constant g and Ω, exact solution satisfies y(t+h) 2y(t)+y(t h) = h 2 ψ(hω) ( g Ω 2 y(t) ), ψ(x) = sinc 2 x 2 (variation-of-constants formula) Gautschi-type exponential integrator y n+1 2y n + y n 1 = h 2 ψ(hω n ) ( g n Ω 2 ) n y n, Ωn = Ω(t n, y n ) choice of g n? p.12
15 2nd order differential equations y n+1 2y n + y n 1 = h 2 ψ(hω n ) ( g n Ω 2 n y n ) obvious choice: g n = g(y n ) Gautschi 61 resonance problems for hω k jπ, ω k eigenvalue of Ω n p.13
16 2nd order differential equations y n+1 2y n + y n 1 = h 2 ψ(hω n ) ( g n Ω 2 n y n ) obvious choice: g n = g(y n ) Gautschi 61 resonance problems for hω k jπ, ω k eigenvalue of Ω n better choice: g n = g(φ(hω n )y n ), φ filter function φ() = 1, φ(kπ) =, k = 1, 2, 3,... p.13
17 2nd order differential equations y n+1 2y n + y n 1 = h 2 ψ(hω n ) ( g n Ω 2 n y n ) obvious choice: g n = g(y n ) Gautschi 61 resonance problems for hω k jπ, ω k eigenvalue of Ω n better choice: g n = g(φ(hω n )y n ), φ filter function φ() = 1, φ(kπ) =, k = 1, 2, 3,... convergence result: (H., Lubich, 99, Grimm 2, 5) Assumptions: g smooth, bounded energy: y n y(t n ) h 2 C(t n ), C(t n ) e t nl p.13
18 Effect of filter function 1 1 without φ with φ p.14
19 Generalization: multiple time-stepping y = f(y) + g(y) given y n and y n compute suitable averaged value y n and solution of u = f(u) + g(y n ), u() = y n, u () = y n. (e.g. by numerical methods with smaller time steps) compute y n+1 and y n+1 from y n+1 2y n + y n 1 = u(h) 2u() + u( h), y n+1 y n 1 = u (h) u ( h), symmetric multiple-time-stepping scheme p.15
20 Results for laser / plasma simulation maximum amplitude error of 1% (phase error ignored) A =.12, Q =.3 dz dt runtime save leap frog leap frog + quasi-envelope % Gautschi % Gautschi + quasi-envelope % PIC multilevel approach not included! 1 day Karle, Schweitzer, H., Laedke, Spatschek, Preprint 25 p.16
21 Quasi-envelope approach x 1 3 A =.12, Q =.3, spectrum k/k replace vectorpotential a by a(z, t) = ã(z, t)e iκz p.17
22 Phase errors at t = t end intensity A =.12, Q =.3 ref PIC KG all tricks z/λ p.18
23 Possible topics for informal sessions Basics on numerical time integration Implementation of exponential integrators (Krylov subspace methods) Numerical methods for time-dependent Schrödinger equations (Magnus integrators) Numerical methods for nonlinear wave equations (Laser / plasma interaction) Numerical methods for parabolic problems p.19
Exponential integrators for oscillatory second-order differential equations
Exponential integrators for oscillatory second-order differential equations Marlis Hochbruck and Volker Grimm Mathematisches Institut Heinrich Heine Universität Düsseldorf Cambridge, March 2007 Outline
More informationNEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING Fernando Casas Universitat Jaume I, Castellón, Spain Fernando.Casas@uji.es (on sabbatical leave at DAMTP, University of Cambridge) Edinburgh,
More informationDynamical Low-Rank Approximation to the Solution of Wave Equations
Dynamical Low-Rank Approximation to the Solution of Wave Equations Julia Schweitzer joint work with Marlis Hochbruck INSTITUT FÜR ANGEWANDTE UND NUMERISCHE MATHEMATIK 1 KIT Universität des Landes Baden-Württemberg
More informationExponential multistep methods of Adams-type
Exponential multistep methods of Adams-type Marlis Hochbruck and Alexander Ostermann KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 KIT University of the State of Baden-Wuerttemberg and National Laboratory
More informationA parallel implementation of a two-dimensional fluid laser-plasma integrator for stratified plasma-vacuum systems
A parallel implementation of a two-dimensional fluid laser-plasma integrator for stratified plasma-vacuum systems Ch. Karle a, J. Schweitzer b, M. Hochbruck b, K.H. Spatschek a, a Institut für Theoretische
More informationTHE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS
THE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS MARIANNA KHANAMIRYAN Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road,
More informationFavourable time integration methods for. Non-autonomous evolution equations
Favourable time integration methods for non-autonomous evolution equations International Conference on Scientific Computation and Differential Equations (SciCADE) Session on Numerical integration of non-autonomous
More informationExponential Runge-Kutta methods for parabolic problems
Exponential Runge-Kutta methods for parabolic problems Marlis Hochbruck a,, Alexander Ostermann b a Mathematisches Institut, Heinrich-Heine Universität Düsseldorf, Universitätsstraße 1, D-4225 Düsseldorf,
More informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More information13.1 Ion Acoustic Soliton and Shock Wave
13 Nonlinear Waves In linear theory, the wave amplitude is assumed to be sufficiently small to ignore contributions of terms of second order and higher (ie, nonlinear terms) in wave amplitude In such a
More informationCOMPOSITION MAGNUS INTEGRATORS. Sergio Blanes. Lie Group Methods and Control Theory Edinburgh 28 June 1 July Joint work with P.C.
COMPOSITION MAGNUS INTEGRATORS Sergio Blanes Departament de Matemàtiques, Universitat Jaume I Castellón) sblanes@mat.uji.es Lie Group Methods and Control Theory Edinburgh 28 June 1 July 24 Joint work with
More informationHigh-order actions and their applications to honor our friend and
High-order actions and their applications to honor our friend and collaborator Siu A. Chin Universitat Politecnica de Catalunya, Barcelona, Spain Lecture: Magnus Expansion and Suzuki s Method Outline of
More informationMulti-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems
Multi-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems Jürgen Geiser and Siu A. Chin geiser@mathematik.hu-berlin.de Abstract. In this paper we discuss the
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationComparative numerical solutions of stiff ordinary differential equations using Magnus series expansion method
NTMSCI, No. 1, 5-45 (215) 5 New Trends in Mathematical Sciences http://www.ntmsci.com Comparative numerical solutions of stiff ordinary differential equations using Magnus series expansion method Mehmet
More information5.1 Classical Harmonic Oscillator
Chapter 5 Harmonic Oscillator 5.1 Classical Harmonic Oscillator m l o l Hooke s Law give the force exerting on the mass as: f = k(l l o ) where l o is the equilibrium length of the spring and k is the
More informationThree approaches for the design of adaptive time-splitting methods
Three approaches for the design of adaptive time-splitting methods Mechthild Thalhammer Leopold Franzens Universität Innsbruck, Austria Workshop on Geometric Integration and Computational Mechanics Organisers:
More informationExponential integrators
Acta Numerica (2), pp. 29 286 c Cambridge University Press, 2 doi:.7/s962492948 Printed in the United Kingdom Exponential integrators Marlis Hochbruck Karlsruher Institut für Technologie, Institut für
More informationS. Sánchez, F. Casas & A. Fernández
New analytic approximations based on the Magnus expansion S. Sánchez, F. Casas & A. Fernández Journal of Mathematical Chemistry ISSN 259-9791 Volume 49 Number 8 J Math Chem (211) 49:1741-1758 DOI 1.17/s191-11-9855-
More informationIntegration of Vlasov-type equations
Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma
More informationEXPONENTIAL ROSENBROCK-TYPE METHODS
EXPONENTIAL ROSENBROCK-TYPE METHODS MARLIS HOCHBRUCK, ALEXANDER OSTERMANN, AND JULIA SCHWEITZER Abstract. We introduce a new class of exponential integrators for the numerical integration of large-scale
More informationTime dependent perturbation theory 1 D. E. Soper 2 University of Oregon 11 May 2012
Time dependent perturbation theory D. E. Soper University of Oregon May 0 offer here some background for Chapter 5 of J. J. Sakurai, Modern Quantum Mechanics. The problem Let the hamiltonian for a system
More informationExponential integrators and functions of the matrix exponential
Exponential integrators and functions of the matrix exponential Paul Matthews, Stephen Cox, Hala Ashi and Linda Cummings School of Mathematical Sciences, University of Nottingham, UK Introduction to exponential
More informationSolutions to Problems in Chapter 4
Solutions to Problems in Chapter 4 Problems with Solutions Problem 4. Fourier Series of the Output Voltage of an Ideal Full-Wave Diode Bridge Rectifier he nonlinear circuit in Figure 4. is a full-wave
More informationAnton ARNOLD. with N. Ben Abdallah (Toulouse), J. Geier (Vienna), C. Negulescu (Marseille) TU Vienna Institute for Analysis and Scientific Computing
TECHNISCHE UNIVERSITÄT WIEN Asymptotically correct finite difference schemes for highly oscillatory ODEs Anton ARNOLD with N. Ben Abdallah (Toulouse, J. Geier (Vienna, C. Negulescu (Marseille TU Vienna
More informationCENTER MANIFOLD AND NORMAL FORM THEORIES
3 rd Sperlonga Summer School on Mechanics and Engineering Sciences 3-7 September 013 SPERLONGA CENTER MANIFOLD AND NORMAL FORM THEORIES ANGELO LUONGO 1 THE CENTER MANIFOLD METHOD Existence of an invariant
More informationApplications of Lie group integrators and exponential schemes
Applications of Lie group integrators and exponential schemes Workshop on Lie group methods and control theory, Edinburgh, 28.6-1.7, 2004, Brynjulf Owren Dept of Math Sci, NTNU Applications of Lie group
More informationHigh-order time-splitting methods for Schrödinger equations
High-order time-splitting methods for Schrödinger equations and M. Thalhammer Department of Mathematics, University of Innsbruck Conference on Scientific Computing 2009 Geneva, Switzerland Nonlinear Schrödinger
More informationReview of Frequency Domain Fourier Series: Continuous periodic frequency components
Today we will review: Review of Frequency Domain Fourier series why we use it trig form & exponential form how to get coefficients for each form Eigenfunctions what they are how they relate to LTI systems
More informationEfficient methods for time-dependence in semiclassical Schrödinger equations
Efficient methods for time-dependence in semiclassical Schrödinger equations Philipp Bader, Arieh Iserles, Karolina Kropielnicka & Pranav Singh November 25, 25 Abstract We build an efficient and unitary
More informationHomework 4. May An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt
Homework 4 May 2017 1. An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt Determine the impulse response of the system. Rewriting as y(t) = t e (t
More informationImplicit-explicit exponential integrators
Implicit-explicit exponential integrators Bawfeh Kingsley Kometa joint work with Elena Celledoni MaGIC 2011 Finse, March 1-4 1 Introduction Motivation 2 semi-lagrangian Runge-Kutta exponential integrators
More informationEE 224 Signals and Systems I Review 1/10
EE 224 Signals and Systems I Review 1/10 Class Contents Signals and Systems Continuous-Time and Discrete-Time Time-Domain and Frequency Domain (all these dimensions are tightly coupled) SIGNALS SYSTEMS
More informationPhysics Reports 470 (2009) Contents lists available at ScienceDirect. Physics Reports. journal homepage:
Physics Reports 47 9 151 38 Contents lists available at ScienceDirect Physics Reports journal homepage: www.elsevier.com/locate/physrep The Magnus expansion and some of its applications S. Blanes a, F.
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationIntroduction to Hamiltonian systems
Introduction to Hamiltonian systems Marlis Hochbruck Heinrich-Heine Universität Düsseldorf Oberwolfach Seminar, November 2008 Examples Mathematical biology: Lotka-Volterra model First numerical methods
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationA pedagogical approach to Magnus expansion
A pedagogical approach to Magnus expansion S Blanes F Casas J A Oteo and J. Ros Instituto de Matemática Multidisciplinar, Universidad Politécnica de Valencia, E-46022 Valencia, Spain Departament de Matemàtiques,
More informationSchrödinger equation
Splitting methods for the time dependent Schrödinger equation S. Blanes, Joint work with F. Casas and A. Murua Instituto de Matemática Multidisciplinar Universidad Politécnica de Valencia, SPAIN Workshop:
More informationarxiv: v1 [math.na] 6 Nov 2017
Efficient boundary corrected Strang splitting Lukas Einkemmer Martina Moccaldi Alexander Ostermann arxiv:1711.02193v1 [math.na] 6 Nov 2017 Version of November 6, 2017 Abstract Strang splitting is a well
More informationA stroboscopic averaging technique for highly-oscillatory Schrödinger equation
A stroboscopic averaging technique for highly-oscillatory Schrödinger equation F. Castella IRMAR and INRIA-Rennes Joint work with P. Chartier, F. Méhats and A. Murua Saint-Malo Workshop, January 26-28,
More informationPHYS-454 The position and momentum representations
PHYS-454 The position and momentum representations 1 Τhe continuous spectrum-a n So far we have seen problems where the involved operators have a discrete spectrum of eigenfunctions and eigenvalues.! n
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationLecture 4: High oscillations
Geometric Numerical Integration TU München Ernst Hairer January February Lecture 4: High oscillations Table of contents A Fermi Pasta Ulam type problem Application of classical integrators 4 3 Exponential
More informationApproximating the matrix exponential of an advection-diffusion operator using the incomplete orthogonalization method
Approximating the matrix exponential of an advection-diffusion operator using the incomplete orthogonalization method Antti Koskela KTH Royal Institute of Technology, Lindstedtvägen 25, 10044 Stockholm,
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationPhysics in the lane. For Ladislav, Bergen December 19th
Physics in the lane For Ladislav, Bergen December 19th A laser pulse and a hydrogen atom ee 0 v mω The story dates back PRL 95, 043601 (005) PRL 97, 043601 (006) i ℏ Ψ =[c α ( p+e A)+V +m c β] Ψ i ℏ Ψ
More informationCommutator-free Magnus Lanczos methods for the linear Schrödinger equation
Commutator-free Magnus Lanczos methods for the linear Schrödinger equation Arieh Iserles, Karolina Kropielnicka & Pranav Singh October 3, 27 AMS Mathematics Subject Classification: Primary 65M7, Secondary
More informationSplittings, commutators and the linear Schrödinger equation
UNIVERSITY OF CAMBRIDGE CENTRE FOR MATHEMATICAL SCIENCES DEPARTMENT OF APPLIED MATHEMATICS & THEORETICAL PHYSICS Splittings, commutators and the linear Schrödinger equation Philipp Bader, Arieh Iserles,
More information2.161 Signal Processing: Continuous and Discrete
MI OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 8 For information about citing these materials or our erms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSES INSIUE
More informationLie algebraic quantum control mechanism analysis
Lie algebraic quantum control mechanism analysis June 4, 2 Existing methods for quantum control mechanism identification do not exploit analytic features of control systems to delineate mechanism (or,
More informationSpatial Statistics with Image Analysis. Lecture L08. Computer exercise 3. Lecture 8. Johan Lindström. November 25, 2016
C3 Repetition Creating Q Spectral Non-grid Spatial Statistics with Image Analysis Lecture 8 Johan Lindström November 25, 216 Johan Lindström - johanl@maths.lth.se FMSN2/MASM25L8 1/39 Lecture L8 C3 Repetition
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 informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationDIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS. Eskil Hansen and Alexander Ostermann
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume X, Number 0X, XX 00X pp. 1 11 DIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS Eskil Hansen and Alexander Ostermann Institut
More informationSignal and systems. Linear Systems. Luigi Palopoli. Signal and systems p. 1/5
Signal and systems p. 1/5 Signal and systems Linear Systems Luigi Palopoli palopoli@dit.unitn.it Wrap-Up Signal and systems p. 2/5 Signal and systems p. 3/5 Fourier Series We have see that is a signal
More informationCalculation of Highly Oscillatory Integrals by Quadrature Method
Calculation of Highly Oscillatory Integrals by Quadrature Methods Krishna Thapa Texas A&M University, College Station, TX May 18, 211 Outline Purpose Why Filon Type Quadrature Methods Our Expectations
More informationA posteriori error estimation for Magnus-type integrators
ASC Report No. /28 A posteriori error estimation for Magnus-type integrators W. Auzinger, H. Hofsta tter, O. Koch, M. Quell and M. Thalhammer Institute for Analysis and Scientific Computing Vienna University
More informationmultistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP):
MATH 351 Fall 217 multistep methods http://www.phys.uconn.edu/ rozman/courses/m351_17f/ Last modified: November 28, 217 Recall that we are interested in the numerical solution of the initial value problem
More informationwe get y 2 5y = x + e x + C: From the initial condition y(0) = 1, we get 1 5 = 0+1+C; so that C = 5. Completing the square to solve y 2 5y = x + e x 5
Math 24 Final Exam Solution 17 December 1999 1. Find the general solution to the differential equation ty 0 +2y = sin t. Solution: Rewriting the equation in the form (for t 6= 0),we find that y 0 + 2 t
More informationSignals and Spectra (1A) Young Won Lim 11/26/12
Signals and Spectra (A) Copyright (c) 202 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version.2 or any later
More information1.17 : Consider a continuous-time system with input x(t) and output y(t) related by y(t) = x( sin(t)).
(Note: here are the solution, only showing you the approach to solve the problems. If you find some typos or calculation error, please post it on Piazza and let us know ).7 : Consider a continuous-time
More informationSOLUTIONS for Homework #2. 1. The given wave function can be normalized to the total probability equal to 1, ψ(x) = Ne λ x.
SOLUTIONS for Homework #. The given wave function can be normalized to the total probability equal to, ψ(x) = Ne λ x. () To get we choose dx ψ(x) = N dx e λx =, () 0 N = λ. (3) This state corresponds to
More informationLinear Algebra and ODEs review
Linear Algebra and ODEs review Ania A Baetica September 9, 015 1 Linear Algebra 11 Eigenvalues and eigenvectors Consider the square matrix A R n n (v, λ are an (eigenvector, eigenvalue pair of matrix A
More informationAlgebraic structures of B-series
Algebraic structures of B-series Philippe Chartier, Ernst Hairer, and Gilles Vilmart November 10, 2009 Abstract B-series are a fundamental tool in practical and theoretical aspects of numerical integrators
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationExponential Integrators
Exponential Integrators John C. Bowman and Malcolm Roberts (University of Alberta) June 11, 2009 www.math.ualberta.ca/ bowman/talks 1 Outline Exponential Integrators Exponential Euler History Generalizations
More informationLINEAR SYSTEMS. J. Elder PSYC 6256 Principles of Neural Coding
LINEAR SYSTEMS Linear Systems 2 Neural coding and cognitive neuroscience in general concerns input-output relationships. Inputs Light intensity Pre-synaptic action potentials Number of items in display
More informationOrdinary Differential Equations
Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,
More informationMini-course on Geometric Numerical Integration
Mini-course on Geometric Numerical Integration Umeå University and University of Innsbruck david.cohen@umu.se http://snovit.math.umu.se/~david/ This mini-course is supported by an initiation grant from
More informationValidated Explicit and Implicit Runge-Kutta Methods
Validated Explicit and Implicit Runge-Kutta Methods Alexandre Chapoutot joint work with Julien Alexandre dit Sandretto and Olivier Mullier U2IS, ENSTA ParisTech 8th Small Workshop on Interval Methods,
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationChapter 4 The Fourier Series and Fourier Transform
Chapter 4 The Fourier Series and Fourier Transform Fourier Series Representation of Periodic Signals Let x(t) be a CT periodic signal with period T, i.e., xt ( + T) = xt ( ), t R Example: the rectangular
More information13. Power Spectrum. For a deterministic signal x(t), the spectrum is well defined: If represents its Fourier transform, i.e., if.
For a deterministic signal x(t), the spectrum is well defined: If represents its Fourier transform, i.e., if jt X ( ) = xte ( ) dt, (3-) then X ( ) represents its energy spectrum. his follows from Parseval
More informationMarlis Hochbruck 1, Michael Hönig 1 and Alexander Ostermann 2
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher REGULARIZATION OF NONLINEAR ILL-POSED PROBLEMS BY EXPONENTIAL INTEGRATORS Marlis
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationCourse and Wavelets and Filter Banks
Course 18.327 and 1.130 Wavelets and Filter Banks Refinement Equation: Iterative and Recursive Solution Techniques; Infinite Product Formula; Filter Bank Approach for Computing Scaling Functions and Wavelets
More informationAttempts at relativistic QM
Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationRunge-Kutta and Collocation Methods Florian Landis
Runge-Kutta and Collocation Methods Florian Landis Geometrical Numetric Integration p.1 Overview Define Runge-Kutta methods. Introduce collocation methods. Identify collocation methods as Runge-Kutta methods.
More informationAspects of Continuous- and Discrete-Time Signals and Systems
Aspects of Continuous- and Discrete-Time Signals and Systems C.S. Ramalingam Department of Electrical Engineering IIT Madras C.S. Ramalingam (EE Dept., IIT Madras) Networks and Systems 1 / 45 Scaling the
More informationGEOMETRIC INTEGRATION METHODS THAT PRESERVE LYAPUNOV FUNCTIONS
BIT Numerical Mathematics 25, Vol. 45, No., pp. 79 723 c Springer GEOMETRIC INTEGRATION METHODS THAT PRESERVE LYAPUNOV FUNCTIONS V. GRIMM and G.R.W. QUISPEL Department of Mathematics, and Centre for Mathematics
More informationNumerical methods for eigenvalue problems in the description of drift instabilities in the plasma edge
Numerical methods for eigenvalue problems in the description of drift instabilities in the plasma edge Dominik Löchel Supervisors: M. Hochbruck and M. Tokar Graduate school Dynamic of hot plasmas Mathematisches
More informationMatrix power converters: spectra and stability
Matrix power converters: spectra and stability Stephen Cox School of Mathematical Sciences, University of Nottingham supported by EPSRC grant number EP/E018580/1 Making It Real Seminar, Bristol 2009 Stephen
More informationStability of Krylov Subspace Spectral Methods
Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot
More informationModelling and Specifying Dispersive Laser Cavities
Modelling and Specifying Dispersive Laser Cavities C. Sean Bohun (UOIT), Yuri Cher (Toronto), Linda J. Cummings (NJIT), Mehran Ebrahimi (UOIT), Peter Howell (Oxford), Laurent Monasse (CERMICS-ENPC), Judith
More information4 The Continuous Time Fourier Transform
96 4 The Continuous Time ourier Transform ourier (or frequency domain) analysis turns out to be a tool of even greater usefulness Extension of ourier series representation to aperiodic signals oundation
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the
More informationModule 4. Related web links and videos. 1. FT and ZT
Module 4 Laplace transforms, ROC, rational systems, Z transform, properties of LT and ZT, rational functions, system properties from ROC, inverse transforms Related web links and videos Sl no Web link
More informationFourier Series. Spectral Analysis of Periodic Signals
Fourier Series. Spectral Analysis of Periodic Signals he response of continuous-time linear invariant systems to the complex exponential with unitary magnitude response of a continuous-time LI system at
More informationTHE UNIVERSITY OF WESTERN ONTARIO. Applied Mathematics 375a Instructor: Matt Davison. Final Examination December 14, :00 12:00 a.m.
THE UNIVERSITY OF WESTERN ONTARIO London Ontario Applied Mathematics 375a Instructor: Matt Davison Final Examination December 4, 22 9: 2: a.m. 3 HOURS Name: Stu. #: Notes: ) There are 8 question worth
More informationQuantum Physics III (8.06) Spring 2016 Assignment 3
Quantum Physics III (8.6) Spring 6 Assignment 3 Readings Griffiths Chapter 9 on time-dependent perturbation theory Shankar Chapter 8 Cohen-Tannoudji, Chapter XIII. Problem Set 3. Semi-classical approximation
More informationNANOSCALE SCIENCE & TECHNOLOGY
. NANOSCALE SCIENCE & TECHNOLOGY V Two-Level Quantum Systems (Qubits) Lecture notes 5 5. Qubit description Quantum bit (qubit) is an elementary unit of a quantum computer. Similar to classical computers,
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More information