陳 希，李麗娜

（華僑大學(xué)機電及自動化學(xué)院，福建 廈門361021）

隨著社會經(jīng)濟的高速發(fā)展、工業(yè)化進程的不斷加快，全球水資源問題日益嚴重。中國水資源占有總量排在世界前列，但由于人口眾多，中國人均水資源占有量僅有世界人均占有量的1/4。保護水資源、防止水污染至關(guān)重要[1]。

水質(zhì)參數(shù)（包括物理性指標(biāo)、化學(xué)性指標(biāo)和生物性指標(biāo)）檢測是保障人類安全用水的前提，通過水質(zhì)檢測可確定水污染分布狀況、預(yù)測污染物的變化趨勢，為水污染防治工作提供技術(shù)支持[2]。

水質(zhì)參數(shù)檢測方法包括常規(guī)化學(xué)分析檢測方法、紫外/可見光譜法、熒光光譜法、原子吸收光譜法、氣相色譜法和可見/近紅外（visible near infrared，VIS-NIR）光譜法等[3-4]?？梢?近紅外光譜分析法具有檢測速度快、無污染、綠色環(huán)保、操作簡單等優(yōu)點，是近年來水質(zhì)參數(shù)檢測的研究熱點之一。

基于可見/近紅外光譜技術(shù)的水質(zhì)參數(shù)檢測方法，在最終實現(xiàn)應(yīng)用前，通常需通過標(biāo)定實驗來確定測量模型[5-7]。測量模型的訓(xùn)練通常在MATLAB軟件上進行，通過反復(fù)多次計算和分析，從而選擇適當(dāng)?shù)念A(yù)處理方法、確定最佳的多元回歸模型參數(shù)，模型訓(xùn)練的計算和分析過程較為復(fù)雜耗時。

為方便研究人員開展水質(zhì)光譜分析實驗數(shù)據(jù)的分析，提高測量模型的建模效率，本設(shè)計研究開發(fā)了一套基于MATLAB的水質(zhì)參數(shù)光譜分析的GUI界面，通過導(dǎo)入數(shù)據(jù)、預(yù)處理、多元校正、數(shù)據(jù)輸出4個模塊，提供常見的預(yù)處理方法和經(jīng)典的多元回歸方法。并通過蒸餾水pH值光譜分析基礎(chǔ)實驗數(shù)據(jù)，對該可視化軟件效果進行了驗證。

1 水質(zhì)參數(shù)光譜分析界面框架設(shè)計

圖形用戶界面，是用戶與計算機程序之間建立聯(lián)系，并進行數(shù)據(jù)交流的一種方式。它包含窗口、圖標(biāo)、菜單、文本、工具欄等圖形對象。用戶借助操作或更改的方式，選擇或激活圖形對象，建立屬于各自研究問題領(lǐng)域的可視化交互界面。通過圖形用戶界面，用戶只要知道如何使用界面的每個組件，就可以與計算機建立聯(lián)系，進行數(shù)據(jù)交流，操作簡單、直觀[8]。

本設(shè)計主要面向水質(zhì)參數(shù)光譜分析領(lǐng)域，基于MATLAB/GUI編輯調(diào)整的可視化分析界面如圖1所示。

圖1 水質(zhì)參數(shù)光譜分析的可視化界面

主要包括導(dǎo)入數(shù)據(jù)模塊、預(yù)處理模塊、多元校正模塊、數(shù)據(jù)輸出模塊。導(dǎo)入數(shù)據(jù)模塊的功能是，導(dǎo)入校正集、驗證集數(shù)據(jù)，并可顯示原始數(shù)據(jù)的光譜圖；預(yù)處理模塊的功能是，提供常見的光譜預(yù)處理方法，包括最大最小歸一化、零均值化、一階導(dǎo)數(shù)、二階導(dǎo)數(shù)、多元散射校正（multiple scatter correction，MSC）和標(biāo)準(zhǔn)正態(tài)變量變換（standard normal variate transformation，SNV），并可顯示預(yù)處理后的光譜圖；多元校正模塊的功能是，提供經(jīng)典的偏最小二乘（partial least squares，PLS）回歸模型，通過輸入主成分個數(shù)，計算該參數(shù)條件下的預(yù)測均方根誤差（root mean square error of prediction，RMSEP），并可顯示預(yù)測值與參考值的比較圖，為工程應(yīng)用人員進行模型訓(xùn)練提供參考。

2 基于MATLAB/GUI的界面開發(fā)方法

MATLAB是美國MathWorks公司開發(fā)的一套高性能數(shù)字計算軟件，MATLAB既能進行科學(xué)計算，又能根據(jù)用戶需求，開發(fā)出便于計算和分析的用戶圖形界面（GUI）。用戶通過選擇、激活這些圖形對象，可實現(xiàn)計算、繪圖，提高工程應(yīng)用人員研究效率[9-10]。本研究設(shè)計的水質(zhì)參數(shù)光譜分析可視化界面，主要開發(fā)步驟和方法如下。

首先，在MATLAB主界面的Command Window中輸入guide，得到用戶圖形界面模板選擇窗口，選擇默認空白（BLANK GUI）選項，如圖2所示，然后可進入編輯界面，如圖3所示。

圖2 GUI模板選擇窗口

圖3 GUI編輯界面

