電子信息工程-數(shù)字譜分析實(shí)踐及誤差討論-文庫(kù)吧
2024-12-24 10:07 本頁(yè)面
【正文】 (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。)。 仿真結(jié)果如下 ： ② 當(dāng)時(shí)窗寬度 L=200時(shí)， 設(shè)定 f1=1kHz， f2= kHz， f3= kHz； fs=5 kHz。 代碼 ：（略） 仿真結(jié)果如下： 仿真結(jié)果分析： 由 物理分辨率計(jì)算公式知： Lfscf ??? （ c=1） ，則當(dāng) L=64， fs=5kHz 時(shí) ?f? = Hz， 又因?yàn)?3 個(gè)信號(hào)的譜峰間隔 f? 為 50Hz，則 ?f? f? .物理分辨率大于 信號(hào)頻率間隔 ，故不能分辨出 3 個(gè)信號(hào) 的頻譜； 當(dāng) L=200， fs=5kHz 時(shí)， ?f? =25Hz，則 ?f? f? 物理分辨率小于信號(hào)頻率間隔，故能分辨出 3 個(gè)信號(hào)的頻譜 Hamming window ① 當(dāng)時(shí)窗寬度 L=100時(shí)，設(shè)定 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。)。 仿真 結(jié)果 如下 ： ② 當(dāng)時(shí)窗寬度 L=300 時(shí)，設(shè)定 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。)。 仿真 結(jié)果 如下 ： 仿真結(jié)果分析： 由物理分辨率計(jì)算公式知 ： Lfscf ??? ，（ c=2）則當(dāng) L=100， fs=5kHz 時(shí) ?f? =100 Hz， 又因?yàn)?3 個(gè)信號(hào)的譜峰間隔 f? 為 50Hz，則 ?f? f? .物理分辨率大于信號(hào)頻率間隔，故不能分辨出 3 個(gè)信號(hào)的頻譜； 當(dāng) L=300， fs=
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1