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

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

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

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

4 Sturm-Liouville Boundary Value Problems

Best Approximation in the 2-norm

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

1 2-D Second Order Equations: Separation of Variables

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

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

MATH 174A: PROBLEM SET 5. Suggested Solution

Green function and Eigenfunctions

Math Fall 2006 Sample problems for the final exam: Solutions

Chapter 28. Fourier Series An Eigenvalue Problem.

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

Linearity, linear operators, and self adjoint eigenvalue problems

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

Numerical Methods I Orthogonal Polynomials

Sturm-Liouville Theory

Hilbert Spaces. Chapter Inner product spaces

Variational Techniques for Sturm-Liouville Eigenvalue Problems

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

Theoretical foundations of Gaussian quadrature

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.

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

1 E3102: a study guide and review, Version 1.0

Math 270A: Numerical Linear Algebra

18 Sturm-Liouville Eigenvalue Problems

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

2 Sturm Liouville Theory

Chapter 3 Polynomials

Abstract inner product spaces

Notes on the Eigenfunction Method for solving differential equations

df dt f () b f () a dt

REGULAR TWO-POINT BOUNDARY VALUE PROBLEMS

Handout 4. Inverse and Implicit Function Theorems.

Eigenfunction Expansions for a Sturm Liouville Problem on Time Scales

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

1 E3102: A study guide and review, Version 1.2

NOTES ON HILBERT SPACE

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

Lecture 19: Continuous Least Squares Approximation

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

STURM-LIOUVILLE PROBLEMS

Functional Analysis I Solutions to Exercises. James C. Robinson

Chapter 3. Vector Spaces

Summary: Method of Separation of Variables

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 2

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

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

STURM-LIOUVILLE PROBLEMS: GENERALIZED FOURIER SERIES

MATRICES AND VECTORS SPACE

21.6 Green Functions for First Order Equations

Review SOLUTIONS: Exam 2

Matrix Eigenvalues and Eigenvectors September 13, 2017

ODE: Existence and Uniqueness of a Solution

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

Partial Differential Equations

Multivariate problems and matrix algebra

Best Approximation. Chapter The General Case

STUDY GUIDE FOR BASIC EXAM

Orthogonal Polynomials

Mapping the delta function and other Radon measures

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

Matrices, Moments and Quadrature, cont d

Lecture 1. Functional series. Pointwise and uniform convergence.

CAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.

Outline. Math Partial Differential Equations. Rayleigh Quotient. Rayleigh Quotient. Sturm-Liouville Problems Part C

International Jour. of Diff. Eq. and Appl., 3, N1, (2001),

Math 554 Integration

Chapter Five - Eigenvalues, Eigenfunctions, and All That

M597K: Solution to Homework Assignment 7

The Regulated and Riemann Integrals

REPRESENTATION THEORY OF PSL 2 (q)

Problem Set 4: Solutions Math 201A: Fall 2016

Advanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015

AMATH 731: Applied Functional Analysis Fall Additional notes on Fréchet derivatives

ECON 331 Lecture Notes: Ch 4 and Ch 5

Sturm Liouville Problems

Lecture 4: Piecewise Cubic Interpolation

Separation of Variables in Linear PDE

MA Handout 2: Notation and Background Concepts from Analysis

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 ),

The Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5

Problem set 5: Solutions Math 207B, Winter r(x)u(x)v(x) dx.

Calculus of Variations

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

Convex Sets and Functions

The Banach algebra of functions of bounded variation and the pointwise Helly selection theorem

Lecture 14: Quadrature

Math 118: Honours Calculus II Winter, 2005 List of Theorems. L(P, f) U(Q, f). f exists for each ǫ > 0 there exists a partition P of [a, b] such that

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

MTH 5102 Linear Algebra Practice Exam 1 - Solutions Feb. 9, 2016

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

Math 115 ( ) Yum-Tong Siu 1. Lagrange Multipliers and Variational Problems with Constraints. F (x,y,y )dx

Math 6455 Oct 10, Differential Geometry I Fall 2006, Georgia Tech

