CHAPTER 10 Fourier Series In the late eighteenth century, it was well-known that complicated functions could sometimes be approximated by a sequence of polynomials. Some of the leading mathematicians at that time, including such luminaries as Daniel Bernoulli, Euler and d’Alembert began studying the possibility of using sequences of trigonometric functions for approximation. In 1807, this idea opened into a huge area of research when Joseph Fourier used series of sines and cosines to solve several outstanding partial di↵erential equations of physics.1 In particular, he used series of the form 1 X an cos nx + bn sin nx n=0 to approximate his solutions. Series of this form are called trigonometric series, and the ones derived from Fourier’s methods are called Fourier series. Much of the mathematical research done in the nineteenth and early twentieth century was devoted to understanding the convergence of Fourier series. This chapter presents nothing more than the tip of that huge iceberg. 1. Trigonometric Polynomials Definition 10.1. A function of the form (111) p(x) = n X ↵k cos kx + k sin kx k=0 is called a trigonometric polynomial. The largest value of k such that |↵k |+ k | = 6 0 is the degree of the polynomial. Denote by T the set of all trigonometric polynomials. Evidently, all functions in T are 2⇡-periodic and T is closed under addition and multiplication by real numbers. Indeed, it is a real vector space, in the sense of linear algebra and the set {sin nx : n 2 N} [ {cos nx : n 2 !} is a basis for T . The following theorem can be proved using integration by parts or trigonometric identities. Theorem 10.2. If m, n 2 Z, then Z ⇡ (112) sin mx cos nx dx = 0, ⇡ 1 Fourier’s methods can be seen in most books on partial di↵erential equations, such as [3]. For example, see solutions of the heat and wave equations using the method of separation of variables. 10-1 10-2 CHAPTER 10. FOURIER SERIES 8 > <0, sin mx sin nx dx = 0, > ⇡ : ⇡ Z (113) ⇡ and Z (114) m 6= n m = 0 or n = 0 m = n 6= 0 8 > <0, cos mx cos nx dx = 2⇡ > ⇡ : ⇡ ⇡ m 6= n m=n=0. m = n 6= 0 If p(x) is as in (111), then Theorem 10.2 shows Z ⇡ 2⇡↵0 = p(x) dx, ⇡ and for n > 0, ⇡↵n = Z ⇡ p(x) cos nx dx, ⇡ ⇡ n = Z ⇡ p(x) cos nx dx. ⇡ Combining these, it follows that if Z Z 1 ⇡ 1 an = p(x) cos nx dx and bn = ⇡ ⇡ ⇡ for n 2 !, then (115) p(x) = ⇡ p(x) sin nx dx ⇡ 1 a0 X + an cos nx + bn sin nx. 2 n=1 (Remember that all but a finite number of the an and bn are 0!) At this point, the logical question is whether this same method can be used to represent a more general 2⇡-periodic function. For any function f , integrable on [ ⇡, ⇡], the coefficients can be defined as above; i.e., for n 2 !, Z Z 1 ⇡ 1 ⇡ (116) an = f (x) cos nx dx and bn = f (x) sin nx dx. ⇡ ⇡ ⇡ ⇡ The numbers an and bn are called the Fourier coefficients of f . The problem is whether and in what sense an equation such as (115) might be true. This turns out to be a very deep and difficult question with no short answer.2 Because we don’t know whether equality in the sense of (115) is true, the usual practice is to write (117) f (x) ⇠ 1 a0 X + an cos nx + bn sin nx, 2 n=1 indicating that the series on the right is calculated from the function on the left using (116). The series is called the Fourier series for f . 2Many people, including me, would argue that the study of Fourier series has been the most important area of mathematical research over the past two centuries. Huge mathematical disciplines, including set theory, measure theory and harmonic analysis trace their lineage back to basic questions about Fourier series. Even after centuries of study, research in this area continues unabated. May 8, 2015 http://math.louisville.edu/⇠lee/ira 2. THE RIEMANN LEBESGUE LEMMA 10-3 Figure 1. This shows f (x) = |x|, s1 (x) and s3 (x), where sn (x) is the nth partial sum of the Fourier series for f . Example 10.1. Let f (x) = |x|. Since f is an even functions and sin nx is odd, Z 1 ⇡ bn = |x| cos nx dx = 0 ⇡ ⇡ for all n 2 N. On the other hand, 1 an = ⇡ for n 2 Z. Therefore, Z ⇡ 4 cos x 2 ⇡ (See Figure 10.1.) |x| ⇠ 8 <⇡, |x| cos nx dx = 2(cos n⇡ : ⇡ n2 ⇡ n=0 ⇡ 4 cos 3x 9⇡ 4 cos 5x 25⇡ 1) 4 cos 7x 49⇡ , n2N 4 cos 9x + ··· 81⇡ There are at least two fundamental questions arising from (117): Does the Fourier series of f converge to f ? Can f be recovered from its Fourier series, even if the Fourier series does not converge to f ? These are often called the convergence and representation questions, respectively. The next few sections will give some partial answers. 2. The Riemann Lebesgue Lemma We learned early in our study of series that the first and simplest convergence test is to check whether the terms go to zero. For Fourier series, this is always the case. Rb Theorem 10.3 (Riemann-Lebesgue Lemma). If f is a function such that a f exists, then Z b Z b lim f (t) cos ↵t dt = 0 and lim f (t) sin ↵t dt = 0. ↵!1 May 8, 2015 a ↵!1 a http://math.louisville.edu/⇠lee/ira 10-4 CHAPTER 10. FOURIER SERIES Proof. Since the two limits have similar proofs, only the first will be proved. Let " > 0 and P be a generic partition of [a, b] satisfying Z b " 0< f D (f, P ) < . 2 a For mi = glb {f (x) : xi 1 < x < xi }, define a function g on [a, b] by g(x) = mi Rb when xi 1 x < xi and g(b) = mn . Note that a g = D (f, P ), so Z b " (118) 0< (f g) < . 2 a Choose n (119) ↵> Since f 4X |mi |. " i=1 g, Z b f (t) cos ↵t dt = a Z Z Z Z b (f (t) g(t)) cos ↵t dt + a Z b (f (t) g(t)) cos ↵t dt + a b g(t) cos ↵t dt a Z b g(t) cos ↵t dt a n b (f g) + (f g) + a 1X mi (sin(↵xi ) ↵ i=1 sin(↵xi 1 )) n b a 2X |mi | ↵ i=1 Use (118) and (119). " " + 2 2 =" < ⇤ Corollary 10.4. If f is integrable on [ ⇡, ⇡] with an and bn the Fourier coefficients of f , then an ! 0 and bn ! 0. 3. The Dirichlet Kernel Suppose f is integrable on [ ⇡, ⇡] and 2⇡-periodic on R, so the Fourier series of f exists. The partial sums of the Fourier series are written as sn (f, x), or more simply sn (x) when there is only one function in sight. To be more precise, n sn (f, x) = a0 X + (ak cos kx + bk sin kx) . 2 k=1 Notice sn is a trigonometric polynomial of degree at most n. We begin with the following calculation. May 8, 2015 http://math.louisville.edu/⇠lee/ira 3. THE DIRICHLET KERNEL 10-5 15 12 9 6 3 -p - p p 2 2 p -3 Figure 2. The Dirichlet kernel Dn (s) for n = 1, 4, 7. n a0 X + (ak cos kx + bk sin kx) 2 k=1 Z ⇡ Z n X 1 1 ⇡ = f (t) dt + (f (t) cos kt cos kx + f (t) sin kt sin kx) dt 2⇡ ⇡ ⇡ ⇡ k=1 ! Z ⇡ n X 1 = f (t) 1 + 2(cos kt cos kx + sin kt sin kx) dt 2⇡ ⇡ k=1 ! Z ⇡ n X 1 = f (t) 1 + 2 cos k(x t) dt 2⇡ ⇡ sn (x) = k=1 Substitute s = x Z 1 (120) = 2⇡ t and use the assumption that f is 2⇡-periodic. ! n ⇡ X f (x s) 1 + 2 cos ks ds ⇡ k=1 The sequence of trigonometric polynomials from within the integral, n X (121) Dn (s) = 1 + 2 cos ks, k=1 is called the Dirichlet kernel. Its properties will prove useful for determining the pointwise convergence of Fourier series. Theorem 10.5. The Dirichlet kernel has the following properties. (a) Dn (s) is an even 2⇡-periodic function for each n 2 N. (b) Dn (0) = 2n + 1 for each n 2 N. (c) |DnZ (s)| 2n + 1 for each n 2 N and all s. ⇡ 1 (d) Dn (s) ds = 1 for each n 2 N. 2⇡ ⇡ May 8, 2015 http://math.louisville.edu/⇠lee/ira 10-6 CHAPTER 10. FOURIER SERIES (e) Dn (s) = sin(n + 1/2)s for each n 2 N and s/2 not an integer multiple sin s/2 of ⇡. Proof. Properties (a)–(d) follow from the definition of the kernel. The proof of property (e) uses some trigonometric manipulation. Suppose n 2 N and s 6= m⇡ for any m 2 Z. Dn (s) = 1 + 2 n X cos ks k=1 Use the facts that the cosine is even and the sine is odd. = n X cos ks + k= n = = 1 sin 2s 1 sin 2s n cos 2s X sin ks sin 2s k= n n ⇣ X k= n n X k= n sin ⌘ s s cos ks + cos sin ks 2 2 1 sin(k + )s 2 This is a telescoping sum. = sin(n + 12 )s sin 2s ⇤ According to (120), sn (f, x) = 1 2⇡ Z ⇡ f (x t)Dn (t) dt. ⇡ Figure 3. This graph shows D50 (t) and the envelope y = ±1/ sin(t/2). As n gets larger, the Dn (t) fills the envelope more completely. May 8, 2015 http://math.louisville.edu/⇠lee/ira 4. DINI’S TEST FOR POINTWISE CONVERGENCE 10-7 This is similar to a situation we’ve seen before within the proof of the Weierstrass approximation theorem, Theorem 9.13. The integral given above is a convolution integral similar to that used in the proof of Theorem 9.13, although the Dirichlet kernel isn’t a convolution kernel in the sense of Lemma 9.14 because it doesn’t satisfy conditions (a) and (c) of that lemma. (See Figure 10.3.) 4. Dini’s Test for Pointwise Convergence Theorem 10.6 (Dini’s Test). Let f : R ! R be a 2⇡-periodic function integrable on [ ⇡, ⇡] with Fourier series given by (117). If there is a > 0 and s 2 R such that Z f (x + t) + f (x t) 2s dt < 1, t 0 then 1 a0 X + (ak cos kx + bk cos kx) = s. 2 k=1 Proof. Since Dn is even, Z ⇡ 1 sn (x) = f (x t)Dn (t) dt 2⇡ ⇡ Z 0 Z ⇡ 1 1 = f (x t)Dn (t) dt + f (x 2⇡ 2⇡ 0 ⇡ Z ⇡ 1 = (f (x + t) + f (x t)) Dn (t) dt. 2⇡ 0 t)Dn (t) dt By Theorem 10.5(d) and (e), Z ⇡ 1 sn (x) s = (f (x + t) + f (x t) 2s) Dn (t) dt 2⇡ 0 Z ⇡ 1 f (x + t) + f (x t) 2s t 1 = · · sin(n + )t dt. 2⇡ 0 t 2 sin 2t Since t/ sin 2t is bounded on (0, ⇡), Theorem 10.3 shows sn (x) Corollary 8.11 to finish the proof. s ! 0. Now use ⇤ Example 10.2. Suppose f (x) = x for ⇡ < x ⇡ and is 2⇡-periodic on R. Since f is odd, an = 0 for all n. Integration by parts gives bn = ( 1)n+1 2/n for n 2 N. Therefore, 1 X 2 f⇠ ( 1)n+1 sin nx. n n=1 For x 2 ( ⇡, ⇡), let 0 < < min{⇡ x, ⇡ + x}. (This is just the distance from x to closest endpoint of ( ⇡, ⇡).) Using Dini’s test, we see Z Z f (x + t) + f (x t) 2x x + t + x t 2x dt = dt = 0 < 1, t t 0 0 so 1 X 2 (122) ( 1)n+1 sin nx = x for ⇡ < x < ⇡. n n=1 May 8, 2015 http://math.louisville.edu/⇠lee/ira 10-8 CHAPTER 10. FOURIER SERIES p p 2 -p - p p 2 2 - p 2p p 2 -p Figure 4. This plot shows the function of Example 10.2 and s8 (x) for that function. In particular, when x = ⇡/2, (122) gives another way to derive (104). When x = ⇡, the series converges to 0, which is the middle of the “jump” for f . This behavior of converging to the middle of a jump discontinuity is typical. To see this, denote the one-sided limits of f at x by f (x ) = lim f (t) and f (x+) = lim f (t), t"x t#x and suppose f has a jump discontinuity at x with f (x ) + f (x+) s= . 2 Guided by Dini’s test, consider Z f (x + t) + f (x t) 2s dt t 0 Z f (x + t) + f (x t) f (x ) f (x+) = dt t 0 Z Z f (x + t) f (x+) f (x t) f (x ) dt + dt t t 0 0 If both of the integrals on the right are finite, then the integral on the left is also finite. This amounts to a proof of the following corollary. Corollary 10.7. Suppose f : R ! R is 2⇡-periodic and integrable on [ ⇡, ⇡]. If both one-sided limits exist at x and there is a > 0 such that both Z Z f (x + t) f (x+) f (x t) f (x ) dt < 1 and dt < 1, t t 0 0 then the Fourier series of f converges to f (x ) + f (x+) . 2 May 8, 2015 http://math.louisville.edu/⇠lee/ira 5. GIBBS PHENOMENON 10-9 Figure 5. This is a plot of s19 (f, x), where f is defined by (123). The Dini test given above provides a powerful condition sufficient to ensure the pointwise convergence of a Fourier series. There is a plethora of ever more abstruse conditions that can be proved in a similar fashion to show pointwise convergence. The problem is complicated by the fact that there are continuous functions with Fourier series divergent at a point and integrable functions with Fourier series diverging everywhere [13]. In Section 6, a continuous function whose Fourier series diverges is constructed. 5. Gibbs Phenomenon For x 2 [ ⇡, ⇡) define (123) f (x) = ( |x| x , 0, 0<x<⇡ x = 0, ⇡ and extend f 2⇡-periodically to all of R. This function is often called a square wave. A straightforward calculation gives f⇠ 1 4 X sin(2k 1)x . ⇡ 2k 1 k=1 Corollary 10.7 shows sn (x) ! f (x) everywhere. This convergence cannot be uniform because all the partial sums are continuous and f is discontinuous at every integer multiple of ⇡. A plot of s19 (x) is shown in Figure 5. Notice the higher peaks in the oscillation of sn (x) just before and after the jump discontinuities of f . This behavior is not unique to f , as it can also be seen in Figure 4. If a function is discontinuous at some point, the partial sums of its Fourier series will always have such higher peaks near that point. This behavior is called Gibbs phenomenon.3 Instead of doing a general analysis of Gibbs phenomenon, we’ll only analyze the simple case shown in the square wave f . It’s basically a calculus exercise. 3It is named after the American mathematical physicist, J. W. Gibbs, who pointed it out in 1899. He was not the first to notice the phenomenon, as the British mathematician Henry Wilbraham had published a little-noticed paper on it it 1848. Gibbs’ interest in the phenomenon was sparked by investigations of the American experimental physicist A. A. Michelson who wrote a letter to Nature disputing the possibility that a Fourier series could converge to a discontinuous function. The ensuing imbroglio is recounted in a marvelous book by Paul J. Nahin [16]. May 8, 2015 http://math.louisville.edu/⇠lee/ira 10-10 CHAPTER 10. FOURIER SERIES Figure 6. This is a plot of s9 (f, x), where f is defined by (123). To locate the peaks in the graph, di↵erentiate the partial sums. n n d 4 X sin(2k 1)x 4X s02n 1 (x) = = cos(2k 1)x dx ⇡ 2k 1 ⇡ k=1 k=1 It is left as Exercise 10.10 to show this has a closed form. 2 sin 2nx s02n 1 (x) = ⇡ sin x Looking at the numerator, we see s02n 1 (x) has critical numbers at x = k⇡/2n for k 2 Z. In the interval (0, ⇡), s2n 1 (k⇡/2n) is a relative maximum for odd k and a relative minimum for even k. (See Figure 6.) The value s2n 1 (⇡/2n) is the height of the left-most peak. What is the behavior of these maxima? From Figure 7 it appears they have an asymptotic limit near 1.179. To prove this, consider the following calculation. n ⇡ ⇣ ⇡ ⌘ 4 X sin (2k 1) 2n s2n 1 f, = 2n ⇡ 2k 1 k=1 ⇣ ⌘ (2k 1)⇡ n 2n 2 X sin ⇡ = (2k 1)⇡ ⇡ n k=1 2n Figure 7. This is a plot of sn (f, ⇡/2n) for n = 1, 2, · · · , 100. The dots come in pairs because s2n May 8, 2015 1 (f, ⇡/2n) = s2n (f, ⇡/2n). http://math.louisville.edu/⇠lee/ira 6. A DIVERGENT FOURIER SERIES 10-11 The last sum is a midpoint Riemann sum for the function sinx x on the interval [0, ⇡] using a regular partition of n subintervals. Example 9.16 shows Z 2 ⇡ sin x dx ⇡ 1.17898. ⇡ 0 x Since f (0+) f (0 ) = 2, this is an overshoot of a bit less than 9%. There is a similar undershoot on the other side. It turns out this is typical behavior at points of discontinuity [21]. 6. A Divergent Fourier Series This is an advanced section that can be omitted. Eighteenth century mathematicians, including Fourier, Cauchy, Euler, Weierstrass, Lagrange and Riemann, had computed the Fourier series for many functions. They believed from these examples that the Fourier series of a continuous function must converge to that function. Fourier went far beyond this claim in his important 1822 book Th´eorie Analytique de la Chaleur. Cauchy even published a flawed proof of this “fact.” Finally, in 1873, Paul du Bois-Reymond, settled the question by giving the construction of a continuous function whose Fourier series diverges at a point [17]. It was finally shown in 1966 by Lennart Carleson [7] that the Fourier series of a continuous function converges to that function everywhere with the exception of a set of measure zero. The problems around the convergence of Fourier series motivated a huge amount of research that is still going on today. In this section, we look at the tip of that iceberg by presenting a continuous function F (t) whose Fourier series fails to converge when t = 0. 6.1. The Conjugate Dirichlet Kernel. Lemma 10.8. If m, n 2 ! and 0 < |t| < ⇡ for k 2 Z, then n X sin kt = cos(m 1 2 )t cos(m + 12 )t 2 sin 2t k=m . Proof. n X k=m n 1 X 1 sin sin kt t 2 sin 2 k=m ✓ ✓ ◆ n 1 X 1 = cos k t 2 sin 2t kt = k=m = cos m 1 2 t cos n + 2 sin 2t 1 2 ✓ ◆ ◆ 1 cos k + t 2 t ⇤ May 8, 2015 http://math.louisville.edu/⇠lee/ira 10-12 CHAPTER 10. FOURIER SERIES Definition 10.9. The conjugate Dirichlet kernel 4 is n X ˜ n (t) = D (124) n 2 N. sin kt, k=1 We’ll not have much use for the conjugate Dirichlet kernel, except as a convenient way to refer to sums of the form (124). Lemma 10.8 immediately gives the following bound. Corollary 10.10. If 0 < |t| < ⇡, then 1 . sin 2t ˜ n (t) D 6.2. A Sawtooth Wave. If the function f (x) = (⇡ x)/2 on [0, 2⇡) is extended 2⇡-periodically to R, then the graph of the resulting function is often referred to as a “sawtooth wave”. It has a particularly nice Fourier series: ⇡ x 2 ⇠ 1 X sin kx k k=1 . According to Corollary 10.7 1 X sin kx k=1 k = ( x = 2n⇡, n 2 Z . otherwise 0, f (x), We’re interested in various partial sums of this series. Lemma 10.11. If m, n 2 ! with m n and 0 < |t| < 2⇡, then n X sin kt 1 t . k m sin 2 k=m Proof. n n ⇣ X X sin kt ˜ k (t) = D k k=m ˜k D 1 (t) k=m ⌘1 k Use summation by parts. = n X k=m n X k=m ˜ k (t) D ˜ k (t) D ✓ ✓ 1 k 1 k+1 1 k 1 k+1 ◆ ◆ + + ˜ n (t) D n+1 ˜n D 1 (t) m ˜ n (t) ˜ n 1 (t) D D + n+1 m 4In this case, the word “conjugate” does not refer to the complex conjugate, but to the ˜ n (t) = 1 + 2 Pn ekit . harmonic conjugate. They are related by Dn (t) + iD k=1 May 8, 2015 http://math.louisville.edu/⇠lee/ira 6. A DIVERGENT FOURIER SERIES 10-13 Apply Corollary 10.10. n ✓ X 1 k 1 2 sin 2t k=m ✓ 1 1 = t 2 sin 2 m 1 = . m sin 2t 1 k+1 ◆ 1 1 + n+1 m ◆ 1 1 1 + + n+1 n+1 m ! ⇤ Proposition 10.12. If n 2 N and 0 < |t| < ⇡, then n X sin kt k=1 k 1 + ⇡. Proof. n X sin kt k=1 k = X sin kt X sin kt + k k 1 1 1k t t <kn X sin kt + k 1 1k t X kt + k 1 1k t X 1k 1t t+ X sin kt k 1 t <kn 1 1 t t sin 2 1 1 2 t t ⇡2 1+⇡ ⇤ 6.3. A Continuous Function with a Divergent Fourier Series. The goal in this subsection is to construct a continuous function F such that lim sup sn (F, 0) = 1. This implies sn (F, 0) does not converge. Example 10.3. The first step in the construction is to define a sequence of trigonometric polynomials that form the building blocks for F . fn (t) = 1 cos t cos 2t cos(n 1)t + + + ··· n n 1 n 2 1 cos(n + 1)t 1 cos(n + 2)t 2 ··· cos 2nt n Note that, (125) May 8, 2015 ( sm (fn , 0) = fn (0) = 0, sm (fn , 0) > 0, m 2n . 1 m < 2n http://math.louisville.edu/⇠lee/ira 10-14 CHAPTER 10. FOURIER SERIES Rearrange the sum in the definition of fn to see fn (t) = n X cos(n k)t k k=1 = 2 sin nt n X k=1 cos(n + k)t sin kt . k This closed form for fn combined with Proposition 10.12 implies the sequence of functions fn is uniformly bounded: (126) |fn (t)| = 2 sin nt n X sin kt k k=1 2 + 2⇡. At last, the main function can be defined. F (t) = 1 X f2n3 (t) n2 n=1 The Weierstrass M-Test along with (126) implies F is uniformly convergent and therefore continuous on R. Using (125), consider Z ⇡ 1 3 s2m (F, 0) = F (t)D2m3 (t) dt 2⇡ ⇡ Z ⇡ X 1 f2n3 (t) 1 = D2m3 (t) dt 2⇡ n2 ⇡ n=1 The uniform convergence allows the sum and integration to be reordered. Z 1 X 1 1 = n2 2⇡ n=1 = ⇡ ⇡ f2n3 (t)D2m3 (t) dt 1 X 1 s m3 f2n3 , 0 n2 2 n=1 Use (125). 1 X 1 = s 3 f 2n3 , 0 2 2m n n=m > 1 s m3 f2m3 , 0 m2 2 m3 2 1 X1 = 2 m k k=1 3 1 > 2 ln 2m m = m ln 2 This implies, lim sup sn (F, 0) May 8, 2015 limm!1 m ln 2 = 1, so sn (F, 0) does not converge. http://math.louisville.edu/⇠lee/ira ´ KERNEL 7. THE FEJER 10-15 7. The Fej´ er Kernel Since pointwise convergence of the partial sums seems complicated, why not change the rules of the game? Instead of looking at the sequence of partial sums, consider a rolling average instead: n n (f, x) = 1 X sn (f, x). n+1 k=0 The trigonometric polynomials n (f, x) are called the Ces`aro means of the partial sums. If limn!1 n (f, x) exists, then the Fourier series for f is said to be (C, 1) summable at x. The idea is that this averaging will “smooth out” the partial sums, making them more nicely behaved. It is not hard to show that if sn (f, x) converges at some x, then n (f, x) will converge to the same thing. But there are sequences for which n (f, x) converges and sn (f, x) does not. (See Exercises 3.19 and 4.26.) As with sn (x), we’ll simply write n (x) instead of n (f, x), when it is clear which function is being considered. We start with a calculation. n n (x) = = (*) = = = 1 X sk (x) n+1 k=0 Z ⇡ n 1 X 1 f (x t)Dk (t) dt n+1 2⇡ ⇡ k=0 Z ⇡ n 1 1 X f (x t) Dk (t) dt 2⇡ n+1 ⇡ k=0 Z ⇡ n 1 1 X sin(k + 1/2)t f (x t) dt 2⇡ n+1 sin t/2 ⇡ k=0 Z ⇡ n X 1 1 f (x t) sin t/2 sin(k + 1/2)t dt 2⇡ (n + 1) sin2 t/2 k=0 ⇡ Use the identity 2 sin A sin B = cos(A 1 = 2⇡ Z The sum telescopes. Z 1 = 2⇡ ⇡ f (x t) f (x t) ⇡ ⇡ ⇡ Use the identity 2 sin2 A = 1 (**) May 8, 2015 1 = 2⇡ Z ⇡ f (x ⇡ B) cos(A + B). n X 1/2 (cos kt (n + 1) sin2 t/2 k=0 1/2 (1 (n + 1) sin2 t/2 cos(k + 1)t) dt cos(n + 1)t) dt cos 2A. 1 t) (n + 1) ✓ sin n+1 2 t sin 2t ◆2 dt http://math.louisville.edu/⇠lee/ira 10-16 CHAPTER 10. FOURIER SERIES 10 8 6 4 2 -p - p 2 p 2 p Figure 8. A plot of K5 (t), K8 (t) and K10 (t). The Fej´er kernel is the sequence of functions highlighted above; i.e., ✓ ◆2 sin n+1 1 2 t (127) Kn (t) = , n 2 N. (n + 1) sin 2t Comparing the lines labeled (*) and (**) in the previous calculation, we see another form for the Fej´er kernel is n (128) Kn (t) = 1 X Dk (t). n+1 k=0 Once again, we’re confronted with a convolution integral containing a kernel: Z ⇡ 1 (x) = f (x t)Kn (t) dt. n 2⇡ ⇡ Theorem 10.13. The Fej´er kernel has the following properties.5 (a) (b) (c) (d) (e) (f) Kn (t) is an even 2⇡-periodic function for each n 2 N. Kn (0) = n + 1 for each n 2 !. Kn (t) R ⇡ 0 for each n 2 N. 1 Kn (t) dt = 1 for each n 2 !. 2⇡ ⇡ If 0 < < ⇡, then Kn ◆ 0 on [ ⇡, ] [ [ , ⇡]. R R⇡ If 0 < < ⇡, then ⇡ Kn (t) dt ! 0 and Kn (t) dt ! 0. Proof. Theorem 10.5 and (128) imply (a), (b) and (d). Equation (127) implies (c). 5Compare this theorem with Lemma 9.14. May 8, 2015 http://math.louisville.edu/⇠lee/ira ´ KERNEL 7. THE FEJER 10-17 Let be as in (e). In light of (a), it suffices to prove (e) for the interval [ , ⇡]. Noting that sin t/2 is decreasing on [ , ⇡], it follows that for t ⇡, ✓ ◆2 sin n+1 1 2 t Kn (t) = (n + 1) sin 2t ✓ ◆2 1 1 (n + 1) sin 2t 1 1 !0 (n + 1) sin2 2 It follows that Kn ◆ 0 on [ , ⇡] and (e) has been proved. Theorem 9.16 and (e) imply (f). ⇤ Theorem 10.14 (Fej´er). If f : R ! R is 2⇡-periodic, integrable on [ ⇡, ⇡] and continuous at x, then n (x) ! f (x). R⇡ R⇡ R Proof. Since f is 2⇡-periodic and ⇡ f (t) dt exists, so does ⇡ (f (x t) f (x)) dt. Theorem 8.3 gives an M > 0 so |f (x t) f (x)| < M for all t. Let " > 0 and choose > 0 such that |f (x) f (y)| < "/3 whenever |x y| < . By Theorem 10.13(f), there is an N 2 N so that whenever n N , Z Z ⇡ 1 " 1 " Kn (t) dt < and Kn (t) dt < . 2⇡ 3M 2⇡ 3M ⇡ We start calculating. | n (x) = Z ⇡ 1 f (x t)Kn (t) dt f (x)Kn (t) dt 2⇡ ⇡ ⇡ Z ⇡ 1 = (f (x t) f (x))Kn (t) dt 2⇡ ⇡ Z (f (x t) f (x))Kn (t) dt + (f (x t) f (x))Kn (t) dt f (x)| = 1 2⇡ Z ⇡ 1 2⇡ Z ⇡ + 1 2⇡ Z (f (x t) M 2⇡ This shows Z Kn (t) dt + ⇡ n (x) ! f (x). ⇡ (f (x t) f (x))Kn (t) dt f (x))Kn (t) dt + ⇡ + < Z 1 2⇡ 1 2⇡ Z Z ⇡ (f (x |f (x t) t) 1 2⇡ Z (f (x t) f (x))Kn (t) dt f (x))Kn (t) dt f (x)|Kn (t) dt + " " 1 <M + 3M 3 2⇡ Z M 2⇡ Z ⇡ Kn (t) dt Kn (t) dt + M " <" 3M ⇤ Theorem 10.14 gives a partial solution to the representation problem. Corollary 10.15. Suppose f and g are continuous and 2⇡-periodic on R. If f and g have the same Fourier coefficients, then they are equal. May 8, 2015 http://math.louisville.edu/⇠lee/ira 10-18 CHAPTER 10. FOURIER SERIES p p p p 2 2 -p - p p 2 2 - 2 p -p p - p p 2 p 2 - 2 -p p 2p p 2 -p Figure 9. These plots illustrate the functions of Example 10.4. On the left are shown f (x), s8 (x) and 8 (x). On the right are shown f (x), 3 (x), 5 (x) and 20 (x). Compare this with Figure 4. Proof. By assumption, 0= n (f, t) = n (f, t) n (g, t) n (g, t) for all n and Theorem 10.14 implies !f g. ⇤ In the case of continuous functions, the convergence is uniform, rather than pointwise. Theorem 10.16 (Fej´er). If f : R ! R is 2⇡-periodic and continuous, then (x) ◆ f (x). n Proof. By Exercise 6.32, f is uniformly continuous. This can be used to show the calculation within the proof of Theorem 10.14 does not depend on x. The details are left as Exercise 10.12. ⇤ A perspicacious reader will have noticed the similarity between Theorem 10.16 and the Weierstrass approximation theorem, Theorem 9.13. In fact, the Weierstrass approximation theorem can be proved from Theorem 10.16 using power series and Theorem 9.24. (See Exercise 10.13.) Example 10.4. As in Example 10.2, let f (x) = x for ⇡ < x ⇡ and extend f to be periodic on R with period 2⇡. Figure 9 shows the di↵erence between the Fej´er and classical methods of summation. Notice that the Fej´er sums remain much more smoothly affixed to the function and do not show Gibbs phenomenon. 8. Exercises 10.1. Prove Theorem 10.2. 10.2. Suppose f : R ! R is 2⇡-periodic and integrable on [ ⇡, ⇡]. If b Rb R⇡ then a f = ⇡ f . a = 2⇡, 10.3. Let f : [ ⇡, ⇡) be a function with Fourier coefficients an and bn . If f is odd, then an = 0 for all n 2 !. If f is even, then bn = 0 for all n 2 N. 10.4. If f (x) = sgn(x) on [ ⇡, ⇡), then find the Fourier series for f . May 8, 2015 http://math.louisville.edu/⇠lee/ira 8. EXERCISES 10.5. Is 1 X 10-19 sin nx the Fourier series of some function? n=1 10.6. Let f (x) = x2 when ⇡ x < ⇡ and extend f to be 2⇡-periodic on all of R. Use Dini’s Test to show sn (f, x) ! f (x) everywhere. 10.7. Use Exercise 10.6 to prove 1 X ⇡2 1 = . 2 6 n n=1 10.8. Suppose f is integrable on [ ⇡, ⇡]. If f is even, then the Fourier series for f has no sine terms. If f is odd, then the Fourier series for f has no cosine terms. 10.9. If f : R ! R is periodic with period p and continuous on [0, p], then f is uniformly continuous. 10.10. Prove n X cos(2k k=1 for t 6= k⇡, k 2 Z and n 2 N. (Hint: 2 1)t = Pn k=1 sin 2nt 2 sin t cos(2k 1)t = D2n (t) Dn (2t).) 10.11. The function g(t) = t/ sin(t/2) is undefined whenever t = 2n⇡ for some n 2 Z. Show that it can be redefined on the set {2n⇡ : n 2 Z} to be periodic and uniformly continuous on R. 10.12. Prove Theorem 10.16. 10.13. Prove the Weierstrass approximation theorem using Fourier series and Taylor’s theorem. 10.14. If f (x) = x for ⇡ x < ⇡, can the Fourier series of f be integrated term-by-term on ( ⇡, ⇡) to obtain the correct answer? 10.15. If f (x) = x for ⇡ x < ⇡, can the Fourier series of f be di↵erentiated term-by-term on ( ⇡, ⇡) to obtain the correct answer? May 8, 2015 http://math.louisville.edu/⇠lee/ira
© Copyright 2025