陸 星 林建華 楊先軍

（中國核動力研究設(shè)計(jì)院第一研究所，四川 成都 610005）

0 引言

在工業(yè)生產(chǎn)過程中， 運(yùn)行人員通常根據(jù)工藝設(shè)備（如容器、管道、換熱器等）內(nèi)介質(zhì)的壓力、流量、液位、溫濕度等一系列過程參數(shù)變化， 對相應(yīng)泵、 閥門或其他現(xiàn)場執(zhí)行機(jī)構(gòu)進(jìn)行實(shí)時控制， 以保證生產(chǎn)過程按既定目標(biāo)持續(xù)、 安全運(yùn)行。 而準(zhǔn)確掌握過程參數(shù)變化的情況，必須依靠工業(yè)測量儀表獲取。 但在測量過程中，因?yàn)楦鳝h(huán)節(jié)測量儀表在信號轉(zhuǎn)換過程中存在誤差，并且每個環(huán)節(jié)誤差會隨測量系統(tǒng)逐級向下傳遞， 為研究這種誤差傳遞過程， 可通過建立各測量環(huán)節(jié)誤差傳遞數(shù)學(xué)模型進(jìn)行定量分析。

MATLAB 由于其功能強(qiáng)大、庫函數(shù)齊全、代碼簡潔、和多種計(jì)算機(jī)語言或軟件具有接口等特點(diǎn)， 在科學(xué)工程計(jì)算、數(shù)據(jù)處理中應(yīng)用廣泛。 同時，配合MATLAB 的GUIDE 圖形界面編輯工具還可方便地為計(jì)算代碼提供友好的圖形界面[1-2]， 并且代碼編譯后可生成可執(zhí)行文件（.exe 文件）可以完全脫離MATLAB 軟件運(yùn)行，具有良好的可移植性， 便可使其他人員也能對相應(yīng)測量系統(tǒng)進(jìn)行定量分析。

1 誤差分析與模型建立

1.1 測量基本原理

過程測量中常見遇到壓力、液位和流量測量。 按最小測量系統(tǒng)原理如圖1 所示。

主要環(huán)節(jié)由變送器（Transducer）、電流電壓轉(zhuǎn)換組件（I/V Component），以下簡稱I/V 組件、開方運(yùn)算組件（Sqrt. Root Component），以下簡稱Sqrt 組件，和記錄儀（Recording meter）或顯示儀（Display meter）組成。 壓力、差壓、 液位測量 （若采用測量液體產(chǎn)生壓力來實(shí)現(xiàn)液位測量的方式） 均可歸為壓力測量。 其原理是將介質(zhì)（如水）通過儀表引壓管引入壓力（或差壓）變送器表腔內(nèi)，變送器將介質(zhì)壓力P（或壓差ΔP）物理量，轉(zhuǎn)換為儀表通用的4 ～20mA 標(biāo)準(zhǔn)直流電流信號， 該電流信號進(jìn)入電流電壓轉(zhuǎn)換組件轉(zhuǎn)換為兩路獨(dú)立的0 ～10V標(biāo)準(zhǔn)電壓信號， 一路輸出給顯示儀， 一路輸出給記錄儀。 信號傳輸原理流程如圖1 中實(shí)線所示。 流量測量若采用介質(zhì)流過節(jié)流裝置 （如孔板、 文丘里管、 噴嘴等），在節(jié)流裝置前后形成差壓的測量方式。 是通過差壓變送器測量裝置產(chǎn)生的壓差ΔP， 然后利用流量Q和ΔP 平方根成正比的函數(shù)關(guān)系計(jì)算得到。 因此，壓力、液位和流量過程參數(shù)測量實(shí)質(zhì)上是都壓力P（或壓差ΔP）測量，流程基本一致，唯一不同的是流量測量系統(tǒng)中，在I/V 組件之后，還需對電壓信號進(jìn)行開方運(yùn)算，即加入了Sqrt 組件，如圖1 中虛線所示。

1.2 組件信號傳遞函數(shù)

為方便系統(tǒng)誤差分析和算法推導(dǎo)， 先分析理想情況下，測量系統(tǒng)信號輸入/輸出的計(jì)算關(guān)系。

以壓力為例： 若所測量壓力的量程范圍為0 ～Pmax，與變送器輸出標(biāo)準(zhǔn)電流信號IT范圍（4 ～20mA）成線性函數(shù)IT=IT(P)=KIP+CI關(guān)系，其中KI，CI分別為比例系數(shù)和常數(shù)， 根據(jù)待定系數(shù)法， 壓力的最小值0 和最大值Pmax分 別 與4 和20mA 相 對 應(yīng)。 于 是KI，CI可 以 由 式（1）的矩陣確定。

