Notes on the Eigenfunction Method for solving differential equations

Similar documents
Chapter Five - Eigenvalues, Eigenfunctions, and All That

LECTURE 3. Orthogonal Functions. n X. It should be noted, however, that the vectors f i need not be orthogonal nor need they have unit length for

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Sturm-Liouville Theory

Homework Problem Set 1 Solutions

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

Sturm-Liouville Eigenvalue problem: Let p(x) > 0, q(x) 0, r(x) 0 in I = (a, b). Here we assume b > a. Let X C 2 1

Math 211A Homework. Edward Burkard. = tan (2x + z)

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

5.4, 6.1, 6.2 Handout. As we ve discussed, the integral is in some way the opposite of taking a derivative. The exact relationship

Problem Set 2 Solutions

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

Best Approximation in the 2-norm

Theoretical foundations of Gaussian quadrature

Green function and Eigenfunctions

g i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH

Math Theory of Partial Differential Equations Lecture 2-9: Sturm-Liouville eigenvalue problems (continued).

Abstract inner product spaces

Orthogonal Polynomials

Best Approximation. Chapter The General Case

The Dirac distribution

M597K: Solution to Homework Assignment 7

Example Sheet 2 Solutions

f a L Most reasonable functions are continuous, as seen in the following theorem:

Hilbert Spaces. Chapter Inner product spaces

Convex Sets and Functions

Problem Set 3 Solutions

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),

21.6 Green Functions for First Order Equations

x ) dx dx x sec x over the interval (, ).

Linearity, linear operators, and self adjoint eigenvalue problems

1 2-D Second Order Equations: Separation of Variables

Matrix & Vector Basic Linear Algebra & Calculus

CHAPTER 9 BASIC CONCEPTS OF DIFFERENTIAL AND INTEGRAL CALCULUS

2 Sturm Liouville Theory

Chapter 28. Fourier Series An Eigenvalue Problem.

STUDY GUIDE FOR BASIC EXAM

Conservation Law. Chapter Goal. 6.2 Theory

3 Mathematics of the Poisson Equation

4 Sturm-Liouville Boundary Value Problems

The Regulated and Riemann Integrals

Practice Problems Solution

Basic Derivative Properties

Do the one-dimensional kinetic energy and momentum operators commute? If not, what operator does their commutator represent?

NOTES ON HILBERT SPACE

Introduction and Review

DEFINITION The inner product of two functions f 1 and f 2 on an interval [a, b] is the number. ( f 1, f 2 ) b DEFINITION 11.1.

Summary: Method of Separation of Variables

1 E3102: a study guide and review, Version 1.0

Introduction to Complex Variables Class Notes Instructor: Louis Block

ENGI 9420 Lecture Notes 7 - Fourier Series Page 7.01

lim P(t a,b) = Differentiate (1) and use the definition of the probability current, j = i (

APPENDIX. Precalculus Review D.1. Real Numbers and the Real Number Line

20 MATHEMATICS POLYNOMIALS

Partial Differential Equations

( x) ( ) takes at the right end of each interval to approximate its value on that

If we have a function f(x) which is well-defined for some a x b, its integral over those two values is defined as

(4.1) D r v(t) ω(t, v(t))

Math Fall 2006 Sample problems for the final exam: Solutions

Course 2BA1 Supplement concerning Integration by Parts

Advanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration.

SUPPLEMENTARY NOTES ON THE CONNECTION FORMULAE FOR THE SEMICLASSICAL APPROXIMATION

1 The Lagrange interpolation formula

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

VII. The Integral. 50. Area under a Graph. y = f(x)

Homework Assignment 5 Solution Set

Improper Integrals, and Differential Equations

Continuous Random Variables

MATH 423 Linear Algebra II Lecture 28: Inner product spaces.

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx

Chapter 3 Polynomials

Linear Differential Equations Physics 129a Solutions to Problems Frank Porter Revision F. Porter

Analytical Methods Exam: Preparatory Exercises

MATH 174A: PROBLEM SET 5. Suggested Solution

Section 17.2 Line Integrals

STURM-LIOUVILLE PROBLEMS

Lecture Notes: Orthogonal Polynomials, Gaussian Quadrature, and Integral Equations

Lecture 1. Functional series. Pointwise and uniform convergence.

sec x over the interval (, ). x ) dx dx x 14. Use a graphing utility to generate some representative integral curves of the function Curve on 5

Fourier series. Preliminary material on inner products. Suppose V is vector space over C and (, )

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature

Math 61CM - Solutions to homework 9

Separation of Variables in Linear PDE

Math& 152 Section Integration by Parts

1 1D heat and wave equations on a finite interval

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

Variational Techniques for Sturm-Liouville Eigenvalue Problems

Inner-product spaces

