Half range Fourier series

In mathematics, a half range Fourier series is a Fourier series defined on an interval ${\displaystyle [0,L]}$ instead of the more common ${\displaystyle [-L,L]}$, with the implication that the analyzed function ${\displaystyle f(x),x\in [0,L]}$ should be extended to ${\displaystyle [-L,0]}$ as either an even (${\displaystyle f(-x)=f(x)}$) or odd function (${\displaystyle f(-x)=-f(x)}$).^[1] This allows the expansion of the function in a series solely of sines (odd) or cosines (even). The choice between odd and even is typically motivated by boundary conditions associated with a differential equation satisfied by ${\displaystyle f(x)}$.^[2]

Given a function ${\displaystyle f(x)}$ defined over the interval ${\displaystyle [0,L]}$, the half-range Fourier series extends ${\displaystyle f}$ to ${\displaystyle [-L,L]}$ with either an even or an odd extension. An even extension satisfies the identity

${\displaystyle f(-x)=f(x)\quad {\text{for}}\quad x\in [0,L]}$

And an odd extension satisfies the identity

${\displaystyle f(-x)=-f(x)\quad {\text{for}}\quad x\in [0,L]}$

An even extension relies on the cosine function since the cosine function is even, and an odd extension uses the sine function since the sine function is odd. The decision between even and odd extensions is usually motivated by the boundary conditions of the problem, especially with partial differential equations such as heat and wave equations.^[2]

Let ${\displaystyle f(x)}$ be defined and piecewise continuous on the interval ${\displaystyle [0,L]}$. The half range Fourier series for sine and cosine are given as follows:^[3]

${\displaystyle {\text{Sine series: }}f(x)\sim \sum _{n=1}^{\infty }b_{n}\sin \left({\frac {n\pi x}{L}}\right)}$

And the coefficients ${\displaystyle b_{n}}$ are given by

${\displaystyle b_{n}={\frac {2}{L}}\int _{0}^{L}f(x)\sin \left({\frac {n\pi x}{L}}\right)dx}$

The cosine series is given as

${\displaystyle {\text{Cosine series: }}f(x)\sim {\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }a_{n}\cos \left({\frac {n\pi x}{L}}\right)}$

And the coefficients ${\displaystyle a_{n}}$ are given by

${\displaystyle a_{0}={\frac {2}{L}}\int _{0}^{L}f(x)dx,\quad a_{n}={\frac {2}{L}}\int _{0}^{L}f(x)\cos \left({\frac {n\pi x}{L}}\right)dx}$

Half range Fourier series are used in solving certain boundary value problems involving heat conduction, vibration, and other physical phenomena modeled by partial differential equations.^[4] The advantage in a half range Fourier series is that the function analyzed need only be defined over ${\displaystyle [0,L]}$, meaning functions only defined over positive numbers can be analyzed with a half range Fourier series. Another advantage is using only sines or only cosines. For example, if ${\displaystyle f(0)=0}$, it can be more practical to use a sine series rather than a full Fourier series. If ${\displaystyle f'(0)=0}$, it can be more practical to use a cosine series instead.^[5]

Example 1.

The half range Fourier sine series expansion for ${\displaystyle f(x)=\cos(x)}$ with ${\displaystyle x\in (0,\pi )}$ is computed as follows. Because we want a sine series, it will take the form:

${\displaystyle f(x)\sim \sum _{n=1}^{\infty }b_{n}\sin(nx)}$

With coefficients given by:

${\displaystyle b_{n}={\frac {2}{\pi }}\int _{0}^{\pi }f(x)\sin(nx)\,dx}$

Noting that because ${\displaystyle \pi =L}$, these terms cancel out giving only ${\displaystyle \sin(nx)}$. We substitute the definition of ${\displaystyle f(x)}$ and obtain

${\displaystyle b_{n}={\frac {2}{\pi }}\int _{0}^{\pi }\cos(x)\sin(nx)\,dx}$

Using the product to sum identity given by ${\displaystyle \sin(A)\cos(B)={\frac {\sin(A+B)+\sin(A-B)}{2}}}$, we convert the integrand into

${\displaystyle {\begin{aligned}b_{n}={\frac {2}{\pi }}\int _{0}^{\pi }{\frac {\sin((n+1)x)+\sin((n-1)x)}{2}}\,dx=\\{\frac {2}{\pi }}\cdot {\frac {1}{2}}\int _{0}^{\pi }{\Big [}\sin((n+1)x)+\sin((n-1)x){\Big ]}\end{aligned}}}$

Applying the fact that for ${\displaystyle m\neq 0}$, ${\displaystyle \int _{0}^{\pi }\sin(mx)\,dx={\frac {1-\cos(m\pi )}{m}}}$, and that ${\displaystyle \cos(m\pi )=(-1)^{m}}$, we have

${\displaystyle \int _{0}^{\pi }\sin(mx)\,dx={\frac {1-(-1)^{m}}{2}}}$

Applying to both terms of the integrand we have

${\displaystyle b_{n}={\frac {1}{\pi }}\left[{\frac {1-(-1)^{n+1}}{n+1}}+{\frac {1-(-1)^{n-1}}{n-1}}\right],\quad n\neq 1}$

If ${\displaystyle n}$ is odd, then both ${\displaystyle n+1}$ and ${\displaystyle n-1}$ are even, meaning ${\displaystyle (-1)^{n+1}=(-1)^{n-1}=1}$, and thus the denominators are equal to ${\displaystyle 0}$. Thus, we have

${\displaystyle b_{n}=0{\text{ when }}n{\text{ is odd}}}$

We apply this to the special case of ${\displaystyle b_{1}}$ as well. Then, if ${\displaystyle n}$ is even, then ${\displaystyle n=2k}$, and thus ${\displaystyle n-1}$ and ${\displaystyle n+1}$ are odd, meaning ${\displaystyle (-1)^{n-1}=(-1)^{n+1}=-1}$. Thus, both denominators become ${\displaystyle 2}$, and we obtain

${\displaystyle b_{2k}={\frac {1}{\pi }}\left[{\frac {2}{2k+1}}+{\frac {2}{2k-1}}\right]={\frac {2}{\pi }}\left[{\frac {1}{2k+1}}+{\frac {1}{2k-1}}\right]}$

Combining into a single fraction, we have

${\displaystyle {\frac {1}{2k+1}}+{\frac {1}{2k-1}}={\frac {2k+1+2k-1}{(2k+1)(2k-1)}}={\frac {4k}{4k^{2}-1}}}$

Therefore:

${\displaystyle b_{2k}={\frac {8k}{\pi (4k^{2}-1)}}}$

Because ${\displaystyle b_{n}=0}$ when ${\displaystyle n}$ is odd, we can substitute ${\displaystyle k}$ directly to get the full summation:

${\displaystyle \cos(x)=\sum _{k=1}^{\infty }{\frac {8k}{\pi (4k^{2}-1)}}\sin(2kx),\quad x\in (0,\pi )}$