M344 - ADVANCED ENGINEERING MATHEMATICS

Orthogonal functions

Chapter 2. Vectors. 2.1 Vectors Scalars and Vectors

Transcription:

Mth 412-501 Theory of Prtil Differentil Equtions Lecture 2-9: Sturm-Liouville eigenvlue problems (continued).

Regulr Sturm-Liouville eigenvlue problem: d ( p dφ ) + qφ + λσφ = 0 ( < x < b), dx dx β 1 φ() + β 2 φ () = 0, β 3 φ(b) + β 4 φ (b) = 0. Here β i R, β 1 + β 2 0, β 3 + β 4 0. Functions p, q, σ re continuous on [, b], p > 0 nd σ > 0 on [, b].

6 properties of regulr Sturm-Liouville problem Eigenvlues re rel. Eigenvlues form n incresing sequence. n-th eigenfunction hs n 1 zeros in (, b). Eigenfunctions re orthogonl with weight σ. Eigenfunctions nd eigenvlues re relted through the Ryleigh quotient. Piecewise smooth functions cn be expnded into generlized Fourier series of eigenfunctions.

Regulr Sturm-Liouville eqution: d ( p dφ ) + qφ + λσφ = 0 ( < x < b). dx dx Consider liner differentil opertor L(f ) = d ( p df ) + qf. dx dx Now the eqution cn be rewritten s L(φ) + λσφ = 0.

Lgrnge s identity: gl(f ) f L(g) = d dx ( ) p(gf fg ) Integrting over [, b], we obtin Green s formul: b ( ) gl(f ) f L(g) dx = p(gf fg ) b Clim If f nd g stisfy the sme regulr boundry conditions, then the right-hnd side in Green s formul vnishes.

Suppose φ n nd φ m re eigenfunctions of the Sturm-Liouville problem corresponding to eigenvlues λ n nd λ m : L(φ n ) + λ n σφ n = 0, L(φ m ) + λ m σφ m = 0. Since φ n nd φ m stisfy the sme regulr boundry conditions, Green s formul implies tht b ( ) φ m L(φ n ) φ n L(φ m ) dx = 0 = b If λ n λ m, then (λ m λ n )φ n (x)φ m (x)σ(x) dx = 0 b φ n (x)φ m (x)σ(x) dx = 0.

Suppose φ is complex-vlued eigenfunction corresponding to complex eigenvlue λ: L(φ) + λσφ = 0, β 1 φ() + β 2 φ () = 0, β 3 φ(b) + β 4 φ (b) = 0. We re going to show tht λ R. Any complex number z = x + iy is ssigned its complex conjugte z = x iy. Let us pply the complex conjugcy to the Sturm-liouville eqution nd the boundry conditions.

L(φ) + λσφ = 0, β 1 φ() + β 2 φ () = β 3 φ(b) + β 4 φ (b) = 0. It is known tht z 1 + z 2 = z 1 + z 2 nd z 1 z 2 = z 1 z 2. L(φ) + λ σ φ = 0, β 1 φ()+β 2 φ () = β 3 φ(b)+β 4 φ (b) = 0. If z is rel then z = z. L(φ) + λσφ = 0, β 1 φ()+β 2 φ () = β 3 φ(b)+β 4 φ (b) = 0.

Let φ denote the complex conjugte function of φ, i.e., φ(x) = φ(x) for x b. We hve tht φ = f + ig, where f nd g re rel-vlued functions. Then φ = f ig. Note tht φ = (f ig) = f ig = f + ig = φ. It follows tht L(φ) = (pφ ) + qφ = (pφ ) + qφ = ( pφ ) + qφ = ( pφ ) + qφ = L(φ).

L ( φ ) + λσφ = 0, β 1 φ() + β 2 φ () = β 3 φ(b) + β 4 φ (b) = 0. If φ is n eigenfunction belonging to n eigenvlue λ, then φ is n eigenfunction belonging to the eigenvlue λ. Assume tht λ λ. Then But b b φ(x)φ(x)σ(x) dx = 0. φ(x)φ(x)σ(x) dx = Thus λ = λ = λ R. b φ(x) 2 σ(x) dx > 0.

