fancy-thinking-about-DSP

Last updated on October 16, 2025 am

解谜周期信号的傅里叶变换:从第一页的”天书”公式到第二页的例题实战

引言

在学习信号与系统或数字信号处理时,当我们翻到”周期信号的傅里叶变换”这一节,常常会遭遇一场”认知风暴”。教材的第一页可能会直接抛出一个充斥着冲激函数和无穷求和的复杂公式,让人摸不着头脑。而紧接着的第二页,例题在使用这个公式时,又似乎”违背”了定义,只计算了”一个周期”的变换。

图片1 图片2

你是否也曾感到困惑:

  1. 第一页的公式是怎么来的?为什么长得这么奇怪?
  2. 为什么例题里用了这个公式,算出来的却不是完整的变换?
  3. 那个神秘的、对所有整数 k 的求和到底是什么意思?

本文就将沿着这条思考路径,为你彻底解开这些谜团。

1. 第一页的公式:为何”凭空出现”,不知所云?

我们看到的第一页公式,通常是这样的:

\[ \widetilde{X}(e^{j\omega}) = \sum_{k=-\infty}^{\infty} \frac{2\pi}{N} \widetilde{X}[k] \delta(\omega - \frac{2\pi k}{N}) \]

初次接触的困惑点:

  • 对象不明: 这个公式处理的是无限长的周期信号 \(\tilde{x}[n]\)
  • 定义突兀: 它没有像非周期信号那样从DTFT的求和式推导过来,而是直接”定义”出来的。因为它知道,对于周期信号,直接套用DTFT求和会得到无穷大,根本不收敛1,2
  • 符号诡异: 求和是对所有整数 k 进行的,但 \(\widetilde{X}[k]\) 是离散傅里叶级数(DFS)系数,它本身是周期为N的序列。一个周期序列的系数,为什么要对所有整数求和?信息不是重复的吗?

这个公式更像是一个“宪法总纲”,它从最高原则上定义了周期信号的傅里叶变换应该是什么样子:即在频率轴上的一串冲激函数。但它并没有立即告诉我们这个”宪法”是如何被”立法”出来的。这正是困惑的源头。

2. 第二页的例题:为何”莫名”变成一个周期的变换?

当我们带着困惑看向例题(例如求 \(\cos(0.1\pi n)\) 的傅里叶变换)时,会发现一个更奇怪的现象。

例题的步骤通常是:

  1. 求出信号的DFS系数 \(\widetilde{X}[k]\)
  2. 将DFS系数代入第一页的公式。
  3. 它并没有真的对k从负无穷到正无穷求和,而是只代入了主值区间(如k=-1, 1)的系数,得到了一个结果:

\[ \widetilde{X}(e^{j\omega}) = \pi[\delta(\omega - 0.1\pi) + \delta(\omega + 0.1\pi)] \]

并称此为”一个周期内的傅里叶变换”。

这里的矛盾感达到了顶峰: 第一页的定义明明要求对全体k求和,凭什么这里只算几个k就说是”一个周期”的结果?这岂不是偷换了概念?

3. 第一页公式的推导:破解”双重求和”的等价性谜题

要解决上述矛盾,我们必须回头补上第一页公式的”立法过程”。推导的核心思想是:将周期信号用DFS表示,然后形式地计算其傅里叶变换1,2

