Math 124B January 17, 2012

Similar documents
Math 124B January 31, 2012

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

14 Separation of Variables Method

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

6 Wave Equation on an Interval: Separation of Variables

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.

MA 201: Partial Differential Equations Lecture - 10

Wave Equation Dirichlet Boundary Conditions

2M2. Fourier Series Prof Bill Lionheart

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

4 1-D Boundary Value Problems Heat Equation

Lecture Notes 4: Fourier Series and PDE s

Assignment 7 Due Tuessday, March 29, 2016

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

Math 220B - Summer 2003 Homework 1 Solutions

4 Separation of Variables

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

Week 6 Lectures, Math 6451, Tanveer

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

1 Heat Equation Dirichlet Boundary Conditions

MA 201: Partial Differential Equations Lecture - 11

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

Week 5 Lectures, Math 6451, Tanveer

221B Lecture Notes Notes on Spherical Bessel Functions

C. Fourier Sine Series Overview

1D Heat Propagation Problems

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

b n n=1 a n cos nx (3) n=1

Higher dimensional PDEs and multidimensional eigenvalue problems

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

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

AMATH 353 Lecture 9. Weston Barger. How to classify PDEs as linear/nonlinear, order, homogeneous or non-homogeneous.

Math 124B January 24, 2012

David Eigen. MA112 Final Paper. May 10, 2002

Chapter 5. Wave equation. 5.1 Physical derivation

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

Mathematical Methods: Fourier Series. Fourier Series: The Basics

22. Periodic Functions and Fourier Series

Suggested Solution to Assignment 5

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

Vibrations of Structures

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

Introduction and preliminaries

Legendre Polynomials - Lecture 8

More Scattering: the Partial Wave Expansion

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

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

An approximate method for solving the inverse scattering problem with fixed-energy data

Method of Separation of Variables

Physics 505 Fall Homework Assignment #4 Solutions

Srednicki Chapter 51

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

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

Transforms and Boundary Value Problems

Separation of Variables and a Spherical Shell with Surface Charge

CS229 Lecture notes. Andrew Ng

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

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].

Problem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.

Integrating Factor Methods as Exponential Integrators

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

Wavelet Galerkin Solution for Boundary Value Problems

XSAT of linear CNF formulas

MAT 167: Advanced Linear Algebra

Lecture notes. May 2, 2017

11.1 One-dimensional Helmholtz Equation

Lecture Note 3: Stationary Iterative Methods

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.

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

V.B The Cluster Expansion

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

arxiv: v1 [math.fa] 23 Aug 2018

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),

Some Measures for Asymmetry of Distributions

A Brief Introduction to Markov Chains and Hidden Markov Models

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Homework 5 Solutions

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

Lecture 6: Moderately Large Deflection Theory of Beams

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

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

V.B The Cluster Expansion

The Wave Equation. Units. Recall that [z] means the units of z. Suppose [t] = s, [x] = m. Then u = c 2 2 u. [u]

Section 6: Magnetostatics

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7

= 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.

Stat 155 Game theory, Yuval Peres Fall Lectures 4,5,6

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons

arxiv: v1 [math.ap] 8 Nov 2014

Eigenvalue Problem. 1 The First (Dirichlet) Eigenvalue Problem

FOURIER SERIES ON ANY INTERVAL

Midterm 2: Sample solutions Math 118A, Fall 2013

Haar Decomposition and Reconstruction Algorithms

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

A Two-Parameter Trigonometric Series

Fourier series. Part - A

B. Brown, M. Griebel, F.Y. Kuo and I.H. Sloan

Overview of Fourier Series (Sect. 6.2). Origins of the Fourier Series.

Partial Differential Equations Summary

