资源描述：
《哈工大数字信号处理上机实验》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、彳丁社络药嫌力热fc摩厂•严严、HarbinInstituteofTechnology实验扳吿课程名称：数字信号处理实验题th用FFT作谱分析院系：电子与信息T程学院班级：姓名:学号：指导教师：实验时间：2013年11月哈尔滨工业大学1.实验目的与原理(1)进一步加深DFT算法原理和基木性质的理解(因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基木性质)。(2)熟悉FFT算法原理和FFT子程序的应用。(3)学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出
2、现的分析误差及其原因，以便在实际中正确应用FFT。2.实验内容与步骤(1)复习DFT的定义、性质和用DFT作谱分析的有关内容。(2)复习FFT算法原理与编程思想，并对照DIT-FFT运算流图和程序框图。(3)编制信号产生程序，并逐个进行谱分析(即画岀幅频特性)。下而给岀针对各信号的FFT变换区间N以及对连续信号x6(t)的采样频率fs,供实验吋参考。Xi(n),x2(n),x3(n),x4(n),x5(n):N=8,16x6(t):fs=64(Hz),N=16,32,64(n=0:l:69)X}
3、(n)=R4(n)0=0&n<=3).*1;subplot(311),ste
4、m(nlx1),xlabel(,n,),ylabel(,x_1(n),),set(gca?,xtick,,0:5);subplot(312),stem(0:7,abs(fft(x1,8))l,r,);set(gca,1xtick*,0:7,1xlinT,[07]),ti11e(1x_1(n)8点DFT*);subplol(313),stem(0:15?abs(fft(x1,16)),*g*);set(gca?1xtickr,0:15,1xlim1,[015]),ti11e(*x_1(n)16点D
5、F「);n=0:10;x2=(n>=0&n<=3).*(n+1)+(n>=4&n<=7).f(8-n);figure,subplot(311),stem(n,x2),xlabel(9n*),ylabel(*x_2(n)*);subplot(312),stem(0:7,abs(fft(x2,8)),'r•),ti11e('x_2(n)8点DFT');subplot(313),stem(0:15,abs(fft(x2,16)),*g*)ti11e(1x_2(n)16点DFT*),set(gca,1x
6、tick*,0:15);x3=(n>=0&n<=3).*(4-n)+(n>=4&n<=7).*(n-3);figure,subplot(311),stem(n,x3),xlabel(9n1),yIabel(9x_2(n)•);subplot(312),stem(0:7,abs(fft(x3,8))1,r*),title(*x_3(n)8点DF「);subplot(313),stem(0:15,abs(fft(x3,16)),1g1)ti11e(1x_3(n)16点DFT1),set(gca,1x
7、tick1,0:15);n=0:19;x4=(n>=0&n<=19).*cos(n*pi/4);figure,subplot(311)lstem(n,x4)lxlabel(ln,),ylabel(lx_4(n)1);subplot(312),stem(0:7,abs(fft(x4(1:8),8)+fft(x4(9:16),8)+fft(x4(17:20),8)))titledx_4(n)8点DFT1);subplot(313),stem(0:15,abs{fft(x4(1:16),16)+fft
8、(x4(17:20),16)))title(*x_4(n)16点DF「),set(gca,'xtick',0:15);x5=(n>=0&n<=19).*sin(n'pi/8);figure,subplot(311),stem(n,x5),xlabel(9n*),ylabel(1x_5(n)1)；subplot(312);stem(0:7,abs(fft(x5(1:8),8)+fft(x5(9:16),8)+fft(x5(17:20),8)))titledx_5(n)8点DFT);subplot(
