趙久強(qiáng) 王震洲

摘要：針對(duì)溫度會(huì)影響紅外CO2傳感器的輸出電壓，造成對(duì)CO2的濃度檢測(cè)誤差較大的問(wèn)題，提出了一種基于L-M貝葉斯正則化BP神經(jīng)網(wǎng)絡(luò)的溫度補(bǔ)償方法。實(shí)驗(yàn)中將傳感器輸出電壓比和溫度作為神經(jīng)網(wǎng)絡(luò)的輸入，CO2濃度作為神經(jīng)網(wǎng)絡(luò)的輸出，并通過(guò)L-M算法和貝葉斯正則化對(duì)神經(jīng)網(wǎng)絡(luò)進(jìn)行優(yōu)化。經(jīng)過(guò)實(shí)驗(yàn)仿真證明，在溫度補(bǔ)償后紅外CO2傳感器測(cè)量輸出的濃度值最大相對(duì)誤差為4.557 8%，具有較高的精確度。因此L-M貝葉斯正則化BP神經(jīng)網(wǎng)絡(luò)能對(duì)紅外CO2傳感器進(jìn)行有效的溫度補(bǔ)償，可為相關(guān)紅外傳感器儀器的改進(jìn)提供參考。

關(guān)鍵詞：計(jì)算機(jī)神經(jīng)網(wǎng)絡(luò)；紅外CO2傳感器；BP神經(jīng)網(wǎng)絡(luò)；L-M算法；貝葉斯正則化；溫度補(bǔ)償

中圖分類號(hào)：TP183文獻(xiàn)標(biāo)志碼：Adoi： 10.7535/hbgykj.2018yx04008

當(dāng)前CO2的排放量不斷增加，加劇了溫室效應(yīng)。能夠?qū)崟r(shí)準(zhǔn)確地檢測(cè)CO2濃度已成為研究熱點(diǎn)。相對(duì)于其他類型的CO2傳感器，紅外吸收型CO2傳感器因靈敏度高、響應(yīng)速度快、選擇性好、無(wú)毒無(wú)害等優(yōu)點(diǎn)，在許多領(lǐng)域得到廣泛應(yīng)用，已成為檢測(cè)CO2氣體濃度的常用傳感器[1-2]。在實(shí)際應(yīng)用中，紅外CO2傳感器測(cè)量的精度易受到溫度的影響[3-5]，因此要對(duì)其進(jìn)行溫度補(bǔ)償。

第4期趙久強(qiáng)，等：L-M貝葉斯正則化BP神經(jīng)網(wǎng)絡(luò)在紅外CO2傳感器的應(yīng)用河北工業(yè)科技第35卷 根據(jù)文獻(xiàn)[6—13]利用BP神經(jīng)網(wǎng)絡(luò)可以不問(wèn)函數(shù)的具體性質(zhì)，只要有足夠的訓(xùn)練樣本，通過(guò)訓(xùn)練和調(diào)整，就可以實(shí)現(xiàn)對(duì)任意函數(shù)逼近，來(lái)建立非線性溫度補(bǔ)償網(wǎng)絡(luò)模型，再結(jié)合L-M算法和貝葉斯正則化算法來(lái)提高神經(jīng)網(wǎng)絡(luò)的推廣能力。通過(guò)實(shí)驗(yàn)證明，該方法對(duì)紅外CO2傳感器能實(shí)現(xiàn)有效的溫度補(bǔ)償。

1紅外CO2傳感器原理

當(dāng)CO2氣體受到紅外光照射時(shí)，CO2氣體分子吸收某些頻率的輻射，使對(duì)應(yīng)于這些吸收區(qū)域的透射光強(qiáng)度減弱，其吸收關(guān)系可由Lambert-Beer表示[14]：

lnI0（λ）I（λ）=klC，（1）

式中：C為CO2濃度；l為介質(zhì)厚度；k為吸收系數(shù)；I0為入射光強(qiáng)；I為透射光強(qiáng)。

