当前位置:首页 > 实用技巧 >

窗函数和频率取样法的区别是什么(窗函数频谱为什么是无限的)

来源:原点资讯(www.yd166.com)时间:2023-11-18 11:02:34作者:YD166手机阅读>>

关于频率分辨率的2种解释

解释一:

频率分辨率可以理解为在使用DFT时,在频率轴上的所能得到的最小频率间隔

窗函数和频率取样法的区别是什么,窗函数频谱为什么是无限的(1)

其中N为采样点数,fs为采样频率,Ts为采样间隔。所以NTs就是采样前模拟信号的时间长度T,信号长度越长,频率分辨率越好。

是不是采样点数越多,频率分辨力提高了呢?其实不是的,因为一段数据拿来就确定了时间T,注意:f0=1/T,而T=NTs,增加N必然减小Ts,因此,增加Nf0是不变的。只有增加点数的同时导致增加了数据长度T才能使分辨率越好。

还有容易搞混的一点,我们在做DFT时,常常在有效数据后面补零达到对频谱做某种改善的目的,我们常常认为这是增加了N,从而使频率分辨率变好了,其实不是这样的,补零并没有增加有效数据的长度,仍然为T

但是补零其实有其他好处:

  • 使数据N为2的整次幂,便于使用FFT。
  • 补零后,其实是对DFT结果做了插值,克服“栅栏”效应,使谱外观平滑化。我把“栅栏”效应形象理解为,就像站在栅栏旁边透过栅栏看外面风景,肯定有被栅栏挡住比较多风景,此时就可能漏掉较大频域分量,但是补零以后,相当于你站远了,改变了栅栏密度,风景就看的越来越清楚了。
  • 由于对时域数据的截短必然造成频谱泄露,因此在频谱中可能出现难以辨认的谱峰,补零在一定程度上能消除这种现象。

那么选择DFT时N参数要注意:

  • 由采样定理:fs>=2fh,
  • 频率分辨率:f0=fs/N,所以一般情况给定了fhf0时也就限制了N范围:N>=fs/f0。

解释二:

频率分辨率也可以理解为某一个算法(比如功率谱估计方法),将原信号中的两个靠得很近的谱峰依然能保持分开的能力。这是用来比较和检验不同算法性能好坏的指标。在信号系统中我们知道,宽度为N的矩形脉冲,它的频域图形为sinc函数,两个一阶零点之间的宽度为4π/N。由于时域信号的截短相当于时域信号乘了一个矩形窗函数,那么该信号的频域就等同卷积了一个sinc函数,也就是频域受到sinc函数的调制了,根据卷积的性质,因此两个信号圆周频率之差W0必须大于4π/N。从这里可以知道,如果增加数据点数N,即增加数据长度,也可以使频率分辨率变好,这一点与第一种解释是一样的。同时,考虑到窗函数截短数据的影响存在,当然窗函数的特性也要考虑,在频率做卷积,如果窗函数的频谱是个冲击函数最好了,那不就是相当于没截断吗?可是那不可能的。

我们考虑窗函数主要是以下几点:

  • 主瓣宽度B最小(相当于矩形窗时的4π/N,频域两个过零点间的宽度)。
  • 最大边瓣峰值A最小(这样旁瓣泄露小,一些高频分量损失少了)。
  • 边瓣谱峰渐近衰减速度D最大(同样是减少旁瓣泄露)。

在此,总结几种很常用的窗函数的优缺点:

  • 矩形窗:B=4π/N A=-13dB D=-6dB/oct
  • 三角窗:B=8π/N A=-27dB D=-12dB/oct
  • 汉宁窗:B=8π/N A=-32dB D=-18dB/oct
  • 海明窗:B=8π/N A=-43dB D=-6dB/oct
  • 布莱克曼窗:B=12π/N A=-58dB D=-18dB/oct

可以看出,矩形窗有最窄的主瓣,但是旁瓣泄露严重。汉宁窗和海明窗虽主瓣较宽,但是旁瓣泄露少,是常选用的窗函数。

采样周期与频率分辨率

fs/N常称作为频率分辨率,它实际是作FFT时谱图中的两条相邻谱线之间的频率间隔,也有称作步长。单位是Hz、Khz等。频率分辨率实际有二重含义,在这里只是其中一种。

1/fs的单位的s、ms、us或分、时...年等。1/fs代表采样周期,是时间域上两个相邻离散数据之间的时间差。因此fs/N用在频率域,只在DFT以后的谱图中使用;而1/fs用时间域,只要数据经采样,离散化后任何其它的应用中都可使用。例如有的数字滤波器中就用到。

窗函数和频率取样法的区别是什么,窗函数频谱为什么是无限的(2)

其中Δf是频率采样间隔,同时也是频率分辨率的重要指标,如果这个值越小,则频率分辨率越高。

1/fs往往用在求时间序列上,如 (0:N-1)*1/fs等等,如果这个不好理解,可以把前面的公式求倒数,这就清楚多了。

采样定理

