，，，

(1.蘇州工業(yè)園區(qū)職業(yè)技術(shù)學院，江蘇 蘇州 215123；2.江蘇大學 電氣信息工程學院，江蘇 鎮(zhèn)江 212013)

0 引言

海洋微生物作為微生物的一類，因其生存的海洋具有特殊的環(huán)境，故其所產(chǎn)生的酶與其他微生物所產(chǎn)生的酶相比具有更加獨特的性質(zhì)，如耐低溫，耐堿性，PH作用范圍寬等，這使得海洋微生物在食品加工、酶工業(yè)、添加劑和醫(yī)藥等發(fā)酵行業(yè)具有極大的開發(fā)潛力和應(yīng)用前景[1-5]。在海洋微生物發(fā)酵過程中，為保證發(fā)酵產(chǎn)物的品質(zhì)和質(zhì)量，需要實時檢測一系列生物參數(shù)，尤其是基質(zhì)濃度、菌體濃度及產(chǎn)物濃度(酶活)。當前，在線測量儀器僅能檢測發(fā)酵過程中某些物理和化學參數(shù)，還沒有成熟實用的儀器來測量這些關(guān)鍵生物參數(shù)[6-8]。在此背景下，許多關(guān)于微生物發(fā)酵過程中生物參數(shù)的軟測量方法應(yīng)運而生，其中，基于神經(jīng)網(wǎng)絡(luò)[9]的預(yù)測方法成為軟測量領(lǐng)域的研究熱點。然而，對于海洋微生物發(fā)酵過程這一類十分復(fù)雜的非線性系統(tǒng)來說，如果沒有先驗知識，就盲目應(yīng)用神經(jīng)網(wǎng)絡(luò)方法，關(guān)鍵生物參數(shù)的測量問題就不能很好的解決。為此，文中將核主元分析法(Kernel Principal Component Analysis，KPCA)[10-12]和動態(tài)模糊神經(jīng)網(wǎng)絡(luò)(Dynamic Fuzzy Neural Network，D-FNN)[13-15]相結(jié)合，提出了一種基于KPCA-DFNN的軟測量方法。

以典型的海洋微生物-海洋蛋白酶發(fā)酵過程為研究對象，首先，確定基質(zhì)濃度、菌體濃度、相對酶活(能夠更好地描述酶活的變化趨勢)這三個參量為軟測量模型的輸出變量，環(huán)境變量為模型的初始輸入變量。然后，利用KPCA對輸入變量進行數(shù)據(jù)壓縮和信息抽取，將所提取的主元作為DFNN的輸入，以上三個變量作為DFNN的輸出，建立了基于KPCA-DFNN的海洋蛋白酶發(fā)酵過程生物參量軟測量模型。仿真結(jié)果表明，該模型較基于DFNN和PCA- DFNN建模具有學習速度快、預(yù)測精度高等優(yōu)勢，有益于海洋蛋白酶的高效、高質(zhì)量生產(chǎn)。

1 核主元分析與動態(tài)模糊神經(jīng)網(wǎng)絡(luò)的原理

1.1 核主元分析

核主元分析是主元分析法(Principal Component Analysis，PCA)的非線性推廣，其具體算法如下：

給定n個樣本，樣本集X={x1,x2,…,xn}，xk∈Rm，由非線性函數(shù)φ(·)將輸入數(shù)據(jù)從原空間映射到高維特征空間F，F(xiàn)中的樣本記為φ(xk)，且滿足：

(1)

在F空間中樣本的協(xié)方差矩陣C為

(2)

對C進行特征值分解

λV=CV

(3)

式中，V是與λ對應(yīng)的特征向量，特征值λ≥0。將(3)式的兩邊左乘以核樣本φ(xk)：

λφ(xk)·V=φ(xk)·CV，k=1,2,…,n

(4)

解方程可得與非零特征值對應(yīng)的特征向量V。其解一定處于φ(x1),φ(x2),···,φ(xn)張成的空間中，所以V可以表示為：

(5)

其中，

φ(x)=[φ(x1),φ(x2),…,φ(xn)]，α表示a1,…,an中的一個列向量。

