Fourier series: Fourier, Dirichlet, Poisson, Sturm, Liouville

Similar documents
An improper integral and an infinite series

Four Mathema,cians Who Shaped Our Understanding of Calculus

Stories from the Development Of Real Analysis

Chapter 2 Joseph Fourier

O1 History of Mathematics Lecture VIII Establishing rigorous thinking in analysis. Monday 31st October 2016 (Week 4)

PHYS 502 Lecture 3: Fourier Series

Group representations and harmonic

Tyn Myint-U Lokenath Debnath. Linear Partial Differential Equations for Scientists and Engineers. Fourth Edition. Birkhauser Boston Basel Berlin

Lecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.

Fourier Series and Recent Developments in Analysis

Beyond Newton and Leibniz: The Making of Modern Calculus. Anthony V. Piccolino, Ed. D. Palm Beach State College Palm Beach Gardens, Florida

Fourier Analysis Partial Differential Equations (PDEs)

MATH-UA 263 Partial Differential Equations Recitation Summary

O1 History of Mathematics Lecture VIII Establishing rigorous thinking in analysis. Monday 30th October 2017 (Week 4)

Mathematics 324 Riemann Zeta Function August 5, 2005

A glimpse of Fourier analysis

ADVANCED ENGINEERING MATHEMATICS MATLAB

Topics in Fourier analysis - Lecture 2.

An Introduction to Partial Differential Equations in the Undergraduate Curriculum

Application of Series in Heat Transfer - transient heat conduction

G: Uniform Convergence of Fourier Series

1 Distributions (due January 22, 2009)

Linear Partial Differential Equations for Scientists and Engineers

1 Solutions in cylindrical coordinates: Bessel functions

5.4 Continuity: Preliminary Notions

Physics 250 Green s functions for ordinary differential equations

Convergence of Fourier series

Ordinary Differential Equations II

On Exponential Decay and the Riemann Hypothesis

Mathematical Methods: Fourier Series. Fourier Series: The Basics

Title Project Summary

Math 46, Applied Math (Spring 2009): Final

Mathematical Modeling using Partial Differential Equations (PDE s)

MATH 104 : Final Exam

PARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS

Primes in arithmetic progressions

INDEX. Bolzano-Weierstrass theorem, for sequences, boundary points, bounded functions, 142 bounded sets, 42 43

9.2 Taylor Polynomials

PART IV Spectral Methods

1.1 Introduction to Limits

The Way of Analysis. Robert S. Strichartz. Jones and Bartlett Publishers. Mathematics Department Cornell University Ithaca, New York

HARMONIC ANALYSIS TERENCE TAO

INTRODUCTION TO REAL ANALYSIS II MATH 4332 BLECHER NOTES

Divergent Series: why = 1/12. Bryden Cais

be the set of complex valued 2π-periodic functions f on R such that

On the computational complexity of mathematical functions

u(e iθ )K(e iθ,z)dθ, z (1). (126) h(e iθ )K(e iθ,z)dθ, z (1),

Complex Analysis: A Round-Up

More Series Convergence

Outline of Fourier Series: Math 201B

Differential Equations

The lecture of 1/23/2013: WHY STURM-LIOUVILLE?

The Prime Number Theorem

Waves on 2 and 3 dimensional domains

Final Exam May 4, 2016

CHAPTER 6 bis. Distributions

= 0. Theorem (Maximum Principle) If a holomorphic function f defined on a domain D C takes the maximum max z D f(z), then f is constant.

ECE Digital Image Processing and Introduction to Computer Vision. Outline

KREYSZIG E Advanced Engineering Mathematics (10th ed., Wiley 2011) Chapter 11 - Fourier analysis

A proof for the full Fourier series on [ π, π] is given here.

The shortest route between two truths in the real domain passes through the complex domain. J. Hadamard

AP Calculus BC. Course Overview. Course Outline and Pacing Guide

Mathematics (MA) Mathematics (MA) 1. MA INTRO TO REAL ANALYSIS Semester Hours: 3

APPLIED PARTIAL DIFFERENTIAL EQUATIONS

Two special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p

JUHA KINNUNEN. Partial Differential Equations

