快速傅里叶变换实验解说.docx

实验二：快速傅里叶变换 一．基本信号（函数）的FFT变换 1. 所利用的matlab代码基本如下： clc;clear; fs=input(选择采样频率：); N=input(选择采样点数：); n=0:N-1;t=n/fs; x=sin(2*pi*t+pi/6)+sin(4*pi*t)+cos(6*pi*t); subplot(3,1,1) stem(t,x) title([采样频率为 num2str(fs) 截断长度为 num2str(N) 的采样信号]); xlabel(t/T); ylabel(x); grid on; y=fft(x,N); y=fftshift(y)/N; mag=abs(y); phase=angle(y)*180/pi; f=n*fs/N-fs/2; subplot(3,1,2) stem(f,mag) title(快速傅里叶变换后的幅频图) xlabel(f/f0); ylabel(S); grid on; subplot(3,1,3) stem(f,phase) title(快速傅里叶变换后的相频图) xlabel(f/f0); ylabel(Φ); grid on; figure(1); （1） 采样频率fs=8f0，截断长度N=16 分析：信号中最高频率fc=3f0，采样频率为fs=8f0?，满足采样定理。采样点数N=16???根据?f=fsN=1T知，分辨率为?f=0.5f0，窗函数宽度T=2T0。 （2） 采样频率fs=8f0，截断长度N=32 分析：信号中最高频率fc=3f0，采样频率为fs=8f0?，满足采样定理。采样点数N=32，根据?f=fsN=1T知，分辨率为?f=0.25f0，窗函数宽度T=4T0。 频谱的分辨率为?f=0.25f0，在输出结果中对应频率为-3f0、-2f0、-f0、f0、2f0、3f0为需要的谱线。 所得图像中并没有混叠现象。 理论上应该存在泄漏效应，但在整周期截取的情况下，旁瓣上采样都约为0，泄漏现象并没有体现。 是整周期截取，并没有栅栏效应 xt=12cos2πf0-π3+12cos-2πf0+π3+12cos?(4πf0-π2)+ 12cos?(-4πf0+π2)+ 12cos6πf0+ 12cos?(-6πf0) 所以f在-3f0、-2f0、-f0、f0、2f0、3f0上的幅值均为12，相角分别为0、π2、π3、-π3、-π2、0。可见实验值与理论值之间并没有误差。 2. （1） 采样频率fs=8f0，截断长度N=16 分析：信号中最高频率fc=11f0，采样频率为fs=8f0?，不满足采样定理。采样点数N=16，根据?f=fsN=1T知，分辨率为?f=0.5f0，窗函数宽度T=2T0。 （2） 采样频率fs=32f0，截断长度N=32 分析：信号中最高频率fc=11f0，采样频率为fs=32f0?，满足采样定理。采样点数N=32，根据?f=fsN=1T知，分辨率为?f=f0，窗函数宽度T=T0。 在（1）中分辨率为?f=0.5f0，（2）中分辨率为?f=f0，在输出结果中-11f0、-f0、f0、11f0是要求的谱线 （1）中有混叠现象，-3f0、3f0处的谱线就是混叠产生的结果，为了消除混叠需要将采样频率提高到要求的最大频率信号的两倍以上。 理论上应该有泄漏现象。但是由于整周期截取，旁瓣上的采样都为0，泄漏现象没有体现出来。 均为整周期截取，没有栅栏效应。 类似1中，在-11f0、-f0、f0、11f0的理论幅值均为0.5，理论相角为π2、π3、-π3、-π2。（1）中在-f0、f0处并没有幅值误差与相角误差。（2）中在-11f0、-f0、f0、11f0也没有幅值误差与相角误差。 3. （1） （2） （1）中的分辨率为?f=0.5f0，（2）中的分辨率为?f=f0，在输出结果中-3f0、3f0所在的谱线是所要求的谱线。 没有混叠现象 从图中可以看出，除了两条幅值较大的谱线外，还有一些幅值较小的谱线，说明发生了泄漏现象。主要原因是截断时加了矩形窗。与（1）相比，（2）的窗宽度减小，主瓣变宽，能量更加分散，而其旁瓣却被压低，幅度A明显减小。泄漏使能量分布变得分散，使要求的谱线能量降低（幅值减小）。为减少泄漏的影响，可以选择性能更好的特殊窗来代替矩形窗，即进行加窗处理。? 由于不是整周期截取，所以出现了栅栏效应。整周期截取就可以避免栅栏效应。 在（1）中，-3f0、3f0的幅值为0.4548，?A=0.0452，其相角为?59.9520 即?φ=59.9520 在（2）中，-3f0、3f0的幅值为0.4910，?A=0.0090，其相角为?27.6898 即?φ=27.6898 4. （1） （1