實(shí)驗(yàn)中采用的紅外光源為IRL715，紅外傳感器為TPS2534。TPS2534為雙通道探測(cè)器，在探測(cè)器前端分別放置著4.26 μm和4.00 μm的濾光片頭，讓相對(duì)應(yīng)波長(zhǎng)的紅外光通過(guò)，其中4.26 μm波長(zhǎng)的紅外光為CO2吸收紅外光，而4.00 μm波長(zhǎng)的紅外光為參考紅外光，CO2對(duì)其沒(méi)有明顯吸收。經(jīng)過(guò)濾波后，紅外光由探測(cè)器轉(zhuǎn)換為相對(duì)應(yīng)的電壓。其輸出電壓與CO2濃度關(guān)系的推導(dǎo)公式如下。

設(shè)吸收紅外光和參考紅外光的入射光強(qiáng)分別為IC0和IR0，二者經(jīng)過(guò)CO2后的光強(qiáng)為IC和IR，由式（1）可知：

IC=IC0e-klC。（2）

而CO2對(duì)參考紅外光沒(méi)有吸收，入射光強(qiáng)和光強(qiáng)可以近似相等：

IR≈IR0。（3）

將式（2）和式（3）相除后可得

lnICIR=-klC+lnIC0IR0。（4）

于是求得CO2的濃度為

C=-1kllnICIR+1kllnIC0IR0。（5）

而吸收紅外光強(qiáng)和參考紅外光強(qiáng)與紅外探測(cè)器輸出的電壓關(guān)系為

UC=KCIC，（6）

UR=KRIR。（7）

式中KC和KR與濾光片的透光系數(shù)和探測(cè)器靈敏度有關(guān)，為與紅外傳感器相關(guān)的常量。則濃度與紅外傳感器輸出電壓比UC/UR的關(guān)系為

C=-1kllnKRUCKCUR+1kllnIC0IR0。（8）

紅外CO2傳感器的光路結(jié)構(gòu)如圖1所示。

式（8）中吸收系數(shù)k是一個(gè)復(fù)雜的變量，它不僅與氣體種類、入射光有關(guān)，還和環(huán)境溫度和氣壓有關(guān)。在一般的應(yīng)用中，環(huán)境壓強(qiáng)基本保持不變，而環(huán)境溫度變化較大，并且環(huán)境溫度不僅影響吸收系數(shù)，還會(huì)影響紅外光源和紅外探測(cè)器，對(duì)測(cè)量精度造成影響。

2測(cè)量系統(tǒng)

測(cè)量系統(tǒng)框圖如圖2所示，單片機(jī)輸出脈沖信號(hào)，通過(guò)光源驅(qū)動(dòng)模塊來(lái)控制紅外光源，實(shí)現(xiàn)對(duì)紅外光源IRL715的調(diào)制。紅外光通過(guò)氣室，在被檢測(cè)氣體CO2吸收后，到達(dá)紅外探測(cè)器TPS2534。紅外探測(cè)器將紅外信號(hào)裝換為電信號(hào)，緊接著經(jīng)過(guò)信號(hào)調(diào)理模塊，進(jìn)行放大濾波得到穩(wěn)定的電信號(hào)。調(diào)理后電信號(hào)經(jīng)過(guò)A/D轉(zhuǎn)換模塊，將模擬信號(hào)轉(zhuǎn)換為數(shù)字信號(hào)，以便進(jìn)行數(shù)據(jù)處理。最后單片機(jī)將采集到的紅外CO2傳感器數(shù)據(jù)和溫度數(shù)據(jù)通過(guò)串口通信傳輸?shù)缴衔粰C(jī)，在上位機(jī)進(jìn)行數(shù)據(jù)的分析處理。

3BP神經(jīng)網(wǎng)絡(luò)優(yōu)化

3.1BP神經(jīng)網(wǎng)絡(luò)

BP神經(jīng)網(wǎng)絡(luò)包含輸入層、隱藏層和輸出層，其中隱藏層可以為一層或多層。各層由神經(jīng)元構(gòu)成，層與層之間的神經(jīng)元通過(guò)權(quán)重全連接，同層神經(jīng)元沒(méi)有連接。BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練過(guò)程中包含信號(hào)正向傳播，而誤差后向傳播，通過(guò)調(diào)節(jié)網(wǎng)絡(luò)權(quán)重，最終使輸出和期望接近。BP網(wǎng)絡(luò)的具體算法參見(jiàn)文獻(xiàn)[15]。

