离散傅里叶变换.pptVIP

图3.4.4 例3.4.1的求解程序运行结果 现在是95页\一共有171页\编辑于星期四 图3.4.5 重叠相加法时域波形 现在是96页\一共有171页\编辑于星期四 3.4.2 用DFT对信号进行谱分析 1． 用DFT对连续信号进行谱分析 　 工程实际中，经常遇到连续信号xa(t)，其频谱函数Xa(jΩ)也是连续函数。为了利用DFT对xa(t)进行频谱分析，先对xa(t)进行时域采样，得到x(n)=xa(nT)，再对x(n)进行DFT，得到的X(k)则是x(n)的傅里叶变换X(ejω)在频率区间[0，2π]上的N点等间隔采样。这里x(n)和X(k)均为有限长序列。 现在是97页\一共有171页\编辑于星期四 　　然而，由傅里叶变换理论知道，若信号持续时间有限长，则其频谱无限宽；若信号的频谱有限宽，则其持续时间必然为无限长。所以严格地讲，持续时间有限的带限信号是不存在的。因此，按采样定理采样时，上述两种情况下的采样序列x(n)=xa(nT)均应为无限长，不满足DFT的变换条件。实际上对频谱很宽的信号，为防止时域采样后产生频谱混叠失真，可用预滤波器滤除幅度较小的高频成分，使连续信号的带宽小于折叠频率。 现在是98页\一共有171页\编辑于星期四 　　对于持续时间很长的信号，采样点数太多, 以致无法存储和计算，只好截取有限点进行DFT。由上述可见，用DFT对连续信号进行频谱分析必然是近似的，其近似程度与信号带宽、采样频率和截取长度有关。实际上从工程角度看，滤除幅度很小的高频成分和截去幅度很小的部分时间信号是允许的。因此，在下面分析中，假设xa(t)是经过预滤波和截取处理的有限长带限信号。 现在是99页\一共有171页\编辑于星期四 设连续信号xa(t)持续时间为Tp，最高频率为fc，如图3.4.6(a)所示。xa(t)的傅里叶变换为 对xa(t)以采样间隔T≤(2fc)-1(即fs=1/T≥2fc)采样得 x(n)=xa(nT)。设共采样N点，并对Xa(jf )作零阶近似(t=nT, dt=T) 得 现在是100页\一共有171页\编辑于星期四 显然， Xa( jf )仍是 f 的连续周期函数，x(n)和X( jf )如图3.4.5(b)所示。 对 X(jf )在区间[0, Fs]上等间隔采样N点， 采样间隔为F， 如图3.4.5(c)所示。 参数Fs 、 Tp、 N和F满足如下关系式： 由于NT=Tp， 所以 (3.4.5) (3.4.6) 将f=kF和式(3.4.5)代入X( jf )中可得Xa( jf )的采样 0≤k≤N-1 现在是101页\一共有171页\编辑于星期四 令 则 (3.4.8) (3.4.7) 同理, 由 可以推出 现在是102页\一共有171页\编辑于星期四 图3.4.6 用DFT分析连续信号谱的原理示意图 现在是103页\一共有171页\编辑于星期四 (3.4.7)式说明，可以通过对连续信号采样并进行DFT再乘以T，近似得到模拟信号频谱的周期延拓函数在第一个周期[0, fs]上的N点等间隔采样 ，如图所示。对满足假设的持续时间有限的带限信号，在满足时域采样定理时， 包含了模拟信号频谱的全部信息（k=0, 1, 2, …, N/2, 表示正频率频谱采样; k=N/2+1，N/2+2，…，N－1, 表示负频率频谱采样）。 现在是104页\一共有171页\编辑于星期四 　　所以，上述分析方法不丢失信息，即可由X(k)恢复Xa(jΩ)或xa(t)，但直接由分析结果X(k)看不到Xa(jΩ)的全部频谱特性，而只能看到N个离散采样点的谱线，这就是所谓的栅栏效应。对实信号，其频谱函数具有共轭对称性，所以分析正频率频谱就足够了。不存在频谱混叠失真时，正频率[0, Fs/2]频谱采样为 （3.4.12） 现在是105页\一共有171页\编辑于星期四 　　值得注意，如果xa(t)持续时间无限长，上述分析中要进行截断处理，所以会产生所谓的截断效应，从而使谱分析产生误差。本节最后将讨论上述误差问题产生的原因及改进措施。  　　下面举例说明截断效应。理想低通滤波器的单位冲激响应ha(t)及其频响函数Ha(f)如图3.4.7(a)、(b)所示（图3.4.7(a)中只画出ha(t)所截取的一段）。图中，  现在是106页\一共有171页\编辑于星期四 现在用DFT来分析ha(t)的频率响应特性。由于ha(t)的持续时间为无穷长，因此要截取一段Tp，假设Tp=8 s，采样间隔T=0.25 s(即采样频率Fs=4 Hz)，采用点数N=