基于LabVIEW的數字變頻FFT設計

發(fā)布時(shí)間:2010-8-17 13:24    發(fā)布者:lavida
關(guān)鍵詞: fft , LabVIEW , 數字變頻
在運動(dòng)目標速度測量中,常利用頻譜分析的方法獲取目標的多普勒頻率,并依據多普勒測速原理來(lái)完成動(dòng)目標速度測量。為達到高精度測速的要求,需進(jìn)一步提高頻率分辨率,在實(shí)際頻譜分析中,要對獲取的試驗數據先進(jìn)行分段處理,在此基礎上再進(jìn)行細化操作,這樣可獲得比常規FFT分析更高的頻率分辨率。近年來(lái),頻譜細化技術(shù)發(fā)展迅速,常見(jiàn)的方法有:HR-FA法,基于多相濾波器的ZFFT法,基于復調制的Zoom-FFT法,自適應Zoom-FFT法,Chirp-Z變換和小波基法等。然而,這些頻譜細化技術(shù)普遍存在運算量大,不易實(shí)現編程的缺點(diǎn)。為此,提出了一種數字變頻FFT的頻譜細化算法,并利用圖形化編程語(yǔ)言LabVIEW進(jìn)行了編程設計。  

1 數字變頻FFT的數學(xué)原理  

在頻譜分析中,頻率分辨率表示頻譜中能夠分辨的兩個(gè)頻率分量的最小間隔,用頻率間隔△f表示為:  


  
要提高FFT的頻率分辨率,可通過(guò)以下兩種途徑實(shí)現:  

(1)降低采樣頻率fs。這會(huì )使頻率分析范圍縮小,其降低的幅度受到采樣定律的限制。  

(2)需要增加分析的采樣點(diǎn)數N。這意味著(zhù)計算機的存儲量和計算量大大增加,由于實(shí)際系統軟、硬件方面的限制,這樣做并不總是可能的。  

可以看出以上兩種方法提高頻率分辨率的能力有限且靈活性差。所用的數字變頻FFT主要指采用移頻特性進(jìn)行頻譜細化的技術(shù),其原理框圖如圖1所示。  


  
設模擬信號為x(t),經(jīng)A/D轉換后得到采樣時(shí)間序列x(n)(n=0,1,…,N-1),設fs為采樣頻率;f1~f3為細化分析頻帶;f0為需要細化的頻帶中心頻率;D為細化倍數;N為FFT分析的點(diǎn)數,算法如下:  

(1)復調制移頻  

所謂復調制移頻就是將頻域坐標向左移或向右移,使得被觀(guān)察的起點(diǎn)為頻域坐標的零頻位置。這里對離散信號x(n)用exp[-2πf0/fs]進(jìn)行復調制,把需要細化的頻帶起點(diǎn)移至頻率軸原點(diǎn),得到:  


  
(2)數字低通濾波  

為保證重新采樣后不發(fā)生頻譜混疊,必須進(jìn)行抗混疊濾波,以濾出所需分析頻段信號。設頻率細化倍數為D,則低通濾波器的截止頻率fc=fs/2D。  

(3)重新采樣  

信號被移頻和低通濾波后,分析信號頻帶變窄,因而可以以較低的采樣頻率fs=fs/D進(jìn)行重采樣,fs比原采樣頻率降低了D倍,即對原采樣點(diǎn)每隔N點(diǎn)再抽樣一次。  

(4)反移頻處理  

實(shí)行反移頻操作就是將頻率中心重新移到需要細化的頻帶起始頻率,使得移頻前后的頻率保持一致。  

(5)FFT運算  

對反移頻后的信號進(jìn)行FFT處理,得到細化后的頻譜,其頻率分辨率提高了D倍。  


  
2 數字變頻FFT的LabVIEW實(shí)現  

2.1 虛擬儀器LabVEW  

LabVIEW是美國NI公司推出的虛擬儀器開(kāi)發(fā)平臺軟件,是一種非常優(yōu)秀的面向對象的圖形化編程語(yǔ)言,用于快速創(chuàng )建測試、測量和控制應用程序。它的優(yōu)點(diǎn)是數據處理速度快,硬件支持等方面功能強大;缺點(diǎn)是在數值處理、分析和算法工具等方面的效率不高。而MathWorks公司開(kāi)發(fā)的Matlab提供了強大的矩陣運算和圖形處理功能,編程效率高,特別擅長(cháng)數值分析和處理,但其界面開(kāi)發(fā)能力較差,不能進(jìn)行實(shí)時(shí)操作和控制。  

鑒于LabVIEW和Matlab兩種語(yǔ)言的優(yōu)點(diǎn),在虛擬儀器開(kāi)發(fā)過(guò)程中,除利用LabVIEW直接編程外,還可結合Matlab進(jìn)行混合編程。通常用LabVIEW設計用戶(hù)圖形界面,負責數據采集和網(wǎng)絡(luò )通信;Matlab在后臺提供大型算法供LabVIEW調用。  

2.2 LabVIEW直接編程實(shí)現數字變頻FFT  

