Vortex Signal Processing Based on Spectrum Analysis

LIN Min

(Medical Instrumentation College, Shanghai Science & Technology University,
Shanghai 200093, China)

　　Key words: sampling in non?complete periods; spectrum analysis; vortex signal processing

1引言
　　渦街流量信號處理的核心問題是擴(kuò)展渦街流量計在低頻段的量程下限。而嘗試運用譜分析方法解決這一問題是目前的研究熱點之一。文獻(xiàn)［1］～［3］采用仿真方法比較了基于FFT的周期圖譜法和基于BURG的最大熵譜法，并得出結(jié)論：FFT算法所需采樣數(shù)據(jù)多，適合抑制低頻確定性噪聲；而BURG算法所需數(shù)據(jù)少，適合抑制隨機(jī)噪聲。但是他們是用信號發(fā)生器產(chǎn)生的理想波形來進(jìn)行譜分析研究的，沒有計算實際渦街信號的頻譜。另外，該頻譜分析是在對信號進(jìn)行整周期采樣的基礎(chǔ)上實現(xiàn)的，而實際系統(tǒng)由于未知信號頻率，不可能對渦街輸出信號實現(xiàn)整周期采樣，因而整周期采樣只是理想情況，在實際系統(tǒng)中幾乎不可能實現(xiàn)。
　　鑒于此，本文在將譜分析運用到實際的流量信號處理前，有必要先通過Matlab軟件來仿真分析非整周期采樣對譜分析處理結(jié)果的影響，以減少研制中的盲目性和縮短研究周期。

　　假設(shè)流量計的輸出信號為　
　　?
信號頻率fsig=20Hz。從t=0時刻開始采樣，總采樣時間保持不變。則當(dāng)φ=360°·k，k=0、1、2、…時為整周期采樣。為了比較非整周期采樣對FFT譜分析和BURG熵譜分析的影響，我們分別取N=1024點、在64個整周期的基礎(chǔ)上進(jìn)行FFT計算以及取N=128點、在16個整周期的基礎(chǔ)上進(jìn)行BURG計算，并令偏移角φ=0°，5°，10°，…，180°，即以5°步長遞增來仿真計算非整周期采樣對兩種算法的功率譜峰值P和頻率估計值f的影響。現(xiàn)將兩組計算結(jié)果分別繪成圖1和圖2，以供更直觀地分析。圖中，橫坐標(biāo)代表信號頻率值f(Hz)；縱坐標(biāo)代表放大1000倍的功率譜密度值P(W/Hz)。

　?、?FFT譜分析在偏角為0°，即整周期采樣時譜峰值P最大，對應(yīng)的頻率值f也最精確，為20Hz；隨著偏角的增大，譜峰值和頻率均減小。而BURG譜分析在偏角為0°時譜峰值P并非最大，但對應(yīng)的頻率值f卻最精確，也為20Hz；其譜峰值P與偏角之間無明顯的規(guī)律關(guān)系。
　　Ⅱ.FFT譜分析隨著偏角的增大，其各偏角所對應(yīng)的頻率值近似于線性地遞減，即偏角差值相等，所對應(yīng)的頻率差值也近似相等。而BURG譜分析不存在這種規(guī)律。
　?、?采樣周期數(shù)越多，非整周期采樣對這兩種譜分析的影響越小。例如，F(xiàn)FT譜分析在采樣64個周期時的頻率誤差是采樣32個周期時的1／2，是采樣16個周期時的1／4，是采樣8個周期時的1／8。即檢測到的頻率誤差隨采樣周期數(shù)線性減少而近似線性增大，F(xiàn)FT譜分析在64、32、16、8個整周期采樣的基礎(chǔ)上偏180°時得到的誤差分別為0.77％、1.54％、3.03％、5.88％。
　?、?FFT譜分析得到的功率譜峰值P準(zhǔn)確；譜峰值對應(yīng)的頻率值準(zhǔn)確。而BURG譜分析得到的功率譜峰值P不準(zhǔn)確，前后無規(guī)律性；但譜峰值所對應(yīng)的頻率值也準(zhǔn)確。
　?、?信號頻率值f與采樣周期數(shù)有關(guān)。只要周期數(shù)取得足夠大時，非整周期采樣帶來的測量誤差就變得很小，可以忽略不計。
　　因此，綜合考慮上面幾點，再參考文獻(xiàn)［1］～［3］的研究結(jié)論，本文選擇基于FFT譜分析來進(jìn)行實際渦街流量信號的處理。?

