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

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