MATHEMATICS. Course Syllabus. Section A: Linear Algebra. Subject Code: MA. Course Structure. Ordinary Differential Equations

Let R be the line parameterized by x. Let f be a complex function on R that is integrable. The Fourier transform ˆf = F f is. e ikx f(x) dx. (1.

Bernoulli Numbers and their Applications

Decoupling course outline Decoupling theory is a recent development in Fourier analysis with applications in partial differential equations and

Chapter 17: Fourier Series

An idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim

Let s Get Series(ous)

Fourier and Partial Differential Equations

On Exponential Decay and the Riemann Hypothesis

SAMPLE FINAL EXAM SOLUTIONS

BASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA

Additional Homework Problems

David M. Bressoud Macalester College, St. Paul, Minnesota Given at Allegheny College, Oct. 23, 2003

Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series

Fourier Series. 1. Review of Linear Algebra

MATH 251 Final Examination December 19, 2012 FORM A. Name: Student Number: Section:

Prime Number Theory and the Riemann Zeta-Function

14 Fourier analysis. Read: Boas Ch. 7.

Advanced Calculus Math 127B, Winter 2005 Solutions: Final. nx2 1 + n 2 x, g n(x) = n2 x

MULTIPOLE EXPANSIONS IN THE PLANE

Primes. Peter Martinson

21 Laplace s Equation and Harmonic Functions

MATH 5640: Fourier Series

CS711008Z Algorithm Design and Analysis

Contents. Preface. Notation

An Examination of Richard Dedekind s Continuity and Irrational Numbers

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.

CONVERGENCE OF THE FOURIER SERIES

n 2 xi = x i. x i 2. r r ; i r 2 + V ( r) V ( r) = 0 r > 0. ( 1 1 ) a r n 1 ( 1 2) V( r) = b ln r + c n = 2 b r n 2 + c n 3 ( 1 3)

Alexander Ostrowski

FOURIER TRANSFORMS. 1. Fourier series 1.1. The trigonometric system. The sequence of functions

Polynomials. p n (x) = a n x n + a n 1 x n 1 + a 1 x + a 0, where

Functional Limits and Continuity

MATH 5400, History of Mathematics

Transcription:

Fourier series: Fourier, Dirichlet, Poisson, Sturm, Liouville Joseph Fourier (1768-1830) upon returning from Egypt in 1801 was appointed by Napoleon Prefect of the Department of Isères (where Grenoble is the main city). He served in this office for 14 years. Beginning in 1804 he investigated the equation of heat propagation, and submitted in 1807 a memoir to the Academy of Sciences with his results, where it was received with skepticism. Laplace, Lagrange, Poisson and Monge reviewed it. Laplace, for example, could not believe cos x could be expressed as a sine series, while Laplace and Poisson thought the derivation of the equation lacked in rigor. The memoir was resubmitted successfully in 1811 (for a competition devised for this purpose; Fourier got the prize), and eventually published in 1822. The main subject is the equation for the temperature u(x, t) in a onedimensional medium: u t = k 2 u t 2 u(x, 0) = f(x) given, either on the whole real line or on an interval 0 < x < L with suitable boundary conditions. The solution is given in terms of Fourier series: u(x, t) = e ( 4π2 k L n2t) (a n cos( 2nπx L ) + b n sin( 2nπx L )) n=0 (for periodic boundary conditions), where: f(x) a 0 2 + n=1 a n 2 cos(2nπx L ) + b n 2 sin(2nπx L ), with the expression for the coefficients (where we assume L = 2π and f is periodic with period 2π): a n = 1 π 2π 0 f(x) cos(nx)dx, b n = 1 π 2π 0 f(x) sin(nx)dx). (Incidentally, the notational idea of attaching the limits of integration to the bottom and top of Leibniz s sign first appeared in Fourier s work.) The question of representability of a continuous function by its Fourier series was considered by Fourier in his treatise, and he believed he had demonstrated it. The key point in the argument was his claim that, for continuous f and infinitely large p one has: sin[p(x y)]f(y)dy = 0. 1

