現(xiàn)代信號處理,信號短時傅里葉變換分析作業(yè)_第1頁
現(xiàn)代信號處理,信號短時傅里葉變換分析作業(yè)_第2頁
現(xiàn)代信號處理,信號短時傅里葉變換分析作業(yè)_第3頁
現(xiàn)代信號處理,信號短時傅里葉變換分析作業(yè)_第4頁
現(xiàn)代信號處理,信號短時傅里葉變換分析作業(yè)_第5頁
已閱讀5頁，還剩2頁未讀， 繼續(xù)免費閱讀

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

1、現(xiàn)代信號處理第一次作業(yè)令x(n)=5ej(0.15n2)+ 6ej(300n-0.15n2)，w(n)為高斯窗函數(shù)。試用matlab軟件，取不同長度的窗函數(shù)，分別求x(n)離散短時傅里葉變換，并進行信號重構。試討論窗函數(shù)長度對時頻分辨率、重構精度的影響。解：matlab代碼見附件。1. 當取高斯窗函數(shù)的長度為7時短時傅里葉變換的時頻圖如下：可以看出，時間分辨率比較高，但頻率分辨率比較低。在此窗函數(shù)作用下，重構圖像與原始圖像的差圖如下：可以看出，兩幅圖最大差值接近4*10-15。2.將高斯窗函數(shù)的長度設為255時所得時頻圖如下：可以看出在窗函數(shù)長度比較大的時候，時間分辨率比較低，但頻率分辨率比較

2、高。此時，進行信號重構，重構后的圖像與原始圖像的差圖如下：此時兩者的最大差值接近2.5*10-14，顯然比窗函數(shù)為7時的要大。所以，stft的時間分辨率由窗函數(shù)的寬度決定 ，stft的頻率分辨率由窗函數(shù)的頻譜寬度決定 。時間分辨率與頻率分辨率的乘積為一個常量。窗函數(shù)寬度越長，時間分辨率越低，頻率分辨率越高，信號重構精度越差。代碼：n=256.75;n=1:0.25:n;m=length(n);%產生線性調頻信號x=5*exp(j*0.15*n.*n)+6*exp(j*(300*n-0.15*n.*n);figure(1);plot(n,x);figure(1);plot(n,x);figure(2);plot(n,x);xlabel(時間 n);ylabel(幅值 a);% 產生高斯窗函數(shù)w = window(gausswin,7);% 計算短時傅里葉變換y=x;tfrstft(y,1:m,m,w);tfr,t,f=tfrstft(y,1:m,m,w);% 進行短時傅里葉反變換g=tfristft(tfr,1:m,w);% 畫出重構信號與原始信號的差圖figure(3)

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間，僅對用戶上傳內容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論