DFT分析信号频谱

利用DFT分析信号频谱一、问题的提出二、四种信号频谱之间的关系三、利用DFT分析连续非周期信号频谱四、混叠现象、泄漏现象、栅栏现象五、DFT参数选取六、工程实际应用1.连续时间非周期信号图1连续非周期信号及其频谱

一、问题的提出ΩX(jΩ)02.连续时间周期信号图2连续周期信号及其频谱

ΩX(nΩ0)03.离散时间非周期信号图3离散非周期信号及其频谱

wX(ejw)0......p2pp-p2-4.离散时间周期信号图4离散周期信号及其频谱

问题提出：如何利用数字方法分析信号的频谱？有限长序列的傅立叶变换DFT.DFT可以直接计算周期序列的DFS可否利用DFT分析以上四种信号的频谱？基本原理

利用信号傅立叶变换具有的信号时域与频域之间的对应关系，建立信号的DFT与四种信号频谱之间的关系。时域的离散化时域的周期化频域周期化频域离散化二、四种信号的时域与频域对应关系

FTFSDTFTDFSΩX(jΩ)0ΩX(nw0)0wX(ejw)0......p2pp-p2-离散傅里叶变换及MATLAB实现用DFT对模拟信号进行谱分析设模拟信号xa(t)的持续时间为Tp,信号最高频率为fc,其傅里叶变换为Xa(jΩ).对其以采样间隔T（采样频率为Fs=1/T）采样N点得到x(n),则有.设x(n)的离散时间傅里叶变换为，则有x(n)的N点DFTX(k)是在上的N点等间隔采样，即令F即为频率采样间隔，将其定义为频谱分辨率。离散傅里叶变换及MATLAB实现显然，频率采样间隔F越小，则频谱分辨率越高，谱分析越接近模拟信号的频谱。从式（3-84）可以总结出提高频谱分辨率的方法如下：★若采样点数N一定，则可以通过减小采样频率Fs，使F减小，但Fs仍需满足奈奎斯特采样定理。★若采样保持采样频率Fs一定，则可以通过增加N使F减小,这种方法包括增加记录时间Tp,或增加DFT点数N.【例3】设某模拟信号由3种频率成分组成，对其进行频谱分析，要求分辨出信号所包含的3种频率成分，抽样频率。即离散傅里叶变换及MATLAB实现由于信号点数取得过少，无法表征信号的全部信息，因此分辨率较差。将信号点数取为80点，分别作80、512和2048点DFT。

离散傅里叶变换及MATLAB实现此时，时域取了信号的一个周期，可以表征信号的全部信息。DFT点数取为80的时候，只分辨出了两种频率成分，增加DFT点数到512和2048时，第3种频率成分隐约能分辨出来。再修改程序，将时域信号取为512点，分别作512、1024和2048点DFT。

离散傅里叶变换及MATLAB实现继续修改程序，将时域信号取为2048点，作2048点DFT。

抽样离散化周期化三、利用DFT分析连续非周期信号的频谱

DFT实现假设连续信号持续时间有限，频带有限讨论：（1）无限长，其频带有限加窗抽样DFTmΩ-Ω)j(ΩX0AmΩp2-p2w)(jweXN0p2-p2w)(jweX0TA（2）有限长，其频带无限抽样DFTΩ)j(ΩX0Ap2-p2w)(jweX0（3）无限长，其频带无限加窗出现三种现象：混叠、泄漏、栅栏抽样DFTp2-p2w)(jweX0p2-p2w)(jweXN0Ω)j(ΩX0A四、混叠现象、泄漏现象、栅栏现象（1）混叠现象：减小抽样间隔T，抗混滤波抗混滤波抽样间隔T抽样DFTΩ)j(ΩX0AΩ)j(ΩX0A0mΩm-Ω0Ap2p2-w)(j0weX（2）泄漏现象：选择合适的窗函数其中：加窗DFT矩形窗窗函数一：时域波形幅度频谱矩形窗：w)(jweWNNNp2Np4pp-Np2-主瓣旁瓣旁瓣0矩形窗：主瓣在处有一个峰值，表示其主要是由直流分量组成。由于矩形窗函数在其两个端点的突然截断，使得频谱中存在许多高频分量。汉宁窗（hanning)：窗函数二：时域波形幅度频谱窗函数三：哈明窗（hamming)时域波形幅度频谱窗函数四：布拉克曼窗（Blackman）时域波形幅度频谱窗函数五：凯塞窗（Kaiser）时域波形幅度频谱窗函数类型主瓣宽度旁瓣峰值衰耗(dB)矩形4p

/N-13Hann8p

/N-31Hamming8p

/N-41Blackman12p

/N-57Kaiser（）10p

/N-57常用窗函数特性

为了说明时域加窗对连续信号频谱分析的影响，现分析一无穷长的余弦信号的频谱。

加窗抽样DFT加窗抽样DFTΩ)j(ΩX00Ω0Ω-(p)(p)Ω)j(32ΩX00Ω0Ω

-频率泄漏)j(64ΩXΩ00Ω0Ω-频率泄漏已知一连续信号为

若以抽样频率Hz对该信号进行抽样，试求由DFT分析其频谱时，能够分辨此两个谱峰所需的最少样本点数。

矩形窗信号样点数N=30信号样点数N=20加矩形窗利用矩形窗计算有限长余弦信号频谱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=25矩形窗N=50海明窗N=25海明窗N=50利用Hamming窗计算有限长余弦信号频谱N=50;%数据的长度L=512;%DFT的点数f1=100;f2=150;fs=600;%抽样频率T=1/fs;%抽样间隔ws=2*pi*fs;t=(0:N-1)*T;f=cos(2*pi*f1*t)+0.15*cos(2*pi*f2*t);wh=(hamming(N))';f=f.*wh;F=fftshift(fft(f,L));w=(-ws/2+(0:L-1)*ws/L)/(2*pi);plot(w,abs(F));ylabel('幅度谱')（3）栅栏现象：序列后补零，ZFFT（3）栅栏现象：序列后补零，ZFFT解：，m=0,1,2,3，m=0,1,2,…,7N=30,N=30,L=64,=600/64N=30,L=128,=600/128N=30,L=256,=600/256五、DFT参数选取1.2.抽样频率:抽样间隔:抽样时间:或记录时间：连续信号被采样的时间区间频率分辨率：DFT算法检测到的存在于待采样连续波形中相邻频率存在的能力。频率分辨率完全由抽样时间或记录时间决定。3.抽样点数:例：某模拟信号最高频率为600Hz，在最低可能的频率下采样，频率分辨率要求为0.5Hz，求采样频率、记录时间和样本数。六、工程实际应用A/DDFTx(t)xN[k]X[m]铁路控制信号识别铁路控制信号的时域波形和频谱铁路控制信号识别铁路控制信号的频谱分析铁路控制信号识别