基于ADSP-TS101的频谱细化

　　在信号处理分析中，最常用到的是时域分析和频域分析。其中信号的频域分析有着广阔的发展和研究空间。随着对信号处理要求的不断提高，用简单的FFT变换来对信号进行时频变换的方法在某些场合可能达不到我们所要求的频率分辨率。如果靠简单的加大FFT点数来增加分辨率无疑也大大地增加了系统的运算量，这在某些实时性要求较高的场合是不允许的。本文详细介绍了如何在快速傅里叶变换(FFT)的基础上辅以傅里叶变换来实现频谱细化。并给出了具体的设计实例加以说明。

　　1 原理介绍

　　对于一个有限长时间的离散信号序列而言，其频谱是连续的。对信号序列作离散傅里叶变化只能得到某些离散点的频率值，而且速度较慢，更多的情况下用的是快速傅里叶变换，可以很快得到信号的频谱。但是离散傅里叶变换有频率分辨率的限制，只能得到fs/N整数倍频率点的幅度和相位，对于他们之间频率的幅度和相位无法得到，这种现象称作栅栏效应。下面以如何准确求得雷达系统中的多普勒频率为例说明如下：

　　对于一个采样频率为fs，采样点数为N，的离散时间序列x(tk)，设：

　　序列的离散傅里叶级数的系数为：

　　n△f处的幅度矢量表达式为an-ibn。

　　快速离散傅里叶变换(FFT)是上述傅里叶变化的特殊情况，即N=2m(m为正整数)的时候，可以很快地计算出信号的频谱。离散傅里叶变换分辨率为△f，和采样长度N成反比，一旦N，确定，频率分辨率无法提高。

　　我们知道，离散信号序列是对连续信号的采样，采样以后的信号的谱以fs为周期重复，因此离散序列x(tk)中已经包含了原始信号中0～fs/2的频率信息，如果把离散傅里叶变换中的n看作是连续的变量，便可得到下面的式子：

　　其中：0≤f≤fs/2。

　　通过以上的式子可以得到离散信号的连续频谱，而且没有分辨率的限制，可以以任意分辨率来分析信号的频谱，下面介绍工程上如何利用上面的公式来估计出信号的准确频率。

　　2 应用实例

　　在雷达系统中需要估计出多普勒频率，以便计算出自己或敌机的飞行速度，这就需要精确地估计出回波信号的多普勒频率是多少，但对于信号的整个频谱不是很关心，因此不需要很细致地去研究信号频潜，所以可以先以粗的分辨率估计出信号的频率大概在那里，然后再用更小的分辨率来细化这段感兴趣的区间，这样就能够更加准确地计算出信号的多普勒频率。ADSP-TSl01是AD公司生产的一款高性能浮点数字信号处理器，他计算1 024点的FFT只需要用32.78 μs，内部有两个计算单元以及大容量的SRAM，能够满足实时信号处理的目的。为了计算出回波信号的多普勒频率是多少，可以这样做：先用FFT计算出信号的整个频谱，找到信号频率所在的大概位置，此时估计的频率比较粗，误差不超过