摘要:對(dui)基于譜分析(xī)的數字式渦(wo)街流量計 在(zài)應用中響應(ying)速度比較慢(màn)及其在小流(liu)量或小口徑(jìng)時流量計輸(shū)出波動兩個(ge)問題進行了(le)深入分析。采(cai)用設置門限(xiàn)值的方法提(tí)高了系統響(xiang)應速度。對頻(pin)譜分析結果(guǒ)進行多次爲(wèi)0判斷,解決了(le)系統輸出斷(duan)續問題。通過(guo)對頻譜分析(xī)結果的限幅(fu),在平均時去(qu)掉超過标準(zhun)差的頻率值(zhi),并去掉最大(dà)值和最小值(zhí),對平均結果(guǒ)進行第二次(cì)平均的方法(fa)減小了系統(tǒng)輸出的波動(dong),提高了流量(liàng)計輸出穩定(dìng)性。
0引言
　　渦街(jiē)流量計具有(yǒu)無可動部件(jiàn)維護方便、壓(ya)力損失小和(hé)輸出頻率量(liàng)等優點,可以(yǐ)測量液體、氣(qi)體和蒸汽流(liú)量,應用相當(dāng)廣泛。目前，渦(wō)街流量計主(zhu)要有兩種處(chu)理方式,即脈(mo)沖整形計數(shù)方式和數字(zì)信号處理方(fang)式。其中,數字(zì)信号處理方(fāng)式具有抗幹(gan)擾能力強、量(liang)程比較寬和(hé)精度高的優(yōu)點。基于譜分(fèn)析的渦街流(liú)量計數字信(xìn)号處理方法(fa),采用DSP(數字信(xin)号處理器)芯(xin)片研制了實(shí)時處理系統(tong),通過标定并(bing)應用于工業(ye)現場,取得了(le)很好的效果(guǒ)。在應用中，針(zhēn)對數字式渦(wo)街流量計響(xiang)應速度慢以(yǐ)及在小流量(liang)或者小口徑(jing)時輸出信号(hào)波動問題，進(jìn)行了深人的(de)分析,并提出(chu)了有效的解(jie)決方案。
1提高(gao)信号響應的(de)速度
　　采用頻(pin)譜估計方法(fǎ)分析渦街信(xin)号的頻率是(shì)基于渦街流(liu)量計信号是(shì)一個平穩随(sui)機過程信号(hao),并且流量在(zai)短時間内是(shi)穩定的或變(bian)化很小。由于(yú)渦街流量計(jì)輸出信号實(shi)時性的要求(qiú),隻能對有限(xiàn)長時間段内(nei)的信号進.行(hang)處理,即隻能(neng)截取一段數(shù)據進行頻譜(pǔ)分析。爲了防(fang)止數據統計(jì)變異性對計(jì)算渦街頻率(lü)的影響,必須(xū)對若幹次頻(pín)譜分析的結(jié)果進行平均(jun),這就需要花(hua)費時間。當平(ping)均次數多時(shi)，數字渦街流(liú)量計的響應(ying)速度就會顯(xian)得比較慢。例(li)如,我們的系(xì)統采用時鍾(zhong)頻率爲40MHz的TMS320LF2407ADSP芯(xin)片，使用彙編(biān)語言編制FFT程(chéng)序計算渦街(jiē)信号的頻譜(pu),對50次計算結(jié)果進行平均(jun1)。每次FFT計算時(shí)間大概是300ms(包(bao)括信号處理(li),人機接口等(deng)對應的其他(tā)軟件執行時(shí)間),總用時大(dà)概爲15s,這個時(shí)間比較長。若(ruo)在這段時間(jiān)裏,流量發生(sheng)突然變化，系(xi)統就需要計(ji)算50次FFT,才能反(fǎn)映流量的變(biàn)化,這将造成(chéng)測量誤差,如(ru)圖1所示。
 
