Sturm-Liouville Theory

Similar documents
Sturm-Liouville Theory

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

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

Notes on the Eigenfunction Method for solving differential equations

Green function and Eigenfunctions

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

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.

Conservation Law. Chapter Goal. 6.2 Theory

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

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

Variational Techniques for Sturm-Liouville Eigenvalue Problems

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

Introduction to Complex Variables Class Notes Instructor: Louis Block

Homework Problem Set 1 Solutions

AMS 212A Applied Mathematical Methods I Lecture 06 Copyright by Hongyun Wang, UCSC. ( ), v (that is, 1 ( ) L i

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH

Theoretical foundations of Gaussian quadrature

Section 6.1 INTRO to LAPLACE TRANSFORMS

Best Approximation in the 2-norm

Chapter 28. Fourier Series An Eigenvalue Problem.

Linearity, linear operators, and self adjoint eigenvalue problems

Problem Set 2 Solutions

Fundamental Theorem of Calculus

MATH 174A: PROBLEM SET 5. Suggested Solution

ax bx c (2) x a x a x a 1! 2!! gives a useful way of approximating a function near to some specific point x a, giving a power-series expansion in x

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O

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

Fourier Series in Complex notation. cos(x) = eix + e ix 2. A n cos + B n sin l. i 2 B n. e inx=l + A n + ib n 2 8. ( 0 m 6= n. C n = 1 2l.

REGULAR TWO-POINT BOUNDARY VALUE PROBLEMS

M344 - ADVANCED ENGINEERING MATHEMATICS

c n φ n (x), 0 < x < L, (1) n=1

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

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

Review of Gaussian Quadrature method

MATH1050 Cauchy-Schwarz Inequality and Triangle Inequality

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

The Regulated and Riemann Integrals

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

Chapter 6 Techniques of Integration

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

dx dt dy = G(t, x, y), dt where the functions are defined on I Ω, and are locally Lipschitz w.r.t. variable (x, y) Ω.

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

Coalgebra, Lecture 15: Equations for Deterministic Automata

STURM-LIOUVILLE PROBLEMS

x dx does exist, what does the answer look like? What does the answer to

The practical version

1 The Riemann Integral

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

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

Abstract inner product spaces

LINEAR ALGEBRA APPLIED

Lecture 2e Orthogonal Complement (pages )

Hilbert Spaces. Chapter Inner product spaces

Introduction and Review

(9) P (x)u + Q(x)u + R(x)u =0

Functional Analysis I Solutions to Exercises. James C. Robinson

Partial Differential Equations

Quadratic Forms. Quadratic Forms

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

Chapter 8: Methods of Integration

1 E3102: a study guide and review, Version 1.0

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

Improper Integrals. The First Fundamental Theorem of Calculus, as we ve discussed in class, goes as follows:

September 13 Homework Solutions

Wave Equation on a Two Dimensional Rectangle

Best Approximation. Chapter The General Case

Section 3.2 Maximum Principle and Uniqueness

Matrices, Moments and Quadrature, cont d

Lecture Solution of a System of Linear Equation

The Evaluation Theorem

Practice Problems Solution

Construction of Gauss Quadrature Rules

1 1D heat and wave equations on a finite interval

Math Lecture 23

1 2-D Second Order Equations: Separation of Variables

Summary: Method of Separation of Variables

Example Sheet 2 Solutions

f(x)dx . Show that there 1, 0 < x 1 does not exist a differentiable function g : [ 1, 1] R such that g (x) = f(x) for all

Inner-product spaces

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

Matrix & Vector Basic Linear Algebra & Calculus

Applied Partial Differential Equations with Fourier Series and Boundary Value Problems 5th Edition Richard Haberman

k and v = v 1 j + u 3 i + v 2

The Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

Farey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University

Things to Memorize: A Partial List. January 27, 2017

38 Riemann sums and existence of the definite integral.

Linear Systems with Constant Coefficients

Stuff You Need to Know From Calculus

ODE: Existence and Uniqueness of a Solution

Chapter 3 Polynomials

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

SUPPLEMENTARY NOTES ON THE CONNECTION FORMULAE FOR THE SEMICLASSICAL APPROXIMATION

The solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr

Matrix Eigenvalues and Eigenvectors September 13, 2017

1. Review. t 2. t 1. v w = vw cos( ) where is the angle between v and w. The above leads to the Schwarz inequality: v w vw.

Math 61CM - Solutions to homework 9

Transcription:

LECTURE 1 Sturm-Liouville Theory In the two preceing lectures I emonstrte the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series re just the tip of the iceerg of the theory n utility of specil functions. Before preceing with the generl theory, let me stte clerly the sic properties of Fourier series we inten to generlize: i The Fourier sine functions { sin nπ L x n = 1,, 3,... } n cosine functions { cos } nπ L x n =, 1,,... re solutions of secon orer liner homogeneous ifferentil eqution 1 y = λy stisfying certin liner homogenous ounry conitions: viz, the Fourier-sine functions stisfy n the Fourier-cosine functions stisfy y = = y L y = = y L. ii The Fourier sine n cosine functions oey certin orthogonlity reltions L nπ mπ { 1 if n = m sin L L x sin L x = if n m L nπ mπ { 1 if n = m cos L L x cos L x = if n m iii Any continuous, ifferentile function on, L] cn e expresse in terms of Fourier-sine or Fourier-cosine expnsion f x = + nπ n cos L x = n=1 nπ n sin L x n=1 Moreover, the orthogonlity reltions ii llow us to etermine the coefficients n n n n = L nπ f x cos L L x n = L nπ f x sin L L x 1. Sturm-Liouville Prolems Definition 1.1. A Sturm-Liouville prolem is secon orer homogeneous liner ifferentil eqution of the form 3 p x y ] q x y + λr x y =, x, ] 55

1. STURM-LIOUVILLE PROBLEMS 56 together with ounry conitions of the form 4 4 y + α y = y + β y = In wht follow, λ is to e regre s constnt prmeter n we shll lwys ssume tht p x n r x re positive functions on the intervl, ] 5 p x > x, ] r x > x, ] Exmple 1.. Tking p x = 1 r x = 1 q x = = = L = = 1 α = β =, we see tht the corresponing Sturm-Liouville prolem y = λy y = y L = hs s its solutions nπ y x = sin L x nπ n λ = n = 1,, 3,... L Note, in prticulr, tht the solution of the Sturm-Liouville prolem only exists for certn vlues of λ; these vlues re eigenvlues of the Sturm-Liouville prolem Nottion 1.3 Liner ifferentil opertor nottion. In wht follows we shll enote y L the liner ifferentil opertor 6 L = p x ] q x which cts on function y x y L y] = p x y ] q x y In terms of L, the ifferentil eqution of Sturm-Liouville prolem cn e expresse L y] = λr x y Theorem 1.4 Lgrnge s ientity. Suppose φ n ψ re two functions stisfying ounry conitions of the form 7 7 then y + α y = y + β y = L φ] ψ φl ψ] =

