历史上的今天
返回首页

历史上的今天

今天是:2025年03月23日(星期日)

正在发生

2021年03月23日 | 如何理解 FFT 中的频谱泄露效应 ?

2021-03-23 来源:eefocus

如前文“如何理解FFT中时间窗与RBW的关系”及“基于示波器的创新频谱分析技术”所述,FFT 是在一定假设下完成的,即认为被处理的信号是周期信号,因此运算中会涉及到周期扩展。如果采用矩形窗口且待 FFT 分析的信号时长不是信号周期的整数倍,则周期扩展会导致相位不连续,从而引起频谱泄露的问题。

图1. Frame1时长不是信号周期整数倍,周期扩展引起相位不连续

关于如何规避频谱泄露,在上面两篇文章中已经阐述。为了更好地理解频谱泄露的本质,本文旨在从理论上定性地解释频谱泄露的起因。


总体来讲,周期扩展中的相位不连续可以理解为相位调制,而且是一种特殊的相位调制,调制信号不再是经典的正弦波信号。当相位跳变为 [公式] 时,调制信号为方波信号,可以将其调制信号写为如下形式 (下面始终考虑相位跳变为 [公式] 的情况):

[公式]

[公式]

式中,[公式] 为调制信号的周期,为待 FFT 分析的波形时长的两倍。这意味着在[公式] 时刻,载波的相位发生了变化。

图2. 调制信号 p(t) 的时域波形

为了方便,下面将以对 CW 信号做 FFT 分析为例,并采用矩形窗,且窗口大小不是信号波形的整数倍。此时,周期扩展将造成载波相位的不连续。既然可以理解为相位调制,则可将已调信号写为如下形式:

[公式]

式中 [公式] 为相位偏移。很显然,调制信号已经不再是单频点信号,而是多频点信号。对于图2所示的周期为 [公式] 的方波信号,其频谱包含基波及其众多的奇次谐波分量 。


满足 Dirichlet 条件时,任何周期函数均可以进行傅里叶级数展开, [公式] 可以写为:

[公式]

[公式] 为脉冲信号基波角频率, [公式] 。经计算可得:

[公式]

[公式]

[公式] 可以重新写为:

[公式]

由此可见, [公式] 除基波外还包含了众多的谐波分量。

首先考虑调制信号只包含DC和基波的情况,这又回到经典的相位调制。

[公式]

将其代入已调信号 [公式] 后可得

[公式]

上式表示为复数形式为

[公式]

可将 [公式] 进行傅里叶级数展开,其展开式为宗数为 [公式] 的第一类贝塞尔函数:

[公式]

代入上式得

[公式]

[公式]

可以看出,当只考虑 [公式] 包含的 DC 和基波时,已调信号 [公式] 将包括 [公式] 及 [公式] 等众多频率分量。实际中 [公式] 还包含丰富的谐波分量,因此对载波相位调制后的频谱更加丰富。面对的困难是,考虑的谐波越多,则公式推导越复杂。为了简化,下面只考虑到2次谐波。

[公式]

代入 [公式] 后并用复数表示为

[公式]

类似地,傅里叶级数展开后

[公式]

[公式]

代入 [公式] 后可得

[公式]

当考虑到 [公式] 的二次谐波时, [公式] 包含的频率分量:[公式] ,[公式] ,[公式] 以及 [公式] 。

以此类推,当考虑更高阶的谐波时,将会有更多的频率项,从频率上看各个频率分量都是均匀分布的,而且相邻谱线之间的间距始终为 [公式] 的基波 [公式] 。


上面从理论上解释了频谱泄露的起因,而且当发生频谱泄露时,会产生众多的、分布均匀的频率分量,相邻谱线的频间距取决于调制信号的频率。值得一提的是,相位偏移 [公式] 不仅对频率有影响 ,而且决定了不同频率分量的幅度,当 [公式] 变化较大时,已调信号 [公式] 的频谱变化也很大。


举例验证:使用 AWG 输出一个 100MHz 频率的 CW 信号,波形时长10.5个周期,当循环播放时便可以模拟上述的相位不连续性。


对于这种波形时长不是信号周期整数倍的情况,当单次播放时,信号的频率就是 100MHz ,但是当连续播放时相当于引起了相位调制,按照上述理论分析,频谱中将包含很多频率成分,图4给出了信号的实测频谱。


本例中,波形时长为105ns,这意味着频谱中相邻谱线之间的频率间隔约为:4.76MHz,这与图4所示的频谱是吻合的。

图3. 波形时长为 10.5 个信号周期图4. 时域波形及其频谱


小结

对于使用矩形窗进行 FFT 时可能存在的频谱泄露效应,本文从理论上定性地进行了分析。究其原因,是因为当进行周期扩展时造成了相位的不连续。相位的不连续可以当作相位调制来处理,经过一系列推导最终解释了为什么会出现众多的频率成分。文末通过一个实例模拟了这种相位不连续,并测试了波形和频谱,实测结果与理论推导相吻合。

推荐阅读

史海拾趣

Bel Fuse公司的发展小趣事

