Introduction to Sturm-Liouville Theory

Transcription

1 Introduction to R. C. Trinity University Partial Differential Equations April 10, 2014

2 Sturm-Liouville problems Definition: A (second order) Sturm-Liouville (S-L) problem consists of A Sturm-Liouville equation on an interval: together with (p(x)y ) +(q(x)+λr(x))y = 0, a < x < b, Boundary conditions, i.e. specified behavior of y at x = a and x = b. Definition: A function y 0 that solves an S-L problem is called an eigenfunction, and the corresponding value of λ is called its eigenvalue.

3 Examples The boundary value problem y +λy = 0, y( p) = y(p), y ( p) = y (p), is an S-L problem on the interval [ p,p]. We have seen that ( ) nπ 2 Eigenvalues: λ = λ n = (n N 0 ) p ( ) ( ) nπx nπx Eigenfunctions: y = y n = a n cos +b n sin p p The boundary value problem y +λy = 0, y(0) = y(l) = 0, is an S-L problem on the interval [0,L]. We have seen that ( nπ ) 2 Eigenvalues: λ = λ n = (n N) L ( nπx ) Eigenfunctions: y = y n = c n sin L

4 The boundary value problem y +λy = 0, y(0) = 0, y (L) = κy(l), is an S-L problem on the interval [0,L]. We have seen that ( Eigenvalues: λ = λ n = µ 2 n tanµ n = µ ) n κ Eigenfunctions: y = y n = c n sin(µ n x) The boundary value problem x 2 y +xy +(λ 2 x 2 p 2 )y = 0, y(0) finite, y(a) = 0 is an S-L problem on the interval [0,a]. We have seen that ( Eigenvalues: λ 2 = λ 2 n = αpn ) 2 (n N) a ( αpn x ) Eigenfunctions: y = y n = c n J p a where α pn is the nth positive zero of J p.

5 Inner products of eigenfunctions Suppose (y j,λ j ), (y k,λ k ) are eigenfunction/eigenvalue pairs of an S-L problem: (py j ) +(q +λ j r)y j = 0, (py k ) +(q +λ k r)y k = 0. Multiply the first by y k and the second by y j, then subtract to get (py j) y k (py k ) y j +(λ j λ k )y j y k r = 0. Moving the λ-terms to one side and adding zero, we get (λ j λ k )y j y k r = (py k ) y j (py j) y k = (py k ) y j +py k y j py j y k (py j ) y k = (py k y j py j y k) = ( p(y k y j y jy k ) ).

6 If λ j λ k, we can divide by λ j λ k to get y j y k r = Now integrate both sides to obtain: Proposition ( ) p(y k y j y j y k). λ j λ k If (y j,λ j ), (y k,λ k ) are eigenfunction/eigenvalue pairs for an S-L problem on the interval [a,b] and λ j λ k, then their inner product with respect to the weight function r(x) is b y j,y k = y j (x)y k (x)r(x)dx a ( ) p(x) y k (x)y j(x) y j (x)y k(x) b =. λ j λ k a

7 Remarks We will be interested in situations when ( ) p(x) y k (x)y j(x) y j (x)y k(x) y j,y k = λ j λ k b a = 0, i.e. when eigenfunctions with distinct eigenvalues are orthogonal. Examples include: Periodic S-L problems, i.e. p(a) = p(b) and y(a) = y(b), y (a) = y (b). S-L problems satisfying p(a) = 0 and y(a) is finite, y(b) = 0.

8 Example Use the preceding results to establish the orthogonality of the trigonometric system { ( ) ( ) ( ) ( πx 2πx πx 1,cos,cos,...,sin,sin 2 πx ) },... p p p p on the interval [ p,p] with respect to the weight function w(x) = 1. ( ) ( ) nπx nπx The functions y = cos and y = sin are p p eigenfunctions of the periodic S-L problem y +λy = 0, y( p) = y(p), y (p) = y ( p) with eigenvalue λ = ( ) nπ 2. p

