傅里叶级数

来自testwiki
跳转到导航 跳转到搜索

Template:NoteTA Template:傅里叶变换

方波的傅里叶级数前四项的部份总和。随着增加更多的谐波,这个部份总和收敛成方波。

数学中,傅里叶级数Template:Lang-enTemplate:IPAc-en)是把类似的函数表示成简单谐波的方式。更正式地说,对于满足狄利克雷定理周期函数,其傅里叶级数是由一组正弦余弦函数加权和表示的方法。傅里叶级数与用来找出无周期函数的频率信息的傅里叶变换有密切的关系。

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

歷史

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

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

定義

上图展示了蓝色无周期函数s(x),它只定义在红色区间[0,P]之上。可以认为傅里叶级数,是对下图中最初函数的“周期延拓”的分析,傅里叶级数总是周期函数,即使最初函数s(x)不是周期函数。

傅里叶级数可以用不同的形式来表达,下面将周期P的一个周期函數s(x), x表达为不同形式的傅里叶级数。

正弦-余弦形式

人们常用sin(x)cos(x)三角級數來表示s(x),就是将所有n諧波sin(2πnxP)cos(2πnxP),乘以其各自在s(x)中的權重,求得它们的總和;这些n階諧波的權重稱爲傅立葉級數係數,它们可以藉由如下積分來獲得:

傅立葉級數係數

Template:NumBlk

符号P表示在选定区间上的积分,典型的选择为[P/2,P/2]或者[0,P]。注意A0是函数s(x)Template:LeTemplate:Efn-ua,这个性质扩展到了类似的变换比如傅里叶变换

通过这些系数定义傅里叶级数为:

傅立葉級數,正弦-余弦形式

Template:NumBlk

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

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

  • 波长等于Pn,并且有着同于x的单位。
  • 频率等于nP,并且有着x的倒数单位。

指数形式

下面藉由歐拉公式 eix=cosx+isinx ,将傅里叶级数系数简化成复数指數形式。

根據定義,我們可以得到:

复数傅里叶级数系数

Template:NumBlk

通过将等式Template:EquationNote代入Template:EquationNote,可以证实[4]

复数傅里叶级数系数

cn=1PPs(x)e2πinxPdxfor n

给定复数傅里叶级数系数,可以用公式复原出AnBn

复数傅里叶级数系数

A0=c0An=cn+cnforn>0Bn=i(cncn)forn>0

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

傅立葉級數,指數形式

Template:NumBlk

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

複數值函數

人们習慣將s(x)值域普遍化到複數上,设s(x)是一個複數值函數,它的實部和虛部,都是實數值函數:

s(x)=Re(s(x))+iIm(s(x)),x

定义cncRn+icIn则:

cn=1PPs(x)ei2πpnx dx=1PPRe(s(x))ei2πpnx dx+i1PPIm(s(x))ei2πpnx dx
cRn=1PPRe(s(x))ei2πpnx dx
cIn=1PPIm(s(x))ei2πpnx dx

对于这个复数值函数,它的傅里叶级数的实部,是它的实部的傅里叶级数;它的傅里叶级数的虚部,是它的虚部的傅里叶级数:

s(x)n=(cRn+icIn)ei2πpnx=n=cRnei2πpnx+in=cInei2πpnx

振幅-相位形式

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

Ancos(2πnPxφn)  Ancos(φn)ancos(2πnPx)+Ansin(φn)bnsin(2πnPx)

然後定義振幅Anan2+bn2相位φnarctan2(bn,an),这里的anbn对应正弦-余弦形式中AnBnA02s(x)Template:Le1PPs(x)dx

傅立葉級數,振幅-相位形式

Template:NumBlk

部份求和算子

在描述傅里叶级数行为的时候,经常会为一个函数f(x)介入部份求和算子SN[5]

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

假設f(x)g(x)是在上的可積函數,f(x)g(x)[π,π]捲積(f*g)(x)為:

(f*g)(x)=ππf(τ)g(xτ)dτ

周期为2π的函數f(x)的傅立葉級數的部份求和,可以经由f(x)狄利克雷核Dn(x)=k=nneikx摺積来表示:

SN(f)(x)=n=NNcneinx=n=NN(12πππf(τ)einτdτ)einx=12πππf(τ)(n=NNein(xτ))dτ=12π(f*DN)(x)

收敛性概要

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

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

其他常用表示法

符号cn在讨论多个不同函数的傅里叶系数时是不够用的。因此习惯上将其替代为函数(这里是函数s)的某种修改形式,即采用函数式符号比如s^[n]S[n],来替代下标式符号:

s(x)=n=s^(n)e2πinx/P常用的数学符号
s(x)=n=S[n]ei2πnx/P常用的工程符号

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

另一个常用频域表示,使用傅里叶级数系数,调制像梳子一样的Template:Le

S(f)  n=S[n]δ(fnP)

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

1{S(f)}=(n=S[n]δ(fnP))ei2πfxdf=n=S[n]δ(fnP)ei2πfxdf=n=S[n]ei2πnx/P   s(x)

