课件数字信号处理实验三.docxVIP

实验三 快速傅里叶变换及应用实验目的在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉相关的MATLAB函数。应用FFT对典型信号进行频谱分析。了解应用FFT进行频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。应用FFT实现序列的线性卷积和相关。实验内容实验中用到的信号序列：高斯序列∶衰减正弦序列：三角波序列：反三角波序列：上机实验内容：观察高斯序列的时域和幅频特性，固定信号中参数 p=8，改变 q 的值，使 q 分别等于2，4，8，观察它们的时域和幅频特性，了解当 q 取不同值时，对信号序列的时域幅频特性的影响；固定 q=8，改变 p，使 p 分别等于 8，13，14，观察参数 p 变化对信号序列的时域及幅频特性的影响，注意 p 等于多少时，会发生明显的泄漏现象，混叠是否 也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。观察衰减正弦序列的时域和幅频特性，a=0.1，f=0.0625，检查谱峰出现位置是否正确，注意频谱的形状，绘出幅频特性曲线，改变 f，使 f 分别等于 0.4375 和 0.5625，观察这两种情况下，频谱的形状和谱峰出现位置，有无混叠和泄漏现象？说明产生现象的原因。观察三角波和反三角波序列的时域和幅频特性，用N=8点FFT分析信号序列：和的幅频特性，观察两者的序列形状和频谱曲线有什么异同？绘出两序列及其幅频特性曲 线。在和末尾补零，用 N=32 点FFT分析这两个信号的幅频特性，观察幅频特性发生了什么变化？两情况的FFT频谱还有相同之处吗？这些变化说明了什么？一个连续信号含两个频率分量，经采样得已知N=16，Δf分别为1/16和1/64，观察其频谱；当N=128时，Δf不变，其结果有何不同。用FFT分别计算(p=8,q=2)和(a=0.1,f=0.0625)的16点循环卷积和线性卷积。用FFT分别计算(p=8,q=2)和(a=0.1,f=0.0625)的16点循环相关和线性相关，问一共有多少种结果，他们之间有何异同点。用FFT分别计算(p=8,q=2)和(a=0.1,f=0.0625)的自相关函数。 各题实验程序与结果分析1、观察高斯序列的时域和幅频特性。实验程序：（1）clc;n=0:1:15;% p=8,q=2p=8;q=2;x1=exp(-(n-p).^2./q);fp1=fft(x1);fp1=abs(fp1);subplot(3,2,1);stem(n,x1);xlabel(n);ylabel(时域特性);title(p=8,q=2);subplot(3,2,2);plot(n,fp1);xlabel(n);ylabel(幅频特性);title(p=8,q=2);% p=8,q=4p=8;q=4;x2=exp(-(n-p).^2./q);fp2=fft(x2);fp2=abs(fp2);subplot(3,2,3);stem(n,x2);xlabel(n);ylabel(时域特性);title(p=8,q=4);subplot(3,2,4);plot(n,fp2);xlabel(n);ylabel(幅频特性);title(p=8,q=4);% p=8,q=8p=8;q=8;x3=exp(-(n-p).^2./q);fp3=fft(x3);fp3=abs(fp3);subplot(3,2,5);stem(n,x3);xlabel(n);ylabel(时域特性);title(p=8,q=8);subplot(3,2,6);plot(n,fp3);xlabel(n);ylabel(幅频特性);title(p=8,q=8);（2）clc;n=0:1:15;% p=8,q=8p=8;q=8;x3=exp(-(n-p).^2./q);fp3=fft(x3);fp3=abs(fp3);subplot(3,2,1);stem(n,x3);xlabel(n);ylabel(时域特性);title(p=8,q=8);subplot(3,2,2);plot(n,fp3);xlabel(n);ylabel(幅频特性);title(p=8,q=8);% p=13,q=8p=13;q=8;x4=exp(-(n-p).^2./q);fp4=fft(x4);fp4=abs(fp4);subplot(3,2,3);stem(n,x4);xlabel(n);ylabel(时域特性);title(p=13,q=8);subplot(3,2,4);plot(n,fp4);xlabel(n);ylabel(幅频特性);title(p=13,q=8);% p=14,q=8p=14;q=8;x5=exp(-(n-p).^2./q);fp5=fft(x5);fp5=abs(fp5);subplot(3,2,