引入核函數(shù)：

Kij=K(xi,xj)=[φ(xi),φ(xj)],i,j=1,2,…,n

(6)

nλα=Kα

(7)

歸一化特征向量V，即(V,V)=1，即可得樣本x在特征空間中的第k(k=1,2,…,n)個主元分量tk(x)：

(8)

(9)

式(9)表示前p個λk的和與總和比值大于E，通常選取E>85%。

(10)

1.2 動態(tài)模糊神經(jīng)網(wǎng)絡(luò)

第1層：輸入層，每個節(jié)點代表一個輸入變量。

第2層：隸屬函數(shù)層，每個節(jié)點代表一個隸屬函數(shù)，該隸屬函數(shù)可用式(11)表示：

(11)

其中，i=1,2,…,r,j=1,2,…,u，r是輸入變量數(shù)，u是隸屬函數(shù)的數(shù)量，即系統(tǒng)的總規(guī)則數(shù)，cij和σj分別是xi的第j個高斯隸屬函數(shù)的中心和寬度，xmin是xi的第j個高斯隸屬函數(shù)。

第3層：T-范數(shù)層，每個節(jié)點代表一個可能的模糊規(guī)則中的IF-部分，即該層的節(jié)點數(shù)反映了模糊規(guī)則數(shù)。第j個規(guī)則Rj的輸出為：

(12)

其中：j=1,2,…,u，X=(x1,x2,...,xr)∈Rr，Cj=(c1j,...,crj)∈Rr是第j個RBF單元的中心。

第4層：歸一化層，這些節(jié)點被稱為N節(jié)點。其數(shù)目等于模糊規(guī)則的節(jié)點數(shù)。第j個節(jié)點Nj的輸出為：

(13)

第5層：輸出層，每個節(jié)點代表一個輸出量，此輸出是所有輸入信號的疊加：

(14)

式中，wk是THEN-部分或者稱為第k個規(guī)則的連接權(quán)，wk=ak0+ak1x1+…akrxr,k=1,2,…u，y是變量的輸出。

圖1 DFNN的結(jié)構(gòu)圖

2 KPCA-DFNN模型的構(gòu)建與驗證

2.1 模型的構(gòu)建

基質(zhì)濃度S、菌體濃度X、相對酶活P對海洋蛋白酶發(fā)酵過程的優(yōu)化控制非常重要，因此，選擇這三個變量作為KPCA-DFNN軟測量模型的輸出變量。通過分析海洋蛋白酶發(fā)酵機理并結(jié)合發(fā)酵過程的實驗數(shù)據(jù)，選取了10個影響生物參數(shù)S、X、P的主要因素作為初始樣本變量，分別是：時間t、溫度T、攪拌速度r、溶解氧濃度DO、空氣流量l、pH值、CO2濃度、基質(zhì)進給速率u、發(fā)酵罐壓力p、反應(yīng)器體積v。

采集了10批發(fā)酵數(shù)據(jù)，前9批作為模型的訓練樣本集，第10批作為模型的測試樣本集，由于采集到的樣本數(shù)據(jù)變化范圍較大，如果直接使用原始測量數(shù)據(jù)進行計算，不僅會夸大量綱數(shù)據(jù)的作用，而且還可能導(dǎo)致信息丟失或引起數(shù)值計算的不穩(wěn)定。需要對樣本數(shù)據(jù)進行歸一化處理。公式如下：

(15)

式中：xmax為樣本數(shù)據(jù)的最大值，xmin為樣本數(shù)據(jù)的最小值，x為原始樣本數(shù)據(jù)，x′為歸一化后的數(shù)據(jù)。歸一化后樣本數(shù)據(jù)在[0，1]之間。

根據(jù)核主元分析法和動態(tài)模糊神經(jīng)網(wǎng)絡(luò)的基本原理，構(gòu)建了基于KPCA-DFNN的海洋蛋白酶發(fā)酵過程生物參量的軟測量模型，建模過程如圖2所示。建模步驟如下。

步驟1：根據(jù)建模對象選取適當?shù)妮斎胼敵鰳颖緮?shù)據(jù)。

