Chapter Five - Eigenvalues, Eigenfunctions, and All That

Similar documents
Sturm-Liouville Theory

Notes on the Eigenfunction Method for solving differential equations

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

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

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

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

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

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.

Green function and Eigenfunctions

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

Differential Equations 2 Homework 5 Solutions to the Assigned Exercises

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

Homework Problem Set 1 Solutions

1 2-D Second Order Equations: Separation of Variables

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

Theoretical foundations of Gaussian quadrature

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

Best Approximation in the 2-norm

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

Homework Assignment 5 Solution Set

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

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

Basic Derivative Properties

Abstract inner product spaces

Matrix & Vector Basic Linear Algebra & Calculus

Wave Equation on a Two Dimensional Rectangle

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

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

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?

Problem Set 2 Solutions

Review of Gaussian Quadrature method

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

Overview of Calculus

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

The Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5

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) Ω.

B Veitch. Calculus I Study Guide

Conservation Law. Chapter Goal. 6.2 Theory

Section 6.1 Definite Integral

Antiderivatives Introduction

(PDE) u t k(u xx + u yy ) = 0 (x, y) in Ω, t > 0, (BC) u(x, y, t) = 0 (x, y) on Γ, t > 0, (IC) u(x, y, 0) = f(x, y) (x, y) in Ω.

Bases for Vector Spaces

Lecture 2e Orthogonal Complement (pages )

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

The Regulated and Riemann Integrals

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

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

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

21.6 Green Functions for First Order Equations

Orthogonal Polynomials

Suppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2.

u(x, y, t) = T(t)Φ(x, y) 0. (THE EQUATIONS FOR PRODUCT SOLUTIONS) Plugging u = T(t)Φ(x, y) in (PDE)-(BC) we see: There is a constant λ such that

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

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

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

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

Generalized Cobb-Douglas function for three inputs and linear elasticity

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

2.4 Linear Inequalities and Interval Notation

= f (c) f (c) the height of the rectangle guaranteed by the MVT for integrals.

Section 7.1 Integration by Substitution

Chapter 0. What is the Lebesgue integral about?

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

Math 142: Final Exam Formulas to Know

Best Approximation. Chapter The General Case

Example Sheet 2 Solutions

The Fundamental Theorem of Calculus Part 2, The Evaluation Part

ODE: Existence and Uniqueness of a Solution

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

1.1 Functions. 0.1 Lines. 1.2 Linear Functions. 1.3 Rates of change. 0.2 Fractions. 0.3 Rules of exponents. 1.4 Applications of Functions to Economics

Course 2BA1 Supplement concerning Integration by Parts

The Wave Equation I. MA 436 Kurt Bryan

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?

Variational Techniques for Sturm-Liouville Eigenvalue Problems

Section 6: Area, Volume, and Average Value

Final Exam Review. Exam 1 Material

Lecture Solution of a System of Linear Equation

CHAPTER 9 BASIC CONCEPTS OF DIFFERENTIAL AND INTEGRAL CALCULUS

Chapter 6 Techniques of Integration

20 MATHEMATICS POLYNOMIALS

Summary: Method of Separation of Variables

Lecture 3. Limits of Functions and Continuity

Introduction to Complex Variables Class Notes Instructor: Louis Block

M344 - ADVANCED ENGINEERING MATHEMATICS

Taylor Polynomial Inequalities

Main topics for the First Midterm

New Expansion and Infinite Series

Year 11 Matrices. A row of seats goes across an auditorium So Rows are horizontal. The columns of the Parthenon stand upright and Columns are vertical

REGULAR TWO-POINT BOUNDARY VALUE PROBLEMS

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

SUPPLEMENTARY NOTES ON THE CONNECTION FORMULAE FOR THE SEMICLASSICAL APPROXIMATION

AP Calculus AB First Semester Final Review

September 13 Homework Solutions

7.2 The Definite Integral

Riemann Sums and Riemann Integrals

10 Vector Integral Calculus

TIME VARYING MAGNETIC FIELDS AND MAXWELL S EQUATIONS

Lecture 6: Coding theory

Quadratic Forms. Quadratic Forms

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

Transcription:

Chpter Five - Eigenvlues, Eigenfunctions, n All Tht The prtil ifferentil eqution methos escrie in the previous chpter is specil cse of more generl setting in which we hve n eqution of the form L 1 xux,tl tux,t Fx,t for x D n t, in which u is specifie on the ounry of D s re initil conitions t t. We ssume tht L 1 x is liner ifferentil opertor in x n L t is liner ifferentil opertor in t. In the previous chpter, L 1 xx x, n we were fortunte tht this opertor h n infinite orthogonl collection of eigenfunctions which pproximte functions in nice wy. We shll see wht hppens in more generl setting. First it of liner lger review. Definition. Suppose V is liner spce together with n inner prouct. A liner opertor L : V V such thtlf,g f,lg for ll f,g V is si to e self-joint. Let V e liner spce of nice functions efine on n intervl x with the inner prouct f,g fxgx. Define the liner opertor L y L px qxx. We ssume tht p n q re rel n continuous n tht p is continuously ifferentile n positive. Such n opertor is clle Sturm-Liouville opertor. Definition. Suppose r is rel continuous n positive function on x.a sclr such tht L r for some nonzero V is clle n eigenvlue of L, n the function is n eigenfunction. For resons tht will soon e cler, we woul very much like to hve our liner opertor L e self-joint. Thus, we wntlf,g f,lg for ll f n g in the spce V. Let s compute: Lf,gf,Lg LfxgxfxLgx px f px g gxqxfxgx fxqxgxfx 1

px f px gx f gx fx g px g fx pgf fg p gf fg. In other wors, L will e self-joint if we choose V to e liner spce of functions such tht p gf fg p gf fg. Exmples. The opertor L is Sturm-Liouville opertor on the intervl x, with px 1 n qx.the self-joint ounry conitions then ecome gf fg gf fg Thus, if we choose V to e the spce of functions such tht, the opertor L is self-joint. So why shoul we cre out self-joint opertors? Let s see. Proposition 1. Every eigenvlue of self-joint opertor is rel. Proof: Suppose L r for some nonzero. Then L,,L ecuse L is self-joint, n this ecomesr,,r, orr,,r. In other wors, rx x rx x. The integrl is not zero, n so we hve. Proposition. Eigenfunctions corresponing to ifferent eigenvlues of self-joint opertor L re orthogonl with respect to the inner prouct f,g r rf,g f,rg rxfxgx. Proof: Suppose L r n L r for. Agin, from the fct tht L is self-joint, we know tht L,,L. Thus, r,,r.from the previous proposition we know

tht n re rel. Thus, r,,r, or, r. But n so it must e true tht, r. Nothing we hve one gurntees the existence of eigenvlues of L. There re, in fct, n infinite numer of them, giving n infinite orthogonl set of eigenfunctions. A proof is too much for these notes. The celerte Sturm-Liouville Theorem sys even more. It sys tht every nice function f cn e expne in series f n n n1 of the eigenfunctions n, n this series converges in the men to f. Exmple. Consier the prolem u xx u t, for x n t u,t n u t,t ux, fx In solving this prolem, I hope it is cler why we egin with the eigenvlue prolem,. Note this is n exmple of the Sturm-Liouville prolem we hve just iscusse, so we know wht to expect. First, we know tht we must in orer to hve nonzero solutions to the prolem. To remin us of tht, let.then x Acosx Bsinx, n the ounry conitions ecome n A 3

Bcos. If or B re zero, then, so it must e true tht cos. Thus must e n o multiple of / : n n1 1 n 1, for n 1,,3, Our eigenvlue re thus n n n 1,n he corresponing eigenfunctions re n x sin n 1 x. Continuing, let ux,t n tsin n x.from here on, things look fmilir. n1 u xx u t n n t n t sin n x. n1 As usul, n n t n t tells us tht n t n e nt n e n t, n so ux,t n e n t sin n x. n1 The constnts n re etermine from the initil conition ux, fx : Thus, ux, fx n sin n x. n1 f,sin n n x sin n x,sin n x fxsin n 4