同 理，I/V 組 價(jià)， 將4 ～20mA 標(biāo) 準(zhǔn) 電 流 信 號 轉(zhuǎn) 換 為0～10V 標(biāo)準(zhǔn)電壓信號線性表示為VI/V=V(IT)=KI/VIT+CI/V，其中比例系數(shù)和常數(shù)KI/V，CI/V由式（2）確定。

圖1 測量原理框圖

Sqrt 組 價(jià)， 將I/V 組 價(jià) 輸 出 的0 ～10V 標(biāo) 準(zhǔn) 電 壓 進(jìn)行開方運(yùn)算，0 ～10V 標(biāo)準(zhǔn)電壓信號也可表示為函數(shù)Vsqrt=Ks，其中常數(shù)Ks 是節(jié)流裝置設(shè)計(jì)時固定的一個常數(shù)，為方便本文取。 數(shù)字顯示儀表或記錄儀要將0 ～10V 的標(biāo)準(zhǔn)電壓信號“還原”成被測物理量壓力， 顯示或記錄。 依次類推， 它們的線性函數(shù)分別為：PDout=P(VI/V)=KDVI/V+CD和PRout=P(VI/V)=KRVI/V+CR，且 變 化 范 圍 為：0 ～P max， 由 此， 比 例 系 數(shù) 和 常 數(shù)KD，CD，KR，CR由式（3），（4）給出：

綜上，所有組件輸入與輸出的傳遞函數(shù)為：

1.3 誤差傳遞模型與算法

上述組件傳遞函數(shù)給出的是無誤差理論上的輸入輸出值的對應(yīng)關(guān)系， 但實(shí)際上， 無論是生產(chǎn)工藝還是環(huán)境影響， 各組件當(dāng)輸入一個值時， 輸出的實(shí)際值和理論值都有偏差， 而這個偏差就是該組件的絕對誤差表達(dá)式為：

式中δX，Xexp，Xth，分別為絕對誤差，組件輸出實(shí)測值和理論值。 實(shí)際工作中， 人們最關(guān)注的就是單個環(huán)節(jié)的誤差。 往往單個環(huán)節(jié)的絕對誤差通過實(shí)驗(yàn)室校準(zhǔn)和檢定測量就能得到。 然而關(guān)鍵的是在整個測量系統(tǒng)中， 不僅僅是各個環(huán)節(jié)的組件各自的誤差， 而且各個環(huán)節(jié)的誤差會因?yàn)樾盘柕淖儞Q隨著測量系統(tǒng)逐級往下傳遞。 舉例說明，由2.1 節(jié)測量原理可知，介質(zhì)的壓力通過變送器變換為電流， 此時的變送器輸出的信號絕對誤差為δIT； 而變送器輸出的電流經(jīng)過I/V 組件時，除了I/V 組件本身存在的輸出信號有絕對誤差δVI/V外，還有δIT成分I/V 組件傳遞到I/V 的輸出中，即I/V 輸出信號的誤差既包含本身信號轉(zhuǎn)換的誤差， 還包含前級變送器輸出誤差引入的誤差，以此類推，當(dāng)I/V組件作為下一個環(huán)節(jié)， 顯示儀、 記錄儀或Sqrt 組件的輸入時， 其輸出誤差又稱為了下級設(shè)備輸出誤差的一個分量。 該誤差傳遞原理如圖2 所示。

圖2 誤差傳遞原理圖

根據(jù)上述分析，以及誤差傳遞公式[3]，該測量系統(tǒng)誤差傳遞算法的數(shù)學(xué)模型建立如下： 設(shè)變送器絕對誤差 為δIT，I/V 組 件 絕 對 誤 差 為δVI/V，數(shù)字顯示儀絕對誤差為δPD，I/V 組價(jià)的電流轉(zhuǎn)電壓后信號總誤差為δVI/Vtotal，數(shù)字顯示儀轉(zhuǎn)換為壓力（壓差）后總誤差為δPtotal，根據(jù)間接量誤差傳遞以及誤差合成公式，可得：

流量測量較壓力測量多了開方組件環(huán)節(jié)， 因此誤差傳遞的同理可得：