　　在實(shí)際處理中,将(jiāng)每次頻譜分(fen)析後得到的(de)渦街頻率放(fàng)人一個循環(huán)隊列中,該隊(duì)列稱爲隊列(liè)Q1,如圖2所示(數(shu)字表示隊列(lie)成員序号)。其(qí)中,最右邊的(de)爲隊首Q1,最左(zuǒ)邊的爲隊尾(wei)Q1。再按先進先(xiān)出的原則對(dui)隊列中最新(xin)的若幹個值(zhí)進行平均,得(de)到渦街輸出(chu)頻率。由于對(dui)隊列中的值(zhí)不加判斷而(ér)直接參與平(ping).均,導緻系統(tǒng)響應速度變(bian)慢。
 
　　爲此,通過(guò)設置門限值(zhí)的方法來解(jie)決這個問題(ti)。根據隊列Q1中(zhōng)保存的過去(qu)的流量信息(xi)來設置一個(ge)頻率突變門(men)限frelimit,如果連續(xu)計算出N個頻(pin)率值Q1[N]~Q1[1]都超過(guò)了fre._limit,則對這N個(gè)頻率值和後(hòu)面進人隊列(lie)的值進行平(ping)均,而前面的(de)值Q1[N+1]~Q1[50]不參加平(píng)均,N一般爲3~5。理(li)論上，門限fre_limit的(de)選擇越小對(dui)微小突變就(jiù)越敏感。但是(shì),由于FFT計算出(chu)來的頻率會(hui)在--定範圍内(nei)波動,所以,門(mén)限值不能選(xuǎn)擇太小。另外(wai)，由于小的突(tu)變造成的誤(wù)差不會太大(dà),所以也沒必(bì)要将門限設(shè)置得很小。
　　一(yi)般來說，渦街(jiē)流量計在小(xiao)流量時(低頻(pín)處)頻率計算(suan)結果波動比(bǐ)較大,而大流(liu)量時(高頻處(chù))比較穩定。所(suǒ)以,門限值可(ke)以根據渦街(jie)流量計所能(neng)測的最低頻(pin)率來設置。由(you)于流量計的(de)最高可測頻(pin)率是已知的(de),而采樣頻率(lü)一般設爲略(luè)大于最高頻(pin)率的2倍,所以(yi),門限值可以(yǐ)根據式(1)來計(jì)算:
 
　　式中:ƒs爲采(cǎi)樣頻率;Ɩ爲儀(yi)表量程比;α(0<lαl<1)爲(wei)系數。調整α可(ke)以調整門限(xiàn)值大小,α<0時爲(wèi)下限;α>0時爲上(shàng)限。一般可以(yǐ)取±4%~±13%。
　　當流量突(tu)變不超過門(mén)限時,在流量(liàng)計調節時間(jian)内累積流量(liang)的最大相對(dui)誤差爲:
 
　　式中(zhong):n(n》1)爲平均次數(shù);ƒp爲突變以前(qián)的頻率;ƒr爲突(tu)變以後的頻(pin)率。而當流量(liang)突變超過fre_limit時(shi),誤差會非常(cháng)小。
　　實際流量(liang)變化是一個(gè)連續變化的(de)過程。但是,隻(zhī)要流量對應(ying)的信号頻率(lü)的變化率超(chao)過fre_limit/tc,其中tc爲FFT計(ji)算的間隔時(shí)間，則系統就(jiù)能快速響應(ying);而流量對應(yīng)的信号頻率(lü)的變化率小(xiǎo)于fre_limit/tc.時,其最大(da)相對誤差小(xiǎo)于式(2)的計算(suàn)值。
　　圖3爲經過(guo)加速處理後(hòu)40口徑渦街流(liu)量計在流量(liàng)突變時的響(xiang)應曲線。可見(jiàn),在流量突變(biàn)後,系統經過(guo)3次FFT計算就能(néng)跟上流量的(de)變化,有效地(dì)提高了流量(liang)計的響應速(su)度。
 
