Math 220B - Summer 2003 Homework 1 Solutions

Similar documents
4 1-D Boundary Value Problems Heat Equation

1 Heat Equation Dirichlet Boundary Conditions

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Math 124B January 31, 2012

6 Wave Equation on an Interval: Separation of Variables

Math 124B January 17, 2012

Strauss PDEs 2e: Section Exercise 1 Page 1 of 7

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form

14 Separation of Variables Method

Week 5 Lectures, Math 6451, Tanveer

Assignment 7 Due Tuessday, March 29, 2016

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations)

Wave Equation Dirichlet Boundary Conditions

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

Lecture Notes 4: Fourier Series and PDE s

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

MA 201: Partial Differential Equations Lecture - 10

4 Separation of Variables

Chapter 5. Wave equation. 5.1 Physical derivation

Lecture Notes for Math 251: ODE and PDE. Lecture 32: 10.2 Fourier Series

C. Fourier Sine Series Overview

Legendre Polynomials - Lecture 8

1D Heat Propagation Problems

Math 220a - Fall 2002 Homework 6 Solutions

David Eigen. MA112 Final Paper. May 10, 2002

MA 201: Partial Differential Equations Lecture - 11

Higher dimensional PDEs and multidimensional eigenvalue problems

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

Suggested Solution to Assignment 5

221B Lecture Notes Notes on Spherical Bessel Functions

CS229 Lecture notes. Andrew Ng

SCHOOL OF MATHEMATICS AND STATISTICS. Mathematics II (Materials) Section A. Find the general solution of the equation

Strauss PDEs 2e: Section Exercise 4 Page 1 of 6

ODE Homework 2. Since M y N x, the equation is not exact. 2. Determine whether the following equation is exact. If it is exact, M y N x 1 x.

Course 2BA1, Section 11: Periodic Functions and Fourier Series

MAS 315 Waves 1 of 8 Answers to Examples Sheet 1. To solve the three problems, we use the methods of 1.3 (with necessary changes in notation).

MATH 251 Final Examination December 16, 2015 FORM A. Name: Student Number: Section:

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons

Method of Separation of Variables

$, (2.1) n="# #. (2.2)

Week 6 Lectures, Math 6451, Tanveer

A SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

SAMPLE FINAL EXAM SOLUTIONS

Strauss PDEs 2e: Section Exercise 4 Page 1 of 5. u tt = c 2 u xx ru t for 0 < x < l u = 0 at both ends u(x, 0) = φ(x) u t (x, 0) = ψ(x),

6 Non-homogeneous Heat Problems

MAT 167: Advanced Linear Algebra

= 1 u 6x 2 4 2x 3 4x + 5. d dv (3v2 + 9v) = 6v v + 9 3v 2 + 9v dx = ln 3v2 + 9v + C. dy dx ay = eax.

WMS. MA250 Introduction to Partial Differential Equations. Revision Guide. Written by Matthew Hutton and David McCormick

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

11.1 One-dimensional Helmholtz Equation

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0

THINKING IN PYRAMIDS

The end conditions on the beam are modelled by assuming the following boundary conditions on the function u: u(0, t) = u x (0, t) = 0;

Numerical methods for PDEs FEM - abstract formulation, the Galerkin method

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27

Problem set 6 The Perron Frobenius theorem.

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection

MAT 372 K.T.D.D. Final Sınavın Çözümleri N. Course. Question 1 (Canonical Forms). Consider the second order partial differential equation

Separation of Variables and a Spherical Shell with Surface Charge

EXISTENCE OF SOLUTIONS FOR ONE-DIMENSIONAL WAVE EQUATIONS WITH NONLOCAL CONDITIONS

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

Lecture notes. May 2, 2017

Lecture 6: Moderately Large Deflection Theory of Beams

Vibrations of Structures

Solution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima

Integrating Factor Methods as Exponential Integrators

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes...

A Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem

