数字信号mtlab实验快速傅里叶变换及其应用.docVIP

实验三 快速傅里叶变换及其应用 上机实验内容 一、观察高斯序列的时域和幅频特性，固定信号xa(n)中参数p=8，改变q的值，分别等于2,4,8，观察他们的时域和幅频特性，了解党q取不同的值时，对信号序列的时域和幅频特性的影响；固定q=8，改变p,使p分别等于8,13,14，观察参数p变化对信号序列的时域及幅频特性的影响，注意p等于多少时，会发生明显的泄露现象，混叠是否也随之出现？ （1）P=8,q=2,4,8时， 程序如下： n=0:15;N=16;p=8;q1=2;q2=4;q3=8; figure(1) x1n=exp(-((n-p1).^2/q)); subplot(3,2,1),stem(n,x1n,.) x2n=exp(-((n-p2).^2/q)); subplot(3,2,3),stem(n,x2n,.) x3n=exp(-((n-p3).^2/q)); subplot(3,2,5),stem(n,x3n,.) x1k=fft(x1n,N); subplot(3,2,2),stem(n,x1k,.) x2k=fft(x2n,N); subplot(3,2,4),stem(n,x2k,.) x3k=fft(x3n,N); subplot(3,2,6),stem(n,x3k,.) 波形如下： p固定时，当p从2变化到8时，时域波形变化缓慢，低频分量增加，频谱泄露和混叠减小。 （2）q=8,p=8,13,14时，程序类同上 波形如下： 固定q=8时，当p从8变化到14，窗口位置偏移，受窗口宽度的影响，波形在高处截断产生严重的频谱泄露。 二、观察衰减正弦序列xb(n)的时域和幅频特性，a=0.1，f=0.0625,检查谱峰出现位置是否正确，注意频谱的形状，会出幅频特性曲线，改变f，使f分别等于0.4375,0.5625,观察着两种情况下，频谱的形状和谱峰出现的位置，有无混叠和泄露现象，说明产生的原因。 程序如下： clear ;close all n=0:15; N=16; a=0.1; f1=0.0625; f2=0.4375; f3=0.5625; x1n=exp(-a.*n).*sin(2*pi*f1*n); subplot(3,2,1),stem(n,x1n,.) x2n=exp(-a.*n).*sin(2*pi*f2*n); subplot(3,2,3),stem(n,x2n,.) x3n=exp(-a.*n).*sin(2*pi*f3*n); subplot(3,2,5),stem(n,x3n,.) x1k=fft(x1n,N); subplot(3,2,2),plot(n,x1k) x2k=fft(x2n,N); subplot(3,2,4),plot(n,x2k) x3k=fft(x3n,N); subplot(3,2,6),plot(n,x3k) 波形如下： 由于DFT的选频性，衰减正弦信号序列的频谱在信号频率，镜像频率和混叠频率点上呈现出峰值。当频率f=0.0625和f=0.4375时，满足采样定理，当频率f=0.5625时，不满足采样定理，频谱发生混叠。 三、观察三角波和反三角波序列的时域和频域特性，用N=8点FFT分析信号序列xc(n)和xd(n)的幅频特性，观察两者的序列性状和频谱曲线有什么异同？在两个序列后面补零，用N=32点FFT分析这两个信号，观察频谱发生了什么变化？两种情况下的FFT频谱还有相同之处么？这些变化说明了什么？ 程序如下： clear,close all ,n = 0:7;x = zeros(size(n));N=8; mask = (n = 0) (n = 3); x1(mask) = n(mask) ; mask = (n =4) (n=7); x1(mask) = 8 - n(mask) mask = (n = 0) (n = 3); x2(mask) = 4-n(mask) ; mask = (n =4) (n=7); x2(mask) = n(mask)-4 subplot(2,2,1),stem(n,x1,.) subplot(2,2,3),stem(n,x2,.) x1k=fft(x1(mask),N); x2k=fft(x2(mask),N); subplot(2,2,2),stem(n,x1k,.) subplot(2,2,4),stem(n,x2k,.) 波形如下：当作N=8点的DFT时 N=32时的波形： 当对8点的序列作8点的DTF变化时，有周期性易知，幅频图像应该是相同的，但是我的程序没有实现这一点，不清楚到底是哪里的问题。当做32点DTF时，幅频特性的曲线看不清楚。 四、一个连续信号含两个频率分量，经采样得x(n)=sin[2*pi*0,