傅里叶级数(三角)

                     

贡献者: addis

预备知识 1 复数,定积分,狄拉克符号

   满足狄利克雷条件1的周期函数 f(x),令周期为 2l,可以使用以下三角级数展开

(1)f(x)=a02+n=1ancos(nπlx)+n=1bnsin(nπlx) .
其中
(2)an=1lllf(x)cos(nπlx)dx ,
(3)bn=1lllf(x)sin(nπlx)dx .

   要特别注意的是,严格来说式 1 并不能用等号。因为若函数 f(x) 在某点不连续,在该点处等式右边的级数未必会收敛到 f(x)。对此本文并不过多讨论,且为了方便仍然使用等号。

   式 1 的所有 an 项为偶函数项,所有 bn 项为奇函数项。若 f(x) 是偶函数,所有 bn 项为零,若是奇函数,则所有 an 项为零2。如果 f(x) 不具有奇偶性,可以表示为偶函数和奇函数之和,分别对应所有 an 项和所有 bn

(4)f(x)=12[f(x)+f(x)]+12[f(x)f(x)] .

复值函数

   若上述傅里叶展开公式对实值函数 f(x) 成立,当 f(x) 的函数值变为复数时是否仍然成立呢?答案是肯定的。所以在以下的例子和推导中,我们都只讨论实函数的情况就足够了。

   推导:我们可以把复值函数拆分成实部函数和虚部函数

(5)f(x)=u(x)+iv(x) .
其中 u(x),v(x) 是和 f(x) 相同周期的实函数,且都满足迪利克雷条件。那么,可以把它们分别做傅里叶展开。令 u(x) 的展开系数为 an,bn,而 v(x) 的系数为 an,bn。分别对 u(x),v(x) 列出式 1 式 3 ,且令 an=an+ianbn=bn+ibn,会发现式 1 式 3 同样成立。

   显然,实函数的展开系数 an,bn 都是实数,而复函数的展开系一定不全是实数。

1. 完备性

   傅里叶级数最奇妙的地方大概就是它能展开任意满足狄利克雷条件的函数,这个性质叫做完备性。我们先不证明完备性,只给出几个例子说明随着求和项数增加,三角函数如何逼近不连续或不光滑的函数。

例 1 方波

  

   首先定义一个方波为