Lecture 19: Continuous Least Squares Approximation

Calculus of Variations

Chapter 5 : Continuous Random Variables

Massachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6

Definite integral. Mathematics FRDIS MENDELU

Differential Equations 2 Homework 5 Solutions to the Assigned Exercises

INTRODUCTION TO INTEGRATION

Review of Gaussian Quadrature method

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Math Solutions to homework 1

Transcription:

Notes on the Eigenfunction Metho for solving ifferentil equtions Reminer: Wereconsieringtheinfinite-imensionlHilbertspceL 2 ([, b] of ll squre-integrble functions over the intervl [, b] (ie, b f(x 2 <1 with the inner prouct efine s n the norm inuce by this inner prouct hf gi : f (xg(x " kfk : p hf fi f(x 2 # 1/2 One cn introuce set of linerly inepenent bsis functions {y n } 1, where y n 2 L 2 ([, b] for n 0, 1,, such tht ny function f 2 L 2 ([, b] cn be represente s liner combintion of these functions f(x c n y n (x Such set is not unique n, in generl, the functions y n re not orthonorml, but they o form complete bsis over the intervl [, b] Furthermore, one cn prouce n orthonorml bsis for the Hilbert spce L 2 ([, b] from {y n } 1 using the Grm-Schmit proceure (ctully, this proceure cn be use to generte n orthonorml bsis for ny inner prouct spce 8 >< >: 0 y 0, ˆ0 0 k 0k 1 y 1 ˆ0h ˆ0 y 1 i, ˆ1 1 k 1k n y n P n 1 k0 ˆkh ˆk y n i, ˆk k k k k Using the set of orthonorml functions { ˆn} 1 one cn represent ny function f 2 L 2 ([, b] s f(x c n ˆn(x with the coefficients c n given by c n h ˆn fi This is true only if the bsis is orthonorml ˆ n (xf(x

1 Hermitin Conjugte of n Opertor Definition The joint of n opertor L,enotebyL,isefineby f (x[l g(x] [L f(x] g(x + bounry terms, (1 where the bounry terms (bt re evlute t the en-points of the intervl [, b] Using the efinition of n inner prouct one my re-write it s hf L gi hl f gi + bt For given liner ifferentil opertor, the metho for formlly fining the joint is to integrte by prts enough times to move the erivtives from g to f Exmple Compute the Hermitin conjugte of the following opertors: 1 L 1 hf gi pple ( f (x g(x u (f (xg(x b pple v 0 u f (x u 0 f (x f (xg(x f(x g(x +(f (xg(x b bt v 0 g(x v g(x 2 L 2 ı hf ı gi pple ( f (x ı g(x u v 0 u f (x u 0 f (x ı(f (xg(x b ı f (xg(x pple ı f(x g(x + ı(f (xg(x b bt v 0 ı g(x v ıg(x ı ı Remrk If L L then it is si tht L is self-joint Definition 1 If in Eq (1 the bt vnishes, then n opertor L is si to be Hermitin over the intervl

[, b] n f (x[l g(x] [L f(x] g(x (2 (equivlently, hf L gi hl f gi 2 Properties of Hermitin opertors From now on we re going to ssume tht the opertor L is Hermitin Relity of the eigenvlues Let y p n y q be two eigenfunctions of n opertor L corresponing to eigenvlues p n q,respectively L y p (x p y p (x L y q (x q y q (x + h y q i + hy p i hl y p y q i h p y p y q i hy p L y q i hy p qy q i hl y p y q i phy p y q i hl y p y q i q hy p y q i + L L hl y p y q i q hy p y q i Thus one hs phy p y q i qhy p y q i 0 p q hy p y q i 0 We get the sme result using the integrl form of the inner prouct Multiplying the eqution L y p (x p y p (x by y q n eqution L y q (x q y q (x by y p,nthenintegrtingovertheintervl[, b] les to The complex conjugte of Eq (3 becomes y q (x[l y p (x] p y p(x[l y q (x] q y q (x y p (x, (3 y p(x y q (x (4 y q (x[l y p (x] p y q (x y p(x Using now the efinition of Hermitin opertor (Eq 2 one gets [L y q (x] y p(x p y q (x y p(x (5

Compring Eq (4 n Eq (5 les to p y q (x y p(x q y p(x y q (x p q yp(x y q (x 0 Now, if p q, then it follows tht p p,ie, p is rel, since b y p(x y p (x hy p y p i n 0, 1,,functionsy n 60 An wht if p 6 q? 0 (for ll Orthogonlity n normliztion of the eigenfunctions We re going to consier two cses: 1 when y p n y q correspon to ifferent eigenvlues, p n q respectively, n 2 when y p n y q correspon to the sme eigenvlue 0 1 If y p n y q correspon to ifferent eigenvlues, then or, in terms of integrls, p q hy p y q i 0 hy p y q i 0 p q yp(x y q (x 0 y p(x y q (x 0 n we get two equivlent sttements of the orthogonlity of y p n y q 2 If 0 is egenerte eigenvlue, ie, L y`(x 0 y` for ` 0, 1,,k 1 then it my hppen tht y` re not orthogonl orthonormliztion proceure one cn obtin k orthogonl vectors However, by performing the Grm-Schmit To normlize the set of bsis functions {y n } 1 it is enough to ivie ech function by its norm ŷ n (x y n(x ky n (xk, obtining in this wy orthonorml set of bsis functions { yˆ n } 1 tht re stisfying the conition hŷ p ŷ q i 8 < ŷp(xŷ 0, if p 6 q q (x : 1, if p q Completeness of the eigenfunctions The infinite set of normlize eigenfunctions of Hermitin opertor form complete bsis set over [, b], ie,nyfunctiony 2 L 2 ([, b] obeying the bounry conitions my be expresse s liner combintion of these functions s y(x c n ŷ n (x,

where, s we know, c n hŷ n yi Thus y(x ŷ n (x On the other hn ŷ n(z y(zz ŷ n(x y(x y(z y(x y(z (x zz, ŷ n (xŷn(z where (x z stns for the Dirc elt function Compring these two equlities les to the completeness (or closure propertyoftheeigenfunctions ŷ n (xŷn(z (x z For continuous spectrum one hs to replce summtion by integrtion 3 Why is it importnt? ŷ n (xŷ n(z (x z We re concerne with the solution to the inhomogeneous eqution of the form L y(x f(x, (6 where f is some known function n the bounry conitions re given The ie for solving this problem is to use the linerity n Hemiticity of n opertor L Howoweoit? Step-1: Fin ll eigenfunctions of n opertor L,ie,solve the eqution L y n (x n y n (x Trnsform eigenfunctions {y n } 1 so tht they form n orthonorml bsis Step-2: Re-write both y n f s superpositions of eigenfunctions of L y(x f(x c n y n (x, f n y n (x where c n n f n stn for some constnt

Step-3: Mke use of the linerity of the opertor L! L y(x L c n y n (x c n L y n (x c n n y n (x f n y n (x f(x Step-4: Tke vntge of the orthonormlity of functions {y n } 1 -tkeninnerproucthy m i of both sies in the bove equlity Since hy m y n i mn,thus c n n hy m y n i f n hy m y n i c m m f m c m f m Step-5: Put your coefficients in to get generl solution of the problem (6 y(x 4 The Sturm-Liouville equtions Let us consier the following eqution p(x 2 y(x 2 + r(x y(x f n n y n (x m + q(xy(x+ (xy(x 0, where r(x p(x, p, q n r re rel functions n (x > 0 is clle weight or ensity function It is cler tht this eqution cn be written s L y(x (xy(x, where L (L is clle the Sturm-Liouville opertor pplep(x 2 2 + r(x + q(x pple p(x + q(x It turns out the with n pproprite bounry conitions the Sturm-Liouville opertor is Hermitin over the intervl [, b] These bounry conitions re [y i (xp(xy 0 j(x] x [y i (xp(xy 0 j(x] xb for ll i, j, where y i n y j stn for two eigenfunctions of the opertor L To check it let us put the Sturm-Liouville

form L (py 0 0 qy into the Eq (2 y i (x[l y j (x] ( ( yi (py 0 j 0 + qy j u yi v 0 (pyj 0 0 u 0 (y i 0 y i pyj 0 b + vnishes u (y i 0 p u 0 ((y i 0 p 0 [(y i 0 py j ] b vnishes v py 0 j (y i 0 py 0 j v 0 y 0 j v y j ((y i 0 p 0 y j [(p(y i 0 0 + qy i ] y j y i (py 0 j 0 y i qy j y i qy j [L y i (x] y j (x y i qy j [((y i 0 p 0 y j + y i qy j ] Why is the Sturm-Liouville eqution importnt? It is becuse ll secon-orer liner orinry ifferentil equtions cn be recst in the form of the LHS of the Sturm-Liouville eqution How? Consier n orinry ifferentil eqution of the form We wnt to rewrite this eqution s P (x 2 y(x 2 + Q(x y(x + R(xy(x 0 pple p(x + q(x y(x (xy(x To o so, one must fin p, q n in the following wy: Exmple 2 Let us consier Bessel s eqution p(x x Q(t e P (t t, (x q(x p(x R(x P (x x 2 y 00 + xy 0 +( 2 x 2 v 2 y 0 Here x t p(x e t 2 t e ln x x (x q(x x ( 2 x 2 v 2 x 2 ( 2 x 2 v 2 x Thus we get (xy 0 0 +( x v 2 y 0 x