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