采样过程所应遵循的规律,又称取样定理、抽样定理。采样定理说明采样频率与信号频谱之间的关系,是连续信号离散化的基本依据。采样定理是1928年由美国电信工程师H.奈奎斯特首先提出来的,因此称为奈奎斯特采样定理。1933年由苏联工程师科捷利尼科夫首次用公式严格地表述这一定理,因此在苏联文献中称为科捷利尼科夫采样定理。1948年信息论的创始人C.E.香农对这一定理加以明确地说明并正式作为定理引用,因此在许多文献中又称为香农采样定理。采样定理有许多表述形式,但最基本的表述方式是时域采样定理和频域采样定理。采样定理在数字式遥测系统、时分制遥测系统、信息处理、数字通信和采样控制理论等领域得到广泛的应用。

时域采样定理

频带为F的连续信号f(t)可用一系列离散的采样值f(t1),f(tΔt),f(t1±2Δt),...来表示,只要这些采样点的时间间隔Δt≤1/2F,便可根据各采样值完全恢复原来的信号f(t)。

窗函数和频率取样法的区别是什么,窗函数频谱为什么是无限的(3)

模拟信号采样示意图

时域采样定理的另一种表述方式是:当时间信号函数f(t)的最高频率分量为Fm时,f(t)的值可由一系列采样间隔小于或等于1/2fM的采样值来确定,即采样点的重复频率f≥2fM。时域采样定理是采样误差理论、随机变量采样理论和多变量采样理论的基础。

频域采样定理

对于时间上受限制的连续信号f(t)(即当│t│>T时,f(t)=0。这里T=T2-T1是信号的持续时间),若其频谱为F(ω),则可在频域上用一系列离散的采样值用下式表示

窗函数和频率取样法的区别是什么,窗函数频谱为什么是无限的(4)

首页 123下一页

栏目热文

eoq模型的8个假设条件(eoq模型的需求和订货数量怎么求)

eoq模型的8个假设条件(eoq模型的需求和订货数量怎么求)

关注泽稷ACCA,每日分享ACCA新鲜干货与财经资讯  在MA中,我们学习过一个很重要的关于存货管理的模型-EOQ (E...

2023-11-18 11:19:31查看全文 >>

aepc指的是什么(aebs功能是什么意思)

aepc指的是什么(aebs功能是什么意思)

阅读此文大概需要2-3分钟,诚邀您点击右上角“关注”按钮,先点后看,养成习惯,您的支持与鼓励便是我们创作的最大动力!近期...

2023-11-18 10:38:50查看全文 >>

队列研究中ar值是指什么(前瞻性回顾性队列研究怎么区分)

队列研究中ar值是指什么(前瞻性回顾性队列研究怎么区分)

公卫执业医师考点总结——队列研究。队列研究,一直是预防专业中流行病学的一个重点章节,其中包含众多考点,如队列研究原理和基...

2023-11-18 11:10:29查看全文 >>

aect定义是什么(AECT定义)

aect定义是什么(AECT定义)

上文结尾说到了94定义和PPT的关系。如果教育技术发展史是一条历史长河的话,那么它们两个谁是主河道、谁是支流呢?怎么看,...

2023-11-18 10:53:43查看全文 >>

aect2005定义中的研究内容是(aect94定义和aect05定义有何异同)

aect2005定义中的研究内容是(aect94定义和aect05定义有何异同)

作者:何克抗摘要 文章首先介绍“教育供给侧改革”的内涵及其产生背景,然后分析当前我国教育技术学是否正在脱离教育创新实践,...

2023-11-18 11:02:22查看全文 >>

eoq模型的基本假设是否合理(eoq模型的基本内容是什么)

eoq模型的基本假设是否合理(eoq模型的基本内容是什么)

在库存管理中,经济订货批量Economic Order Quantity(简称EOQ)是最经典的生产计划和库存管理模型之...

2023-11-18 10:41:03查看全文 >>

通信中的cs域和ps域是什么意思(通信中的cs和ps都是什么意思)

通信中的cs域和ps域是什么意思(通信中的cs和ps都是什么意思)

对从事通信行业的工程师来说,不管是无线还是核心网,要想拿高薪走上人生巅峰,读协议是必经之路。但3GPP规范涉及到的东西太...

2023-11-18 10:52:44查看全文 >>

墨镜严重划痕怎么修复(墨镜划痕修复最佳方法)

墨镜严重划痕怎么修复(墨镜划痕修复最佳方法)

修复眼镜上的划痕可能有一些方法,但请注意,这些方法可能并不适用于所有类型的眼镜或严重的划痕。以下是一些尝试修复眼镜划痕的...

2023-11-18 11:22:07查看全文 >>

墨镜镜片有划痕可以修复吗

墨镜镜片有划痕可以修复吗

墨镜上出现划痕。(来源:受访者)信网6月2日讯(记者 郑倩楠) 王女士在天猫COACH官方outlet旗舰店购买墨镜后,...

2023-11-18 11:21:44查看全文 >>

塑料镜片有划痕怎么修复(镜片有划痕怎么修复小妙招)

塑料镜片有划痕怎么修复(镜片有划痕怎么修复小妙招)

你知道吗,根据数据统计超过6亿中国人近视,近70%的学生佩戴眼镜?近视问题是否有如一层阴影笼罩在你们的头顶?同时,大多数...

2023-11-18 10:55:01查看全文 >>

文档排行