周安然，楊光永，2，吳海鋒，趙錦劍，朱 芹

(1．云南民族大學(xué)電氣信息工程學(xué)院，云南昆明 650500；2．華南理工大學(xué)機(jī)械與汽車(chē)工程學(xué)院，廣東廣州 510640)

0 引言

近20年來(lái)，離子濃度測(cè)量的研究熱點(diǎn)主要是圍繞電極過(guò)程動(dòng)力學(xué)[1]、提高離子電極性能[2]和離子濃度信號(hào)通道特性[3]等方面進(jìn)行，主要分為兩大類(lèi)：一類(lèi)著眼于改進(jìn)離子濃度儀的硬件設(shè)計(jì)、電氣特性[4]，構(gòu)建測(cè)量平臺(tái)的最優(yōu)系統(tǒng)方案，如選擇性、免校準(zhǔn)離子電極[5]，采用具有超高輸入阻抗、高共模抑制比和極低輸入偏置電流的儀表放大器[6]，采用可配置模擬前端(AFE)專(zhuān)用器件(如LMP91200)，選擇16位甚至24位Δ-Σ型AD轉(zhuǎn)換器以提高信號(hào)的轉(zhuǎn)換精度，或采用帶AD轉(zhuǎn)換器的增強(qiáng)型混合單片機(jī)[7]；另一類(lèi)側(cè)重于噪聲抑制和非線性標(biāo)定技術(shù)，采用現(xiàn)代信號(hào)方法，包括遞歸最小均方誤差濾波器(RLS)[8]、最大熵譜方法[9]、盲源分離[10]等，抑制超高斯和亞高斯噪聲，以低階統(tǒng)計(jì)量信號(hào)的穩(wěn)態(tài)過(guò)程看待離子活度測(cè)量系統(tǒng)。

無(wú)論是在線水質(zhì)監(jiān)測(cè)，還是離線水質(zhì)樣品分析，由于選擇性離子電極的暫態(tài)響應(yīng)特性，使離子活度信號(hào)不可避免地混入高斯噪聲和非高斯噪聲[11]，而基于低階統(tǒng)計(jì)量的RLS、LMS濾波器很難抑制或消除高階統(tǒng)計(jì)噪聲的影響，從而降低離子濃度儀的測(cè)量精度和信噪比。

為了提高離子濃度儀的測(cè)量性能，文中利用高階統(tǒng)計(jì)量均方誤差準(zhǔn)則對(duì)非高斯噪聲的強(qiáng)一致估計(jì)特性，在測(cè)量系統(tǒng)的模擬前端與微處理器單元之間加入高階統(tǒng)計(jì)量自適應(yīng)濾波器(HoSRLS)，消除離子濃度信號(hào)中的超高斯和亞高斯噪聲。

1 高階統(tǒng)計(jì)量自適應(yīng)濾波器

設(shè)離子濃度測(cè)量系統(tǒng)的輸出過(guò)程為：

(1)

式中：誤差e(n-k)是均值為零、與沖激響應(yīng)h(n)相互獨(dú)立、同分布的隨機(jī)信號(hào)。

x(n)經(jīng)過(guò)濾波器的輸出εx，θ(n) 為：

(2)

在四階統(tǒng)計(jì)量的均方誤差代價(jià)函數(shù)：

(3)

式中：

B=E[y(n+k1)εy，w(n)]E[y(n+k2)εy，θ(n)]；

(4)

(5)

式中：γ2c和γ4e分別是e(n)的二階和四階累積量；β為標(biāo)量因子。

(6)

按高階統(tǒng)計(jì)量均方誤差準(zhǔn)則建立新的強(qiáng)一致估計(jì)器：

Φw=φ

(7)

滿足遞歸形式：

(8)

窗口寬度為2M+1，每次有2M+1個(gè)數(shù)據(jù)被分析。式(7)中w是濾波器在n時(shí)刻的抽頭權(quán)向量：

w(n)=[w0(n)，w1(n)，…，wM-1(n)]T

(9)

使y(n)-y(n+1)達(dá)到最大值：

