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

Similar documents
AMS 212A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang, UCSC. ( ) cos2

Green function and Eigenfunctions

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

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

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

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

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH

Sturm-Liouville Theory

Chapter Five - Eigenvalues, Eigenfunctions, and All That

The Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5

u(t)dt + i a f(t)dt f(t) dt b f(t) dt (2) With this preliminary step in place, we are ready to define integration on a general curve in C.

Variational Techniques for Sturm-Liouville Eigenvalue Problems

REGULAR TWO-POINT BOUNDARY VALUE PROBLEMS

Differential Equations 2 Homework 5 Solutions to the Assigned Exercises

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

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

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

Wave Equation on a Two Dimensional Rectangle

1 E3102: a study guide and review, Version 1.0

Chapter 28. Fourier Series An Eigenvalue Problem.

Section 6.1 Definite Integral

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

Improper Integrals. Introduction. Type 1: Improper Integrals on Infinite Intervals. When we defined the definite integral.

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

21.6 Green Functions for First Order Equations

10 Elliptic equations

Module 9: The Method of Green s Functions

ODE: Existence and Uniqueness of a Solution

p-adic Egyptian Fractions

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.

Topics Covered AP Calculus AB

1 2-D Second Order Equations: Separation of Variables

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

Review of Gaussian Quadrature method

Section 3.3: Fredholm Integral Equations

Bases for Vector Spaces

M344 - ADVANCED ENGINEERING MATHEMATICS

Lecture 24: Laplace s Equation

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

Improper Integrals, and Differential Equations

Polynomials and Division Theory

Linearity, linear operators, and self adjoint eigenvalue problems

Thomas Whitham Sixth Form

Math Solutions to homework 1

4 Sturm-Liouville Boundary Value Problems

Definite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +

Linear Systems with Constant Coefficients

Lecture 2e Orthogonal Complement (pages )

STURM-LIOUVILLE PROBLEMS: GENERALIZED FOURIER SERIES

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

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

Name Solutions to Test 3 November 8, 2017

September 13 Homework Solutions

Best Approximation in the 2-norm

3 Mathematics of the Poisson Equation

Math Fall 2006 Sample problems for the final exam: Solutions

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

Simplifying Algebra. Simplifying Algebra. Curriculum Ready.

1B40 Practical Skills

