摘要:針對氣體超聲流量(liàng)計
在(zài)測量中存(cún)在回(huí)波信号衰(shuāi)減大(dà)、波形(xíng)易受(shòu)工況(kuàng)影響(xiǎng)的問題,提(tí)出了(le)一種基于(yú)精度高時(shí)差的(de)氣體(tǐ)超(chāo)聲流量(liàng)測量(liàng)方法。該方(fāng)法首先(xiān)通過相(xiàng)似度評估回波(bō)信号(hào),對回波信号特征(zhēng)點(diǎn)進行(háng)準确定(dìng)位(wèi),進而(ér)獲(huò)取飛行時間(jiān)差的粗(cū)測(cè)量值,其次(cì)選取(qǔ)特定(dìng)回波波(bō)形(xíng)進行(háng)互相關法計算(suàn)獲得時差(chà)的細(xì)測量(liàng)值,最後對(duì)兩次測量結果(guǒ)相加得(dé)到精度(dù)高時差(chà),從(cóng)而實現精度高的流(liú)量測量。不同壓力下的(de)聲速測(cè)量實驗(yàn)表明該(gāi)方法(fǎ)在100kpa至500kpa範圍(wéi)内可準确測量(liàng)飛行(háng)時間(jiān)和時(shí)差。氣體(tǐ)流(liú)量計(jì)樣機的流量測量誤差小(xiǎo)于1%,重(zhòng)複性優于0.2%,并在(zài)大流量下(xià)與傳(chuán)統阈(yù)值法相比(bǐ)具有更高(gāo)的正(zhèng)确率和(hé)更優的(de)重複性。
氣體超(chāo)聲流(liú)量計以結(jié)構簡單、壓(yā)損低、精度(dù)高、量(liàng)程寬(kuān)及易(yì)于維(wéi)護等(děng)優(yōu)點(diǎn),成爲天然(rán)氣貿易(yì)中的重(zhòng)要流量儀表。目(mù)前氣體(tǐ)超(chāo)聲流(liú)量計(jì)測量方法(fǎ)多(duō)采用時差法(fǎ),該方法通過超(chāo)聲波(bō)在(zài)管道内順、逆(nì)流傳播的(de)飛行時間及聲速來(lái)計算流量(liàng)。飛行(háng)時間測量(liàng)常采用曲線拟(nǐ)合法、互(hù)相(xiàng)關法和阈(yù)值法,其中(zhōng)曲線拟合法計(jì)算過(guò)程比較複(fú)雜,而阈值法與(yǔ)互相(xiàng)關法(fǎ)的理(lǐ)論成熟,被(bèi)廣泛(fàn)應用于氣體超(chāo)聲流(liú)量測(cè)量中。阈值法通(tōng)過定(dìng)位回波(bō)信号的(de)特征(zhēng)點來測量飛(fēi)行(háng)時間(jiān),然而測量工況的變(biàn)化會使特(tè)征點(diǎn)定位錯誤,導緻(zhì)測量(liàng)結果誤差(chà)偏(piān)大。對此基于(yú)回波信号峰值的比例阈值法(fǎ),該方法(fǎ)通(tōng)過回波峰值(zhí)調整阈(yù)值來定位(wèi)特征(zhēng)點以求(qiú)得時差(chà)。基于分段(duàn)流速的可變(biàn)阈值法(fǎ),通過在(zài)不(bú)同流(liú)速(sù)下設置不同(tóng)阈值對特征點(diǎn)進行定(dìng)位(wèi),進而求得時(shí)差。基于回波極值點的幅(fú)值對阈(yù)值進(jìn)行調(diào)節的自适(shì)應阈值法,該方法使用(yòng)當前工(gōng)況下的回波極值點對阈(yù)值進(jìn)行修正,進而準(zhǔn)确定(dìng)位特征點(diǎn)位置并得(dé)到時(shí)差。上(shàng)述方法均根據不同(tóng)的工(gōng)況對阈值進行調整,以(yǐ)提(tí)高時差測(cè)量(liàng)精(jīng)度,但(dàn)面對(duì)複雜(zá)的測(cè)量環境仍(réng)存在(zài)局限性。互相關法通過将(jiāng)順、逆流回波(bō)信号進(jìn)行互(hù)相關計算(suàn)以得到(dào)時(shí)差值(zhí),可以解(jiě)決由于(yú)波形變化(huà)引起(qǐ)的特征點定位(wèi)錯誤問題。通過(guò)選取各(gè)換能器(qì)靜态下的(de)回波(bō)信号均值作爲互相關(guān)計算的(de)參考信号以提(tí)高(gāo)測量(liàng)結(jié)果的(de)抗幹(gàn)擾能力。提(tí)出了使(shǐ)用實時(shí)動态參考波形(xíng).進行互相(xiàng)關計算的方法,有效(xiào)解決了由環(huán)境因素(sù)導緻相關(guān)性降低(dī)的(de)問題,然而(ér)以上方法(fǎ)均存在(zài)計(jì)算量(liàng)較大(dà)的問(wèn)題(tí)。
針對阈值法(fǎ)與互(hù)相關(guān)法存在的(de)問題,本(běn)文(wén)提出了基(jī)于相似度(dù)和互(hù)相關法的(de)精度(dù)高時差測(cè)量方法(timedifferencemeasurementmethodbasedonsimi-larityandcross-correlation,tdm-sc)。該(gāi)方法(fǎ)通過(guò)回波相似(sì)度評(píng)估,實現特征點的正确定(dìng)位,并結合(hé)傳輸(shū)時差法與(yǔ)互相(xiàng)關法(fǎ)分(fèn)别對時(shí)差進行粗(cū)、細兩(liǎng)次測量,以提高(gāo)其測(cè)量(liàng)精度。
測量原理
1.1時(shí)差法基本(běn)原(yuán)理(lǐ)
時差(chà)法超聲流(liú)量計(jì)的測量原(yuán)理如(rú)圖1所示。超(chāo)聲換能器a,b分(fèn)别(bié)安裝(zhuāng)在流量(liàng)計管道的上(shàng)下遊位置(zhì),超聲波從a傳(chuán)播(bō)到b爲(wèi)順(shùn)流時間(jiān),超聲波從b傳播(bō)到a爲逆(nì)流時間(jiān),流體流(liú)量與順、逆流時間差,的關系如式(shì)(1)所示:
式中:q是管道中(zhōng)氣體瞬時流量(liàng),d爲管(guǎn)道直徑,△t是(shì)時差,c爲(wèi)聲(shēng)速(sù),α是(shì)信号傳播路徑與管道(dào)軸(zhóu)線的(de)夾角。由式(1)可知,時差(chà)測量精度(dù)将直接影響氣(qì)體超聲流(liú)量計流量(liàng)計(jì)算的精(jīng)度。
1.2時差測量方案
時差測量方(fāng)法的原(yuán)理如(rú)圖2所示(shì)。首先(xiān)通過(guò)對采(cǎi)集的回波(bō)信号(hào)與(yǔ)參(cān)考信(xìn)号進行(háng)相(xiàng)似度(dù)計算,獲得回波特征(zhēng)點。其次通過特征(zhēng)點結合(hé)采樣頻率(lǜ)得到“粗”時(shí)差值(zhí);同時(shí)以特(tè)征(zhēng)點作爲起始(shǐ)點來(lái)選取(qǔ)特定的波形數(shù)據,并(bìng)将(jiāng)選取波(bō)形進(jìn)行上(shàng)采樣(yàng)處理,通(tōng)過互相(xiàng)關運算得(dé)到“細(xì)”時差(chà)值(zhí),最終獲(huò)得精度高時差測(cè)量結果。
2基于(yú)相似度的(de)特征點定位(wèi)
由于噪(zào)聲幹擾和測量環境會(huì)使回波信号的幅(fú)值發生變(biàn)化,最終導(dǎo)緻回波信号起(qǐ)始點(diǎn)定位錯誤(wù)。因此需(xū)在回波(bō)信号上找到一(yī)個穩定的特征(zhēng)點,如圖3所示。該(gāi)特征點(diǎn)與(yǔ)回波(bō)起始點之間時(shí)間(jiān)恒定(dìng),通過特(tè)征點結合(hé)采樣(yàng)頻率計算(suàn)得到(dào)順、逆(nì)流的兩個(gè)特征(zhēng)飛行(háng)時間(jiān)tcharacter,将兩(liǎng)者(zhě)相減可(kě)抵消(xiāo)固定(dìng)時延,從而(ér)得(dé)到傳播(bō)時間差值(zhí)。
目前廣泛(fàn)使用的特(tè)征點定位方法(fǎ)是雙阈值法,其(qí)原理(lǐ)如圖4所示(shì)。第一阈(yù)值線電壓值約(yuē)爲(wèi)0.35v,0v幅值(zhí)線作(zuò)爲(wèi)第二阈(yù)值用于(yú)定(dìng)位到(dào)過零(líng)采樣(yàng)點,即(jí)回波特(tè)征點。在(zài)不同(tóng)流量或工(gōng)況下(xià),回波(bō)信号(hào)的幅(fú)值(zhí)特性會發生變(biàn)化(huà)。此時(shí)若采用固(gù)定阈(yù)值來确(què)定回波信号(hào)特征點,會(huì)造成飛行時(shí)間測量存在(zài)數個(gè)周期的誤差。如圖4所(suǒ)示,當環境(jìng)壓力從500kpa變(biàn)化爲(wèi)101kpa時,原(yuán)本通過第一阈值定(dìng)位的第三(sān)個波形會錯誤地定位在第四個波形,上(shàng)引起.測(cè)量(liàng)誤差(chà)。改進的阈值法結合不(bú)同的工況來(lái)對阈(yù)值進(jìn)行(háng)調整(zhěng),然而在(zài)複雜的測(cè)量環(huán)境下(xià),這些(xiē)方法依然存在一定(dìng)局限性。
針對以上問題,本文采(cǎi)用基于(yú)相似(sì)度的回(huí)波特征點(diǎn)定位(wèi)方法來獲(huò)取特征(zhēng)點。首先(xiān)通過(guò)0v幅值線獲(huò)得回(huí)波信(xìn)号的(de)多個過零采樣點,以作爲“備選”特征點,即圖5方(fāng)框内采樣點。随(suí)後将采集得(dé)到(dào)的回(huí)波信(xìn)号峰(fēng)值電壓與标準(zhǔn)工況(kuàng)下的峰值(zhí)電壓進(jìn)行相似(sì)度計(jì)算,從而正(zhèng)确定位到(dào)回波特征(zhēng)點。
回波信号相似度評估選(xuǎn)擇标準工況(kuàng)下的(de)參考回波信号i和實測回(huí)波信号j作(zuò)爲相(xiàng)似估(gū)計對象。參(cān)考信号第2,3,4峰值(zhí)電壓值(zhí)與實測(cè)信号(hào)各個峰值(zhí)電壓(yā)值xi,xj爲特征(zhēng)參數,參數數量(liàng)n取3。通過計算,選(xuǎn)取與參考回波信号歐氏(shì)距離最小(xiǎo)的(de)一組實(shí)測回(huí)波信号峰(fēng)值,即(jí)實(shí)際回波(bō)信号(hào)的第2,3,4波峰(fēng)值,将特征(zhēng)點準确定位到實際(jì)回波(bō)信号第2波後的過零點,即(jí)圖5中點p2。
3精度高(gāo)時差的測量
3.1粗(cū)時差測(cè)量
激勵(lì)信号(hào)驅(qū)動(dòng)超聲(shēng)換能器發射聲(shēng)波後(hòu),采樣電(diàn)路(lù)開始進行(háng)回波(bō)信号采(cǎi)集。超聲(shēng)波順、逆流(liú)傳播(bō)的(de)飛(fēi)行時間tui和(hé)tdi通過其對應采樣點(diǎn)數(shù)n與采樣(yàng)間隔t的(de)乘積表(biǎo)示,求得粗時差(chà)值,計算如式(3)所(suǒ)示:
式中:n1和n2分别(bié)爲順、逆(nì)流(liú)下回(huí)波信号特征點(diǎn)對應(yīng)的采樣點(diǎn)數。
3.2細(xì)時差測量(liàng)
3.2.1波形(xíng)選取與(yǔ)上(shàng)采樣處理(lǐ)
針對(duì)互相關計(jì)算過(guò)程中(zhōng)運(yùn)算量較(jiào)大的問題(tí),選擇(zé)回波(bō)信(xìn)号特(tè)征(zhēng)點後三個(gè)周期(qī)的采樣點(diǎn)作爲(wèi)待處理數據以降(jiàng)低運算量。具(jù)體信号(hào)如(rú)圖6虛(xū)線方框内點所(suǒ)示。
對選取信号(hào)進行上采樣處(chù)理來提高采樣(yàng)率。上采樣處理(lǐ)包括數據的插值和(hé)低通濾波(bō)兩個步驟(zhòu)。首先将采集(jí)到的數(shù)據量爲n的原始(shǐ)信号(hào)x[n]中每兩個(gè)采樣點之(zhī)間插(chā)人l-1個(gè)零值得到(dào)信号(hào)xu[n],如式(4)所示(shì):
爲了更好(hǎo)地觀(guān)察信(xìn)号處(chù)理前(qián)後的頻(pín)率(lǜ)特性(xìng),通過式(5)、式(shì)(6)将信号x[n]、xu[n]轉移到頻域(yù),如式(7)所示(shì),并得到幅度譜(pǔ)圖,如(rú)圖7(a)、圖7(b)所示。
對于(yú)因子(zǐ)爲l的插零(líng)擴展,相較于圖(tú)7(a),插值後的(de)信号(hào)在基(jī)帶(dài)_上(shàng)有l-1個(gè)額外的原信号(hào)譜鏡像産生。随(suí)後通過低(dī)通濾(lǜ)波濾除這(zhè)l-1個鏡(jìng)像,等同(tóng)于(yú)将内(nèi)插樣本(běn)值“填入(rù)”到xu[n]中(zhōng)的零樣本(běn),上,實現原(yuán)采集(jí)信号(hào)x[n]的(de)上(shàng)采樣(yàng)處理。
設計(jì)的低通濾(lǜ)波器的頻(pín)域表達爲(wèi)式(8):
當(dāng)c=l時滿(mǎn)足零初始(shǐ)條件,濾波器的(de)頻域表達如式(10)所示:
采樣(yàng)信号x[n]與經(jīng)過l=20進(jìn)行上(shàng)采樣處理後信号xu[n]的(de)數據與幅(fú)度譜圖如圖8(a)和(hé)圖8(b)所(suǒ)示(shì),結果表(biǎo)明上(shàng)采樣處理後的信(xìn)号采樣(yàng)率增大(dà)了(le)20倍,同(tóng)時處理後(hòu)的數據曲線(xiàn)光滑,證(zhèng)明上采樣處理(lǐ)符合預期(qī)效果(guǒ)。
3.2.2互相關(guān)計算
将(jiāng)順、逆(nì)流(liú)回(huí)波信号的原(yuán)始(shǐ)采樣(yàng)數據(jù)進行(háng)上采樣處理得到xu(n)、yu(n)後(hòu),通過離散互(hù)相關運算式(11)得到(dào)互相關(guān)函數rxy(m):
式中:m=(-n+1,n-1),n爲回波數據(jù)的信号長度。如圖9所示,互相關(guān)函數rxy(m)的峰(fēng)值b所(suǒ)對應(yīng)的時間值(zhí)即爲兩信(xìn)号時差。爲(wèi)了進一步(bù)提高時差精(jīng)度(dù),選取互相(xiàng)關函(hán)數rxy(m)中峰(fēng)值(zhí)處的三個(gè)最高(gāo)點a、b、c進(jìn)行(háng)曲線(xiàn)拟(nǐ)合以(yǐ)得到(dào)更(gèng)精确的(de)峰值d(max,ymax)。
通過式(12)得到xmax對應的細時(shí)差值(zhí)△tcorr,其(qí)中t爲采(cǎi)樣間隔。
4系統實(shí)現
4.1硬件設計
采用msp430f6638芯(xīn)片作爲核(hé)心控制(zhì)單(dān)元,負(fù)責整(zhěng)個(gè)測量過(guò)程(chéng)中(zhōng)時序(xù)和所(suǒ)屬電路(lù)的控制。fpga模塊用以産生(shēng)驅動電路(lù)的觸(chù)發脈沖以(yǐ)及對采樣數據進行(háng)實時獲取與存(cún)儲,如圖10所示。包(bāo)括兩路(lù)激勵電(diàn)路、切(qiē)換(huàn)接收電(diàn)路、回波(bō)信号處(chù)理電路(濾波(bō)放(fàng)大電路、回(huí)波到(dào)達電(diàn)路、峰值檢測電路)和信号采樣(yàng)電路(lù)等。激(jī)勵電路将(jiāng)觸發脈沖(chòng)進行推挽放大後(hòu)輸人到超聲波(bō)換(huàn)能器(qì)并(bìng)使其發射超(chāo)聲波(bō)。回波信号(hào)接收(shōu)後經過回(huí)波到(dào)達探測電路産(chǎn)生一個回波到(dào)達(dá)信号再(zài)輸入到(dào)單(dān)片機(jī)。msp430單片機通(tōng)過内(nèi)部ad對經(jīng)過峰值(zhí)檢測(cè)電路(lù)的回波(bō)信号(hào)進行采(cǎi)集,獲得回波(bō)的(de)最大峰值。放大(dà)後的(de)回波信号(hào)由fpga配合高速ad以(yǐ)及(jí)ram進(jìn)行(háng)模(mó)數轉(zhuǎn)換和數據存(cún)儲(chǔ),采集(jí)到的(de)數據(jù)通過485通(tōng)信電路(lù)傳輸(shū)到計算機(jī)進(jìn)行數據(jù)處理(lǐ)。
電路(lù)采用的超聲換能(néng)器中心(xīn)頻率爲200khz,驅動信(xìn)号幅值爲20v。采樣電路中高速采(cǎi)集芯片選用ad9237-40,采樣頻率設定爲(wèi)5mhz。
4.2軟件設計(jì)
軟件設計包含msp430程序(xù)和matlab程序兩(liǎng)個部分,如圖(tú)11所示(shì)。
msp430程序流程如下所述。系(xì)統(tǒng)初上電(diàn)後,msp430f6638将(jiāng)對(duì)i0口、定時器及(jí)fpga模(mó)塊等各(gè).個參數進行初(chū)始化并(bìng)進人低(dī)功耗模式(shì)。定時(shí)器達(dá)到0.5s時(shí),微處(chù)理器(qì)控制fpga芯片(piàn)産生(shēng)激勵信号(hào)輸人到指(zhǐ)定的發射(shè)換(huàn)能器中。當單片(piàn)機(jī)接收(shōu)到回波到達信(xìn)号後,,微控制(zhì)器(qì)使能fpga對處理後(hòu)的回波信(xìn)号進(jìn)行采(cǎi)樣并(bìng)存儲在fpga的ram中(zhōng),同時開(kāi)啓單片(piàn)機内部(bù)ad對回波最(zuì)大峰(fēng)值電壓進行采(cǎi)集。随(suí)後(hòu),通(tōng)過上(shàng)位(wèi)機通訊(xùn)将采(cǎi)集到的回波數據傳輸到(dào)matlab程序(xù)。matlab程序(xù)首先根據(jù)回波相似度計(jì)算定位(wèi)到回波信号(hào)的(de)特(tè)征(zhēng)點(diǎn),其次以特征點(diǎn)爲基(jī)礎(chǔ)結合采(cǎi)樣頻率和互相關法(fǎ)得到精度(dù)高的(de)飛行時間差以及實(shí)時聲(shēng)速值,利用時差(chà)法計(jì)算式(1)得到(dào)瞬時流量值。
5實(shí)驗驗(yàn)證
爲評估(gū)方法(fǎ)法的(de)有效(xiào)性,采用(yòng)壓力實(shí)驗驗證時間差(chà)測量(liàng)的(de)穩定性(xìng),進而通過流量(liàng)實驗驗證整體(tǐ)算法的(de)精度。
5.1壓(yā)力實驗研究
裝(zhuāng)置如(rú)圖12所(suǒ)示,包括氮(dàn)氣(qì)鋼瓶和(hé)密封管路裝置(zhì)等。選擇(zé)101kpa、200kpa、300kpa、400kpa及500kpa五個(gè)壓力(lì)點進(jìn)行相(xiàng)關的壓(yā)力。
采用(yòng)本文的(de)信号(hào)處(chù)理方法(fǎ)和(hé)基于tdc-gp22測量(liàng)模塊的傳(chuán)統雙阈值法時(shí)差測量(liàng)方(fāng)法(timedifferencemeasurementmethodbasedontdc-gp22moduleofdoublethresholdmethod,tdm-dt)進行對比(bǐ)。由于在(zài)測量過(guò)程中時差值難(nán)以直觀(guān)表示,而聲(shēng)速測(cè)量與(yǔ)時差測量(liàng)均以飛(fēi)行(háng)時間爲基礎,因此在(zài)各個壓(yā)力(lì)下比(bǐ)較兩種方法測(cè)量得到的(de)聲速(sù)值與理論聲速(sù)值來(lái)間接(jiē)驗證(zhèng)測量(liàng)的(de)穩定性(xìng),結果如表(biǎo)1所示(shì)。
由表1可知,使用(yòng)基于回波(bō)相似(sì)度進(jìn)行(háng)特征點(diǎn)定位的(de)方(fāng)法測(cè)得的5個(gè)壓力試(shì)驗點下聲(shēng)速值(zhí)均與理論聲(shēng)速(sù)吻合,最大誤差(chà)僅爲-0.13m/s。而傳統雙(shuāng)阈值法計(jì)算得(dé)到的(de)聲速(sù)在101.9kpa、203.2kpa及(jí)305.5kpa下與(yǔ)理(lǐ)論聲速(sù)吻合,但在(zài)405.2kpa壓力(lì)下與理論聲速(sù)産生明顯偏差,與此同時(shí)壓力(lì)越大,偏(piān)差數值(zhí)越大。而在509.5kpa下,聲(shēng)速測量值與理(lǐ)論聲速差值高達7.89m/s。實驗結果證(zhèng)明基于回波相(xiàng)似度的特征點(diǎn)定位(wèi)信号處理方法(fǎ)能在不同(tóng)壓力下實現飛行時(shí)間(jiān)差(chà)測量(liàng)的正(zhèng)确率(lǜ)。
5.2流量實驗研(yán)究(jiū)
選用(yòng)圖13所示精(jīng)度等級爲(wèi)0.25級的lqb-1000臨界(jiè)流文丘裏音速噴嘴校準裝置(zhì),采用管(guǎn)徑(jìng)爲50mm的(de)氣體超聲流量測量系統樣機(jī),流量(liàng)範圍爲(wèi)2m'/h~160m'/h。根(gēn)據超聲流(liú)量計(jì)檢定(dìng)規(guī)程(chéng)《jjg1030-2007超聲(shēng)流量計》,選(xuǎn)擇分界流量(liàng)點爲16m2/h。各(gè)個流量(liàng)檢(jiǎn)定(dìng)點(diǎn)爲qmin、qt、0.25qmax、0.4qmax、0.7qmax,和qmax,每個流量(liàng)點測量90s。将(jiāng)測量得到(dào)的流量值(zhí)和(hé)标準裝(zhuāng)置的(de)平均(jun1)流量值進(jìn)行比較,計(jì)算誤(wù)差并進(jìn)行(háng)三次(cì)實驗來得到重複性。基于tdm-sc與tdm-dt兩(liǎng)種方法的測量結果如(rú)表(biǎo)2所示(shì)。
表2數據(jù)表(biǎo)明,基(jī)于tdm-sc的(de)氣體超聲(shēng)流量測量(liàng)系統(tǒng)測量誤差小于1%,重複性(xìng)優(yōu)于0.2%,符合一級(jí)表的要(yào)求。同時(shí)在大流量下(xià),方(fāng)法依然(rán)能保(bǎo)持(chí)低于1%的測量(liàng)誤差和良好的重複(fú)性。
6結論
提(tí)出(chū)了基于(yú)精度高時差的(de)氣體(tǐ)超聲流量測量方法(fǎ),該方法通過回波相(xiàng)似度(dù)評估對回波特(tè)征點進行(háng)準确定位(wèi),在特征點基礎(chǔ)上結(jié)合傳(chuán)輸時(shí)間法(fǎ)與互(hù)相關法對時差(chà)進行粗、細兩次(cì)測量以得(dé)到準确的(de)時差值,最終實現精(jīng)度高(gāo)的(de)流(liú)量測(cè)量。 結(jié)果(guǒ)表明,該(gāi)方法(fǎ)在100kpa至500kpa的(de)壓力下能對時差(chà)進行準确測量(liàng)。系統(tǒng)樣機(jī)的流量測(cè)量精(jīng)度滿(mǎn)足1級(jí)精(jīng)度(dù)的要(yào)求,并在大流量下測(cè)量誤差和(hé)重複性(xìng)優(yōu)于傳(chuán)統雙阈(yù)值(zhí)法。
本(běn)文來源于網絡(luò),如有侵權聯系(xì)即删除!