傅里叶级数

的傅里叶级数前四项的部份总和。随着增加更多的谐波,这个部份总和收敛成方波。]]
在数学中,傅里叶级数(,)是把类似波的函数表示成简单谐波的方式。更正式地说,对于满足狄利克雷定理的周期函数,其傅里叶级数是由一组正弦与余弦函数的加权和表示的方法。傅里叶级数与用来找出无周期函数的频率信息的傅里叶变换有密切的关系。

傅里叶级数是傅里叶分析的一个研究分支,也是采样定理原始证明的核心。傅里叶级数在数论、组合数学、信号处理、概率论、统计学、密码学、声学、光学等领域都有着广泛的应用。

歷史
傅里叶级数得名于法国数学家约瑟夫·傅里叶(1768年–1830年),他提出任何函数都可以展开为三角级数。此前数学家欧拉、达朗贝尔和克莱罗,已发现在认定一個函数有三角级数展开后,通过积分方法计算其系数的公式,而拉格朗日等人已经找到了一些非周期函数的三角级数展开。将周期函数分解为简单振荡函数的总和的最早想法,可以追溯至公元前3世紀古代天文學家的均輪和本輪學說。

傅里叶的工作得到了丹尼尔·伯努利的赞助,傅里叶介入三角级数用來解热传导方程,其最初论文雖經、加斯帕尔·蒙日同意,但在1807年经拉格朗日、拉普拉斯和勒讓德评審后被拒绝出版,他的现在被称为的理论后来发表于1822年出版的《热的解析理论》。

定義
[[Image:Fourier series illustration.svg|thumb|right|360px|上图展示了蓝色无周期函数s(x),它只定义在红色区间[0,P]之上。可以认为傅里叶级数,是对下图中最初函数的“周期延拓”的分析,傅里叶级数总是周期函数,即使最初函数s(x)不是周期函数。]]
傅里叶级数可以用不同的形式来表达,下面将周期为P的一个周期函數s(x), \ x\in \R 表达为不同形式的傅里叶级数。

正弦-余弦形式
人们常用\sin(x)與\cos(x)的三角級數來表示s(x),就是将所有n階諧波\sin(\frac{2\pi nx}{P})與\cos(\frac{2\pi nx}{P}),乘以其各自在s(x)中的權重,求得它们的總和;这些n階諧波的權重稱爲傅立葉級數係數,它们可以藉由如下積分來獲得:

{{Equation box 1
|indent =:
|title=傅立葉級數係數
|equation = {{NumBlk||
\begin{aligned}
A_0 &= \frac{1}{P}\int_{P}s(x)\,dx \\
A_n &= \frac{2}{P}\int_{P}s(x)\cos\left(\frac{2\pi nx}{P}\right)\,dx\qquad\text{for } n\geq 1\qquad\\
B_n &= \frac{2}{P}\int_{P}s(x)\sin\left(\frac{2\pi nx}{P}\right)\,dx\qquad\text{for } n\geq 1
\end{aligned}
|}}
|cellpadding=
|border
|border colour = #0073CF
|background colour=#F5FFFA}}

符号\int_{P}表示在选定区间上的积分,典型的选择为[-P/2, P/2]或者[0, P]。注意A_0是函数s(x)的平均值,这个性质扩展到了类似的变换比如傅里叶变换。

通过这些系数定义傅里叶级数为:
{{Equation box 1
|indent =:
|title=傅立葉級數,正弦-余弦形式
|equation = {{NumBlk||
s(x) \sim A_0 + \sum_{n=1}^\infty \left(A_n \cos\left(\frac{2\pi nx}{P} \right) + B_n \sin\left(\frac{2\pi nx}{P}\right)\right)
|}}
|cellpadding=
|border
|border colour = #0073CF
|background colour=#F5FFFA}}

这里使用符号\sim,表示傅里叶级数的求和不一定总是等于s(x)。普遍來說n是理論上趨近於無限大的,但是就算趨近於無限大,對所有的x(例如在某一點上不連續),傅立葉級數也不一定收斂到s(x) 。尽管不收斂的可能性始终存在,在科学和工程领域中经常将中的\sim直接替代为=。

在傅里叶级数系数中的整数索引n,是级数中相应的\cos或\sin,在这个函数的周期P中,形成的圆周(cycle)的数目。因此对应于A_n和B_n的项有着:

  • 波长等于\tfrac{P}{n},并且有着同于x的单位。
  • 频率等于\tfrac{n}{P},并且有着x的倒数单位。

指数形式
下面藉由歐拉公式\ e^{ix} = \cos x+i\sin x\ ,将傅里叶级数系数简化成复数指數形式。

根據定義,我們可以得到:
{{Equation box 1
|border
|indent=:
|title=复数傅里叶级数系数
|equation={{NumBlk||\begin{aligned} c_0 &= A_0 &\\
c_n &= (A_n - i B_n)/2 \qquad &\text{for } n > 0 \\
c_n &= (A_{-n}+iB_{-n})/2 \qquad &\text{for } n |}}
|cellpadding=6|border colour=#0073CF|background colour=#F5FFFA}}

通过将等式代入,可以证实:

{{Equation box 1
|border
|indent=:
|title=复数傅里叶级数系数
|equation=c_n = \frac{1}{P}\int_{P} s(x)e^{-\frac{2\pi inx}{P}}\,dx\qquad\text{for}\ n \in \mathbb{Z}
|cellpadding=6|border colour=#0073CF|background colour=#F5FFFA}}

给定复数傅里叶级数系数,可以用公式复原出A_n和B_n:
{{Equation box 1
|border
|indent=:
|title=复数傅里叶级数系数
|equation=\begin{aligned} A_0 &= c_0 &\\
A_n &= c_n+c_{-n} \qquad &\textrm{for}~ n > 0 \\
B_n &= i(c_n-c_{-n}) \qquad &\textrm{for}~ n > 0\end{aligned}
|cellpadding=6|border colour=#0073CF|background colour=#F5FFFA}}

通过这些定义,傅里叶级数可以写为:

{{Equation box 1
|indent =:
|title=傅立葉級數,指數形式
|equation = {{NumBlk||
s(x) \sim \sum_{n=-\infty}^\infty c_n\cdot e^{\frac{2\pi i nx}{P}}
|}}
|cellpadding= 6
|border
|border colour = #0073CF
|background colour=#F5FFFA}}

这是可推广到复数值域函数的惯用形式。n的负数值对应于负频率。

複數值函數
人们習慣將s(x)的值域普遍化到複數上,设s(x)是一個複數值函數,它的實部和虛部,都是實數值函數:
:s(x) = \operatorname{Re}(s(x)) +i\cdot \operatorname{Im}(s(x)) , \quad x\in \R
定义c_n \triangleq c_{_{Rn}}+i\cdot c_{_{In}}则:
:c_n = \frac{1}{P}\int_{P} s(x)\cdot e^{-i\tfrac{2\pi}{p}nx}\ dx = \frac{1}{P}\int_{P} \operatorname{Re}(s(x))\cdot e^{-i\tfrac{2\pi}{p}nx}\ dx + i\cdot \frac{1}{P} \int_{P} \operatorname{Im}(s(x))\cdot e^{-i\tfrac{2\pi}{p}nx}\ dx
:c_{_{Rn}} = \frac{1}{P}\int_P \operatorname{Re}(s(x))\cdot e^{-i\tfrac{2\pi}{p}nx}\ dx
:c_{_{In}} = \frac{1}{P}\int_P \operatorname{Im}(s(x))\cdot e^{-i\tfrac{2\pi}{p}nx}\ dx
对于这个复数值函数,它的傅里叶级数的实部,是它的实部的傅里叶级数;它的傅里叶级数的虚部,是它的虚部的傅里叶级数:
:s(x) \sim \sum_{n=-\infty}^\infty \left(c_{_{Rn}}+i\cdot c_{_{In}}\right) \cdot e^{i\tfrac{2\pi}{p}nx} = \sum_{n=\infty}^\infty c_{_{Rn}}\cdot e^{i\tfrac{2\pi}{p}nx} + i\cdot \sum_{n=-\infty}^\infty c_{_{In}}\cdot e^{i\tfrac{2\pi}{p}nx}

振幅-相位形式
還可以利用三角恆等式\ \cos(\alpha - \beta)=\cos\alpha\cos\beta + \sin\alpha\sin\beta\,,把正弦-余弦形式中後面的正弦函數跟餘弦函數合併起來:

:A_n\cdot \cos\left(2\pi\tfrac{n}{P}x-\varphi_n\right) \ \equiv \ \underbrace{A_n \cos(\varphi_n)}_{a_n}\cdot \cos\left(2\pi\tfrac{n}{P}x\right) + \underbrace{A_n \sin(\varphi_n)}_{b_n}\cdot \sin\left(2\pi\tfrac{n}{P}x\right)

然後定義振幅A_n \triangleq \sqrt{a_n^2+b_n^2},相位\varphi_n \triangleq \operatorname{arctan2}(b_n,a_n),这里的a_n和b_n对应正弦-余弦形式中A_n和B_n。\tfrac {A_0}{2}是s(x)的 \frac{1}{P}\int_{P}s(x)\,dx。

{{Equation box 1
|indent =:
|title=傅立葉級數,振幅-相位形式
|equation = {{NumBlk||s(x) \sim \frac{A_0}{2} + \sum_{n=1}^\infty A_n\cdot \cos\left(2\pi\tfrac{n}{P}x - \varphi_n \right)|}}
|cellpadding= 6
|border
|border colour = #0073CF
|background colour=#F5FFFA}}

部份求和算子
在描述傅里叶级数行为的时候,经常会为一个函数f(x)介入部份求和算子S_N:
{{Equation box 1
|indent =:
|title=
|equation = {{NumBlk||S_N(f) = \sum_{n=-N}^N c_n e^{\frac{2\pi i n x}{P}}|}}
|cellpadding= 6
|border
|border colour = #0073CF
|background colour=#F5FFFA}}

这里的c_n是f的傅里叶系数。不同于微积分中的级数,傅里叶级数的部份求和必须采用对称形式,否则收敛结果可能不成立。

假設f(x)與g(x)是在\R 上的可積函數,f(x)與g(x)在[-\pi,\pi] 的捲積(f*g)(x)為:

:(f*g)(x)=\int_{- \pi}^{\pi}f(\tau)g(x-\tau)d\tau

周期为2\pi的函數f(x)的傅立葉級數的部份求和,可以经由f(x)与狄利克雷核D_n(x)=\sum_{k=-n}^n e^{ikx}的摺積来表示:
:\begin{align}
S_{N}(f)(x)&=\sum_{n=-N}^{N}c_n e^{inx}\\
&=\sum_{n=-N}^{N} \left( \frac{1}{2 \pi}\int_{-\pi}^{\pi}f(\tau)e^{-in\tau}d\tau \right) \cdot e^{inx} =\frac{1}{2 \pi}\int_{-\pi}^{\pi}f(\tau)\left(\sum_{n=-N}^{N}e^{in(x-\tau)}\right)d\tau \\
&=\frac{1}{2 \pi}(f*D_{N})(x)
\end{align}

收敛性概要
s_N(x)在[x_0,\ x_0+P]近似了s(x),该近似程度会随着N \rightarrow \infty 逐渐改善。这个无穷和s_{\infty}(x)叫做 s的傅里叶级数表示。傅里叶级数的收敛性取决于函数有限数量的极大值和极小值,这就是通常称为傅里叶级数的狄利克雷条件。参见之一。对于广义函数或分布也可以用范数或定义傅里叶系数。在s(x)的不可导点上,如果我们只取无穷级数中的有限项求和,那么在这些点上会有幅度不随N增大而持续变小的起伏,这叫做吉布斯现象,一个简单的例子是方波信号。

在工程应用中,一般假定傅里叶级数除了在不连续点以外处处收敛,原因是工程上遇到的函数比数学家提供的这个假定的反例表现更加良好。特别地,傅里叶级数绝对收敛且一致收敛于s(x),只要在s(x)的导数(或许不会处处存在)是平方可积的。如果一个函数在区间[x_0,x_0+P]上是平方可积的,那么此傅里叶级数在几乎处处的点都收敛于该函数。

File:Fourier series sawtooth wave circles animation.gif|一个相同幅度和频率的锯齿波的近似的可视化
File:Fourier series square wave circles animation.gif|另一个分别采用傅里叶级数的前 1, 2, 3, 4 项近似方波的可视化。(可以在这里看到一个交互式的动画)
File:Example_of_Fourier_Convergence.gif |收敛于某个任意函数的例子。注意其中的吉布斯现象

其他常用表示法
符号c_n在讨论多个不同函数的傅里叶系数时是不够用的。因此习惯上将其替代为函数(这里是函数s)的某种修改形式,即采用函数式符号比如\hat{s}[n]或S[n],来替代下标式符号:
:s(x)= \sum_{n=-\infty}^\infty \hat{s}(n)\cdot e^{2\pi i nx/P} \quad常用的数学符号
:s(x) = \sum_{n=-\infty}^\infty S[n]\cdot e^{i\,2\pi nx/P} \quad常用的工程符号

在工程上,特别是在变量x表示时间的时候,系数序列叫做频域表示。经常使用方括号来强调这个函数的定义域是频率的离散集合。

另一个常用频域表示,使用傅里叶级数系数,调制像梳子一样的:
:S(f) \ \triangleq \ \sum_{n=-\infty}^\infty S[n]\cdot \delta \left(f-\frac{n}{P}\right)

这里的f表示连续频域。在变量x以秒为单位的时候,f以赫兹为单位。采样的间隔为基本频率\tfrac{1}{P}的n倍(即为谐波)。 s_{\infty}(x)可以通过从这种表示恢复出来:

:\begin{align}
\mathcal{F}^{-1}\{S(f)\} &= \int_{-\infty}^\infty \left( \sum_{n=-\infty}^\infty S[n]\cdot \delta \left(f-\frac{n}{P}\right)\right) e^{i 2 \pi f x}\,df \\
&= \sum_{n=-\infty}^\infty S[n]\cdot \int_{-\infty}^\infty \delta\left(f-\frac{n}{P}\right) e^{i 2 \pi f x}\,df \\
&= \sum_{n=-\infty}^\infty S[n]\cdot e^{i\,2\pi nx/P} \ \ \triangleq \ s_\infty(x)
\end{align}

构造出的函数S(f),因而通常称为“傅里叶变换”,即使一个周期函数的傅里叶积分在这个谐波频率上不收敛{{efn-ua|因为周期函数的傅里叶变换的积分定义不是收敛的,需要将周期函数和它的变换视为分布。在这种意义上,\mathcal{F} \{ e^{i \frac{2\pi nx}{P} } \}是一个狄拉克δ函数,它是分布的是例子。}}。

常用的傅里叶级数
下表列出常用的周期函数及其傅里叶级数系数。
*s(x)指示周期P的周期函数。
*A_0、A_n和B_n指示周期函数s(x)的傅里叶级数系数(正弦-余弦形式)。

}
当x= \pi时,傅里叶级数收敛于0,为在x= \pi 处s的左极限和右极限之和的一半。这是傅里叶级数的狄利克雷定理的特例。

这个例子为我们引出了巴塞尔问题的一种解法。

傅里叶诱导
在上例中我们的函数的傅里叶级数展开式看起来不比s(x)=\tfrac{x}{\pi}简单,因此人们需要傅里叶级数的原因也就不会立即显现出来。但还有很多应用,我们举用傅里叶诱导解热方程的例子。考虑边长为\pi米的方形金属版,坐标为(x,y) \in [0,\pi] \times [0,\pi]。如果板内没有热源,并且四个边中三个都保持在0摄氏度,而第四条边y=\pi,对于x \in (0,\pi),保持在温度梯度T(x,\pi) = x摄氏度。在这种情况下,稳态(或者说很长时间过后的)热分布函数T(x,y)不能得出解析解,但却可以证明:
: T(x,y) = 2\sum_{n=1}^\infty \frac{(-1)^{n+1}}{n} \sin(nx) {\sinh(ny) \over \sinh(n\pi)}
这里的\sinh是双曲正弦函数。热方程的这个解是通过将\pi s(x)的傅里叶级数的每一项乘以\tfrac{\sinh(ny)}{\sinh(n\pi)}得到的。尽管示例的函数s(x)的傅里叶级数似乎很复杂,用傅里叶的方法却可以求解这个热分布问题。

其他例子
我們也可以應用傅立葉級數去證明等周不等式,或是構造處處連續而處處不可微的函數。

收斂性
至今还没有判断傅里叶级数的收敛性充分必要条件,但是对于实际问题中出现的函数,有很多种判别条件可用于判断收敛性。比如x(t)的可微性或级数的一致收敛性。在闭区间上满足狄利克雷条件的函数表示成的傅里叶级数都收敛。狄利克雷条件如下:

在定义区间上,x(t)须绝对可积;

在任一有限区间中,x(t)只能取有限个极值点;

在任何有限区间上,x(t)只能有有限个第一类间断点。

满足以上条件的x(t)傅里叶级数都收敛,且:
:1.当t是x(t)的连续点时,级数收敛于x(t);
:2.当t是x(t)的间断点时,级数收敛于\frac{1}{2}[x(t^-)+x(t^+)]。

1966年,里纳特·卡尔松证明了勒贝格二次可积函数的傅立叶级数一定是几乎处处收敛的,即级数在除了一个勒贝格零测集外均收敛。

傅立葉級數收斂证明
假設一個函數在f(x)在[0,2 \pi]上是平方可積,則會有:

:\frac{1}{2 \pi}\int_{0}^{2\pi}|f(x)-S_{N}(f)(x)|^{2}dx \rightarrow 0 當N \rightarrow \infty

證明的第一步:

考慮一系列正交基底,\{e_{n}\}_{n \in \Z},其中e_{n}(x)=e^{-inx},且有

:(e_{n},e_{m})= \begin{cases} 1, & \text{if }n=m \\ 0, & \text{if }n \neq m \end{cases}

然後有(f,e_{n})=\frac{1}{2 \pi} \int_{0}^{2\pi}f(x)e^{-inx} dx= \hat{f}(n)

特別的有,f(x)的傅立葉級數的部分和S_{N}(f)(x)=\sum_{|n|\leq N}\hat{f}(n)e_{n}

然後根據f = f-\sum_{|n|\leq N}\hat{f}(n)e_{n}+\sum_{|n|\leq N}\hat{f}(n)e_{n} 以及畢氏定理,可以有:

:||f ||^{2}= ||f-\sum_{|n|\leq N}\hat{f}(n)e_{n}||^{2}+||\sum_{|n|\leq N}\hat{f}(n)e_{n}||^{2} 替換一下後有 ||f ||^{2}= ||f-S_{N}(f)(x)||^{2}+||\sum_{|n|\leq N}\hat{f}(n)e_{n}||^{2}

如果右邊第一項收斂到0,再根據正交的性質,可以看出上述式子中的右手邊第二項:

:||\sum_{|n|\leq N}\hat{f}(n)e_{n}||^{2}= \sum_{|n|\leq N}| \hat{f}(n)|^{2},這就證明了帕塞瓦尔定理。

证明的第二步:

回到證明右邊第一項,因為函數f(x)可積,找到一個連續函數g(x),然後根據最佳逼近引理,可以找到一個三角多項式p(x),使得

:|f-S_{N}(f)(x)| \leq |f(x)-g(x)|+|g(x)-S_{N}(f)(x)|

故當N \rightarrow \infty ,函數f(x)跟S_{N}(f)(x)的差為0。

其他性质
傅立葉級數的唯一性
如果有一個定義在[-\pi,\pi] 的函數f(x)和g(x),其中函數f(x)和g(x)的傅立葉係數\hat{f}(n)還有\hat{g}(n)相同,且傅立葉級數都收斂到函數本身,那麼可以證明此傅立葉級數具有唯一性,也就是f(x)=g(x)。換句話說,如果函數f(x)在[-\pi,\pi] 上可積,傅立葉係數\hat{f}(n)為0,對所有的n \in \N,那麼函數f(x)=0

卷积定理
给定周期为P的函数s_{_P}和r_{_P},它们具有傅里叶级数系数S[n]和R[n],这里的n \in \mathbb{Z}。

逐点乘积h_{_P}(x) \triangleq s_{_P}(x)\cdot r_{_P}(x),也是周期为P,并且它的傅里叶级数系数是序列S和R的离散卷积:H[n] = (SR)[n]。
周期卷积h_{_P}(x) \triangleq (s_{_P} r)(x) = (s * r_{_P})(x) = \int_{P} s_{_P}(\tau)\cdot r_{_P}(x-\tau)\, d\tau,也是周期为P,它具有傅里叶级数系数:H[n] = P \cdot S[n]\cdot R[n]。
*在c_0(\mathbb{Z})中的双无限序列\left \{c_n \right \}_{n \in Z},是在L^1([0,2\pi])中的傅里叶系数的序列,当且仅当它是在\ell^2(\mathbb{Z})中的两个序列的卷积。

微分性質
我們說f(x)屬於在C^k(\mathbb{T}) \Rightarrow 如果f(x)是一個在實數上以2\pi為週期的函數,且k次可微而且k階連續。

  • 如果f(x)屬於在C^1(\mathbb{T}),那麼f'(x)傅立葉係數\hat{f'}(n)可以被用f(x)傅立葉係數\hat{f}(n)的表示,藉由公式\hat{f'}(n)=in \hat{f}(n)
  • 如果f(x)屬於在C^k(\mathbb{T}),\hat{f^{(k)}}(n)=(in)^k \hat{f}(n)。特別的,當固定k\geq 1,我們有\hat{f^{(k)}}(n)趨近於0當n \rightarrow \infty,且有\hat{f^{(k)}}(n)=O(1/n^k)。

黎曼-勒贝格定理
如果S是可积函数,则\lim_{|n| \to \infty} S[n]=0,\lim_{n \to +\infty} a_n=0而 \lim_{n \to +\infty} b_n=0。

帕塞瓦尔定理
如果函數f(x)屬於在L^2([- \pi, \pi])之中,那麼便有\sum_{- \infty}^{\infty}|\hat{f}(n)|^2= \frac{1}{2 \pi}\int_{- \pi}^{\pi} |f(x)|^2 dx =||f||。

普朗歇爾定理
如果c_0,\, c_{\pm 1},\, c_{\pm 2}, \ldots是系数,并且\sum_{n=-\infty}^\infty |c_n|^2 ,则有一个唯一的函数s\in L^2(P)使得对于所有n有着S[n] = c_n。

延伸
希尔伯特空间的解读
所谓的两个不同向量正交是指它们的内积为0,这也就意味着这两个向量之间没有任何相关性,例如,在三维欧氏空间中,互相垂直的向量之间是正交的。事实上,正交是垂直在数学上的一种抽象化和一般化。一组n个互相正交的向量必然是线性无关的,所以必然可以张成一个n维空间,也就是说,空间中的任何一个向量可以用它们来线性表出。

在希爾伯特空間釋義下,函數的集合{en = einx; nZ}是[−π, π]平方可積函數L2([−π, π])的正交基。這個空間實際上是一個希爾伯特空間,有著針對任何兩個的元素fg的如下內積:
:\langle f,\, g \rangle \;\stackrel{\mathrm{def}}{=} \; \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)\overline{g(x)}\,dx.
三角函数族的正交性用公式表示出来就是:
:\int_{-\pi}^{\pi} \cos(mx)\, \cos(nx)\, dx = \pi \delta_{mn}, \quad m, n \ge 1, \,
:\int_{-\pi}^{\pi} \sin(mx)\, \sin(nx)\, dx = \pi \delta_{mn}, \quad m, n \ge 1
(這裡的δmn是克羅內克函數),而
:\int_{-\pi}^{\pi} \cos(mx)\, \sin(nx)\, dx = 0;\,

参阅

  • 离散时间傅里叶级数
  • 傅里叶变换
  • 维尔斯特拉斯逼近定理

注释
引用
延伸阅读
*

  • 2003 unabridged republication of the 1878 English translation by Alexander Freeman of Fourier's work Théorie Analytique de la Chaleur, originally published in 1822.

*
*
Felix Klein, Development of mathematics in the 19th century. Mathsci Press Brookline, Mass, 1979. Translated by M. Ackerman from Vorlesungen über die Entwicklung der Mathematik im 19 Jahrhundert*, Springer, Berlin, 1928.
*

  • The first edition was published in 1935.
  • 電機電子類科《工程數學》,ISBN 978-957-584-377-9,作者 陳錫冠、曾致煌,高立出版社。

外部链接
*
*
*

评论 (0)

  • 还没有评论,来抢沙发吧。