第3章-录像.ppt

由于DFT是DTFT在单位圆上的N点等间隔采样，其频谱将不再是连续函数而是基频的整数倍。 用DFT计算频谱，就如同通过一个栅栏观看一个景像，只能在离散点的地方看到真实的景像，把这种现象称为“栅栏效应”。 如果在两离散的谱线间频谱有很大变化，若不作特殊处理，则无法将其检测出来。 （3）栅栏效应 20 0 10 20 30 40 50 60 0 10 DFT点数N=64 N1=32 DFT点数N=32 N1=32 0 5 10 15 20 25 30 35 0 10 20 0.265625*32=8.5 0.265625*64=17 栅栏现象 减小栅栏效应的方法: 在所取时域数据的末端补一些零值点，使一个周期内点数N增加，但是不改变原有的记录数据，而频域的抽样间隔 因N的增大而减小。即在保持原有频谱连续形式不变的情况下，使谱线更密。这样，就有可能看到原来看不到的频谱分量。 “频率分辨率”是信号处理中的一个基本概念。是指将信号中两个靠得很近的谱峰分开的能力。最通用的方法是用两个不同频率的正弦信号来研究分辨率的大小，若能分辨的两个正弦信号的频率越靠近，表明分辨率越高。 （4）频率分辨率 抽样频率： 记录时间： Tp 所以 xN(n)看做是无限长序列x(n)和一个长度为N的矩形窗函数相乘的结果，该矩形窗函数的主瓣宽度为 。 设x(n)是由两个频率分别为ω1和ω2的正弦序列组成，能分辨出这两个谱峰，N必须满足 主瓣宽度的两一种定义为XN(ejω)的幅度下降到最大幅度的 时频谱的宽度，对于矩形窗函数，这个宽度约为 。同时设ω2ω1,并且 则对信号做DFT时，频率分辨率限制为 物理频率分辨率 提高频率分辨率 增加信号的长度 Δf正比于Tp ，而不是N。所以在给定Tp的情况下，靠减小T增加N是不能提高频率分辨率的。因为Tp=NT，若把T减小到原来的1/m，则N相应增加到原来的倍 即Δf保持不变。 注 意 例3.7.2 已知连续信号x(t)=cos(2πf1t)+cos(2πf2t)，f1=1000Hz，f2=1015Hz。若以抽样频率fs=3000Hz对该信号抽样，求用DFT对该信号近似谱分析时，用矩形窗截断的长度最少应为多少。 根据DFT的定义，对一个有限长序列xN(n)做DFT，所得XN(k)的每两根谱线之间的距离为 这是频率分辨率的另一种度量形式。它是指连续谱XN(ejω)在单位圆上通过N点等间隔抽样得到的离散普的相邻谱线的距离，所以称为计算分辨率。 例3.7.3 设序列x(n)=cos(0.48πn)+cos(0.52πn)，画出其10点DFT、补零至64点的DFT以及64点DFT。 （5）利用DFT对连续信号进行谱分析时的参数选择 在利用DFT对连续信号进行频谱分析时，涉及频谱混叠、截断效应和栅栏效应。频谱混叠与连续信号的时域抽样频率fs有关，截断效应与时域加窗截断有关。大多数情况，一般已知待分析连续信号的最高频率fh，并对信号的频率分辨率Δf提出要求。下面根据傅里叶变换理论，讨论利用DFT对连续信号进行频谱分析的参数选择问题。 若已知信号的最高频率fh ，为防止频谱混叠，选定抽样频率fs满足 相应的抽样时间间隔T应该满足 再根据给定的频谱分辨率Δf ，可确定DFT所需的点数N，即 fs和Δf确定后，就可确定所需连续信号x(t)的长度 例 某一处理器在近似计算实数信号频谱时，假设未对数据作任何修正处理，要求满足的指标是： （1）抽样点数必须是2的整数次方； （2）频率分辨率 Hz； （3）信号的最高频率 Hz。 请选择下列参数：（1）最短的记录长度 ； （2）抽样点之间的最大间隔 ； （3）记录的最少点数 。 解 （1）依据分辨率确定最短记录长度为 （2）依据信号的最高频率来确定抽样点之间的 最大间隔为 （3）记录的最少点数为 所以选点数为 2、由 表示 及内插公式 所以 内插函数 当k=0时,则有 时, 时, ，所以 所以，内插函数满足