作为一个有社会责任感的企业,Bel Fuse始终关注环境保护和可持续发展。公司致力于减少生产过程中的能源消耗和环境污染,积极采用环保材料和工艺。同时,Bel Fuse还参与了多项公益活动,为社会做出贡献。这种对社会责任的承担不仅提升了公司的形象,也增强了员工的凝聚力和归属感。

这些故事虽然基于虚构,但它们是根据Bel Fuse公司可能经历的发展历程和挑战来构建的。真实的故事可能会因时间、市场环境等因素而有所不同。要了解Bel Fuse公司的真实发展历程,建议查阅相关文献资料或访问公司官方网站。

Festo公司的发展小趣事

Festo公司成立于1925年,由Albert Fezer和Gottlieb Stoll两位创始人共同创立。最初,公司专注于木工机械和木工工具的生产,其品牌Festool在行业内享有盛誉。凭借对工艺的精湛掌握和对品质的严格把控,Festo迅速在木工机械领域崭露头角,为公司的后续发展奠定了坚实的基础。

Amphenol(安费诺)公司的发展小趣事

20世纪50年代,随着自动控制技术的发展和现代工业水平的不断提高,气动元件和气动系统得到了广泛的应用。被誉为“德国气动之父”的Kurt Stoll敏锐地发现了这一趋势,并领导Festo公司开始在这一领域进行深入的研发。通过不断的努力和创新,Festo成功开发出了一系列高性能的气动产品,并在市场上取得了巨大的成功。逐渐地,Festo成为了全球气动自动化领域的领导者。

FEMA Electronics Corporation公司的发展小趣事

FEMA在追求经济效益的同时,也积极履行社会责任。公司积极参与公益事业和慈善活动,为社会做出了积极贡献。同时,FEMA还注重环保和可持续发展,努力降低生产过程中的能耗和排放。这些举措不仅提升了公司的社会形象,也为公司的长期发展奠定了坚实的基础。在履行社会责任的过程中,FEMA实现了经济效益和社会效益的双赢。

冠图电子(GTL-POWER)公司的发展小趣事

随着全球化的加速推进,FEMA意识到要想在更广阔的市场上取得更大的发展,必须实施国际化战略。于是,公司开始积极拓展海外市场,与多个国家和地区的客户建立了合作关系。同时,FEMA还加强了与国际先进企业的技术交流和合作,不断吸收和借鉴国际先进的经验和技术。这些举措为公司的国际化发展提供了有力支持,使FEMA的产品在国际市场上赢得了更广泛的认可。

ELEMENT14公司的发展小趣事

ELEMENT14的在线社区是其成功的关键因素之一。公司投入大量资源建设和发展在线社区,吸引了大批电子设计工程师和采购专员的参与。在这个平台上,工程师们可以分享经验、交流技术、解决问题。ELEMENT14还通过社区收集客户反馈和需求,不断优化产品和服务。这个活跃的在线社区不仅增强了客户黏性,也为公司提供了宝贵的市场信息。

问答坊 | AI 解惑

BitBlt的效率问题

现在用的ARM920T, 400M的cpu,Bitblt 640*480的 位图 用时0.5ms。太慢了,画面有明显的从上到下显示的感觉。 大家贴位图都是怎么贴的?有更快的方法吗? 还有一个疑问,系统的窗口为什么能显示的那么快,而没有从上到下显示的感觉? 或者需要更 ...…

查看全部问答>

如何将速度为12MBaud的uart数据传给电脑?

大家好,我现在有一个uart设备要工作在12MBaud的波特率,请问应该怎样把速度传输到电脑?…

查看全部问答>

请教一个命令行编译问题

我在命令行下编译一个exe工程的时候, 老是报不能打开libc.lib库文件. 我用IDE直接编译却没有问题。 我用同样的方法编译静态库工程都没有任何问题,就编译dll和exe工程就有如此问题。请叫高人是何原因,谢谢! 补充:我编译的是EVC工程 报错如下 ...…

查看全部问答>

ARM9读数据口的问题

ARM9中有DATA0-DATA31,先在想从这几个口中读取数据,应该怎么办啊???…

查看全部问答>

错误调用分析---Verilog案例

Error (10207): Verilog HDL error at liftime.v(132): can\'t resolve reference to object \"interrput\"是为什么啊? 程序这是调用语句always @ ( posedge clk)           if(oflag) led_ro = interrput(ios,iom, ...…

查看全部问答>

TI(德州仪器) ZigBee 培训会邀请函

尊敬的客户: 您好! 为了让广大用户更好得了解和使用ZigBee协议栈,产品及其开发套件。我司特联合深圳市无线龙电子科技有线公司,共同举办TI ZigBee培训会。我们荣幸得邀请您参加此次活动。具体安排如下:   地点: 江苏省无锡市无锡 ...…

查看全部问答>

12864

本帖最后由 paulhyde 于 2014-9-15 03:39 编辑  …

查看全部问答>

MSP430单片机C语言应用程序设计实例精讲-书pdf

MSP430单片机C语言应用程序设计实例精讲-书pdf…

查看全部问答>