然而在BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練過(guò)程中會(huì)出現(xiàn)過(guò)擬合的問(wèn)題，導(dǎo)致對(duì)訓(xùn)練集以外的新鮮樣本不能得到適當(dāng)?shù)妮敵?，降低泛化能力?/p>

而正則化是改善BP神經(jīng)網(wǎng)絡(luò)泛化能力的有效方法。正則化是在網(wǎng)絡(luò)均方誤差性能函數(shù)的基礎(chǔ)上，引入網(wǎng)絡(luò)權(quán)值的均方誤差，實(shí)現(xiàn)對(duì)性能函數(shù)的修正，調(diào)高泛化能力。網(wǎng)絡(luò)均方誤差函數(shù)：ED=1n∑ni=1（ti-xi）2。（9）網(wǎng)絡(luò)所有權(quán)值均方誤差：EW=1m∑mj=1w2j。（10）修正后的性能函數(shù)：f（W）=αEW+βED。（11）α，β為性能函數(shù)參數(shù)，參數(shù)的大小決定了網(wǎng)絡(luò)訓(xùn)練的效果。當(dāng)α越大于β，網(wǎng)絡(luò)的輸出會(huì)更加平滑，但訓(xùn)練集可能欠擬合；當(dāng)α越小于β，網(wǎng)絡(luò)誤差會(huì)更小，但過(guò)擬合可能還會(huì)存在。因此需要保證α，β為最優(yōu)參數(shù)，才能使網(wǎng)絡(luò)性能最佳。常規(guī)的正則化方法很難確定α，β，而貝葉斯正則化在網(wǎng)絡(luò)訓(xùn)練過(guò)程中可以完成對(duì)α，β自適應(yīng)，使其到達(dá)最佳。

3.2貝葉斯正則化和L-M 算法

貝葉斯正則化的目的是獲取性能最優(yōu)、參數(shù)最小的神經(jīng)網(wǎng)絡(luò)，從而使過(guò)擬合的問(wèn)題得到有效控制。根據(jù)文獻(xiàn)[16—20]，假設(shè)樣本集D中的噪聲和權(quán)重向量W服從高斯分布，則權(quán)值的后驗(yàn)概率密度函數(shù)與優(yōu)化后的性能函數(shù)關(guān)系為

P（w|D，α，β，M）=

1ZD（β）ZW（α）×exp[-（αEW+βED）]P（D|α，β，M）=

exp（-f（W））Z（α，β） 。（12）

式中：P（D|α，β，M）為數(shù)據(jù)集的似然函數(shù)；

ZD（β）=（πβ）n2；ZW（α）=（πα）n2。

最優(yōu)的權(quán)重向量的后驗(yàn)概率密度P（w|D，α，β，M）應(yīng)為最大，則在Z（α，β）一定的情況下，性能函數(shù)f（W）應(yīng)為最小。

對(duì)性能函數(shù)f（W）在梯度為零的最小值點(diǎn)WMP附近進(jìn)行泰勒級(jí)數(shù)展開(kāi)，忽略高次項(xiàng)，并根據(jù)貝葉斯規(guī)則，在WMP優(yōu)化的正則化參數(shù)可寫為α=γ2EW（WMP），（13）

β=N-γ2ED（WMP），（14）

γ=n-α·Trα（H-1），（15）

H=βΔ2ED+αΔ2EW。（16）式中：γ為神經(jīng)網(wǎng)絡(luò)有效參數(shù)個(gè)數(shù)；n為神經(jīng)網(wǎng)絡(luò)參數(shù)總個(gè)數(shù)；H為f（W）最小點(diǎn)處的Hessian矩陣。

在貝葉斯正則化過(guò)程中，參數(shù)的確定需要計(jì)算f（W）在最小點(diǎn)WMP的Hessian矩陣，這里通過(guò)使用L-M算法來(lái)對(duì)Hessian矩陣近似求解。H≈2βJTJ+2αIn。（17）式中：J為雅可比矩陣；In為單位矩陣。