2.1消除輸出(chū)信号的斷續(xu)現象
　　頻譜分(fèn)析法是按能(néng)量最大原則(zé),即把FFT計算出(chū)來的功率譜(pu)中能量最大(da)的頻率值作(zuo)爲渦街頻率(lǜ)。但是,由于幹(gàn)擾的存在,不(bú)得不設置一(yī)個最小截止(zhǐ)幅值(min_amp)來限制(zhi)小幅值幹擾(rǎo)信号,即當FFT計(ji)算幅值低于(yu)最小截止幅(fú)值時,流量計(ji)輸出爲0,從而(er)提高抗幹擾(rǎo)能力;而由于(yú)信号質量、對(dui)采樣信号加(jia)窗FFT非整周期(qi)采樣等原因(yin),導緻FFT計算幅(fú)值會有波動(dong)。當流量在流(liu)量計下限附(fu)近時,會出現(xiàn)FFT計算幅值低(dī)于最小截止(zhǐ)幅值的情況(kuàng),這将導緻流(liu)量計輸出出(chu)現斷續現象(xiàng)。
　　流量計斷續(xù)問題主要是(shi)由于偶爾幾(jǐ)次的FFT計算幅(fú)值低于最小(xiao)截止幅值所(suǒ)緻。通過實驗(yan)發現，在流量(liang)計可測流量(liàng)範圍内-般不(bu)會出現連續(xù)幾次FFT計算幅(fú)值均小于最(zuì)小截止幅值(zhi)的情況，所以(yǐ)在碰到FFT計算(suan)幅值低于最(zuì)小截止幅值(zhí)時不是直接(jiē)輸出0,而是将(jiāng)0放入隊列Q1。在(zài)計算平均頻(pín)率以前先對(duì)隊列Q1中最近(jìn)的N(N可取3~4)個值(zhí)進行判斷,如(rú)果N個值，全部(bù)爲0,則流量計(ji)輸出爲0,否則(zé)繼續做平均(jun)計算,将平均(jun)值作爲流量(liang)計的輸出頻(pin)率,這樣就可(kě)以解決輸出(chu)斷續的問題(ti)。
2.2解決輸出信(xìn)号的波動問(wèn)題
　　在小流量(liang)的時候渦街(jiē)流量計信号(hao)微弱,信噪比(bi)降低,信号質(zhì)量比較差,頻(pin)譜分析得到(dao)的渦街頻率(lǜ)波動比較大(dà)。特别是小口(kǒu)徑渦街流量(liang)計，由于傳感(gan)器安裝等原(yuán)因,這種現象(xiàng)更爲明顯。圖(tú)4中,曲線1爲15口(kou)徑渦街流量(liàng)計在下限頻(pin)率附近,每次(cì)FFT計算得到的(de)渦街頻率;曲(qǔ)線2爲40口徑離(lí)下限頻率比(bi)較遠時每次(ci)FFT計算後得到(dao)的渦街頻率(lǜ)。可見，曲線1波(bo)動較大，數值(zhí)約爲30Hz。此時,即(jí)便采取平均(jun)措施,也無法(fa)使輸出頻率(lǜ)值穩定。
 
　　波動(dòng)問題主要是(shi)由于信号較(jiao)差,每次FFT計算(suan)出的頻率值(zhi)偏差太大所(suǒ)造成的。所以(yi),首先必須保(bǎo)證隊列Q1中的(de)所有頻率值(zhí)在可信範圍(wéi)内,即在流量(liang)沒有突變的(de)情況下,隊列(lie)Q1中的頻率值(zhí)不能與前次(ci)輸出頻率差(chà)得太多。爲此(ci)，可以設置一(yi)個頻率允許(xǔ)跳動範圍的(de)限幅值。當計(jì)算的頻率值(zhí)超出範圍時(shi),不是簡單地(dì)将這個值去(qù)掉,而是用最(zuì)近一次的輸(shū)出頻率加上(shang)限幅值來代(dai)替這個值。這(zhè)在某種程度(du)上能保留反(fan)映流量變化(huà)趨勢的部分(fèn)信息。通過實(shi)驗發現,可以(yi)根據fre__limit來設置(zhi)這個限幅值(zhi)。
處理波動問(wèn)題的具體步(bu)驟爲:
①對隊列(lie)Q,的第N個元素(sù)Q1[N]用式(3)處理,得(de)到數值ƒin爲:
 
