摘要:爲(wèi)了改善(shan)渦街流(liú)量計 的(de)性能,提(ti)出了一(yi)種基于(yú)卡爾曼(màn)濾波的(de)渦街信(xìn)号處理(lǐ)方法。根(gen)據渦街(jiē)信号的(de)特點,設(shè)計了一(yi)個線性(xìng)渦街信(xìn)号模型(xing)。結合模(mó)糊搜索(suo)和叠代(dài)算法,通(tōng)過分析(xī)卡爾曼(màn)濾波器(qì)算法的(de)原理和(hé)關鍵參(cān)數,改進(jìn)卡爾曼(màn)濾波器(qì)算法。通(tōng)過仿真(zhen)模拟和(he)實際流(liu)量實驗(yan)驗證了(le)所提出(chu)的方法(fǎ),并與其(qi)他方法(fǎ)進行了(le)比較。實(shi)驗結果(guo)表明,所(suǒ)提出的(de)方法具(jù)有自适(shi)應濾波(bo)、抗幹擾(rao)能力和(he)濾波速(sù)度的優(you)點。
0引言(yán)
　　渦街流(liu)量計作(zuo)爲一種(zhǒng)振動型(xíng)流量計(jì),具有應(ying)用範圍(wei)廣、測量(liàng)介質多(duo)耐高溫(wen)、耐高壓(yā)等優點(diǎn)，具有廣(guang)闊的發(fā)展前景(jǐng)。渦街流(liú)量計利(li)用流體(tǐ)振動的(de)原理來(lái)測量流(liu)量。當流(liu)體通過(guò)一個垂(chui)直放置(zhì)的非流(liu)線型旋(xuan)渦發生(shēng)體時,發(fa)生體兩(liang)側會産(chǎn)生兩排(pái)交錯排(pai)列的旋(xuan)渦,被稱(chēng)爲卡門(men)渦街,如(rú)圖1所示(shi)。發生體(tǐ)後方的(de)應力型(xing)壓電傳(chuan)感器将(jiāng)旋渦産(chan)生的壓(ya)力差轉(zhuǎn)換爲電(diàn)荷信号(hao)。電荷信(xìn)号的變(biàn)化頻率(lǜ)與旋渦(wō)的産生(shēng)頻率-緻(zhì)。
 
