13.07.2015 Views

Mathematical Methods for Physics and Engineering - Matematica.NET

Mathematical Methods for Physics and Engineering - Matematica.NET

Mathematical Methods for Physics and Engineering - Matematica.NET

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

CALCULUS OF VARIATIONSUsing (22.13) <strong>and</strong> the fact that y does not appear explicitly, we obtain(∂ρ ∂y )− ∂ (τ ∂y )=0.∂t ∂t ∂x ∂xIf, in addition, ρ <strong>and</strong> τ do not depend on x or t then∂ 2 y∂x = 1 ∂ 2 y2 c 2 ∂t , 2where c 2 = τ/ρ. This is the wave equation <strong>for</strong> small transverse oscillations of a tautuni<strong>for</strong>m string. ◭22.6 General eigenvalue problemsWe have seen in this chapter that the problem of finding a curve that makes thevalue of a given integral stationary when the integral is taken along the curveresults, in each case, in a differential equation <strong>for</strong> the curve. It is not a greatextension to ask whether this may be used to solve differential equations, bysetting up a suitable variational problem <strong>and</strong> then seeking ways other than theEuler equation of finding or estimating stationary solutions.We shall be concerned with differential equations of the <strong>for</strong>m Ly = λρ(x)y,where the differential operator L is self-adjoint, so that L = L † (with appropriateboundary conditions on the solution y) <strong>and</strong>ρ(x) is some weight function, asdiscussed in chapter 17. In particular, we will concentrate on the Sturm–Liouvilleequation as an explicit example, but much of what follows can be applied toother equations of this type.We have already discussed the solution of equations of the Sturm–Liouvilletype in chapter 17 <strong>and</strong> the same notation will be used here. In this section,however, we will adopt a variational approach to estimating the eigenvalues ofsuch equations.Suppose we search <strong>for</strong> stationary values of the integral∫ b []I = p(x)y ′2 (x) − q(x)y 2 (x) dx, (22.22)awith y(a) =y(b) =0<strong>and</strong>p <strong>and</strong> q any sufficiently smooth <strong>and</strong> differentiablefunctions of x. However, in addition we impose a normalisation conditionJ =∫ baρ(x)y 2 (x) dx = constant. (22.23)Here ρ(x) is a positive weight function defined in the interval a ≤ x ≤ b, butwhich may in particular cases be a constant.Then, as in section 22.4, we use undetermined Lagrange multipliers, § <strong>and</strong>§ We use −λ, rather than λ, so that the final equation (22.24) appears in the conventional Sturm–Liouville <strong>for</strong>m.790

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!