傅里叶级数是一个专门的工具,它允许任何周期信号(满足某些条件)分解为无限个永恒正弦波的总和。这对许多人来说可能并不明显,但它在数学和图形上都是可以证明的。实际上,这允许傅里叶级数的用户将周期信号理解为各种频率分量的总和。(或者)将信号在特定时间段内的表示形式表示为正交函数线性组合称为傅里叶级数。
矩形级数将信号表示为正弦和余弦项的总和。周期信号可以分解成的正弦波类型完全取决于周期信号的特性。
如果我们有一个周期为 2L 的函数 f(x),我们可以将其分解为正弦和余弦函数的总和,如下所示
![{\displaystyle f(x)={\frac {1}{2}}a_{0}+\sum _{n=1}^{\infty }\left[a_{n}\cos \left({\frac {n\pi x}{L}}\right)+b_{n}\sin \left({\frac {n\pi x}{L}}\right)\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/558f6f57e13337a5d98bcc4582bb2dcd4b897f0d)
系数 a 和 b 可以使用以下积分求出


"n" 是一个整数变量。它可以取正整数(1、2、3 等)。n 的每个值都对应于 A 和 B 的值。幅度为 A 和 B 的正弦波称为谐波。使用傅里叶表示,谐波是信号的原子(不可分割)成分,据说它是正交的。
当我们将 n 设置为 1 时,从上述方程得到的正弦频率值称为基频。给定信号的基频是信号中最强大的正弦分量,也是最需要忠实传输的分量。由于 n 取整数,因此信号的所有其他频率分量都是基频的整数倍。
如果我们考虑时间上的信号,周期T0类似于上述定义中的 2L。那么基频由下式给出

那么基频角频率为

因此,我们可以用更简洁的
替换每个
项。
不同的信号特性会转化为傅里叶级数的特定特性。如果我们能够提前识别这些特性,可以避免不必要的计算。
如果周期信号存在直流偏移,则信号的傅里叶级数将包含一个零频分量,称为直流分量。如果信号不存在直流偏移,则直流分量的幅值为 0。由于傅里叶级数过程的线性性,如果去除直流偏移,我们可以进一步分析信号(例如对称性),并在最后添加回直流偏移。
如果信号为偶函数(关于参考垂直轴对称),则它由余弦波组成。如果信号为奇函数(关于参考垂直轴反对称),则它由正弦波组成。如果信号既不是偶函数也不是奇函数,则它由正弦波和余弦波组成。
如果信号不连续(即存在“跳跃”),则每个谐波n的幅值将按 1/n 的比例下降。
如果信号是连续的,但信号的导数不连续,则每个谐波n的幅值将按 1/n2 的比例下降。
如果信号具有半波对称性,则不存在直流偏移,并且信号由仅位于奇次谐波(1、3、5 等)上的正弦波组成。这很重要,因为具有半波对称性的信号,要传输相同数量的谐波,需要比没有半波对称性的信号两倍的带宽。

如果一个 2L 周期信号具有四分之一波对称性,那么它也必须是半波对称的,所以没有偶次谐波。如果信号是偶函数且具有四分之一波对称性,我们只需要对第一个四分之一周期进行积分。

我们也知道,因为信号是半波对称的,所以没有直流偏移。

因为信号是偶函数,所以没有正弦项。

如果信号是奇函数且具有四分之一波对称性,那么我们可以说
因为信号是奇函数,所以没有余弦项。


由于半波对称性,没有偶次正弦项,并且由于四分之一波对称性,我们只需要对第一个四分之一周期进行积分。

按照惯例,余弦分量的系数用“a”标记,正弦分量的系数用“b”标记。然后可以提到一些重要的结论。
- 如果函数有直流偏移,a0 将不为零。没有 B0 项。
- 如果信号是偶函数,所有 b 项为 0(没有正弦分量)。
- 如果信号是奇函数,所有 a 项为 0(没有余弦分量)。
- 如果函数具有半波对称性,那么所有偶次系数(正弦项和余弦项的系数)都为零,我们只需要对信号进行一半积分。
- 如果函数具有四分之一波对称性,我们只需要对信号的四分之一进行积分。
- 正弦波或余弦波的傅里叶级数只包含一个谐波,因为正弦波或余弦波不能分解成其他正弦波或余弦波。
- 我们可以通过观察信号或信号导数的不连续性来检查级数。如果存在不连续性,谐波将以 1/n 的速度衰减,如果导数不连续,谐波将以 1/n2 的速度衰减。
傅里叶级数也可以用极坐标形式表示,这种形式更紧凑,更易于操作。
如果我们有矩形傅里叶级数的系数 a 和 b,我们可以定义一个系数 x 和一个相位角 φ,它们可以通过以下方式计算:



然后,我们可以使用余弦基函数来根据我们新的傅里叶表示来定义 f(x):

使用余弦基而不是正弦基是一个任意区别,但很重要。如果我们想使用正弦基而不是余弦基,则需要修改上面 φ 的方程。
我们可以明确地证明极坐标余弦基函数等效于具有正弦项和余弦项的“笛卡尔”形式。

根据余弦的倍角公式
![{\displaystyle f(x)=x_{0}+\sum _{n=1}^{\infty }x_{n}\left[\cos(n\omega x)\cos(-\phi _{n})-\sin(n\omega x)\sin(-\phi _{n})\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/af757a1ae198f1d66b3db3dddcf810c92447085b)
根据余弦和正弦的奇偶性质
![{\displaystyle f(x)=x_{0}+\sum _{n=1}^{\infty }x_{n}\left[\cos(n\omega x)\cos(\phi _{n})+\sin(n\omega x)\sin(\phi _{n})\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8df8ef4d9dc97fff848bf1d1aed24bbeb8a2cb7f)
将系数分组
![{\displaystyle f(x)=x_{0}+\sum _{n=1}^{\infty }[x_{n}\cos(\phi _{n})\cos(n\omega x)+x_{n}\sin(\phi _{n})\sin(n\omega x)]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/44055f758ce5411c6d1f4fbf3939a285a95d9834)
这等效于给定以下条件的矩形级数


除以,我们得到


平方并相加,我们得到
![{\displaystyle a_{n}^{2}+b_{n}^{2}=x_{n}^{2}\left[\cos ^{2}(\phi _{n})+sin^{2}(\phi _{n})\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/24d33857ef1a40bafb557f47efbdaeb7bc599ec7)


因此,鉴于上述对 xn 和 φn 的定义,两者是等价的。对于正弦基函数,只需使用正弦倍角公式。其余过程非常相似。
使用欧拉公式,以及一些技巧,我们可以将标准的矩形傅里叶级数转换为指数形式。即使复数更难理解,我们出于多种原因使用这种形式
- 只需要执行一次积分
- 单个指数比正弦波的总和更容易操作
- 它为进一步讨论傅里叶变换提供了逻辑过渡。
我们可以使用欧拉公式从矩形级数构建指数级数

矩形级数由下式给出
![{\displaystyle f(x)=a_{0}+\sum _{n=1}^{\infty }\left[a_{n}\cos(n\omega x)+b_{n}\sin(n\omega x)\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/839af9beab5f76e2f27dd4098f281c1cd6fa73df)
将欧拉公式代入
![{\displaystyle f(x)=a_{0}+\sum _{n=1}^{\infty }\left[{\frac {a_{n}}{2}}e^{in\omega x}+{\frac {a_{n}}{2}}e^{-in\omega x}-{\frac {ib_{n}}{2}}e^{in\omega x}+{\frac {ib_{n}}{2}}e^{-in\omega x}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f565197cfdc5bff016c2f4027c5d987088ef228e)
将表达式分成“正n”和“负n”两部分,得到
![{\displaystyle f(x)=a_{0}+\sum _{n=1}^{\infty }\left[{\frac {a_{n}}{2}}e^{in\omega x}-{\frac {ib_{n}}{2}}e^{in\omega x}\right]+\sum _{n=-\infty }^{-1}\left[{\frac {a_{-n}}{2}}e^{in\omega x}+{\frac {ib_{-n}}{2}}e^{in\omega x}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eeb283e29ae4ebcf1d8d23091ca6803a6a2931ac)

现在将这个表达式合并成一个

我们可以将cn与矩形级数中的an和bn联系起来

这是f(x)的指数傅里叶级数。注意,cn一般来说是复数。还要注意


我们可以直接计算2L周期函数的cn

这可以与使用欧拉公式的矩形形式中的an和bn定义相关联:
.
傅里叶级数的指数形式做了一些与级数的矩形和极坐标形式相比非常有趣的事情:它允许存在负频率分量。为此,指数级数通常被称为“双边傅里叶级数”,因为频谱同时具有正负两侧。当然,这引出了一个问题,“什么是负频率?”
负频率似乎违反直觉,许多人会很快将其视为无稽之谈。然而,进一步研究电气工程(超出了本书的范围)将提供许多负频率在对某些系统建模和理解中起着非常重要作用的示例。虽然一开始可能不太容易理解,但在研究傅里叶域时需要考虑负频率。
负频率遵循对称性的重要规则:对于真实信号,负频率分量始终是正频率分量的镜像。一旦掌握了这个规则,一旦绘制了正侧频谱,绘制负侧频谱就变得微不足道。
然而,在查看双边频谱时,需要考虑负频率的影响。如果负频率是正频率的镜像,并且如果负频率类似于正频率,那么将负分量添加到信号中的效果与将正分量加倍相同。这是计算中指数傅里叶级数系数乘以二分之一的主要原因之一:因为一半的系数位于负频率。
注意:负频率的概念实际上是不切实际的。负频率仅在我们使用傅里叶级数的指数形式时出现在频谱中。为了表示余弦函数,欧拉的关系告诉我们,需要正负指数。为什么?因为要表示一个实函数,如余弦,指数表示法中存在的虚部必须消失。因此,欧拉公式中的负指数使它看起来存在负频率,而实际上并不存在。
理解负频率的另一种方法是使用它们来描述物理世界的数学完整性。假设我们要将头顶正上方旋转的吊扇的旋转情况描述给坐在附近的人。我们会说“它以每分钟 60 转的速度逆时针旋转”。但是,如果我们要将它的旋转描述给从上面观察风扇的人,我们会说“它以每分钟 60 转的速度顺时针旋转”。如果我们习惯上用负号表示顺时针旋转,那么我们会用正号表示逆时针旋转。我们使用正负号来描述相同的过程,具体取决于我们选择的参考系。
带宽是信号传输所需的频率范围的名称,也是特定传输介质的频率容量的名称。例如,如果给定信号的带宽为 10kHz,则它需要带宽至少为 10kHz 的传输介质才能在不衰减的情况下进行传输。
带宽可以用赫兹或弧度每秒来衡量。带宽仅测量正频率分量。所有真实信号都具有负频率分量,但由于它们只是正频率分量的镜像,因此不包括在带宽计算中。
需要注意的是,大多数周期信号是由无限个正弦波的总和组成的,因此需要无限的带宽才能在没有失真地传输。不幸的是,任何可用的通信介质(线路、光纤、无线)都没有无限的带宽可用。这意味着某些谐波将通过介质,而信号的其他谐波将被衰减。
工程学都是权衡。这里的问题是“我需要传输多少个谐波,我可以安全地去除多少个?”使用更少的谐波会导致带宽要求降低,但也导致信号失真增加。这些主题将在未来更详细地考虑。
利用周期和频率之间的关系,我们可以看到一个重要的结论

随着信号周期的减小,基频增加。这意味着每个额外的谐波将被进一步隔开,现在传输相同数量的谐波将需要更多的带宽!一般来说,在考虑周期信号时,必须遵循一条规则:时域中较短的周期需要频域中更多的带宽。在频域中使用较少带宽的信号将需要时域中较长的周期。
让我们考虑基于三次多项式的重复模式

并且f(x)是 2π 周期的
通过观察,我们可以确定傅里叶级数的一些特征
- 该函数是奇函数,因此余弦系数 (an) 将全部为零。
- 该函数没有直流偏移,因此将没有常数项 (a0)。
- 存在间断,因此我们预计会有 1/n 的下降。
因此,我们只需要计算bn项。这些可以通过以下公式找到

将目标函数代入,得

利用分部积分法求解,
![{\displaystyle b_{n}={1 \over \pi }\left({\left[{-x^{3}{{\cos \left({nx}\right)} \over n}}\right]_{-\pi }^{\pi }+\int \limits _{-\pi }^{\pi }{3x^{2}{{\cos \left({nx}\right)} \over n}dx}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/970d6717dfbf1d685ad8a8e93e15112290557d19)
提取公因子
![{\displaystyle b_{n}={1 \over {n\pi }}\left({\left[{-x^{3}\cos \left({nx}\right)}\right]_{-\pi }^{\pi }+3\int \limits _{-\pi }^{\pi }{x^{2}\cos \left({nx}\right)dx}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4f8e74ebf2adbc00c7316b7ecc062a919dafcbd1)
将积分上限和下限代入方括号中,并再次使用分部积分法求解
![{\displaystyle b_{n}={1 \over {n\pi }}\left({-\left({\pi ^{3}\cos \left({n\pi }\right)+\pi ^{3}\cos \left({-n\pi }\right)}\right)+3\left({\left[{x^{2}{{\sin \left({nx}\right)} \over n}}\right]_{-\pi }^{\pi }-\int \limits _{-\pi }^{\pi }{2x{{\sin \left({nx}\right)} \over n}dx}}\right)}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2413f5ba00a6b0982d4a8f582afd97b0880eee98)
因为 cos(x) 是偶函数,所以 cos(-nπ) = cos(nπ)。另外,将积分中的 1/n 提取出来
![{\displaystyle b_{n}={1 \over {n\pi }}\left({-\left({\pi ^{3}\cos \left({n\pi }\right)+\pi ^{3}\cos \left({n\pi }\right)}\right)+{3 \over n}\left({\left[{x^{2}\sin \left({nx}\right)}\right]_{-\pi }^{\pi }-2\int \limits _{-\pi }^{\pi }{x\sin \left({nx}\right)dx}}\right)}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/09884cb9f813005cfc7d492a8e11cc087c6826db)
简化左边部分,并将积分上限和下限代入方括号中,

回想一下,对于整数 n,sin(nπ) 始终等于零。

提出因子并使用分部积分法
![{\displaystyle b_{n}={1 \over {n\pi }}\left({-2\pi ^{3}\cos \left({n\pi }\right)-{6 \over n}\left[{-\left[{x{{\cos \left({nx}\right)} \over n}}\right]_{-\pi }^{\pi }+\int \limits _{-\pi }^{\pi }{{{\cos \left({nx}\right)} \over n}dx}}\right]}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1588db7d8ec2601e8c34270ff8ec650d881fd42f)
![{\displaystyle b_{n}={1 \over {n\pi }}\left({-2\pi ^{3}\cos \left({n\pi }\right)-{6 \over {n^{2}}}\left[-{\left[{x\cos \left({nx}\right)}\right]_{-\pi }^{\pi }+\int \limits _{-\pi }^{\pi }{\cos \left({nx}\right)dx}}\right]}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d8235f94f8de61502b8b90b07844a8d293dff3d5)
求解现在简单的积分,并将极限代入方括号中
![{\displaystyle b_{n}={1 \over {n\pi }}\left({-2\pi ^{3}\cos \left({n\pi }\right)-{6 \over {n^{2}}}\left[{-\left({\pi \cos \left({n\pi }\right)+\pi \cos \left({-n\pi }\right)}\right)+\left[{\sin \left({nx}\right)}\right]_{-\pi }^{\pi }}\right]}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5a0aed3dca748e03ce97d83766ccdf7021bd856a)
由于正弦波一个周期内的面积为零,因此我们可以消除积分。我们再次使用 cos(x) 是偶函数这一事实来简化
![{\displaystyle b_{n}={1 \over {n\pi }}\left({-2\pi ^{3}\cos \left({n\pi }\right)-{6 \over {n^{2}}}\left[{-2\pi \cos \left({n\pi }\right)+0}\right]}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/671904296dbe2732485166f3f51e54c06f3ea31b)
简化




现在,利用cos(nπ)=(-1)n

这是我们最终的bn。我们看到,我们有一个近似1/n的关系(常数“6”随着n的增长变得无关紧要),正如我们预期的那样。现在,我们可以根据以下公式找到傅里叶逼近
![{\displaystyle f\left(x\right)={1 \over 2}a_{0}+\sum \limits _{n=1}^{\infty }{\left[{a_{n}\cos \left({nx}\right)+b_{n}\sin \left({nx}\right)}\right]}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8428da57855f184a157026bcb9c3b57a2a0dbc6e)
由于所有a项均为零,

所以,f(x) = x3的傅里叶级数逼近是
![{\displaystyle f\left(x\right)=\sum \limits _{n=1}^{\infty }{\left[{{{-2\left({-1}\right)^{n}} \over {n^{3}}}\left({\pi ^{2}n^{2}-6}\right)\sin \left({nx}\right)}\right]}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b65e9b06cfbf8eeaaf24c3f4255146db1be03498)
下图显示了前 7 项(红色)和前 15 项(蓝色)的逼近。原始函数以黑色显示。
我们有以下方波信号,作为电压的函数,通过通信介质传播
我们将设置以下值:A = 4 伏特,T = 1 秒。此外,还给定单个脉冲的宽度为T/2。
找到该信号的矩形傅里叶级数。
首先,我们可以清楚地看到该信号确实具有直流值:该信号完全存在于水平轴之上。直流值意味着我们将不得不计算a0项。接下来,我们可以看到,如果我们移除直流分量(将信号向下移动,直到它以水平轴为中心),我们的信号就是一个奇信号。这意味着我们将有bn项,但没有an项。我们还可以看到,该函数具有不连续性和半波对称性。让我们回顾一下
- 直流值(必须计算a0)
- 奇函数(an = 0 对于 n > 0)
- 不连续性(项随着 1/n 衰减)
- 半波对称性(没有偶次谐波)
现在,我们可以按如下方式计算这些值


![{\displaystyle a_{0}={\frac {1}{T}}\left[4t\right]_{0}^{T/2}={\frac {4\times T}{T\times 2}}=2}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8388a5c4d9c99821c9d264dac4b8e1d73b34ea9e)
直观地看,由于信号的占空比为 50%,因此平均值为最大值的一半。
由于函数的奇性,没有余弦项
.
由于半波对称性,只有奇数的正弦项,由以下公式给出:


![{\displaystyle b_{n}=-2\left[{\frac {4}{2n\pi }}\cos \left(2n\pi t\right)\right]_{0}^{1/2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f93b851bd127b01bce4f76bca7721dfd9b774af5)
![{\displaystyle b_{n}=-2\left[{\frac {4}{2n\pi }}\cos \left({\frac {2n\pi }{2}}\right)-{\frac {4}{2n\pi }}\cos \left(2n\pi \times 0\right)\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f39fc22957f9daa15c5d71b989e43fa49c3b7ecb)
![{\displaystyle b_{n}=-{\frac {4}{n\pi }}\left[\cos \left(n\pi \right)-1\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/37a19cc5728594aac73678449dd2ef0be7d54202)
考虑到 cos(nπ)=(-1)n

对于任何偶数 n,这个值等于零,这与我们根据半波对称性做出的预测相符。由于存在不连续性,它也按 1/n 的速率衰减,正如我们预期的那样。
最后,我们可以将傅里叶级数组合在一起,如下所示:


这与以下公式相同:

我们可以看到傅里叶级数与原始函数非常吻合。
维基百科有一篇关于傅里叶级数的文章,但文章的数学推导非常严格。