09CnvSer.mws

Partial Differential Equations PowerTool

by Dr. Jim Herod

Section 2.3: Convergence of Fourier Series

Maple Packages for Section 2.3

>    restart:

     In the Section 2.1, we discussed three general types of convergence in C([0, 1]): normed, pointwise, and uniform. Here, we apply these ideas to the special case that we have a function f  and we have the Fourier series

                      S[n](x) = Sum(`<,>`(f,phi[p])*phi[p](x),p = 0 .. n) .

where phi[1], phi[2], phi[3]   ...   is the usual orthonormal sequence. These results are presented here as Theorems, without proofs. The results may be found in many texts in texts on Fourier Series and Boundary Value Problems.

There are several terms we will use in this module.

A function is sectionally continuous  on an interval [a, b] if it is continuous on that interval except for possibly a finite number of jumps and removable discontinuities.

A function is sectionally smooth  on an interval [a, b] if both f   and   f  ' are sectionally continuous on the interval [a, b].

Examples:  The function f(x) = signum(x) is sectionally continuous on [-1, 1], but the function g(x) = 1/x is not sectionally continuous on that interval. Both are continuous on that interval except at zero. The signum function has a jump at x = 0, but 1/x has a more serious discontinuity.

>    plot([signum(x),1/x+1/10],x=-1..1,y=-3..3,discont=true,
   color=[BLACK, RED]);

[Maple Plot]

>   

Example: The function sqrt(abs(x))  is continuous, but not sectionally smooth on [-1, 1]. It is not sectionally smooth because the derivative goes to infinity  as x approaches zero.

>    Diff(sqrt(abs(x)),x)=diff(sqrt(abs(x)),x);
plot([rhs(%),sqrt(abs(x))],x=-1..1,y=-3..3,discont=true,color=[red,black]);

Diff(abs(x)^(1/2),x) = 1/2*1/abs(x)^(1/2)*abs(1,x)

[Maple Plot]

>   

Here is the first  result for convergence of series.

THEOREM:  If the function f  is sectionally smooth and periodic with period 2  c , then at each point x the Fourier series for f  converges and

        a[0]+sum(a[n]*cos(n*pi*x/c)+b[n]*sin(n*pi*x/c),n = 1 .. infinity)  = [ f(x+) + f(x-)] / 2.

Example:  The function signum(x) is sectionally smooth. Therefore the Fourier series for this function converges to 1 for 0 < x < 1, to  -1 for -1 < x < 0, and to 0 for x = -1, 0, or 1. The Fourier series for the function has period 2. Here is a plot for 5 terms of the series.

>    c:=1;
f:=x->signum(x);
a[0]:=int(f(x),x=-c..c)/int(1^2,x=-c..c);
for n from 1 to 5 do
   a[n]:=int(f(x)*cos((n*Pi*x)/c),x=-c..c)/
          int(cos(n*Pi*x/c)^2,x=-c..c);
   b[n]:=int(f(x)*sin((n*Pi*x)/c),x=-c..c)/
          int(sin(n*Pi*x/c)^2,x=-c..c);
od;
n:='n':

c := 1

f := signum

a[0] := 0

a[1] := 0

b[1] := 4/Pi

a[2] := 0

b[2] := 0

a[3] := 0

b[3] := 4/3/Pi

a[4] := 0

b[4] := 0

a[5] := 0

b[5] := 4/5/Pi

>    s:=x->a[0]+sum(a[n]*cos((n*Pi*x)/c)+b[n]*sin((n*Pi*x)/c),n=1..5);

s := proc (x) options operator, arrow; a[0]+sum(a[n]*cos(n*Pi*x/c)+b[n]*sin(n*Pi*x/c),n = 1 .. 5) end proc

>    plot([[x,f(x),x=-c..c],[x,s(x),x=-2*c..2*c]]);

[Maple Plot]

>   

In this previous example, the function was not continuous. We will see that a discontinuity always leads to the overshoot that can be observed, no matter how many terms of the series are taken. Is it not clear that the function would be getting close to f ( x ) if more terms were used?

In this next Theorem, we see a more powerful result that will insure uniform convergence for all x .

THEOREM:  If the series sum(abs(a[n])+abs(b[n]),n = 1 .. infinity)  converges, then the Fourier series for f  converges uniformly in the interval [-c, c].

Example:  Take the series a[n] = 2*(-1+(-1)^n)/(Pi*n^2)   , with b[n] = 0 . Since the number series sum(1/(n^2),n = 1 .. infinity)  converges, the series of absolute values converges, and so a Fourier Series with these coefficients converges. We draw graphs for this series on [- Pi , Pi  ].

>    c:=Pi;
a[0]:=Pi/2;
for n from 1 to 5 do
    a[n]:=2/Pi*(-1+(-1)^n)/n^2;
    b[n]:=0;
od;
n:='n':

c := Pi

a[0] := 1/2*Pi

a[1] := -4/Pi

b[1] := 0

a[2] := 0

b[2] := 0

a[3] := -4/9/Pi

b[3] := 0

a[4] := 0

b[4] := 0

a[5] := -4/25/Pi