步驟2：利用式(15)對輸入輸出數(shù)據(jù)進行預(yù)處理。

步驟3：根據(jù)KPCA算法對輸入變量進行數(shù)據(jù)壓縮和信息抽取，消除輸入變量之間的相關(guān)性，進行特征選取。文中按累積方差百分比大于95% ，KPCA選定了2 個特征主元，PCA選定了6個特征主元。

步驟4：將提取非線性主元作為DFNN的輸入，X、S、P作為模型的輸出變量，利用訓練樣本集對DFNN模型進行訓練，選取最佳模型構(gòu)建參數(shù)。

步驟5：利用測試樣本集對建好的KPCA-DFNN軟測量模型進行驗證。

圖2 海洋蛋白酶發(fā)酵過程預(yù)測模型框圖

2.2 模型驗證

用測試樣本集對建好的KPCA-DFNN模型進行仿真驗證。仿真結(jié)果如圖3、圖4和表2所示。

圖3顯示出了海洋蛋白酶發(fā)酵過程X、S、P的離線化驗值(真實值)和軟測量值(預(yù)測值)對比結(jié)果。

圖3 KPCA-DFNN生物參量預(yù)測曲線

圖4 DFNN、PCA-DFNN、KPCA-DFNN生物參量預(yù)測值誤差曲線

雖然試驗過程中采集到的樣本值分散性很大和重復(fù)性很小，但從圖3中可以看出，對于X、S、P，基于KPCA-DFNN的軟測量模型，輸出的軟測量值都能夠很好的追蹤離線化驗值，這說明，KPCA-DFNN具有較好的逼近能力。

采用均方根誤差(RMSE)和平均絕對誤差(MAD)更直觀的反映DFNN、PCA-DFNN、KPCA-DFNN建模方式的預(yù)測效果，如表1所示。

表1 DFNN、PCA-DFNN、KPCA-DFNN預(yù)測值誤差對比

(16)

(17)

其中:n是數(shù)據(jù)對的數(shù)目，t(i)和y(i)是第i個期望輸出與實際輸出。

另外，為了更清晰地描述DFNN、PCA-DFNN、KPCA-DFNN的預(yù)測能力，圖4中給出了這三種建模方法預(yù)測誤差。

從表1和圖4中可以看出，基于KPCA-DFNN的建模所得到基質(zhì)濃度、菌體濃度和相對酶活的預(yù)測誤差要比DFNN和PCA-DFNN小得多，這進一步表明，即使在同樣的訓練數(shù)據(jù)和測試數(shù)據(jù)的條件下，KPCA-DFNN的預(yù)測效果比DFNN和PCA-DFNN好。既基于KPCA-DFNN在發(fā)酵樣本預(yù)測中具有顯著的作用。

3 結(jié)論

為解決海洋微生物發(fā)酵過程中生物參量難以實時在線測量的問題，文中以海洋蛋白酶發(fā)酵過程為例，將核主元分析法(KPCA)與動態(tài)模糊神經(jīng)網(wǎng)絡(luò)(DFNN)相結(jié)合，提出一種基于KPCA-DFNN的建模方法。首先確定了X、S、P作為海洋蛋白酶預(yù)測模型的輸出變量，同時，通過對海洋蛋白酶發(fā)酵過程進行機理分析，選定了10個影響輸出變量的環(huán)境變量作為KPCA-DFNN模型的輸入變量。用訓練樣本集對這些變量進行預(yù)處理后，按照累積方差百分比大于95%，KPCA提取了2 個特征主元。根據(jù)DFNN算法對提取后的數(shù)據(jù)進行學習訓練，建立了基于KPCA-DFNN的海洋蛋白酶發(fā)酵過程生物參量的軟測量模型，用測試樣本集對模型進行仿真驗證，試驗結(jié)果表明與DFNN、PCA-DFNN建模方法相比，所建立的KPCA-DFNN模型具有良好的預(yù)測精度，所得生物參量預(yù)測值的最大RMSE為0.582 9，最大MAD為0.454 2，滿足發(fā)酵過程中生物參量的在線測量要求。