数字信号处理第三章-FFT.ppt

用DFT对模拟信号作频谱分析 信号的频谱分析：计算信号的傅里叶变换 模拟信号进行频谱分析时，有几个重要的参数要选择， 采样频率 信号的最高截止频率 频率分辨率F0 能够分辨的两个频率分量最 小的间隔 频谱分析范围 采样点N 与对信号的观测时间有关 频率响应的混叠失真及参数的选择 参数确定 (1)由采样定理， （2）由频谱分辨率确定N （3）确定最小记录长度 用DFT（FFT）对周期信号进行谱分析的误差来源 1.频谱混叠 采样频率不满足采样定理 一般模拟信号的频谱函数并不是锐截止的 信号中总含有或多或少的干扰或噪声 2 频谱泄漏 改善方法： 比较截断前、后的幅度谱的差别： 泄露 定义：原来的离散谱线向两边展宽，这种将谱线展宽的现象称为频谱泄漏。 约束因素：矩形窗的长度越长，展宽的宽度就越窄。 影响：泄漏会使频谱模糊，谱的分辨率降低。 谱间干扰 出现原因：频谱卷积以后存在着的旁瓣，引起不同频率分量干扰。旁瓣淹没信号主谱线或者旁瓣误认为另一信号谱线。 影响：降低谱分辨率 泄漏和谱间干扰统称为信号的截断效应。 减轻截断效应的方法---提高谱分析的分辨率及精确性 。 1:增加窗时域的宽度，运算量及存贮量增大 2：其他缓变的窗，比如三角形窗等 3. 栅栏效应 定义：N点DFT（FFT）得到的只是N个采样点上的频谱值，两点之间的频谱值是不知道的，就好像被栅栏遮住一样 应对策略： 增多DFT（FFT）的变换点数 在信号尾部加零的方法加大DFT（FFT）的变换点数。谱线更密。 频率分辨率 提高频率分辨率方法： 增加信号实际记录长度 提高采样频率不能提高频率分辨率 提高了采样频率，虽然在相同的观察时长那的点数增多了，但与此同时采样频率也变大了，点数增加几倍采样频率增加几倍，所以不改变观察时长而仅仅提高采样频率并不能提高DFT谱的频率分辨率。 补零并不能提高频率分辨率 补零信号的谱，是通过截短信号的谱进行了推测（插值算法）得来的，它并不能反映原信号的谱，所以虽然补零信号的谱线间隔变小了，但是除了从截短信号的谱中取出来的谱线以外，其余的谱线都是无效的。去掉这些无效的谱线，采样频率不变，有效的谱线数不变，所以其物理频率分辨率自然没有改变。 3.8 DFT(FFT)应用举例 DFT因为具有快速算法FFT，应用非常广泛，限于篇幅，这里 仅介绍DFT在线性卷积和频谱分析两方面的应用。 本节主要论述： 3.8.1 用DFT(FFT)计算两个有限长序列的线性卷积 当h(n)或x(n)序列较长时，直接计算线性卷积的时间会很 长，满足不了实时处理的要求。可用FFT将时域卷积变换为 频域相乘来计算线性卷积，达到快速实时要求。 3.8.2 用DFT计算有限长序列与无限长序列的线性卷积 用FFT计算有限长序列与无限长序列的线性卷积 问题： ·h(n) 为某滤波器的单位脉冲响应,长度有限 ·输入信号 x(n)很长，h(n) 要补许多零再进行计算, 计算量有很大的浪费 解决方法： 重叠相加法 重叠保留法 重叠相加法： 将x(n)分段，每段长度为L，然后依次计算各段与h(n) 的卷积，再由各段的卷积结果得到y(n)。 假设h(n)和x(n)是因果序列，对x(n)进行分段，每段长为L， 则x(n)可以表示为如下形式: 的长度为N=L+M-1。用N点FFT和IFFT计算 ,再按上式求和， 得到。 计算步骤: (1)计算H(k)=DFT[h(n)] ，N=L+M-1，i=0； (2)读入 并计算 ； (3) ； (4) ，n=0, 1, 2, …, N-1； (5) 返回 * * * * 例：设x1(n)和x2(n)都是N点的实数序列，试用一次N点DFT运算来计算它们各自的DFT: 公式推导及参数选择 信号最高频率与频率分辨率之间的矛盾 例：有一调幅信号  用DFT做频谱分析，要求能分辨 的所有频率分量，问 (1)抽样频率应为多少赫兹（Hz）? (2)抽样时间间隔应为多少秒（Sec）? (3)抽样点数应为多少点？ （1）抽样频率应为 解： （2）抽样时间间隔应为 对时域截短，使频谱变宽拖尾，称为泄漏 1）增加x(n