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

Similar documents
Sturm-Liouville Theory

Notes on the Eigenfunction Method for solving differential equations

Chapter Five - Eigenvalues, Eigenfunctions, and All That

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

Theoretical foundations of Gaussian quadrature

Abstract inner product spaces

Best Approximation in the 2-norm

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

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

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.

Chapter 14. Matrix Representations of Linear Transformations

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

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

The Regulated and Riemann Integrals

Lecture 19: Continuous Least Squares Approximation

Best Approximation. Chapter The General Case

ODE: Existence and Uniqueness of a Solution

Chapter 3 Polynomials

The final exam will take place on Friday May 11th from 8am 11am in Evans room 60.

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.

Orthogonal Polynomials

Inner-product spaces

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

MATRICES AND VECTORS SPACE

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

Conservation Law. Chapter Goal. 6.2 Theory

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

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

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH

4. Calculus of Variations

AM1 Mathematical Analysis 1 Oct Feb Exercises Lecture 3. sin(x + h) sin x h cos(x + h) cos x h

ENGI 9420 Lecture Notes 7 - Fourier Series Page 7.01

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

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

Homework Assignment 5 Solution Set

Chapter 28. Fourier Series An Eigenvalue Problem.

Math& 152 Section Integration by Parts

Review SOLUTIONS: Exam 2

Hilbert Spaces. Chapter Inner product spaces

B Veitch. Calculus I Study Guide

20 MATHEMATICS POLYNOMIALS

Matrix & Vector Basic Linear Algebra & Calculus

1 2-D Second Order Equations: Separation of Variables

Green function and Eigenfunctions

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

8 Laplace s Method and Local Limit Theorems

Lecture 1. Functional series. Pointwise and uniform convergence.

STUDY GUIDE FOR BASIC EXAM

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

AP Calculus Multiple Choice: BC Edition Solutions

Introduction to Complex Variables Class Notes Instructor: Louis Block

Review of Calculus, cont d

Lecture 4: Piecewise Cubic Interpolation

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

Analytical Methods Exam: Preparatory Exercises

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

The Dirac distribution

REPRESENTATION THEORY OF PSL 2 (q)

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

DEFINITION OF ASSOCIATIVE OR DIRECT PRODUCT AND ROTATION OF VECTORS

ODE: Existence and Uniqueness of a Solution

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

Homework Problem Set 1 Solutions

STURM-LIOUVILLE PROBLEMS

Math Solutions to homework 1

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

Necessary and sufficient conditions for some two variable orthogonal designs in order 44

MATH 409 Advanced Calculus I Lecture 19: Riemann sums. Properties of integrals.

Stuff You Need to Know From Calculus

Elementary Linear Algebra

Basic Derivative Properties

set is not closed under matrix [ multiplication, ] and does not form a group.

1 Parametric Bessel Equation and Bessel-Fourier Series

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

Bases for Vector Spaces

Practice Problems Solution

Lecture 14: Quadrature

1 The Riemann Integral

Review of Riemann Integral

MATH 174A: PROBLEM SET 5. Suggested Solution

1 E3102: a study guide and review, Version 1.0

Summary: Method of Separation of Variables

Sections 5.2: The Definite Integral

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

Math 61CM - Solutions to homework 9

Math 554 Integration

Review of basic calculus

Mapping the delta function and other Radon measures

1 Linear Least Squares

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Lecture 3. Limits of Functions and Continuity

Linearity, linear operators, and self adjoint eigenvalue problems

18.06 Problem Set 4 Due Wednesday, Oct. 11, 2006 at 4:00 p.m. in 2-106

Calculus in R. Chapter Di erentiation

NOTES ON HILBERT SPACE

1 1D heat and wave equations on a finite interval

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

4 Sturm-Liouville Boundary Value Problems

Bernoulli Numbers Jeff Morton

Construction of Gauss Quadrature Rules

Transcription:

ECTURE 3 Orthogonl Functions 1. Orthogonl Bses The pproprite setting for our iscussion of orthogonl functions is tht of liner lgebr. So let me recll some relevnt fcts bout nite imensionl vector spces. Every vector in R n cn be represente s sum of the form 3.1 v = v i e i i=1 where e i is the unit vector long the i th coorinte xis. However if ff i g is ny set of n linerly inepenent vectors then we lso hve unique representtion of v s 3.2 v = c i f i : i=1 It shoul be note however tht the vectors f i nee not be orthogonl nor nee they hve unit length for the expnsion 3.2 to work. If however 3.3 f i f j = ij ; tht is if the f i re orthonorml then the coecients c i cn esily be compute. For if we tke the ot prouct of v with bsis vector f i we get tht is to sy f i v = f i @ X n c j f j 1 A = n X c j f i f j = c j ij = c i ; 3.4 c i = f i v : 2. Bses of Orthogonl Functions The relevnce of these remrks now comes from the observtion tht the set CëR; Rë of continous rel-vlue functions on the rel line is lso vector spce; for the opertions of ition n sclr multipliction of functions re well-ene: ëf + gëx = fx + gx ëc fëx = cfx 8 c 2 R: Note however tht CëR; Rë is n innite imensionl vector spce. Inee the Tylor expnsion 2.16 of cn be thought of s n expnsion of with respect to bsis of monomil functions ff mn x; t =x m t n g. Unfortuntely this bsis is not orthonorml t lest not with respect to ny obvious inner prouct. Fourier series re better exmples since they constitute expnsions of functions with respect to n orthonorml bsis of CëR;Rë. 1