max[y(n)-y(n+1)]

(10)

即式(10)成立時(shí)所取得的n即為AR(p)的階數(shù)。

2 仿真和模型驗(yàn)證

為了驗(yàn)證高階統(tǒng)計(jì)量自適應(yīng)濾波器的收斂速度和噪聲抑制性能，文中以不同噪聲功率的加性噪聲、不同幅值的窗函數(shù)設(shè)計(jì)仿真數(shù)學(xué)模型，以獲得仿真數(shù)據(jù)。

加性噪聲由三類(lèi)噪聲組成：

(1)高斯白噪聲。X：N(μ，σ2)，其中，均值μ=0，方差σ2=0.1，噪聲功率Ρ=5 dBW，用于模擬正態(tài)分布的加性噪聲分量；

(2)瑞利分布噪聲。X：Reyleigh(B)，分布參數(shù)B取值25，用于模擬亞高斯噪聲分量；

(3)指數(shù)分布噪聲。X：Exponential(λ)，率參數(shù)λ=0．1，模擬超高斯噪聲分量。

設(shè)采樣序列為k∈[1，200]，以氫離子電極動(dòng)力特性為仿真對(duì)象，將參考信號(hào)d劃分為不同幅值的Blackman-Harris窗函數(shù)d1、d2和d3：

系統(tǒng)仿真數(shù)學(xué)模型為：

x=d+N(P)+Reyleigh(В)+Exponential(λ)

在采樣序列k1=160、k2=100、k3=40，分別獲得3個(gè)極值點(diǎn)，以測(cè)試在極值和不同峭度條件下HOSRLS濾波器的收斂性能。

HOSRLS和RLS濾波器的仿真結(jié)果如圖1所示，從圖中可看出，當(dāng)?shù)螖?shù)較少(i=5)，在峭度和偏度較大的區(qū)域(如A、B和C)，RLS濾波器仍然含有超高斯和亞高斯噪聲分量Reyleigh(B)和Exponential(λ)，且超調(diào)量高于10%，而HOSRLS則完全濾除超高斯和亞高斯噪聲，其收斂速度優(yōu)于RLS.

圖1 HOSRLS和RLS濾波器的仿真結(jié)果

另一方面，在同樣的迭代次數(shù)(i=5)條件下(如圖2所示的D、E和F點(diǎn))，RLS和HOSRLS抑制高斯白噪聲的能力相當(dāng)，但對(duì)于疊加非高斯噪聲分量的系統(tǒng)仿真數(shù)學(xué)模型來(lái)說(shuō)，HOSRLS濾波器在采樣序列k1=160、k2=100、k3=40分別獲得45 dB、28 dB和55 dB，而RLS分別獲得25 dB、25 dB和22 dB，可看出HOSRLS比RLS濾波器獲得較高的信噪比。

圖2 HOSRLS和RLS濾波器的信噪比比較

3 實(shí)驗(yàn)及結(jié)果分析

實(shí)驗(yàn)平臺(tái)的總體結(jié)構(gòu)框圖如圖3所示，在Cyclone III系列EP3C10F256上創(chuàng)建系統(tǒng)，以標(biāo)準(zhǔn)型NIOS II微處理器為核心，將HOSRLS濾波器模塊掛接到Avalon總線上，以硬件加速方式設(shè)計(jì)HOSRLS模塊，NIOS II微處理器訪問(wèn)HOSRLS模塊，并實(shí)時(shí)計(jì)算離子濃度，經(jīng)IEEE 802．15．4收發(fā)片上系統(tǒng)CC2430，將離子濃度值發(fā)送到遠(yuǎn)程數(shù)據(jù)管理中心。

圖3 無(wú)線傳感器網(wǎng)絡(luò)離子濃度監(jiān)測(cè)節(jié)點(diǎn)結(jié)構(gòu)框圖

