06maxort.mws

Partial Differential Equations PowerTool

by Dr. Jim Herod

Section 1.6: A Maximal Orthonormal Family

Maple Packages for Section 1.6

>    restart;

>   

In the previous sections, we have developed the notions of an orthogonal family, or of an orthonormal family. We emphasized the importance of these in Fourier Analysis. There is one final notion that we introduce in this section. It is that of a maximal  orthonormal family. It should become apparent what is the importance of this notion in this section.

We review that important Fourier Inequality. Suppose that { E, <x, y> } is a linear space on which there is defined an inner product. If { phi[1] , phi[2] , phi[3] , ... } is an orthonormal  family and f  is in the space, then

      0 ` ` <= ` `    abs(f-sum(a[p]*phi[p],p))^2  = abs(f)^2  + sum(abs(`<,>`(f,phi[p])-a[p])^2,p)   - sum(abs(`<,>`(f,phi[p]))^2,p) .

Several important facts follow from this inequality.

Fact 1:  Suppose that n is a positive integer, that S  is the span of { phi[1] , phi[2] , phi[3] , ... , phi[n]  }, and that f  is in E. The closest element in S  is given by the Fourier Expansion

                                                 approximation for f = ~    sum(`<`*f*`,`*phi[p]*`>`*phi[p],p = 1 .. n)  .

We have argued that this is true in Section 1.3.

Fact 2:  If { phi[1] , phi[2] , phi[3] , ... } is an infinite sequence, then the number series   sum(abs(`<`*f*`,`*phi[p]*`>`)^2,p)  converges.

We argue this as follows: By choosing a[p]  = < f  , phi[p]  > we get Bessel's inequality. From this we see that the number series is increasing and bounded above. That is, it is a series of non-negative numbers and bounded above. Thus, it has to converge.

Fact 3:  If { phi[1] , phi[2] , phi[3] , ... } is an infinite sequence, then the series in E given by

                                                    sum(`<`*f*`,`*phi[p]*`>`*phi[p],p = 1 .. infinity)  

converges in E.

This follows from Fact 2 through a series of steps: Suppose m > n.

      sum(`<`*f*`,`*phi[p]*`>`*phi[p],p = 1 .. m)     -     sum(`<`*f*`,`*phi[p]*`>`*phi[p],p = 1 .. n)     =    sum(`<`*f*`,`*phi[p]*`>`*phi[p],p = n .. m)

Also,

      abs(sum(`<`*f*`,`*phi[p]*`>`*phi[p],p = n .. m))^2   =   sum(` `,p = n .. m) | < f , phi[p]  > | ` `^2   =   sum(` `,p = 1 .. m) | < f , phi[p]  > | ` `^2   -    sum(` `,p = 1 .. n) | < f , phi[p]  > | ` `^2

Using Fact 2, we have the convergence of the vector series.

Remark: The naive novice might think that having an infinite orthonormal family and having that the series of Fact 3 converges implies that it converges to f . Only a minute or two of consideration, however, will yield the folly of this mistake. Here is a specific example. Take the orthonormal sequence to be the infinite sequence consisting of

                           {2 sin(2 Pi  x), 2 sin(4 Pi  x), 2 sin(6 Pi  x), ... }.

Take f  to be

                               f ( x ) = x^2*(1-x) .

Look at how badly the approximation misses.

>    f:=x->x^2*(1-x);

f := proc (x) options operator, arrow; x^2*(1-x) end proc

>    for n from 1 to 10 do
   a[n]:=int(f(x)*sqrt(2)*sin(2*n*Pi*x),x=0..1);
od;

a[1] := -3/4*2^(1/2)/Pi^3

a[2] := -3/32*2^(1/2)/Pi^3

a[3] := -1/36*2^(1/2)/Pi^3

a[4] := -3/256*2^(1/2)/Pi^3

a[5] := -3/500*2^(1/2)/Pi^3

a[6] := -1/288*2^(1/2)/Pi^3

a[7] := -3/1372*2^(1/2)/Pi^3

a[8] := -3/2048*2^(1/2)/Pi^3

a[9] := -1/972*2^(1/2)/Pi^3

a[10] := -3/4000*2^(1/2)/Pi^3

>    approx:=x->sum(a[p]*sqrt(2)*sin(2*p*Pi*x),p=1..10);

approx := proc (x) options operator, arrow; sum(a[p]*sqrt(2)*sin(2*p*Pi*x),p = 1 .. 10) end proc

>    plot([f(x),approx(x)],x=0..1,color=[black,red]);

[Maple Plot]

>   

We missed all the odd terms. There were not enough orthonormal vectors. One needs a maximal  family to be able to guarantee convergence.

Definition:  An orthogonal sequence is maximal , or complete , if the only vector in the space that is orthogonal to every element in the sequence is the zero vector.

Fact 4 . Suppose that { phi[1] , phi[2] , ..., phi[n]  , ...} is a  maximal orthonormal sequence. Then

                                   sum(`<,>`(f,phi[p])*phi[p],p = 1 .. infinity)  

converges to f in norm.

We argue this by showing that the element   g  defined as  f -    sum(`<,>`(f,phi[p])*phi[p],p = 1 .. infinity)   is orthogonal to each phi[p] . Just do the computations:

   < g , phi[n]  > = < f , phi[n]  > -  < f , phi[n]  > = 0.

Here, we have used the fact that the dot product of phi[p]  and phi[n]  is zero, unless p = n. In that case, the dot product is 1. Because the orthonormal sequence was maximal and g is orthogonal to every phi[p] , it must be that g is zero, and the series converged to f .

We need to repair the previous example by taking all the sine terms:

                            {2 sin( Pi  x), 2 sin(2 Pi  x), 2 sin(3 Pi  x), ... }.

>    f:=x->x^2*(1-x);

f := proc (x) options operator, arrow; x^2*(1-x) end proc

>    for n from 1 to 10 do
   a[n]:=int(f(x)*sqrt(2)*sin(n*Pi*x),x=0..1);
od;

a[1] := 2*2^(1/2)/Pi^3

a[2] := -3/4*2^(1/2)/Pi^3

a[3] := 2/27*2^(1/2)/Pi^3

a[4] := -3/32*2^(1/2)/Pi^3

a[5] := 2/125*2^(1/2)/Pi^3

a[6] := -1/36*2^(1/2)/Pi^3

a[7] := 2/343*2^(1/2)/Pi^3

a[8] := -3/256*2^(1/2)/Pi^3

a[9] := 2/729*2^(1/2)/Pi^3

a[10] := -3/500*2^(1/2)/Pi^3

>    approx:=x->sum(a[p]*sqrt(2)*sin(p*Pi*x),p=1..10);

approx := proc (x) options operator, arrow; sum(a[p]*sqrt(2)*sin(p*Pi*x),p = 1 .. 10) end proc

We off set the graph of f just a little so that we can see the difference.

>    plot([f(x)+0.003,approx(x)],x=0..1,color=[black,red]);

[Maple Plot]

>   

This was the culmination of the review of concepts from Linear Algebra that are needed to understand some of the methods for solving partial differential equations in this set of lecture notes. The new and important idea of this Section is the notion of a maximal orthonormal family. With such a family, we can approximate in norm all vectors in the space having the associated dot product.

>   

>   

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