2. BASES OF ORTHOGONA FUNCTIONS 11 Theorem 3.1. Fourier Theorem If f is function such tht fx n f x tht re piecewise continuous on the intervl ë;ë R then the series 1X nx 1X nx 2 + n cos + b n sin where the coecients n b n re etermine by converges pointwise to fx for ll x 2 ë;ë. R nx cos n = 2 b n = 2 R sin 2x fx fx In view of the formuls 2 mx sin 2 mx sin 2 mx cos sin cos cos nx nx nx x = m;n x = x = m;n the Fourier expnsion of f mounts to expnsion of f with respect to the bsis sin mx which is orthonorml with respect to the inner prouct f g = 2 fxgx : ; cos nx Sturm-iouville theory is generliztion of Fourier theory. It provies mens of constructing other sets of orthonoml bses for spces of functions. Theorem 3.2. Sturm-iouville Theorem Consier bounry vlue problem of the form 3.5 px y + qx + rx y = c1y +c2y = 1yb +2y b = where px n rx re smooth positive functions on the intervl ; b. Then i For ll but iscrete set S of choices of there re no solutions to 3.5. There exists miniml 2 S. n one cn rrnge the set S of missible so tht with One then hs S = f;1;2;::: ;g é1 é2 é lim n!1 n =+1 : ii To ech eigenvlue n 2 S there correspons exctly one solution n x of 3.5. iii If n ; m 2 S n n x m x re the corresponing solutions then n 6= m n x m xrx = :

2. BASES OF ORTHOGONA FUNCTIONS 12 iv The set of functions 8 é é: m = m x hr b m x 2 rxs i 1=2 m =; 1; 2; 9 é= é; form complete orthonorml bsis for spce of piecewise continuous functions on the intervl ë; bë. Tht is to sy every piecewise continous function f :ë; bë! R cn be expne s 3.6 1X fx n n t where the coecients n re etermine by 3.7 3.8 The series 3.6 converges in the sense tht lim N!1 n = fx n xrx : fx NX 2 n n x! = : Remrks: i Existence of Eigenvlues n Eigenfunctions. Sttement i is me plusible by consiering the simple exmple y + 2 y = y = y = 9 = ; = n nx yx = A sin n inee in pplictions before one cn mke use of the expnsion 3.6 one hs to rst n the eigenvlues i n so prt i my be regre s prove constructively. However n bstrct proof lso exists. ii Uniqueness of Eigenfunctions. Sttement ii follows from the existence n uniqueness theorem for secon orer liner ODE's. iii Orthogonlity Property of Eigenfunctions. et 1 n 2 be solutions of 3.5 for = 1 n = 2 respectively. Assume 1 6= 2. Then px 1 +qx +1 rx 1 = px 2 +qx +2 rx 2 = Multiplying the rst eqution by 2 n the secon by 1 n subtrcting the two equtions yiels 1 2 rx 1 2 = 2 px 1 1 px 2

3. EXAMPES 13 Integrting both sies of the expression bove between x = n x = b yiels 1 2 1 x 2 xrx = 2 x px 1 x 1 x px 2 x = 1 x 2 xpx b 1 x 2 xpx b + 1 x 2 xpx 1 x 2 xpx = rb 1 b 2 b 1 b 2 b r 1 2 1 2 Now in orer for the bounry conitions t x = c1 1 +c2 1 = c1 2 +c2 2 = to hve solutions with c1;c2 not both zero we must hve 1 2 1 2 = : This cn be seen s follows: Recll from liner lgebr tht if M is n n mtrix n v is n-imensionl vector then Mv = hs non-trivil solutions if n only if et M =. The sttement bove then follows by consiering 1 M = 2 1 2 c1 ; v = c2 Similrly in orer for the bounry conitions t x = b to be stise for c3;c4 not both zero we must hve 1 b 2 b 1 b 2 b = : : Thus we hve So if 1 6= 2. 1 2 1 x 2 xrx = 1 x 2 xrx = iv Completeness of Eigenfunctions The hr thing to unerstn is the remrkble completeness property expresse in sttement iv. The proof of this sttement is not terribly icult - however it oes require moerte igression into the Clculus of Vritions. At the en of the course time-permitting we will evelop the Clculus of Vritions n then prove iv s smple ppliction. 3. Exmples Exmple 3.3. Fourier Sine Series: y + y = y = y =

3. EXAMPES 14 This is Sturm-iouville type problem with px = rx = 1qx =. The generl solution of the ODE is given by p p yx =Asin x + B cos x however such solutions will stisfy the bounry conitions y = y = if n only if B = n p = n ; n =1; 2; 3;::: : The Sturm-iouville inner prouct is n the functions f; g = y n = fxgx r 2 nx sin constitute complete orthonorml bsis for the set of piecewise continuous functions on the intervl ;. Exmple 3.4. Fourier Cosine Series: y + y = y = y = This is Sturm-iouville type problem with px = rx = 1qx =. The generl solution of the ODE is given by p p yx =Asin x + B cos x however such solutions will stisfy the bounry conitions y = y = if n only if A = n p = 2n ; n =1; 2; 3;::: : The Sturm-iouville inner prouct is n the functions f; g = y n = fxgx r 2 2nx cos constitute complete orthonorml bsis for the set of piecewise continuous functions on the intervl ;. Exmple 3.5. Bessel Functions Exmple 3.6. egenre Functions 3.9 Exmple 3.7. Hermite Functions 3.1 x J n; n2 2 x J n; + x + J n; = 1 x 2 l e x2 H +1+ll + 1 l = +2e x2 H = In ech of the lst three exmples just s in the cse of the rst two exmples there exist ierent sets of orthogonl functions epening on the bounry contions impose.