数字信号处理--matlab滤波器实验.docVIP

广 西 工 学 院 实 验 报 告 用 纸 实验名称 快速Fourier变换（FFT）及其应用 实验成绩 指导老师 陈艳 系（院） 计算机工程系 班级 通信042 学号 学生姓名 一、实验目的 1．在理论学习的基础上，通过本实验，加深对FFT的理解。 2．熟悉应用FFT对典型信号进行频谱分析的方法。 3．了解应用FFT进行信号频谱分析过程中可能出现的问题以便正确应用。 二、实验原理 快速离散Fourier变换在计算机上的实现，当序列x（n）的长度为N时，它的DFT定义为： 反变换为： 在运用DFT进行频谱分析的过程可能产生三种误差： 混叠 序列的频谱在被采样信号的周期延拓，当采样速率不满足Nyquist定理时，就会发生频谱混叠，使采样 后的信号序列不能真实反映原信号频谱。 泄露 我们往往用截短的序列来代替很长的序列，这样可以用较短的DFT来对信号进行频谱分析，这种截短等价于给原信号序列乘以一个矩形窗函数，所得的频谱是原序列频谱的扩展。 （3） 栅栏效应 DFT是对单位圆上Z变换的均匀采样，所以它不可能将频谱视为一个连续函数，从一定意义上看，用DFT 来观察频谱就好像通过一个栅栏来观看一个图景一样，只能在离散点看到真实的频谱。 三、实验内容 1、观察高斯序列的时域和幅频特性，固定信号xa（n）中的参数p=8,改变q的值，使q分别等于2，4，8，观察时域和幅频特性，了解q取不同值时的信号特性；固定q=8，改变p=8，13，14，观察信号特性，p=多少时，发生明显泄露现象？ 2、观察衰减正弦序列xb（n）的时域和幅频特性，a=0.1，f=0.0625，检查谱峰出现的位置是否正确，改变f，使f分别等于0.04375，0.5625，观察频谱形状和谱峰出现位置，有无混叠和泄露现象？ 3、观察三角波和反三角波序列和幅频特性，用N=8点FFT变换分析信号xc（n）和xd（n）的幅频特性，有何不同？在xc（n）和xd（n）末尾补零，用N=16点FFT分析频谱特性，观察发生了什么变化？还有相同之处吗？共 4 张 第 1 张 广 西 工 学 院 实 验 报 告 用 纸 实验名称 实验成绩 指导老师 系（院） 班级 学号 学生姓名4、一个连续信号含两个频率分量，经采样得x(n)=sin2π*0.125n+cos2π*(0.125+Δf)n n=0,1……,N-1已知N=16，Δf分别为1/16和1/64，观察其频谱；当N=128时，Δf不变，其结果有何不同？ ●实验程序如下： 1、 n=0:15; m=16; p=input(‘p=’); q=input(‘q=’); xa=exp(-(n-p).^2/q); w=linspace(-3*pi,3*pi,512); h=freqz(xa,[1],w).*exp(j*m*w); subplot(2,1,1);stem(n,xa); subplot(2,1,2);plot(w/pi,abs(h));title(‘magnitude response’); p=8,q=8p=8,q=4p=8,q=2 p=13,q=8p=14,q=8 2、 n=0:15; a=0.1; f=input(‘f=’); xb=exp(-a*n).*sin(2*pi*f*n); h=fft(xb); subplot(2,1,1);stem(n,xb); subplot(2,1,2);stem(n,abs(h)); 共 4 张 第 2 张 广 西 工 学 院 实 验 报 告 用 纸 实验名称 实验成绩 指导老师 系（院） 班级 学号 学生姓名 f=0.0625f=0.4375f=0.5625 3、 n=0:7; xc=[1 2 3 4 4 3 2 1]; xd=[4 3 2 1 1 2 3 4]; figure(1) hc=fft(xc); hd=fft(xd); subplot(2,1,1);stem(n,abs(hc));title(8 point fft transform); subplot(2,1,2);stem(n,abs(hd));title(8 point fft transform); figure(2) n1=0:15; xc1=[1 2 3 4 4 3 2 1 0 0 0 0 0 0 0 0]; xd1=[4 3 2 1 1 2 3 4 0 0 0 0 0 0 0 0]; hc1=fft(xc1); hd1=fft(xd1); subplot(2,1,1);stem(n1,abs(hc1));title(16 point fft transform); subplot(2,1,2);stem(n1,abs(hd1));title(