為了提高信號(hào)整定單元的性能，減少溫度漂移，取消高輸入阻抗運(yùn)算放大器ICL7621和電壓溫度漂移補(bǔ)償電路，精簡(jiǎn)傳感器及其模擬前端的電路結(jié)構(gòu)，采用化學(xué)傳感器可配置模擬前端LMP91200，NIOS II處理器經(jīng)SPI總線，配置離子電極的共模電壓、對(duì)溫度補(bǔ)償傳感器的電流源編程配置、選擇標(biāo)定/測(cè)試模式等。

無(wú)線傳感器網(wǎng)絡(luò)離子濃度監(jiān)測(cè)節(jié)點(diǎn)在1．68、4．01、6．86、7．00和10．00等不同pH標(biāo)準(zhǔn)緩沖溶液、25 ℃標(biāo)定實(shí)驗(yàn)條件下，測(cè)量并計(jì)算其絕對(duì)誤差和獨(dú)立線性度，其結(jié)果見(jiàn)表1?？煽闯?，加入HOSRLS濾波器模塊之后，離子濃度監(jiān)測(cè)節(jié)點(diǎn)的線性度有明顯提高。

表1 絕對(duì)誤差和線性度測(cè)試結(jié)果

4 結(jié)論

在文中采用基于四階累積量的均方誤差準(zhǔn)則來(lái)設(shè)計(jì)自適應(yīng)濾波器，用于離子敏傳感器來(lái)測(cè)量溶液中氫離子濃度。仿真和實(shí)驗(yàn)結(jié)果表明高階統(tǒng)計(jì)量自適應(yīng)濾波器方法比低階自適應(yīng)濾波器方法在氫離子濃度的測(cè)量上，極大地提高了信噪比，濾波效果明顯優(yōu)于低階自適應(yīng)濾波器，收斂速度更快，這也驗(yàn)證了高階統(tǒng)計(jì)量能有效地抑制高斯噪聲和非高斯噪聲。

參考文獻(xiàn)：

[1] 查全性．電極過(guò)程動(dòng)力學(xué)導(dǎo)論．3版 北京：科學(xué)出版社，2002．

[2] 李芬．鋅空氣電池之氣體擴(kuò)散電極性能研究：[學(xué)位論文].合肥：中國(guó)科學(xué)技術(shù)大學(xué)，2010．

[3] 黃運(yùn)生，毛先柏，陳學(xué)．基于微控制器的智能鈷離子濃度分析儀．儀表技術(shù)與傳感器，2005(10)：6-8．

[4] 楊光永．KA0324A 離子濃度計(jì)的研究與設(shè)計(jì)：[學(xué)位論文].昆明：昆明理工大學(xué)，2005．

[5] 馬鈺華．金屬離子選擇性電極的研究及應(yīng)用：[學(xué)位論文].重慶：西南大學(xué)，2010．

[6] 黃勤，周婷，劉益良，等．高共模抑制比腦電放大電路設(shè)計(jì)．儀器儀表與檢測(cè)技術(shù)，2009，28(4)：77-79．

[7] 盧立蘋(píng)．離子濃度計(jì)的研制：[學(xué)位論文].沈陽(yáng)：東北大學(xué)，2008．

[8] 江清潘，常太華，朱紅路，等．自適應(yīng)噪聲消除算法的性能比較與仿真．科學(xué)技術(shù)與工程，2009，9(19)：5835- 5839．

[9] 聶士忠，王玉泰．最大熵譜分析方法和MATLAB中對(duì)短記錄資料的譜分析．山東師范大學(xué)學(xué)報(bào)(自然科學(xué)版)，2005，20(3)：40-41．

[10] CARDOSO J F．Blind signal separation：statistical principles．Proceedings of the IEEE，1998，86(10)：2009-2025．

[11] YANG G Y,HU G Q,GUO B．Multiple Reflection Cancellation Using High Order Statistics Adaptive Filter for Scattered Triangulation Laser Displacement Sensor．Advanced Materials Research，2013，774-776：1613-1616．

[12] DELOPOULOS A N,GIANNAKIS G B．Strongly Consistent Identification Algorithms and Noise Insensitive MSE Criteria．IEEE Transactions on Signal Processing，1992，40(8)：1955-1970．

[13] 田昭武．電化學(xué)研究方法．北京：科學(xué)出版社，1984．