The idea is that the contribution between two zeros of the sine function is cancelled by the contributions between the next two zeros, since the distance between zeros is infinitely small and the function f (being continuous) is essentially constant over such an interval. This argument is not rigorous, and indeed the result is false: Du-Bois Raymond in 1876 gave an example of a continuous 2π-periodic function with the property that a sequence of partial sums S n (x) at zero of its Fourier series tends to infinity: lim S n i (0) =, n i. i The first (and still the most useful) satisfactory result on the representation of periodic functions by Fourier series is due to P.G.Lejeune Dirichlet (1805-1859), in 1829: Theorem: Assume f(x) is 2π-periodic, bounded and continuous with bounded, continuous first derivative except at a finite number of points in [0, 2π], where the one-sided limits exist. Then S n (x) f(x), except where f is discontinuous (where it converges to the average of the one-sided limits.) Dirichlet had spent the years 1822-1827 in Paris (supporting himself by tutoring), and during this time he had close contact with Fourier s school. Even as a student he must have made an impression, for on his return he received teaching appointments at Breslau (1827) and Berlin (1829), on Alexander v. Humboldt s suggestion. In the years 1847-1849, in Berlin, he taught Bernhard Riemann (1826-1866), a visiting doctoral student from Göttingen. In 1855, Dirichlet was appointed to succeed Gauss in Göttingen. In addition to his fundamental contributions to analysis, Dirichlet also proved a major result in analytic number theory: in a paper of 1840 he solved a fifty-year old conjecture of Legendre: if a and d are co-prime, the arithmetic progression with first term a and differences d contains infinitely many primes; the proof combines a generalization of Euler s product formula with a generalized zeta function, henceforth known as Dirichlet L-function. In 1833, George Green (1793-1841) first used what came to be known as the Dirichlet principle to give a variational solution to a problem in potential theory: for suitable boundary conditions, the variational problem in a domain D R n : minimize u 2 dx has a solution. This principle was used by Dirichlet in potential theory, and later by Riemann for existence proofs in complex function theory; it was 2 D

later attacked by an example of Weierstrass (1870), and eventually rescued by Hilbert (1899). Dirichlet was Riemann s teacher, and Riemann was very interested in the question of necessary conditions for representability by Fourier series. Riemann s investigations were published in the 1853 Habilitation paper On the representability of a function by means of a trigonometric series, and led naturally to Riemann s formulation of what we now know as the Riemann integral (see [Dunham-Gallery]), as follows. The key issue is whether: 2π lim f(x) sin(nx)dx = 0. n 0 (now known as Riemann-Lebesgue lemma.) Riemann argues we can break up the integral into intervals of length 2π/n, and observe that the sine is positive in the first half of such an interval, negative in the second half; so that if in the k th. interval: m k f(x) M k and hence M k f(x) m k, the integral over the two halves is less than 2 M k m k n, that is: less than the length of the interval times the oscillation of the function over all intervals. Thus if the sum over all intervals of these products tends to zero as n (the Riemann integrability condition), we have that the desired limit is also zero (that is, a n, b n 0.) Note that this is a rigorization (in the spirit of Dirichlet, Cauchy and Weierstrass) of Fourier s intuitive idea. Also in Fourier s work one saw for the first time the Fourier transform : F[f](s) = ˆf(s) = f(t)e ist dt, (defined, say, for f continuous with f integrable over the whole real line.) The main application given by Fourier was to solve the heat equation over the whole real line, with given initial condition f(x). The solution formula is: u(x, t) = 1 4πkt f(x y)e y2 4kt dt; the right-hand side is the convolution integral of f and the heat kernel : E 1/ 2kt (x) = 1 4πkt e x2 4kt, EK (x) = 3 1 K 2π e K 2 x 2 2

