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

本(ben)文來源(yuán)于網絡(luò),如有侵(qīn)權聯系(xi)即删除(chú)！