3.7利用FFT分析时域连续信号频谱.ppt

* 第3章 快速傅里叶变换 3.7 利用FFT分析时域连续信号频谱 3.7.1 基本步骤 图 3-21 时域连续信号离散傅里叶分析的处理步骤 在图3-21中，前置低通滤波器LPF（预滤波器）的引入，是为了消除或减少时域连续信号转换成序列时可能出现的频谱混叠的影响。 在实际工作中，时域离散信号x(n)的时宽是很长的， 甚至是无限长的（例如语音或音乐信号）。由于DFT的需要(实际应用FFT计算)，必须把x(n)限制在一定的时间区间之内，即进行数据截断。 数据的截断相当于加窗处理（窗函数见6.2节）。 因此， 在计算FFT之前，用一个时域有限的窗函数w(n)加到x(n)上是非常必要的。 xc(t)通过A/D变换器转换(忽略其幅度量化误差)成采样序列x(n)，其频谱用X(ejω)表示，它是频率ω的周期函数，即 (3-47) 式中，Xc(jΩ)或 为xc(t)的频谱。 在实际应用中，前置低通滤波器的阻带不可能是无限衰减的， 故由Xc(jΩ)周期延拓得到的X(ejω)有非零重叠，即出现频谱混叠现象。  由于进行FFT的需要，必须对序列x(n)进行加窗处理，即v(n)=x(n)w(n)，加窗对频域的影响，用卷积表示如下: (3-48) 最后是进行FFT运算。 加窗后的DFT是 0≤k≤N-1 (3-49) 式中，假设窗函数长度L小于或等于DFT长度N，为进行FFT运算，这里选择N为2的整数幂次即N=2m。 有限长序列v(n)=x(n)w(n)的DFT相当于v(n)傅里叶变换的等间隔采样。 (3-50) V(k)便是sc(t)的离散频率函数。因为DFT对应的数字域频率间隔为Δω=2π/N，且模拟频率Ω和数字频率ω间的关系为ω=ΩΤ， 其中Ω=2πf。所以，离散的频率函数第k点对应的模拟频率为: (3-51) (3-52) 由式（3-52）很明显可看出，数字域频率间隔Δω=2π/N对应的模拟域谱线间距为 (3-53) 谱线间距，又称频谱分辨率（单位：Hz）。所谓频谱分辨率是指可分辨两频率的最小间距。它的意思是，如设某频谱分析的F=5Hz，那么信号中频率相差小于5 Hz的两个频率分量在此频谱图中就分辨不出来。 长度N=16 的时间信号v(n)=(1.1)nR16(n)的图形如图 3-22(a)所示, 其16点的DFT V(k)的示例如图 3 - 22(b)所示。 其中T为采样时间间隔（单位：s）；fs为采样频率（单位：Hz）；tp为截取连续时间信号的样本长度（又称记录长度，单位：s）；F为谱线间距，又称频谱分辨率（单位：Hz）。 注意：V(k)示例图给出的频率间距F及N个频率点之间的频率fs为对应的模拟域频率(单位： Hz)。 图 3 - 22 由图可知： (3-54) (3-55) 在实际应用中, 要根据信号最高频率fh和频谱分辨率F的要求， 来确定T、tp和N的大小。  （1）首先，由采样定理，为保证采样信号不失真，fs≥2fh(fh为信号频率的最高频率分量，也就是前置低通滤波器阻带的截止频率)， 即应使采样周期T满足 (3-56) （2） 由频谱分辨率F和T确定N, (3-57) 为了使用FFT运算，这里选择N为2的幂次即N=2m，由式(3-55)可知，N大，分辨率好，但会增加样本记录时间tp。 （3） 最后由N, T确定最小记录长度，tp=NT。 例 3-3 有一频谱分析用的FFT处理器，其采样点数必须是2的整数幂， 假定没有采用任何特殊的数据处理措施，已给条件为: ① 频率分辨率≤10 Hz;　② 信号最高频率≤4kHz。试确定以下参量：  (1) 最小记录长度tp； (2) 最大采样间隔T（即最小采样频率）； (3) 在一个记录中的最少点数N。 解 (1) 由分辨率的要求确定最小长度tp 所以记录长度为 (2) 从信号的最高频率确定最大可能的采样间隔T（即最小采样频率fs=1/T）。 按采样定理 即 (3) 最小记录点数N应满足 取 如果我们事先不知道信号的最高频率，可以根据信号的时域波形图来估计它。例如， 某信号的波形如图 3-23 所示。 先找出相邻的波峰与波谷之间的距离，如图中t1，t2，t3，t4。 然后，选出其中最小的一个如t4。这里, t4可能就是由信号的最高频率分量形成的。 峰与谷之间的距离就是周