The solution to this problem had previously been given by Laplace (following work of Poisson), and Fourier showed how it would also follow from his methods. If f is bounded and uniformly continuous, one indeed has u(x, t) f(x) as t 0+, uniformly in x. The use of such uniformity conditions in analysis was emphasized by Karl Weierstrass (1815-1897), who used Fourier s solution to prove an important theorem: Theorem (Weierstrass density theorem (1885)). Any continuous function f(x) on a closed interval [a, b] can be approximated arbitrarily closely by polynomials; that is, given ɛ > 0 one may find a polynomial p(x) so that: max x [a,b] f(x) p(x) < ɛ. Here is the general idea of the proof: using the standard power series for the exponential, one sees that, for n large, the heat kernel E K (x) is arbitrarily close to the polynomial (uniformly on a large interval [ R, R]): Q 2n (x) = (1/k 2π) n ( k 2 x 2 /2) r /r!. Now, as just seen, f(x) is uniformly close (in [a, b]) to the solution u(x, t), for t sufficiently close to zero, and this solution is the convolution of (a suitable extension to the real line of) f with E 1/ 2kt. Finally, it is easy to see that the convolution of the extended f with a polynomial is itself a polynomial. One of the problems Fourier studies in his memoir is the cooling of a sphere of radius R, with heat loss at the boundary proportional to the excess temperature: u r r=r = hu(r, t). This gives rise to the rotationally symmetric heat equation: r=0 1 u k t = 2 u r 2 + 2 u r r, and by separation of variables, eventually to an eigenvalue problem : ρ + 2 r ρ + λ 2 ρ = 0, with suitable boundary conditions. It turns out the possible values of λ are solutions of the equation: Rλ tan Rλ = 1 Rh. 4

Fourier was not able to show that there are no complex solutions λ to this equation, as pointed out by Siméon-Denis Poisson (1781-1840), in a paper of 1823 including an alternative derivation of the heat equation, along lines suggested by Laplace (1809). A few years later (1826), Poisson found a proof that the eigenfunctions (of this and similar problems) are orthogonal, and the eigenvalues are real- exactly the same proof presented in courses today. Poisson wrote his own Théorie analytique de la chaleur (1835), where he introduced the multidimensional version of these orthogonality results. The general theory of eigenvalue problems of this type (for second-order ordinary differential equations) was developed by Joseph Liouville (1809-1882) and Jacques-Charles-Francois Sturm (1803-1855), in a series of papers in the 1830s. A little more precisely, two major memoirs by Sturm (1836) developed completely the oscillation theorems and spectral theory for the whole class of Sturm-Liouville eigenvalue problems, generalizing the work of Fourier and Poisson. They were followed by three papers by Liouville (in 1836 and 1837), establishing the convergence of the associated generalized Fourier series, and by a joint paper in 1837 (Sturm and Liouville were friends; the latter supported the election of the former, who was Swiss, to the French Academy.) Poisson made many important contributions to the use of trigonometric series in problems of mathematical physics (and also to analytical mechanics and probability.) For example, the method of Poisson summation obtains the representation of a 2π- periodic function by its (complex) Fourier series: f(x) n= c n e inx as the boundary limit of functions defined in the unit disk. Poisson showed that: g(r, x) := c n r n e inx n= converges uniformly in x for each 0 < r < 1; and: g(r, x) f(x) as r 1, uniformly in x. In the process, he discovered the famous Poisson kernel : P r (x) = r n e inx = n= 5 1 r 2 1 2r cos x + r 2,

which is used to solve the Dirichlet problem (find a harmonic function with given boundary values) for the unit disk. Poisson s result shows that: lim lim N r 1 N n= N c n r n e inx exists, with value f(x); this is not the same as saying that lim lim N N r 1 n= N c n r n e inx exists, for arbitrary continuous f; indeed this is false, by the du Bois Raymond example mentioned earlier. Poisson failed to appreciate this subtle point, and believed he had proved the convergence of all Fourier series (something that also Dirichlet expected to be true.) In the period between Cauchy and Weierstrass, the idea of rigor in Analysis was in a state of flux, as evidenced by this quote, from a letter of Jacobi to Alexander v. Humboldt, December 1846 (quoted in [Laugwitz, p.63]): Only Dirichlet-not I, not Cauchy, not Gauss- knows what a perfectly rigorous proof is; we learn it only from him. When Gauss says he has proved something, I think it very likely; when Cauchy says it, it is a fifty-fifty bet; when Dirichlet says it, it is certain. I d rather not go into these delicate matters... SOURCES: 1. T.W.Körner, Fourier Analysis, Cambridge U.P.1988; 2. D.Laugwitz, Bernhard Riemann- Turning Points in the Conception of Mathematics, Birkhäuser Boston 1999 (Abe Shenitzer, transl.)3. J.Lützen, Joseph Liouville-Master of Pure and Applied Mathematics, Springer-Verlag 1990. 6