在LabVIEW直接編程中,通過(guò)調用Ramp Pattern.vi節點(diǎn)和Exponential.vi節點(diǎn)完成信號的復調制和反移頻操作;通過(guò)調用Decimate.vi節點(diǎn)對復調制后的信號進(jìn)行濾波,壓縮頻帶,然后進(jìn)行整數倍抽取。其中,對信號進(jìn)行數字低通濾波的濾波器為切比雪夫I型濾波器,其主要參數:階數為8階,截止頻率為0.8·(fs/2)/D;通過(guò)調用FFT.vi和Array Size.vi節點(diǎn)完成FFT運算,輔以其他相應的計算處理節點(diǎn),根據數字變頻的原理圖最終實(shí)現數字變頻FFT,其程序框圖如圖2所示。  


  
2.3 LabVIEW與Matlab混合編程實(shí)現數字變頻FFT  

LabVIEW與Matlab的混合編程,可通過(guò)調用LabVIEW中的Matlab Script節點(diǎn)實(shí)現。其中,MatlabScript節點(diǎn)本身具有多輸入、多輸出的特點(diǎn),一次處理的信息量可以很大。在Matlab中,根據數字變頻FFT的數學(xué)原理,編程實(shí)現復調制移頻、低通濾波、重采樣、反移頻和FFT操作處理,經(jīng)調試無(wú)誤后,導入到Matlab Script節點(diǎn)中;然后在LabVIEW中,通過(guò)調用SineWave.vi產(chǎn)生仿真信號,或從文件中讀取信號數據,同時(shí)添加采樣頻率,細化倍數等控制節點(diǎn);最后連接各圖標,實(shí)現數字變頻FFT,其程序框圖如圖3所示。  


  
3 仿真分析  

利用LabVIEW編程實(shí)現數字變頻FFT軟件處理平臺,調用Functions\Analyze\Signal Processing\Signal Generation子模板中的Sine Wave.vi創(chuàng )建正弦信號發(fā)生器,構造仿真信號:  

x(t)=sin(2πf1t)+2sin(2πf2t)+5sin(2πf3t)  

式中:f1=2 002 Hz;f2=2 004 Hz;f3=2 006 Hz;采樣點(diǎn)數N=5 120;采樣頻率為51 200 Hz。根據式(1) 可知,此時(shí)頻率分辨率為10 Hz,在頻域內分辨不出這3個(gè)信號。  

若要把分辨率提到1 Hz,即細化10倍,就要采樣51 200個(gè)點(diǎn),然后把分析頻帶(2 000~2 010 Hz)的起始頻率f=2 000 Hz點(diǎn)移到原點(diǎn),當細化倍數D=10時(shí),即51 200個(gè)采樣點(diǎn)每隔10個(gè)點(diǎn)進(jìn)行抽取,完成對復調制移頻、濾波后的信號重采樣。新的采樣頻率即為5 120 Hz,降低了10倍,抽取得到5 120個(gè)點(diǎn)。為了使細化后的頻率與細化前的一致,在作FFT前應該實(shí)行反移頻,這樣就可以得到分析頻帶上的細化頻譜。如果采用LabVIEW直接編程處理,其細化頻譜如圖4所示;如果采用LabVIEW和Matlab混合編程處理,其細化頻譜如圖5所示。  





  
從圖4和圖5中可以看出,基于LabVIEW的兩種編程方法都實(shí)現了頻譜細化的功能,2 002 Hz,2 004 Hz和2 006 Hz三個(gè)頻率點(diǎn)對應的幅值譜清晰可見(jiàn),且幅值相差2.5倍,分辨率為1 Hz。  

4 結語(yǔ)  

數字變頻FFT是頻譜分析中一種約束條件少,可操作性強的方法。在此借助功能強大的LabVIEW軟件編程處理方法,使得數字變頻實(shí)現簡(jiǎn)便,能夠滿(mǎn)足提高頻率分辨率的要求,并具有很高的實(shí)時(shí)性。  

顯然,在動(dòng)目標速度測量中,利用數字變頻FFT進(jìn)行頻譜細化處理,可獲得更高的頻率測量精度。根據多普勒原理,也可獲得更高的測速精度。因此,研究中基于LabVIEW的數字變頻FFT頻譜細化方法在工程實(shí)踐中有著(zhù)重要的應用價(jià)值,可廣泛應用于汽車(chē)、飛機等運動(dòng)目標的速度測量中。
本文地址:http://selenalain.com/thread-21861-1-1.html     【打印本頁(yè)】

本站部分文章為轉載或網(wǎng)友發(fā)布,目的在于傳遞和分享信息,并不代表本網(wǎng)贊同其觀(guān)點(diǎn)和對其真實(shí)性負責;文章版權歸原作者及原出處所有,如涉及作品內容、版權和其它問(wèn)題,我們將根據著(zhù)作權人的要求,第一時(shí)間更正或刪除。
您需要登錄后才可以發(fā)表評論 登錄 | 立即注冊

相關(guān)視頻

關(guān)于我們  -  服務(wù)條款  -  使用指南  -  站點(diǎn)地圖  -  友情鏈接  -  聯(lián)系我們
電子工程網(wǎng) © 版權所有   京ICP備16069177號 | 京公網(wǎng)安備11010502021702
快速回復 返回頂部 返回列表
午夜高清国产拍精品福利|亚洲色精品88色婷婷七月丁香|91久久精品无码一区|99久久国语露脸精品|动漫卡通亚洲综合专区48页