數(shù)字信號(hào)處理實(shí)驗(yàn)二.doc
《數(shù)字信號(hào)處理實(shí)驗(yàn)二.doc》由會(huì)員分享，可在線閱讀，更多相關(guān)《數(shù)字信號(hào)處理實(shí)驗(yàn)二.doc（11頁珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
實(shí)驗(yàn)報(bào)告（本科） 學(xué) 號(hào) 2015141443002 姓 名 柏沖 專 業(yè) 通信工程 日 期 2017/12/4 實(shí)驗(yàn)題目 時(shí)域采樣和頻域采樣 一、實(shí)驗(yàn)?zāi)康? 時(shí)域采樣理論與頻域采樣理論是數(shù)字信號(hào)處理中重要的理論。要求掌握模擬信號(hào)采樣前后頻譜的變化，以及如何選擇采樣頻率才能使得采樣后的信號(hào)不丟失信息；要求掌握頻率采樣會(huì)引起時(shí)域周期化的概念，以及頻域采樣定理及其對(duì)頻域采樣點(diǎn)數(shù)選擇的指導(dǎo)作用。 2、 實(shí)驗(yàn)過程 附：源程序 (1)時(shí)域采樣 Tp=64/1000; %觀察時(shí)間Tp=64毫秒 %產(chǎn)生M長采樣序列x(n) Fs=1000; T=1/Fs; M=Tp*Fs; n=0:M-1; A=444.128; a=pi*50*2^0.5; omega=pi*50*2^0.5; xnt=A*exp(-a*n*T).*sin(omega*n*T); Xk=T*fft(xnt,M); %M點(diǎn)FFT[(xnt)] subplot(3,2,1); stem(xnt,.); %調(diào)用編繪圖函數(shù)stem繪制序列圖 box on;title((a) Fs=1000Hz); k=0:M-1;fk=k/Tp; subplot(3,2,2);stem(fk,abs(Xk),.);title((a) T*FT[xa(nT)],Fs=1000Hz); xlabel(f(Hz));ylabel(幅度);axis([0,Fs,0,1.2*max(abs(Xk))]); % Fs=300Hz和 Fs=200Hz的程序與上面Fs=1000Hz完全相同。 Tp=64/1000; %觀察時(shí)間Tp=64毫秒 %產(chǎn)生M長采樣序列x(n) Fs=300; T=1/Fs; M=Tp*Fs; n=0:M-1; A=444.128; a=pi*50*2^0.5; omega=pi*50*2^0.5; xnt=A*exp(-a*n*T).*sin(omega*n*T); M1=fix(M); Xk=T*fft(xnt,M1); %M點(diǎn)FFT[(xnt)] subplot(3,2,3); stem(xnt,.); %調(diào)用自編繪圖函數(shù)stem繪制序列圖 box on;title((b) Fs=300Hz); k=0:M-1;fk=k/Tp; subplot(3,2,4);stem(fk,abs(Xk),.);title((b) T*FT[xa(nT)],Fs=300Hz); xlabel(f(Hz));ylabel(幅度);axis([0,Fs,0,1.2*max(abs(Xk))]); Tp=64/1000; %觀察時(shí)間Tp=64毫秒 %產(chǎn)生M長采樣序列x(n) Fs=200; T=1/Fs; M=Tp*Fs; n=0:M-1; A=444.128; a=pi*50*2^0.5; omega=pi*50*2^0.5; xnt=A*exp(-a*n*T).*sin(omega*n*T); M2=fix(M); Xk=T*fft(xnt,M2); %M點(diǎn)FFT[(xnt)] yn=xa(nT);subplot(3,2,5); stem(xnt,.); %調(diào)用自編繪圖函數(shù)stem繪制序列圖 box on;title((c) Fs=200Hz); k=0:M-1;fk=k/Tp; subplot(3,2,6);stem(fk,abs(Xk),.);title((c) T*FT[xa(nT)],Fs=200Hz); xlabel(f(Hz));ylabel(幅度);axis([0,Fs,0,1.2*max(abs(Xk))]); 程序運(yùn)行結(jié)果 （2） 頻域采樣 M=27;N=32;n=0:M-1; %產(chǎn)生M長三角波序列x(n) xa=1:ceil(M/2); xb= fix(M/2):-1:1; xn=[xa,xb]; Xk=fft(xn,1024); %1024點(diǎn)FFT[x(n)], 用于近似序列x(n)的Ft X32k=fft(xn,32); %32點(diǎn)FFT[x(n)] x32n=ifft(X32k); %32點(diǎn)IFFT[X32(k)]得到x32(n) X16k=X32k(1:2:N); %隔點(diǎn)抽取X32k得到X16(K) x16n=ifft(X16k,N/2); %16點(diǎn)IFFT[X16(k)]得到x16(n) subplot(3,2,2);stem(n,xn,.);box on title((b) 三角波序列x(n));xlabel(n);ylabel(x(n));axis([0,32,0,20]) k=0:1023;wk=2*k/1024; % subplot(3,2,1);plot(wk,abs(Xk));title((a)FT[x(n)]); xlabel(\omega/\pi);ylabel(|X(e^j^\omega)|);axis([0,1,0,200]) k=0:N/2-1; subplot(3,2,3);stem(k,abs(X16k),.);box on title((c) 16點(diǎn)頻域采樣);xlabel(k);ylabel(|X_1_6(k)|);axis([0,8,0,200]) n1=0:N/2-1; subplot(3,2,4);stem(n1,x16n,.);box on;title((d)16點(diǎn)IDFT[X_1_6(k)]); xlabel(n);ylabel(x_1_6(n));axis([0,32,0,20]) k=0:N-1; subplot(3,2,5);stem(k,abs(X32k),.);box on title((e) 32點(diǎn)頻域采樣);xlabel(k);ylabel(|X_3_2(k)|);axis([0,16,0,200]) n1=0:N-1; subplot(3,2,6);stem(n1,x32n,.);box on title((f)32點(diǎn)IDFT[X_3_2(k)]); xlabel(n);ylabel(x_3_2(n));axis([0,32,0,20]) 程序運(yùn)行結(jié)果 3、 問題分析（思考題以及老師提問） （1）如果序列x(n)的長度為M，希望得到其頻譜X(ejω)在]2,0[p上的N點(diǎn)等間隔采樣，當(dāng)N下載提示(請(qǐng)認(rèn)真閱讀)
- 1.請(qǐng)仔細(xì)閱讀文檔，確保文檔完整性，對(duì)于不預(yù)覽、不比對(duì)內(nèi)容而直接下載帶來的問題本站不予受理。
- 2.下載的文檔，不會(huì)出現(xiàn)我們的網(wǎng)址水印。
- 3、該文檔所得收入（下載+內(nèi)容+預(yù)覽）歸上傳者、原創(chuàng)作者；如果您是本文檔原作者，請(qǐng)點(diǎn)此認(rèn)領(lǐng)！既往收益都?xì)w您。
文檔包含非法信息？點(diǎn)此舉報(bào)后獲取現(xiàn)金獎(jiǎng)勵(lì)！
下載文檔到電腦，查找使用更方便
9.9 積分
| 下載 |
還剩頁未讀，繼續(xù)閱讀
- 配套講稿：
如PPT文件的首頁顯示word圖標(biāo)，表示該P(yáng)PT已包含配套word講稿。雙擊word圖標(biāo)可打開word文檔。
- 特殊限制：
部分文檔作品中含有的國旗、國徽等圖片，僅作為作品整體效果示例展示，禁止商用。設(shè)計(jì)者僅對(duì)作品中獨(dú)創(chuàng)性部分享有著作權(quán)。
- 關(guān) 鍵 詞：
- 數(shù)字信號(hào) 處理 實(shí)驗(yàn)
裝配圖網(wǎng)所有資源均是用戶自行上傳分享，僅供網(wǎng)友學(xué)習(xí)交流，未經(jīng)上傳用戶書面授權(quán)，請(qǐng)勿作他用。
相關(guān)資源
更多
正為您匹配相似的精品文檔
相關(guān)搜索
鏈接地址：http://m.jqnhouse.com/p-6623406.html