Theory and applications of time-frequency analysis
|
|
- Joel Tucker
- 6 years ago
- Views:
Transcription
1 Theory and applications of time-frequency analysis Ville Turunen Abstract: When and how often something happens in a signal? By properly quantizing these questions, we obtain the Born-Jordan time-frequency transform, defining a sharp phase-space energy density. We study properties of different time-frequency transforms, and also present computed examples from acoustic signal processing, quantum mechanics and medical sciences. Mathematical optics, image modelling and algorithms , Hannover 1 / 26
2 Waveform of speech... 1 Signal (word Why ), sampling rate 4000 Hz / 26
3 Wigner for speech (compare to next Born Jordan...) 700 "Why?": Discrete time Wigner distribution / 26
4 ... Born Jordan for speech (compare to next spectrogram...) 700 "Why?": Discrete time Born Jordan distribution / 26
5 ... a spectrogram (compare to previous Born Jordan) 700 Spectrogram with a medium wide Gaussian window / 26
6 ... a spectrogram with wide Gaussian window Spectrogram with a wide Gaussian window / 26
7 ... a spectrogram with narrow Gaussian window 700 Spectrogram with a narrow Gaussian window / 26
8 Some history Tools: Fourier analysis, FFT! 1920s foundations of quantum mechanics (Heisenberg, Born, Jordan, Schrödinger; Dirac, Wigner, Weyl, von Neumann) Spectrograms (early 1940s Bell Labs, Gabor) 1966 Leon Cohen s class of time-frequency distributions (physics, signal processing) 1960s pseudodifferential operators (Hörmander et al) 2010 Born-Jordan L 2 -continuity (Boggiatto De Donno Oliaro) 2011 Born-Jordan uncertainty (Boggiatto Oliaro Carypis) 8 / 26
9 Order on R characterizes Born Jordan [V.T.] Studying a quantum particle on the real line R, we ask: (A) Is particle on the right? (B) Is particle moving right? So, we just ask for the directions of location and of movement (Separate right from left, up from down, future from past...) In other words, the order relation on R is essential here. The uncertainty in (A, B) characterizes the Born Jordan transform, leading to sharp time-frequency (phase-space) analysis. Other time-frequency transforms do not properly answer to (A, B). This has important consequences in signal processing and quantum mechanics. 9 / 26
10 Direction of position Wavefunction ψ H = L 2 (R) describing a quantum particle on R. By Max Born, probability of finding position right of x R is x ψ(y) 2 dy [0, 1]. (1) Localization to [x, ) R is given by projection A x : H H, { u(y) when y x, A x u(y) := (2) 0 otherwise. Observable A x : Is particle having position right of x? Yes (eigenvalue 1) with probability x ψ(y) 2 dy, and we find the updated wavefunction u/ u with u = A x ψ. No (eigenvalue 0) with probability x ψ(y) 2 dy, and the wavefunction becomes v/ v with v = ψ A x ψ. 10 / 26
11 Direction of momentum Change from position x to momentum η by Fourier transform. By Max Born, probability of finding momentum above η R is η ψ(ξ) 2 dξ [0, 1]. (3) Localization to [η, ) R is given by projection B η : H H, B η u(ξ) := {û(ξ) when ξ η, 0 otherwise. (4) Observable B η : Is particle having momentum above η? Yes (eigenvalue 1) with probability η ψ(ξ) 2 dξ, and we find the updated wavefunction u/ u with u = B η ψ. No (eigenvalue 0) with probability η ψ(ξ) 2 dξ. and the wavefunction becomes v/ v with v = ψ B η ψ. 11 / 26
12 Expectation of directional uncertainty Uncertainty observable of the observable pair (A, B) is i2π [A, B] = i2π (AB BA). (5) An application of the Cauchy Schwarz inequality yields i2π [A, B]u, u = 4π Im Au, Bu 4π Au Bu. (6) This gives the Heisenberg uncertainty inequality i2π [A, B] 4π ( A) ( B), (7) where in state ψ observable M has expectation M := Mψ, ψ and uncertainty M := Mψ M ψ. For u, v H, define Q(u, v)(x, η) := i2π [A x, B η ]u, v. (8) We call Q(u, v) : R R C the time-frequency transform, and Q[u] = Q(u, u) : R R R the time-frequency distribution. Q[ψ](x, η) is the expectation of uncertainty of (A x, B η ) in state ψ. Especially: Q[ψ](0, 0) is the expectation of uncertainty in directional location and movement, 12 / 26
13 Born Jordan transform... and a brief calculation yields the Born Jordan transform 1/2 Q(u, v)(x, η) = e i2πy η u(x+(τ+1/2)y) v(x+(τ 1/2)y) dτdy. R 1/2 Q[ψ] = Q(ψ, ψ) is a quasi-probability distribution of ψ, or Q[u] = Q(u, u) is an energy density of u. Alternatively, the Born Jordan transform is given by (9) FQ(u, v)(ξ, y) = sinc(ξ y) FW (u, v)(ξ, y), (10) where F is the symplectic Fourier transform, and the Wigner transform W (u, v) is defined by W (u, v)(x, η) = e i2πy η u(x + y/2) v(x y/2) dy. (11) R Unfortunately, Wigner has a very bad property: it is extremely sensitive to noise (unlike the non-unitary Q. Also, Q is not causal nor positive.) 13 / 26
14 Wigner distribution is sensitive to noise: speech example Wigner distribution / 26
15 ... and corresponding Born Jordan distribution Speech: Discrete time Born Jordan distribution / 26
16 Bound for Born Jordan energy density [V.T. 2016] There is the optimal Born Jordan bound Q(u, v)(x, η) π u v (12) for all (x, η) R R. Especially, Q[u](x, η) π u 2. Proof. Point localization at the origin : Lu, v = Q(u, v)(0, 0) 1 z/2 = u(t + z/2) v(t z/2) dt dz z z/2 = K L (x, y) u(y) v(x) dx dy, with the Schwartz kernel K L (x, y) = { x y 1 if xy < 0, 0 if xy / 26
17 Proof of the bound for energy density From Lu, v = = u(y) v(x) y x u(y)v( x) y + x dy dx + dy dx u(y) v(x) dy dx x y u( y)v(x) dy dx, x + y denoting u = Pu + Nu, where Pu(x) = 1 R +(x) u(x), we get Lu, v Hilbert π Pu Nv + π Nu Pv = π ( Pu, Nu ) ( Nv, Pv ) Cauchy Schwarz π = π u v. Pu 2 + Nu 2 Nv 2 + Pv 2 Especially, Lu, u π u / 26
18 EKG (Electro-KardioGram) 0.6 Channel 1: Waveform. Press any key to continue / 26
19 EKG: Born Jordan distribution (absolute value) Channel 1: Time index interval [56001:72000]. Frequency unit Hz. Comparison un / 26
20 Another characterization for Born Jordan [V.T. 2016] Cohen class transforms: C = k W for some k S (R R). Necessary and sufficient (i,ii,iii) for C = Q: (i) C is scale-invariant. This means that if v(x) = λ 1/2 u(λx) for λ > 0 then C[v](x, η) = C[u](λx, η/λ). (ii) C is position-local. This means that if u(x) = 0 whenever x [a, b] R then C[u](x, η) = 0 whenever x [a, b] R. (iii) C maps Dirac delta comb to Dirac delta grid. This means C[δ Z ](x, η) = δ Z (x) + δ Z (η) 1, where the Dirac delta comb is δ Z (x) = k Z δ k (x) = κ Z e i2πx κ, with δ k being the Dirac delta distribution at k Z. Think δ Z as a ticking-of-a-clock. Notice that F (δ Z ) = δ Z. Wigner distribution has properties (i,ii) but not (iii). Spectrograms satisfy none of the properties (i,ii,iii). 20 / 26
21 Proof idea of Born Jordan characterization Let FC[u] = φ FW [u]. We must show φ(ξ, y) = sinc(ξ y). (1) φ(ξ, y) = ϕ(ξ y) for some ϕ S (R). (2,1) ϕ(x) = 0 for almost all x > 1/2. So we have analytic ϕ Thm. {: C C by Paley Wiener Schwartz (3,2,1) ϕ(k) = 1 if k = 0, 0 if k Z \ {0}. So û = ϕ/sinc : C C is analytic. Then ϕ = u 1 [ 1/2,1/2] gives u(x + k y) dy [ 1/2,1/2] ϕ(x + k) = k Z k Z = u(x) dx = û(0) = ϕ(0) sinc(0) = 1, so ϕ(x) = 1 for almost all x [ 1/2, 1/2]. Thus we obtain 1/2 ϕ(ξ y) = e i2πxy ξ ϕ(x) dx = e i2πxy ξ dx = sinc(ξ y). 1/2 21 / 26
22 Inversion of Born Jordan distribution [V.T.] Mapping u Q[u] is not invertible: Q[λu] = λ 2 Q[u]. But mapping [u] Q[u] is invertible, where the equivalence class [u] = {u : λ C, λ = 1}. It is easy to see invertibility of Wigner [u] W [u], and by (10) also [u] Q[u] can be inverted. Notice the marginals: Q[u](x, η) dη = u(x) 2, Q[u](x, η) dx = û(η) 2. (13) R R Yet these marginals are not enough to find [u]. Suppose u(x) 0 for Schwartz function u. Then we have inversion u(x + h) = h u(x) 1 R(x kh, h) (14) k=0 for all h 0, where 1 is the partial derivative in the first variable, R(x, y) := e i2πy η Q[u](x, η) dη. (15) R 22 / 26
23 Quantization a a C Cohen class transform C gives the corresponding quantization a a C, a C = a C (X, D) from symbols a : R R C to pseudo-differential operators a C. Such linear operator a C = a C (X, D) is defined by the duality u, a C v = C(u, v), a. (16) Properties of C are reflected in a a C. E.g. marginal conditions C[u](x, η) dη = u(x) 2, C[u](x, η) dx = û(η) 2, (17) R mean that if a(x, η) = f (x) and b(x, η) = ĝ(η) then a C u(x) = f (x) u(x) (multiplication) and b C u(x) = g u(x) = R g(x y) u(y) dy (convolution). Especially, if (17) and a(x, η) = x and b(x, η) = η, then a C = X and b C = D, where Xu(x) = x u(x), Du(x) = 1 i2π u (x). R 23 / 26
24 Quantization examples For example, the Wigner transform C = W from formula (11) gives rise to the Weyl quantization a a W : a W u(x) = e i2π(x y) η a( x + y, η) u(y) dy dη. (18) 2 R R Born Jordan transform Q gives Born Jordan quantization a a Q, 1 a Q u(x) = e i2π(x y) η 2 a((τ+ 1 R R 1 2 )x (τ 1 )y, η) dτ u(y) dy dη. 2 2 (19) What makes the Born Jordan quantization unique among quantizations is that if a(x, η) = f (x) and b(x, η) = ĝ(η) then {a, b} Q = i2π [a Q, b Q ], (20) where the Poisson bracket {a, b} is reduced to ( a b {a, b}(x, η) = x η a ) b (x, η) = f (x) ĝ (η). (21) η x (20)-remark: no-go-theorems [Groenewold 1946, van Hove 1951]. 24 / 26
25 MRI + patient (frequency unit 1 Hz) Restore/Delete/+/ /Zoom/Compute. Frequency unit Hz / 26
26 MRI + patient after noise filtering Restore/Delete/+/ /Zoom/Compute. Frequency unit Hz / 26
QUANTUM MECHANICS LIVES AND WORKS IN PHASE SPACE
Two slit experiment The Wigner phase-space quasi-probability distribution function QUANTUM MECHANICS LIVES AND WORKS IN PHASE SPACE A complete, autonomous formulation of QM based on the standard c- number
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationSharp Gårding inequality on compact Lie groups.
15-19.10.2012, ESI, Wien, Phase space methods for pseudo-differential operators Ville Turunen, Aalto University, Finland (ville.turunen@aalto.fi) M. Ruzhansky, V. Turunen: Sharp Gårding inequality on compact
More informationLecture 4 (Sep. 18, 2017)
Lecture 4 8.3 Quantum Theory I, Fall 07 Lecture 4 (Sep. 8, 07) 4. Measurement 4.. Spin- Systems Last time, we said that a general state in a spin- system can be written as ψ = c + + + c, (4.) where +,
More informationFRAMES AND TIME-FREQUENCY ANALYSIS
FRAMES AND TIME-FREQUENCY ANALYSIS LECTURE 5: MODULATION SPACES AND APPLICATIONS Christopher Heil Georgia Tech heil@math.gatech.edu http://www.math.gatech.edu/ heil READING For background on Banach spaces,
More informationQuantum Theory and Group Representations
Quantum Theory and Group Representations Peter Woit Columbia University LaGuardia Community College, November 1, 2017 Queensborough Community College, November 15, 2017 Peter Woit (Columbia University)
More informationThe Calderon-Vaillancourt Theorem
The Calderon-Vaillancourt Theorem What follows is a completely self contained proof of the Calderon-Vaillancourt Theorem on the L 2 boundedness of pseudo-differential operators. 1 The result Definition
More informationElliptic Problems for Pseudo Differential Equations in a Polyhedral Cone
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 9, Number 2, pp. 227 237 (2014) http://campus.mst.edu/adsa Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone
More informationMicrolocal Analysis : a short introduction
Microlocal Analysis : a short introduction Plamen Stefanov Purdue University Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Analysis : a short introduction 1 / 25 Introduction
More informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationLecture 1 Some Time-Frequency Transformations
Lecture 1 Some Time-Frequency Transformations David Walnut Department of Mathematical Sciences George Mason University Fairfax, VA USA Chapman Lectures, Chapman University, Orange, CA 6-10 November 2017
More informationUncertainty and Spectrogram Geometry
From uncertainty...... to localization Spectrogram geometry CNRS & École Normale Supérieure de Lyon, France Erwin Schrödinger Institute, December 2012 * based on joint work with François Auger and Éric
More informationMATH FALL 2014
MATH 126 - FALL 2014 JASON MURPHY Abstract. These notes are meant to supplement the lectures for Math 126 (Introduction to PDE) in the Fall of 2014 at the University of California, Berkeley. Contents 1.
More informationMICROLOCAL ANALYSIS METHODS
MICROLOCAL ANALYSIS METHODS PLAMEN STEFANOV One of the fundamental ideas of classical analysis is a thorough study of functions near a point, i.e., locally. Microlocal analysis, loosely speaking, is analysis
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationPage 684. Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02
Page 684 Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02 Time Transformations Section 12.5 Symmetries: Time Transformations Page 685 Time Translation
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationLecture 3 Kernel properties and design in Cohen s class time-frequency distributions
Lecture 3 Kernel properties and design in Cohen s class time-frequency distributions Time-frequency analysis, adaptive filtering and source separation José Biurrun Manresa 22.02.2011 Time-Frequency representations
More informationWave Packet Representation of Semiclassical Time Evolution in Quantum Mechanics
Wave Packet Representation of Semiclassical Time Evolution in Quantum Mechanics Panos Karageorge, George Makrakis ACMAC Dept. of Applied Mathematics, University of Crete Semiclassical & Multiscale Aspects
More informationPDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation
PDEs, Homework #3 Solutions. Use Hölder s inequality to show that the solution of the heat equation u t = ku xx, u(x, = φ(x (HE goes to zero as t, if φ is continuous and bounded with φ L p for some p.
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More information4 Sobolev spaces, trace theorem and normal derivative
4 Sobolev spaces, trace theorem and normal derivative Throughout, n will be a sufficiently smooth, bounded domain. We use the standard Sobolev spaces H 0 ( n ) := L 2 ( n ), H 0 () := L 2 (), H k ( n ),
More informationA new class of pseudodifferential operators with mixed homogenities
A new class of pseudodifferential operators with mixed homogenities Po-Lam Yung University of Oxford Jan 20, 2014 Introduction Given a smooth distribution of hyperplanes on R N (or more generally on a
More informationDangerous and Illegal Operations in Calculus Do we avoid differentiating discontinuous functions because it s impossible, unwise, or simply out of
Dangerous and Illegal Operations in Calculus Do we avoid differentiating discontinuous functions because it s impossible, unwise, or simply out of ignorance and fear? Despite the risks, many natural phenomena
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationLecture 6. Four postulates of quantum mechanics. The eigenvalue equation. Momentum and energy operators. Dirac delta function. Expectation values
Lecture 6 Four postulates of quantum mechanics The eigenvalue equation Momentum and energy operators Dirac delta function Expectation values Objectives Learn about eigenvalue equations and operators. Learn
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationLINEAR SECOND-ORDER EQUATIONS
LINEAR SECON-ORER EQUATIONS Classification In two independent variables x and y, the general form is Au xx + 2Bu xy + Cu yy + u x + Eu y + Fu + G = 0. The coefficients are continuous functions of (x, y)
More informationPhase Space Formulation of Quantum Mechanics
Phase Space Formulation of Quantum Mechanics Tony Bracken Centre for Mathematical Physics and Department of Mathematics University of Queensland SMFT07, Melbourne, January-February 2007 Lecture 1 Introduction:
More informationOscillatory integrals
Chapter Oscillatory integrals. Fourier transform on S The Fourier transform is a fundamental tool in microlocal analysis and its application to the theory of PDEs and inverse problems. In this first section
More information2 Resolvents and Green s Functions
Course Notes Solving the Schrödinger Equation: Resolvents 057 F. Porter Revision 09 F. Porter Introduction Once a system is well-specified, the problem posed in non-relativistic quantum mechanics is to
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationContinuous quantum states, Particle on a line and Uncertainty relations
Continuous quantum states, Particle on a line and Uncertainty relations So far we have considered k-level (discrete) quantum systems. Now we turn our attention to continuous quantum systems, such as a
More informationApplied and Computational Harmonic Analysis
Appl. Comput. Harmon. Anal. 8 010) 313 319 Contents lists available at ScienceDirect Applied and Computational Harmonic Analysis www.elsevier.com/locate/acha Smoothed affine Wigner transform A. Athanassoulis
More informationELEMENTARY APPLICATIONS OF FOURIER ANALYSIS
ELEMENTARY APPLICATIONS OF FOURIER ANALYSIS COURTOIS Abstract. This paper is intended as a brief introduction to one of the very first applications of Fourier analysis: the study of heat conduction. We
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationA review: The Laplacian and the d Alembertian. j=1
Chapter One A review: The Laplacian and the d Alembertian 1.1 THE LAPLACIAN One of the main goals of this course is to understand well the solution of wave equation both in Euclidean space and on manifolds
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationHHT: the theory, implementation and application. Yetmen Wang AnCAD, Inc. 2008/5/24
HHT: the theory, implementation and application Yetmen Wang AnCAD, Inc. 2008/5/24 What is frequency? Frequency definition Fourier glass Instantaneous frequency Signal composition: trend, periodical, stochastic,
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 6
Strauss PDEs e: Setion.1 - Exerise 1 Page 1 of 6 Exerise 1 Solve u tt = u xx, u(x, 0) = e x, u t (x, 0) = sin x. Solution Solution by Operator Fatorization By fatoring the wave equation and making a substitution,
More informationA Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets
A Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets George A. Hagedorn Happy 60 th birthday, Mr. Fritz! Abstract. Although real, normalized Gaussian wave packets minimize the product
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationOn a class of pseudodifferential operators with mixed homogeneities
On a class of pseudodifferential operators with mixed homogeneities Po-Lam Yung University of Oxford July 25, 2014 Introduction Joint work with E. Stein (and an outgrowth of work of Nagel-Ricci-Stein-Wainger,
More informationNotes. 1 Fourier transform and L p spaces. March 9, For a function in f L 1 (R n ) define the Fourier transform. ˆf(ξ) = f(x)e 2πi x,ξ dx.
Notes March 9, 27 1 Fourier transform and L p spaces For a function in f L 1 (R n ) define the Fourier transform ˆf(ξ) = f(x)e 2πi x,ξ dx. Properties R n 1. f g = ˆfĝ 2. δλ (f)(ξ) = ˆf(λξ), where δ λ f(x)
More information5.4 Given the basis e 1, e 2 write the matrices that represent the unitary transformations corresponding to the following changes of basis:
5 Representations 5.3 Given a three-dimensional Hilbert space, consider the two observables ξ and η that, with respect to the basis 1, 2, 3, arerepresentedby the matrices: ξ ξ 1 0 0 0 ξ 1 0 0 0 ξ 3, ξ
More informationIntroduction to Pseudodifferential Operators
Introduction to Pseudodifferential Operators Mengxuan Yang Directed by Prof. Dean Baskin May, 206. Introductions and Motivations Classical Mechanics & Quantum Mechanics In classical mechanics, the status
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationNOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS
NOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS JARED WUNSCH Note that these lectures overlap with Alex s to a degree, to ensure a smooth handoff between lecturers! Our notation is mostly, but not completely,
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationFiltering and Edge Detection
Filtering and Edge Detection Local Neighborhoods Hard to tell anything from a single pixel Example: you see a reddish pixel. Is this the object s color? Illumination? Noise? The next step in order of complexity
More informationFractional order operators on bounded domains
on bounded domains Gerd Grubb Copenhagen University Geometry seminar, Stanford University September 23, 2015 1. Fractional-order pseudodifferential operators A prominent example of a fractional-order pseudodifferential
More informationMathematical Methods for Physics
Mathematical Methods for Physics Peter S. Riseborough June 8, 8 Contents Mathematics and Physics 5 Vector Analysis 6. Vectors................................ 6. Scalar Products............................
More information4.1. Dirichlet, Poisson and Neumann boundary value problems
Chapter 4 Laplace s equation 4.. Dirichlet, Poisson and Neumann boundary value problems The most commonly occurring form of problem that is associated with Laplace s equation is a boundary value problem,
More informationFourier and Laplace Transforms, Convolution
Fourier and Laplace Transforms, Convolution Contents: 1. The Fourier transform in L 1. 2. Extensions of the Fourier transform. 3. The Laplace transform. 4. Convolution. 5. Fourier Transform and Differential
More informationA non-local problem with integral conditions for hyperbolic equations
Electronic Journal of Differential Equations, Vol. 1999(1999), No. 45, pp. 1 6. ISSN: 17-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu ftp ejde.math.unt.edu (login:
More informationUNIVERSITY OF SURREY FACULTY OF ENGINEERING AND PHYSICAL SCIENCES DEPARTMENT OF PHYSICS. BSc and MPhys Undergraduate Programmes in Physics LEVEL HE2
Phys/Level /1/9/Semester, 009-10 (1 handout) UNIVERSITY OF SURREY FACULTY OF ENGINEERING AND PHYSICAL SCIENCES DEPARTMENT OF PHYSICS BSc and MPhys Undergraduate Programmes in Physics LEVEL HE PAPER 1 MATHEMATICAL,
More informationQuantum Electrodynamics Test
MSc in Quantum Fields and Fundamental Forces Quantum Electrodynamics Test Monday, 11th January 2010 Please answer all three questions. All questions are worth 20 marks. Use a separate booklet for each
More informationHermite functions. Jordan Bell Department of Mathematics, University of Toronto. September 9, 2015
Hermite functions Jordan Bell jordan.bell@gmail.com Department of Mathematics, University of Toronto September 9, 25 Locally convex spaces If V is a vector space and {p α : α A} is a separating family
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationFOURIER TRANSFORMS. 1. Fourier series 1.1. The trigonometric system. The sequence of functions
FOURIER TRANSFORMS. Fourier series.. The trigonometric system. The sequence of functions, cos x, sin x,..., cos nx, sin nx,... is called the trigonometric system. These functions have period π. The trigonometric
More informationKirchhoff, Fresnel, Fraunhofer, Born approximation and more
Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Oberseminar, May 2008 Maxwell equations Or: X-ray wave fields X-rays are electromagnetic waves with wave length from 10 nm to 1 pm, i.e., 10
More informationUNCERTAINTY PRINCIPLES FOR RADAR AMBIGUITY FUNCTIONS AND MOMENTS
UNCTAINTY PINCIPLS FO ADA AMBIGUITY FUNCTIONS AND MOMNTS PHILIPP JAMING Abstract. In this article, we give an estimate of the zero-free region around the origin of the ambiguity function of a one-dimensional
More informationIntroduction to the Mathematics of Medical Imaging
Introduction to the Mathematics of Medical Imaging Second Edition Charles L. Epstein University of Pennsylvania Philadelphia, Pennsylvania EiaJTL Society for Industrial and Applied Mathematics Philadelphia
More informationDRAFT: April 28, 2008
1 DRAFT: April 28, 2008 1 Lectures on PDEs- APMA- E6301 DRAFT IN PROGRESS Michael I. Weinstein April 28, 2008 Contents 1 First order partial differential equations and the method of characteristics 4 2
More informationarxiv: v1 [eess.sp] 4 Nov 2018
Estimating the Signal Reconstruction Error from Threshold-Based Sampling Without Knowing the Original Signal arxiv:1811.01447v1 [eess.sp] 4 Nov 2018 Bernhard A. Moser SCCH, Austria Email: bernhard.moser@scch.at
More informationIntroduction to pseudodifferential operators. 1. Kohn-Nirenberg calculus, simplest symbol classes S m
February 8, 17) Introduction to pseudodifferential operators Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/fun/notes 16-17/psidos.pdf]
More informationSEMICLASSICAL LAGRANGIAN DISTRIBUTIONS
SEMICLASSICAL LAGRANGIAN DISTRIBUTIONS SEMYON DYATLOV Abstract. These are (rather sloppy) notes for the talk given in UC Berkeley in March 2012, attempting to explain the global theory of semiclassical
More informationPROBLEMS IN UNBOUNDED CYLINDRICAL DOMAINS
PROBLEMS IN UNBOUNDED CYLINDRICAL DOMAINS PATRICK GUIDOTTI Mathematics Department, University of California, Patrick Guidotti, 103 Multipurpose Science and Technology Bldg, Irvine, CA 92697, USA 1. Introduction
More informationWIGNER FUNCTION OF A QUBIT
GEOMETRY OF JETS AND FIELDS BANACH CENTER PUBLICATIONS, VOLUME 110 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 016 WIGNER FUNCTION OF A QUBIT JERZY KIJOWSKI, PIOTR WALUK, KATARZYNA SENGER
More informationTIME-FREQUENCY ANALYSIS: TUTORIAL. Werner Kozek & Götz Pfander
TIME-FREQUENCY ANALYSIS: TUTORIAL Werner Kozek & Götz Pfander Overview TF-Analysis: Spectral Visualization of nonstationary signals (speech, audio,...) Spectrogram (time-varying spectrum estimation) TF-methods
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP
Dedicated to Professor Gheorghe Bucur on the occasion of his 7th birthday ON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP EMIL POPESCU Starting from the usual Cauchy problem, we give
More informationConvolution and Linear Systems
CS 450: Introduction to Digital Signal and Image Processing Bryan Morse BYU Computer Science Introduction Analyzing Systems Goal: analyze a device that turns one signal into another. Notation: f (t) g(t)
More informationON THE PRANDTL EQUATION. 0. Introduction The purpose of this paper is to investigate the integro-differential equation
GEORGIAN MATHEMATICAL JOURNAL: Vol. 6, No. 6, 1999, 525-536 ON THE PRANDTL EQUATION R. DUDUCHAVA AND D. KAPANADZE Abstract. The unique solvability of the airfoil (Prandtl) integrodifferential equation
More informationHölder continuity of the solution to the 3-dimensional stochastic wave equation
Hölder continuity of the solution to the 3-dimensional stochastic wave equation (joint work with Yaozhong Hu and Jingyu Huang) Department of Mathematics Kansas University CBMS Conference: Analysis of Stochastic
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science Quantum Optical Communication
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.453 Quantum Optical Communication Date: Tuesday, September 13, 216 Dirac-notation Quantum Mechanics. Introduction
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 informationHermite functions and uncertainty principles for the Fourier and the windowed Fourier transforms
[BDJam] Revista Matemática Iberoamericana 19 (003) 3 55. Hermite functions and uncertainty principles for the Fourier and the windowed Fourier transforms Aline BONAMI, Bruno DEMANGE & Philippe JAMING Abstract
More informationTopics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality
Topics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality Po-Lam Yung The Chinese University of Hong Kong Introduction While multiplier operators are very useful in studying
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationLaplace equation. In this chapter we consider Laplace equation in d-dimensions given by. + u x2 x u xd x d. u x1 x 1
Chapter 6 Laplace equation In this chapter we consider Laplace equation in d-dimensions given by u x1 x 1 + u x2 x 2 + + u xd x d =. (6.1) We study Laplace equation in d = 2 throughout this chapter (excepting
More informationWave packet decompositions adapted to (non-self-adjoint) operators
1 / 42 Wave packet decompositions adapted to (non-self-adjoint) operators Joe Viola Laboratoire de Mathématiques Jean Leray Université de Nantes Joseph.Viola@univ-nantes.fr 26 March 2018 2 / 42 Outline
More informationMIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain
10/31/05 wk9-a-1 The spatial frequency domain Recall: plane wave propagation x path delay increases linearly with x λ z=0 θ E 0 x exp i2π sinθ + λ z i2π cosθ λ z plane of observation 10/31/05 wk9-a-2 Spatial
More informationVariations on Quantum Ergodic Theorems. Michael Taylor
Notes available on my website, under Downloadable Lecture Notes 8. Seminar talks and AMS talks See also 4. Spectral theory 7. Quantum mechanics connections Basic quantization: a function on phase space
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More informationIntroduction to time-frequency analysis. From linear to energy-based representations
Introduction to time-frequency analysis. From linear to energy-based representations Rosario Ceravolo Politecnico di Torino Dep. Structural Engineering UNIVERSITA DI TRENTO Course on «Identification and
More informationEstimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.
Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations
More informationEvolution equations and microlocal analysis
Evolution equations and microlocal analysis Jean-Michel Bony Introduction Given a self-adjoint unbounded operator A on L 2 (R n ), we know from abstract functional analysis that there exists a one parameter
More informationGreen Function of a Time-Dependent Linear PDE
Green Function of a Time-Dependent Linear PDE Consider one-dimensional linear PDE of the form γu t = u xx, (1) where u = u(x, t), x (, ), u t = u/ t, u xx = 2 u/ x 2. We want to find u(x, t), provided
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 14: The Power Function and Native Space Error Estimates Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationQuantification of Gaussian quantum steering. Gerardo Adesso
Quantification of Gaussian quantum steering Gerardo Adesso Outline Quantum steering Continuous variable systems Gaussian entanglement Gaussian steering Applications Steering timeline EPR paradox (1935)
More informationNumerical methods for a fractional diffusion/anti-diffusion equation
Numerical methods for a fractional diffusion/anti-diffusion equation Afaf Bouharguane Institut de Mathématiques de Bordeaux (IMB), Université Bordeaux 1, France Berlin, November 2012 Afaf Bouharguane Numerical
More informationOn Time-Frequency Sparsity and Uncertainty
CNRS & E cole Normale Supe rieure de Lyon, France * partially based on joint works with Franc ois Auger, Pierre Borgnat and E ric Chassande-Mottin Heisenberg (classical) from or to and Heisenberg refined
More informationMath 5440 Problem Set 6 Solutions
Math 544 Math 544 Problem Set 6 Solutions Aaron Fogelson Fall, 5 : (Logan,.6 # ) Solve the following problem using Laplace transforms. u tt c u xx g, x >, t >, u(, t), t >, u(x, ), u t (x, ), x >. The
More informationGEOMETRIC QUANTIZATION
GEOMETRIC QUANTIZATION 1. The basic idea The setting of the Hamiltonian version of classical (Newtonian) mechanics is the phase space (position and momentum), which is a symplectic manifold. The typical
More information