1. STURM-LIOUVILLE PROBLEMS 57 Proof. We hve * L φ] ψ = = = p φ ψ = p φ ψ = p φ ψ p φ p φ ] ψ qφψ ] ψ p φ φ φp ψ p ψ ψ + qφψ φ = p φ ψ φψ + φl ψ] v where we hve twice utilize the integrtion y prts formul Now consier the first term in the lst line f g = fg qφψ qφψ p ψ f g qφψ p φ ψ φψ Since φ n ψ stisfy 7 n 7 we hve ssuming, α,, β n Hence φ = α φ, ψ = α ψ φ = β φ, ψ = β ψ φ ψ φ ψ = α φ ψ φ α φ ψ φ ψ = β φ ψ φ ψ β ψ p u v uv 1 = p 1 u 1 v 1 u 1 v 1 p u v u v = p 1 p = We cn hence conclue from * tht n the theorem follows. L φ] ψ = + φl ψ] = = Recll tht, for exmple, the Fourier-sine functions sin λx not only stisfy certin secon orer liner ifferentil eqution f = λ f ut they lso stisfy certin liner homogeneous ounry conitions f = = f L

1. STURM-LIOUVILLE PROBLEMS 58 exctly when the prmeter λ is tune to the ounry conitions { = f = λ f f = f L = λ = nπ L f x = sin nπ L x Wht my seem little surprising t first is tht the fct tht the Fourier-sine functions re solutions to Sturm-Liouville prolem is lso responsile for their orthogonlity properties L nπ mπ { 1 if n = m sin L L x sin L x = if n m Theorem 1.5. If φ n ψ re two solutions of Sturm-Liouville prolem corresponing to ifferent vlues of the prmeter λ then Proof. Suppose φ x ψ x r x =. L φ] = λ 1 r x φ L ψ] = λ r x ψ with λ 1 λ. Then since φ, ψ stisfy the Strum-Liouville ounry conitions, we hve y Theorem 11.1 L φ] ψ = φ L ψ] n ecuse they stisfy Sturm-Liouville-type ifferentil equtions or λ 1 r x φ x ψ x = λ 1 λ By hypothesis, λ 1 λ n so we must conclue φ x λ r x ψ x φ x ψ x r x = φ x ψ x r x = As remrke ove, it is common to write in shorthn the ifferentil eqution of Sturm-Liouville prolem s 8 L y] = λr x y n refer to the prmeter λ corresponing to prticulr solution y x s n eigenvlue of the ifferentil opertor L. This lnguge, of course, is erive from the nomenclture of liner lger where if A is n n n mtrix, v is n n 1 column vector n λ is numer such tht Av = λv then v is clle n eigenvector of A n λ is the corresponing eigenvlue of A. In fct, it is common to use such liner-lgeric-like nomenclture throughout Sturm-Liouville theory; n henceforth we shll refer to the functions y x tht stisfy 8 s eigenfunctions of L n the numers λ s the corresponing eigenvlues of L. Theorem 1.6. All the eigenvlues of Sturm-Liouville prolem re rel.

