數(shù)字信號處理實驗二.doc
《數(shù)字信號處理實驗二.doc》由會員分享,可在線閱讀,更多相關(guān)《數(shù)字信號處理實驗二.doc(11頁珍藏版)》請在裝配圖網(wǎng)上搜索。
實驗報告(本科) 學(xué) 號 2015141443002 姓 名 柏沖 專 業(yè) 通信工程 日 期 2017/12/4 實驗題目 時域采樣和頻域采樣 一、實驗?zāi)康? 時域采樣理論與頻域采樣理論是數(shù)字信號處理中重要的理論。要求掌握模擬信號采樣前后頻譜的變化,以及如何選擇采樣頻率才能使得采樣后的信號不丟失信息;要求掌握頻率采樣會引起時域周期化的概念,以及頻域采樣定理及其對頻域采樣點數(shù)選擇的指導(dǎo)作用。 2、 實驗過程 附:源程序 (1)時域采樣 Tp=64/1000; %觀察時間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點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; %觀察時間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點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; %觀察時間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點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))]); 程序運行結(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點FFT[x(n)], 用于近似序列x(n)的Ft X32k=fft(xn,32); %32點FFT[x(n)] x32n=ifft(X32k); %32點IFFT[X32(k)]得到x32(n) X16k=X32k(1:2:N); %隔點抽取X32k得到X16(K) x16n=ifft(X16k,N/2); %16點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點頻域采樣);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點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點頻域采樣);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點IDFT[X_3_2(k)]); xlabel(n);ylabel(x_3_2(n));axis([0,32,0,20]) 程序運行結(jié)果 3、 問題分析(思考題以及老師提問) (1)如果序列x(n)的長度為M,希望得到其頻譜X(ejω)在]2,0[p上的N點等間隔采樣,當(dāng)N- 1.請仔細閱讀文檔,確保文檔完整性,對于不預(yù)覽、不比對內(nèi)容而直接下載帶來的問題本站不予受理。
- 2.下載的文檔,不會出現(xiàn)我們的網(wǎng)址水印。
- 3、該文檔所得收入(下載+內(nèi)容+預(yù)覽)歸上傳者、原創(chuàng)作者;如果您是本文檔原作者,請點此認(rèn)領(lǐng)!既往收益都歸您。
下載文檔到電腦,查找使用更方便
9.9 積分
下載 |
- 配套講稿:
如PPT文件的首頁顯示word圖標(biāo),表示該PPT已包含配套word講稿。雙擊word圖標(biāo)可打開word文檔。
- 特殊限制:
部分文檔作品中含有的國旗、國徽等圖片,僅作為作品整體效果示例展示,禁止商用。設(shè)計者僅對作品中獨創(chuàng)性部分享有著作權(quán)。
- 關(guān) 鍵 詞:
- 數(shù)字信號 處理 實驗
鏈接地址:http://www.hcyjhs8.com/p-6623406.html