【正文】 ) plot((0:N21)/N2*fs,abs(X2),39。 vct_w1=0: 2*pi/N1: 2*pi2*pi/N1。 %w為頻率數(shù)組 DFT/DTFT仿真 clear all clf j=sqrt(1)。L=20039。f1=500。 yL=dtft(xL,w)。 由物理分辨率計算公式 Lfscf ???可知 ， 在采樣頻率 fs 一定時， 增加時窗寬度 L，兩種類型的時窗的物理分辨率均提高。 w1=2*pi*f1/fs。 plot(w/2/pi*fs,abs(yL),39。 代碼為： clear。 fs=5000。得出相應結論，給出提高物理分辨率的方法。 5. 分組討論，獨立撰寫項目報告，提交 Word 電子文檔。 w=0:2*pi/300:2*pi。f1=1000。L=10039。 wh=*cos(2*pi*n/(L1))。 L=100。) gtext(39。w3=2*pi*f3/fs。我們可以通過 增加采樣的長度 或者 加 海明 窗 來 降低頻譜泄露量 。 w1=2*pi*f1/fs。 X2=dtft(x,vct_w2)。DFT, L=10 N=6439。 n=0:L21。DFT, L=100 N=1639。b:39。,39。 subplot(3,1,1) axis([0 fs/2 0 60]) hold on [NN HN1]=size(abs(X1))。DTFT39。DTFT39。 x=sin(w1*n)+sin(w2*n)+sin(w3*n)。*w)。) gtext(39。 fs=5000。 w=0:2*pi/500:2*pi。 仿真 結果 如下 ： 仿真結果分析： 由物理分辨率計算公式知 ： Lfscf ??? ，（ c=2）則當 L=100， fs=5kHz 時 ?f? =100 Hz， 又因為 3 個信號的譜峰間隔 f? 為 50Hz，則 ?f? f? .物理分辨率大于信號頻率間隔，故不能分辨出 3 個信號的頻譜； 當 L=300， fs=5kHz 時， ?f? =，則 ?f? f? 物理分辨率小于信號頻率間隔，故能分辨出 3 個信號的頻譜 結論 ：取相同的時窗寬度 L，則矩形窗的物理分辨率比海明窗高。f3=1100。 yL=dtft(xL_H,w)。 代碼 ：（略） 仿真結果如下： 仿真結果分析： 由 物理分辨率計算公式知： Lfscf ??? （ c=1） ，則當 L=64， fs=5kHz 時 ?f? = Hz， 又因為 3 個信號的譜峰間隔 f? 為 50Hz，則 ?f? f? .物理分辨率大于 信號頻率間隔 ，故不能分辨出 3 個信號 的頻譜； 當 L=200， fs=5kHz 時， ?f? =25Hz，則 ?f? f? 物理分辨率小于信號頻率間隔，故能分辨出 3 個信號的頻譜 Hamming window ① 當時窗寬度 L=100時，設定 f1=1kHz， f2= kHz， f3= kHz； fs=5 kHz。f3=1100。 二 項目內(nèi)容 1. 舉例說明時窗寬度、類型與物理分辨率的關系。 三 項目報告正文 1. 數(shù)字譜分析的意義 信號特征提取，語音、圖象處理，通