数字信号处理第二版(吴镇扬)第三,四次实验(山师) .pdfVIP

实验三快速Fourier变换(FFT)及其应用

一、实验目的

1．在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB

中的有关函数。

2．应用FFT对典型信号进行频谱分析。

3.了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正

确应用FFT。

4．应用FFT实现两个序列的线性卷积和方法。

二、实验原理与方法

上机实验内容：

(1)、观察高斯序列的时域和幅频特性，固定信号xa(n)中参数p=8，改变q的

值，使q分别等于2，4，8，观察它们的时域和幅频特性，了解当q取不同值时，

对信号序列的时域幅频特性的影响；固定q=8，改变p，使p分别等于8，13，

14，观察参数p变化对信号序列的时域及幅频特性的影响，观察p等于多少时，

会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出

相应的时域序列和幅频特性曲线。

程序：

function[x,F]=gauss(p,q);

n=0:15;

x(n+1)=exp(-(n+1-p).^2/q);

F=fft(x);

end

clearall;

figure(1)

[x1,F1]=gauss(8,2);

n=0:15;

subplot(3,2,1);

plot(n,x1);text(6,0.2,p=8,q=2);

gridon;

xlabel(n);

ylabel(时域);

subplot(3,2,2);

plot(abs(F1));text(7.5,2,p=8,q=2);

gridon;

xlabel(k);

ylabel(频域);

[x2,F2]=gauss(8,4);

n=0:15;

subplot(3,2,3);

plot(n,x2);text(6,0.2,p=8,q=4);

gridon;

xlabel(n);

ylabel(时域);

subplot(3,2,4);

plot(abs(F2));text(7.5,2,p=8,q=4);

gridon;

xlabel(k);

ylabel(频域);

[x3,F3]=gauss(8,8);

n=0:15;

subplot(3,2,5);

plot(n,x3);text(6,0.2,p=8,q=8);

gridon;

xlabel(n);

ylabel(时域);

subplot(3,2,6);

plot(abs(F3));text(7.5,2,p=8,q=8);

gridon;

xlabel(k);

ylabel(频域);

figure(2)

[x4,F4]=gauss(8,8);

n=0:15;

subplot(3,2,1);

plot(n,x4);text(6,0.2,p=8,q=8);

gridon;

xlabel(n);

ylabel(时域);

subplot(3,2,2);

plot(abs(F4));text(7.5,2,p=8,q=8);

gridon;

xlabel(k);

ylabel(频域);

[x5,F5]=gauss(13,8);

n=0:15;

subplot(3,2,3);

plot(n,x5);text(6,0.2,p=13,q=8);

gridon;

xlabel(n);

ylabel(时域);

subplot(3,2,4);

plot(abs(F5));text(7.5,2,p=13,q=8);

gridon;

xlabel(k);

ylabel(频域);

[x6,F6]=gauss(14,8);

n=0:15;

subplot(3,2,5);

plot(n,x6);text(6,0.2,p=14,q=8);

gridon;

xlabel(n);

ylabe