Discrete Orthogonal Harmonic Transforms
|
|
- Erik Carter
- 6 years ago
- Views:
Transcription
1 Discrete Orthogonal Harmonic Transforms Speaker: Chun-Lin Liu, Advisor: Soo-Chang Pei Ph. D Image Processing Laboratory, EEII 530, Graduate Institute of Communication Engineering, National Taiwan University. May 26, 2012 CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
2 Outline 1 Introduction 2 Computation of 1D Discrete Orthogonal Functions Proposed method to the 1D discrete orthogonal functions Simulation results 3 Discrete implementation of 2D Fourier Transform Eigenfunctions The general form of 2D FT eigenfunctions and the discrete implementation Simulation Results 4 Conclusions CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
3 Outline 1 Introduction 2 Computation of 1D Discrete Orthogonal Functions Proposed method to the 1D discrete orthogonal functions Simulation results 3 Discrete implementation of 2D Fourier Transform Eigenfunctions The general form of 2D FT eigenfunctions and the discrete implementation Simulation Results 4 Conclusions CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
4 Introduction Some well-known discrete orthogonal transforms are Discrete Fourier transform (DFT), Discrete Sine/Cosine/Hartley transforms, Discrete wavelet transforms. Properties of orthogonal transforms More coefficients lead to smaller reconstruction error. Perfect reconstruction. Problem Explicit definitions of the transform kernels. For instance, the DFT is defined as F k = N 1 n=0 2πnk j f n e N, n, k = 0, 1, 2,..., N 1. (1) CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
5 A simple method to other discrete transforms Choose a continuous orthogonal transform and take samples of the transform kernel. Pros: Good approximation to the continuous transforms. Cons: Not orthogonal. For instance, assume that the continuous transform kernel {ψ n (x)} n=0 satisfies the orthogonal relation: ψ m(x)ψ n (x) d x = δ m,n. (2) However, taking the samples x = l x does not ensure orthogonality: ψm (l x ) ψ n (l x ) x δ m,n. (3) l CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
6 Problem formulation The 1D discrete orthogonal transform is characterized by the transform matrix Ψ, as the following relationship F 0 F =.. = [ f 0 ] H ψ 0... ψ. N 1. = Ψ H f, (4) F N 1 f N 1 f is the input signal, F is the transformed signal, {ψ n } N 1 n=0 are the discrete transform kernels, f, F, and ψ n are all N-by-1 column vectors. Problem statement We want to the discrete orthogonal transform to be General ( there are general ways to find entries of Ψ). Good approximation (ψ n is close to the samples of ψ n (x)). Orthogonal (Ψ is an unitary matrix). CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
7 Outline 1 Introduction 2 Computation of 1D Discrete Orthogonal Functions Proposed method to the 1D discrete orthogonal functions Simulation results 3 Discrete implementation of 2D Fourier Transform Eigenfunctions The general form of 2D FT eigenfunctions and the discrete implementation Simulation Results 4 Conclusions CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
8 Hermite Gaussian functions and Hermite transforms The Hermite Gaussian functions (HGFs) are defined by ( ) 1 1/2 h n (x) = 2 n n! H n (x)e x2 /2, (5) π where H n (x) is the Hermite polynomial. Differential equation (D x = d / d x is the differential operator): ( D 2 x x 2) h n (x) = (2n + 1)h n (x). (6) Properties: hn (x) is the eigenfunction of the Fourier transform with eigenvalue ( j) n. Complete and orthonormal basis for L 2 (R). Hermite transforms are defined by a n = f(x)h n (x) d x. (7) CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
9 Outline 1 Introduction 2 Computation of 1D Discrete Orthogonal Functions Proposed method to the 1D discrete orthogonal functions Simulation results 3 Discrete implementation of 2D Fourier Transform Eigenfunctions The general form of 2D FT eigenfunctions and the discrete implementation Simulation Results 4 Conclusions CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
10 Implement discrete HGFs by differential equations Implement the differential equation as the matrix eigen-problem. 1 Convert L into L. 2 Solve the eigenvectors and eigenvalues of L numerically. Continuous: L ({}}{ D 2 x x 2) h n (x) = (2n + 1) h n (x),... Discrete: = λ n.. Three basic operations }... {{ } L }{{} h n 1 Multiplied by a constant c. (ci) f(x) = cf(x). 2 Multiplied by x. X f(x) = xf(x). 3 Differentiate with respect to x, D x = d / d x. }{{} h n CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
11 Discrete version of ci and X 1 ci N is the discrete version of ci. I N is the identity matrix of size N. Continuous: ci f(x) = cf(x), c f 1 cf 1 0 c... 0 f 2 Discrete: = cf c f N 1 cf N 1 2 Assume that the functions are sampled uniformly on x = n x, n = [ N 1 2, N 3 2,..., N 1 2 ]T. X = diag(x) represents the discrete version of X. X f(x) = xf(x), N f 1 N 1 0 N 3 2 xf f 2 x = N 3 2 xf 2.. N f N 1 2 N 1 2 xf N 1 CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
12 Discrete version of differential operators In numerical analysis, there are forward difference, backward difference, and central difference. However, they are simple but inaccurate. We derive the differential operator from Fourier transforms. D x f(x) = F 1 FD x f(x) = F 1 (jx ) Ff(x), (8) where F denotes the Fourier transform operator. Proposed discrete differential operator 3 The discrete differential operator is obtained by D = F 1 (jx) F, (9) where F is the DFT matrix and X is the discrete multiplied-by-x operator. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
13 Outline 1 Introduction 2 Computation of 1D Discrete Orthogonal Functions Proposed method to the 1D discrete orthogonal functions Simulation results 3 Discrete implementation of 2D Fourier Transform Eigenfunctions The general form of 2D FT eigenfunctions and the discrete implementation Simulation Results 4 Conclusions CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
14 Simulation on the Discrete HGFs The number of discrete points N = 101. The sample interval x = 2π/N The eigenvalues are all negative and sorted in the descent order. We want to verify whether The eigenvectors are close to the continuous samples. Note that the results are identical to those of the n 2 -matrix 1. 1 S. C. Pei, J. J. Ding, W. L. Hsue and K. W. Chang, Generalized commuting matrices and their eigenvectors for DFTs, offset DFTs, and other periodic operations, IEEE Trans. on Signal Processing, vol.56, No.8, pp , Aug CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
15 First four Discrete HGFs Discrete Continuous Discrete Continuous ψ 0 (x) ψ 1 (x) x Discrete Continuous x Discrete Continuous ψ 2 (x) ψ 3 (x) x x CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
16 Outline 1 Introduction 2 Computation of 1D Discrete Orthogonal Functions Proposed method to the 1D discrete orthogonal functions Simulation results 3 Discrete implementation of 2D Fourier Transform Eigenfunctions The general form of 2D FT eigenfunctions and the discrete implementation Simulation Results 4 Conclusions CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
17 2D discrete orthogonal functions Another question is raised: How do we implement high dimensional discrete orthogonal functions while the orthogonality is kept? Possible approach Take the samples of the continuous functions. Very simple but not orthogonal. Start from the differential equations. Partial differential equations require a very huge matrix. Large scale numerical eigen-decomposition is very difficult. Different coordinate systems lead to multiple solutions. Due to these difficulties, we narrow our discussion to the eigenfunctions of 2D Fourier transforms 2. 2 S. C. Pei and C. L. Liu, A general form of 2D Fourier transform eigenfunctions, 2012 IEEE International Conference on Acoustics, Speech and Signal Processing, Kyoto, Japan,, March CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
18 Known eigenfunctions of 2D FT There are three known eigenfunctions for 2D FT. Separable HGFs (SHGFs) in Cartesian coordinates h m,n (x, y) = h m (x)h n (y). (10) Rotated HGFs (RHGFs) in rotated Cartesian coordinates h m,n (α; x, y) = h m,n (x cos α + y sin α, x sin α + y cos α). (11) Laguerre-Gaussian functions (LGFs) in polar coordinates l m,n (r, θ) = N p,l r l L l ( p r 2 ) e r2 /2 e jlθ, (12) where p = min {m, n}, l = m n, N p,l is the normalization factor, and L l p ( ) is the associated Laguerre polynomial. All of them form a complete and orthonormal basis for L 2 (R 2 ). CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
19 Observations on the 2D Fourier transform eigenfunctions Their differential equations have identical operators, identical eigenvalues, but different eigenfunctions. ( ) 2 x y 2 x2 y 2 h m,n (x, y) = 2 (m + n + 1) h m,n (x, y), ( ) 2 x y 2 x2 y 2 h m,n (α; x, y) = 2 (m + n + 1) h m,n (α; x, y), ( 1 r r r r ) r 2 θ 2 r2 l m,n (r, θ) = 2 (m + n + 1) l m,n (r, θ). Eigenvalues of the 2D FT are the same, ( j) m+n, for each solution. Our question Is there a general form to represent these eigenfunctions? CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
20 Outline 1 Introduction 2 Computation of 1D Discrete Orthogonal Functions Proposed method to the 1D discrete orthogonal functions Simulation results 3 Discrete implementation of 2D Fourier Transform Eigenfunctions The general form of 2D FT eigenfunctions and the discrete implementation Simulation Results 4 Conclusions CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
21 The general form of the 2D FT eigenfunctions The 2D Fourier transform eigenfunctions ψ m,n (x, y) is ψ m,n (x, y) = ψ 0,m+n(x, y) ψ 1,m+n 1(x, y) =. ψ m+n,0(x, y) m+n c m,n p h p,m+n p (x, y), (13) p=0 c 0,m+n 0 c 0,m+n 1... c 0,m+n m+n c 1,m+n 1 0 c 1,m+n c 1,m+n 1 m+n } c m+n,0 0 c m+n,0 1 {{... c m+n,0 m+n } T h 0,m+n(x, y) h 1,m+n 1(x, y). h m+n,0(x, y) m, n = 0, 1, 2,.... c m,n p are the combination coefficients. ψ m,n (x, y) corresponds eigenvalue ( j) m+n of the 2D FT. Mathematical interpretation: linear combination of the eigenfunctions with the same eigenvalue/eigenspace gives another eigenfunctions. ψ m,n (x, y) are orthogonal if and only if T is unitary. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
22 Known combination coefficients The SHGFs, RHGFs, and LGFs all satisfy the general for of 2D Fourier transform eigenfunctions. SHGFs h m,n (x, y) and c m,n p = δ[p m]. RHGFs h m,n (α; x, y) 3, c m,n p = p!(m + n p)! m!n! where P n (α,β) ( ) is the Jacobi polynomial. LGFs l m,n (r, θ) 4, (sin α) m p (cos α) n p P (m p,n p) p (cos 2α), c m,n p (LGF) = j p c m,n p (RHGF) α=π/4. (14) 3 A. Wünsche, Hermite and Laguerre 2D polynomials, Journal of Computational and Applied Mathematics, vol. 133, pp , M. W. Beijersbergen, L. Allen, H. van der Veen, and J. P. Woerdman, Astigmatic laser mode converters and transform of orbital angular momentum, Opt. Comm., vol. 96, pp , CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
23 Discrete 2D Fourier transform eigenfunctions Implementation of discrete 2D Fourier transform eigenfunctions 1 Solve the 1D discrete Hermite Gaussian functions. 2 Construct discrete separable HGFs by multiplying 1D discrete HGFs in each dimension. 3 Compute combination coefficients c m,n p. 4 Obtain discrete 2D Fourier transform eigenfunctions according to ψ m,n (x, y) = m+n p=0 c m,n p h p,m+n p (x, y). Advantages: More general than the explicit expression. Different combination coefficients yield different eigenfunctions. Samples are always on the Cartesian grids. Orthogonal, as long as the matrix composed of c m,n p is unitary. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
24 Problem on very high order discrete eigenfunctions An illustrative example The number of discrete points in each dimension N = 6. We want to get the discrete function of h 1,5 (α; x, y). The implementation steps indicates 1 We have only 6 1D discrete HGFs: h 0, h 1, h 2, h 3, h 4, h 5. 2 Construct separable HGFs. Then compute the combination coefficients. 3 We have a problem in the last step. Write out the general form c 2,4 0 h 6h T 0 + c 2,4 1 h 5h T 1 + c 2,4 2 h 4h T 2 + c 2,4 3 h 3h T 3 +c 2,4 4 h 2h T 4 + c 2,4 5 h 1h T 5 + c 2,4 6 h 0h T 6 The problem is the absence of h 6 in the 1D discrete HGFs. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
25 Two possible solutions to high-order discrete eigenfunctions 1 Truncate the combination coefficients. h 1,5 (α; x, y) c 2,4 0 h 6h T 0 + c 2,4 1 h 5h T c 2,4 5 h 1h T 5 + c 2,4 6 h 0h T 6 Not orthogonal High accuracy 2 Replace the combination coefficients with those of lower-order eigenfunctions. The concept is to mirror the coefficients. h 0,4 (α; x, y) d 2,4 0 h 4h T d 2,4 4 h 0h T 4 h 1,5 (α; x, y) c 2,4 0 h 6h T 0 +c 2,4 1 h 5h T c 2,4 5 h 1h T 5 + c 2,4 6 h 0h T 6 Orthogonal Low accuracy CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
26 The discrete Laguerre Gaussian transform The Laguerre Gaussian transform is L m,n = f(x, y)lm,n(r, θ) d 2 r. (15) R 2 According to the general form, (13), we have m+n ( ) L m,n = c m,n p f(x, y)h p,m+n p (x, y) d 2 r. (16) p=0 R } 2 {{} 2D separable Hermite transform The Laguerre Gaussian transform is composed of 1 2D separable Hermite transform 2 Linear transformation associated with c m,n p. In terms of discrete implementation: 1 Separable transforms and the linear transformation are faster than non-separable transforms. 2 High-order L m,n are adjusted either by truncating or by mirroring the transformation coefficients c m,n p. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
27 Outline 1 Introduction 2 Computation of 1D Discrete Orthogonal Functions Proposed method to the 1D discrete orthogonal functions Simulation results 3 Discrete implementation of 2D Fourier Transform Eigenfunctions The general form of 2D FT eigenfunctions and the discrete implementation Simulation Results 4 Conclusions CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
28 Simulation results 1 Compare the continuous-sampled Laguerre Gaussian functions with the discrete Laguerre Gaussian functions. 2 Compare the high-order continuous-sampled LGFs with the discrete LGFs with either truncated coefficients or mirrored coefficients. 3 Image expansion and reconstruction using discrete Laguerre Gaussian transforms. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
29 Simulation 1: Accuracy of the discrete eigenfunctions l 8,3 (r, θ) Magnitude Phase Cont Number of points in each dimension N = 101. Sampling interval x = y = 2π/N. Disc Error comes from very small values ( ) on the boundaries CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
30 Simulation 2: High-order discrete LGFs Continuous samples l 63,1 (r, θ) Truncate Mirror Number of points in each dimension N = 64. Sampling interval x = y = 2π/N. Compare the magnitudes of l m,n (r, θ). The number of rings should be min {m, n} + 1. For the mirrored discrete function, coefficients are borrowed from those of l 62,0 (r, θ). Error comes from Inaccurate combination terms. Inaccurate high-order Hermite Gaussian functions. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
31 Simulation 3: Image expansion and reconstruction The input image I Reconstruction schemes n L m L m,n K Image expansion: Convert input image I into its LGT L m,n. Image reconstruction: partial L m,n and the inverse LGT yields the reconstructed image Î. N = 90, x = y = 2π/N. The input image I is the 2D squared image. We have two type of reconstruction schemes, shown on the left: Fix L, increase K gradually. Fix K, increase L gradually. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
32 Image reconstruction along K (details in θ-direction) The input image I discrete LGT = Reconstruction scheme L m,n n m K = 0, 2, 4, 6 K K = 8, 10, 12, 14 CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
33 Image reconstruction along L (details in r-direction) The input image I discrete LGT = Reconstruction scheme L m n L m,n L = 0, 2, 4, 6 L = 8, 10, 12, 14 CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
34 Outline 1 Introduction 2 Computation of 1D Discrete Orthogonal Functions Proposed method to the 1D discrete orthogonal functions Simulation results 3 Discrete implementation of 2D Fourier Transform Eigenfunctions The general form of 2D FT eigenfunctions and the discrete implementation Simulation Results 4 Conclusions CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
35 Conclusion on 1D discrete orthogonal functions In the first part, we proposed an alternative method to the discrete orthogonal functions. As long as the differential equations are known, we can replace the continuous operators L with discrete operator L. This approach is more general and the explicit form is not required. We can apply this procedure to solve any linear differential equations. Classical orthogonal polynomials Schrödinger equations in quantum mechanics Scale transforms Fractional Fourier transforms Linear canonical transforms and their eigenfunctions We utilize the 1D results to implement the 2D case. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
36 Conclusion on 2D Fourier transform eigenfunctions In the section part, we implement the 2D Fourier transform eigenfunctions along with the transform derived from these eigenfunctions. The general form of 2D Fourier transform eigenfunctions. Deal with high-order discrete eigenfunctions. We can truncate or mirror the combination coefficients. Convert the non-separable transform (LGT) into 2D separable Hermite transforms and then a linear transformation (fast). The concept can be generalized to three-dimensional eigenfunctions. CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
37 Thank you for your attention! Q & A Time CL-Liu (GICE, NTU) Discrete Orthogonal Harmonic Transforms May 26, / 37
A new method of image encryption with. multiple-parameter discrete fractional Fourier transform
2012 International Conference on Information and Computer Networks (ICICN 2012) IPCSIT vol. 27 (2012) (2011) IACSIT Press, Singapore A new method of image encryption with multiple-parameter discrete fractional
More informationAngular momentum. Quantum mechanics. Orbital angular momentum
Angular momentum 1 Orbital angular momentum Consider a particle described by the Cartesian coordinates (x, y, z r and their conjugate momenta (p x, p y, p z p. The classical definition of the orbital angular
More informationMathematical Physics Homework 10
Georgia Institute of Technology Mathematical Physics Homework Conner Herndon November, 5 Several types of orthogonal polynomials frequently occur in various physics problems. For instance, Hermite polynomials
More informationSpin Dynamics Basic Theory Operators. Richard Green SBD Research Group Department of Chemistry
Spin Dynamics Basic Theory Operators Richard Green SBD Research Group Department of Chemistry Objective of this session Introduce you to operators used in quantum mechanics Achieve this by looking at:
More informationPage 712. Lecture 42: Rotations and Orbital Angular Momentum in Two Dimensions Date Revised: 2009/02/04 Date Given: 2009/02/04
Page 71 Lecture 4: Rotations and Orbital Angular Momentum in Two Dimensions Date Revised: 009/0/04 Date Given: 009/0/04 Plan of Attack Section 14.1 Rotations and Orbital Angular Momentum: Plan of Attack
More informationSample Quantum Chemistry Exam 2 Solutions
Chemistry 46 Fall 7 Dr. Jean M. Standard Name SAMPE EXAM Sample Quantum Chemistry Exam Solutions.) ( points) Answer the following questions by selecting the correct answer from the choices provided. a.)
More informationNIH Public Access Author Manuscript Proc Soc Photo Opt Instrum Eng. Author manuscript; available in PMC 2013 December 17.
NIH Public Access Author Manuscript Published in final edited form as: Proc Soc Photo Opt Instrum Eng. 2008 March 18; 6913:. doi:10.1117/12.769604. Tomographic Reconstruction of Band-limited Hermite Expansions
More informationSection 9 Variational Method. Page 492
Section 9 Variational Method Page 492 Page 493 Lecture 27: The Variational Method Date Given: 2008/12/03 Date Revised: 2008/12/03 Derivation Section 9.1 Variational Method: Derivation Page 494 Motivation
More informationProblem 1: Spin 1 2. particles (10 points)
Problem 1: Spin 1 particles 1 points 1 Consider a system made up of spin 1/ particles. If one measures the spin of the particles, one can only measure spin up or spin down. The general spin state of a
More informationRepresentation of the quantum and classical states of light carrying orbital angular momentum
Representation of the quantum and classical states of light carrying orbital angular momentum Humairah Bassa and Thomas Konrad Quantum Research Group, University of KwaZulu-Natal, Durban 4001, South Africa
More informationLecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor
Lecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor It turns out that the boundary condition of the wavefunction going to zero at infinity is sufficient to quantize the value of energy that
More informationPhysics 137A Quantum Mechanics Fall 2012 Midterm II - Solutions
Physics 37A Quantum Mechanics Fall 0 Midterm II - Solutions These are the solutions to the exam given to Lecture Problem [5 points] Consider a particle with mass m charge q in a simple harmonic oscillator
More information8.04 Spring 2013 March 12, 2013 Problem 1. (10 points) The Probability Current
Prolem Set 5 Solutions 8.04 Spring 03 March, 03 Prolem. (0 points) The Proaility Current We wish to prove that dp a = J(a, t) J(, t). () dt Since P a (t) is the proaility of finding the particle in the
More informationSimple Harmonic Oscillator
Classical harmonic oscillator Linear force acting on a particle (Hooke s law): F =!kx From Newton s law: F = ma = m d x dt =!kx " d x dt + # x = 0, # = k / m Position and momentum solutions oscillate in
More informationMath 118, Handout 4: Hermite functions and the Fourier transform. n! where D = d/dx, are a basis of eigenfunctions for the Fourier transform
The Hermite functions defined by h n (x) = ( )n e x2 /2 D n e x2 where D = d/dx, are a basis of eigenfunctions for the Fourier transform f(k) = f(x)e ikx dx on L 2 (R). Since h 0 (x) = e x2 /2 h (x) =
More informationG : Quantum Mechanics II
G5.666: Quantum Mechanics II Notes for Lecture 5 I. REPRESENTING STATES IN THE FULL HILBERT SPACE Given a representation of the states that span the spin Hilbert space, we now need to consider the problem
More informationOne-electron Atom. (in spherical coordinates), where Y lm. are spherical harmonics, we arrive at the following Schrödinger equation:
One-electron Atom The atomic orbitals of hydrogen-like atoms are solutions to the Schrödinger equation in a spherically symmetric potential. In this case, the potential term is the potential given by Coulomb's
More informationRepresentation theory and quantum mechanics tutorial Spin and the hydrogen atom
Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition
More informationNon Adiabatic Transitions in a Simple Born Oppenheimer Scattering System
1 Non Adiabatic Transitions in a Simple Born Oppenheimer Scattering System George A. Hagedorn Department of Mathematics and Center for Statistical Mechanics, Mathematical Physics, and Theoretical Chemistry
More information1. Matrix multiplication and Pauli Matrices: Pauli matrices are the 2 2 matrices. 1 0 i 0. 0 i
Problems in basic linear algebra Science Academies Lecture Workshop at PSGRK College Coimbatore, June 22-24, 2016 Govind S. Krishnaswami, Chennai Mathematical Institute http://www.cmi.ac.in/~govind/teaching,
More informationSolutions: Problem Set 3 Math 201B, Winter 2007
Solutions: Problem Set 3 Math 201B, Winter 2007 Problem 1. Prove that an infinite-dimensional Hilbert space is a separable metric space if and only if it has a countable orthonormal basis. Solution. If
More informationLinear Algebra in Hilbert Space
Physics 342 Lecture 16 Linear Algebra in Hilbert Space Lecture 16 Physics 342 Quantum Mechanics I Monday, March 1st, 2010 We have seen the importance of the plane wave solutions to the potentialfree Schrödinger
More informationHilbert Space Problems
Hilbert Space Problems Prescribed books for problems. ) Hilbert Spaces, Wavelets, Generalized Functions and Modern Quantum Mechanics by Willi-Hans Steeb Kluwer Academic Publishers, 998 ISBN -7923-523-9
More informationSpecial Functions of Mathematical Physics
Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More informationTotal Angular Momentum for Hydrogen
Physics 4 Lecture 7 Total Angular Momentum for Hydrogen Lecture 7 Physics 4 Quantum Mechanics I Friday, April th, 008 We have the Hydrogen Hamiltonian for central potential φ(r), we can write: H r = p
More informationwhich implies that we can take solutions which are simultaneous eigen functions of
Module 1 : Quantum Mechanics Chapter 6 : Quantum mechanics in 3-D Quantum mechanics in 3-D For most physical systems, the dynamics is in 3-D. The solutions to the general 3-d problem are quite complicated,
More information3.024 Electrical, Optical, and Magnetic Properties of Materials Spring 2012 Recitation 1. Office Hours: MWF 9am-10am or by appointment
Adam Floyd Hannon Office Hours: MWF 9am-10am or by e-mail appointment Topic Outline 1. a. Fourier Transform & b. Fourier Series 2. Linear Algebra Review 3. Eigenvalue/Eigenvector Problems 1. a. Fourier
More informationQuantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras. Lecture - 21 Square-Integrable Functions
Quantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 21 Square-Integrable Functions (Refer Slide Time: 00:06) (Refer Slide Time: 00:14) We
More informationPHYSICS 234 HOMEWORK 2 SOLUTIONS. So the eigenvalues are 1, 2, 4. To find the eigenvector for ω = 1 we have b c.
PHYSICS 34 HOMEWORK SOUTIONS.8.. The matrix we have to diagonalize is Ω 3 ( 4 So the characteristic equation is ( ( ω(( ω(4 ω. So the eigenvalues are,, 4. To find the eigenvector for ω we have 3 a (3 b
More informationarxiv:math/ v2 [math.ap] 3 Oct 2006
THE TAYLOR SERIES OF THE GAUSSIAN KERNEL arxiv:math/0606035v2 [math.ap] 3 Oct 2006 L. ESCAURIAZA From some people one can learn more than mathematics Abstract. We describe a formula for the Taylor series
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 informationQuantum Mechanics: Postulates
Quantum Mechanics: Postulates 25th March 2008 I. Physical meaning of the Wavefunction Postulate 1: The wavefunction attempts to describe a quantum mechanical entity (photon, electron, x-ray, etc.) through
More informationReading: Mathchapters F and G, MQ - Ch. 7-8, Lecture notes on hydrogen atom.
Chemistry 356 017: Problem set No. 6; Reading: Mathchapters F and G, MQ - Ch. 7-8, Lecture notes on hydrogen atom. The H atom involves spherical coordinates and angular momentum, which leads to the shapes
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationThe quantum state as a vector
The quantum state as a vector February 6, 27 Wave mechanics In our review of the development of wave mechanics, we have established several basic properties of the quantum description of nature:. A particle
More informationLecture 11 Spin, orbital, and total angular momentum Mechanics. 1 Very brief background. 2 General properties of angular momentum operators
Lecture Spin, orbital, and total angular momentum 70.00 Mechanics Very brief background MATH-GA In 9, a famous experiment conducted by Otto Stern and Walther Gerlach, involving particles subject to a nonuniform
More informationPHY 407 QUANTUM MECHANICS Fall 05 Problem set 1 Due Sep
Problem set 1 Due Sep 15 2005 1. Let V be the set of all complex valued functions of a real variable θ, that are periodic with period 2π. That is u(θ + 2π) = u(θ), for all u V. (1) (i) Show that this V
More information5. Orthogonal matrices
L Vandenberghe EE133A (Spring 2017) 5 Orthogonal matrices matrices with orthonormal columns orthogonal matrices tall matrices with orthonormal columns complex matrices with orthonormal columns 5-1 Orthonormal
More information18.06 Problem Set 9 - Solutions Due Wednesday, 21 November 2007 at 4 pm in
8.6 Problem Set 9 - Solutions Due Wednesday, 2 November 27 at 4 pm in 2-6. Problem : (5) When A = SΛS is a real-symmetric (or Hermitian) matrix, its eigenvectors can be chosen orthonormal and hence S =
More information235 Final exam review questions
5 Final exam review questions Paul Hacking December 4, 0 () Let A be an n n matrix and T : R n R n, T (x) = Ax the linear transformation with matrix A. What does it mean to say that a vector v R n is an
More information1 Infinite-Dimensional Vector Spaces
Theoretical Physics Notes 4: Linear Operators In this installment of the notes, we move from linear operators in a finitedimensional vector space (which can be represented as matrices) to linear operators
More informationLecture 6: Vector Spaces II - Matrix Representations
1 Key points Lecture 6: Vector Spaces II - Matrix Representations Linear Operators Matrix representation of vectors and operators Hermite conjugate (adjoint) operator Hermitian operator (self-adjoint)
More informationProblems and Multiple Choice Questions
Problems and Multiple Choice Questions 1. A momentum operator in one dimension is 2. A position operator in 3 dimensions is 3. A kinetic energy operator in 1 dimension is 4. If two operator commute, a)
More informationLecture 4.6: Some special orthogonal functions
Lecture 4.6: Some special orthogonal functions Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics
More informationPrincipal Components Analysis
Principal Components Analysis Santiago Paternain, Aryan Mokhtari and Alejandro Ribeiro March 29, 2018 At this point we have already seen how the Discrete Fourier Transform and the Discrete Cosine Transform
More informationA Quantum Mechanical Model for the Vibration and Rotation of Molecules. Rigid Rotor
A Quantum Mechanical Model for the Vibration and Rotation of Molecules Harmonic Oscillator Rigid Rotor Degrees of Freedom Translation: quantum mechanical model is particle in box or free particle. A molecule
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationarxiv: v3 [physics.optics] 24 Jul 2013
Non-integer OAM beam shifts of Hermite-Laguerre-Gaussian beams arxiv:1303.495v3 [physics.optics] 4 Jul 013 Abstract A.M. Nugrowati, J.P. Woerdman Huygens Laboratory, Leiden University P.O. Box 9504, 300
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
More information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
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 informationarxiv: v1 [math-ph] 3 Nov 2011
Formalism of operators for Laguerre-Gauss modes A. L. F. da Silva (α), A. T. B. Celeste (β), M. Pazetti (γ), C. E. F. Lopes (δ) (α,β) Instituto Federal do Sertão Pernambucano, Petrolina - PE, Brazil (γ)
More information11.D.2. Collision Operators
11.D.. Collision Operators (11.94) can be written as + p t m h+ r +p h+ p = C + h + (11.96) where C + is the Boltzmann collision operator defined by [see (11.86a) for convention of notations] C + g(p)
More informationOutline. Introduction Definition of fractional fourier transform Linear canonical transform Implementation of FRFT/LCT
04//5 Outline Introduction Definition of fractional fourier transform Linear canonical transform Implementation of FRFT/LCT The Direct Computation DFT-like Method Chirp Convolution Method Discrete fractional
More informationTime part of the equation can be separated by substituting independent equation
Lecture 9 Schrödinger Equation in 3D and Angular Momentum Operator In this section we will construct 3D Schrödinger equation and we give some simple examples. In this course we will consider problems where
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More information1.6: 16, 20, 24, 27, 28
.6: 6, 2, 24, 27, 28 6) If A is positive definite, then A is positive definite. The proof of the above statement can easily be shown for the following 2 2 matrix, a b A = b c If that matrix is positive
More informationClosed-Form Discrete Fractional and Affine Fourier Transforms
1338 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 48, NO. 5, MAY 2000 Closed-Form Discrete Fractional and Affine Fourier Transforms Soo-Chang Pei, Fellow, IEEE, and Jian-Jiun Ding Abstract The discrete
More informationSchrödinger equation for central potentials
Chapter 2 Schrödinger equation for central potentials In this chapter we will extend the concepts and methods introduced in the previous chapter for a one-dimensional problem to a specific and very important
More informationECE 595, Section 10 Numerical Simulations Lecture 7: Optimization and Eigenvalues. Prof. Peter Bermel January 23, 2013
ECE 595, Section 10 Numerical Simulations Lecture 7: Optimization and Eigenvalues Prof. Peter Bermel January 23, 2013 Outline Recap from Friday Optimization Methods Brent s Method Golden Section Search
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationAngular Momentum. Classically the orbital angular momentum with respect to a fixed origin is. L = r p. = yp z. L x. zp y L y. = zp x. xpz L z.
Angular momentum is an important concept in quantum theory, necessary for analyzing motion in 3D as well as intrinsic properties such as spin Classically the orbital angular momentum with respect to a
More informationPY 351 Modern Physics - Lecture notes, 3
PY 351 Modern Physics - Lecture notes, 3 Copyright by Claudio Rebbi, Boston University, October 2016. These notes cannot be duplicated and distributed without explicit permission of the author. Time dependence
More informationKernel Principal Component Analysis
Kernel Principal Component Analysis Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr
More informationThe Schrodinger Equation and Postulates Common operators in QM: Potential Energy. Often depends on position operator: Kinetic Energy 1-D case:
The Schrodinger Equation and Postulates Common operators in QM: Potential Energy Often depends on position operator: Kinetic Energy 1-D case: 3-D case Time Total energy = Hamiltonian To find out about
More informationCollection of formulae Quantum mechanics. Basic Formulas Division of Material Science Hans Weber. Operators
Basic Formulas 17-1-1 Division of Material Science Hans Weer The de Broglie wave length λ = h p The Schrödinger equation Hψr,t = i h t ψr,t Stationary states Hψr,t = Eψr,t Collection of formulae Quantum
More informationPhysics 505 Homework No. 1 Solutions S1-1
Physics 505 Homework No s S- Some Preliminaries Assume A and B are Hermitian operators (a) Show that (AB) B A dx φ ABψ dx (A φ) Bψ dx (B (A φ)) ψ dx (B A φ) ψ End (b) Show that AB [A, B]/2+{A, B}/2 where
More informationMassachusetts Institute of Technology Physics Department
Massachusetts Institute of Technology Physics Department Physics 8.32 Fall 2006 Quantum Theory I October 9, 2006 Assignment 6 Due October 20, 2006 Announcements There will be a makeup lecture on Friday,
More informationSchrödinger equation for central potentials
Chapter 2 Schrödinger equation for central potentials In this chapter we will extend the concepts and methods introduced in the previous chapter ifor a one-dimenional problem to a specific and very important
More informationSEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationDiatomic Molecules. 7th May Hydrogen Molecule: Born-Oppenheimer Approximation
Diatomic Molecules 7th May 2009 1 Hydrogen Molecule: Born-Oppenheimer Approximation In this discussion, we consider the formulation of the Schrodinger equation for diatomic molecules; this can be extended
More informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
More informationThe Time{Dependent Born{Oppenheimer Approximation and Non{Adiabatic Transitions
The Time{Dependent Born{Oppenheimer Approximation and Non{Adiabatic Transitions George A. Hagedorn Department of Mathematics, and Center for Statistical Mechanics, Mathematical Physics, and Theoretical
More informationPolynomial Chaos and Karhunen-Loeve Expansion
Polynomial Chaos and Karhunen-Loeve Expansion 1) Random Variables Consider a system that is modeled by R = M(x, t, X) where X is a random variable. We are interested in determining the probability of the
More information(1) ... (2) /$ IEEE
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: REGULAR PAPERS, VOL 55, NO 11, DECEMBER 2008 3469 Closed-Form Orthogonal DFT Eigenvectors Generated by Complete Generalized Legendre Sequence Soo-Chang Pei,
More informationSummary: angular momentum derivation
Summary: angular momentum derivation L = r p L x = yp z zp y, etc. [x, p y ] = 0, etc. (-) (-) (-3) Angular momentum commutation relations [L x, L y ] = i hl z (-4) [L i, L j ] = i hɛ ijk L k (-5) Levi-Civita
More informationAddition of Angular Momenta
Addition of Angular Momenta What we have so far considered to be an exact solution for the many electron problem, should really be called exact non-relativistic solution. A relativistic treatment is needed
More informationFunctional Analysis. James Emery. Edit: 8/7/15
Functional Analysis James Emery Edit: 8/7/15 Contents 0.1 Green s functions in Ordinary Differential Equations...... 2 0.2 Integral Equations........................ 2 0.2.1 Fredholm Equations...................
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
More informationSupplementary information I Hilbert Space, Dirac Notation, and Matrix Mechanics. EE270 Fall 2017
Supplementary information I Hilbert Space, Dirac Notation, and Matrix Mechanics Properties of Vector Spaces Unit vectors ~xi form a basis which spans the space and which are orthonormal ( if i = j ~xi
More informationCourse Outline. Date Lecture Topic Reading
Course Outline Date Lecture Topic Reading Graduate Mathematical Physics Tue 24 Aug Linear Algebra: Theory 744 756 Vectors, bases and components Linear maps and dual vectors Inner products and adjoint operators
More informationmsqm 2011/8/14 21:35 page 189 #197
msqm 2011/8/14 21:35 page 189 #197 Bibliography Dirac, P. A. M., The Principles of Quantum Mechanics, 4th Edition, (Oxford University Press, London, 1958). Feynman, R. P. and A. P. Hibbs, Quantum Mechanics
More informationApproximation of Functions by Multivariable Hermite Basis: A Hybrid Method
Approximation of Functions by Multivariable Hermite Basis: A Hybrid Method Bartlomiej Beliczynski Warsaw University of Technology, Institute of Control and Industrial Electronics, ul. Koszykowa 75, -66
More informationarxiv: v1 [math-ph] 10 Sep 2013
arxiv:1309.2544v1 [math-ph] 10 Sep 2013 An Overview of the Relationship between Group Theory and Representation Theory to the Special Functions in Mathematical Physics Ryan D. Wasson Advisor: Robert Gilmore,
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationQuantum Mechanics: Vibration and Rotation of Molecules
Quantum Mechanics: Vibration and Rotation of Molecules 8th April 2008 I. 1-Dimensional Classical Harmonic Oscillator The classical picture for motion under a harmonic potential (mass attached to spring
More informationarxiv: v1 [quant-ph] 14 Dec 2009
Collocation method for fractional quantum mechanics Paolo Amore Facultad de Ciencias, CUICBAS,Universidad de Colima, Bernal Díaz del Castillo 340, Colima, Colima, Mexico Francisco M. Fernández INIFTA Conicet,
More informationNumerical Simulation of Spin Dynamics
Numerical Simulation of Spin Dynamics Marie Kubinova MATH 789R: Advanced Numerical Linear Algebra Methods with Applications November 18, 2014 Introduction Discretization in time Computing the subpropagators
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 informationIntroduction to Mathematical Physics
Introduction to Mathematical Physics Methods and Concepts Second Edition Chun Wa Wong Department of Physics and Astronomy University of California Los Angeles OXFORD UNIVERSITY PRESS Contents 1 Vectors
More informationLecture #1. Review. Postulates of quantum mechanics (1-3) Postulate 1
L1.P1 Lecture #1 Review Postulates of quantum mechanics (1-3) Postulate 1 The state of a system at any instant of time may be represented by a wave function which is continuous and differentiable. Specifically,
More informationCHEM 301: Homework assignment #5
CHEM 30: Homework assignment #5 Solutions. A point mass rotates in a circle with l =. Calculate the magnitude of its angular momentum and all possible projections of the angular momentum on the z-axis.
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 informationM.Sc. DEGREE EXAMINATION, DECEMBER First Year. Physics. Paper I MATHEMATICAL PHYSICS. Answer any FIVE questions.
(DPHY 01) M.Sc. DEGREE EXAMINATION, DECEMBER 009. Paper I MATHEMATICAL PHYSICS (5 0 = 100) 1. (a) Obtain the series solution for Legendre s differential equation. (b) Evaluate : J ( ) 1.. (a) Starting
More informationQuantum mechanics. Chapter The quantum mechanical formalism
Chapter 5 Quantum mechanics 5.1 The quantum mechanical formalism The realisation, by Heisenberg, that the position and momentum of a quantum mechanical particle cannot be measured simultaneously renders
More informationElliptic Fourier Transform
NCUR, ART GRIGORYAN Frequency Analysis of the signals on ellipses: Elliptic Fourier Transform Joao Donacien N. Nsingui Department of Electrical and Computer Engineering The University of Texas at San Antonio
More information