Fourier Series

Fourier Series
Part 2: Convergence Issues
So far we have proved the following theorem:
∞
Theorem: If 𝑓(π‘₯) = π‘Ž0 + βˆ‘βˆž
𝑛=1 π‘Žπ‘› cos(𝑛π‘₯) + βˆ‘π‘›=1 𝑏𝑛 sin(𝑛π‘₯) for π‘₯ ∈ [βˆ’πœ‹, πœ‹] then
π‘Ž0 =
1 πœ‹
∫ 𝑓(π‘₯)𝑑π‘₯ ,
2πœ‹ βˆ’πœ‹
π‘Žπ‘› =
1 πœ‹
∫ 𝑓(π‘₯) cos(𝑛π‘₯) 𝑑π‘₯ ,
πœ‹ βˆ’πœ‹
𝑏𝑛 =
1 πœ‹
∫ 𝑓(π‘₯) sin(𝑛π‘₯) 𝑑π‘₯
πœ‹ βˆ’πœ‹
In the previous segment we used these formulas to compute a number of Fourier series. In this segment we want to
switch to a more theoretical discussion and we want to figure out the convergence behavior of a Fourier series. There
are two parts to this question:
1. Find out for which functions a Fourier series converges
2. Determine for which functions a Fourier series converges to the original function
Let’s check the first question first (of course, or I would not have called question 1): Let’s start with some function 𝑓(π‘₯):
ο‚·
if π‘Žπ‘› =
ο‚·
if 𝑏𝑛 =
1 πœ‹
∫ 𝑓(π‘₯) cos(𝑛π‘₯) 𝑑π‘₯ , does βˆ‘βˆž
𝑛=1 π‘Žπ‘› cos(𝑛π‘₯) converge
πœ‹ βˆ’πœ‹
1 πœ‹
∫ 𝑓(π‘₯) sin(𝑛π‘₯) 𝑑π‘₯ , does βˆ‘βˆž
𝑛=1 𝑏𝑛 sin(𝑛π‘₯) converge
πœ‹ βˆ’πœ‹
The series are, on one level, just some infinite series, so we could apply the various convergence tests we know. Let’s
start with the divergence test:
Proposition 1: If 𝑓(π‘₯) is a 𝐢 1 ([βˆ’πœ‹, πœ‹]) function (i.e. f is differentiable and f’ is continuous), then the divergence test fails.
That means, of course, that lim π‘Žπ‘› cos(𝑛π‘₯) = 0 and lim 𝑏𝑛 sin(𝑛π‘₯) = 0, and tells us nothing about convergence of
π‘›β†’βˆž
π‘›β†’βˆž
the series.
1
πœ‹
To show this, let’s use integration by parts on the coefficients π‘Žπ‘› = πœ‹ βˆ«βˆ’πœ‹ 𝑓(π‘₯)cos(𝑛π‘₯) 𝑑π‘₯. Let 𝑒(π‘₯) = f(x) and
1
𝑣 β€² (π‘₯) = cos(𝑛π‘₯). Then 𝑒′ (π‘₯) = 𝑓′(π‘₯) and 𝑣(π‘₯) = 𝑛 sin(𝑛π‘₯) so that:
πœ‹
∫ 𝑓(π‘₯)cos(𝑛π‘₯) 𝑑π‘₯ =
βˆ’πœ‹
πœ‹
1
1 πœ‹
1 πœ‹
𝑓(π‘₯) sin(𝑛π‘₯)| βˆ’ ∫ 𝑓 β€² (π‘₯) sin(𝑛π‘₯) 𝑑π‘₯ = 0 βˆ’ ∫ 𝑓 β€² (π‘₯) sin(𝑛π‘₯) 𝑑π‘₯
𝑛
𝑛 βˆ’πœ‹
𝑛 βˆ’πœ‹
βˆ’πœ‹
Since 𝑓′ is continuous on a closed, bounded interval, it is bounded by, say, 𝑀. But then
|π‘Žπ‘› | ≀
11
2𝑀
𝑀 2πœ‹ =
π‘›πœ‹
𝑛
so that lim π‘Žπ‘› = 0. But then lim π‘Žπ‘› cos(𝑛π‘₯) = 0 for all x because |cos(𝑛 π‘₯)| ≀ 1, so that the divergence test
π‘›β†’βˆž
π‘›β†’βˆž
fails.
The proof for the other series is similar – but not quite: make sure to try it!
Based on this proof it is now easy to show that:
Proposition 2: Suppose 𝑓 is a 𝐢 2 [βˆ’πœ‹, πœ‹] function (i.e. 𝑓 is twice differentiable and 𝑓′′ is continuous). Then the Fourier
∞
series π‘Ž0 + βˆ‘βˆž
𝑛=1 π‘Žπ‘› cos(𝑛π‘₯) + βˆ‘π‘›=1 𝑏𝑛 sin(𝑛π‘₯) with π‘Ž0 =
1 πœ‹
∫ 𝑓(π‘₯) sin(𝑛π‘₯) 𝑑π‘₯
πœ‹ βˆ’πœ‹
1 πœ‹
∫ 𝑓(π‘₯)𝑑π‘₯ ,
2πœ‹ βˆ’πœ‹
π‘Žπ‘› =
1
πœ‹
πœ‹
βˆ«βˆ’πœ‹ 𝑓(π‘₯) cos(𝑛π‘₯) 𝑑π‘₯ , and 𝑏𝑛 =
converges uniformly.
Note that this theorem will be easy to prove but it is not satisfying: for one thing, it does not show whether the series
converges to the original function, and for another we know from our examples that the Fourier series for a function 𝑓
seems to converge (to the original function) for functions that are not even continuous, while this theorem assumes 𝑓 is
twice continuously differentiable. But, it is a start, so let’s go prove it.
Proof: We already know that
πœ‹
1 πœ‹
∫ 𝑓(π‘₯)cos(𝑛π‘₯) 𝑑π‘₯ = βˆ’ ∫ 𝑓 β€² (π‘₯) sin(𝑛π‘₯) 𝑑π‘₯
𝑛 βˆ’πœ‹
βˆ’πœ‹
Integrate by parts again, with 𝑒(π‘₯) = 𝑓 β€² (π‘₯), 𝑣 β€² (π‘₯) = sin(𝑛π‘₯), and, consequently, 𝑒′ (π‘₯) = 𝑓′′(π‘₯), 𝑣(π‘₯) =
1
𝑛
βˆ’ cos(𝑛π‘₯):
πœ‹
πœ‹
1 πœ‹
1
1
1 πœ‹
∫ 𝑓(π‘₯) cos(𝑛π‘₯) 𝑑π‘₯ = βˆ’ ∫ 𝑓 β€² (π‘₯) sin(𝑛π‘₯) 𝑑π‘₯ = βˆ’ (βˆ’ 𝑓 β€² (π‘₯) cos(𝑛π‘₯)| + ∫ 𝑓 β€²β€² (π‘₯) cos(𝑛π‘₯) 𝑑π‘₯) =
𝑛 βˆ’πœ‹
𝑛
𝑛
𝑛 βˆ’πœ‹
βˆ’πœ‹
βˆ’πœ‹
πœ‹
1
1
= 2 (𝑓 β€² (πœ‹) βˆ’ 𝑓 β€² (βˆ’πœ‹)) cos(𝑛 πœ‹) βˆ’ 2 ∫ 𝑓 β€²β€² (π‘₯) cos(𝑛π‘₯) 𝑑π‘₯
𝑛
𝑛 βˆ’πœ‹
Therefore:
𝐾
|π‘Žπ‘› | ≀ 2
𝑛
𝐾
for some constant K, so that βˆ‘βˆž
𝑛=1 π‘Žπ‘› cos(𝑛π‘₯) converges by the comparision test with βˆ‘ 𝑛2 , a convergent p-series
with 𝑝 = 2, and the convergence is uniform by Weierstrass Convergence theorem.
The proof that βˆ‘βˆž
𝑛=1 𝑏𝑛 sin(𝑛π‘₯) converges is similar. But then the original Fourier series must converge.
The above theorem is nice and easy to prove, but using Abel’s convergence test for series we can do slightly better:
Proposition 3: If 𝑓(π‘₯) is a 𝐢 1 ([βˆ’πœ‹, πœ‹]) function (i.e. f is differentiable and f’ is continuous), then the Fourier series for 𝑓
converges.
Apply Abel’s convergence test (see IRA 4.2 for details). The details are left as an exercise.
But we can do a little better still. Here is a β€˜sharp’ result for the convergence of a Fourier series. The proof is a little
detailed so we are going to skip it.
Theorem: (Fourier Series)
Suppose 𝑓 is a piecewise 𝐢 1 [βˆ’πœ‹, πœ‹] function, i.e. 𝑓 consists of finitely many parts, each of which is continuously
differentiable, and all discontinuities are jump discontinuities, then the Fourier series of 𝑓 converges to 𝑓(π‘₯) at
all points where 𝑓 is continuous, and to 1/2(𝑓(π‘₯ + ) βˆ’ 𝑓(π‘₯ βˆ’ )) at each jump.
Example: Check the convergence behavior of the Fourier series of 𝑓(π‘₯) = {
πœ‹, 0 < π‘₯ ≀ πœ‹
0, βˆ’πœ‹ ≀ π‘₯ ≀ 0
Pay special attention to the limit of the Fourier series at the jump at π‘₯ = 0
πœ‹
2
In the previous segment we worked out the Fourier series for this function to be 𝐹(π‘₯) = 2 + βˆ‘βˆž
𝑛=0 2𝑛+1 sin((2𝑛 + 1) π‘₯).
πœ‹
1
1
πœ‹
Clearly 𝐹(0) = 2 = 2 (𝑓(0+ ) βˆ’ 𝑓(0βˆ’ )) = 2 (πœ‹ βˆ’ 0) = 2 , so that the series does converge at least at the jump as
advertised. And it seems to converge to 𝑓(π‘₯) elsewhere as well, judging by several graphs of the N-th partial Fourier
series with 𝑁 = 3, 𝑁 = 10, and 𝑁 = 40, respectively.
3
2
1
3
3
3
2
2
2
1
1
1
1
2
3
3
2
1
1
2
3
3
2
1
1
2
3
Note how the Fourier series β€˜overshoots’ near the jump but then quickly dampens down to hug the original function
closely. That behavior is typical for Fourier series and has a number of interesting consequences in the physics of signal
transmission.
Anyway, we have achieved our goal: the Fourier Theorem above tells us for which functions the Fourier Series converges
back to the original function.
We do have one pesky limitation left: each of the function we considered so far was defined on the interval [βˆ’πœ‹, πœ‹].
While we can extend the function by periodicity to all real numbers and the Fourier series would converge to that
extended periodic function for all π‘₯, what about functions defined on an arbitray interval [a, b]? We will check that out
in our next segment.
Exercises:
1. Prove the second part of Proposition 1
2. Prove the second part of Proposition 2
1
3. Prove Proposition 3. Hint: Prove and then use the fact that | βˆ‘π‘˜π‘›=1 sin(𝑛π‘₯) | ≀ |sin(π‘₯/2)| together with Abel’s
convergence test.
πœ‹
4
4. Use the Fourier series for 𝑓(π‘₯) = π‘₯ on [βˆ’πœ‹, πœ‹] to show that = βˆ‘βˆž
𝑛=0
(βˆ’1)𝑛
2𝑛+1
(we already knew that by looking at
the Taylor series for π‘Žπ‘Ÿπ‘π‘‘π‘Žπ‘›(π‘₯) but our reasoning is different this time around)
5. Use the Fourier series for 𝑓(π‘₯) = |π‘₯| on [βˆ’πœ‹, πœ‹] to show that
πœ‹2
8
1
= βˆ‘βˆž
𝑛=0 (2𝑛+1)2
1
6. Show that βˆ‘βˆž
𝑛=1 𝑛2 =
πœ‹2
6
1, π‘₯ < 0
7. Check the Fourier Theorem for 𝑓(π‘₯) = {
, especially at the jump π‘₯ = 0. Verify the convergence behavior
π‘₯, π‘₯ β‰₯ 0
graphically, using Mathematica.