常见波形的傅里叶级数展开式

引言

近来，在开展课题时遇到了需要将梯形波进行傅里叶级数展开的问题，查询了一些资料（惭愧，一开始就没想着自己动手积分），然后没有找到自己想要的结果（其实有相近的，只不过不是任意周期的，当时没有转变过来），最后还是动手算出来了，在这里做一个小小的记录，算是回顾以前的知识吧，捂脸。

由于像三角波，矩形波，梯形波这种波形不连续，因此在仿真软件中很容易出现计算不收敛的情况。所以，在这种情况下，利用一系列谐波叠加的形式来等价于原来的波形，可以很好的优化模型。

方波拟合

预备知识

公式

给定一个周期为 $T$ 的函数 $x(t)$ ，那么它可以表示为无穷级数：

$f(x) \sim \frac{a_0}{2} + \sum_{n=1}^{\infty} \left [a_n \cos \left (\frac{2 {\pi} nx}{T} \right ) + b_n \sin \left (\frac{2 {\pi} nx}{T} \right ) \right ] = \sum_{n=-\infty}^{\infty} c_n e^{i \frac{2 {\pi} nx}{T} }$

其中傅里叶系数为：

$\left \{ \begin{aligned} a_n = &\frac{2}{T} \int_{t_0}^{t_0+T} f(t) \cdot \cos \left (\frac{2 {\pi} nt}{T} \right )dt \qquad &n=0, 1, 2, \cdots \\[2ex] b_n = &\frac{2}{T} \int_{t_0}^{t_0+T} f(t) \cdot \sin \left (\frac{2 {\pi} nt}{T} \right )dt &n=1, 2,3, \cdots \\[2ex] c_n = &\frac{1}{T} \int_{t_0}^{t_0+T} f(t) \cdot e^{-i \frac{2 {\pi} nt}{T} }dt &n=0, \pm 1, \pm 2, \cdots \end{aligned} \right.$

性质

收敛性

在闭区间上满足狄利克雷条件的函数表示成的傅里叶级数都收敛。狄利克雷条件如下：

在定义区间上， $x(t)$ 需绝对可积；

在任一有限区间中， $x(t)$ 只能取有限个极值点；

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

满足上述条件的 $x(t)$ 傅里叶级数都收敛，且：

当 $t$ 是 $x(t)$ 的连续点时，级数收敛于 $x(t)$

当 $t$ 是 $x(t)$ 的间断点时，级数收敛于 $\frac{1}{2} \left [x(t^-)+x(t^+) \right ]$

正交性

所谓的两个不同向量正交是指它们的内积为0，这也就意味着这两个向量之间没有任何相关性，例如，在三维欧式空间中，互相垂直的向量之间是正交的。三角函数族的正交性用公式表示出来就是：

$\left\{ \begin{aligned} &\int_0^{2 \pi} \cos(mx) \cdot \cos(nx) dx =0 \qquad (m \ne n) \\[2ex] &\int_0^{2 \pi} \sin(nx) \cdot \sin(nx) dx = \pi \\[2ex] &\int_0^{2 \pi} \cos(nx) \cdot \cos(nx) dx = \pi \end{aligned} \right.$

奇偶性

奇函数 $f_o(x)$ 可以表示为正弦级数，而偶函数 $f_e(x)$ 则可以表示成余弦级数：

$\begin{aligned} f_o(x) &\sim \sum_{n=1}^{\infty} b_n \sin \left (\frac{2 {\pi} nx}{T} \right ) \\ f_e(x) &\sim \frac{a_0}{2} + \sum_{n=1}^{\infty} a_n \cos \left (\frac{2 {\pi} nx}{T} \right )\end{aligned}$

几种常见波形的傅里叶级数展开式

梯形波

如上图所示，该梯形波是一个周期为T的奇函数，幅值为 $A_{max}$ ，上升沿时间为 $d$ ，在区间 $\left [0, \frac{T}{2} \right ]$ 的函数表达式为：

$f(t) = \begin{cases} \frac {A_{max} }{d} t, \qquad \qquad & 0 \le t \le d \\[2ex] A_{max}, & d \le t \le \frac{T}{2} - d \\[2ex] \frac {A_{max} }{d} \left (\frac{T}{2} - t \right ), & \frac{T}{2} - d \le t \le \frac{T}{2} \end{cases}$

由奇偶性可知，该波形在区间 $\left [-\frac{T}{2}, \frac{T}{2} \right ]$ 的傅里叶级数展开式为：

$f(t) \sim \sum_{n=1}^{\infty} b_n \sin \left (\frac{2 {\pi} nt}{T} \right )$

其中傅里叶系数为：

$b_n = \frac{2}{T} \int_{-\frac{T}{2} }^{\frac{T}{2} } f(t) \cdot \sin \left (\frac{2 {\pi} nt}{T} \right )dt \quad \quad n=1, 2,3, \cdots$

将 $f(t)$ 函数代入傅里叶系数表达式中，可得：

$\begin{aligned} b_n &= \frac{2}{T} \int_{-\frac{T}{2} }^{\frac{T}{2} } f(t) \cdot \sin \left (\frac{2 {\pi} nt}{T} \right )dt = \frac{4}{T} \int_{0}^{\frac{T}{2} } f(t) \cdot \sin \left (\frac{2 {\pi} nt}{T} \right )dt \\[2ex] &= \frac{4}{T} \left [ \int_0^d \frac{A_{max} }{d} t \cdot \sin \left (\frac{2 {\pi} nt}{T} \right )dt + \int_d^{\frac{T}{2} - d} A_{max} \cdot \sin \left (\frac{2 {\pi} nt}{T} \right )dt + \int_{\frac{T}{2} - d}^{\frac{T}{2} } \frac{A_{max} }{d} \left (\frac{T}{2}-t \right ) \cdot \sin \left (\frac{2 {\pi} nt}{T} \right )dt \right ] \\[2ex] &= \left . \left . \frac{4}{T} \left [ -\frac{A_{max} }{d} {T \over {2 \pi n} } \cdot t \cdot \cos \left ( \frac{2\pi nt}{T} \right ) \right | _0^d + {A_{max} \over d } {T^2 \over 4 \pi^2 n^2} \cdot \sin \left( {2 \pi n t \over T} \right ) \right | _0^d \right ] + {4 \over T} \left [ \left . - {A_{max}T \over 2 \pi n} \cdot \cos \left( {2\pi nt \over T}\right) \right |_d^{ {T \over 2}-d} \right ] + {} \\[2ex] &{} + \left. {4 \over T} \left [ -{A_{max} \over d} {T \over 2 \pi n} \cdot \left ( {T \over 2} - t \right ) \cdot \cos \left ( \frac{2\pi nt}{T} \right ) \right | _{ {T \over 2} - d}^{T \over 2} - \left . {A_{max} \over d}{T^2 \over 4 \pi^2 n^2} \cdot \sin \left( {2 \pi n t \over T} \right ) \right | _{ {T \over 2} - d}^{T \over 2} \right ] \\[2ex] &= {4 \over T} \left \{ {A_{max} T^2 \over 4d \pi ^2 n^2} \left [ \sin \left( {2 \pi n \over T} \cdot d \right) + \sin \left( {2 \pi n \over T} \cdot \left( {T \over 2} - d \right ) \right) \right ] \right \} ={A_{max} T \over d \pi ^2 n^2} \left[ \sin\left( {2 \pi n d \over T} \right) + \sin \left( n\pi - {2 \pi n d \over T} \right) \right] \end{aligned}$

由
$\sin \left( n\pi - {2 \pi n d \over T} \right) = \begin{cases} \sin\left( {2 \pi n d \over T} \right) \qquad \qquad &n=2N-1 \\[2ex] -\sin\left( {2 \pi n d \over T} \right) &n=2N \end{cases}$
可得：

$b_n = \begin{cases} {2 A_{max} T \over d \pi ^2 n^2} \left[ \sin\left( {2 \pi n d \over T} \right) \right ] \qquad \qquad &n=2N-1 \\[2ex] 0 &n=2N \end{cases}$

综上所述，可以得到该梯形波在区间 $\left[ -{T \over 2}, {T \over 2} \right ]$ 的傅里叶级数展开式为：

$f(t)\sim\frac{4A_{max}}{\pi\omega d}\sum_{n=1}^{\infty}{\sin ( { (2n-1) \omega d }) \over (2n-1)^2} \cdot \sin((2n-1) \omega t) \qquad n=1,2,3,\cdots$
其中： $\omega = {2 \pi \over T}$

脉冲波（偶函数）

脉冲波

如上图所示，该脉冲波是一个周期为T的偶函数，幅值为 $A_{max}$ ，脉冲宽度为 $\alpha T$ ，在区间 $\left[ -{T \over 2}, {T \over 2} \right ]$ 的函数表达式为：

$f(t)= \begin{cases} A_{max}, \quad &|t| \le { {\alpha T} \over 2} \\[2ex] 0, &|t| \gt { { {\alpha T} } \over 2} \end{cases}, \qquad - {T \over 2} \le t \le {T \over 2}$

由奇偶性可知，该波形在区间 $\left[ -{T \over 2}, {T \over 2} \right ]$ 的傅里叶级数展开式为：

$f(t) \sim \frac{a_0}{2} + \sum_{n=1}^{\infty} a_n \cos \left (\frac{2 {\pi} n t}{T} \right )$

其中傅里叶系数为：

$a_n =\left\{\begin{aligned} &{2 \over T} \int_{-{T \over 2} }^{T \over 2} f(t) dt \qquad \qquad \qquad \qquad & n=0 \\[2ex] &\frac{2}{T} \int_{-\frac{T}{2} }^{\frac{T}{2} } f(t) \cdot \cos \left (\frac{2 {\pi} nt}{T} \right )dt & n=1, 2,3, \cdots \end{aligned}\right.$

将 $f(t)$ 函数代入傅里叶系数表达式中，可得：

$\begin{aligned} a_0 &= {2 \alpha A_{max} } \\[2ex] a_n &={2 \over T} \int_{-{\alpha T \over 2} }^{\alpha T \over 2} A_{max} \cos \left (\frac{2 {\pi} nt}{T} \right )dt = {4 \over T} \int_{0}^{\alpha T \over 2} A_{max} \cos \left (\frac{2 {\pi} nt}{T} \right )dt\\[2ex] &= \left . {4 \over T} {A_{max} T \over 2 \pi n } \sin \left( {2 \pi n t \over T} \right) \right|_0^{\alpha T \over 2} ={2A_{max} \over n \pi} \sin (\alpha n \pi) \qquad \qquad \qquad \qquad n=1,2,3,\cdots \end{aligned}$

因此，可以得到该梯形波在区间 $\left [ -{T \over 2}, {T \over 2} \right ]$ 的傅里叶级数展开式为：

$f(t) \sim \alpha A_{max} + {2A_{max} \over \pi} \sum _{n=1}^{\infty}{\sin (\alpha n \pi) \over n} \cos(n\omega t) \qquad n=1,2,3, \cdots$

其中： $\omega = {2 \pi \over T}$

方波（奇函数）

方波

同理，该方波在区间 $\left [ -{T \over 2}, {T \over 2} \right ]$ 的傅里叶级数展开式为：

$f(t) \sim {4A_{max} \over \pi} \sum_{n=1}^{\infty}{\sin((2n-1)\omega t) \over 2n-1} \qquad \qquad n=1,2,3,\cdots$

其中： $\omega = {2 \pi \over T}$

三角波（奇函数）

三角波

同理，该三角波在区间 $\left [ -{T \over 2}, {T \over 2} \right ]$ 的傅里叶级数展开式为：

$f(t) \sim {8A_{max} \over \pi^2}\sum_{n=1}^{\infty}{(-1)^{n-1}\sin((2n-1)\omega t) \over (2n-1)^2} \qquad n=1,2,3,\cdots$

锯齿波（非奇非偶函数）

锯齿波

该锯齿波如上图所示，在区间 $[0, T]$ 的函数表达式为：

$f(t)={A_{max} \over T}t \qquad \qquad 0 \le t \le T$

由于该函数为非奇非偶函数，因此，该波形在区间 $[0, T]$ 的傅里叶级数展开式为：

$f(x) \sim \frac{a_0}{2} + \sum_{n=1}^{\infty} \left [a_n \cos \left (\frac{2 {\pi} nx}{T} \right ) + b_n \sin \left (\frac{2 {\pi} nx}{T} \right ) \right ]$

其中傅里叶系数为：

$\begin{aligned} a_n &= \frac{2}{T} \int_{0}^{T} f(t) \cdot \cos \left (\frac{2 {\pi} nt}{T} \right )dt \qquad \qquad & n=0, 1, 2, \cdots \\[2ex] b_n &= \frac{2}{T} \int_{0}^{T} f(t) \cdot \sin \left (\frac{2 {\pi} nt}{T} \right )dt & n=1, 2,3, \cdots \end{aligned}$

将 $f(t)$ 函数代入傅里叶系数表达式中，可得：

$\begin{aligned} a_0 &=A_{max} \\[2ex] a_n &= {2 \over T} \int_0^T{A_{max} \over T}t \cdot \cos \left( {2 \pi nt\over T} \right) dt \\[2ex] & = {2 \over T} {A_{max} \over T}{T \over 2 \pi n} \left . \left . \left[ t\sin \left({2\pi nt \over T}\right) \right|_0^T + {T \over 2 \pi n} \cos \left({2\pi nt \over T}\right) \right|_0^T \right] =0 \qquad &n=1,2,3,\cdots\\[2ex] b_n &= {2 \over T} \int_0^T{A_{max} \over T}t \cdot \sin \left( {2 \pi nt\over T} \right) dt \\[2ex] &= -{2 \over T} {A_{max} \over T}{T \over 2 \pi n} \left . \left . \left[ t\cos \left({2\pi nt \over T}\right) \right|_0^T - {T \over 2 \pi n}\sin \left({2\pi nt \over T}\right) \right|_0^T \right] =-{A_{max} \over n\pi} &n=1,2,3,\cdots \end{aligned}$

因此，可以得到该锯齿波在区间 $[0,T]$ 的傅里叶级数展开式为：

$f(t) \sim {A_{max} \over 2}-{A_{max} \over \pi} \sum_{n=1}^{\infty} {\sin(n \omega t) \over n} \qquad \qquad n=1,2,3,\cdots$

结语

这里仅仅列出了极小部分的波形的傅里叶级数展开式，对于其它波形，类似代入计算即可，给出公式之后，更多的是考验数学积分计算了。

参考文献

维基百科编者. 傅里叶级数
百度百科编者. 傅里叶级数
Fourier Series Examples

最后编辑于：2020.05.05 15:35:41

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 213,616评论 6赞 492
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 91,020评论 3赞 387
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 159,078评论 0赞 349
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 57,040评论 1赞 285
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 66,154评论 6赞 385
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 50,265评论 1赞 292
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 39,298评论 3赞 412
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 38,072评论 0赞 268
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 44,491评论 1赞 306
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 36,795评论 2赞 328
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 38,970评论 1赞 341
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 34,654评论 4赞 337
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 40,272评论 3赞 318
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 30,985评论 0赞 21
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 32,223评论 1赞 267
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 46,815评论 2赞 365
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 43,852评论 2赞 351

常见波形的傅里叶级数展开式

引言

预备知识

公式

性质

几种常见波形的傅里叶级数展开式

脉冲波（偶函数）

方波（奇函数）

三角波（奇函数）

锯齿波（非奇非偶函数）

结语

参考文献

推荐阅读更多精彩内容