(6)f(x)={1(2kπ<x(2k+1)π)1((2k1)π<x2kπ) ,
注意 f(x)x=kπ 处存在间断点。进行傅里叶级数展开得
(7)f(x)=k=04π(2k+1)sin(2k+1)x .
注意由于 f(x) 是奇函数,求和只有正弦项。取级数的前 m 项求和并画图如图 1

图
图 1:有限项傅里叶级数逼近方波

   理论上虽然随着项数变多三角级数会收敛于方波,但是容易发现在函数不连续点的附近,收敛会比较慢,而且会出现高频振荡。这被称为龙格现象,在多项式插值中同样会遇到。

例 2 正弦函数的绝对值

   偶函数 f(x)=|sinx| 存在不光滑的点,展开成傅里叶级数为

(8)f(x)=|sinx|=2π4πk=1+14k21cos2kx ,
m 项和画图如图 2

图
图 2:有限项傅里叶级数逼近 |sin(x)|

习题 1 

   根据傅里叶级数的定义用定积分计算式 7 式 8

2. 系数公式推导

预备知识 2 正交归一基底

   在泰勒级数中,我们通过求 n 阶导数的方式来 “过滤” 出第 n 阶系数。这里我们用积分的方式来 “ 过滤” 系数。我们在推导中还要介绍一个重要的类比来讲解,即把级数中的 sincos 的性质对比几何矢量在正交(不归一)基底上的展开(式 2 )。

   给出一组无穷多个函数

(9)12,sin(πlx),cos(πlx),sin(2πlx),cos(2πlx),sin(nπlx),cos(nπlx) 
其中 n 是正整数。我们把任意满足狄利克雷条件的函数(以下简称任意函数)比作几何矢量,把上面这组函数(式 9 )比作矢量基底(称为函数基底) ,任意函数都可以表示成这组基底的线性组合(即式 1 )。现在把两个任意函数(矢量)f(x)g(x)内积定义为它们的乘积在 [l,l] 内积分
(10)f|g=llf(x)g(x)dx .
可以证明这组基底正交(即任意两个不同的基底内积为 0,证明见下文)但不归一3。与 “几何矢量在正交但不归一的基底上展开” 一样,我们只要把函数分别与各个基底内积,再除以基底的模长平方(模方) 即可获得线性组合的系数 anbn。可以证明所有基底的模方为 l,这样我们就得到了系数公式(式 2 式 3 )。

3. 正弦基底

   若我们只需要在一个区间 [0,l] 上展开函数 f(x) 而不在意其他地方,那么我们可以假想 f(x) 是以 2l 为周期的奇函数,这样,我们只需要用正弦基底展开 f(x) 即可4。于是我们可以说,式 9 中给出的所有正弦基底在区间 [0,l] 具有完备性。

例 3 用正弦基底展开闭区间内的函数

  

图
图 3:三角形函数

   定义区间 [0,a] 内的一个三角形函数如下(图 3 ),先把函数归一化,再用正弦基底做傅里叶展开。

(11)f(x)={2ax(0xa2)22ax(a2<xa) 

   首先假想 f(x) 为奇函数,原则上可以直接使用式 3 计算展开系数,但由于被积函数 f(x)sin(x) 为偶函数,可以先把式 3 化简为

(12)bn=2a0af(x)sin(nπax)dx .

   又注意 f(x) 关于区间中点的对称性,我们可以进一步判断出 n 为偶数时 bn 项都为零,而 n 为奇数项的 bn[0,a/2] 的积分等于在 [a/2,a] 的积分,所以

(13)bn=4a0a/2f(x)sin(nπax)dx=8a20a/2xsin(nπax)dx .
这个积分可以用 Wolfram Alpha 或 Mathematica 完成。结果是
(14)bn=(1)n128π2n2(n 为奇数) .

   进一步说,我们可以用正弦基底在任意区间 [x0,x0+l] 上展开任意函数。要这样做,我们只需要把所有正弦基底平移 x0 即可。

(15)sinπl(xx0),sin2πl(xx0),sinnπl(xx0) 

4. 任意区间的展开

   有时候我们希望对函数 f(x) 在某个任意区间 [a,b] 内做类似式 1 的三角级数,却并不关心 f(x) 是否有周期性,周期是多少,甚至不关心 f(x) 在该区间外是否有定义。那么此时我们将有多种可能的展开方式。

   第一种方法是令 g(x) 是周期为 2l 的周期函数,且在区间 [a,b]g(x)=f(x),然后对 g 傅里叶展开。令式 1 l=(ba)/2,那么容易证明只需要把系数公式式 2 式 3 的积分范围改为 [a,b] 即可(留做习题)。

   第二种方法式 1 l=(ba)/2,把 f(x) 先向左平移 x0=(a+b)/2,也就是使原来 [a,b] 的移动到 [l,l],展开成式 1 的傅里叶级数,再向右平移 x0 即可。这样傅里叶展开变为(a<x<b,下同)

(16)f(x)=a02+n=1ancos[nπl(xx0)]+n=1bnsin[nπl(xx0)] ,
(17)an=1labf(x)cos[nπl(xx0)]dx ,
(18)bn=1labf(x)sin[nπl(xx0)]dx .

   第三种方法l=ba,然后只用正弦或余弦函数来展开(例 3 ):

(19)f(x)=n=1bnsin[nπl(xa)] ,
(20)bn=2labf(x)sin[nπl(xa)]dx .
把两式中的 sin 替换为 cos 也成立,以下不再赘述。从图像上来看,式 19 n=1 的项表示横跨 [a,b] 的一个拱形(sin 的半个周期),n=2 代表两个拱形(sin 的一个周期),等。

   若令式 19 的右边为 g(x),会发现在区间 [al,a+l] 中, g(x) 是关于点 x=a 的奇函数,而在整个实轴上,g(x) 是周期为 2l 的函数。

   为什么这里只需要用 sin 基底就可以实现完备?这是因为若把 g(x) 向左平移 a 后,在区间 [l,l] 上是一个奇函数,而奇函数的傅里叶展开是不需要 cos 项的。如果对比式 19 的右边和式 16 的右边,它们在 [a,b] 内的的部分是相同的,而在 [a,b] 外的部分却不同。

   在实际操作中要如何选择各种方法呢?前两种方法实际上是等价的,把式 16 利用两两角和公式(式 4 )拆开即可证明。另一个主要的考虑因素是例 1 提到的龙格效应。若 f(a)=f(b),那么用前两种方法把 f(x) 拓展为周期函数后,[a,b] 的两个端点不是间断点,这就可以避免此处的龙格效应。特殊地,若 f(a)=f(b)=0,那么用第三种方法中的 sin 展开也可以避免龙格效应,若区间两端的导数 f(a)=f(b)=0,那么当然最好只用 cos 来展开。

   事实上,使用任何 l>(ba)/2 都可以在区间 [a,b] 内展开函数,所以存在无穷多种展开的方法。

5. 证明函数基底正交

   现在证明式 9 中任意两个不同的基底内积等于 0。当 mn 时,有(使用积化和差式 19 式 21

(21)ππsin(mx)sin(nx)dx=12ππcos[(mn)x]dx12ππcos[(m+n)x]dx=0 ,
(22)ππcos(mx)cos(nx)dx=12ππcos[(mn)x]dx+12ππcos[(m+n)x]dx=0 ,
(23)ππsin(mx)cos(nx)dx=12ππsin[(m+n)x]dx12ππsin[(mn)x]dx=0 ,

   其中任意一个函数与自己内积都等于 π(除了常函数 1 积分为 2π)。

(24)ππsin2(nx)dx=π ,
(25)ππcos2(nx)dx=π ,
(26)ππ12dx=2π .

6. 使用归一化的基底

   我们也可以将式 9 的每个基底乘以一个常数进行归一化,即让每个基底的模长为 1。我们将归一化后的正弦项和余弦项分别记为 |Sn|Cnn=1,2,),将常数项记为 |C0

(27)Sn|Sn=ll[Bnsin(nπlx)]2dx=1 ,
(28)Cn|Cn=ll[Ancos(nπlx)]2dx=1 .
解得 An=Bn=1/l 。最后将常数项归一化变为 1/2l 即可。所以正交归一得傅里叶级数基底为
(29)12l,1l sin(πlx),1l cos(πlx),1l sin(nπlx),1l cos(nπlx) 
易证将基底乘以常数不影响它们的正交性,我们可以将正交归一基底的性质用克罗内克 δ 函数(式 6 )表示(常数基底可以记为 |C0
(30)Sm|Sn=δm,n ,Cm|Cn=δm,n ,Sm|Cn=0 .

   如果将任意函数 f(x) 记为 |f,傅里叶级数的公式可以记为(下面默认 |Cn 求和时 n=0,1,,对 |Sn 求和时 n=1,2,

(31)|f=n(an|Cn+bn|Sn) ,
(32)an=f|Cn ,bn=f|Sn .
这比上文的定义要更简洁。我们还可以证明两函数 f(x)g(x) 的内积为
(33)f|g=m(amCm|+bmSm|)n(cn|Cn+dn|Sn)=m,namcnδm,n+bmdnδm,n=n(ancn+bndn) .
该式可以类比几何矢量的坐标计算(式 3 ),即把两个函数在所有基底上的展开系数分别相乘再相加。

   特殊地,一个任意函数的模方为

(34)f|f=n(an2+bn2) .
即将所有展开系数分别平方再相加。这与几何矢量模方公式也相同。这叫做 Parseval 定理


1. ^ 函数值有限,存在有限个间断点和有限个极值点。如无特殊说明,本书中做傅里叶级数的展开的函数都满足狄利克雷条件。
2. ^ 证明:如果 f(x) 是偶函数,那么式 3 中的被积函数就是奇函数,所以在区间 [l,l] 的积分为零。奇函数的证明类似。
3. ^ 归一化的矢量也叫单位矢量,与自身内积等于 1,详见 “点乘
4. ^ 显然我们也可以假想 f(x) 是以 2l 为周期的偶函数,用余弦基底展开 f(x)。选择正弦或余弦还应该考虑到函数在该区间的边界条件:正弦可以保证两端函数值为 0,余弦可以保证区间的导数值为 0。


致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利