实验二-快速傅里叶变换(FFT)及其应用

《数字信号处理》课程（2010-2011学年第1学期）成绩：实验二快速傅里叶变换（FFT）及其应用学生姓名：闫春遐所在院系：电子信息工程学院自动化系年级专业：2008级自动化系学号：00824049指导教师：王亮完成日期：2010年9月27日实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第1页共19页实验二快速傅里叶变换（FFT）及其应用一、实验目的（1）在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB中的有关函数。（2）应用FFT对典型信号进行频谱分析。（3）了解应用FFT进行信号频谱分析过程可能出现的问题，以便在实际中正确应用FFT。（4）应用FFT实现序列的线性卷积和相关。二、实验内容实验中用到的信号序列：a）高斯序列2()015()0npqaenxn其他b）衰减正弦序列sin(2)015()0anbefnnxn其他c）三角波序列03()8470cnnxnnn其他d）反三角波序列403()4470dnnxnnn其他上机实验内容：（1）观察高斯序列的时域和幅频特性，固定信号()axn中参数8p，改变q的值，使q分别等于2、4、8，观察他们的时域和幅频特性，了解当q取不同值时，对信号的时域和幅频特性的影响；固定8q，改变p，使p分别等于8、13、实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第2页共19页14，观察参数p变化对信号序列的时域及幅频特性的影响，注意p等于多少时，会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。解答:n=0:1:15;xn=exp(-(n-8).^2/2);subplot(1,2,1);stem(n,xn);xlabel('t/T');ylabel('x(n)');xk1=fft(xn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');xn=exp(-(n-8).^2/4);subplot(1,2,1);stem(n,xn);xlabel('t/T');ylabel('x(n)');xk1=fft(xn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第3页共19页xn=exp(-(n-8).^2/8);subplot(1,2,1);stem(n,xn);xlabel('t/T');ylabel('x(n)');xk1=fft(xn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');xn=exp(-(n-13).^2/8);subplot(1,2,1);stem(n,xn);xlabel('t/T');ylabel('x(n)');xk1=fft(xn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');xn=exp(-(n-14).^2/8);subplot(1,2,1);stem(n,xn);xlabel('t/T');ylabel('x(n)');xk1=fft(xn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第4页共19页随着q值的增大，时域信号幅值变化缓慢，频域信号频谱泄露程度减小。随着p的增大，时域信号幅值不变，会在时间轴移位。（2）观察衰减正弦序列()bxn的时域和幅频特性，0.1a，0.0625f，检查普峰出现的位置是否正确，注意频谱的形状，绘出幅频特性曲线，改变f，使f分别等于0.4375和0.5625，观察这两种情况下，频谱的形状和普峰出现的位置，有无混叠和泄漏现象？说明产生现象的原因。解答：n=0:1:15;xn=exp(-0.1*n).*sin(2*pi*0.0625*n);subplot(1,2,1);stem(n,xn);xlabel('t/T');ylabel('x(n)');xk1=fft(xn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第5页共19页xn=exp(-0.1*n).*sin(2*pi*0.4375*n);subplot(1,2,1);stem(n,xn);xlabel('t/T');ylabel('x(n)');xk1=fft(xn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');xn=exp(-0.1*n).*sin(2*pi*0.5625*n);subplot(1,2,1);stem(n,xn);xlabel('t/T');ylabel('x(n)');xk1=fft(xn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');（3）观察三角波和反三角波的时域和幅频特性，用8N点FFT分析信号序列()cxn和()dxn的幅频特性，观察两者的序列形状和频谱曲线有什么异同？绘出两序列及其幅频特性曲线。在()cxn和()dxn末尾补零，用32N点FFT分析这两个信号的幅频特性，观察幅频特性发生了什么变化？两种情况下的FFT频谱还有相同之处吗？这些变实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第6页共19页化说明了什么？解答：forn=0:1:3xcn(n+1)=n;end;forn=4:1:7xcn(n+1)=8-n;end;xcnxcn=01234321n=0:1:7;subplot(1,2,1);stem(n,xcn);xlabel('t/T');ylabel('x(n)');xk1=fft(xcn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');forn=0:1:3xdn(n+1)=4-n;end;实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第7页共19页forn=4:1:7xdn(n+1)=n-4;end;xdnxdn=43210123n=0:1:7;subplot(1,2,1);stem(n,xdn);xlabel('t/T');ylabel('x(n)');xk1=fft(xdn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');xcn=[xcn,zeros(1,24)];n=0:1:31;subplot(1,2,1);stem(n,xcn);xlabel('t/T');ylabel('x(n)');xk1=fft(xcn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第8页共19页xdn=[xdn,zeros(1,24)];n=0:1:31;subplot(1,2,1);stem(n,xdn);xlabel('t/T');ylabel('x(n)');xk1=fft(xdn);xk1=abs(xk1);subplot(1,2,2);stem(n,xk1);xlabel('k');ylabel('X(k)');8N时，()cxn和()dxn的幅频特性相同，在()cxn和()dxn末尾补零，用32N点FFT分析这两个信号的幅频特性时，它们还有相同之处，即当k取4的整数倍时对应幅值相等。分析：8N点FFT分析信号的幅频特性：1121()0()()*NNjnkNnXkxne32N点FFT分析信号的幅频特性：224222411()()4400()()*()*NNNjnkjnkNNnnXkxnexne实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第9页共19页由上两式可知，当k2=4k1时，两个信号的对应频率幅值相等，即对信号末尾补零加长整数个周期可以对原信号达到细化频谱的作用。（4）一个连续时间信号含两个频率分量，经采样得()sin[20.125]cos[2(0.125)]0,1,,1xnnfnnN已知16N，f分别为1/16和1/64，观察其频谱；当128N时，f不变，其结果有何不同，为什么？解答：n=0:1:15;x1n=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1/16)*n);xk1=fft(x1n);xk1=abs(xk1);subplot(1,2,1);stem(n,xk1);xlabel('k');ylabel('X(k)');legend('f=1/16');x2n=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1/64)*n);xk2=fft(x2n);xk2=abs(xk2);subplot(1,2,2);stem(n,xk2);xlabel('k');ylabel('X(k)');legend('f=1/64');n=0:1:127;x1n=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1/16)*n);xk1=fft(x1n);xk1=abs(xk1);stem(n,xk1);xlabel('k');ylabel('X(k)');legend('f=1/16');实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第10页共19页x2n=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1/64)*n);xk2=fft(x2n);xk2=abs(xk2);stem(n,xk2);xlabel('k');ylabel('X(k)');legend('f=1/64');分析：由于离散傅里叶变换的选频性质：()2/ojqwnoxnewN2()2()/1()[()]01jqkjqkNNqkeXkDFTxnqke当q不等于整数时，则信号频谱会发生泄漏。（5）用FFT分别计算()axn（8,2pq）和()bxn（0.1,0.0625af）的16点循环卷积和线性卷积。解答：n=0:1:15;xan=exp(-(n-8).^2/2);xbn=exp(-0.1*n).*sin(2*pi*0.0625*n);实验二快速傅里叶变换（FFT）及其应用闫春遐00824049第11页共19页subplot(4,1,1);stem(n,xan);xlabel('n');ylabel('xa(n)');subplot(4,1,2);stem(n,xbn);xlabel('n');ylabel('xb(n)');xak=fft(xan);xbk=fft(xbn);x1k=xak.*xbk;x1n=ifft(x1