哈工大數(shù)字信號處理實驗一源程序

實驗?zāi)康呐c原理(1)進一步加深DFT算法原理和基本性質(zhì)的理解(因為FFT只是DFT的一種快速算法，所以FFT的運算結(jié)果必然滿足DFT的基本性質(zhì))。(2)熟悉FFT算法原理和FFT子程序的應(yīng)用。(3)學習用FFT對連續(xù)信號和時域離散信號進行譜分析的方法，了解可能出現(xiàn)的分析誤差及其原因，以便在實際中正確應(yīng)用FFT。實驗內(nèi)容與步驟(1)復習DFT的定義、性質(zhì)和用DFT作譜分析的有關(guān)內(nèi)容。(2)復習FFT算法原理與編程思想，并對照DIT-FFT運算流圖和程序框圖。(3)編制信號產(chǎn)生程序，并逐個進行譜分析（即畫出幅頻特性）。下面給出針對各信號的FFT變換區(qū)間N以及對連續(xù)信號x6(t)的采樣頻率fs，供實驗時參考。x1(n),x2(n),x3(n),x4(n),x5(n):N=8,16x6(t):fs=64(Hz),N=16,32,64（n=0:1:69）(4)令x(n)=x4(n)+x5(n),用FFT計算8點和16點離散傅里葉變換，X(k)=DFT［x(n)］(5)令x(n)=x4(n)+jx5(n)，重復(2)。實驗原程序n1=0:3;x1=[1,1,1,1];%產(chǎn)生信號figure,subplot(311),stem(n1,x1),xlabel('n'),ylabel('x_1(n)'),title('原始序列');subplot(312),stem(0:7,abs(fft(x1,8)),'r'),title('x_1(n)8點DFT');%8點FFTsubplot(313),stem(0:15,abs(fft(x1,16)),'g'),title('x_1(n)16點DFT');%16點FFTn2=0:7,x2=[1,2,3,4,4,3,2,1];%產(chǎn)生信號figure,subplot(311),stem(n2,x2),xlabel('n'),ylabel('x_2(n)'),title('原始序列');subplot(312),stem(0:7,abs(fft(x2,8)),'r'),title('x_1(n)8點DFT');%8點FFTsubplot(313),stem(0:15,abs(fft(x2,16)),'g'),title('x_1(n)16點DFT');%16點FFTn3=0:7,x3=[4,3,2,1,1,2,3,4];%產(chǎn)生信號figure,subplot(311),stem(n3,x3),xlabel('n'),ylabel('x_3(n)'),title('原始序列');subplot(312),stem(0:7,abs(fft(x3,8)),'r'),title('x_1(n)8點DFT');%8點FFTsubplot(313),stem(0:15,abs(fft(x3,16)),'g'),title('x_1(n)16點DFT');%16點FFTn4=0:19;x4=cos(n4*pi/4);figure,subplot(311),stem(n4,x4),xlabel('n'),ylabel('x_4(n)'),title('原始序列');subplot(312),stem(0:7,abs(fft(x4,8)),'r'),title('x_4(n)8點DFT');%8點FFTsubplot(313),stem(0:15,abs(fft(x4,16)),'g'),title('x_4(n)16點DFT');%16點FFTn5=0:19;x5=sin(n5*pi/8);%產(chǎn)生信號figure,subplot(311),stem(n5,x5),xlabel('n'),ylabel('x_5(n)'),title('原始序列');subplot(312),stem(0:7,abs(fft(x5,8)),'r'),title('x_5(n)8點DFT');%8點FFTsubplot(313),stem(0:15,abs(fft(x5,16)),'g'),title('x_5(n)16點DFT');%16點FFTn6=0:69;x6=cos(pi/8*n6)+cos(pi/4*n6)+cos(pi/16*5*n6);%產(chǎn)生信號figure,subplot(221);stem(n6,abs(x6)),xlabel('n'),ylabel('x_6(n)'),title('原始序列');subplot(222),stem(0:15,abs(fft(x6,16)),'r'),title('x_6(n)16點DFT');%16點FFTsubplot(223),stem(0:31,abs(fft(x6,32)),'r'),title('x_6(n)32點DFT');%32點FFTsubplot(224),stem(0:63,abs(fft(x6,64)),'r'),title('x_6(n)64點DFT');%64點FFTn7=0:19;x=x4+x5;%產(chǎn)生信號figure,subplot(311),stem(n7,x),xlabel('n'),ylabel('x'),title('原始序列');subplot(312),stem(0:7,abs(fft(x,8)),'r'),title('x8點DFT');%8點FFTsubplot(313),stem(0:15,abs(fft(x,16)),'g'),title('x16點DFT');%16點aFFTn8=0:19;xk=x4+j*x5;%產(chǎn)生信號figure,subplot(311),stem(n8,x),xlabel('n'),ylabel('xk'),title('原始序列');subplot(31