哈工大数字信号处理上机实验.docVIP

Harbin Institute of Technology 实验报告 课程名称： 数字信号处理I 实验题目：一：利用FFT对信号进行频谱分析 二：窗函数法设计FIR数字滤波器 院 系： 班 级： 姓 名： 学 号： 指导教师： 实验时间： 2015 年 6 月 哈尔滨工业大学 实验一：利用FFT对信号进行频谱分析 1.实验目的 (1) 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质)。 (2) 熟悉FFT算法原理和FFT子程序的应用。 (3) 学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 2.实验步骤 (1) 复习DFT的定义、 性质和用DFT作谱分析的有关内容。 (2) 复习FFT算法原理与编程思想， 并对照DIT-FFT运算流图和程序框图。 (3) 编制信号产生子程序， 产生以下典型信号供谱分析用： (4) 编写主程序。 (5) 按实验内容要求， 上机实验， 并写出实验报告。 3.实验程序与结果 （1） 程序n=0:7; x=(n0).*0+(n3).*0+((n=0)(n=3)).*1; N1=8; n1=0:N1-1; y1=fft(x,N1); subplot(2,2,1);stem(n,x); title(x_1); subplot(2,2,2); stem(n1,abs(y1)); xlabel(k); ylabel(|X_1(k)|); title(x_1的8点FFT); N2=16; n2=0:N2-1; y2=fft(x,N2); subplot(2,2,3); stem(n2,abs(y2)); title(x_1的16点FFT); xlabel(k); ylabel(|X_1(k)|); 结果如图： 实验误差分析：FFT频谱为一个抽样信号，由图形看出，随着采样率的提高，得到的FFT频谱分辨率就越高，当N趋于无限时频谱包络接近理论的抽样函数。 程序： clear all N=8; n=0:N-1; x=(n0)*0+(n7)*0+((n=0)(n=3)).*(n+1)+((n=4)(n=7)).*(8-n); subplot(2,2,1); stem(n,x); title(x_2); xlabel(n); ylabel(x_2); y=fft(x,N); subplot(2,2,2); stem(n,abs(y)) xlabel(k); ylabel(|X_2(k)|); title(x_2(n)的8点FFT); N2=16; n2=0:15; y2=fft(x,N2); subplot(2,2,3); stem(n2,abs(y2)); xlabel(k); ylabel(|X_2(k)|); title(x_2(n)的8点FFT); 结果： 实验误差分析：?理论FFT频谱为一个抽样函数平方的信号，由图形看出，随着采样率的提高，得到的FFT频谱分辨率就越高，当N趋于无限时频谱包络接近理论的抽样函数。 程序： clear all N=8; n=0:N-1; x=(n0)*0+(n7)*0+((n=0)(n=3)).*(4-n)+((n=4)(n=7)).*(n-3); subplot(2,2,1); stem(n,x); title(x_3); xlabel(n); ylabel(x_3); y=fft(x,N); subplot(2,2,2); stem(n,abs(y)) xlabel(k); ylabel(|X_3(k)|); title(x_3(n)的8点FFT); N2=16; n2=0:15; y2=fft(x,N2); subplot(2,2,3); stem(n2,abs(y2)); xlabel(k); ylabel(|X_3(k)|); title(x_3(n)的16点FFT); 结果： 实验误差分析： 与的8点DFT的模相等，因为。但是，当N=16时，与不满足循环移位关系，此时 与的16点DFT的模不相等。 程序： clear all N=8; n=0:19; x=(n0)*0+(n19)*0+((n=0