G: Uniform Convergence of Fourier Series
Total Page:16
File Type:pdf, Size:1020Kb
G: Uniform Convergence of Fourier Series From previous work on the prototypical problem (and other problems) 8 < ut = Duxx 0 < x < l ; t > 0 u(0; t) = 0 = u(l; t) t > 0 (1) : u(x; 0) = f(x) 0 < x < l we developed a (formal) series solution 1 1 X X 2 2 2 nπx u(x; t) = u (x; t) = b e−n π Dt=l sin( ) ; (2) n n l n=1 n=1 2 R l nπy with bn = l 0 f(y) sin( l )dy. These are the Fourier sine coefficients for the initial data function f(x) on [0; l]. We have no real way to check that the series representation (2) is a solution to (1) because we do not know we can interchange differentiation and infinite summation. We have only assumed that up to now. In actuality, (2) makes sense as a solution to (1) if the series is uniformly convergent on [0; l] (and its derivatives also converges uniformly1). So we first discuss conditions for an infinite series to be differentiated (and integrated) term-by-term. This can be done if the infinite series and its derivatives converge uniformly. We list some results here that will establish this, but you should consult Appendix B on calculus facts, and review definitions of convergence of a series of numbers, absolute convergence of such a series, and uniform convergence of sequences and series of functions. Proofs of the following results can be found in any reasonable real analysis or advanced calculus textbook. 0.1 Differentiation and integration of infinite series Let I = [a; b] be any real interval. Theorem 1 (term-by-term differentiation): If, on the interval I, P1 a) f(x) = n=1 fn(x) converges uniformly; P1 0 b) n=1 fn(x) converges uniformly; 1We can get away with less restrictive conditions, but the analysis is beyond the level of these Notes. 1 0 c) fn(x) is continuous. P1 0 Then the series n=1 fn(x) can be differentiated term-by-term; that is, f (x) = P1 0 n=1 fn(x). Theorem 2 (term-by-term integration): If, on the interval I, P1 a) f(x) = n=1 fn(x) converges uniformly; 0 b) each fn(x) is integrable. P1 R b Then the series n=1 fn(x) can be integrated term-by-term; that is, a f(x)dx = P1 R b n=1 a fn(x)dx. Remark: Note the weaker conditions needed for integration versus differ- entiation. Theorem 3 (Weirstrass M-test for uniform convergence of an infi- nite series of functions): Suppose ffngn≥1 is a sequence of functions defined on I, and suppose jfn(x)j ≤ P1 Mn for all x 2 I, n ≥ 1. Then n=1 fn(x) converges uniformly in I if the P1 series of numbers n=1 Mn converges absolutely. Theorem 4 (Ratio test): P1 A series of numbers n=1 Mn converges absolutely if, for some N ≥ 1, and some constant 0 < ρ < 1, jMn+1=Mnj ≤ ρ < 1, for all n ≥ N. Now fix t0 > 0 and consider series (2) for all t ≥ t0, 0 ≤ x ≤ l, and also R l assume f is absolutely integrable (i.e. 0 jf(x)jdx < 1); f being piecewise continuous on 0 ≤ x ≤ l is more than sufficient. Then we have 2 2 2 nπx 2 2 2 ju (x; t)j = jb je−n π Dt=l j sin( )j ≤ jb je−n π Dt0=l ; n n l n 2 2 2 R l R l −n π Dt0=l where jbnj ≤ (2=l) 0 jf(y)jdy. Define Mn := ((2=l) 0 jf(y)jdy)e . Then, for all t ≥ t0, jun(x; t)j ≤ Mn. Also, 2 2 2 −(n+1) π Dt0=l Mn+1 e 2 2 2 2 2 2 2 2 = = e[n −(n+1) ]π Dt0=l = e−(2n+1)π Dt0=l ≤ e−π Dt0=l < 1 : −n2π2Dt =l2 Mn e 0 2 2 2 −π Dt0=l P1 Thus, with ρ := e , n=1 Mn converges absolutely, by Theorem 4, P1 and by Theorem 3, n=1 un converges uniformly on [0; l], t ≥ t0 > 0. A similar argument holds for the convergence of derivatives ut; uxx. So, for t ≥ t0, the series (2) for u can be differentiated term-by-term. Since each term, un(x; t), satisfies the pde and b.c.s, so does u. After our discussion of the properties of the Fourier series, and the uniform convergence result on the Fourier series, the convergence of u holds all the way down to t = 0 (given the appropriate conditions on u(x; 0) = f(x)). For the latter argument we turn to the uniform convergence of the trig series. 0.2 Uniform convergence of classical Fourier series Let2 f be piecewise smooth on (−1; 1), continuous on [−1; 1], with f(−1) = f(1). Thus, when f is considered extended to the whole real line, it is continuous everywhere, and is a 2-periodic function on R. Then its Fourier series converges everywhere (pointwise) to f. For purposes below we will also 0 R 1 0 2 assume f is square-integrable; that is, −1(f (x)) dx < 1. Let 1 a0 X S (x) = + fa cos(nπx) + b sin(nπx)g ; N 2 n n n=1 where fan; bng are the Fourier coefficients for f : 2 3 2 3 an Z 1 cos(nπx) 4 5 = f(x) 4 5 dx : −1 bn sin(nπx) The purpose here is to show the sequence of partial sums fSN gN≥1 is a Cauchy sequence, which is sufficient for the Fourier series for f to converge uniformly to f. Recall that fSN gN≥1 is a Cauchy sequence if, for any " > 0, ∗ ∗ there is an N such that for any N; M > N , jSN (x) − SM (x)j < " for all x 2 (−1; 1). Assume without loss of generality that N > M; from the definition of SN , N X SN (x) − SM (x) = fan cos(nπx) + bn sin(nπx)g : (3) n=M+1 2For convenience only we have scaled the interval to have semi-length unity by redefin- ingx ~ = x=l, then dropping the tilde notation. 3 The essential \trick" is to write the right side of (3) as a dot product of two vectors: SN (x) − SM (x) = (4) (π(M + 1)aM+1; π(M + 1)bM+1; : : : ; πNaM+1; πNaM+1)· cos((M + 1)πx) sin((M + 1)πx) cos(Nπx) sin(Nπx) ( ; ;:::; ; ) : π(M + 1) π(M + 1) πN πN Now we need three pieces of information, namely, P1 1 2 A. n=1 n2 = π =6 : N N B. (Schwarz's inequality): given sequences fcngn=1; fdngn=1, then N N N X X 2 1=2 X 2 1=2 j cndnj ≤ ( cn) ( dn) : n=1 n=1 n=1 0 C. (Bessel's inequality): if fan; bng are the Fourier coefficients of f, where f is square-integrable on (−1; 1), then 1 Z 1 2 X 2 2 2 0 2 π n (an + bn) ≤ (f (x)) dx : n=1 −1 Remark: In the last section of this appendix are proofs of these three state- ments. Now apply Schwarz's inequality B to (4): N N X X 1 jS (x) − S (x)j ≤ ( (πn)2(a2 + b2 ))1=2( [cos2(nπx) + sin2(nπx)])1=2 N M n n (nπ)2 n=M+1 n=M+1 N N X X 1 = ( n2(a2 + b2 ))1=2( )1=2 n n n2 n=M+1 n=M+1 N 1 Z 1 X 1 ≤ ( jf 0(x)j2dx)1=2( )1=2 π2 n2 −1 n=M+1 using Bessel's inequality C in the last step. Since, from fact A, N M N π2 X 1 π2 X 1 X 1 0 < − = − − ; 6 n2 6 n2 n2 n=1 n=1 n=M+1 4 then 1=2 M !1=2 1 Z 1 π2 X 1 jS (x) − S (x)j ≤ jf 0(x)j2dx − ; (5) N M π2 6 n2 −1 n=1 where the first factor on the right side of the inequality is just a constant, and the second factor can be made arbitrarily small if M is taken sufficient large. Hence, fSN g is a Cauchy sequence. Remark: We can let N ! 1 in (5) to obtain 1=2 M !1=2 1 Z 1 π2 X 1 jf(x) − S (x)j ≤ jf 0(x)j2dx − : M π2 6 n2 −1 n=1 Example: f(x) = x2 on [−1; 1]; computing the Fourier coefficients, we have 1 1 4 X (−1)n f(x) = + cos(nπx) : 3 π2 n2 n=1 (Note: the equality is fully justified in this case.) Now, Z 1 Z 1 1=2 0 2 8 1 0 2 ∼ (f (x)) dx = ; so 2 (f (x)) dx = 0:519798 : −1 3 π −1 Then for π2 PM 1 1=2 M (0:519798)( 6 − n=1 n2 ) 2 0.32667 10 0.16035 50 0.07314 100 0.05185 1000 0.01643 100000 0.00164 Remark: In actuality the theory gives very conservative estimates. You only need a small number of terms usually to get a high order accurate approxi- mation for series that converge uniformly. 5 Exercise: Compute estimates of jf(x) − SM (x)j for f(x) = 1 − jxj on [−1; 1], for M = 10; 100; 1000. 0.3 Proofs of statements A, B, and C 1.