Applications of Fourier Series to Differential Equations

Trigonometry

Trigonometry Logo

Applications of Fourier Series to Differential Equations

Fourier theory was initially invented to solve certain differential equations. Therefore, it is of no surprise that Fourier series are widely used for seeking solutions to various ordinary differential equations (ODEs) and partial differential equations (PDEs).

In this section, we consider applications of Fourier series to the solution of ODEs and the most well-known PDEs:

  • the heat equation
    \[\frac{{\partial u}}{{\partial t}} = k \frac{{{\partial ^2}u}}{{\partial {x^2}}};\]
  • the wave equation
    \[\frac{{{\partial ^2}u}}{{\partial {t^2}}} = {a^2} \frac{{{\partial ^2}u}}{{\partial {x^2}}};\]
  • Laplace's equation
    \[\frac{{{\partial ^2}u}}{{\partial {x^2}}} + \frac{{{\partial ^2}u}}{{\partial {y^2}}} = 0.\]

Solved Problems

Click or tap a problem to see the solution.

Example 1

Find the Fourier series solution to the differential equation \[y^{\prime\prime} + 2y = 3x\] with the boundary conditions \(y\left( 0 \right) = y\left( 1 \right) = 0.\)

Example 2

Find the periodic solutions of the differential equation \[y' + ky = f\left( x \right),\] where \(k\) is a constant and \(f\left( x \right)\) is a \(2\pi\)-periodic function.

Example 1.

Find the Fourier series solution to the differential equation \[y^{\prime\prime} + 2y = 3x\] with the boundary conditions \(y\left( 0 \right) = y\left( 1 \right) = 0.\)

Solution.

We will use the Fourier sine series for representation of the nonhomogeneous solution to satisfy the boundary conditions. Using the results of Example 3 on the page Definition of Fourier Series and Typical Examples, we can write the right side of the equation as the series

\[3x = \frac{6}{\pi }\sum\limits_{n = 1}^\infty {\frac{{{{\left( { - 1} \right)}^{n + 1}}}}{n}\sin n\pi x} .\]

We assume that the solution has the form

\[y\left( x \right) = \sum\limits_{n = 1}^\infty {{b_n}\sin n\pi x} .\]

Substituting this into the differential equation, we get

\[\sum\limits_{n = 1}^\infty {\left( { - {n^2}{\pi ^2}} \right){b_n}\sin n\pi x} + 2\sum\limits_{n = 1}^\infty {{b_n}\sin n\pi x} = \frac{6}{\pi }\sum\limits_{n = 1}^\infty {\frac{{{{\left( { - 1} \right)}^{n + 1}}}}{n}\sin n\pi x} .\]

Since the coefficients of each sine mode must be equal to each other, we obtain the algebraic equation

\[\left( {2 - {n^2}{\pi ^2}} \right){b_n} = \frac{{6{{\left( { - 1} \right)}^{n + 1}}}}{{n\pi }}\;\; \text{or}\;\;{b_n} = \frac{{6{{\left( { - 1} \right)}^{n + 1}}}}{{n\pi \left( {2 - {n^2}{\pi ^2}} \right)}}.\]

Hence, the solution of the given differential equation is described by the series

\[y\left( x \right) = \frac{6}{\pi }\sum\limits_{n = 1}^\infty {\frac{{{{\left( { - 1} \right)}^{n + 1}}}}{{n\left( {2 - {n^2}{\pi ^2}} \right)}}\sin n\pi x} .\]

Example 2.

Find the periodic solutions of the differential equation \[y' + ky = f\left( x \right),\] where \(k\) is a constant and \(f\left( x \right)\) is a \(2\pi\)-periodic function.

Solution.

We represent the function \(f\left( x \right)\) on the right-hand side of the equation as a Fourier series:

\[f\left( x \right) = \sum\limits_{n = - \infty }^\infty {{c_n}{e^{inx}}} .\]

The complex Fourier coefficients are defined by the formula

\[{c_n} = \frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {f\left( x \right){e^{ - inx}}dx} .\]

Assuming that the solution can be represented as a Fourier series expansion

\[y = \sum\limits_{n = - \infty }^\infty {{y_n}{e^{inx}}} ,\]

we find the expression for the derivative:

\[y' = \sum\limits_{n = - \infty }^\infty {in{y_n}{e^{inx}}} .\]

Substituting this into the differential equation, we get

\[\sum\limits_{n = - \infty }^\infty {in{y_n}{e^{inx}}} + k\sum\limits_{n = - \infty }^\infty {{y_n}{e^{inx}}} = \sum\limits_{n = - \infty }^\infty {{c_n}{e^{inx}}} .\]

As this equation is valid for all \(n,\) we obtain

\[in{y_n} + k{y_n} = {c_n}\;\; \text{or}\;\;{y_n} = \frac{{{c_n}}}{{in + k}}.\]

Here \({c_n}\) and \(k\) are known numbers. Consequently, the solution is given by

\[y\left( x \right) = \sum\limits_{n = - \infty }^\infty {\frac{{{c_n}}}{{in + k}}{e^{inx}}} .\]