Lecture Note 3: Stationary Iterative Methods

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Physics 506 Winter 2006 Homework Assignment #6 Solutions

Problem (p.613) Determine all solutions, if any, to the boundary value problem. y + 9y = 0; 0 < x < π, y(0) = 0, y (π) = 6,

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

LECTURE 19: SEPARATION OF VARIABLES, HEAT CONDUCTION IN A ROD

Math 1600 Lecture 5, Section 2, 15 Sep 2014

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE

Mat 1501 lecture notes, penultimate installment

There are five problems. Solve four of the five problems. Each problem is worth 25 points. A sheet of convenient formulae is provided.

Physics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions

2.4 Eigenvalue problems

BALANCING REGULAR MATRIX PENCILS

v(x, 0) = g(x) where g(x) = f(x) U(x). The solution is where b n = 2 g(x) sin(nπx) dx. (c) As t, we have v(x, t) 0 and u(x, t) U(x).

BASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set

Wave Equation With Homogeneous Boundary Conditions

LIMITS AT INFINITY MR. VELAZQUEZ AP CALCULUS

SVM: Terminology 1(6) SVM: Terminology 2(6)

MATH-UA 263 Partial Differential Equations Recitation Summary

2M2. Fourier Series Prof Bill Lionheart

XSAT of linear CNF formulas

MATH 251 Final Examination May 3, 2017 FORM A. Name: Student Number: Section:

A Brief Introduction to Markov Chains and Hidden Markov Models

Transcription:

Math 0B - Summer 003 Homework Soutions Consider the eigenvaue probem { X = λx 0 < x < X satisfies symmetric BCs x = 0, Suppose f(x)f (x) x=b x=a 0 for a rea-vaued functions f(x) which satisfy the boundary conditions Show there are no negative eigenvaues Answer: Let X be an eigenfunction with eigenvaue λ λ X, X = X, X using the assumption that X X x= x=0 0 = X, X X X x= x=0 0 Consider the eigenvaue probem, { X = λx a < x < b X satisfies certain BC s Suppose µ is an eigenvaue of mutipicity m > Let X,, X m denote m ineary independent eigenfunctions (which may or may not be orthogona) associated with the eigenvaue µ Use these eigenfunctions to construct m eigenfunctions Y,, Y m which are necessariy orthogona Answer: First, define Y = X X Next, we define a function orthogona to Y by subtracting the piece of X which is parae to Y Therefore, et Z = X X, Y Y We normaize, etting Y = Z / Z Continuing in this way, we et Z 3 = X 3 X 3, Y Y X 3, Y Y and then define Y 3 = Z 3 / Z 3 In genera, we define Z m = X m X m, Y Y X m, Y m Y m, and et Y m = Z m / Z m

3 Consider the eigenvaue probem X = λx X (0) + X(0) = 0 X() = 0 0 < x < (a) Find an equation for the positive eigenvaues Answer: If λ = β > 0, then X + β X = 0 = X(x) = A cos(βx) + B sin(βx) X (0) + X(0) = Bβ + A = 0 X() = 0 = A cos(β) + B sin(β) = 0 Combining these two equations, we concude that λ is a positive eigenvaue if and ony if β cos(β) + sin(β) = 0 or tan(β) = β (b) Show graphicay that there are an infinite number of positive eigenvaues Answer: Let f(β) = β g(β) = tan(β) We know that g(β) has vertica asymptotes at β = ( n + ) π/ and g(β) + as β ( n + ) π/ By graphing both functions, we see that f(β) and g(β) have an infinite number of intersections (c) Show that λ = 0 is an eigenvaue if and ony if = Find a corresponding eigenfunction in this case Answer: λ = 0 = X = 0 = X(x) = A + Bx X (0) + X(0) = 0 = B + A = 0 X() = 0 = A + B = 0

Combining these two equations impies B( ) = 0 Therefore, zero is an eigenvaue if and ony if = In this case, the corresponding eigenfunction is X 0 (x) = B 0 ( x) for any B 0 (d) Show that if, then there are no negative eigenvaues, but if >, then there is one negative eigenvaue Find the corresponding eigenfunction Answer: λ = γ < 0 = X(x) = A cosh(γx) + B sinh(γx) X (0) + X(0) = 0 = Bγ + A = 0 X() = 0 = A cosh(γ) + B sinh(γ) = 0 Combining these two equations, we concude that Bγ cosh(γ) + B sinh(γ) = 0 If B = 0, then X(x) 0 Therefore, we must have or γ cosh(γ) + sinh(γ) = 0, tanh(γ) = γ To ook for roots of this equation, we consider the graphs of the functions f(γ) = γ g(γ) = tanh(γ) We note that f (γ) = and g (γ) = sech (γ) Therefore, f (0) =, g (0) = Further, we note that g (γ) = sech (γ) tanh(γ) < 0 Aso, we note that g(γ) as γ + From this information, we concude that f(γ) and g(γ) ony intersect at γ = 0 if, and, consequenty, there are no negative eigenvaues if On the other hand, if >, f(γ) and g(γ) wi intersect once, impying there is exacty one negative eigenvaue if > 3

4 Use separation of variabes to sove u t ku xx = 0 0 < x <, t > 0 u(x, 0) = φ(x) u(0, t) = 0 u x (, t) = 0 Answer: We use separation of variabes, etting u(x, t) = X(x)T (t) Pugging this function into our PDE, we have kt = X X = λ We are ed to the eigenvaue probem { X = λx 0 < x < X(0) = 0 = X () Looking for positive eigenvaues, Therefore, T λ = β > 0 = X(x) = A cos(βx) + B sin(βx) X(0) = 0 = A = 0 X () = 0 = Bβ cos(β) = 0 = β n = ( n + ) π (( ) (( ) n + λ n = ) π n + X n (x) = B n sin ) π x By straightforward cacuation, we see these are a the eigenvaues Soving the equation for T n, Therefore, we et T = kλ n T = T n (t) = C n e kλnt u(x, t) = = T n (t) = C n e k (n+ (( ) n + C n sin ) π x e n= )π k (n+ t ) π! t In order for the initia condition u(x, 0) = φ(x) to be satisfied, we need ( ) (n+ sin )π x, φ(x) C n = ( ) ( ) (n+ (n+ sin x, sin x )π )π 4

5 Use separation of variabes to sove u t ku xx + u x = 0 0 < x <, t > 0 u(x, 0) = φ(x) u(0, t) = 0 u x (0, t) = 0 (Hint: Introduce a function f(x) such that v(x, t) = f(x)u(x, t) wi satisfy a PDE of the form v t kv xx + av = 0 with new initia and boundary conditions soution, sove for u) Answer: We wi give two methods Sove the equation for v, and from this Method We want to find a function f such that defining v = fu, then v satisfies a PDE of the form v t kv xx + av = 0 Let g = /f Therefore, u = gv Pugging this function into our PDE, we have gv t kgv xx + [ g + g]v x + [ kg + g ]v = 0 In order to eiminate the v x term, we want to choose g such that g + g = 0, which impies g(x) = e x/ Then et v = fu = (/g)u = e x/ u Soving for kg + g, we see that v wi satisfy the foowing initia/boundary-vaue probem v t kv xx + 4k v = 0 0 < x <, t > 0 v(x, 0) = e x/ φ(x) v(0, t) = 0 v x (, t) + v(, t) = 0 Now use separation of variabes on this initia/boundary-vaue probem Pugging in a function v(x, t) = X(x)T (t) impies We need to sove the ODE T kt = X X 4k = λ X X = λ + 4k Let µ = λ + We are ed to the eigenvaue probem 4k X = µx X(0) = 0 X () + X() = 0 0 < x < If µ = β > 0, then X(x) = A cos(βx) + B sin(βx) X(0) = 0 = A = 0 5

X () + X() = 0 = Bβ cos(β) + B sin(β) = 0 If B = 0, then X(x) 0 Therefore, we need or The corresponding eigenfunctions are β cos(β) + sin(β) = 0, tan(β) = β X n (x) = B n sin(β n x) If µ = 0, then X(x) = A + Bx X(0) = 0 = A = 0 X () + X() = 0 = B + B = 0 = B = 0 Therefore, zero is not an eigenvaue Using the resut from probem (), we can concude that there are no negative eigenvaues Recaing the definition of µ, we see that where Soving our equation for T n, we have λ n = β n + 4k tan(β n ) = β n T n (t) = C n e kλnt Therefore, we et v(x, t) = C n sin(β n x)e kλ nt n= where λ n and β n are defined above We want v(x, 0) = e x/ φ(x) Therefore, et C n = sin(β nx), e x/ φ(x) sin(β n x), sin(β n x) 6

Finay, we reca that u(x, t) = e x/ v(x, t) for v defined above Method Here is an aternate method, which aows one to use separation of variabes directy, but requires knowedge about Sturm-Liouvie probems (You are not responsibe for this materia, but I thought I d write up a soution using this method because some students used this technique) Let u = XT Pugging u into our equation, we have T kt = X k X X = λ We need to consider the eigenvaue probem { X + k X λx = 0 0 < x < X(0) = 0 = X () Mutipying by the integrating factor e x/k, we see this ODE can be rewritten as ( e x/k X ) λe x/k X = 0, an exampe of a Sturm-Liouvie probem We note that eigenfunctions of this Sturm- Liouvie probem are orthogona with respect to the weight function e x/k To find the soutions of our second-order ODE, X k X + λx = 0, we ook at the characteristic poynomia p(r) = r r + λ = 0 The roots of this k poynomia are given by ± 4λ k k Let µ = 4k λ Then r = The first case is µ > 0 In this case, = ± 4k λ r = ± µ r = ± i µ, which impies the associated eigenfunctions are given by X(x) = e x/ [A cos( µx) + B sin( µx)] 7

X(0) = 0 = A = 0 X () = 0 = [ sin( µ) + µ cos( µ)] = 0 = tan( µ) = µ Thus far, we have eigenvaues where and the associated eigenfunctions λ n = 4k + µ n tan( µ n ) = µ n X n (x) = e x/ sin( µ n x) Next, we consider µ = 0 In this case, r = which impies X(x) = Ae x/ + Bxe x/ X () = 0 = B X(0) = 0 = A = 0 [ e x/ + ] e/ = 0 But this can ony happen if B = 0 which impies X(x) 0 Therefore, µ 0 Next, we consider µ < 0 In this case, r = ± µ, which impies [ X(x) = e x/ Ae ] µx + Be µx X(0) = 0 = A + B = 0 X () = 0 = ( + ) µ e µ = µ or B = 0 Since both terms in the expression above are greater than, we must have B = 0 Therefore, µ 0 We concude that λ = 4k + µ 8

where µ > 0 such that and tan( µ) = µ X n (x) = e x/ sin( µ n x) The soution for our equation for T is given by T n (t) = C n e kλnt Therefore, et u(x, t) = C n e x/ sin( µ n x)e kλnt n= We want u(x, 0) = φ(x) That is, we want to choose C n such that C n e x/ sin( µ n x) = φ(x) n= Using the fact that eigenfunctions are orthogona with respect to the weight function e x/, we can mutipy both sides by sin( µ m x) and integrate over [0, ] with respect to the weight function e x/ In particuar, we concude that C n = 0 φ(x)ex/ sin( µ n x)e x/k dx 0 (ex/ sin( µ n x)) e x/k dx Simpifying, we have C n = 0 φ(x)e x/ sin( µ n x) dx 0 sin ( µ n x) dx 9