0
  • 聊天消息
  • 系统消息
  • 评论与回复
登录后你可以
  • 下载海量资料
  • 学习在线课程
  • 观看技术视频
  • 写文章/发帖/加入社区
会员中心
创作中心

完善资料让更多小伙伴认识你,还能领取20积分哦,立即完善>

3天内不再提示

数字信号处理 FFT分析信号频谱

信号与系统和数字信号处理 来源:fqj 2019-05-07 14:09 次阅读

离散傅里叶变换的应用之一——用FFT分析信号频谱

题目如下:

数字信号处理 FFT分析信号频谱

欲解此题,关键掌握以下几点:

第一,分清“截取数据长度”(即窗函数长度)与“DFT点数”二者的不同;

第二,能够根据模拟频率推断出DFT谱峰处对应的序号k的数值,方法如下:

首先,由模拟频率转换为数字域频率:

数字信号处理 FFT分析信号频谱

图1

然后,数字域频率对应到DFT的序号k

数字信号处理 FFT分析信号频谱

图2

综合以上两式,得到:

数字信号处理 FFT分析信号频谱

图3

【题目分析与解答】

我们按照DFT分析信号频谱的三个步骤来分别求解:

第一步:采样。所以首先,我们写出采样后离散时间信号的表达式:

数字信号处理 FFT分析信号频谱

其周期为10。

第二步:时域加窗(即截取)

截取10点长,相当于将该周期信号x(n)与10点长的矩形窗相乘,得到v(n),所以我们求v(n)的DTFT。

数字信号处理 FFT分析信号频谱

图4

数字信号处理 FFT分析信号频谱

图5

注意,上图中只画出了[-Π,Π] 区间的图形,实际上DTFT是以2Π为周期的(所以吓人的公式中有西格玛求和符号)。

第三步:频域抽样,也就是对V(e^jw)在 [0,2Π] 区间抽取N个点(N为DFT点数,而非第二步中截取的长度)

再次把狐假虎威的公式摆出来(因为有些同学有强迫症,非要看看公式长什么样,注意,要先把图4中的V(e^jw)的西格玛求和符号去掉,只取0~2Π区间的,所以,后面一项是w-2Π/5+2Π,也就是w+8Π/5,然后再把w变成2Πk/N)(在4月28日发的文章中,下面这个公式写错了,没有人给我提出来。看来你们都没仔细看公式。)

数字信号处理 FFT分析信号频谱

图6

其实V(k) 就是下图中的红点点啦。

数字信号处理 FFT分析信号频谱

图7

好了,那么最后的问题就是,N取不同值时(也就是做不同点数的DFT)这些红点点显然也不同。对于此题来说,这三种N的取值(10、20、128)得到的结果到底是什么呢?

你看到或者听到这里的话,暂停一下,自己算算呗。

我直接把matlab画图的结果给出来。

用前面的公式算一下:f=kfs/N,最后那个128点DFT的图,最高的谱峰序号k是多少?

数字信号处理 FFT分析信号频谱

图8

最后附上matlab程序。

clc;clear all;

f0=1;fs=5;%单位:Hz

n=0:1000;L=10;

xn=cos(2*pi*f0*n/fs);%时域离散时间信号

Xk1=fft(xn(1:L),10);

Xk2=fft(xn(1:L),20);

Xk3=fft(xn(1:L),128);

subplot(311);stem((0:length(Xk1)-1),abs(Xk1));title('10点DFT');

subplot(312);stem((0:length(Xk2)-1),abs(Xk2));title('20点DFT');

subplot(313);stem((0:length(Xk3)-1),abs(Xk3));title('128点DFT');

下面分析一下:

首先看N=10时,此时的结果看似最为干净清爽,只有干干净净两根线。但有的同学要问了,单频信号,只有一个频率成分,应该只有一根谱线呀?为什么会有两根?