Some fcts bout Eucliden spce Eucliden spce R 3. Let v = (v 1, v 2, v 3 ), u = (u 1, u 2, u 3 ) be two vectors. v u = v 1 u 1 + v 2 u 2 + v 3 u 3 is the dot product. v nd u re orthogonl if v u = 0. v = v v. Vectors e 1 = (1, 0, 0), e 2 = (0, 1, 0), e 3 = (0, 0, 1) form n orthonorml bsis. v = v 1 e 1 + v 2 e 2 + v 3 e 3 = (v e 1 )e 1 + (v e 2 )e 2 + (v e 3 )e 3.

Let v 1,v 2,v 3 be orthogonl nonzero vectors. They form bsis in R 3 so tht for ny u R 3 we hve u = c 1 v 1 + c 2 v 2 + c 3 v 3. Note tht u v n = c n v n v n so tht c n = u v n v n v n. Pythgoren theorem implies tht u 2 = c 1 v 1 2 + c 2 v 2 2 + c 3 v 3 2. Observe tht c n v n 2 = c n 2 v n v n. Hence u u = u v 1 2 v 1 v 1 + u v 2 2 v 2 v 2 + u v 3 2 v 3 v 3 (Prsevl s equlity)

Let v 1,v 2 be orthogonl nonzero vectors. Given vector u R 3, let u 0 = u (c 1 v 1 + c 2 v 2 ), where c n = u v n v n v n. It is esy to check tht u 0 v 1 = u 0 v 2 = 0 so tht u 0 (u u 0 ) = 0. Pythgoren theorem implies tht u 2 = c 1 v 1 2 + c 2 v 2 2 + u 0 2 c 1 v 1 2 + c 2 v 2 2. Since c n v n 2 = c n 2 v n v n, we get u u u v 1 2 v 1 v 1 + u v 2 2 v 2 v 2 (Bessel s inequlity)

Suppose A nd B re liner opertors in R 3. We sy tht B is djoint to A (denoted B = A ) if Au v = u Bv for ll u,v R 3. Let A = ( ij ) 1 i,j 3, B = (b ij ) 1 i,j 3. Then Ae j = 1j e 1 + 2j e 2 + 3j e 3, hence ij = Ae j e i. Similrly, b ij = Be j e i = e i Be j. It follows tht ij = b ji, i.e., B is the trnspose of A. A is clled self-djoint if A = A. Self-djoint opertors hve only rel eigenvlues. Suppose v 1, v 2 re eigenvectors of A belonging to eigenvlues λ 1, λ 2. Then λ 1 v 1 v 2 = Av 1 v 2 = v 1 Av 2 = λ 2 v 1 v 2. If λ 1 λ 2 then v 1 v 2 = 0.

From Eucliden spce to Hilbert spce Hilbert spce is n infinite-dimensionl nlogue of Eucliden spce. One reliztion is L 2 [, b] = {f : b f (x) 2 dx < }. Inner product of functions: f, g = b f (x)g(x) dx. Since fg 1 2 ( f 2 + g 2 ), the inner product is well defined for ny f, g L 2 [, b]. Norm of function: f = f, f. Convergence: we sy tht f n f in the men if f f n 0 s n.

Functions f, g L 2 [, b] re clled orthogonl if f, g = 0. Alterntive inner product: f, g w = b where w is the weight function. f (x)g(x)w(x) dx, Functions f nd g re clled orthogonl with weight w if f, g w = 0.

A set f 1, f 2,... of pirwise orthogonl nonzero functions is clled complete if it is mximl, i.e., there is no nonzero function g such tht g, f n = 0, n = 1, 2,.... A complete set forms bsis of the Hilbert spce, tht is, ech function g L 2 [, b] cn be expnded into series g = n=1 c nf n tht converges in the men. The expnsion is unique: c n = g, f n f n, f n.