在每次迭代中，神經(jīng)網(wǎng)絡(luò)有效參數(shù)個(gè)數(shù)γ由當(dāng)前的正則參數(shù)α，β所確定，而網(wǎng)絡(luò)的性能函數(shù)通過(guò)L-M算法達(dá)到最小后，通過(guò)式（10）和式（11）得到新的正則參數(shù)α，β。隨著性能函數(shù)參數(shù)的每次近似，性能函數(shù)將會(huì)改變，最小點(diǎn)將不斷移動(dòng)，從而自適應(yīng)調(diào)節(jié)α，β。

在訓(xùn)練過(guò)程中，網(wǎng)絡(luò)權(quán)值的迭代公式為

Wh+1=Wh-[J（Wh）TJ（Wh）+μI]-1J（Wh）Te。（18）

式中：Wh為第h次迭代的網(wǎng)絡(luò)權(quán)值向量；μ為自適應(yīng)的標(biāo)量；e為誤差矩陣。當(dāng)μ很小時(shí)，L-M算法變?yōu)榕ｎD法，μ很大時(shí)，則為梯度下降法。

4實(shí)驗(yàn)和仿真

4.1紅外CO2傳感器標(biāo)定

首先對(duì)紅外CO2傳感器進(jìn)行數(shù)據(jù)標(biāo)定，利用靜態(tài)配氣法，在標(biāo)準(zhǔn)大氣壓，室溫為25 ℃的環(huán)境下，通過(guò)配置14種不同濃度的CO2氣體樣品，對(duì)紅外傳感器進(jìn)行多濃度標(biāo)定，其中標(biāo)準(zhǔn)CO2氣體濃度的測(cè)定采用比較先進(jìn)的CO2氣體濃度測(cè)量?jī)x器Model 400。數(shù)據(jù)的擬合采用origin來(lái)完成，擬合的結(jié)果如圖3所示。

4.2溫度補(bǔ)償實(shí)驗(yàn)

實(shí)驗(yàn)將紅外CO2測(cè)量裝置置于溫度箱內(nèi)，溫度變化范圍為15～35 ℃，每5 ℃為一個(gè)溫度實(shí)驗(yàn)點(diǎn)，往箱內(nèi)分別通入N2和250×10-6，500×10-6，750×10-6，1 000×10-6，1 500×10-6，2 000×10-6，2 500×10-6，3 000×10-6的標(biāo)準(zhǔn)CO2氣體，待氣體均勻后，測(cè)量計(jì)算傳感器相應(yīng)的輸出電壓比UC/UR，電壓值采樣10次，這樣得到450組樣本，將405組用作訓(xùn)練樣本，45組用作測(cè)試樣本。圖4為測(cè)試樣本數(shù)據(jù)，在不同溫度下，濃度與電壓比的分布圖。從圖4可以明顯看出，在不同溫度、相同濃度下，UCUR值各不相同。

將不同溫度下的UC/UR代入擬合的曲線中，得到測(cè)量值如表1所示。其中最大的相對(duì)誤差為49.6%。

實(shí)驗(yàn)中BP神經(jīng)網(wǎng)絡(luò)采用三層結(jié)構(gòu)，輸入層為2個(gè)節(jié)點(diǎn)，隱藏層的節(jié)點(diǎn)通過(guò)訓(xùn)練比較選用9個(gè)，輸出層節(jié)點(diǎn)為1個(gè)。神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)精度為10-4，最大訓(xùn)練次數(shù)為500次。

表2為測(cè)試樣本的輸出數(shù)據(jù)，圖5為測(cè)試樣本的輸出與實(shí)際值的比較。從圖中可以看出，神經(jīng)網(wǎng)絡(luò)的輸出值能很好地逼近實(shí)際值。在不同溫度下，對(duì)CO2濃度測(cè)量的相對(duì)誤差都在5%以內(nèi)，其中最大的相對(duì)誤差為4.557 8%。

5結(jié)語(yǔ)