接著，在GUI編輯界面上建立需要的控件，編輯控件的基本屬性，將界面各個控件擺放整齊、美觀、方便。

然后，根據(jù)各控件需要完成的功能，編輯回調(diào)函數(shù)。

最后，進行可視化界面的結(jié)構(gòu)優(yōu)化設(shè)計和程序調(diào)試。

對于導(dǎo)入數(shù)據(jù)模塊，由八個靜態(tài)文本框、六個按鈕組成。設(shè)置靜態(tài)文本框的屬性作為提醒使用者的提示，導(dǎo)入數(shù)據(jù)按鈕按下后運行對應(yīng)的回調(diào)函數(shù)，打開文件夾并且從文件夾中選擇需要的文件加載到程序中。

MATLAB函數(shù)global可用于定義全局變量。

MATLAB函數(shù)[filename,filepath]=uigetfile('*.*','請選擇文件')可實現(xiàn)從文件夾中選擇數(shù)據(jù)。

MATLAB函數(shù)set(handles.text6,'String',[filename])可實現(xiàn)將文件名顯示在界面上。

MATLAB函數(shù)importdata(strcat([filepath,filename]))可實現(xiàn)文件的導(dǎo)入。

MATLAB函數(shù)axes(handles.axes1)、plot，可實現(xiàn)在axes1坐標(biāo)上繪圖，點擊“生成光譜圖”按鈕，就會在界面繪圖區(qū)顯示光譜圖。

類似的，對于預(yù)處理模塊、多元校正模塊、數(shù)據(jù)輸出模塊，基于MATLAB編輯好對應(yīng)的預(yù)處理函數(shù)、多元回歸函數(shù)及其回調(diào)函數(shù)，可實現(xiàn)預(yù)處理方法的選擇、PLS模型參數(shù)“主成分個數(shù)”設(shè)置、PLS模型的運行、RMSEP的計算以及預(yù)測值與參考值比較的圖形繪制。

MATLAB函數(shù)k=str2num(get(handles.edit1,'String'))可實現(xiàn)從可輸入編輯框中獲得參數(shù)k。

MATLAB函數(shù)guidata(hObject,handles)可實現(xiàn)從可輸入編輯框中獲得輸入的主成分個數(shù)。

3 蒸餾水pH值光譜分析實驗數(shù)據(jù)處理

采用可見/近紅外分光光度計V-1800PC（上海美譜達儀器公司）測量樣品光譜，光譜波長間隔2 nm。1 cm的石英池用于裝載水樣。pH值的參考值采用精度為±0.1個單位的酸度計（上海力辰儀器公司）測量。為了獲得不同酸堿度的水樣，在蒸餾水中隨機加入0.1 mmol/L的HCL或0.1 mmol/L的NaOH，用酸度計測量獲得pH參考值。采用透射式測量獲取每個樣品的光譜數(shù)據(jù)。采集的光譜波長范圍是400～1 100 nm，共351個波長。

實驗中共獲得34個樣本，其中24個樣本作為校正集、11個樣本作為驗證集，用于PLS模型訓(xùn)練，來測試本設(shè)計開發(fā)GUI界面的性能。圖4為應(yīng)用本設(shè)計開發(fā)的可視化界面繪制的圖像，圖4（a）為原始光譜圖，圖4（b）為經(jīng)標(biāo)準(zhǔn)正態(tài)變量變換預(yù)處理后的光譜圖。通過設(shè)置不同的主成分個數(shù)，然后運行PLS，可以得到不同的模型精度參數(shù)值RMSEP，根據(jù)RMSEP來確定最優(yōu)的主成分個數(shù)，應(yīng)用本研究設(shè)計開發(fā)的可視化界面，可以方便的調(diào)整參數(shù)，確定最佳的PLS模型，提高建模效率。例如，經(jīng)標(biāo)準(zhǔn)正態(tài)變量變換預(yù)處理后，設(shè)置主成分個數(shù)為10，PLS模型的預(yù)測參數(shù)RMSEP顯示為0.335 221。PLS模型的預(yù)測值與參考值的對比圖如圖5所示。

圖4 GUI輸出光譜圖界面

圖5 GUI輸出PLS模型的預(yù)測值與參考值的對比圖頁面

4 結(jié)語

水質(zhì)檢測用于確保用水安全，對居民生活、工業(yè)生產(chǎn)和環(huán)境保護都非常重要。可見/近紅外光譜分析技術(shù)具有簡便、快速、經(jīng)濟、無污染等優(yōu)點，是目前極具發(fā)展前景的水質(zhì)參數(shù)檢測方法之一。為方便研究人員進行光譜分析，本研究設(shè)計開發(fā)了一套基于MATLAB/GUI的水質(zhì)參數(shù)光譜分析的可視化界面。應(yīng)用該水質(zhì)參數(shù)光譜分析界面，可快速切換預(yù)處理方法和主成分個數(shù)，根據(jù)PLS模型精度參數(shù)、以及預(yù)測值和參考值的圖像分析，達到快速確定最佳預(yù)處理方法和多元校正模型，提高建模效率的目的。基于MATLAB/GUI的可視化圖形界面，設(shè)計界面簡潔清晰，可極大方便研究人員開展水質(zhì)參數(shù)光譜分析工作，提高研究效率。