渦街信(xin)号具有(you)以下兩(liang)個特點(dian)。
1)壓電傳(chuan)感器輸(shu)出的電(dian)荷信号(hào)QH滿足正(zheng)弦規律(lü),如式(1)所(suo)示
QH=γρƒ2sin(2πƒt+φa)(1)
　　式中(zhong):γ爲壓電(dian)傳感器(qì)的系數(shu),ρ爲流體(ti)密度,ƒ爲(wei)渦街頻(pín)率,φa爲初(chu)始相位(wei)。
2)在渦街(jie)流量計(jì)中,電荷(he)放大器(qì)一般用(yòng)于将傳(chuán)感器輸(shu)出的電(dian)荷信号(hào)轉換成(cheng)電壓信(xìn)号。電壓(yā)信号是(shi)待處理(li)的原始(shǐ)信号,和(hé)電荷信(xin)号一樣(yàng)爲正弦(xián)波形式(shì)。當流體(ti)密度和(hé)發生體(tǐ)的寬度(du)爲固定(dìng)值時，渦(wō)街信号(hào)的幅值(zhi)與頻率(lǜ)的平方(fang)成正比(bi),如式(2)所(suo)示
α∞ƒ2(2)
式中(zhong):α爲渦街(jie)信号的(de)幅值。
　　由(you)于渦街(jiē)流量計(ji)的傳感(gan)器屬于(yú)振動型(xing)傳感器(qì),很容易(yi)受到振(zhèn)動幹擾(rao)，這對渦(wo)街流量(liang)計的測(ce)量精度(dù)和測量(liàng)範圍有(you)很大的(de)影響。由(you)于渦街(jie)信号的(de)幅值與(yu)頻率的(de)平方成(chéng)正比，在(zai)高流速(su)下信号(hao)的信噪(zao)比大,信(xìn)号受噪(zao)聲影響(xiǎng)小;在低(dī)流速下(xià)信号的(de)信噪比(bǐ)小,信号(hào)受噪聲(sheng)影響大(dà),增加了(le)信号檢(jiǎn)測的難(nan)度。國内(nei)外衆多(duo)學者對(dui)這一問(wèn)題進行(háng)了研究(jiu)(2-3]。近年來(lái)，卡爾曼(man)濾波方(fāng)法也開(kai)始被應(yīng)用在渦(wo)街流量(liang)計的信(xin)号處理(lǐ)中。
　　卡爾(er)曼濾波(bo)是科學(xué)家R.E.Kalman等在(zai)1960年提出(chū)的一種(zhǒng)适用于(yú)離散随(sui)機非平(ping)穩系統(tǒng)的最優(yōu)估計算(suàn)法。它基(jī)于線性(xing)離散系(xi)統,将最(zuì)優濾波(bō)理論與(yu)狀态空(kong)間思想(xiǎng)相結合(hé)。宋開臣(chen)等[4]針對(duì)壓電式(shì)渦街流(liú)量計抗(kang)幹擾性(xìng)差的缺(que)點,提出(chū)了基于(yu)多傳感(gan)器融合(he)的渦街(jie)信号檢(jian)測方法(fǎ)。該方法(fǎ)通過無(wú)迹卡爾(er)曼濾波(bō)算法将(jiang)壓差傳(chuan)感器測(ce)量的鈍(dun)體前後(hòu)壓差和(he)其下遊(you)的渦街(jie)信号頻(pín)率進行(háng)融合,增(zēng)強了壓(ya)電式渦(wo)街流量(liang)計的抗(kàng)振能力(li),有效提(tí)高了數(shu)字帶通(tōng)濾波器(qi)的測量(liang)精度。Shao等(deng)[5]針對渦(wō)街流量(liang)計提出(chū)了一種(zhǒng)基于分(fen)段卡爾(er)曼濾波(bo)的數字(zì)信号處(chu)理方法(fa)。該方法(fa)首先分(fèn)析瞬态(tai)沖擊的(de)特性,建(jian)立數學(xue)模型,然(ran)後在流(liú)量信号(hào)數據中(zhōng)找到含(han)有強瞬(shùn)态沖擊(ji)的數據(jù)段,并對(dui)數據段(duàn)進行卡(kǎ)爾曼濾(lǜ)波,以降(jiang)低瞬态(tai)沖擊的(de)功率。
　　本(běn)文提出(chu)了一種(zhong)基于渦(wō)街信号(hào)模型的(de)卡爾曼(man)濾波的(de)信号處(chu)理方法(fa)(以下簡(jiǎn)稱“本方(fang)法")。首先(xiān),根據微(wēi)分原理(li)和線性(xìng)矩陣對(dui)非線性(xing)的渦街(jie)信号進(jin)行線性(xing)化處理(li)并建立(li)模型;其(qi)次,根據(jù)渦街信(xin)号的幅(fú)頻關系(xì),初始化(huà)渦街信(xin)号模型(xing)頻率;再(zài)次,将濾(lǜ)波後的(de)輸出頻(pín)率作爲(wei)下一個(ge)循環的(de)初始渦(wō)街模型(xíng)頻率進(jin)行叠代(dai)和模糊(hú)搜索,直(zhi)到輸出(chū)頻率與(yǔ)模型頻(pin)率的誤(wu)差在渦(wō)街流量(liang)計允許(xǔ)的誤差(cha)範圍内(nei)爲止;最(zuì)後,通過(guo)仿真實(shí)驗和實(shi)流實驗(yan)對該方(fang)法進行(hang)了驗證(zhèng)。
1卡爾曼(man)濾波原(yuán)理和渦(wō)街系統(tong)模型
1.1卡(kǎ)爾曼濾(lü)波原理(lǐ)
　　卡爾曼(man)濾波是(shì)一種利(li)用線性(xìng)系統狀(zhuang)态方程(cheng),通過系(xì)統輸人(rén)的觀測(ce)數據對(duì)系統狀(zhuang)态進行(hang)最優估(gū)計的算(suan)法。其基(jī)本原理(li)如下:假(jiǎ)設有一(yī)個離散(san)的線性(xing)系統Xk,通(tong)過k-1時刻(ke)的最優(you)估計Xk-1得(de)到k時刻(kè)的預測(ce)值Xk|k-1，并用(yòng)k時刻的(de)觀測值(zhí)Zk修正預(yù)測值，從(cong)而得到(dào)h時刻的(de)最優估(gū)計Xk。圖2顯(xian)示了卡(kǎ)爾曼濾(lǜ)波原理(lǐ)。
 
　　對于沒(mei)有控制(zhì)輸人的(de)系統,卡(kǎ)爾曼濾(lǜ)波算法(fǎ)的狀态(tài)方程和(he)觀測方(fang)程可用(yòng)式(3).式(4)表(biǎo)示
Xk+1=AXk+Bwk(3)
Yk+1=HXk+1+Dʋk+1(4)
　　式中(zhōng):Xk爲n維的(de)狀态變(bian)量在k時(shi)刻的值(zhí),wk爲p維的(de)過程噪(zao)聲,Yk+1爲m維(wéi)的觀測(ce)變量,ʋk+1爲(wèi)m維的觀(guān)測噪聲(shēng),A爲變量(liang)Xk的狀态(tài)轉移矩(ju)陣,H爲系(xì)統參數(shù)矩陣,B爲(wèi)過程噪(zào)聲的系(xì)數矩陣(zhen),D爲觀測(cè)噪聲的(de)系數.矩(ju)陣。wk和ʋk+1是(shi)均值爲(wèi)0且互不(bu)相關的(de)高斯白(bai)噪聲。不(bú)難看出(chu),由式(3)和(he)式(4)構建(jiàn)的系統(tong)模型不(bu)包含渦(wō)街信号(hao)的特征(zhēng)。
1.2渦街系(xi)統模型(xing)
　　由于渦(wo)街信号(hao)是非線(xiàn)性正弦(xián)波信号(hao),不滿足(zú)卡爾曼(màn)濾波系(xì)統模型(xíng)的線性(xing)要求,不(bu)能成爲(wèi)卡爾曼(man)濾波算(suan)法的系(xi)統模型(xíng),需要利(li)用微分(fen)原理和(hé)線性矩(jǔ)陣對渦(wo)旋信号(hao)進行線(xian)性化處(chu)理,具體(ti)推導過(guo)程如下(xia)。
　　假設,渦(wo)街信号(hao)的數學(xué)模型如(rú)式(5)所示(shi):
s(t)=asin(2πƒt)(5)
那麽，其(qí)二階導(dǎo)數可以(yǐ)用式(6)表(biǎo)示:
s"(t)=-4π2aƒ2sin(2πƒt)(6)
将式(shi)(5)代人式(shì)(6),得到式(shì)(7):
s"(1)=-4π2ƒ2s(t)(7)
　　根據導(dǎo)數的定(dìng)義,當t>△t且(qiě)△t→0時,可得(dé)式(8):
 
2算法(fǎ)實現
2.1算(suan)法設計(ji)
　　卡爾曼(màn)濾波是(shì)用觀測(cè)量(實際(jì)信号)對(duì)預測變(biàn)量(模型(xing)信号)進(jìn)行修正(zheng)，濾波結(jié)果介于(yu)實際信(xin)号和模(mo)型信号(hào)之間。同(tóng)樣,濾波(bō)後的信(xìn)号頻率(lü)也介于(yú)實際信(xin)号頻率(lü)和模型(xing)信号頻(pín)率之間(jian)。爲此,設(shè)計了一(yī)種基于(yu)渦街模(mo)型的卡(kǎ)爾曼濾(lü)波算法(fǎ),以叠代(dài)的方法(fǎ)搜索渦(wo)街信号(hào)的頻率(lǜ)。
　　首先,根(gēn)據渦街(jiē)信号的(de)幅頻關(guān)系設置(zhi)初始系(xi)統模型(xing)頻率。按(an)照式(2)對(dui)液體介(jie)質管道(dào)上采集(ji)到的渦(wō)流信号(hào)的幅值(zhi)和頻率(lǜ)進行二(er)次多項(xiàng)式拟合(hé),得到在(zai)液體介(jie)質中50mm口(kǒu)徑渦街(jiē)流量計(ji)信号的(de)幅值和(he)頻率的(de)關系,如(ru)式(15)所示(shì)。
α=1.789x10-5ƒ2(15)
　　同理,對(duì)采集到(dao)的氣體(ti)數據進(jìn)行拟合(he),可以得(dé)到氣體(tǐ)信号的(de)幅值與(yǔ)頻率的(de)關系，如(ru)式(16)所示(shì)
α=2.622x10-8ƒ2(16)
　　用ƒm表示(shì)狀态模(mó)型頻率(lü),用ƒmax,表示(shi)渦街信(xin)号的最(zui)大頻率(lǜ),令ƒm=ƒmax。這樣(yàng)設置的(de)目的是(shì)減少叠(die)代次數(shu)和計算(suan)量。于是(shì),式(3)中的(de)系數矩(jǔ)陣A可用(yòng)式(17)表示(shì)。
 
　　當流速(sù)低時,渦(wō)街信号(hào)能量弱(ruò),噪聲較(jiao)大,因而(er)噪聲系(xi)數D較大(da);反之,噪(zao)聲系數(shu)D較小。由(yóu)此可見(jiàn),噪聲系(xì)數D與渦(wo)流頻率(lǜ)ƒ成反比(bi)。多次實(shi)驗數據(jù)分析表(biǎo)明,當D爲(wèi)觀測信(xìn)号αmax與渦(wo)街信号(hao)模型幅(fú)值α之比(bǐ)時,得到(dao)了理想(xiang)的良好(hao)濾波效(xiao)果，如式(shì)(18)所示。
 
　　式(shi)中:ϒ爲不(bú)同介質(zhi)中幅頻(pín)關系的(de)系數。
　　以(yǐ)50mm.口徑管(guǎn)道的液(ye)體介質(zhi)爲例,對(duì)不同流(liu)量點采(cǎi)集的實(shí)驗數據(jù)進行噪(zao)聲系數(shu)D和信号(hao)頻率ƒ的(de)曲線拟(nǐ)合，拟合(he)得到的(de)關系式(shi)如式(20)所(suǒ)示。
 
 
　　最後(hou),對最優(yōu)估計Xk的(de)周期進(jìn)行統計(jì)分析,去(qù)除組内(nei)雜散數(shu)據後,取(qu)平均值(zhí)的倒數(shù)作爲渦(wo)街信号(hao)的頻率(lǜ),以得到(dào)的頻率(lǜ)爲渦街(jiē)信号的(de)新系統(tong)模型頻(pin)率,對原(yuan)始信号(hào)進行卡(ka)爾曼濾(lǜ)波。由于(yú)原始信(xìn)号中渦(wō)街信号(hào)的頻率(lǜ)保持不(bu)變,濾波(bō)器輸出(chū)頻率介(jie)于渦街(jiē)信号頻(pin)率和模(mó)型頻率(lǜ)之間，濾(lü)波器輸(shu)出頻率(lǜ)和模型(xing)頻率在(zai)叠代中(zhong)逐漸收(shou)斂到渦(wo)旋信号(hao)頻率。當(dang)輸出頻(pin)率與模(mo)型頻率(lü)的相對(duì)誤差在(zài)預設值(zhí)以内時(shi),停止叠(die)代，最終(zhong)輸出渦(wo)街頻率(lü)。
2.2算法流(liú)程
具體(ti)的算法(fa)步驟整(zheng)理如下(xià)。
步驟一(yī):采集一(yī)組觀測(cè)信号序(xu)列Yk(k=1,2,3,，,N),對卡(ka)爾曼濾(lǜ)波參數(shù)B、H、Q、D初始化(huà)，并拟合(he)出R與ƒ的(de)關系式(shì)。
步驟二(èr):首先,根(gēn)據在叠(dié)代中不(bu)斷變化(huà)的狀态(tài)模型頻(pin)率ƒm.對轉(zhuǎn)移矩陣(zhèn)A和觀測(ce)噪聲協(xié)方差R進(jin)行參數(shu)更新;然(rán)後,對觀(guan)測信号(hao)Yj進行狀(zhuàng)态預測(ce),并輸出(chū)最優估(gu)計信号(hao)序列Xk(h=1,2,3,,N)。
步(bù)驟三:通(tōng)過脈沖(chòng)翻轉整(zhěng)形方法(fa)對最優(yōu)估計信(xìn)号序列(lie)進行頻(pin)率計算(suan)。設置翻(fan)轉上阈(yu)值Athr和翻(fān)轉下阈(yu)值-Ar,當信(xìn)号由低(di)向高上(shang)升到Athr時(shí),将輸出(chu)的信号(hào)電平置(zhì)高。當信(xìn)号由高(gao)向低下(xia)降到-Athr時(shí)，将輸出(chū)的信号(hào)電平置(zhì)低,最終(zhōng)輸出脈(mò)沖信号(hao)序列Z(h=1,2,3,,N)。通(tōng)過脈沖(chòng)計數方(fang)法直接(jiē)求出脈(mo)沖信号(hào)Zk的周期(qī)序列T;(i=1,2,3,.,M),計(ji)算出周(zhou)期序列(lie)Ti,的平均(jun)值Tavg,得到(dao)濾波輸(shū)出信号(hào)的平均(jun)頻率ƒout=1/Tavg。
步(bù)驟四:輸(shu)出信号(hào)頻率ƒout和(he)狀态模(mó)型頻率(lü)ƒm若滿足(zú)Iƒout-ƒmI≤ƒmx1%,則跳轉(zhuǎn)到步驟(zhòu)五。若lƒout-ƒm|>ƒmx1%,且(qiě)ƒout≥ƒmin,則令fm=fe,并(bìng)跳轉到(dào)步驟二(èr);否則,應(yīng)停止搜(sou)索并保(bao)持輸出(chu)上一輪(lun)信号處(chù)理得到(dào)的渦街(jiē)信号頻(pin)率,跳轉(zhuan)到步驟(zhou)一。
步驟(zhou)五:輸出(chū)信号幅(fú)值Aout和拟(nǐ)合的渦(wo)街信号(hào)幅值α的(de)關系若(ruo)滿足|Aout-αl<αx10%,則(ze)判斷爲(wèi)渦街信(xìn)号頻率(lü)輸出頻(pin)率ƒout,并跳(tiao)轉到步(bu)驟一;若(ruò)|Aout-α|≥αX10%,則認爲(wèi)是周期(qi)振動噪(zao)聲頻率(lü),跳轉到(dao)步驟六(liu)。
步驟六(liu):令ƒm=ƒmin9時,繼(ji)續向下(xia)搜索渦(wo)街信号(hào)頻率。當(dang)ƒm≥ƒmin時,跳轉(zhuǎn)到步驟(zhòu)二。若輸(shu)出頻率(lǜ)ƒout,仍等于(yu)噪聲頻(pín)率,則重(zhòng)複步驟(zhou)六;否則(zé)跳轉到(dao)步驟四(si)。當ƒm<ƒmin時,應(yīng)停止搜(sou)索并輸(shū)出上一(yī)次正确(que)的渦街(jiē)信号頻(pín)率,并跳(tiao)轉到步(bu)驟一。
3實(shi)驗驗證(zhèng)
　　爲驗證(zhèng)本方法(fa)的有效(xiào)性、測量(liang)精度和(he)抗幹擾(rǎo)性,采用(yòng)仿真信(xin)号和實(shí)流信号(hào)在不同(tong)管徑、不(bu)同介質(zhì)下進行(hang)實驗測(cè)試。
3.1仿真(zhēn)實驗
　　本(běn)文帶有(yǒu)管道噪(zào)聲的渦(wo)街信号(hào)模型是(shi)基于牛(niú)津大學(xue)獲得的(de)渦街信(xìn)号功率(lü)譜密度(dù),Shao等在此(ci)基礎.上(shang)加人時(shi)域波形(xíng)規律和(hé)幅度衰(shuai)減現象(xiang)建立的(de)渦街信(xin)号進行(háng)仿真模(mó)型分析(xī)。數學模(mó)型表達(da)式如式(shi)(26)所示。
 
　　式(shi)中:α0爲渦(wo)街信号(hào)幅值;f爲(wèi)渦街信(xin)号頻率(lü);Kƒ、Kα分别爲(wèi)調頻靈(líng)敏度和(he)調幅靈(líng)敏度,K,爲(wei)渦街信(xìn)号頻率(lǜ)與采樣(yang)頻率的(de)比值,設(shè)Kα=1;δα(t)和δƒ(t)分别(bie)爲高斯(si)白噪聲(shēng)和渦街(jiē)信号幅(fu)度和頻(pín)率的波(bō)動偏差(chà);n(t)爲其他(ta)噪聲幹(gan)擾,包括(kuò)低頻振(zhèn)蕩幹擾(rao)、工頻幹(gan)擾、周期(qi)振動幹(gan)擾和随(suí)機幹擾(rǎo)。
　　在上述(shù)模型中(zhong),加人具(ju)有多個(ge)單自由(yóu)度阻尼(ní)彈性系(xi)統線性(xìng)組合特(te)性的瞬(shùn)态沖擊(ji)振動幹(gan)擾模型(xing),如式(27)所(suǒ)示。
 
　　式中(zhōng):n爲系統(tǒng)的自由(yóu)度,取n=6;ne(t)爲(wèi)高斯白(bái)噪聲;ξi爲(wei)阻尼系(xi)數;ƒi爲振(zhen)動頻率(lǜ);Φi爲初始(shǐ)相位;αi、bi、ξi爲(wèi)常數,取(qǔ)值參考(kao)相關文(wén)獻。
3.1.1本方(fang)法的仿(pang)真驗證(zhèng)
　　首先,驗(yàn)證本方(fang)法對瞬(shùn)态沖擊(jī)的濾波(bō)效果。渦(wo)街信号(hào)仿真模(mo)型的采(cai)樣頻率(lü)爲10kHz、采樣(yàng)時間爲(wèi)6s,加入兩(liang)次瞬态(tài)振動幹(gàn)擾,管道(dao)直徑分(fen)别爲25mm;和(hé)50mm,流體介(jie)質爲氣(qì)體和液(ye)體。以管(guǎn)徑爲25mm、頻(pín)率爲9.54Hz的(de)液體介(jie)質信号(hào)爲例,含(han)有瞬态(tài)沖擊幹(gàn)擾的渦(wō)街流量(liang)信号波(bō)形及其(qí)頻譜圖(tú)如圖3所(suǒ)示經過(guo)本方法(fǎ)處理後(hòu)的波形(xíng)及頻譜(pu)圖如圖(tu)4所示。從(cóng)圖3和圖(tu)4中可以(yi)看出,瞬(shùn)态沖擊(ji)幹擾被(bèi)有效濾(lü)除了。
 
3.1.2本(běn)方法對(duì)比仿真(zhēn)實驗
　　将(jiang)本方法(fa)與傳統(tǒng)卡爾曼(man)濾波方(fāng)法、經驗(yan)模态分(fen)解方法(fǎ)(EMD方法)進(jin)行仿真(zhen)實驗對(dui)比，在液(yè)體介質(zhì)中的仿(pang)真實驗(yan)結果列(lie)于表1,在(zai)氣體介(jie)質中的(de)仿真實(shi)驗結果(guo)列于表(biao)2。
 
 
　　在表1、表(biǎo)2中，實際(jì)頻率是(shì)指模拟(nǐ)渦街信(xin)号的頻(pin)率,相對(duì)誤差是(shì)指實測(ce)頻率與(yǔ)實際頻(pin)率的誤(wu)差絕對(duì)值與實(shi)際頻率(lǜ)的比值(zhí),按式(28)計(jì)算。從中(zhong)可以看(kan)出,本方(fāng)法的測(ce)量相對(dui)誤差小(xiǎo)于傳統(tong)卡爾曼(man)濾波方(fāng)法和EMD方(fāng)法的測(ce)量相對(duì)誤差,在(zai)低流量(liàng)的情況(kuang)下，其測(cè)量低誤(wù)差優勢(shì)更爲明(míng)顯。
 
式中(zhong):Er爲相對(duì)誤差,ƒmea爲(wèi)實測頻(pín)率,ƒa爲實(shí)際頻率(lǜ)。
3.2實流實(shi)驗
　　本文(wen)采用由(yóu)上海質(zhì)量監督(du)檢驗技(ji)術研究(jiū)院提供(gòng)的移動(dòng)式氣體(ti)流量标(biāo)定裝置(zhi)進行氣(qi)體介質(zhi)下的仿(pang)真實驗(yan)該裝置(zhì)由被檢(jiǎn)儀表、标(biāo)準儀表(biǎo)、風機、工(gong)控機、穩(wen)壓箱和(hé)變頻器(qi)組成,其(qí)标定流(liu)量範圍(wei)爲0.5~270m3/h,測量(liàng)相對擴(kuò)展不确(que)定度不(bu)大于0.63%,穩(wen)定性和(he)重複性(xing)均不超(chao)過0.3%。
　　本文(wen)采用由(you)上海質(zhì)量監督(dū)檢驗技(ji)術研究(jiu)院提供(gong)的移動(dong)式液體(tǐ)流量标(biāo)定裝置(zhì)進行液(ye)體介質(zhì)下的仿(páng)真實驗(yàn)。該裝置(zhì)由被檢(jian)儀表、标(biāo)準儀表(biǎo)、水泵、工(gōng)控機、穩(wěn)壓罐和(he)變頻器(qi)組成。标(biao)定裝置(zhi)可提供(gòng)近似穩(wen)定的流(liú)量,通過(guò)标定時(shi)間内的(de)累計流(liú)量可驗(yan)證裝置(zhi)的精度(dù)可達0.001m3/h。
　　實(shi)流實驗(yan)的管道(dào)口徑爲(wei)50mm,流體介(jiè)質爲氣(qi)體和液(ye)體,采樣(yàng)頻率爲(wei)10kHz,采樣時(shi)間爲6s。每(měi)組實驗(yan)選取10個(gè)流量點(dian),主要是(shì)受噪聲(shēng)影響較(jiao)大的低(dī)流速信(xìn)号。表3和(he)表4分别(bie)爲管徑(jìng)爲50mm液體(ti)和50mm氣體(ti)的3種方(fāng)法的處(chu)理結果(guo)，其中實(shi)際頻率(lü)爲标定(dìng)裝置上(shang)标準表(biao)的信号(hào)頻率。
 
 
　　實(shí)流實驗(yan)結果表(biao)明,相比(bi)于其他(tā)兩種方(fang)法,本方(fāng)法.具有(you)更小的(de)誤差。
4結(jie)語
　　本文(wén)提出了(le)一-種基(ji)于渦街(jie)信号模(mó)型的卡(ka)爾曼濾(lü)波的渦(wo)街流量(liang)計信号(hào)處理方(fang)法。首先(xian)分析了(le)卡爾曼(man)濾波算(suan)法的原(yuán)理,利用(yong)微分原(yuan)理和線(xian)性矩陣(zhèn)建立渦(wo)街信号(hao)的線性(xìng)系統模(mó)型。模型(xíng)的初始(shǐ)頻率由(yóu)渦街信(xìn)号的最(zui)大頻率(lǜ)決定,提(tí)高了算(suan)法的計(jì)算效率(lǜ)。而後結(jié)合模糊(hu)搜索和(hé)叠代算(suan)法對卡(ka)爾曼濾(lü)波算法(fǎ)進行改(gǎi)進,通過(guo)叠代搜(sōu)索使濾(lǜ)波結果(guǒ)逐漸接(jie)近渦街(jiē)信号。經(jīng)驗證,循(xun)環叠代(dai)次數一(yi)般在3~10次(cì)之間,複(fu)雜度低(dī)，響應速(su)度快。接(jie)着爲叠(die)代循環(huan)設置終(zhōng)止條件(jiàn)，判斷是(shi)否找到(dào)渦街信(xin)号,并通(tōng)過渦街(jiē)信号的(de)特性設(she)置邊界(jie)條件,防(fang)止叠代(dài)過程發(fā)散。實現(xiàn)了卡爾(ěr)曼濾波(bō)器的自(zi)适應濾(lü)波功能(neng)。最後通(tōng)過仿真(zhen)實驗和(hé)實流實(shí)驗計算(suan)信号頻(pín)率和相(xiang)對誤差(chà),并與傳(chuán)統的卡(ka)爾曼濾(lǜ)波方法(fa)和EMD方法(fa)進行比(bi)較。實驗(yàn)結果表(biǎo)明,與其(qí)他兩種(zhong)方法相(xiang)比,所提(tí)方法具(ju)有測量(liang)精度、抗(kàng)振性。渦(wō)街信号(hao)的幅值(zhí)與頻率(lǜ)的關系(xi)是本文(wén)算法初(chū)始參數(shù)和輸出(chu)條件的(de)設計依(yi)據,其系(xì)數易受(shòu)流體溫(wen)度和探(tàn)頭損耗(hào)的影響(xiǎng)，從而影(yǐng)響算法(fǎ)精度。因(yīn)此,本文(wen)設計的(de)算法适(shi)用于低(di)流體密(mi)度、低腐(fǔ)蝕、低溫(wēn)波動的(de)場合。

本(běn)文來源(yuán)于網絡(luò),如有侵(qīn)權聯系(xì)即删除(chu)！