需要指出的是： 各環(huán)節(jié)自身的絕對誤差， 即δIT，δVI/V，δPD一般由有資質(zhì)的實(shí)驗(yàn)室校準(zhǔn)或檢定得到一般作為已知量，而δVI/Vtotal，δPtotal和δQtotal作 為 未 知 量，需 要通過上述方式得到， 才是整個系統(tǒng)的測量誤差。 若系統(tǒng)更為復(fù)雜情況， 亦可按照上述分析方式逐級推演獲得。

因本文采用 MATLAB 進(jìn)行計(jì)算與程序設(shè)計(jì)，MATLAB 對于矩陣運(yùn)算上具有極大優(yōu)勢，因此，將上述表達(dá)式進(jìn)行變換，可得：

2 計(jì)算程序設(shè)計(jì)

2.1 GUIDE 設(shè)計(jì)

為使其他人員在不需要對代碼進(jìn)行編輯修改的情況下， 就能快速利用各環(huán)節(jié)組件絕對誤差計(jì)算整個系統(tǒng)的誤差。根據(jù)上述數(shù)學(xué)模型的算法，采用MATLAB 的GUIDE 工具編寫了具有GUI 圖形界面的誤差分析程序的，如圖3 所示。

圖3 GUI 圖形界面

2.2 程序功能

“Browse”按鈕，原始數(shù)據(jù)文件瀏覽選擇鍵，利用該鍵可以方便地將事先錄入的原始數(shù)據(jù)EXCEL 表格（支持格式.doc 和.docx） 讀入并顯示讀入進(jìn)度， 數(shù)據(jù)讀入成功后， 將在程序圖形界面的預(yù)設(shè)位置顯示讀入數(shù)據(jù)和整個EXCEL 表格中數(shù)據(jù)的行、列數(shù)，以此判斷初始數(shù)據(jù)輸入的正確性；

“P Transmitter”或“Q Transmitter”為壓力變送器和流量變送器模式， 分別對應(yīng)程序的壓力和流量的算法。 使用人員可根據(jù)所測物理量，適當(dāng)選擇，壓力、差壓和液位選擇 “P Transmitter” 模式，流量選擇 “Q Transmitter”模式；

“Calculate”鍵為程序計(jì)算開始，點(diǎn)選后程序計(jì)算開始，顯示計(jì)算進(jìn)度，計(jì)算結(jié)束后，程序?qū)⑤敵鲇?jì)算結(jié)果至之前EXCEL 表格預(yù)留的結(jié)果欄，并繪制相應(yīng)的函數(shù)圖像。

3 驗(yàn)證舉例

圖4 誤差傳遞輸出曲線

為驗(yàn)證計(jì)算程序功能實(shí)現(xiàn)情況， 壓力和流量測量誤差傳遞分析的對象為一套0 ～5Mpa 的壓力測量系統(tǒng)，和一套0-40t/h-1的差壓流量測量系統(tǒng)，將各環(huán)節(jié)誤差數(shù)據(jù)和系統(tǒng)參數(shù)輸入程序， 輸出結(jié)果如圖4 所示。

4 分析與討論

兩套測量系統(tǒng)輸出曲線自變量和應(yīng)變量分別對應(yīng)線性和平方根關(guān)系， 證明此系統(tǒng)中的誤差傳遞過程未改變測量對應(yīng)關(guān)系。 從壓力曲線可以看出 0Mpa、3.5Mpa 和4.5Mpa 位置處測得的壓力誤差較為明顯，且線性度變差。 而流量測量曲線因?yàn)榱阄挥忻黠@誤差，導(dǎo)致曲線在0-20kpa 范圍內(nèi)有所偏移，曲線其余部分誤差比較小， 分析原因?yàn)樽兯推髁泓c(diǎn)存在較大漂移。

5 結(jié)論

（1） 對誤差傳遞的探討有助于詳細(xì)掌握系統(tǒng)各環(huán)節(jié)誤差分量來源、貢獻(xiàn)以及系統(tǒng)整體誤差的影響；

（2）本文建立誤差傳遞的理論模型，可詳細(xì)掌握有效評估整個測量系統(tǒng)最大誤差變化范圍， 對過程測量系統(tǒng)設(shè)計(jì)維修評估提供有效理論支撐；

（3） 基于誤差傳遞的理論所編寫的 MATLAB GUIDE 程序操作方便， 為其他誤差分析人員在不要編寫代碼情況下， 將參數(shù)輸入程序， 一鍵操作就能給出各測點(diǎn)的誤差， 避免了大量繁瑣的人工計(jì)算， 極大提高了工作效率。