步骤详解:

  1. 起点: 周期信号 \(\tilde{x}[n]\) 可以用其DFS系数精确重建: \[ \tilde{x}[n] = \frac{1}{N} \sum_{k=\langle N \rangle} \widetilde{X}[k] e^{j \frac{2\pi}{N} k n} \] 这里的 \(\sum_{k=\langle N \rangle}\) 表示在任意一个长度为N的区间内求和,例如k=0到N-1。这是准确且有限的求和。

  2. 代入与变换: 我们将上述表达式代入(广义)傅里叶变换的框架。我们知道,一个复指数序列 \(e^{j\omega_0 n}\) 的傅里叶变换是 \(2\pi \sum_{l=-\infty}^{\infty} \delta(\omega - \omega_0 - 2\pi l)\)。因此,对于DFS和式中的每一项 \(\frac{1}{N}\widetilde{X}[k] e^{j \frac{2\pi}{N} k n}\),其傅里叶变换是 \(\frac{2\pi}{N}\widetilde{X}[k] \sum_{l=-\infty}^{\infty} \delta(\omega - \frac{2\pi k}{N} - 2\pi l)\)

  3. 得到初始形式: 根据线性性,整个周期信号的傅里叶变换为: \[ \widetilde{X}(e^{j\omega}) = \sum_{k=\langle N \rangle} \frac{2\pi}{N} \widetilde{X}[k] \sum_{l=-\infty}^{\infty} \delta(\omega - \frac{2\pi k}{N} - 2\pi l) \] 这个公式已经是正确的,但它是一个双重求和

  4. 最关键的一步:重新索引(indexing) 现在,我们进行一个巧妙的变量代换,来理解为什么双重求和可以简化为单重求和。

    • \(m = k + lN\),其中 \(m\) 是任意整数。
    • 对于任意一个 \(m\),都存在唯一的一组 \(k\) (在0到N-1之间) 和整数 \(l\),使得 \(m = k + lN\)
    • 由于 \(\widetilde{X}[k]\) 的周期性,\(\widetilde{X}[m] = \widetilde{X}[k]\)
    • 再看冲激函数的参数:\(\omega - \frac{2\pi m}{N} = \omega - \frac{2\pi (k + lN)}{N} = \omega - (\frac{2\pi k}{N} + 2\pi l)\)
    • 因此,双重求和 \(\sum_{k=\langle N \rangle} \sum_{l=-\infty}^{\infty} \delta(\omega - \frac{2\pi k}{N} - 2\pi l)\) 实际上等价于遍历了所有整数 \(m\) 的单重求和 \(\sum_{m=-\infty}^{\infty} \delta(\omega - \frac{2\pi m}{N})\)

    于是,公式被简化为: \[ \widetilde{X}(e^{j\omega}) = \sum_{m=-\infty}^{\infty} \frac{2\pi}{N} \widetilde{X}[m] \delta(\omega - \frac{2\pi m}{N}) \] 将索引 \(m\) 换回 \(k\),就得到了第一页那个”天书”公式。

至此,谜题一解开: 第一页的公式并非凭空出现,它是通过严谨的数学推导,将DFS表示代入傅里叶变换定义后,经过重新索引化简得到的最终形式。那个对全体 \(k\) 的求和,其物理意义是遍历频率轴上所有(无限多个)等间隔的频率点

4. 第二页例题的合理性:为何能只算”一个周期”?

在理解了第一页公式的由来后,第二页例题的做法就变得非常自然和合理了3

  • 频谱的周期性: 离散时间傅里叶变换 \(X(e^{j\omega})\) 本身就是一个以 \(2\pi\) 为周期的周期函数。所以,要描述它,我们只需要写出它在一个周期(通常是 \([-\pi, \pi)\)\([0, 2\pi)\))内的样子就够了,这称为主值区间

  • 例题的聪明之处: 例题知道,对于周期信号,其频谱是离散的冲激。在频率主值区间 \([-\pi, \pi)\) 内,只会出现有限个冲激。例题通过计算DFS系数发现,只有在 \(k=1\)\(k=-1\) (对应 \(k=19\)) 时,系数不为零。这两个 \(k\) 值对应的频率点 \(\omega = \pm 0.1\pi\) 正好落在主值区间 \([-\pi, \pi)\) 内。

  • 因此,例题的 \[ \pi[\delta(\omega - 0.1\pi) + \delta(\omega + 0.1\pi)] \] 正是整个傅里叶变换在主值区间内的表达式。而最后一步加上对 \(k\) 的求和 \(\sum_{k=-\infty}^{+\infty}\),只是将这个主值区间内的频谱模式进行周期延拓,从而得到在整个频率轴上都成立的完整表达式。这两者在数学上是完全等价的。

结论

第一页的公式是”宪法”,定义了周期信号傅里叶变换的完整形式。第二页的例题是”司法解释”,展示了如何高效地运用这部”宪法”:先找到主值区间内的有效信息(一个周期内的变换),再利用周期性宣告其在整个定义域内生效。

希望这篇”破案”笔记,能帮助你彻底理解周期信号傅里叶变换的来龙去脉,让那两页纸的内容不再神秘。


fancy-thinking-about-DSP
http://example.com/2025/10/14/fancy-thinking-about-DSP/
Author
Mingtao Nie
Posted on
October 14, 2025
Licensed under