3.1渦街流量信號采集系統(tǒng)
　　采集系統(tǒng)由兩部分組成：渦街流量信號采集裝置和可安裝于PC機(jī)中的PC-6333多功能模入模出接口卡。采集裝置由水管、天津儀表廠的LUGB型渦街流量計、ADMACSE電磁流量計、滬東電機(jī)公司的Y90S-2三相異步電動機(jī)及水泵、西門子公司的MICROMASTER420變頻器、水箱等部分組成(圖3)。圖中，箭頭的方向表示水流的流向。
　　采集裝置的工作流程為：啟動變頻器，設(shè)置變頻器的顯示頻率，電動機(jī)在變頻器的控制下按一定的速度旋轉(zhuǎn)，帶動水泵工作，將水從水箱抽上來。水流經(jīng)過電磁流量計，顯示瞬時的流量百分比，再經(jīng)渦街流量計檢測后返回到水箱。渦街流量計的輸出模擬信號由PC-6333模入模出接口卡轉(zhuǎn)換成數(shù)字信號后送入PC機(jī)中再作進(jìn)一步處理。
3.2基于譜分析的渦街流量信號處理與分析
　　通過變頻器設(shè)置水泵工作頻率為50Hz、47.5Hz、45Hz、……，以2.5Hz的步長逐漸降低，當(dāng)?shù)陀?5Hz時，水泵電動機(jī)就自動逐漸地停止旋轉(zhuǎn)。對應(yīng)每一個頻率點，分別采樣2048點數(shù)據(jù)。將2048點采樣數(shù)據(jù)載入Matlab程序中，通過sptool工具箱對其進(jìn)行近似整周期FFT計算。為了便于比較和提高計算精度，在每個泵頻率點上，我們?nèi)×藘啥尾煌秶牟蓸訑?shù)據(jù)，每一段都近似整周期采樣，再求取兩段的FFT平均值，如圖4所示。橫坐標(biāo)表示譜分析得到的頻率平均值f(Hz)；縱坐標(biāo)表示各泵頻率點上電磁流量計顯示的流量平均值(％)。

　　(1)各泵頻率點上計算得到的整周期FFT平均值隨著“泵頻率／流量”值的逐漸降低而減小。
　　(2)根據(jù)用戶手冊知道LUGB型渦街流量計所能檢測到的流體頻率范圍為13.191Hz～131.91Hz，這是該流量計用傳統(tǒng)的電路閾值方法處理渦街信號時的頻率檢測范圍。而用頻譜分析方法處理該流量計的渦街信號，在“泵頻率／流量”為35Hz／60.58％時得到的渦街信號頻率為11.0726Hz，比用傳統(tǒng)的電路閾值方法所能檢測到的下限頻率低2Hz左右。從這點看，頻譜分析方法處理渦街信號要比傳統(tǒng)的電路閾值方法優(yōu)越，特別是在測量低流速段時。
　　(3)本實驗選取的采樣頻率為1kHz，采樣點數(shù)為2048點，因而頻率分辨率
　　
　　則系統(tǒng)在低流速時的測量相對誤差為
　　
　　由此可見，若希望測量精度為0.45％，則在信號頻率不變的情況下，分辨率應(yīng)該在0.05Hz以下，在采樣點數(shù)不變的情況下，則要求采樣頻率降低到100Hz以下，這就不能滿足香農(nóng)采樣定理。若采樣頻率保持不變，要使分辨率在0.05Hz以下，則采樣點數(shù)需要增加到20480點以上?？墒屈c數(shù)增加則增大數(shù)據(jù)存儲量，同時增大計算量，增加計算時間，會降低系統(tǒng)的實時性。而在采樣頻率和采樣點數(shù)一定的情況下，信號頻率越低，測量誤差越大。對此，用分段設(shè)置采樣頻率的辦法以達(dá)到同時滿足計算精度和系統(tǒng)實時性的要求。進(jìn)行頻率分段，要進(jìn)行采樣頻率的頻繁切換，這就很難實現(xiàn)在線的信號采集，也很難真正滿足系統(tǒng)的實時性，因此該方法不能從根本上解決問題。

　　在低流速時，特別是當(dāng)信號頻率在10Hz以下時，渦街信號和噪聲信號幾乎重疊在一起，甚至噪聲的幅值還略大于渦街信號幅值。此時，用功率譜分析方法來處理渦街流量信號，很可能得到的噪聲頻譜峰值要高于信號頻譜幅值，這樣就會將噪聲頻率錯認(rèn)為是渦街信號頻率。由此看來，單純地用頻譜分析方法要達(dá)到擴(kuò)展渦街流量計在低流速時的量程下限非常困難。
　　但是，利用頻譜分析能很好地展現(xiàn)信號的頻率分布特征，能初步提供渦街信號的頻率，為進(jìn)一步進(jìn)行渦街信號頻率的準(zhǔn)確檢測做好了準(zhǔn)備，提供了研究基礎(chǔ)和方向。在頻譜分析的基礎(chǔ)上，我們設(shè)計了渦街信號閾值處理專家系統(tǒng)，并已取得了一定的進(jìn)展。下一步，我們將繼續(xù)結(jié)合譜分析，運用人工智能和專家系統(tǒng)的方法來處理渦街流量信號，以期能對擴(kuò)展渦街流量計量程下限的問題有所突破。