我是靠谱客的博主 长情啤酒,这篇文章主要介绍【DSP】计算DFT频谱的参数选择,现在分享给大家,希望可以做个参考。

文章目录

    • 时域离散化(满足采样定理)
    • 频域离散化
    • 参数选择
    • 例:音乐频谱分析参数选择

利用DFT在计算连续时间信号时,由于时域的采样往往会出现一些问题,下面就DFT的一些参数做一些探讨。

时域离散化(满足采样定理)

若信号的最高频率为 f h f_h fh​, 根据奈奎斯特采样定理,抽样频率(采样率)应满足:
f s > 2 f h f_s > 2f_h fs​>2fh​
也就是抽样间隔 T T T应满足:
T = 1 f s < 1 2 f h T = frac{1}{f_s}<frac{1}{2f_h} T=fs​1​<2fh​1​
一般取 f s = (  from  2.5  to  3.0 ) f h f_s=( text{ from }2.5 text{ to } 3.0)f_h fs​=( from 2.5 to 3.0)fh​

频域离散化

对于DFT来说,时域和频域都是离散的。
因此频率函数也要抽样,变成离散的序列,其抽样间隔为 F 0 F_0 F0​,这就是我们所说的频率分辨率。
由它可引出时间函数的周期,也就是所取的记录长度 T 0 = 1 F 0 T_0=frac{1}{F_0} T0​=F0​1​


参数选择

由以上公式可见,要想增加信号的最高频率分量 f h f_h fh​,则时域的抽样间隔就一定减小,而 f s f_s fs​增加,由于抽样点数满足:
N = f s F 0 = T 0 T N= frac{f_s}{F_0} = frac{T_0}{T} N=F0​fs​​=TT0​​
若是在 N N N固定的情况下, F 0 F_0 F0​必然增加,即频率分辨力下降。
反之要提高频率分辨力(减小 F 0 F_0 F0​),就要增加记录长度 T 0 T_0 T0​( f s f_s fs​减小)。
若要保证不发生频谱混叠,则必然会减小高频容量 f h f_h fh​。

要想兼顾高频容量(大 f h f_h fh​)和频谱分辨率(小或不变 F 0 F_0 F0​),唯一方法就是增加记录长度的点数N,即满足:

N = f s F 0 > 2 f h F 0 N= frac{f_s}{F_0} > frac{2f_h}{F_0} N=F0​fs​​>F0​2fh​​
重要!!!

这个公式是未采用任何特殊数据处理(e.g. 加窗)的情况下,为实现基本DFT算法必须满足的最低条件。
如果加窗处理,相当于时域与窗函数相乘,则频域周期卷积,则必然加宽频谱分量。频率分辨力就可能变坏,为了保证频率分辨力不变,则必须增加记录长度,也就是增加数据长度 T 0 T_0 T0​。

例:音乐频谱分析参数选择

用FFT进行针对音质的频谱分析,88键钢琴的最高音是四千多赫兹,很多中频到高频的乐器,2、3次泛音就在4-12kHz,是主要泛音区域,是音色质感的主要表现区域。要求抽样点数必须是2的整数幂,假定没有任何特殊的数据处理措施。已给条件为频率分辨力<=2Hz, 信号最高频率<=12kHz。试确定以下参量:

1. 最小记录长度 T 0 T_0 T0​
最小记录长度由频率分辨力确定,
T 0 ≥ 1 F 0 = 1 2 = 0.5 ( s ) T_0 ge frac{1}{F_0}=frac{1}{2}=0.5(s) T0​≥F0​1​=21​=0.5(s)
2. 抽样点间的最大抽样时间间隔T(即最小采样率)
f s > 2 f h = 24 k H z f_s > 2f_h=24kHz fs​>2fh​=24kHz
T < 1 2 f h = 1 2 × 12 × 1 0 3 = 4.17 × 1 0 − 5 ( s ) T< frac{1}{2f_h}=frac{1}{2 times 12 times10^{3}}=4.17times10^{-5}(s) T<2fh​1​=2×12×1031​=4.17×10−5(s)
3. 在一个记录中最少点数N
N > 2 f h F 0 = 12000 N > frac{2f_h}{F_0}=12000 N>F0​2fh​​=12000
N = 2 14 = 16384 N = 2^{14} = 16384 N=214=16384

参考书:数字信号处理 程佩青

最后

以上就是长情啤酒最近收集整理的关于【DSP】计算DFT频谱的参数选择的全部内容,更多相关【DSP】计算DFT频谱内容请搜索靠谱客的其他文章。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(298)

评论列表共有 0 条评论

立即
投稿
返回
顶部