Fourier Series Intro - Fourier Series


http://mathworld.wolfram.com/FourierSeries.html


Fourier Series

DOWNLOAD Mathematica NotebookEXPLORE THIS TOPIC IN the MathWorld ClassroomFourierSeriesExamples

A Fourier series is an expansion of a periodic function f(x) in terms of an infinite sum ofsines and cosines. Fourier series make use of the orthogonality relationships of the sine and cosine functions. The computation and study of Fourier series is known asharmonic analysis and is extremely useful as a way to break up anarbitrary periodic function into a set of simple terms that can be plugged in, solved individually, and then recombined to obtain the solution to the original problem or an approximation to it to whatever accuracy is desired or practical. Examples of successive approximations to common functions using Fourier series are illustrated above.

In particular, since the superposition principle holds for solutions of a linear homogeneous ordinary differential equation, if such an equation can be solved in the case of a single sinusoid, the solution for an arbitrary function is immediately available by expressing the original function as a Fourier series and then plugging in the solution for each sinusoidal component. In some special cases where the Fourier series can be summed in closed form, this technique can even yield analytic solutions.

Any set of functions that form a complete orthogonal system have a corresponding generalized Fourier series analogous to the Fourier series. For example, using orthogonality of the roots of aBessel function of the first kind gives a so-called Fourier-Bessel series.

The computation of the (usual) Fourier series is based on the integral identities

int_(-pi)^pisin(mx)sin(nx)dx=pidelta_(mn)
(1)
int_(-pi)^picos(mx)cos(nx)dx=pidelta_(mn)
(2)
int_(-pi)^pisin(mx)cos(nx)dx=0
(3)
int_(-pi)^pisin(mx)dx=0
(4)
int_(-pi)^picos(mx)dx=0
(5)

for m,n!=0, wheredelta_(mn) is theKronecker delta.

Using the method for a generalized Fourier series, the usual Fourier series involving sines and cosines is obtained by takingf_1(x)=cosx andf_2(x)=sinx. Since these functions form acomplete orthogonal system over [-pi,pi], the Fourier series of a functionf(x) is given by

 f(x)=1/2a_0+sum_(n=1)^inftya_ncos(nx)+sum_(n=1)^inftyb_nsin(nx),
(6)

where

a_0=1/piint_(-pi)^pif(x)dx
(7)
a_n=1/piint_(-pi)^pif(x)cos(nx)dx
(8)
b_n=1/piint_(-pi)^pif(x)sin(nx)dx
(9)

and n=1, 2, 3, .... Note that the coefficient of the constant terma_0 has been written in a special form compared to the general form for ageneralized Fourier series in order to preserve symmetry with the definitions ofa_n andb_n.

The Fourier cosine coefficient a_n and sine coefficientb_n are implemented in theWolfram Language asFourierCosCoefficient[expr,t,n] and FourierSinCoefficient[expr,t,n], respectively.

