Differential Equations Project

Computer Lab Modules

 

Fourier Series Expansions

 

 

Background. Fourier series are used to expand periodic functions in the trigonometric form.

Theorem 1. (Fourier Expansion) If f(t) has period [Graphics:e24.txtgr1.gif] and f '(t) is piecewise continuous, the Fourier expansion is
[Graphics:e24.txtgr2.gif]

Theorem 2. (Fourier Cosine Series) Assume that f(t) is an even function and has period [Graphics:e24.txtgr3.gif]. If f(t) and f '(t) are piecewise continuous, the Fourier series for f(t) involves only the cosine terms, (i.e. [Graphics:e24.txtgr4.gif]):
[Graphics:e24.txtgr5.gif]

Theorem 3. (Fourier Sine Series) Assume that f(t) is an odd function and has period [Graphics:e24.txtgr6.gif]. If f(t) and f '(t) are piecewise continuous, the Fourier series for f(t) involves only the sine terms, (i.e. [Graphics:e24.txtgr7.gif]):
[Graphics:e24.txtgr8.gif]

Load Mathematica's FourierTransform package.

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr9.gif]
 

Computer Lab Work.

Exercise 1. Find the Fourier series expansion for [Graphics:e24.txtgr11.gif] extended periodically with period [Graphics:e24.txtgr12.gif].

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr13.gif]

Notice that the function f[t] is odd, so that the coefficients [Graphics:e24.txtgr14.gif] are all zero. The following computations assist with the computation of the coefficients [Graphics:e24.txtgr15.gif].

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr16.gif]

The first six coefficients are:

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr17.gif]

Plot the graph of [Graphics:e24.txtgr18.gif].

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr19.gif]

We can use Mathematica's built in Fourier series procedure to perform our computations, for Exercise:

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr20.gif]

The pattern we have for the Fourier trigonometric polynomials are given by the following summation which we can check out for the case where 5 terms are added in the sum.

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr21.gif]

Mathematica can sum the infinite Fourier series and obtain a closed form for the "periodic extension" g(t) of f(t).

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr22.gif]


Exercise 2. Find the Fourier series expansion for [Graphics:e24.txtgr23.gif] extended periodically with period [Graphics:e24.txtgr24.gif].

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr25.gif]

Notice that the function f[t] is even, so that the coefficients [Graphics:e24.txtgr26.gif] are all zero. The following computations assist with the computation of the coefficients [Graphics:e24.txtgr27.gif].

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr28.gif]

The first six coefficients are:

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr29.gif]

Since [Graphics:e24.txtgr30.gif] we will plot the Fourier trigonometric polynomials corresponding to n=1,3,5.

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr31.gif]

The pattern we have for the Fourier trigonometric polynomials are given by the following summation which we can check out for the case where 3 terms are added in the sum.

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr32.gif]

Mathematica can sum the infinite Fourier series and obtain a closed form for the "periodic extension" g(t) of f(t).

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr33.gif]


Exercise 3. Find the Fourier series expansion for [Graphics:e24.txtgr34.gif] extended periodically with period [Graphics:e24.txtgr35.gif].

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr36.gif]

Notice that the function f[t] is odd, so that the coefficients [Graphics:e24.txtgr37.gif] are all zero. The following computations assist with the computation of the coefficients [Graphics:e24.txtgr38.gif].

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr39.gif]

The first seven coefficients are:

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr40.gif]

Since [Graphics:e24.txtgr41.gif] we will plot the Fourier trigonometric polynomials corresponding to n=1,3,5,7.

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr42.gif]

The pattern we have for the Fourier trigonometric polynomials are given by the following summation which we can check out for the case where 3 terms are added in the sum.

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr43.gif]

Mathematica can sum the infinite Fourier series and obtain a closed form for the "periodic extension" g(t) of f(t).

[Graphics:e24.txtgr10.gif][Graphics:e24.txtgr44.gif]
 

Solutions.

 

 

Return to the Differential Equations Project

Return to the Numerical Analysis Project

Return to the Complex Analysis Project

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews, 1998