我们先看第一根,k=2那个谱线,对应频率为2*fs/10=2*5/10=1Hz,与题设cos(2Πt)完全吻合。k=8那根谱线是怎么一回事呢?是负频率周期延拓过去的,本来在-2,-2+10就等于8了。所以,(敲黑板,以下结论很重要)

对实信号做N点DFT,我们只需要看前N/2根谱线就行了,不用关注N/2~N-1之间的。

再看N=20和N=128的DFT结果,怎么出来那么多根谱线呢?

回过头去看一看,做DFT之前的截取L点长的序列cos(2Πn/5)的频谱到底是什么样子呢?是图7中的虚线所示。而N点DFT,是对V(e^jw)在 [0,2Π] 区间抽取N个点。相当于把连续的频谱图(V(e^jw),如图7中的虚线所示),用一张不透明的纸盖住,纸上以2Π/N为间隔开了一些缝,露出来的点才是我们得到的DFT的结果。这就是频域抽样产生的”栅栏效应“。

所以,不管是10点DFT干干净净的两根线,还是128点DFT密密麻麻的那么多根线,背后隐藏的,都是连续的频谱函数。之所以10点DFT的结果看起来更顺眼,无非是因为因为2Π/N也恰好是旁瓣的宽度(因为信号的周期和截取长度也是10)

在DFT谱分析中,当DFT点数N大于数据本身的实际点数L时,相当于在数据后面补上了L-N个0再做DFT,称为“补零DFT“。补零经常是必要的,补零相当于对信号频谱以更小的间隔采样,得到更多的频谱的信息。而且有时数据长度不是2的整数次幂,如果我们想采用基2FFT算法,就必须进行补零。

此题以单频周期信号为例,展示了不同点数DFT时结果的不同。给人一种错觉:补零,似乎没带来任何好处,反倒是点数最少的10点DFT的结果最好看。

是不是这样呢?当然了,如果你知道这个周期信号的周期是多少,毫无疑问,就截取一个周期的数据,做同样点数的DFT,结果最好看。但在分析实际问题时,要么信号根本没有周期性;要么虽然有周期性,但你不知道;而且实际信号也不会是简单的单频信号,会包含多个频率分量。这个时候,在系统性能、实时性、存储量等等容许的范围内,截取尽量长的数据(即获取更多的信息,得到更高的频率分辨率,下一篇会专门讲频率分辨率),做尽量多点数的DFT(即对频谱进行更为精细的采样)。

而且,实际应用中,由于DFT的点数一般都比较大,我们一般不会以离散的形式画频谱图,而是直接将频谱图化成连续的曲线。例如,上例中,我们截取128点长的数据,做128点DFT,用连续曲线形式画图(matlab中为plot函数),并且只画出前一半(即0~N/2-1)的点,并且把横轴直接转换为Hz,如下图所示:

数字信号处理 FFT分析信号频谱

Matlab代码如下:

数字信号处理 FFT分析信号频谱

上图中,横轴单位为:Hz。采用如下公式,将序号k转换为模拟频率Hz:

数字信号处理 FFT分析信号频谱

DFT分析信号频谱,是实际中应用最广泛的数字信号处理算法,还有很多种题目可以出。还是那句话,题目无穷无尽,而原理就那么多,大家只有掌握了其真正含义,才能以不变应万变。

声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
  • FFT
    FFT
    +关注

    关注

    15

    文章

    434

    浏览量

    59367
  • 离散傅里叶变换

    关注

    0

    文章

    27

    浏览量

    7159

原文标题:数字信号处理系列串讲第12篇(离散信号的频域分析之五)——FFT分析信号频谱题目解答(补充修订版)

文章出处:【微信号:SignalAndSystem_DSP,微信公众号:信号与系统和数字信号处理】欢迎添加关注!文章转载请注明出处。

