電子信息工程-數字譜分析實踐及誤差討論-文庫吧
2024-12-24 10:07 本頁面
【正文】 (w3*n)。 w=0:2*pi/300:2*pi。 yL=dtft(xL,w)。 plot(w/2/pi*fs,abs(yL),39。k39。) axis([500,2500,0,L/]) xlabel(39。f/Hz39。) gtext(39。L=6439。)。 仿真結果如下 ： ② 當時窗寬度 L=200時， 設定 f1=1kHz， f2= kHz， f3= kHz； fs=5 kHz。 代碼 ：（略） 仿真結果如下： 仿真結果分析： 由 物理分辨率計算公式知： 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。 代碼為： clear。 close all。 L=100。 n=0:L1。 fs=5000。f1=1000。f2=1050。f3=1100。 w1=2*pi*f1/fs。w2=2*pi*f2/fs。w3=2*pi*f3/fs。 xL=sin(w1*n)+sin(w2*n)+sin(w3*n)。 w=0:2*pi/500:2*pi。 wh=*cos(2*pi*n/(L1))。 xL_H=xL.*wh。 yL=dtft(xL_H,w)。 plot(w/2/pi*fs,abs(yL),39。k39。) axis([500,1500,0,L/3]) xlabel(39。f/Hz39。) gtext(39。L=10039。)。 仿真 結果 如下 ： ② 當時窗寬度 L=300 時，設定 f1=1kHz， f2= kHz， f3= kHz； fs=5 kHz。 代碼： clear。 close all。 L=300。 n=0:L1。 fs=5000。f1=1000。f2=1050。f3=1100。 w1=2*pi*f1/fs。w2=2*pi*f2/fs。w3=2*pi*f3/fs。 xL=sin(w1*n)+sin(w2*n)+sin(w3*n)。 w=0:2*pi/1000:2*pi。 wh=*cos(2*pi*n/(L1))。 xL_H=xL.*wh。 yL=dtft(xL_H,w)。 plot(w/2/pi*fs,abs(yL),39。k39。) axis([500,1500,0,L/3]) xlabel(39。f/Hz39。) gtext(39。L=30039。)。 仿真 結果 如下 ： 仿真結果分析： 由物理分辨率計算公式知 ： Lfscf ??? ，（ c=2）則當 L=100， fs=5kHz 時 ?f? =100 Hz， 又因為 3 個信號的譜峰間隔 f? 為 50Hz，則 ?f? f? .物理分辨率大于信號頻率間隔，故不能分辨出 3 個信號的頻譜； 當 L=300， fs=
公司管理相關推薦
鄂ICP備17016276號-1