　　式(shì)中:ƒin爲處理後(hou)隊列Q1[N]中的頻(pín)率值;f。爲處理(li)前隊列Q1N]中的(de)頻率值,即原(yuán)先FFT計算出來(lái)的頻率值;f,爲(wèi)流量計最近(jìn)--次輸出的頻(pín)率值;c爲常數(shu),可以根據實(shi)際情況選取(qu),實驗中取1。需(xū)要注意的是(shi),這裏是對Q1[N]處(chù)理,而不是對(duì)隊首元素Q1[1]處(chu)理。原因是,Q1[1]~Q1[N-1]中(zhōng)的信息會反(fan)映流量的突(tu)變,如果對其(qí)進行限幅處(chù)理,會影響系(xi)統的響應速(sù)度。
②對隊列Q1中(zhōng)的值進行去(qu)N-1次最大、最小(xiǎo)值處理,算出(chu)均值μ,根據均(jun1)值算出标準(zhun)差σ,然後,再利(lì)用标準差對(dui)隊列Q1中的值(zhi)進行平均。而(ér)超出标準差(chà)的數據不參(cān)加平均。
③經過(guò)以上處理後(hou)得到的頻率(lǜ)已經比較穩(wen)定,但仍有一(yī)定範圍的波(bo)動。如果流量(liang)是穩定的話(hua),流量計輸出(chu)應該接近穩(wen)定,波動比較(jiào)小。因此,可以(yǐ)考慮對計算(suàn)頻率值進行(hang)第二次平均(jun),相當于進行(hang)二次濾波。爲(wei)此,建立另一(yī)個隊列Q2。隊列(liè)Q2中的元素是(shì)隊列Q;經過上(shang)述步驟處理(li)後的結果;再(zài)對隊列Q2進行(hang)平均計算,得(dé)到的均值作(zuò)爲流量計的(de)輸出。注意,這(zhè)裏對隊列Q2也(yě)要進行加速(su)處理,方法與(yǔ)處理Q1的類似(si)。圖5所示爲經(jing)過處理以後(hòu)渦街流量計(ji)的輸出，從圖(tú)中可以看出(chū)流量計輸出(chu)已經比較平(ping)穩。
 
3結束語
　　根(gēn)據頻譜分析(xi)後的結果,采(cai)用設置門限(xiàn)值的方法對(dui)流量計輸出(chu)進行加速處(chù)理,并分析了(le)加速處理方(fang)法的誤差,有(yǒu)效地提高了(le)流量計的響(xiang)應速度降低(dī)了調節時間(jiān)内由于響應(ying)慢造成的誤(wu)差。
　　對頻譜分(fen)析後的結果(guǒ)進行0值判斷(duan),連續幾次爲(wei)0,流量計輸出(chū)才爲0,消除了(le)流量計在可(ke)測流量下限(xian).時輸出斷續(xù)現象。
　　在對頻(pin)譜分析後的(de)結果平均前(qián)采取了限幅(fu)、去掉超過标(biāo)準差以外的(de)值和去掉最(zui)大值最小值(zhi)的操作,并對(dui)平均後的結(jie)果進行第二(èr)次平均,減小(xiao)了流量計輸(shu)出信号的波(bo)動,提高了流(liú)量計輸出的(de)穩定性。

文章(zhāng)來源于網絡(luò),如有侵權聯(lián)系即删除！