b[5] := 0

>    s:=x->a[0]+sum(a[n]*cos((n*Pi*x)/c)+b[n]*sin((n*Pi*x)/c),n=1..5);

s := proc (x) options operator, arrow; a[0]+sum(a[n]*cos(n*Pi*x/c)+b[n]*sin(n*Pi*x/c),n = 1 .. 5) end proc

>    plot(s(x),x=-2*c..2*c);

[Maple Plot]

>   

Example:  If { phi[1], phi[2], phi[3] , ... } is an orthogonal sequence and we construct the Fourier Series

                       sum(a[p]*phi[p],p)

then this series will converge in norm if and only if

                       sum(abs(a[p])^2,p)  

converges. We established this earlier with the Fourier Inequality. This result and the following example should be contrasted with the previous example. In this example, we have normed convergence, but not uniform convergence.

The example is sum(sin(n*x)/n,n = 1 .. infinity) . I try to convince you and me that the convergence is not uniform by seeing that the series converges in norm, but not to a continuous function.

>    f2:=x->sum(sin(n*x)/n,n=1..15):
g2:=x->sum(sin(n*x)/n,n=1..25):
h2:=x->sum(sin(n*x)/n,n=1..35):

>    plot({f2(x),g2(x),h2(x)},x=-Pi..Pi,y=-2..2,
   color=[red,blue,green]);

[Maple Plot]

>   

Here is the last Theorem that we state in this section. It gives a condition for uniform convergence based on properties of the function, not the series  as the last theorem did. After all, it is usually the function we know at the outset, and not properties of the series.

THEOREM:  If f ( x ) is periodic, continuous, and has a sectionally continuous derivative, then the Fourier Series corresponding to f  converges uniformly to f(x)  for the entire real line.

As an example of this Theorem, take the function f ( x ) = | x | on the interval [ -pi , pi  ]. It is continuous on that interval and its periodic extension is also continuous. While its derivative is not continuous, it is sectionally continuous. Further, the coefficients are the ones used in the previous example. Check that.

On the other hand, the function f ( x ) = x  on the interval [-1, 1] is continuous there. While the periodic extension is sectionally  continuous, the periodic extension is not continuous . Thus, you cannot expect the Fourier Series to converge uniformly. At these discontinuities, you should expect the series to converge to

          ( f ( c -) + f (- c +) )/2,

the average of the jump from the left and the jump from the right at the points of discontinuities. Recall the first Theorem in this Section

Test for Failure of a Fourier Series to Converge Uniformly

A commonly known result in undergraduate mathematics is that if a sequence of continuous functions converges uniformly, the the limit must be continuous. A Fourier Series

          S ( N ) = sum(A[p]*sin(p*x),p = 1 .. N)

is a sequence of continuous functions. Indeed, S ( N ) is simply a sum of trigonometric function (and, hence, infinitely differentiable).

 

Thus, if the series converges uniformly, the limit must be continuous. This result is invoked most often if the periodic extension of a function is not continuous. The test for failure to have uniform convergence goes like this: If the periodic extension of a function is not continuous, then the Fourier Series for the function cannot converge uniformly.

Here is an illustration: we take f ( x ) = x  over [-1,1]. Think about the periodic extension, with period 2. This periodic extension is NOT continuous. Thus, the Fourier series cannot converge uniformly.

>    c:=1;
f:=x->x;
a[0]:=int(f(x),x=-c..c)/int(1^2,x=-c..c);
for n from 1 to 5 do
a[n]:=int(f(x)*cos((n*Pi*x)/c),x=-c..c)/
          int(cos(n*Pi*x/c)^2,x=-c..c);
   b[n]:=int(f(x)*sin((n*Pi*x)/c),x=-c..c)/
          int(sin(n*Pi*x/c)^2,x=-c..c);
od;
n:='n':

c := 1

f := proc (x) options operator, arrow; x end proc

a[0] := 0

a[1] := 0

b[1] := 2/Pi

a[2] := 0

b[2] := -1/Pi

a[3] := 0

b[3] := 2/3/Pi

a[4] := 0

b[4] := -1/(2*Pi)

a[5] := 0

b[5] := 2/5/Pi

>    s:=x->a[0]+sum(a[n]*cos((n*Pi*x)/c)+b[n]*sin((n*Pi*x)/c),n=1..5);

s := proc (x) options operator, arrow; a[0]+sum(a[n]*cos(n*Pi*x/c)+b[n]*sin(n*Pi*x/c),n = 1 .. 5) end proc

>    plot([[x,f(x),x=-c..c],[x,s(x),x=-2*c..2*c]]);

[Maple Plot]

>   

>   

We have presented two type results in this Section. One type result was information about  convergence based on the character of the coefficients. The other result was based on the character of the function to which the series was to converge. The most frequent application is the last result: if the function is periodic, piecewise smooth, and continuous, the Fourier Trigonometric Series converges uniformly. If the function is not continuous, the series will not converge uniformly.

EMAIL: herod@math.gatech.edu   or    jherod@tds.net

URL: http://www.math.gatech.edu/~herod

Copyright ©  2003  by James V. Herod

All rights reserved