資源描述:
《哈工大數(shù)字信號(hào)的處理實(shí)驗(yàn)報(bào)告材料》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、實(shí)用標(biāo)準(zhǔn)文案HarbinInstituteofTechnology實(shí)驗(yàn)報(bào)告課程名稱:數(shù)字信號(hào)處理院系:電子與信息工程學(xué)院班級:姓名:學(xué)號(hào):指導(dǎo)教師:冀振元實(shí)驗(yàn)時(shí)間:2014年11月哈爾濱工業(yè)大學(xué)精彩文檔實(shí)用標(biāo)準(zhǔn)文案實(shí)驗(yàn)一:用FFT作譜分析實(shí)驗(yàn)?zāi)康模?.進(jìn)一步加深DFT算法原理和基本性質(zhì)的理解(因?yàn)镕FT只是DFT的一種快速算法,所以FFT的運(yùn)算結(jié)果必然滿足DFT的基本性質(zhì))。2.熟悉FFT算法原理和FFT子程序的應(yīng)用。3.學(xué)習(xí)用FFT對連續(xù)信號(hào)和時(shí)域離散信號(hào)進(jìn)行譜分析的方法,了解可能出現(xiàn)的分析誤差及其原因
2、,以便在實(shí)際中正確應(yīng)用FFT。實(shí)驗(yàn)步驟(1)復(fù)習(xí)DFT的定義、性質(zhì)和用DFT作譜分析的有關(guān)內(nèi)容。(2)復(fù)習(xí)FFT算法原理與編程思想,并對照DIT-FFT運(yùn)算流圖和程序框圖,讀懂本實(shí)驗(yàn)提供的FFT子程序。(3)編制信號(hào)產(chǎn)生子程序,產(chǎn)生以下典型信號(hào)供譜分析用:精彩文檔實(shí)用標(biāo)準(zhǔn)文案(4)編寫主程序。(5)按實(shí)驗(yàn)內(nèi)容要求,上機(jī)實(shí)驗(yàn),并寫出實(shí)驗(yàn)報(bào)告。實(shí)驗(yàn)程序與結(jié)果:實(shí)驗(yàn)程序:%%整個(gè)程序通過switch語句切換得到x1~x6的波形產(chǎn)生和FFT變換,波形。clc;clear;n0=-10:0.1:20;n1=1:4;
3、x1=[1111];n2=0:7;x2=(n2>=0&n2<=3).*(n2+1)+(n2>=4&n2<=7).*(8-n2);n3=0:7;x3=(n3>=0&n3<=3).*(4-n3)+(n3>=4&n3<=7).*(n3-3);n4=0:19;x4=cos(pi/4*n4);n5=0:19;x5=sin(pi/8*n5);switchv精彩文檔實(shí)用標(biāo)準(zhǔn)文案case1subplot(311);stem(n1,x1);xlabel('n');ylabel('x1(n)');title('原時(shí)間序列')
4、;N=8;n1=0:N-1;X=fft(x1,N);X=abs(X);subplot(312);stem(n1,fftshift(X));title('N=8');N=16;n1=0:N-1;X=fft(x1,N);X=abs(X);subplot(313);stem(n1,fftshift(X));title('N=16');case2subplot(311);stem(n2,x2);xlabel('n');ylabel('x2(n)');title('原時(shí)間序列');N=8;n2=0:N-1;X=ff
5、t(x2,N);X=abs(X);subplot(312);stem(n2,fftshift(X));title('N=8');N=16;n2=0:N-1;X=fft(x2,N);X=abs(X);subplot(313);精彩文檔實(shí)用標(biāo)準(zhǔn)文案stem(n2,fftshift(X));title('N=16');case3subplot(311);stem(n3,x3);xlabel('n');ylabel('x3(n)');title('原時(shí)間序列');N=8;n3=0:N-1;X=fft(x3,N);
6、X=abs(X);subplot(312);stem(n3,fftshift(X));title('N=8');N=16;n3=0:N-1;X=fft(x3,N);X=abs(X);subplot(313);stem(n3,fftshift(X));title('N=16');case4subplot(311);stem(n4,x4);xlabel('n');ylabel('x4(n)');title('原時(shí)間序列');subplot(312);N=8;n4=0:N-1;X=fft(x4,N);X=abs
7、(X);stem(n4,fftshift(X));title('N=8');精彩文檔實(shí)用標(biāo)準(zhǔn)文案N=16;n4=0:N-1;X=fft(x4,N);X=abs(X);subplot(313);stem(n4,fftshift(X));title('N=16');case5subplot(311);stem(n5,x5);xlabel('n');ylabel('x5(n)');title('原時(shí)間序列');subplot(312);N=8;n5=0:N-1;X=fft(x5,N);X=abs(X);stem
8、(n5,fftshift(X));title('N=8');N=16;n5=0:N-1;X=fft(x5,N);X=abs(X);subplot(313);stem(n5,fftshift(X));title('N=16');case6t=0:0.1:69;x6=cos(8*pi*t)+cos(16*pi*t)+cos(20*pi*t);subplot(411);plot(t,x6);xlabel('n');ylabel('x