Module for Fourier Series for O. D. E.

Numerical Methods for O. D. E.  using Mathematica. (c) John H. Mathews, 2003

Background.

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

Definition  (Fourier Series).  If [Graphics:Images/FourierSeriesDEMod_gr_1.gif] is periodic with period [Graphics:Images/FourierSeriesDEMod_gr_2.gif] and is piecewise continuous on[Graphics:Images/FourierSeriesDEMod_gr_3.gif], then the Fourier Series [Graphics:Images/FourierSeriesDEMod_gr_4.gif] for  [Graphics:Images/FourierSeriesDEMod_gr_5.gif] is

    [Graphics:Images/FourierSeriesDEMod_gr_6.gif],

where the coefficients [Graphics:Images/FourierSeriesDEMod_gr_7.gif]  are given by the so-called Euler's formulae:  

    [Graphics:Images/FourierSeriesDEMod_gr_8.gif],  
    and  
    [Graphics:Images/FourierSeriesDEMod_gr_9.gif].  

Theorem (Fourier Expansion).  Assume that [Graphics:Images/FourierSeriesDEMod_gr_10.gif] is the Fourier Series for [Graphics:Images/FourierSeriesDEMod_gr_11.gif].  If [Graphics:Images/FourierSeriesDEMod_gr_12.gif] are piecewise continuous on [Graphics:Images/FourierSeriesDEMod_gr_13.gif], then [Graphics:Images/FourierSeriesDEMod_gr_14.gif] is convergent for all [Graphics:Images/FourierSeriesDEMod_gr_15.gif].  The relation [Graphics:Images/FourierSeriesDEMod_gr_16.gif] holds for all [Graphics:Images/FourierSeriesDEMod_gr_17.gif]where f is continuous.  If [Graphics:Images/FourierSeriesDEMod_gr_18.gif] is a point of discontinuity of f, then

     [Graphics:Images/FourierSeriesDEMod_gr_19.gif],  

where [Graphics:Images/FourierSeriesDEMod_gr_20.gif] denote the left-hand and right-hand limits, respectively.  With this understanding, we have the Fourier Series expansion:

    [Graphics:Images/FourierSeriesDEMod_gr_21.gif] .  

Theorem (Fourier Cosine Series).  Assume that [Graphics:Images/FourierSeriesDEMod_gr_22.gif] is an even function and has period [Graphics:Images/FourierSeriesDEMod_gr_23.gif].  If [Graphics:Images/FourierSeriesDEMod_gr_24.gif] are piecewise continuous, the Fourier series for   [Graphics:Images/FourierSeriesDEMod_gr_25.gif]  involves only the cosine terms,  [Graphics:Images/FourierSeriesDEMod_gr_26.gif],  and we write  

    [Graphics:Images/FourierSeriesDEMod_gr_27.gif],  

where [Graphics:Images/FourierSeriesDEMod_gr_28.gif].  

Theorem (Fourier Sine Series).  Assume that [Graphics:Images/FourierSeriesDEMod_gr_29.gif] is an odd function and has period [Graphics:Images/FourierSeriesDEMod_gr_30.gif].  If [Graphics:Images/FourierSeriesDEMod_gr_31.gif] are piecewise continuous, the Fourier series for   [Graphics:Images/FourierSeriesDEMod_gr_32.gif]  involves only the sine terms,  [Graphics:Images/FourierSeriesDEMod_gr_33.gif],  and we write  

    [Graphics:Images/FourierSeriesDEMod_gr_34.gif],  

where  [Graphics:Images/FourierSeriesDEMod_gr_35.gif].  

Method I. Fourier Series - Trigonometric Polynomials.  Execute the following sells to set up the procedure S[n,t] for finding the Fourier Series - Trigonometric Polynomial of degree n.

[Graphics:Images/FourierSeriesDEMod_gr_36.gif]
[Graphics:Images/FourierSeriesDEMod_gr_37.gif]
[Graphics:Images/FourierSeriesDEMod_gr_38.gif]
[Graphics:Images/FourierSeriesDEMod_gr_39.gif]

Example 1.  Find the Fourier Series expansion of  [Graphics:Images/FourierSeriesDEMod_gr_40.gif].  

Solution 1.

Method II.  Fourier Series - Trigonometric Polynomials.  

    Modify procedure S[n,t] for finding the Fourier Series - Trigonometric Polynomial of degree n to the case of a piecewise continuous function on the subintervals[Graphics:Images/FourierSeriesDEMod_gr_72.gif]

[Graphics:Images/FourierSeriesDEMod_gr_73.gif]
[Graphics:Images/FourierSeriesDEMod_gr_74.gif]
[Graphics:Images/FourierSeriesDEMod_gr_75.gif]
[Graphics:Images/FourierSeriesDEMod_gr_76.gif]

Example 2.  The function  [Graphics:Images/FourierSeriesDEMod_gr_77.gif]  has the Fourier series representation  

    [Graphics:Images/FourierSeriesDEMod_gr_78.gif].  

Solution 2.

Example 3.  Find the Fourier series expansion for  [Graphics:Images/FourierSeriesDEMod_gr_110.gif]  over [Graphics:Images/FourierSeriesDEMod_gr_111.gif],  extended periodically with period  [Graphics:Images/FourierSeriesDEMod_gr_112.gif].  

Solution 3.

Research Experience for Undergraduates


Converted by Mathematica      January 28, 2003