Exercises 1. Fin ll eigenvlues n eigenfunctions for,.. Fin ll eigenvlues n eigenfunctions:, Let n e non-negtive integer, n consier the Sturm-Liouville prolem x n x x xx, x c. Note tht this opertor oes not quite fit our efinitions. The coefficient qx n x of is not nice t the left en point of the intervl uner consiertion n px x is zero t the left en point of the intervl. This is wht is known s singulr Sturm-Liouville prolem. If, however, we consier functions tht, together with their erivtives, re resonly well-ehve t the left en point, then our self-joint ounry conitions ecome p gf fg p gf fg c gf fg, n re stisfie y functions in the spce of those nice functions tht vnish t x. The eqution x n x x xx, or x x n x is the worl-fmous Bessel s Eqution. We shll seek solutions y mens of the Froenius metho. Tht is, ssume solutionof the form 5

Then x c k x k, c. k k k k c k x k, kc k x k1, kk1c k x k, n we hve x x n x k k Writing c k x k c k x k gives us k k kk1kn c k x k c k x k x x n x k n c k x k k k c k x k n c x 1 n c 1 x 1 k n c k c kx k. k In orer for the ifferentil eqution to e stisfie, ech of the coefficients of the powers of x must e zero: n c 1 n c 1 k n c k c k, for k,3,4,. In the first eqution, c n so it must e true tht n, giving us n, or n. Let s strt with n. The reminer of the equtions ecome 6

n1c 1 knkc k c k, for k,3,4,. Thus, c 1 c 3 c 5, n c k knk c k, for k,3,4,. Letting k m results in Then c m mnm c m1 mnm c m1. c m mnm c m1 1 m m n! m m!nm! c where c is ny constnt. For the ske of netness, choose c 1 n! n. Then c m 1 m m n! m m!nm! c n 1 m m!nm! m, n our solutions, t lst, x n m x x n m 1 m m!nm! 1 m m!nm! m x mn,or x m. The function J n x x n m 1 m m!nm! x m is the Bessel function of the first kin of orer n. Our solution is thus x J n x. 7

A secon linerly inepenent solution is foun reucing the orer of the originl eqution. It turns out not to e nice t x, n so we hve solutionsx AJ n x.the ounry conition is thus J n. It is known tht J n hs n infinite numer of non-negtive zeros. These n ll sorts of other informtion out Bessel functions cn e foun in most hnooks or y mens of most computer lger systems; e.g., Mple or Mthemtic. Here re some pictures. First, J : 1.8.6.4. -. 4 6 8 1 1 14 16 18 x -.4 Now here re J 1,J,n J 3 :.6.4. -. 4 6 8 1 1 14 16 18 x -.4 -.6 It is cler then tht our eigenvlues re nm, with nm z nm /, where z nm is the m th zero of J n. Here is short tle of some of these vlues: m n n 1 n n 3 1. 45 3. 83 5. 135 6. 379 5. 5 7. 16 8. 417 9. 76 3 8. 654 1. 173 11. 6 13. 17 4 11. 79 13. 33 14. 796 16. 4 The corresponing eigenfunctions re, of course, J n nm x. Note there is non-constnt weight function here, so our orthogonlity ecomes 8

xj n nk xj n nl x, for k l. Exmple. Consier the prolem x xu x u tt, x L, t ul,t ux, fx, u tt x, gx. The solution u gives the isplcement of hnging chin of length L : It shoul e cler y now why we re intereste in the eigenvlue prolem x L, n nice. This is tntlizingly close to Bessel s eqution of orer zero. We coul use the metho of Froenius to fin solutions, etc. Inste, let s mke chnge of vrile. Let z x, or x z /4.Define z z /4. Then, Our ifferentil eqution thus ecomes z z, or z z. z x z z z z. z z 4 z z, or 9

Now we hve Bessel s eqution with n. The solution is J z J x.the ounry conitionl is J L, n so we hve the eigenvlues m z m L, where z m is the m th zero of J. We set ux,t m tj m m1 x. Then x xu x u tt m m t m tj m x. m1 This gives us m t m cos m t m sin m t. Thus, ux,t m cos m t m sin m tj m m1 x Use the initil conitions to fin the constnts m, m : ux, m J m m1 x fx. Hence, m L fxj m L J m x. x Also, n so, u t x, m m J m m1 x, 1

m L gxj m L m J m x x Exercise 3. Fin the first four terms of the series for u in cse L 1, fx x1x, n gx. Drw some pictures of your pproximtion for ifferent vlues of t. 11