Transcription:

Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia data. In the same way the periodic boundary conditions u(, t) = u(, t), u x (, t) = u x (, t), (1) ead to the fu Fourier expansion, i.e. the Fourier series that contains both sines and cosines. Such boundary conditions arise naturay in appications. For exampe heat conduction phenomena in a circuar rod wi be described by the heat equation subject to boundary conditions (1). To see how the fu Fourier series comes about, et us consider the eigenvaue probem associated with the boundary conditions (1), X = λx X( ) = X(), (2) X ( ) = X (). To find the eigenvaues and eigenfunctions for (2), we consider the quaitativey distinct cases λ <, λ = and λ > separatey. We first assume λ = γ 2 <. In this case the equation takes the form X = γ 2 X, and the soutions are thus X(x) = Ce γx + De γx. The boundary conditions then impy Ce γ + De γ = Ce γ + De γ, Cγe γ Dγe γ = Cγe γ Dγe γ. Mutipying the first equation by γ and adding to the second equation gives 2Cγe γ = 2Cγe γ 2Cγe γ (1 e 2γ ) = C =, since γ,. Simiary, D =, which impies that there are no negative eigenvaues. We next assume that λ =, which resuts in the equation X =. The soution is then X(x) = C+Dx, and the boundary conditions impy C D = C + D 2D = D =. D = D So X(x) 1 is an eigenfunction corresponding to the eigenvaue λ =. Finay, we consider the case λ = β 2 >. The equation then takes the form X = β 2 X, and hence, has the soution X(x) = C cos βx + D sin βx. Checking the boundary conditions gives C cos β D sin β = C cos β + D sin β Cβ sin β + Dβ cos β = Cβ sin β + Dβ cos β 2D sin β = 2Cβ sin β = Since C and D cannot both be equa to zero, and β, we must have sin β = β = nπ, for n = 1, 2,.... Then the positive eigenvaues and the corresponding eigenfunctions are λ n = ( nπ ) 2, Xn (x) = C + D sin nπx, for n = 1, 2,.... Notice that for each of the eigenvaues λ n, n = 1, 2,... we have two ineary independent eigenfunctions, cos(nπx/) and sin(nπx/). This differs from the case of the Dirichet and Neumann boundary conditions, where we had ony one ineary independent eigenfunction for each of the same eigenvaues, 1

namey sin(nπx/) for Dirichet, and cos(nπx/) for Neumann. We wi see in the next ecture that in the cases where there are more than one ineary independent eigenfunctions corresponding to the same eigenvaue, one can aways choose a pairwise orthogona set of eigenfunctions, which is necessary for the method of computing the Fourier coefficients to go through. Having soved the eigenvaue probem (2), it is now cear that the series soutions of PDEs subject to the periodic boundary conditions (1) wi have both sines and cosines in their expansions. For exampe the soution to the heat equation satisfying the boundary conditions (1) wi be given by the series u(x, t) = A e (nπ/)2kt (A n + B n sin nπx ), where the coefficients A n, B n come from the expansion of the initia data φ(x) = A A n + B n sin nπx. (3) The series (3) is caed the fu Fourier series of the function φ(x) on the interva (, ). It is now cear that to sove boundary vaue probems with periodic boundary conditions (1) via separation of variabes, one needs to find the coefficients in the expansion (3). This procedure is simiar to finding the coefficients of the Fourier cosine and sine series, and reies on the pairwise orthogonaity of the eigenfunctions. Notice that (3) is a inear combination of the eements of the set 1,, sin nπx : n = 1, 2,... }. (4) Now, if we can prove that the eements of this set are pairwise orthogona, then the coefficients in (3) can be found by taking the dot product of the series (3) with the corresponding eigenfunction. Since the interva in this case is (, ), we define the dot product of two functions to be (f, g) = ˆ f(x)g(x) dx. Let us now check the orthogonaity of the set (4). We need to show that ˆ ˆ ˆ ˆ ˆ sin nπx 1 1 sin mπx sin mπx sin mπx dx =, for a n, m = 1, 2,..., (5) dx =, for n m, (6) dx =, for n m, (7) dx =, for a m = 1, 2,..., (8) dx =, for a m = 1, 2,.... (9) Notice that, since sin(nπx/) and cos(nπx/) sin(mπx/) are odd functions, (5) and (9) are trivia, since the integration is performed over a symmetric interva. The integrands in (6), (7) and (8) are even functions, so the integras are equa to twice the corresponding integras over the interva (, ). Then, say for (6), we have ˆ ˆ dx = 2 dx =, for n m, 2

as was computed in the ast ecture, when considering Fourier cosine series over the interva (, ). Proofs of(7) and (8) are simiar. Thus, the set (4) is indeed orthogona. To compute the coefficients in (3), we take the dot product of the series with the corresponding eigenfunctions. For exampe, ˆ φ(x) dx = ˆ ( A 2 + A n + B n sin nπx ) and simiary for the eigenfunctions 1 = cos(πx/), and sin(mπx/). But ˆ ˆ ˆ sin 2 mπx 1 2 dx = 2. So the coefficients can be computed as ˆ dx = 2 dx = 2 ˆ sin 2 mπx dx = 2 2 =, dx = 2 2 =, = A m ˆ dx, A m = 1 B m = 1 ˆ ˆ φ(x) φ(x) sin mπx dx, n =, 1, 2,..., dx, n = 1, 2,.... (1) Exampe 3.1. Find the fu Fourier series of the function φ(x) = x on the interva (, ). We compute the coefficients in the expansion using formuas (1). x = A A n A n = 1 ˆ x + B n sin nπx dx =, since the function x cos(nπx/) is odd, and the integration interva is symmetric. Using evenness of the function x sin(nπx/), we have B n = 1 ˆ x sin nπx dx = 2 ˆ x sin nπx n+1 2 dx = ( 1) nπ, where the ast integra is exacty the coefficient of the Fourier sine series of function φ(x) = x, computed in the ast ecture. Thus the fu Fourier series of x on the interva (, ) is x = ( 1) n+1 2 nπ nπx sin, which exacty coincides with the Fourier sine series of the function x on the interva (, ). 3

