FFT中频率和实际频率的关系.docxVIP

FFT中频率和实际频率的关系 一，四个名词：实际物理频率，角频率，圆周频率，归一化频率 ? 实际物理频率表示ad采集物理信号的频率，FS表示采样频率，可由奈奎斯特采样定理确定 知道，fs必须≥信号最高频率的2倍才不会发生信号混叠，因此fs能采样到的信号最高频率为fs/2。 ? 角频率是物理频率的2*PI倍，也就是模拟频率。（Lu注：由于信号周期（例如 流电）是360度，即2pi。故角频率就是转了多少个2pi。设置角频率纯粹为了便于计算。） ? 归一化频率是根据FS对实际物理频率进行归一化的结果，最高信号频率为FS/2 一化频率0.5，这也就是为什么在matlab的fdtool工具中归一化频率为什么最大只到0.5的原因。 ? 周向频率是归一化频率的2*PI倍，也就是数字频率。也就是说，标准化的角频率。 二，有关fft频率与实际物理频率的分析 进行N点FFT意味着在时域中从原始信号中提取N点进行频谱分析，N点FFT变换的结果仍然是N点。 换句话说，就是将2pi数字频率w分成n份，而整个数字频率w的范围覆盖了从0-2pi*fs的模拟频率范围。这里的fs是采样频率。而我们通常只关心0-pi中的频谱，因为根据奈科斯特定律，只有f=fs/2范围内的信号才是被采样到的有效信号。那么，在w的范围内，得到的频谱肯定是关于n/2对称的。 例如，如果对16个点进行FFT分析，原始模拟信号的最高频率为f=32kHz，采样频率为64khz，n的范围为0,1,2 15。（Lu注：这意味着原始模拟信号已采样8次。）此时，64khz的模拟频率被分成16个点，每个点都是4kHz，这称为频率分辨率（Lu注：用于FFT的点越多，频率分辨率越高）。在横坐标中，当n=1时，f为4kHz，当n=2时，f为8kHz，当n=15时，f为60KHZ。你的光谱在n=8左右是对称的。你只需要关心n=0到7之间的频谱，因为原始信号的最高模拟频率是32kHz。有两个结论。 ? 首先，我们必须知道原始信号的采样频率FS，然后才能知道每个n对应的实际频率 是多少，第k个点的实际频率的计算为f(k)=k*(fs/n) ? 其次，在64khz下进行16点FFT后，因为频率分辨率是4kHz，如果原始信号 在5khz或者63khz有分量，你在频谱上是看不见的，这就表示你越想频谱画得逼真，就必须取越多的点数来做fft，n就越大，你在时域上就必须取更长的信号样本来做分析。但是无论如何，由于离散采样的原理，你不可能完全准确地画出原来连续时间信号的真实频谱，只能无限接近（就是n无限大的时候），这个就叫做频率泄露。在采样频率fs不变得情况下，频率泄漏可以通过取更多的点来改善，也可以通过做fft前加窗来改善，这就是另外一个话题了。 离散信号傅里叶变换周期性的讨论 （卢注：从下图可以看出： S平面，相当于直角坐标系，其实轴为复数的实部，虚轴为复数的虚部。这里，它可以理解为该频率下信号的振幅； z平面，相当于极坐标系，与re轴的夹角相当于频率，向量长度相当于幅值。）要分析这个，我们先从laplace变换与z变换之间的关系谈起。由 得到了z平面与s平面的关系图 图中的关系有以下几点： ? s平面的虚轴是否映射到z平面的单位圆？s平面的负半轴映射到z平面的单位圆？s平面的正半轴映射到z平面单位圆的外侧 laplace变换是用于连续信号的变换，相对应的z变换是应用到z平面的变换。因此从另一个角度，上面谈到的角频率（模拟频率）对应的是s平面，圆周频率对应的是z平面（也是为什么称为圆周频率的原因）。 现在让我们看看模拟频率在s平面虚轴上的变换如何导致z平面单位圆的变化：？当模拟频率在s平面的虚轴上从0变为FS时，数字频率在z平面的单位圆上从0变为2pi。 ?当模拟频率在s平面的虚轴上从2fs变到4fs时，数字频率在z平面单位圆上仍然从0变到 2便士 。。。。。。z平面如此循环重复 我们知道离散信号的傅里叶变换对应于单位圆上的z变换，因此上述结论验证了离散信号的傅里叶变换为什么是周期性的：根本原因是单位圆上的周期性。 考虑到我们实际应用中可选择一个周期，这也能够解释：因为实际信号的频率总是在fs/2以下，这就对应到z平面单位圆上的0~pi，在一个周期范围内就可以进行信号分析了。