1. STURM-LIOUVILLE PROBLEMS 59 Proof. Consier the following inner prouct on the spce of complex-vlue solutions of Sturm-Liouville prolem 8 for vrious λ Suppose u, v = u x v x u x = v x = R x + ii x tht is, suppose we set u equl to v n split it up into its rel n imginry prts. Then u, u = = R x + ii x R x ii x R x + I x ] Note now tht the integrn is lwys non-negtive. By well-known result from Clculus, the integrl of continuous non-negtive function over finite intervl is lwys non-negtive, n zero only if the function is ienticlly equl to zero. We conclue tht for ny continuous, non-zero, complex-vlue function u x u x u x > Now consier eqution otine from the conclusion of Theorem 11.4 y setting v = u, 9 L u] v = n suppose u is Sturm-Liouville eigenfunction: ul u] 1 L u] = λru for some λ C Tking the complex conjugte of this eqution, we hve, ecuse p x, q x, n r x re ll ssume to e rel-vlue functions, or L u] = λrv = = 11 L u] = λru = L u] = λru p x u ] q x u = λru p x u ] q x u = λru An so plugging the right hn sies of 1 n 11 into 9 we cn conclue tht or 1 λ ru u = λu ru λ λ ruu = We hve rgue ove tht u x u x is non-negtive function, n y ssumption r x is positive function on, 1]. Thus, the integrn is non-negtive function n < ruu so long s u x is not ienticlly zero. We cn therefore conclue from 1 tht λ λ = = λ R.

1. STURM-LIOUVILLE PROBLEMS 6 In Liner Alger, we sy tht n eigenvlue λ of mtrix A hs multiplicity m if the imension of the corresponing eigenspce is m; tht is to sy, m = im NullSp A λi In Sturm-Liouville theory, we sy tht the multiplicity of n eigenvlue λ of Sturm-Liouville prolem 1 φ + φ = 1 φ 1 + φ 1 = L φ] = λr x φ x if there re exctly m linerly inepenent solutions for tht vlue of λ. Theorem 1.7. The eigenvlues of Sturm-Liouville prolem re ll of multiplicity one. Moreover, the eigenvlues form n infinite sequence n cn e orere ccoring to incresing mgnitue: {eigenvlues} = {λ 1, λ, λ 3,...}, λ 1 < λ < λ 3 < The orering of the eigenvlues λ, n the fct tht they multiplicity free, estlishes certin cnonicl orering of the corresponing Sturm-Liouville eigenfunctions. We enote y φ n the solution of the Sturm- Liouville prolem with λ = λ n. Actully, this only etermines φ n x up to constnt multiple ecuse if φ x stisfies 3, 4, 4 then so oes ny constnt multiple of φ x. It is common prctise to remove this miguity y emning in ition which fixes φ x up to sclr fctor of the form e iα. φ n x φ n xr x = 1 Moreover, the fcts tht the eigenvlues re ll rel n multiplicity free, lso implies tht we cn choose the φ n x to e rel-vlue functions of x ecuse the complex spn of the functions φ n x must coincie with tht of φ n x. Thus, we cn choose the eigenfunctions of Sturm-Liouville prolem to e rel-vlue functions φ n x φ n xr x = 1 Theorem 1.8. Let φ 1, φ, φ 3,... e the normlize eigenfunctions of Sturm-Liouville prolem n suppose f is piecewise continuous function on, 1]. Then if the series converges to t ech point on, 1. c n := f x φ n x r x c n φ n x n=1 f x + + f x Here f x ± := lim f x ± ε ε + Note tht t ny point x where f is continuous f x + + f x = f x.

1. STURM-LIOUVILLE PROBLEMS 61 In summry, ny time you hve ifferentil eqution of the form p x y ] q x y = λr x y with homogeneous ounry conitions of the form Then: 1 y + y = = 1 y 1 + y 1 Solutions will exist for only iscrete ut otherwise infinite set of vlues for the prmeter λ, ll of which re rel numers. There is unique lowest eigenvlue λ n the other eigenvlues cn e totlly orere λ < λ 1 < λ < λ 3 < For ech λ n, there is solution φ n of p x y ] q x y = λ n r x y tht is unique up to sclr fctor. The solutions φ n, n =, 1,,... cn e normlize such tht { 1 if n = m φ n x φ m x r x = otherwise Any continuous function f x on the intervl, 1] cn e expne in terms of the Sturm-Liouville eigenfunctions φ n, n =, 1,,,... f x = n φ n x n= with the coefficients n eing etermine y n = f x φ n x r x If f x is merely piece-wise continuous on, 1], one hs inste f x + + f x = n φ n x n= with the coefficients n etermine y the sme formul.