MA123, Chapter 10: Formulas for integrals: integrals, antiderivatives, and the Fundamental Theorem of Calculus (pp.

Summary: Method of Separation of Variables

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

APM346H1 Differential Equations. = u x, u = u. y, and u x, y =?. = 2 u t and u xx= 2 u. x,t, where u t. x, y, z,t u zz. x, y, z,t u yy.

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

Section 6.1 INTRO to LAPLACE TRANSFORMS

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

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

MTH 505: Number Theory Spring 2017

Lecture Solution of a System of Linear Equation

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Lecture 4 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell

1 E3102: A study guide and review, Version 1.2

PDE Notes. Paul Carnig. January ODE s vs PDE s 1

Numerical Methods for Partial Differential Equations

4.1. Probability Density Functions

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

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

Math 5440 Problem Set 3 Solutions

Example Sheet 2 Solutions

Lecture 1. Functional series. Pointwise and uniform convergence.

Section 4: Integration ECO4112F 2011

Elliptic Equations. Laplace equation on bounded domains Circular Domains

Math 61CM - Solutions to homework 9

10 Vector Integral Calculus

Partial Differential Equations

Math 360: A primitive integral and elementary functions

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?

2. VECTORS AND MATRICES IN 3 DIMENSIONS

AP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

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

MATH1050 Cauchy-Schwarz Inequality and Triangle Inequality

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

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

1 1D heat and wave equations on a finite interval

Determinants Chapter 3

DYNAMICAL SYSTEMS SUPPLEMENT 2007 pp Natalija Sergejeva. Department of Mathematics and Natural Sciences Parades 1 LV-5400 Daugavpils, Latvia

Lecture 2: January 27

Math 5440 Problem Set 3 Solutions

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

Transcription:

AMS A Applied Mthemticl Methods I Lecture 6 Copyright y Hongyun Wng, UCSC Recp of Lecture 5 Clssifiction of oundry conditions Dirichlet eumnn Mixed Adjoint opertor, self-djoint opertor Sturm-Liouville prolem in self-djoint form: Lu ( ) r( x)u, x (, ) u( )+ u( ) u ( )+ u( ) with respect to the inner product u, v u( x)v( x)r( x)dx where the liner differentil opertor is Lu ( ) p( x)u x ( )u ( ) x + q x Results of Sturm-Liouville theory: u, r( x) L( v) r( x) Lu All eigenvlues re rel. ( ), v (tht is, r x ( ) L i For ech eigenvlue, we cn find rel eigenfunction. () is self-djoint). Eigenfunctions for different eigenvlues re orthogonl to ech other. Results of Sturm-Liouville theory (continued) Since eigenvlues nd eigenfunctions re oth rel, we only need to work with rel functions. All eigenvlues re simple. (i.e., one corresponding eigenfunction for ech eigenvlue.) - -

AMS A Applied Mthemticl Methods I Suppose Clim: Proof: u(x) nd v(x) re eigenfunctions corresponding to eigenvlue. u(x) nd v(x) re proportionl to ech other. Lu ( ) r( x)u L( v) r( x)v > ul( v) vl( u) u( r( x)v) v( r( x)u) On the other hnd, using the differentil form of L( ), we write ul( v) vl( u) u( ( pv x ) x + qv) v( ( pu x ) x + qu) u( pv x ) x vpu ( x ) x (Using Lemm of Lecture 5) ( ) x puv ( x vu x ) Comining these two results, we hve ( puv ( x vu x )) x > puv ( x vu x ) const Since oth u(x) nd v(x) stisfy the oundry conditions, we hve ( uv vu) x (Lemm of Lecture 5) [ ] > pu ( v vu), x, > uv vu, > x [, ] v u, x, v > u const End of proof [ ] Eigenvlue sequence Eigenvlues form n unounded strictly incresing infinite sequence < < < < n < - -

AMS A Applied Mthemticl Methods I lim n + n+ We will prove this fter we descrie the completeness of eigenfunctions. Completeness of eigenfunctions The set of eigenfunctions v ( x), v ( x), v ( x),, ( x), piecewise continuous functions on [, ]. Tht is, ny piecewise continuous function ƒ(x) cn e expnded s f( x) x ( ) where coefficient is given y f,, { } is complete sis for ll f( x) ( x)r( x)dx ( ( x) ) r x ( )dx Below we do 3 things: prove the eigenvlue sequence, introduce the Ryleigh quotient, nd then prove the completeness Proof of eigenvlue sequence Eigenvlue nd eigenfunction u(x) stisfy the Sturm-Liouville eqution ( p( x)u x ) x + q( x)u r( x)u, x (, ) nd oundry conditions u( )+ u( ) u ( )+ u( ) The key step of the proof is Prufer sustitution: ( ) p( x)u x ( x) ( x)cos ( x) u( x) ( x)sin( ( x) ) Before the Prufer sustitution, we hve second order liner ODE for u(x). After the Prufer sustitution, (s we will see) we hve - 3 -

AMS A Applied Mthemticl Methods I first order non-liner ODE system for (x) nd (x). Let us look t the dvntges of Prufer sustitution. Advntge #: [ ] ( x), x, Suppose ( x ) t some point x in [, ]. > p(x ) p( x )u x ( x ) u( x ) > u x ( x ) u( x ) Tht is, u(x) stisfies the zero initil conditions t x. The Sturm-Liouville eqution is second order ODE of u(x) > u( x) (due to the zero initil conditions t x ) In the Sturm-Liouville prolem, we re only interested in non-trivil solutions. Advntge #: Boundry conditions ffect only Boundry condition t x : u( )+ u( ) ( ) > ( )sin( ( ) )+ p ( ) > p ( )sin( ( ) )+ cos( ( ) ) > sin( ( ) ) cos ( ( )) where stisfies cos( ) ( ) ( ) p, sin( ) ( p) + The oundry condition ecomes ( ) n ( p ( )) + > ( ) + n We cn shift the whole function (x) to write the oundry condition t x s - 4 -

AMS A Applied Mthemticl Methods I BC t x : ( ), < After the shifting, we write the oundry condition t x s BC t x : ( ) μ + n, < μ ote: we use shifting nd pick suitle vlue of n to mke < nd < μ, which is importnt in the proof of Lemm nd proof of eigenvlue sequence. Advntge #3: The evolution eqution of (x) is independent of (x). (x) stisfies the first order non-liner ODE x ( r( x)+ q( x) )sin + p x ( ) cos Key steps in derivtion: equtions for (x) nd (x): Eqution : ( p( x)u x ) x + ( r( x)+ q( x) )u (the Sturm-Liouville eqution) > ( cos) x + ( r( x)+ q( x) )sin Eqution : cos p( x)u x p( x) ( sin ) x (See Exercise #6 for derivtion). Lemm : Let (x, ) e the solution of x ( r( x)+ q( x) )sin + p x ( ), < ( ) cos where p(x) > nd r(x) > for x in [, ] (from the Sturm-Liouville prolem). Then, (x, ) hs the 3 properties elow. For x >, ) (x, ) is strictly incresing function of. ) lim x, + c) lim x, (For proof of Lemm, see Appendices of Lecture 6 in seprte PDF file). - 5 -

AMS A Applied Mthemticl Methods I ow we use the results of Lemm to prove the eigenvlue sequence. (, ) < μ (since < μ ) lim > μ < (, )< μ for negtive nd lrge. As increses to infinity, (, ) lso increses to infinity. At some vlue,, we will hve (, ) μ. > is the smllest vlue of for which (, ) μ + n is stisfies. > is the smllest eigenvlue. As increses further, t some vlue,, we will hve (, ) μ + ; > is the smllest eigenvlue fter. At some vlue,, we will hve (, ) μ + ; > is the smllest eigenvlue fter. At some vlue, 3, we will hve (, 3 ) μ + 3; > 3 is the smllest eigenvlue fter. Thus, we hve strictly incresing sequence of eigenvlues < < < < n < ext we prove lim n +. n+ Since { n } increses monotoniclly, we hve either lim n + or lim n finite. n+ n+ Suppose lim n ( c) finite. We hve n+ ( c) > n for ny n > (, ( c) )> (, n ) μ + n for ny n > (, ( c) )+, which contrdicts with ( c) finite. End of proof of eigenvlue sequence Ryleigh quotient Definition: - 6 -

AMS A Applied Mthemticl Methods I Ryleigh quotient is defined s Ru ( ) u, r( x) Lu ( ) u, u ul( u)dx u r( x)dx Ryleigh quotient on eigenfunctions: Suppose < < < < n < is the eigenvlue sequence nd { v ( x), v ( x), v ( x),, ( x), } is the sequence of corresponding eigenfunction. It is strightforwrd to verify tht Ryleigh quotient stisfies Rv ( n ) n, n,,, R min n { } R,, (See Exercise #6 for derivtion). This lst result motivtes the Ryleigh principle elow. Ryleigh principle: Let C BC u( x)u C nd stisfies u ( )+ u ( ) u ( )+ u( ) Consider the minimiztion of R(u) over C BC. Let u rg min Ru ( ) uc BC, u Clim : u is n eigenfunction corresponding to eigenvlue nd min Ru ( ) uc BC, u Proof of Clim : Consider function of sclr vrile s - 7 -

AMS A Applied Mthemticl Methods I gs ( ) Ru ( + su) where u C BC By definition of u, we hve g( ) for ny u C BC The derivtives of numertor nd denomintor of R(u + s u) re d ds u + su, u + su s u, u d ds u + su, r x ( ) Lu ( + su) s (See Exercise #6 for derivtion). u, r( x) Lu ( ) The derivtive of R(u + s u) hs the expression d ds Ru ( + su) s u, u, u r( x) Lu ( ) + Ru ( ) u, u u, u u, (See Exercise #6 for derivtion). r( x) Lu ( )+ Ru ( )u The derivtive of R(u + s u) is zero for ny u C BC. d ds Ru ( + su) for ny u C BC s > u, r( x) Lu ( )+ Ru ( )u for ny u C BC > r( x) ( )+ Ru ( )u > Lu ( ) R( u )r( x)u > R(u ) is n eigenvlue nd u is corresponding eigenfunction. ext we show Ru ( ). By definition of u, we hve Ru ( ) min Ru ( ) uc BC > Ru ( ) Rv ( ) (v is n eigenfunction corresponding to ). On the other hnd, R(u ) is n eigenvlue nd is the smllest eigenvlue. Thus, we hve - 8 -

AMS A Applied Mthemticl Methods I Ru ( ) Comining these two results, we conclude Ru ( ) min Ru ( ) End of proof of Clim uc BC ext we look t the minimum of Ryleigh quotient over suset of C BC. Let W spn{ v,, v } { nd u W } W uuc BC Consider the minimiztion of R(u) over W. Let u + rg min Ru ( ) uw, u Clim : u + is n eigenfunction corresponding to eigenvlue + nd + min Ru ( ) uw, u (For proof of Clim, see Appendices of Lecture 6 in seprte PDF file) Ryleigh principle provides wy of pproximting the lowest few eigenvlues, especilly the smllest eigenvlue. Exmple: u u u( ), u() This is Sturm Liouville prolem with p( x), r( x), q( x), Lu ( ) u We solved it previously. The eigenvlues re n ( n + ), n,,, ow we use Ryleigh principle to pproximte the smllest eigenvlue. (Skip the detils in Lecture) We try w( x) x( x) (it hs to stisfy the oundry conditions) - 9 -

AMS A Applied Mthemticl Methods I The numertor nd denomintor of Ryleigh quotient re w, r( x) L( w) w( x) w ( x)dx x( x)dx 3 ( ( )) dx w, w w x x x 3 ( ) dx > Rw ( ) w, r x w, w ( ) L( w) > min Ru ( ) Rw ( ) uc BC The true vlue is 9.87. 3 3 Proof of the completeness of eigenfunctions For piecewise continuous function ƒ(x), we show tht where lim f n f,, u, v u( x)v( x)r( x)dx u u, u The proof consists of results -4 elow. Result : f v j for j This follows directly from the definition of. Result : f min { c n } f c n - -

AMS A Applied Mthemticl Methods I This follows directly from Result ove f + v n n, f + n f, f v n n + n, n > f c n f Result 3: We only need to show tht for piecewise continuous function ƒ(x), we hve Result 4: for ny > there exists nd {c n } such tht f c n < (***) Since C BC is dense in the set of piecewise continuous functions, we only need to show tht for f( x)c BC, we hve for ny > there exists nd {c n } such tht f c n < (***) For proof of (***), see Appendices of Lecture 6 in seprte PDF file. - -