收藏 人收藏

    评论

    相关推荐

    数字信号处理fft的verilog应用程序

    在学习了数字信号处理之后,在底部使用verilog实现了FFT的应用
    发表于 09-28 06:49

    基于Matlab的系统信号FFT频谱分析与显示pdf

    基于Matlab的系统信号FFT频谱分析与显示
    发表于 06-19 14:55 63次下载

    数字信号处理课件下载

    数字信号处理课件是PPT电子教案,主要内容有离散时间信号和离散时间,离散傅里叶级数及其性质,频率取样,N为合数的FFT算法,快速傅里叶变换(FFT
    发表于 10-30 12:38 0次下载
    <b class='flag-5'>数字信号</b><b class='flag-5'>处理</b>课件下载

    应用FFT信号进行频谱分析

    实验 应用FFT信号进行频谱分析一、 实验目的1、在理论学习的基础上,通过本次实验,加深对快速傅里叶变换的理解,熟悉FFT算法及其程序的编写。2、熟悉
    发表于 05-10 10:00 64次下载

    数字信号处理第4章-快速傅里叶变换(FFT)

    数字信号处理第4章-快速傅里叶变换(FFT)
    发表于 12-28 14:23 0次下载

    数字信号处理技术FFT算法与FPGA的FFT变换设计

    随着集成电路的飞速发展,在图像处理,通信和多媒体等很多领域中,数字信号处理技术已经被广泛应用。快速傅立叶变换(FFT)算法的提出,使得数字信号
    发表于 10-15 10:54 20次下载

    傅里叶变换的应用 FFT分析信号频谱

    数字信号处理中,FFT分析信号频谱是一个重要内容。而频谱分析
    的头像 发表于 05-07 09:49 1w次阅读
    傅里叶变换的应用 <b class='flag-5'>FFT</b><b class='flag-5'>分析</b><b class='flag-5'>信号</b><b class='flag-5'>频谱</b>

    数字信号处理的PPT课件合集免费下载

    本文档的主要内容详细介绍的是数字信号处理的PPT课件合集免费下载包括了:数字信号处理概述,模数转换和数模转换,数字信号
    发表于 11-11 08:00 0次下载

    数字信号处理的资料概述

    本文档的主要内容详细介绍的是数字信号处理的资料概述包括了:1 信号与系统,2 模/数和数/模转换,3 数字信号及其频谱,4
    发表于 11-13 11:49 20次下载
    <b class='flag-5'>数字信号</b><b class='flag-5'>处理</b>的资料概述

    数字信号处理教程之数字信号频谱的详细资料说明

    数字信号处理教程之数字信号频谱的详细资料说明包括了:1 频谱的意义,2 非周期数字信号
    发表于 11-13 17:08 24次下载

    数字信号处理FFT的Verilog工程文件和程序免费下载

    本文档的主要内容详细介绍的是数字信号处理FFT的Verilog工程文件和程序免费下载。
    发表于 11-29 17:13 33次下载
    <b class='flag-5'>数字信号</b><b class='flag-5'>处理</b><b class='flag-5'>FFT</b>的Verilog工程文件和程序免费下载

    基于数字信号处理器实现雷达信号处理系统的设计

    FFT是雷达信号处理的重要工具。DSP内部的硬件乘法器、地址产生器(反转寻址)和多处理内核,保证DSP在相同条件下,完成FFT算法的速度比通
    发表于 12-31 10:26 2738次阅读

    数字信号处理器的常见故障

    数字信号处理器是一种专门用于数字信号处理的微处理器。它能够高效地执行数字信号
    的头像 发表于 05-31 11:53 3138次阅读

    实时频谱分析FFT功能如何采集信号

    实时频谱分析FFT功能如何采集信号? 实时频谱分析仪是一种用于分析信号
    的头像 发表于 01-19 15:01 1187次阅读

    频谱分析仪怎么测量数字信号

    在现代通信、数据传输和信号处理等领域,数字信号的应用日益广泛。数字信号以其高可靠性、抗干扰能力强和易于处理等优点,成为现代通信系统的核心。然
    的头像 发表于 05-17 18:10 991次阅读