《离散傅里叶变换-第三章》.ppt

3.4.2 用DFT对信号进行谱分析 信号的谱分析：就是计算信号的傅里叶变换。 连续信号与系统的傅里叶分析不便于直接用计算机进行计算， 应用受到限制。 DFT是一种时域和频域均离散化的变换，适合数值运算， 成为分析离散信号和系统的有力工具。 1. 用DFT对连续信号进行谱分析 工程中经遇到的连续信号xa(t)，其频谱函数Xa(jΩ)也是连续函数。 先对xa(t)进行时域采样，得到时域离散信号x(n)=xa(nT)； 对x(n)进行DFT，得到的X(k)是x(n)的傅里叶变换X(ejw)在区间[0, 2?]上的N点等间隔采样； x(n)和X(k)均是有限长序列； DFT对xa(t)进行频谱分析 傅里叶变换理论 信号持续时间有限长，其频谱是无限宽。 信号的频谱有限长，在时域中，该信号的持续时间无限长。 上述两种情况，在时域或频域中进行采样，得到的序列都是无限长序列，不满足DFT的变换条件。 采用的处理方法：在频域中用滤波器滤除高于折叠频率的高频分量，在时域中则是截取有限点进行DFT。 结论：用DFT对连续信号进行谱分析是一种近似的分析，近似程度与信号带宽、采样频率和截取的长度有关。 3.4 DFT的应用举例 设连续信号xa(t)持续时间为Tp，最高频率为fc， 如下图(a)所示。 则xa(t)的傅里叶变换为： 3.4 DFT的应用举例 对xa(t)以采样频率fs= 1/T≥2fc进行采样得：x(n)= Xa(nT)。 设共采样N点，并对Xa(jf)作零阶近似(t=nT, dt=T)得: 对 x(jf)在区间[0, fs]上等间隔采样N点，采样间隔为F，参数fs 、 Tp、 N和F满足如下关系式： 令f=KF，频域N点采样得 : 令X(jkF)=Xa(k)，xa(nT)=x(n)，代入得 3.4 DFT的应用举例 结论: (1)连续信号的频谱特性可以通过对连续信号采样，并进行DFT再乘以T的近似方法得到。 (2)连续信号的时域采样信号可以通过对其频谱函数进行采样，并进行IDFT再乘以1/T的近似方法得到。 误差现象： (1)分析的结果看不到xa(jf)的全部特性，只能看到N个离散采样点的谱特性，这就是栅栏效应。 (2)如果持续时间无限长，分析时要进行截断处理，这样会产生频谱混叠和泄漏现象，使谱分析产生误差。 3.4 DFT的应用举例 【例】理想低通滤波器的单位冲激响应ha(t)及其频响函数Ha(if)如图所示。 3.4 DFT的应用举例 对连续信号进行谱分析主要关心的两个问题： 谱分析的范围fc ：受采样频率fs的限制，fc fs / 2，确保不产生频谱混叠失真。 频率分辨率：用采样间隔F来描述，表示谱分析中能够分辩两个频谱分量的最小间隔。F越小，频谱分辨率越高。 用DFT对连续信号进行谱分析的参数选择原则： 采样频率fs： fs 2 fc 谱分辩率: F= fs / N 采样点数N的选择： N ? 2fc/F 信号观察时间Tp的选择： Tp ? 1/F 例：对实信号进行谱分析，要求谱分辨率F≤10 Hz，信号最高频率fc=2.5kHz，试确定最小记录时间TPmin，最大的采样间隔Tmax，最少的采样点数Nmin。如果fc不变，要求谱分辨率增加一倍，最少的采样点数和最小的记录时间是多少? 解：根据信号观察时间TP的选择原则：TP ≥1/F=1/10=0.1s 因为要求: fs≥2fc，最小的采样频率为2fc ，所以: 频率分辨率提高一倍， 即：F=5 Hz TPmin = 1/ F = 1/5 = 0.2s 2. 用DFT对序列进行谱分析 单位圆上的Z变换就是序列傅里叶变换。 X(ejw)是w的连续周期函数，对序列x(n)进行N点DFT，得到X(k) ，X(k) 是在区间[0, 2?]上的N点等间隔采样。 序列x(n)的傅里叶变换可利用DFT来计算。 (2)对周期序列的频谱分析 设序列 (n)=x(n+rN)是周期为N的周期序列，则其傅立叶变换为： 周期序列的频谱结构可以用离散傅里叶级数系数 表示 取 的主值序列 进行N点DFT，得到 周期序列的频谱结构也可以用其主值序列的离散傅里叶变换X(k)来表示(分析) 3.4 DFT的应用举例 周期序列的频谱结构也可以用xM(k)表示 分析: (1)只有在k=rm时，XM(rm)=m ，表示 (n)的r次谐波谱线，幅度扩大了m倍，在其它k值， XM(