哈工大數(shù)字信號(hào)處理實(shí)驗(yàn)一源程序

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