数字信号处理_8利用DFT分析连续信号频谱.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 解： 例： ，m=0,1,2,…,7 解： 例： N=30, N=30, L=64, = 600/64 N=30, L=128, = 600/128 N=30, L=256, = 600/256 k 0 1 N-1 m 0 1 N-1 DFT参数选取 1. 2. 3. 抽样频率: 抽样间隔: 抽样时间: 抽样点数: 例： 试利用DFT分析一连续信号，已知其最高频率=1000Hz，要求频率分辨率Dfd?2Hz，DFT的点数必须为2的整数幂次，确定以下参数：最大的抽样间隔，最少的信号持续时间，最少的DFT点数。 解：(1)最大的抽样间隔Tmax为 (2)最少的信号持续时间Tpmin为 (3) 由最大的抽样间隔Tmax与最少信号持续时间Tpmin，可得最少DFT点数N为 选择DFT的点数N=1024，以满足其为2的整数幂次。 思考题 在利用DFT分析连续非周期信号频谱过程中， (2) 如何合理选择窗函数？ (1) 如果由于截短信号造成泄漏而导致频谱分辨率下降，可否通过在截短后序列补零得到改善？ 本章小结 四种信号傅里叶变换的特点及数学公式 有限长序列DFT的定义与物理含义 DFT的性质(循环位移、对称性和循环卷积) 序列DFT、DTFT及Z变换的相互关系 有限长线性卷积与循环卷积关系。 分段卷积常用方法(重叠法相加法与重叠保留法) DFT分析信号频谱出现三种现象，以及改进方法 * * * * * * * * * * * * * * * * * * * * * * * * 利用DFT分析连续非周期信号的频谱 3. 无限长，其频带无限 加窗 出现三种现象：混叠、泄漏、栅栏 抽样 DFT 混叠现象、泄漏现象、栅栏现象 （1）. 混叠现象：减小抽样间隔T，抗混滤波 抗混滤波 抽样间隔T 抽样 DFT 混叠现象、泄漏现象、栅栏现象 （2）. 泄漏现象：选择合适的窗函数 加窗 DFT 其中： ? ? ? ? ? í ì = 凯塞窗 布拉克曼窗 哈明窗 汉宁窗 矩形窗 ] [ k w N 频谱泄漏 在实际应用中，通常将所观测的信号 x[k]限制在一定的时间间隔内，也就是说，在时域对信号进行截短操作。或称作加时间窗，亦即用时间窗函数乘以信号，由卷积定理可知，时域相乘，频域为卷积，这就造成拖尾现象，称之为频谱泄漏。 改善方法： 1）增加x[k]长度 2）缓慢截短 混叠现象、泄漏现象、栅栏现象 2. 泄漏现象： 矩形窗 窗函数一： 时域波形 幅度频谱 矩形窗： 矩形窗： 主瓣在?=0处有一个峰值，表示其主要是由直流分量组成。由于矩形窗函数在其两个端点的突然截断，使得频谱中存在许多高频分量。 例：为了说明时域加窗对连续信号频谱分析的影响，现分析一无穷长的余弦信号的频谱。 加窗 抽样 DFT 加窗 抽样 DFT 结 论 加窗处理对信号频谱分析影响： 加窗后的频谱存在更多高频成份。 谱线变成了具有一定宽度的谱峰，降低了频率分辨率。 书81公式(2-79)理解 结 论 窗函数主瓣的有效宽度 采样周期 矩形窗主瓣的有效宽度 例：已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析 其频谱时，能够分辨此两个谱峰所需的最少样本点数。 矩形窗 例：已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析 其频谱时，能够分辨此两个谱峰所需的最少样本点数。 利用矩形窗计算有限长余弦信号频谱 N=30; %数据的长度 L=512; %DFT的点数 f1=100; f2=120; fs=600; %抽样频率 T=1/fs; %抽样间隔 ws=2*pi*fs; t=(0:N-1)*T; f=cos(2*pi*f1*t)+cos(2*pi*f2*t); F=fftshift(fft(f,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); hd=plot(w,abs(F)); ylabel(幅度谱) 例：已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析 其频谱时，能够分辨此两个谱峰所需的最少样本点数。 信号样点数 N=30 信号样点数 N=20 加矩形窗 在选择窗宽度时，所要主要问题： 增加N可以减少窗函数频谱的主瓣宽度