時域采樣和頻域采樣

時域采樣和頻域采樣

ID:27485557

大?。?66.50 KB

頁數(shù):8頁

時間:2018-12-04

時域采樣和頻域采樣_第1頁
時域采樣和頻域采樣_第2頁
時域采樣和頻域采樣_第3頁
時域采樣和頻域采樣_第4頁
時域采樣和頻域采樣_第5頁
資源描述:

《時域采樣和頻域采樣》由會員上傳分享,免費在線閱讀,更多相關內容在行業(yè)資料-天天文庫。

1、完美WORD格式.整理實驗二:時域采樣與頻域采樣一、實驗目的:時域采樣理論與頻域采樣理論是數(shù)字信號處理中的重要理論。要求掌握模擬信號采樣前后頻譜的變化,以及如何選擇采樣頻率才能使采樣后的信號不丟失信息;要求掌握頻率域采樣會引起時域周期化的概念,以及頻率域采樣定理及其對頻域采樣點數(shù)選擇的指導作用。二、實驗原理與方法:1、時域采樣定理的要點:1)對模擬信號以間隔T進行時域等間隔理想采樣,形成的采樣信號的頻譜是原模擬信號頻譜以采樣角頻率()為周期進行周期延拓。公式為:2)采樣頻率必須大于等于模擬信號最高頻率的兩倍以上,才能使采樣信號的

2、頻譜不產生頻譜混疊。利用計算機計算上式并不方便,下面我們導出另外一個公式,以便用計算機上進行實驗。理想采樣信號和模擬信號之間的關系為對上式進行傅立葉變換,得到:.專業(yè)資料分享.完美WORD格式.整理在上式的積分號內只有當時,才有非零值,因此上式中,在數(shù)值上=,再將代入,得到:上式的右邊就是序列的傅立葉變換,即上式說明理想采樣信號的傅立葉變換可用相應的采樣序列的傅立葉變換得到,只要將自變量ω用代替即可。2、頻域采樣定理的要點:a)對信號x(n)的頻譜函數(shù)X(ejω)在[0,2π]上等間隔采樣N點,得到則N點IDFT[]得到的序列就

3、是原序列x(n)以N為周期進行周期延拓后的主值區(qū)序列,公式為:b)由上式可知,頻域采樣點數(shù)N必須大于等于時域離散信號的長度M(即N≥M),才能使時域不產生混疊,則N點IDFT[]得到的序列就是原序列x(n),即=x(n)。如果N>M,比原序列尾部多N-M零點;如果N

4、到一個有用的結論,這兩個采樣理論具有對偶性:“時域采樣頻譜周期延拓,頻域采樣時域信號周期延拓”。因此放在一起進行實驗。三、實驗內容及步驟:1、時域采樣理論的驗證:給定模擬信號,式中A=444.128,=50π,=50πrad/s,它的幅頻特性曲線如圖2.1圖2.1的幅頻特性曲線現(xiàn)用DFT(FFT)求該模擬信號的幅頻特性,以驗證時域采樣理論。安照的幅頻特性曲線,選取三種采樣頻率,即=1kHz,300Hz,200Hz。觀測時間選。為使用DFT,首先用下面公式產生時域離散信號,對三種采樣頻率,采樣序列按順序用,,表示。因為采樣頻率不同

5、,得到的,,的長度不同,長度(點數(shù))用公式計算。選FFT的變換點數(shù)為M=64,序列長度不夠64的尾部加零。X(k)=FFT[x(n)],k=0,1,2,3,-----,M-1式中k代表的頻率為。要求:編寫實驗程序,計算、和.專業(yè)資料分享.完美WORD格式.整理的幅度特性,并繪圖顯示。觀察分析頻譜混疊失真。2、頻域采樣理論的驗證給定信號如下:編寫程序分別對頻譜函數(shù)在區(qū)間上等間隔采樣32和16點,得到:再分別對進行32點和16點IFFT,得到:分別畫出、的幅度譜,并繪圖顯示x(n)、的波形,進行對比和分析,驗證總結頻域采樣理論。提示

6、:頻域采樣用以下方法容易變程序實現(xiàn)。①直接調用MATLAB函數(shù)fft計算就得到在的32點頻率域采樣②抽取的偶數(shù)點即可得到在的16點頻率域采樣,即。當然也可以按照頻域采樣理論,先將信號x(n)以16為周期進行周期延拓,取其主值區(qū)(16點),再對其進行16點DFT(FFT),得到的就是在的16點頻率域采樣。.專業(yè)資料分享.完美WORD格式.整理四.實驗程序:實驗一:%========================closeall;clearall;clc;Tp=64/1000;%觀察時間Tp=64微秒%產生M長采樣序列x(n)%F

7、s=1000;T=1/Fs;Fs=1000;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=A*exp(-alph*n*T).*sin(omega*n*T);Xk=T*fft(xnt,M);%M點FFT[xnt)]yn='xa(nT)';subplot(3,2,1);tstem(xnt,yn);%調用自編繪圖函數(shù)tstem繪制序列圖boxon;title('(a)Fs=1000Hz');k=0:M-1;fk=k/Tp;subplot(

8、3,2,2);plot(fk,abs(Xk));title('(a)T*FT[xa(nT)],Fs=1000Hz');xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))])%===========

當前文檔最多預覽五頁,下載文檔查看全文

此文檔下載收益歸作者所有

當前文檔最多預覽五頁,下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學公式或PPT動畫的文件,查看預覽時可能會顯示錯亂或異常,文件下載后無此問題,請放心下載。
2. 本文檔由用戶上傳,版權歸屬用戶,天天文庫負責整理代發(fā)布。如果您對本文檔版權有爭議請及時聯(lián)系客服。
3. 下載前請仔細閱讀文檔內容,確認文檔內容符合您的需求后進行下載,若出現(xiàn)內容與標題不符可向本站投訴處理。
4. 下載文檔時可能由于網絡波動等原因無法下載或下載錯誤,付費完成后未能成功下載的用戶請聯(lián)系客服處理。