构造出的函数S(f),因而通常称为“傅里叶变换”,即使一个周期函数的傅里叶积分在这个谐波频率上不收敛Template:Efn-ua

常用的傅里叶级数

下表列出常用的周期函数及其傅里叶级数系数。

  • s(x)指示周期P的周期函数。
  • A0AnBn指示周期函数s(x)的傅里叶级数系数(正弦-余弦形式)。
时域
s(x)
绘图 频域(正弦-余弦形式)
A0Anfor n1Bnfor n1
注释 引用
s(x)=A|sin(2πPx)|for 0x<P
A0=2AπAn={4Aπ1n21n even0n oddBn=0 全波整流正弦 [8]Template:Rp
s(x)={Asin(2πPx)for 0x<P/20for P/2x<P
A0=AπAn={2Aπ1n21n even0n oddBn={A2n=10n>1 半波整流正弦 [8]Template:Rp
s(x)={Afor 0x<DP0for DPx<P
A0=ADAn=Anπsin(2πnD)Bn=2Anπ(sin(πnD))2 0D1
s(x)=AxPfor 0x<P
A0=A2An=0Bn=Anπ 锯齿函数 [8]Template:Rp
s(x)=AAxPfor 0x<P
A0=A2An=0Bn=Anπ 反锯齿函数 [8]Template:Rp
s(x)=4AP2(xP2)2for 0x<P
A0=A3An=4Aπ2n2Bn=0 反全波整流 [8]Template:Rp

基本性質

下表展示在时域中的一些数学运算及其对应的在傅里叶级数系数上的效果。

  • 复数共轭指示为上标星号 * 
  • s(x)r(x)指示周期为P的函数或只定义在x[0,P]中的函数。
  • S[n]R[n]指示sr的傅里叶级数系数(指数形式)。
性质 时域 频域(指数形式) 注释 引用
线性 as(x)+br(x) aS[n]+bR[n] a,b
时间反转/频率反转 s(x) S[n] [9]Template:Rp
时间共轭 s*(x) S*[n] [9]Template:Rp
时间反转且共轭 s*(x) S*[n]
时间实部 Re(s(x)) 12(S[n]+S*[n])
时间虚部 Im(s(x)) 12i(S[n]S*[n])
频率实部 12(s(x)+s*(x)) Re(S[n])
频率虚部 12i(s(x)s*(x)) Im(S[n])
时间移位/频率调制 s(xx0) S[n]ei2πPnx0 x0 [9]Template:Rp
频率移位/时间调制 s(x)ei2πPn0x S[nn0] n0 [9]Template:Rp

对称性质

所有的函数都可以分解成唯一性的Template:Lef(x)=fe(x)+fo(x),这里的fe(x)=f(x)+f(x)2fo(x)=f(x)f(x)2。实数参数的复数值函数f:,对于所有x,如果f(x)=f(x)则称其为“偶对称”,如果f(x)=f(x)则称其为“奇对称”,这里z的上顶横线指示复数共轭

一个复数值函数的实部和虚部,分解成各自的偶部和奇部,就有了四个分量,分别用下标标明为RE、RO、IE和IO。一个复数值时间参数函数的四个分量,与它的复数频率变换的四个分量之间,有着一一映射[10]

时 域s=sRE+sRO+i sIE+i sIO      频 域S=SRE+i SIO+i SIE+SRO

由此可见,各种关系是显而易见的,例如:

範例

一个简单的傅里叶级数

锯齿波周期函数的图
前五个部分傅里叶级数的动态图

我们现在用上面的公式给出一个简单函数的傅里叶级数展开式。考虑一个锯齿波:

s(x)=xπ,forπ<x<π
s(x+2πk)=s(x),for<x< and k

在这种情况下,傅里叶级数为:

An=1πππs(x)cos(nx)dx=0,n0Bn=1πππs(x)sin(nx)dx=2πncos(nπ)+2π2n2sin(nπ)=2(1)n+1πn,n1

可以证明,当s可微时,傅立叶级数在每个点x都收敛于s(x),于是: Template:NumBlkx=π时,傅里叶级数收敛于0,为在x=πs的左极限和右极限之和的一半。这是傅里叶级数的狄利克雷定理的特例。

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

傅里叶诱导

金属板内的热分布,使用傅里叶方法求解

在上例中我们的函数的傅里叶级数展开式看起来不比s(x)=xπ简单,因此人们需要傅里叶级数的原因也就不会立即显现出来。但还有很多应用,我们举用傅里叶诱导解热方程的例子。考虑边长为π米的方形金属版,坐标为(x,y)[0,π]×[0,π]。如果板内没有热源,并且四个边中三个都保持在0摄氏度,而第四条边y=π,对于x(0,π),保持在温度梯度T(x,π)=x摄氏度。在这种情况下,稳态(或者说很长时间过后的)热分布函数T(x,y)不能得出解析解,但却可以证明:

T(x,y)=2n=1(1)n+1nsin(nx)sinh(ny)sinh(nπ)

这里的sinh双曲正弦函数。热方程的这个解是通过将πs(x)的傅里叶级数的每一项乘以sinh(ny)sinh(nπ)得到的。尽管示例的函数s(x)的傅里叶级数似乎很复杂,用傅里叶的方法却可以求解这个热分布问题。

其他例子

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

收斂性

Template:Main

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

  1. 在定义区间上,x(t)绝对可积
  2. 在任一有限区间中,x(t)只能取有限个极值点;
  3. 在任何有限区间上,x(t)只能有有限个第一类间断点

满足以上条件的x(t)傅里叶级数都收敛,且:

1.当tx(t)的连续点时,级数收敛于x(t)
2.当tx(t)的间断点时,级数收敛于12[x(t)+x(t+)]

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

傅立葉級數收斂证明

假設一個函數在f(x)[0,2π]上是平方可積,則會有:

12π02π|f(x)SN(f)(x)|2dx0N

證明的第一步:

考慮一系列正交基底,{en}n,其中en(x)=einx,且有

(en,em)={1,if n=m0,if nm

然後有(f,en)=12π02πf(x)einxdx=f^(n)

特別的有,f(x)的傅立葉級數的部分和SN(f)(x)=|n|Nf^(n)en

然後根據f=f|n|Nf^(n)en+|n|Nf^(n)en 以及畢氏定理,可以有:

||f||2=||f|n|Nf^(n)en||2+|||n|Nf^(n)en||2 替換一下後有 ||f||2=||fSN(f)(x)||2+|||n|Nf^(n)en||2

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

|||n|Nf^(n)en||2=|n|N|f^(n)|2,這就證明了帕塞瓦尔定理

证明的第二步:

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

|fSN(f)(x)||f(x)g(x)|+|g(x)SN(f)(x)|

故當N,函數f(x)SN(f)(x)的差為0。

其他性质

傅立葉級數的唯一性

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

卷积定理

Template:Main

给定周期为P的函数sPrP,它们具有傅里叶级数系数S[n]R[n],这里的n

  • 逐点乘积hP(x)sP(x)rP(x),也是周期为P,并且它的傅里叶级数系数是序列SR离散卷积H[n]=(S*R)[n]
  • 周期卷积hP(x)(sP*r)(x)=(s*rP)(x)=PsP(τ)rP(xτ)dτ,也是周期为P,它具有傅里叶级数系数:H[n]=PS[n]R[n]
  • c0()中的双无限序列{cn}nZ,是在L1([0,2π])中的傅里叶系数的序列,当且仅当它是在2()中的两个序列的卷积[11]

微分性質

我們說f(x)屬於在Ck(𝕋) 如果f(x)是一個在實數上以2π為週期的函數,且k次可微而且k階連續。

  • 如果f(x)屬於在C1(𝕋),那麼f(x)傅立葉係數f^(n)可以被用f(x)傅立葉係數f^(n)的表示,藉由公式f^(n)=inf^(n)
  • 如果f(x)屬於在Ck(𝕋)f(k)^(n)=(in)kf^(n)。特別的,當固定k1,我們有f(k)^(n)趨近於0當n,且有f(k)^(n)=O(1/nk)

黎曼-勒贝格定理

Template:Main 如果S可积函数,则lim|n|S[n]=0limn+an=0limn+bn=0

帕塞瓦尔定理

Template:Main 如果函數f(x)屬於在L2([π,π])之中,那麼便有|f^(n)|2=12πππ|f(x)|2dx=||f||

普朗歇爾定理

Template:Main 如果c0,c±1,c±2,是系数,并且n=|cn|2<,则有一个唯一的函数sL2(P)使得对于所有n有着S[n]=cn

延伸

希尔伯特空间的解读

Template:Main

正弦和餘弦形成了正交集合。正弦、餘弦及其乘積的積分,當mn不同或二函數不同時是0(綠色和紅色區域相等抵消),僅當mn相等並且函數相同時為π。

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

希爾伯特空間釋義下,函數的集合{en = einx; nZ}是[−π, π]平方可積函數L2([−π, π])的正交基。這個空間實際上是一個希爾伯特空間,有著針對任何兩個的元素fg的如下內積:

f,g=def12πππf(x)g(x)dx.

三角函数族的正交性用公式表示出来就是:

ππcos(mx)cos(nx)dx=πδmn,m,n1,
ππsin(mx)sin(nx)dx=πδmn,m,n1

(這裡的δmn克羅內克函數),而

ππcos(mx)sin(nx)dx=0;

参阅

注释

Template:Div col Template:Notelist-ua Template:Div col end

引用

Template:Reflist

延伸阅读

Template:Refbegin

  • Template:Cite book
  • Template:Cite book 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.
  • Template:Cite journal
  • Template:Cite book
  • 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.
  • Template:Cite book
  • Template:Cite book The first edition was published in 1935.
  • 電機電子類科《工程數學》,ISBN 978-957-584-377-9,作者 陳錫冠、曾致煌,高立出版社。

Template:Refend

外部链接

Template:PlanetMath attribution Template:级数 Template:- Template:Authority control