9 Since r(x) 1 in this case, we automatically find that p ( ) ( ) mπx nπx p ( ) ( ) mπx nπx cos cos dx = cos sin dx = p p p p p p p ( ) ( ) mπx nπx sin sin dx = 0 for m n. p p p ( ) ( ) nπx nπx We must verify orthogonality of cos and sin p p manually, since they have the same eigenvalue: p ( ) ( ) nπx nπx cos sin dx = p ( ) nπx p p p 2nπ sin2 = 0. p p That covers every case, so we re done. p

10 Example Use the preceding results to establish the orthogonality of the Bessel function system { ( αp1 x ) ( αp2 x ) ( αp3 x ) } J p,j p,j p,..., a a a on the interval [0,a] with respect to the weight function w(x) = x. ( αpnx) The function y = J p a is an eigenfunction of the S-L problem x 2 y +xy +(λ 2 x 2 p 2 )y = 0, y(0) finite, y(a) = 0 ) 2. The S-L form of this equation is with eigenvalue λ 2 = ( α pn a (xy ) + ) ( p2 x +λ2 x y = 0, which shows that this problem is of the second type mentioned above. Since r(x) = x, we re done.

11 Regular Sturm-Liouville problems A regular Sturm-Liouville problem on [a, b] has the form where: (p(x)y ) +(q(x)+λr(x))y = 0, a < x < b, c 1 y(a)+c 2 y (a) = 0, (1) d 1 y(b)+d 2 y (b) = 0, (2) (c 1,c 2 ) (0,0) and (d 1,d 2 ) (0,0); p, p, q and r are continuous on [a,b]; p and r are positive on [a,b]. Remark: We will focus on the boundary conditions (1) and (2), since they yield orthogonality of eigenfunctions.

12 Examples The boundary value problem y +λy = 0, 0 < x < L, y(0) = y(l) = 0, is a regular S-L problem (p(x) = 1, r(x) = 1, q(x) = 0). The boundary value problem ((x 2 +1)y ) +(x +λ)y = 0, 1 < x < 1, y( 1) = y (1) = 0, is a regular S-L problem (p(x) = x 2 +1, q(x) = x, r(x) = 1).

13 Non-example Although important, the boundary value problem x 2 y +xy +(λ 2 x 2 p 2 )y = 0, 0 < x < a, y(0) finite, y(a) = 0, is not a regular Sturm-Liouville problem. In Sturm-Liouville form we had p(x) = r(x) = x, q(x) = p 2 /x. p and r are not positive when x = 0. q is not continuous when x = 0. The boundary condition at x = 0 is irregular. This is an example of a singular Sturm-Liouville problem.

14 Orthogonality in regular S-L problems Suppose y j and y k are eigenfunctions of a regular S-L problem with distinct eigenvalues. The boundary condition at x = a gives c 1 y j (a)+c 2 y j (a) = 0, ( yj (a) y c 1 y k (a)+c 2 y k (a) = 0, j (a) ) ( ) ( ) c1 0 y k (a) y k (a) =. c 2 0 Since (c 1,c 2 ) (0,0) the determinant must be zero: y j (a)y k (a) y k(a)y j(a) = 0. Likewise, the boundary condition at x = b gives y j (b)y k (b) y k(b)y j(b) = 0. This means that ( ) p(x) y k (x)y j(x) y j (x)y k(x) y j,y k = λ j λ k b a = 0.

15 Eigenfunctions and eigenvalues of regular S-L problems Theorem The eigenvalues of a regular S-L problem on [a,b] form an increasing sequence of real numbers with lim n λ n =. λ 1 < λ 2 < λ 3 < Eigenfunctions corresponding to distinct eigenvalues are orthogonal on [a,b] with respect to the weight r(x). Moreover, the eigenfunction y n corresponding to λ n is unique (up to a scalar multiple), and has exactly n 1 zeros in the interval a < x < b.

16 Remarks Aside from orthogonality, the proof of this result is beyond the scope of our class. Up to the moreover statement, this result holds for many irregular S-L problems (as we have seen). Orthogonality will help us extract coefficients from eigenfunction expansions of the form a n y n. n=1 The results of S-L theory unify and explain all of the ODE boundary value problems we have encountered throughout the semester!