針對(duì)溫度影響紅外CO2傳感器測(cè)量的問(wèn)題，提出了一種基于L-M貝葉斯正則化BP神經(jīng)網(wǎng)絡(luò)來(lái)對(duì)其進(jìn)行溫度補(bǔ)償?shù)姆椒?。通過(guò)實(shí)驗(yàn)仿真，在未經(jīng)過(guò)溫度補(bǔ)償?shù)那闆r下，傳感器測(cè)量出的濃度值的最大相對(duì)誤差為49.6%，經(jīng)過(guò)溫度補(bǔ)償后，測(cè)量出的濃度值最大相對(duì)誤差為4.557 8%，證明此方法具有較高的精確度和一定的實(shí)用性。雖然目前關(guān)于神經(jīng)網(wǎng)絡(luò)的理論還不完善，并且獲取樣本和訓(xùn)練神經(jīng)網(wǎng)絡(luò)的過(guò)程還比較繁瑣，但隨著相關(guān)技術(shù)的發(fā)展，神經(jīng)網(wǎng)絡(luò)會(huì)有廣泛的應(yīng)用前景。

參考文獻(xiàn)/References：

[1]張廣軍，武曉利.新型高性能紅外二氧化碳傳感器[J].紅外與激光工程，2002，31（6）：540-544.

ZHANG Guangjun， WU Xiaoli. Novel IR carbon dioxide sensor with high performance[J]. Infrared and Laser Engineering， 2002， 31（6）：540-544.

[2]涂志華，趙陽(yáng)，鄭力文，等.基于非分散紅外（NDIR）技術(shù)的土壤剖面二氧化碳濃度的測(cè)定[J].光譜學(xué)與光譜分析，2015，35（4）：997-1000.

TU Zhihua， ZHAO Yang， ZHENG Liwen， et al. Detection of carbon dioxide concentration in soil profile based on non-dispersive infrared spectroscopy technique[J]. Spectroscopy and Spectral Analysis，2015，35（4）：997-1000.

[3]閆玲，姚毅，昌學(xué)年，等. 紅外吸收法測(cè)定煙氣中CO2濃度的研究[J].儀器儀表與分析監(jiān)測(cè)，2009，100（4）：27-29.

YAN Ling， YAO Yi， CHANG Xuenian， et al. Study on infrared absorption measurement for the concentration of carbon dioxide in flue gas[J]. Instrumentation Analysis Monitoring， 2009， 100（4）：27-29.

[4]琚雪梅，張巍，畢東云，等. 紅外吸收型 CO2氣體傳感器的設(shè)計(jì)[J].傳感器技術(shù)，2005，24（8）：62-64.

JU Xuemei， ZHANG Wei， BI Dongyun， et al. Design of infrared absorption CO2 gas sensor[J]. Journal of Transducer Technology， 2005， 24（8）：62-64.

[5]王莉，陳小平.基于紅外傳感器的 CO2測(cè)量裝置[J].傳感器與微系統(tǒng)，2011，30（11）：115-117.

WANG Li， CHEN Xiaoping. CO2 measuring device based on infrared sensor[J]. Transducer and Microsystem Technologies， 2011，30（11）：115-117.

[6]崔穎，宋紹樓，彭繼慎，等.關(guān)于瓦斯傳感器的非線性校正的研究[J].微計(jì)算機(jī)信息，2010，26（25）：138-139.

CUI Ying， SONG Shaolou， PENG Jishen， et al. The study for non-linear correction based on gas sensor[J]. Microcomputer Information， 2010， 26（25）：138-139.

[7]王建梅，覃文忠. 基于L-M算法的BP神經(jīng)網(wǎng)絡(luò)分類器[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2005，30（10）：928-931.

WANG Jianmei， QIN Wenzhong. BP neural network classifier based on Levenberg-Marquardt algorithm[J]. Geomatics and Information Science of Wuhan University， 2005， 30（10）：928-931.

[8]趙映，陳小平.BP神經(jīng)網(wǎng)絡(luò)在紅外CO2體積分?jǐn)?shù)測(cè)量中的應(yīng)用[J].傳感器與微系統(tǒng)，2015，34（3）：151-153.