3.1 Even, odd and periodic functions In the previous exampe we coud take any odd function φ(x), and the coefficients of the cosine terms in the fu Fourier series woud vanish for exacty the same reason, eading to the Fourier sine series. Thus the fu Fourier series of an odd function on the interva (, ) coincides with its Fourier sine series on the interva (, ). Anaogousy, the fu Fourier series of an even function on the interva (, ) coincides with its Fourier cosine series on the interva (, ). This is not surprising, taking into account the fact that the sine and cosine functions are themseves respectivey odd and even. So how can the same function φ(x), be it odd, even, or neither, have both a nonzero Fourier cosine and a nonzero Fourier sine expansion on the interva (, )? The answer to this question is trivia, if one recas the refection method used to sove boundary vaue probems on the finite interva (, ). In the refection method one takes either an odd (Dirichet BCs), or an even (Neumann BCs) extension of the data. Then the resuting extended data is either odd or even on the symmetric interva (, ), and is extended to be periodic with period 2 to the rest of the number ine. But then one can find the fu Fourier series of this extended data, which wi contain ony sines or cosines, depending on the evenness or the oddness of the extension, thus giving either the Fourier cosine, or Fourier sine series of the function φ(x) on the interva (, ). In short, one has the foowing correspondence between the boundary conditions, extensions, and Fourier series. Dirichet: u(, t) = u(, t) = odd extension Fourier sine series. Neumann: u x (, t) = u x (, t) = even extension Fourier cosine series. Periodic: u(, t) = u(, t), u x (, t) = u x (, t) periodic extension fu Fourier series. 3.2 The compex form of the fu Fourier series Instead of writing the Fourier series in terms of cosines and sines, one can use compex exponentia functions. Reca the Euer s formua, e iθ = cos θ + i sin θ, (11) from which one can get the formua for e iθ in terms of cosine and sine by pugging in θ. Soving from these two equations, we have the expressions of sine and cosine in terms of the compex exponentias cos θ = eiθ + e iθ 2, sin θ = eiθ e iθ. (12) 2i So instead of (3), one can write the Fourier expansion in terms of the functions e inπx/ and e inπx/, themseves eigenfunctions corresponding to the eigenvaue λ n = n 2 π 2 / 2. The convenience of the compex form is in that the Fourier series can be written compacty as φ(x) = n= C n e inπx/, (13) which is a inear combination in terms of the eements of the set 1, e iπx/, e iπx/, e i2πx/, e i2πx/,... } = 1, e ±inπx/ : n = 1, 2,... }, (14) repacing the set (4). In order to find the coefficients in the expansion (13), we first need to show that the set (14) is orthogona, after which the coefficients can be found by taking the dot product of the series with the corresponding eigenfunction. Notice that the functions in the set (14) are compex vaued. For compex vaued functions the dot product is defined as (f, g) = ˆ f(x)g(x) dx, 4

where the bar stands for the compex conjugate. Then the dot product of two distinct eements of the set (14) is ˆ e inπx/ e imπx/ dx = 1 i(n m)π ei(n m)πx/ = e i(n m)π i(n m)π (ei2(n m)π 1) =, for n m. So taking the dot product of the series (13) with one of the eigenfunctions in (14) gives ˆ φ(x)e imπx/ dx = ˆ n= Thus, the coefficients can be computed as 3.3 Concusion ˆ C n e inπx/ e imπx/ dx = C m e imπx/ e imπx/ dx = C m 2. C n = 1 2 ˆ φ(x)e inπx/ dx. (15) We saw how the separation of variabes in the case of periodic boundary conditions eads to the fu Fourier series (3) in the same way as the Dirichet and Neumann boundary conditions ead to the Fourier sine and cosine series respectivey. The coefficients in the fu Fourier expansion were computed by making use of the pairwise orthogonaity of eigenfunctions, resuting in the coefficients formuas (1). We aso expored the connection between the Fourier cosine, respectivey sine, expansions on the interva (, ) and the fu Fourier expansion on the interva (, ), using the idea of even and odd extensions. Finay using Euer s formua (11), we rewrote the fu Fourier series in the compex form (13), the coefficients of which were again computed by taking the dot product of the series with the eigenfunctions, and reying on the fact that the distinct eigenfunctions are orthogona. We wi see in the next ecture that one aways has orthogona eigenfunctions, provided the boundary conditions are symmetric in an appropriate sense. This wi ead to the notion of generaized Fourier series. 5