实验四__信号处理中FFT的应用析.doc

实验名称：实验四 信号处理中FFT的应用析一、实验目的1、理解用FFT对周期序列进行频谱分析时所面临的问题并掌握其解决方法。2、掌握用时域窗函数加权处理的技术。3、理解用FFT对非周期信号进行频谱分析所面临的问题并掌握其解决方法。二、实验原理与计算方法、对周期序列进行频谱分析应注意的问题对时间序列作FFT时，实际上要作周期延拓（如果取长序列的一段进行计算还要先作截断）。周期序列是无限长时间序列，如果截断区间刚好就是该序列周期的整数倍，那么在进行周期延拓后，将还原出原来的周期序列，由此可以较精确地计算出的该周期序列的频谱。反之，如果截断区间并不是该序列周期的整数倍，那么在进行周期延拓后，就不可能还原出原来的周期序列，由此计算出的频谱与该周期序列的频谱存在误差，而且误差的大小与截断区间的选取直接相关，如图4-1所示，其中幅度频谱的量值表示为，以dB(分贝)为单位。这种误差是由于周期序列与矩形截断序列相乘在频域产生二者的频谱卷积形成的。矩形窗的频谱是抽样函数序列。除了k = 0处主瓣内集中了大部分的能量外，两旁的较小峰值处的旁瓣也分散了一部分能量，它与周期序列频谱卷积的结果使原来集中的频谱展宽，称为频率泄漏。如果对已知周期的信号作频谱分析，在进行时域截断时，完全可以选取其周期的整数倍裁取，从而可以避免这种频率泄漏的发生。不过，通常需要进行频谱分析的信号是周期未知的信号，或随机信号，无法判断它的周期值，为了尽量避免频率泄漏对结果的影响，在作时间截断时，就应选取其频谱的旁瓣较小的截断函数，以减轻泄漏问题。2、时域窗函数的应用作为截断函数，矩形窗在作时间截断时，对所截取区间内的信号不加以任何影响，而其它的窗函数都将对所截取区间内的信号作加权处理。除了在实验二中已经介绍过三角窗、Hanning窗和Hamming窗外，常用的窗函数还有很多，例如Parzen窗、Kaiser窗、Chebyshev窗、Tukey窗、Poisson窗、Caushy窗、Gaussian窗和Blackman窗等等。本次实验仍是采用实验二中的几种窗函数，但是利用窗函数作时域加权截断。3、对非周期序列进行频谱分析应注意的问题（）混叠一般非周期信号作FFT之前要进行时域采样和周期延拓(无限长时间信号还应先截断再延拓)。根据Fourier变换理论，经等周期的冲激采样后，离散序列的频谱是原信号频谱以为周期的周期延拓。按照Nyquist采样定理，由采样引起周期延拓后频谱之间不发生混叠的条件是：(1)原信号应该是有限带宽信号，设其频带宽度为fm；(2)频谱的周期，即采样周期应满足Nyquist 条件。由于实际上有限长时间信号的FT是频域的无限函数，因此采样所得的离散序列的频谱必定产生混叠，减小采样周期只能减小而不能消除混叠。对于时间有限函数，当采样周期较大时，也会在FFT得到的频域出现混叠，形成频谱失真，造成频谱分析结果与原信号的实际频谱的差异，也无法恢复出原信号。当然，实际工作中只要采样和截断产生的误差在许可的范围内就行了，但应该认识到混叠是引起频谱分析误差的一个主要原因。还应该注意的是，离散Fourier变换的频域也是周期化的，区间内的样点实际上是负频率区的量值，因此如果出现混叠，就将在一个周期内出现，并发生在附近的区域，如图4-4所示。要减少混叠，就要尽量减小采样周期。（）泄漏周期函数截断引起的频率泄漏问题，在非周期函数截断处理后同样存在，这种误差是由于采样后的离散序列与矩形截断序列相乘在频域造成二者的频谱卷积形成的。矩形窗的频谱是抽样函数序列，它与离散序列频谱卷积的结果使原来集中于每一个样点处的频谱展宽，其影响在高频区(接近N/2的样点)特别明显，如图4-5所示。同样，为了尽量避免频率泄漏对结果的影响，在对非周期函数作时间截断时，除尽量增加截断序列的宽度外，也应选取其频谱的旁瓣较小的截断函数，以减轻泄漏问题。在选取了适当的窗函数后，应当使窗函数的宽度与被处理的序列长度相同，如果作变换前还需要补零（例如为了作卷积运算或避免栅栏效应），则应将原序列与窗函数相乘后再补零，即补零的样点不用窗函数加权处理。（）栅栏效应非周期信号应具有连续的频谱，在对作抽样后进行DFT，得到的是离散的频谱。如果排除混叠和泄漏等误差的影响，所得的结果也只是的连续频谱上的个样值。这就象通过栅栏上的等间距缝隙观看到的另一边的景象，故此称栅栏效应。被栅栏遮住的景象中有可能存在与显现出的频谱差异较大的变化，即显示信号特征的频谱分量。为了使被栅栏遮住的部分能尽可能地显现出来，可以采用增加频域样点密度的方法，即在不增加信号采样点的情况下，用时域补零加宽变换尺度N来实现，称为补零重构。例如原来信号采样得到12个样点，在其后面再加上4个零，使序列的总长度为16个样点。这样处理的结果原来信号的采样间隔和频率都没有改变，设采样频率为，经补零重构之后，采样频率仍然为，但是原来频域样点间宽度为/12，经补零重构之后频域样点间宽度为/16。这就使补零重构之后频域样点密度增加，而且显示出原来没有显露的一些频率位置的频谱。三、实验内容（1）将数字频率f=1/53，样点数为256的余弦序列作FFT，画出频谱曲线，观察并记录频率泄漏现象，然后用Hamming窗和三角窗作加权截断，观察并记录泄漏的衰减。所用程序： n=0:255;N=256;Ts=1./53;t=n.*Ts;xn=cos(2.*pi.*t);Xn=fft(xn);magXn=abs(Xn);k=0:1:255;w1=hamming(N);w2=bartlett(N);H=xn.*w1;B=xn.*w2;Xh=fft(H);magXh=abs(Xh);Xb=fft(B);magXb=abs(Xb);subplot(6,1,1);stem(n,xn,.);title(xn的离散函数图形)subplot(6,1,2);stem(k,magXn,.);title(xn频谱曲线的图形)subplot(6,1,3);stem(n,H,.);title(Hamming窗加权)subplot(6,1,4);stem(k,magXh,.);title(Hamming窗加权频谱曲线的图形)subplot(6,1,5);stem(n,B,.);title(三角窗加权)subplot(6,1,6);stem(k,magXb,.);title(三角窗加权频谱曲线)产生的图形：（2）将幅度为1，周期为2的方波信号，按Ts=1/37 s的间距抽样，取样点数N=256作FFT，画出频谱曲线，然后用Hamming窗和三角窗作加权截断，观察并记录作不同的加权截断引起的频谱差异。所用的程序： n=0:255;N=256;Ts=1./37;t=n.*Ts;xn=square(1.*pi.*t,50);Xn=fft(xn);magXn=abs(Xn);k=0:1:255;w1=hamming(N);w2=bartlett(N);H=xn.*w1;B=xn.*w2;Xh=fft(H);magXh=abs(Xh);Xb=fft(B);magXb=abs(Xb);subplot(6,1,1);stem(n,xn,r.);title(xn的离散函数图形)subplot(6,1,2);stem(k,magXn,.);title(xn频谱曲线的图形)subplot(6,1,3);stem(n,H,m.);title(Hamming窗加权)subplot(6,1,4);stem(k,magXh,.);title(Hamming窗加权频谱曲线图形)subplot(6,1,5);stem(n,B,m.);title(三角窗加权)subplot(6,1,6);stem(k,magXb,.);title(三角窗加权频谱曲线) 产生的图形：（3）将单边指数函数抽样截断后作FFT，首先选取不同的抽样周期s，并取N = 256，观察频谱混叠。然后作不同宽度的截断，选取矩形窗宽为4，8，16，32等，并保持N = 128，观察频谱泄漏。说明在什么情况下混叠和泄漏基本消失。所用的程序:Ts1=input(输入周期值Ts1=)Ts2=input(输入周期值Ts2=)Ts3=input(输入周期值Ts3=)Ts4=input(输入周期值Ts4=)n1=0:3;N=256;t1=n1.*Ts1; xn1=exp(-t1).*(t1=0);B1=boxcar(4);H1=xn1.*B1;Xn1=fft(xn1);magXn1=abs(Xn1);magHw1=abs(fft(H1,N);subplot(4,2,1);stem(magXn1,.);title(Ts1=0.1抽样函数)subplot(4,2,2);stem(magHw1,.);title(矩形窗宽4加权频谱曲线)n2=0:7;N=256;t2=n2.*Ts2;xn2=exp(-t2).*(t2=0);B2=boxcar(8);H2=xn2.*B2;magXn2=abs(Xn2);magHw2=abs(fft(H2,N);subplot(4,2,3);stem(magXn2,.);title(Ts2=0.5抽样函数)subplot(4,2,4);stem(magHw2,.);title(矩形窗宽8加权频谱曲线)n3=0:15;N=256;t3=n3.*Ts3xn3=exp(-t3).*(t3=0);B3=boxcar(16);H3=xn3.*B3;Xn3=fft(xn3);magXn3=abs(Xn3);magHw3=abs(fft(H3,N);subplot(4,2,5);stem(magXn3,.);title(Ts3=1抽样函数)subplot(4,2,6);stem(magHw3,.);title(矩形窗宽16加权频谱曲线)n4=0:31;N=256;t4=n4.*Ts4;xn4=exp(-t4).*(t4=0);B4=boxcar(32);H4=xn4.*B4;Xn4=fft(xn4);magXn4=abs(Xn4);magHw4=abs(fft(H4,N);subplot(4,2,7);stem(magXn4,.);title(Ts4=2抽样函数)subplot(4,2,8);stem(magHw4,.);title(矩形窗宽32加权频谱曲线) 产生的图形:（4）计算下面三个正弦函数的组合的频谱其中频率，令，抽样周期。分别取N=32,64,128将其截断后作FFT，观察和记录混叠和泄漏形态。分别采取补零加宽和增加截取时间宽度的方法作出频谱图，了解栅栏效应和频率分辨力的意义。所用程序:n=0:31;N=32;Ts=1./64;t=n.*Ts;xn=(sin(2.*pi.*6.3.*t)+sin(2.*pi.*9.7.*t)+sin(2.*pi.*15.3.*t).*(t=0);magXn1=abs(fft(xn,N);subplot(3,2,1);stem(magXn1,.);title(N=32截断图形)subplot(3,2,2);stem(n,abs(fft(xn),.)n=0:63;N=64;Ts=1./64;t=n.*Ts;xn=(sin(2.*pi.*6.3.*t)+sin(2.*pi.*9.7.*t)+sin(2.*pi.*15.3.*t).*(t=0);magXn2=abs(fft(xn,N);subplot(3,2,3);stem(magXn2,.);title(N=64截断图