A Fourier series converges to the function f^_ (equal to the original function at points of continuity or to the average of the two limits at points of discontinuity)

 f^_={1/2[lim_(x->x_0^-)f(x)+lim_(x->x_0^+)f(x)]   for -pi<x_0<pi; 1/2[lim_(x->-pi^+)f(x)+lim_(x->pi_-)f(x)]   for x_0=-pi,pi
(10)

if the function satisfies so-called Dirichlet boundary conditions. Dini's test gives a condition for the convergence of Fourier series.

FourierSeriesSquareWave

As a result, near points of discontinuity, a "ringing" known as the Gibbs phenomenon, illustrated above, can occur.

For a function f(x) periodic on an interval[-L,L] instead of[-pi,pi], a simple change of variables can be used to transform the interval of integration from[-pi,pi] to[-L,L]. Let

x=(pix^')/L
(11)
dx=(pidx^')/L.
(12)

Solving for x^' givesx^'=Lx/pi, and plugging this in gives

 f(x^')=1/2a_0+sum_(n=1)^inftya_ncos((npix^')/L)+sum_(n=1)^inftyb_nsin((npix^')/L).
(13)

Therefore,

a_0=1/Lint_(-L)^Lf(x^')dx^'
(14)
a_n=1/Lint_(-L)^Lf(x^')cos((npix^')/L)dx^'
(15)
b_n=1/Lint_(-L)^Lf(x^')sin((npix^')/L)dx^'.
(16)

Similarly, the function is instead defined on the interval [0,2L], the above equations simply become

a_0=1/Lint_0^(2L)f(x^')dx^'
(17)
a_n=1/Lint_0^(2L)f(x^')cos((npix^')/L)dx^'
(18)
b_n=1/Lint_0^(2L)f(x^')sin((npix^')/L)dx^'.
(19)

In fact, for f(x) periodic with period2L,any interval (x_0,x_0+2L) can be used, with the choice being one of convenience or personal preference (Arfken 1985, p. 769).

The coefficients for Fourier series expansions of a few common functions are given in Beyer (1987, pp. 411-412) and Byerly (1959, p. 51). One of the most common functions usually analyzed by this technique is thesquare wave. The Fourier series for a few common functions are summarized in the table below.

functionf(x)Fourier series
Fourier series--sawtooth wavex/(2L)1/2-1/pisum_(n=1)^(infty)1/nsin((npix)/L)
Fourier series--square wave2[H(x/L)-H(x/L-1)]-14/pisum_(n=1,3,5,...)^(infty)1/nsin((npix)/L)
Fourier series--triangle waveT(x)8/(pi^2)sum_(n=1,3,5,...)^(infty)((-1)^((n-1)/2))/(n^2)sin((npix)/L)

If a function is even so that f(x)=f(-x), thenf(x)sin(nx) isodd. (This follows sincesin(nx) isodd and aneven function times anodd function is anodd function.) Therefore,b_n=0 for alln. Similarly, if a function isodd so thatf(x)=-f(-x), thenf(x)cos(nx) isodd. (This follows sincecos(nx) iseven and aneven function times anodd function is anodd function.) Therefore,a_n=0 for alln.

The notion of a Fourier series can also be extended to complex coefficients. Consider a real-valued function f(x). Write

 f(x)=sum_(n=-infty)^inftyA_ne^(inx).
(20)

Now examine

int_(-pi)^pif(x)e^(-imx)dx=int_(-pi)^pi(sum_(n=-infty)^(infty)A_ne^(inx))e^(-imx)dx
(21)
=sum_(n=-infty)^(infty)A_nint_(-pi)^pie^(i(n-m)x)dx
(22)
=sum_(n=-infty)^(infty)A_nint_(-pi)^pi{cos[(n-m)x]+isin[(n-m)x]}dx
(23)
=sum_(n=-infty)^(infty)A_n2pidelta_(mn)
(24)
=2piA_m,
(25)

so

 A_n=1/(2pi)int_(-pi)^pif(x)e^(-inx)dx.
(26)

The coefficients can be expressed in terms of those in the Fourier series

A_n=1/(2pi)int_(-pi)^pif(x)[cos(nx)-isin(nx)]dx
(27)
={1/(2pi)int_(-pi)^pif(x)[cos(nx)+isin(|n|x)]dx n<0; 1/(2pi)int_(-pi)^pif(x)dx n=0; 1/(2pi)int_(-pi)^pif(x)[cos(nx)-isin(nx)]dx n>0
(28)
={1/2(a_n+ib_n) for n<0; 1/2a_0 for n=0; 1/2(a_n-ib_n) for n>0.
(29)

For a function periodic in [-L/2,L/2], these become

f(x)=sum_(n=-infty)^(infty)A_ne^(i(2pinx/L))
(30)
A_n=1/Lint_(-L/2)^(L/2)f(x)e^(-i(2pinx/L))dx.
(31)

These equations are the basis for the extremely important Fourier transform, which is obtained by transforming A_n from a discrete variable to a continuous one as the length L->infty.

The complex Fourier coefficient is implemented in the Wolfram Language as FourierCoefficient[expr,t,n].



  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值