数字信号处理实验报告2j.docVIP

南京信息工程大学 实验（实习）报告 实验（实习）名称 实验二、用FFT作普分析 日期 2011年5月16日 得分 指导老师 行鸿彦 系 电子与信息工程 专业 电子信息工程 年级 班次 姓名 学号 一、实验目的 （1）进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质)。 （2）熟悉FFT算法原理和FFT子程序的应用。 （3）学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 二、实验内容 (1) 复习DFT的定义、 性质和用DFT作谱分析的有关内容。 (2) 复习FFT算法原理与编程思想， 并对照DIT-FFT运算流图和程序框图， 分析和编写FFT子程序。 (3) 产生指数序列、正弦序列、单位脉冲序列、矩形序列、三角序列等并绘图说明。 1）并对所给出的信号逐个进行谱分析，用FFT计算 8 点和 16 点离散傅里叶变换进行计算。 ? 2）对任意两个信号进行叠加， 进行16 点离散傅里叶变换计算，并观察信号频谱的变化。 （4）编写不同频率正弦信号的叠加，并进行频谱分析，观察信号不同参数变化时，FFT的计算变化，并绘图说明。 （5）编写正弦信号与带噪信号叠加，并进行频谱分析，观察信号和噪声不同参数变化时，FFT的计算变化，并绘图说明。 （6）结合实验中所分析的典型序列幅频特性曲线，与理论结果比较， 并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法。 （7）验证卷积定理 三、实验结果 n=0:9:100; x1=exp(-0.1*n); subplot(521);plot(x1); subplot(522);stem(fft(x1,8)); x2=sin(n); subplot(523);plot(x2); subplot(524);stem(fft(x2,8)); x3 = [1,zeros(1,15)]; subplot(525); plot(x3); subplot(526);stem(fft(x3,8)); x4 = ones(1,16); subplot(527);plot(x4); subplot(528);stem(fft(x4,8)); for i = 1:3 x(i) = i + 1; end for i = 4:7 x(i) = 8 - i; end for i = 8:15 x(i) = 0; end x5 = [1,x]; subplot(529);plot(x5); subplot(5,2,10);stem(fft(x5,8)); n=0:9:100; x1=exp(-0.1*n); subplot(521);plot(x1); subplot(522);stem(fft(x1,16)); x2=sin(n); subplot(523);plot(x2); subplot(524);stem(fft(x2,16)); x3 = [1,zeros(1,15)]; subplot(525); plot(x3); subplot(526);stem(fft(x3,16)); x4 = ones(1,16); subplot(527);plot(x4); subplot(528);stem(fft(x4,16)); for i = 1:3 x(i) = i + 1; end for i = 4:7 x(i) = 8 - i; end for i = 8:15 x(i) = 0; end x5 = [1,x]; subplot(529);plot(x5); subplot(5,2,10);stem(fft(x5,16)); n=0:9:100; x1=exp(-0.1*n); x2=sin(n); x3=exp(-0.1*n)+sin(n); plot(x3); stem(fft(x3,16)); n=0:9:100; x1=sin(n); x2=sin(2*n); x=x1+x2; subplot(231);stem(x1); subplot(232);stem(x2); subplot(233);stem(x); subplot(234);stem(fft(x1)); subplot(235);stem(fft(x2)); subplot(236