第三篇 测试信号分析3.ppt

3.3-五 离散傅立叶变换DFT 频域：对 进行周期延拓，使 之在时域上成为周期信号，周期信号的频谱是离散的；对以上信号采样前进行滤波，保证带宽有限（频域有限长度）； 在时域上对N点的采样信号进行周期延拓，就成为离散的周期信号，在频域其频谱也成为离散的周期频谱。 （时域离散-频域周期，时域周期-频域离散） 频域上的离散化实际上就是将频域上的横轴进行了离散化， ，T为时域采集周期。 DFT： 在一个周期内求其逆变换 多对信号进行周期延拓时，相当于对频谱进行采样，使频谱变成原来的1/T，故应乘以T倍。 离散傅立叶变换的特性： ★离散性：在时域具有重复周期 ，故频谱是离散的。 ★谐波性：谐波频率 是基波频率的整数倍 ★周期性：采样导致其周期性，共有N条谱线有效，前N/2条是主分量，后N/2条是高频分量。 DFT的计算次数与FFT 4个采样值 根据DFT变换可知，计算一条谱线需进行N次复数乘法和（N-1）次复数加法，N条谱线需进行N2次复数乘法和N(N-1)次复数加法，当N=4，总计需28次运算。当N=1024时，需进行2096068次运算。早期不太可能在实际工程中应用。 FFT：将长序列的离散傅立叶变换逐次分解成长度较短的序列来计算傅立叶变换，然后按一定的规律组合起来，得到原来序列的傅立叶变换。 ★简化了计算，计算量大大减少，减小了计算机字长对乘积结果舍入误差，使得计算精度提高了。 4 DFT的相关问题 f1=2Hz、f2=2.02Hz、f3=2.07Hz三个频率，fs=10Hz Δf1=0.0390625Hz, Δf2=0.009765Hz；f1-f2=-0.02,f1-f3=-0.07,f2-f3=-0.05 2、栅栏效应与补零 计算机采样根据采样频率fs和采样点数N对原信号进行了截短，信号长度为tp= fs* N，并以此长度转化成周期信号。相应的原来的连续频谱变成离散频谱，在离散频谱之间的频率分量就被挡住了而无法得到，这种现象称为栅栏效应。 但是补零并没改善频率分辨率，并没有增加数据的有效长度，分辨率是根据采样信号的实际长度来决定的，与DFT计算的点数无关。补零后的长度是DFT的长度，增加了频域谱线的点数，但谱线之间还是同样存在一定的重叠。 补零还经常用于把采样点数凑成2的整数次幂，以便于用FFT算法。谱线外观平滑，克服栅栏效应。 由于截短造成频谱泄露，频谱中出现一些难以确认的谱峰，补零可以消除这种现象。 细化FFT ZOOM-FFT 由于窗函数对原信号的截短，导致截短后的信号频谱中包含了窗函数旁瓣的频谱，这就是频谱的泄露。 谱间干扰 因截短而使在主谱线两边形成许多旁瓣，引起不同分量之间的干扰，称为谱间干扰。这不仅影响谱分辨率，严重时强信号谱的旁瓣可能淹没弱信号的主谱线，或者把强信号的旁瓣误认为另外的一个信号的谱线，形成假信号。 矩形窗谱线过零点： 当某一频率分量的输出谱线中除主分量外其他输出谱线均为0，则在与其他频率分量输出谱线叠加时不会引起其失真。即过零点谱线的频率等于输出谱线的频率。即 为了防止产生泄漏，应使窗的宽度为信号周期T的整数倍。 ★如果对信号分析中选择的采样频率和截取的信号长度不合适，就会导致分析结果不是实际的频谱，而是泄露的频谱。 频谱的泄漏主要是矩形窗边界的突变特性造成的，这种突变在频域就会产生高频成分。如果用其它类型的边缘较平缓的窗函数（如海宁窗、海明窗、三角窗等）来代替矩形窗，可以使高频分量减弱，从而减少由于泄漏产生的误差。(矩形窗主瓣与第一旁瓣的幅度比为100:25，海宁窗为100:2.5,但是海宁窗的主瓣宽度是矩形窗的两倍) 但在N一定时，旁瓣越小的窗主瓣越宽，频谱变得“模糊”，这就使减小谱间干扰与提高谱分辨率矛盾，只能折衷。 例：设余弦信号的频率为f0＝100Hz，以下采用不同的采样频率fs（或采样周期Ts）和截取不同的信号长度T（或采样点数N）对其进行分析。 ⑴ 设Ts＝78.125μs,N=256，则： 窗的宽度 ＝N Ts ＝20ms； 谱线间距 ； 谱线的过零点为：0、50、150、200…（Hz）； 输出的谱线：0、50、100、150、200，…（Hz）共计256条； 100Hz余弦信号的谱线在输出谱线之列，故2次谐波的输出即是所分析的信号的频谱，其余输出谱线均为0，此时频谱的泄漏得以消除。 3