ZHAO Ying， CHEN Xiaoping. Application of BP neural networks in IR CO2 volume fraction measuring[J]. Transducer and Microsystem Technologies， 2015， 34（3）： 151-153.

[9]張鵬，楊興銳，嚴(yán)翔，等.基于BP神經(jīng)網(wǎng)絡(luò)的風(fēng)洞傳感器非線性誤差修正方法[J].兵工自動(dòng)化，2016，35（3）：31-34.

ZHANG Peng， YANG Xingrui， YANG Xiang， et al. Nonlinear error correction method of wind tunnel sensor based on BP neural network[J]. Ordnance Industry Automation， 2016， 35（3）：31-34.

[10]楊春生，牛紅濤，隋良紅，等.基于貝葉斯正則化算法BP神經(jīng)網(wǎng)絡(luò)釩電池SOC預(yù)測(cè)[J].現(xiàn)代電子技術(shù)，2016，39（8）：158-161.

YANG Chunsheng，NIU Hongtao，SUI Lianghong，et al. Application of BP neural network improved by Bayesian regularization algorithm in VRB SOC prediction[J]. Modern Electronics Technique， 2016， 39（8）：158-161.

[11]劉剛，劉學(xué)仁，嵇英華，等. 基于改進(jìn)型 BP 神經(jīng)網(wǎng)絡(luò)的瓦斯傳感器的非線性校正[J].傳感器與微系統(tǒng)，2007，26（1）：15-17.

LIU Gang，LIU Xueren， JI Yinghua，et al. Non-linear Correction of methane Sensor based on improved BP neural network[J]. Transducer and Microsystem Technologies， 2007， 26（1）：15-17.

[12]ROBERT F， THOMAS T. An automotive Bi-Source spectroscopic carbon dioxide sensor with pressure compensation[J]. Sensors and Actuators B， 2007， 127（1）： 82-88.

[13]KING T L， WEI M G， BREDA K， et al. Non-linear carbon dioxide determination using infrared gas sensors and neural network s with Bayesian regularization[J]. Sensors and Actuators B， 2009， 136（1）：242-247.

[14]汪獻(xiàn)忠，王曉東，赫樹(shù)開(kāi).基于紅外吸收型智能CO2濃度分析儀[J].儀表技術(shù)與傳感器，2010，325（2）：28-30.

WANG Xianzhong， WANG Xiaodong， HE Shukai. Design of intelligent CO2 concentration analyzer based on infrared absorption[J].Instrument Technique and Sensor， 2010， 325（2）： 28-30.

[15]陳明. Matlab神經(jīng)網(wǎng)絡(luò)原理與實(shí)例精解[M].北京： 清華大學(xué)出版社，2013.

[16]LAMPINEN J， VEHTARI A. Bayesian approach for neural networks：Review and case studies[J]. Neural Networks the Official Journal of the International Neural Network Society， 2001， 14（3）：257-274.

[17]MACKAY D J C. Probable networks and plausible prediction：A review of practical Bayesian methods for supervised neural networks[J]. Network Computation in Neural Systems， 1995， 6（3）：469-505.

[18]吳方良，石仲堃，楊向暉，等.基于L-M貝葉斯正則化方法的BP神經(jīng)網(wǎng)絡(luò)在潛艇聲納部位自噪聲預(yù)報(bào)中的應(yīng)用[J]. 船舶力學(xué)，2007，11（1）：136-142.

WU Fangliang， SHI Zhongkun，YANG Xianghui， et al.Submarine sonar self- noise forecast base on BP neural network and L-M Bayesian regularization algorithm[J]. Journal of Ship Mechanics， 2007，11（1）：136-142.

[19]OKUT H， WU X L， ROSA G J， et al. Predicting expected progeny difference for marbling score in Angus cattle using artificial neural networks and Bayesian regression models[J]. Genetics Selection Evolution， 2013， 45（1）：34.

[20]VITO S D， PIGA M， MARTINOTTO L， et al. CO， NO2， and NOx， urban pollution monitoring with on-field calibrated electronic nose by automatic bayesian regularization[J]. Sensors & Actuators B Chemical， 2009， 143（1）：182-191.