1804 年,傅里叶首次提出“在有限区间上由任意图形定义的任意函数都可以表示为单纯的正弦与余弦之和”,但并没有给出严格的证明。1829 年,由德国数学家狄利克雷 (Dirichlet) 证明的狄利克雷条件为傅里叶级数奠定了理论基础。
傅立叶级数的三角形式
周期为 T 的函数 f(x) 可以使用傅立叶级数表示为:
f(x)=2a0+n=1∑∞(ancosT2πnx+bnsinT2πnx)
其中,an 和 bn 是傅立叶系数,可以通过原函数积分计算得到:
an=T2∫0Tf(x)cosT2πnxdx(n=0,1,2,…)
bn=T2∫0Tf(x)sinT2πnxdx(n=0,1,2,…)
周期窗口(积分区间)的选取是任意的,只要长度为一个完整周期即可。
使用 ω=T2π 的写法多用于工科,上述写法多用于数学教材。
傅立叶级数的复指数形式
下面借由欧拉公式 eix=cosx+isinx 将傅里叶级数系数简化成复数指数形式。
f(x)=2a0+n=1∑∞(ancosT2πnx+bnsinT2πnx)
使用 eix=cosx+isinx 将正弦和余弦表示为复指数形式
cosT2πnx=21[eiT2πnx+e−iT2πnx]
sinT2πnx=2i1[eiT2πnx−e−iT2πnx]
代入傅立叶级数得到:
f(x)=2a0+n=1∑∞[an(2eiT2πnx+e−iT2πnx)+bn(2ieiT2πnx−e−iT2πnx)]
重新整理,我们可以将上面的表达式转化为:
f(x)=n=−∞∑∞cneiT2πnx
其中,复数系数 cn 定义为:
cn=T1∫0Tf(x)e−iT2πnxdx
复数系数与原傅立叶系数的关系
对于 n>0:
cn=2an−ibn
对于 n<0:
cn=2a∣n∣+ib∣n∣
对于 n=0:
c0=2a0
这样,原来的傅立叶级数形式就可以简化为复数指数形式。
例题:求以
T 为周期的函数的离散频谱和它的傅里叶级数的复指数形式. $$f_{T} (t) = \left{ \begin{array}{l l} 0, & - T / 2 < t < 0, \ 2, & 0 < t < T / 2 \end{array} \right.$$
解:令 ω0=2π/T ,当 n=0 时,
c0=F(0)=T1∫−T/2T/2fT(t)dt=T1∫0T/22dt=1;
当 n=0 时,
cn=F(nω0)=T1∫−T/2T/2fT(t)e−jnω0tdt=T2∫0T/2e−jnω0tdt=nπj(e−jn2ω0T−1)=nπj(e−jnπ−1)={0,−nπ2j,当n为偶数,当n为奇数.
fT(t) 的傅里叶级数的复指数形式为$$f_{T} (t) = 1 + \sum_{n = - \infty}^{+ \infty} \frac {- 2 j}{(2 n - 1) \pi} \mathrm{e}^{\mathrm{j} (2 n - 1) \omega_{0} t}.$$
振幅谱为$$\mid F (n \omega_{0}) \mid = \left{ \begin{array}{l l} 1, & n = 0, \ 0, & n = \pm 2, \pm 4, \dots , \ \frac {2}{| n | \pi}, & n = \pm 1, \pm 3, \dots . \end{array} \right.$$
相位谱为$$\arg F (n \omega_{0}) = \left{ \begin{array}{l l} 0, & n = 0, \pm 2, \pm 4, \dots , \ \frac {- \pi}{2}, & n = 1, 3, 5, \dots , \ \frac {\pi}{2}, & n = - 1, - 3, - 5, \dots . \end{array} \right.$$
其图形如图所示.
(a)
(b)
(c)
总结
利用欧拉公式可以将傅立叶级数的三角形式转换为复数指数形式。这种形式在许多应用中更为方便,尤其在涉及到对称性或频谱分析时。傅立叶系数 cn 可以通过直接对函数 f(x) 进行复数指数形式的积分来计算,并且与原三角形式中的 an 和 bn 有直接的关系。
奇偶性
有些函数的傅立叶级数只包含正弦或者余弦项
三角级数
1. 周期函数
正弦函数是一种常见而简单的周期函数, 例如描述简谐振动的函数
y=Asin(ωt+ϕ)
实际问题中还有非正弦函数外的周期函数, 如矩形波
周期函数反映了客观世界的周期运动
2. 矩形波展开为正弦函数:
前面介绍过使用函数的幂级数展开函数, 周期函数可以展开为三角函数
f(t)=A0+n=1∑∞Ansin(nωt+ϕn)
- A0称为直流分量
- A1sin(ωt+ϕ)称为一次谐波或基波
- Ansin(ωt+ϕ) 称为n次谐波
傅里叶级数
使用三角公式变形级数项
Ansin(nωt+ϕn)=Ansin(ϕ)cos(nωt)+Ancos(ϕn)sin(nωt)
为了方便, 将上式符号改写为
2a0+n=1∑∞(ancos(nx)+bnsin(nx))
可以得到,其中an bn均为常数。
这样一来我们把周期为ω2π的三角函数变成了周期为2π的三角级数
下面讨论以T=2π为周期的三角级数
三角函数系
例如: 1, cos(x), sin(x), sin(2x), cos(2x)...cos(nx), sin(nx)
三角函数系有正交性: 任意两个三角函数系的乘积在区间[−π,π]上的定积分=0
证明使用积化和差
周期函数展开为三角函数系
如果周期函数可以展开为
f(x)=2a0+n=1∑∞(ancos(nx)+bnsin(nx))
我们不禁要问an bn 之间存在什么关系,换句话说,如何使用已知信息f(x)将an bn表示出来。
为此,我们进一步假设等式两端可以积分。
1. 计算a0
将傅立叶级数的表达式在 [−π,π] 区间上积分:
∫−ππf(x)dx=∫−ππ[2a0+n=1∑∞(ancos(nx)+bnsin(nx))]dx
由于余弦函数和正弦函数在 [−π,π] 内的积分为零,即:
∫−ππcos(nx)dx=0(对于 n=0)
∫−ππsin(nx)dx=0(对于所有 n)
所以,上述积分式简化为:
∫−ππf(x)dx=∫−ππ2a0dx=2a0⋅2π=a0⋅π
因此,
a0=π1∫−ππf(x)dx
2. 计算 an
为了得到 an,将原级数表达式乘以 cos(mx) (其中 m 是任意正整数)后在 [−π,π] 区间上积分:
∫−ππf(x)cos(mx)dx=∫−ππ[2a0cos(mx)+n=1∑∞(ancos(nx)cos(mx)+bnsin(nx)cos(mx))]dx
根据三角函数的正交性:
∫−ππcos(nx)cos(mx)dx={π,0,当 n=m当 n=m
∫−ππsin(nx)cos(mx)dx=0(对于所有 n 和 m)
由此可得:
∫−ππf(x)cos(mx)dx=am⋅π
所以:
am=π1∫−ππf(x)cos(mx)dx(其中 m=1,2,3,…)
3. 计算 bn
类似地,为了得到 bn,将原级数表达式乘以 sin(mx) 后在 [−π,π] 区间上积分:
∫−ππf(x)sin(mx)dx=∫−ππ(2a0sin(mx)+n=1∑∞(ancos(nx)sin(mx)+bnsin(nx)sin(mx)))dx
根据三角函数的正交性:
∫−ππsin(nx)sin(mx)dx={π,0,当 n=m当 n=m
∫−ππcos(nx)sin(mx)dx=0(对于所有 n 和 m)
由此可得:
∫−ππf(x)sin(mx)dx=bm⋅π
所以:
bm=π1∫−ππf(x)sin(mx)dx(其中 m=1,2,3,…)
于是我们得到了傅立叶级数展开中的所有系数:
a0=2π1∫−ππf(x)dx
an=π1∫−ππf(x)cos(nx)dx,(n=1,2,…)
bn=π1∫−ππf(x)sin(nx)dx,(n=1,2,…)
根据假设的前提 T=2π/ω,可以将上式推广至任意周期函数