我是靠谱客的博主 长情啤酒,最近开发中收集的这篇文章主要介绍【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=fs1<2fh1
一般取 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=F01


参数选择

由以上公式可见,要想增加信号的最高频率分量 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=F0fs=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=F0fs>F02fh
重要!!!

这个公式是未采用任何特殊数据处理(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) T0F01=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<2fh1=2×12×1031=4.17×105(s)
3. 在一个记录中最少点数N
N > 2 f h F 0 = 12000 N > frac{2f_h}{F_0}=12000 N>F02fh=12000
N = 2 14 = 16384 N = 2^{14} = 16384 N=214=16384

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

最后

以上就是长情啤酒为你收集整理的【DSP】计算DFT频谱的参数选择的全部内容,希望文章能够帮你解决【DSP】计算DFT频谱的参数选择所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

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

评论列表共有 0 条评论

立即
投稿
返回
顶部