周期信号的傅里叶变换：从级数系数到冲激串频谱的完整推导

周期信号的傅里叶变换这个话题在信号与系统课程里卡住过不少人。很多教材一上来就把非周期信号的傅里叶变换讲得明明白白到了周期信号这里突然说“周期信号不存在严格意义上的傅里叶变换”然后转头又要你用傅里叶变换去分析周期信号这个弯子绕得人头晕。这篇把这个问题掰开揉碎从数学本质讲到实际做题的套路。1. 为什么周期信号也需要“强行”做傅里叶变换先说个反直觉的结论周期信号在传统定义下根本不存在傅里叶变换但我们在实际工程和考试里依然天天对它做傅里叶变换。这不是数学错了而是我们引入了工具把砖墙凿开了一道门。1.1 傅里叶变换和傅里叶级数的本质区别傅里叶级数处理的是周期信号把信号分解成基波和各次谐波的正弦叠加频率是离散的间隔是 ( 2\pi/T )。傅里叶变换处理的是非周期信号频谱是连续的。两者本来各管各的井水不犯河水。但问题是周期的冲激串、周期矩形波、周期三角波这类信号在工程里随处可见。比如你拿示波器看一个时钟信号就是周期方波电力系统里的正弦交流电也是周期信号。如果只能说“这个信号是周期的我用傅里叶级数去分析”那很多系统分析工具就用不了——比如系统的频率响应、滤波器的设计全部建立在傅里叶变换的框架上。1.2 矛盾焦点不满足绝对可积条件教科书上给傅里叶变换定义了一个存在条件信号必须绝对可积也就是[ \int_{-\infty}^{\infty} |x(t)| , dt \infty ]周期信号永远不满足这个条件因为它从负无穷到正无穷一直在振荡或者重复积分结果是无穷大。所以严格按定义周期信号没有傅里叶变换。但如果就此打住那通信里所有等幅振荡、雷达用的周期脉冲串都不需要频谱分析了这显然说不过去。解决办法就是引入冲激函数——用冲激的“面积”概念来突破绝对可积的限制。这里的本质思想是周期信号的频谱虽然无限大但它集中在离散的频点上每个频点上的“面积”是有限的我们就把这个有限的面积定义为冲激强度用它来刻画周期信号的频谱结构。2. 从傅里叶级数出发推导周期信号的傅里叶变换理解了为什么要做这件事再看怎么推。其实道理不复杂关键就藏在傅里叶级数和傅里叶变换的数学关系里。2.1 用冲激函数表示复指数信号的变换傅里叶变换里有一个重要的对应关系[ e^{j\omega_0 t} \leftrightarrow 2\pi \delta(\omega - \omega_0) ]这个式子的意思是一个无限长的复指数振荡它的频谱只在 ( \omega \omega_0 ) 这一个频率点上有值强度是 ( 2\pi )。注意是 ( 2\pi ) 而不是 1这个系数的由来是因为傅里叶变换定义中积分核里有指数项而在逆变换里有一个 ( 1/2\pi ) 的系数两个约定凑在一起正变换的冲激就必须带一个 ( 2\pi ) 才能保证逆变换能还原信号。很多人第一次看到这个 ( 2\pi ) 会困惑其实可以这样想两个函数做卷积如果一个是 ( x(t) )一个是冲激串时域里就是复制频域里就是乘法。而傅里叶级数系数相当于在频域给一个个离散频点标了权重要把这些权重变成傅里叶变换频谱上的幅值就得让每个频点对应一个面积为 ( 2\pi ) 倍权重的冲激。后面推导会看得更清楚。2.2 把周期信号表达成级数再逐项变换任何一个周期信号 ( x(t) )周期为 ( T )都可以展开成傅里叶级数[ x(t) \sum_{n-\infty}^{\infty} a_n e^{jn\omega_0 t} ]其中 ( \omega_0 2\pi/T )( a_n ) 是傅里叶级数系数。现在对等式两边同时做傅里叶变换[ X(\omega) \mathcal{F}\left{ \sum_{n-\infty}^{\infty} a_n e^{jn\omega_0 t} \right} ]因为傅里叶变换是线性的把求和号提出来逐项变换[ X(\omega) \sum_{n-\infty}^{\infty} a_n \cdot 2\pi \delta(\omega - n\omega_0) ]这样就得出了关键结论周期信号的傅里叶变换是一串离散的冲激每个冲激位于 ( \omega n\omega_0 ) 处强度是 ( 2\pi a_n )。这个式子非常优雅傅里叶级数系数 ( a_n ) 描述的是每个频率分量的相对大小而傅里叶变换把每个频率分量变成一根冲激线线的强度就是去掉系数 ( 2\pi ) 后的结果。也就是说周期信号在频域里是离散的谱线只不过每根谱线窄到面积为 0、高度无穷大只有“强度”是有限的。2.3 系数 ( 2\pi ) 的物理意义和记忆方法很多同学背公式时容易把 ( 2\pi ) 写在奇怪的位置。这里给一个永远忘不掉的记忆锚点做傅里叶逆变换时时域要除以 ( 2\pi )做正变换时频域就只能在每个冲激点上带一个 ( 2\pi ) 保持平衡。本质上这是傅里叶变换对之间保持“能量守恒”的体现。你可以把这个关系类比成单位换算公斤和斤换算时有个系数 2傅里叶级数域和傅里叶变换域之间的换算系数就是 ( 2\pi )。级数系数相当于“斤”变换出来的冲激强度相当于“公斤”乘个 ( 2\pi ) 就能对上账。3. 常见周期信号的傅里叶变换结果光看懂推导还不行得动手算几个典型信号。这部分是考试和工程分析的常客。3.1 周期冲激串周期冲激串定义为[ p(t) \sum_{n-\infty}^{\infty} \delta(t - nT) ]这是一个以 ( T ) 为周期的冲激串每个冲激强度为 1。它的傅里叶级数系数有一个极简的结果[ a_n \frac{1}{T} ]也就是说所有谐波分量的幅度是相等的都是 ( 1/T )。将这个系数代入上面的结论[ P(\omega) \frac{2\pi}{T} \sum_{n-\infty}^{\infty} \delta(\omega - n\omega_0) \omega_0 \sum_{n-\infty}^{\infty} \delta(\omega - n\omega_0) ]这是傅里叶变换里最经典的一个结论周期冲激串的傅里叶变换仍是一个周期冲激串只不过冲激的间隔从时域的 ( T ) 变成了频域的 ( \omega_0 )而强度乘了一个 ( \omega_0 )。这个结果在采样定理里是核心工具。信号经过理想采样就是时域乘以一个周期冲激串频域变成频谱的周期延拓。如果理解的底子够深就能一眼看出采样为什么会造成频谱混叠——因为频域里全是等间隔的冲激副本靠靠挤挤带宽一宽就重叠了。3.2 周期矩形脉冲串设周期为 ( T )脉冲宽度为 ( \tau )幅度为 1 的周期矩形脉冲串。它的傅里叶级数系数[ a_n \frac{\tau}{T} \cdot \frac{\sin(n\omega_0 \tau/2)}{n\omega_0 \tau/2} \frac{\tau}{T} \operatorname{Sa}\left( \frac{n\omega_0 \tau}{2} \right) ]其中 ( \omega_0 2\pi/T )。于是傅里叶变换是[ X(\omega) 2\pi \sum_{n-\infty}^{\infty} \frac{\tau}{T} \operatorname{Sa}\left( \frac{n\omega_0 \tau}{2} \right) \delta(\omega - n\omega_0) ]说人话就是频域里是一根一根离散的谱线谱线的位置是基频的整数倍谱线的强度按照 (\operatorname{Sa}) 函数的包络变化。这个包络就是单个矩形脉冲的傅里叶变换形状——一个 (\sin x / x) 的形态。这里有个值得注意的细节很多初学者以为周期矩形脉冲的频谱就是连续的那个 (\operatorname{Sa}) 函数形状这是错的。连续的 (\operatorname{Sa}) 函数是单个矩形脉冲非周期的频谱而周期的矩形脉冲串频谱是挨个取 (\operatorname{Sa}) 函数在某些点上的值然后在这些点上放冲激。一个是连续曲线一个是离散抽样的谱线两者有血缘关系但不是一回事。3.3 周期三角脉冲和它的记忆套路热词里特别提到了三角脉冲的傅里叶变换记忆方法。单个三角脉冲底宽为 ( 2\tau )高度为 1的傅里叶变换是[ X(\omega) \tau \operatorname{Sa}^2\left( \frac{\omega\tau}{2} \right) ]而周期三角脉冲串的频谱就是把上面这个包络按周期抽样。所以记忆三角脉冲的关键是盯住“平方”和“底宽折半”这两个特征。为什么是平方因为三角脉冲可以看成两个矩形脉冲的卷积时域卷积对应频域相乘矩形脉冲频谱是 (\operatorname{Sa})两个一乘就是平方。这是一个非常漂亮的记忆角度一旦你看穿这一层根本不需要背公式。从卷积的角度来记单个矩形脉冲的频谱( \tau \operatorname{Sa}(\omega\tau/2) )两个宽度相同的矩形脉冲卷积得到三角脉冲时域卷积 频域相乘所以三角脉冲频谱就是 ( \tau \operatorname{Sa}^2(\omega\tau/2) )我自己记的时候还会加上一个检查动作看零点的位置。矩形脉冲频谱的第一个零点在 ( \omega 2\pi/\tau )三角脉冲因为平方了零点位置一样但零点的“平滑度”更高靠近零点时衰减更快。画图的时候三角脉冲频谱的旁瓣比矩形脉冲的更矮更胖一看就知道是平方的效果。3.4 周期信号的三大变换对对比为了帮助记忆整理一个对比表格。这类的对比适合贴在笔记本上反复对照。信号类型时域形式傅里叶级数系数傅里叶变换频谱记忆要点周期冲激串( \sum \delta(t-nT) )( 1/T )( \omega_0 \sum \delta(\omega-n\omega_0) )时域频域都是冲激串周期矩形串矩形脉冲周期延拓( \frac{\tau}{T}\operatorname{Sa}(\frac{n\omega_0\tau}{2}) )离散谱线包络为(\operatorname{Sa})谱线落在基频整数倍上周期三角串三角脉冲周期延拓( \frac{\tau}{T}\operatorname{Sa}^2(\frac{n\omega_0\tau}{2}) )离散谱线包络为(\operatorname{Sa}^2)看成矩形卷积频域平方单矩形脉冲非周期矩形不适用( \tau\operatorname{Sa}(\frac{\omega\tau}{2}) )连续谱零点在 (2\pi/\tau)这个表格集中体现了一个思想周期化和非周期化的关系就是频域的离散化和连续化。时域做周期延拓频域必然从连续曲线变成离散谱线谱线的包络就是原非周期信号的频谱。4. 典型考研题型从频谱反推信号参数结合热词里提到的“大理大学考研真题信号与系统”很多学校的考研题喜欢考一类题给一个周期信号的傅里叶变换一串冲激让你反推时域信号的形式和参数。这类题考察的正是对“冲激强度对应级数系数”这一关系的理解深度。4.1 一道经典的例题拆解假设某信号的傅里叶变换为[ X(\omega) 2\pi \sum_{n-\infty}^{\infty} e^{-|n|} \delta(\omega - 3n) ]第一问判断这个信号是否是周期的如果是求周期。从表达式看冲激出现在 ( \omega 3n ) 的位置( n ) 为整数间隔是 3 rad/s即 ( \omega_0 3 )。这说明时域信号的基波角频率是 3周期 ( T 2\pi/3 )。因为频域是离散等间隔的冲激串时域必然是周期信号。第二问求时域表达式。既然 ( X(\omega) 2\pi \sum a_n \delta(\omega - n\omega_0) )对比系数发现[ a_n e^{-|n|} ]所以时域信号是[ x(t) \sum_{n-\infty}^{\infty} e^{-|n|} e^{j3nt} ]这个结果可以化简成实函数形式。把正负 ( n ) 合并利用欧拉公式 ( e^{j\theta} e^{-j\theta} 2\cos\theta )得到[ x(t) 1 2\sum_{n1}^{\infty} e^{-n}\cos(3nt) ]实际上这就是一个由无数谐波组成的周期信号基频为 3 rad/s。虽然谐波次数趋于无穷但系数按指数级衰减实际能量主要集中在低次谐波上。这个衰减特性就是我们做题时能“截断”求和的理论依据工程上只取前几项就能近似出信号的主要形态。4.2 这类题的通用解题步骤做题时我的习惯步骤是四步走看冲激间隔冲激在频域等间隔排列间距就是 ( \omega_0 )反解出周期 ( T 2\pi/\omega_0 )看系数对应每个冲激前面的系数除以 ( 2\pi )就得到傅里叶级数系数 ( a_n )写级数( x(t) \sum a_n e^{jn\omega_0 t} )化简实函数把正负 n 的项配对合并用欧拉公式转成余弦形式第 4 步是很多人丢分的地方。只要题目给的是实信号级数系数就有共轭对称性 ( a_{-n} a_n^* )正负项合并后一定可以变成实余弦叠加。如果化简不出来基本可以肯定前面系数或者频点位置弄错了。顺带一个避坑建议有些题目给的冲激强度不是 ( 2\pi a_n )而是直接写成 ( \pi a_n ) 之类的形式这时候不要慌除回 ( 2\pi ) 再对比就行。关键是以 F 变换的标准形式为基准不要被出题人牵着鼻子走。5. 三角脉冲与矩形脉冲的频谱对比技巧三角脉冲的记忆问题是考研和期末的热门值得单独拿出来说透。5.1 从时域卷积到频域平方的完整推导设 ( g(t) ) 是高度 1、宽度 ( \tau ) 的矩形脉冲。两个这样的矩形脉冲做卷积[ g(t) * g(t) ]卷积结果是一个高度 ( \tau )、底宽 ( 2\tau ) 的三角脉冲。画图时能看到当两个矩形完全重叠时卷积值最大等于矩形的面积 ( \tau )因为一个矩形在另一个矩形上的积分随着平移从 0 线性增长到 ( \tau ) 再线性减到 0整体呈等腰三角形。由傅里叶变换的卷积定理[ \mathcal{F}{ g(t) * g(t) } G(\omega) \cdot G(\omega) ]而 ( G(\omega) \tau \operatorname{Sa}(\omega\tau/2) )所以[ \mathcal{F}{ 三角脉冲 } \tau^2 \operatorname{Sa}^2(\omega\tau/2) ]如果题目里说的三角脉冲底宽是 ( \tau )而不是 ( 2\tau )那它就是两个宽度 ( \tau/2 ) 的矩形卷积的结果这时候公式里所有 ( \tau ) 都要缩一半。这也是我最早踩过的坑——拿到题不仔细看底宽直接套公式结果整个频谱幅度和零点全错。后来养成习惯先看定义域到底宽度是多少再决定公式里的 ( \tau ) 填什么。5.2 用零点位置快速验证矩形脉冲频谱零点( \omega 2\pi/\tau, 4\pi/\tau, 6\pi/\tau, \cdots )三角脉冲频谱零点因为平方了零点还是在 ( \omega 2\pi/\tau, 4\pi/\tau, \cdots )但是它是二阶零点。意思是频谱在这些点的附近不会穿过横轴产生负值而是轻轻触碰横轴后弹回去。画图时这会让旁瓣全部是正的不会出现矩形脉冲旁瓣那种正负交替的情况。这是判断你算对没算对一个很有用的几何特征。我第一次认真对比这两种函数图形时觉得它们像一对性格迥异的兄弟矩形脉冲频谱旁瓣有正有负性格“尖锐”三角脉冲频谱整个压在零轴上方性格“温和”。之所以旁瓣全正就是因为平方运算把所有负值都变成了正值。5.3 课堂之外的记忆口诀送大家一组自编口诀读起来也不拗口矩形脉冲一阶 Sa零点正负交替画三角脉冲 Sa 平方旁瓣全正零触碰时域卷积频域乘三角就是矩形摞底宽变窄频域胖底宽变宽频域瘦口诀不如理解可靠但当你在考场上有 5 秒钟的暂时性遗忘这个口诀能帮你快速定位到正确的公式方向。6. 周期矩形脉冲频谱随参数变化的行为周期矩形脉冲是最有工程感觉的信号因为它直接联系到数字电路里的时钟、脉冲调制和雷达信号。搞懂它随参数变化的规律比多背十个公式都管用。6.1 脉宽 ( \tau ) 变化对频谱的影响固定周期 ( T )只改变脉宽 ( \tau )。当 ( \tau ) 变小脉冲变窄时(\operatorname{Sa}) 函数的第一个零点右移因为零点在 ( 2\pi/\tau )(\tau) 越小零点越远。于是谱线的包络变得“胖而矮”。时域脉冲越窄频域主瓣越宽这和傅里叶变换的“时宽-带宽积为常数”正好对上。物理上的直觉是脉冲越窄意味着时域波形变化越剧烈需要更多的高频分量去构造这个陡峭的跳变所以频谱自然就展宽了。6.2 周期 ( T ) 变化对频谱的影响固定脉宽 ( \tau )只改变周期 ( T )。( T ) 变大基频 ( \omega_0 2\pi/T ) 变小频域谱线变密。同时每个谱线的强度 ( 2\pi a_n 2\pi(\tau/T)\operatorname{Sa}(\cdot) ) 中 ( \tau/T ) 变小谱线幅度整体降低。这个变化在采样和频谱分析里非常关键周期越大脉冲串越稀疏频域谱线越密集。极端情况是 ( T \to \infty )信号退化成单个脉冲频域谱线密到连成一片就从离散谱过渡到连续谱了。可以说从周期信号到非周期信号的频谱演化本质上就是谱线间隔从有限值压缩到无穷小的过程。我自己在上课的时候喜欢画一个动态图T 从 4τ 慢慢拉到 20τ学生能看见谱线从稀疏变密集、从高耸变低平。这个动态印象比任何公式都刻骨铭心。6.3 工程上对谱线的取舍经验实际做信号分析时我们不可能处理无穷多根谱线。通常取主瓣范围内的谱线作为信号的“有效带宽”。对周期矩形脉冲主瓣宽度为 ( 4\pi/\tau )从负零点到正零点在这个范围内集中了信号的绝大部分能量。如果传输系统带宽太窄只让基波通过波形就会畸变成接近正弦要让波形保持矩形形状至少要让前几根谱线通过。这也是为什么数字电路里对时钟信号的上升沿有要求的原因之一——上升沿越陡等效的高频分量越多对传输线的带宽要求就越高。7. 周期信号傅里叶变换的三大常见误区知识点讲到这里顺手把最常见的三个误区一次性说透。7.1 误区一认为周期信号没有傅里叶变换严格按黎曼积分定义确实没有但引入冲激函数分布意义下的傅里叶变换之后周期信号有了广义变换。关键是不要把“不存在”和“无法分析”画上等号。理解这个东西需要一点“广义函数”的概念但对信号与系统这门课你只需要知道冲激函数让这个问题从“无解”变成了“有解”而且解是唯一的。遇到学得较真的同学问“冲激函数不是普通函数凭什么这么算”我一般会这样解释就像实数轴上没有 ( \sqrt{-1} )但引入虚数单位后整个复分析大厦拔地而起。冲激函数的地位在信号分析里跟虚数在代数里的地位类似——它是一个扩展工具虽然不是普通函数但用它的性质办事结论自洽且能返回工程应用。7.2 误区二混淆级数系数和变换冲激强度这是最容易丢分的地方。看到一个周期信号的傅里叶变换表达式[ X(\omega) \sum_{n-\infty}^{\infty} C_n \delta(\omega - n\omega_0) ]写答案时脑子里一定要有一根弦级数系数 ( a_n C_n / 2\pi )。很多同学的错题本里都留着“忘记除 ( 2\pi ) 导致频谱幅度差 6 倍”的血泪记录。做题时建议先写下 ( X(\omega) 2\pi \sum a_n \delta(\omega - n\omega_0) ) 这个标准形式再跟题目对比就不会错位。7.3 误区三把周期信号频谱的包络当成连续频谱画周期矩形脉冲频谱时谱线是按 (\operatorname{Sa}) 包络排列的离散线不是连续的曲线。不少同学画着画着就描成了一条平滑曲线甚至在记号上把离散谱线和连续谱混为一谈。考试时这种画法通常会被扣分因为它模糊了“周期信号频谱离散”这一核心认识。怎么避免画图的时候先画包络虚线再在虚线对应的频率点上画竖线。只要虚线只起“串起谱线”的作用图形语言就不会骗人。8. 周期信号傅里叶变换的工程应用思路课程内容学到最后还是得回到为什么学它。周期信号的傅里叶变换在三个常见工程领域里是地基。8.1 信号采样与频谱混叠理想采样在时域就是原信号乘以周期冲激串。从频域看这个乘法对应的是原信号频谱与冲激串频谱的卷积——注意这里用到的是周期冲激串傅里叶变换的结论。因为冲激串频谱本身是离散冲激的间隔 ( \omega_s )所以卷积的结果就是原信号频谱以采样角频率 ( \omega_s ) 为周期进行搬移。当采样频率低于奈奎斯特率搬移过去的频谱会和原频谱重叠产生混叠恢复信号就成了不可能完成的任务。整个过程只要理解了周期冲激串的傅里叶变换全部顺理成章。8.2 频闪测速与周期结构分析旋转机械的转速测量常常用频闪仪。频闪灯以某个频率闪烁当频闪频率与旋转周期吻合机械看起来静止不动。这个现象在数学上就是周期信号的傅里叶变换中基波分量的匹配问题。你看到的“静止”本质是采样频率恰好等于信号基频的整数倍。做频谱分析的时候旋转机械的振动信号往往是周期信号叠加随机噪声。用傅里叶变换取出离散谱线对应的基频和倍频就能诊断齿轮是否磨损、转轴是否偏心。这些应用背后的数学工具就是我们刚刚推理的 ( X(\omega) 2\pi \sum a_n \delta(\omega - n\omega_0) )。8.3 通信系统中的载波与调制通信系统中最基本的调幅信号就是把基带信号 ( m(t) ) 乘以高频载波 ( \cos(\omega_c t) )。如果把基带信号看成慢变的周期信号频域里就是原频谱搬到载频 ( \omega_c ) 附近。周期信号的离散谱线在这里扮演了“载波泄漏”的角色工程师要专门设计电路来抑制不需要的谱线。理解了周期信号的傅里叶变换再看频谱仪上那些竖线就不会再觉得神秘了。9. 写在最后的几个学习建议每门课都有它的“任督二脉”信号与系统的任督二脉我觉得就是频谱思维。周期信号的傅里叶变换是连接傅里叶级数和傅里叶变换的一座桥——你在时域里看到的每一个重复波形在频域里都是离散的谱线你在频域里看到的每一根谱线都能还原成时域里对应的谐波振荡。自己学这段内容的时候花过不少时间死磕为什么周期信号的变换结果是一串冲激。后来想通了其实是在跟自己做和解主动接受冲激函数这个抽象工具要知道数学工具的扩展就是为了能描述更广的现象。想通这一点以后后面学采样定理、学拉普拉斯变换、学 Z 变换都没有再被类似的“扩张概念”卡住过。如果你也在备考或者正在信号与系统的苦海里扑腾我的建议是有两条第一画图比背公式重要。把矩形脉冲、三角脉冲、冲激串的时域波形和对应的频域谱线都画一遍画熟了之后很多题目其实不用动笔算看一眼频谱的疏密、包络形状、谱线位置答案就大概有了。第二做题时先写标准形式再代入。所有周期信号的傅里叶变换题先写 ( X(\omega) 2\pi \sum a_n \delta(\omega - n\omega_0) )再根据题目给的信息往里填东西。这一条能帮你避开八成以上的低级错误。最后分享一个小习惯我备课时常把“周期冲激串的变换仍是冲激串”这句话贴在草稿纸上因为它实在是太重要太容易忘又太好用了。建议你也找一个这样的“锚点”把它焊死在记忆里。