基于改進(jìn)免疫遺傳算法的近紅外光譜變量選擇方法

陶煥明，高美鳳

（江南大學(xué) 輕工過(guò)程先進(jìn)控制教育部重點(diǎn)實(shí)驗(yàn)室，物聯(lián)網(wǎng)工程學(xué)院，江蘇 無(wú)錫 214122）

近紅外光譜分析檢測(cè)技術(shù)具有快速、無(wú)損以及低成本的優(yōu)勢(shì)，已經(jīng)廣泛運(yùn)用于眾多領(lǐng)域。如何快速有效地提取出有用光譜信息，即光譜變量選擇，是近紅外光譜技術(shù)研究的重點(diǎn)之一。目前，許多國(guó)內(nèi)外學(xué)者提出了大量基于不同原理、策略的變量選擇方法，如通過(guò)無(wú)信息變量消除法（UVE）對(duì)偏最小二乘（PLS）回歸系數(shù)設(shè)定閾值限制來(lái)選擇有效變量［1］；區(qū)間偏最小二乘（iPLS）［2］、移動(dòng)窗口偏最小二乘（MWPLS）［3］、向前和向后間隔偏最小二乘（FB－iPLS）［4］、區(qū)間隨機(jī)蛙跳（iRF）［5］等則均為以光譜區(qū)間為篩選對(duì)象的方法。近年來(lái)，以智能算法及其改進(jìn)算法為核心進(jìn)行波長(zhǎng)組合優(yōu)化的算法研究較多，主要有遺傳算法（GA）［6－8］、粒子群算法（PSO）［9］、模擬退火算法（SAA）［10］、蟻群算法（ACO）［11］等。GA算法具有群體搜索特性，廣泛應(yīng)用于多種領(lǐng)域，但易過(guò)早收斂。劉鑫等［12］將改進(jìn)的GA 算法用于土壤養(yǎng)分預(yù)測(cè)，采用一種改進(jìn)的實(shí)數(shù)編碼差分變異算子，擴(kuò)大了全局最優(yōu)解搜索空間，加快了收斂速度，提高了預(yù)測(cè)精度。免疫遺傳算法（IGA）在GA 算法的基礎(chǔ)上引入了生物免疫機(jī)制，本質(zhì)也是對(duì)GA 算法的一種改進(jìn)［13－15］。相較于GA算法及其改進(jìn)算法，IGA算法避免了容易過(guò)早收斂的缺點(diǎn)，維持了種群中個(gè)體的多樣性，可有效防止算法陷入局部最優(yōu)，但其精度性能仍存在提升空間。國(guó)內(nèi)外對(duì)IGA 算法的研究也從未停止腳步［16－17］。

本文提出了一種改進(jìn)免疫遺傳算法（iIGA）用于近紅外光譜變量選擇。IGA 算法結(jié)合了適應(yīng)度（親和度）和抗體濃度形成聚合適應(yīng)度進(jìn)行后續(xù)的遺傳操作來(lái)維持種群的多樣性以及免疫平衡，抗體濃度越高，越受抑制，反之亦然。其中抗體濃度是指與抗體相似的個(gè)體總數(shù)和群體總數(shù)的比值，其中抗體之間相似度達(dá)到相應(yīng)的閾值判定抗體相似。原算法中固定了抗體相似度閾值，可能會(huì)存在閾值過(guò)大或者過(guò)小而導(dǎo)致抗體濃度不具備差異性，使得引入的抗體濃度失去意義。本文提出自適應(yīng)相似度閾值，所選取的閾值可使群體抗體濃度的標(biāo)準(zhǔn)差達(dá)到最大，完全體現(xiàn)出抗體濃度的差異性；同時(shí)采用精英保留策略，將親和度高的若干個(gè)抗體存入到記憶細(xì)胞中，直接復(fù)制到下一代；另外在算法中引入貪心算法思想，將算法由局部最優(yōu)擴(kuò)展到全局最優(yōu)。將每一次迭代中的親和度最優(yōu)基因隨機(jī)反轉(zhuǎn)一位基因位進(jìn)行局部性探優(yōu)，并且保證探優(yōu)之后的基因親和度值大于探優(yōu)之前；然后以探優(yōu)成功后的基因替換本次迭代中親和度最差基因，使得算法朝著正確的方向進(jìn)行局部性探優(yōu)，最終達(dá)到全局最優(yōu)，提升模型預(yù)測(cè)精度。

1 算法原理

1.1 IGA算法實(shí)現(xiàn)原理

GA算法中通過(guò)模擬種群個(gè)體對(duì)環(huán)境適應(yīng)度大小的不同進(jìn)行優(yōu)勝劣汰，經(jīng)過(guò)不斷的迭代更新，使種群個(gè)體對(duì)環(huán)境的適應(yīng)度不斷提升，最終達(dá)到群體適應(yīng)度的最佳表現(xiàn)。但在算法迭代后期，種群個(gè)體之間的差異性越來(lái)越小，導(dǎo)致算法陷入局部最優(yōu)。因此，免疫遺傳算法在個(gè)體對(duì)環(huán)境的適應(yīng)度中引入抗體濃度，促使抗體濃度越高，適應(yīng)度反遭抑制，維持了種群中個(gè)體的多樣性，可有效防止算法陷入局部最優(yōu)。

1.1.1 產(chǎn)生初始解 在種群中設(shè)置N個(gè)個(gè)體，作為算法模型的抗體角色。將需篩選的近紅外光譜波段均分為L(zhǎng)個(gè)波段，每個(gè)波段包括多個(gè)波長(zhǎng)點(diǎn)，對(duì)應(yīng)抗體基因上的基因位，按照0－1 形式對(duì)其進(jìn)行隨機(jī)編碼，1表示選中該波段，0表示未選中。

1.1.2 計(jì)算親和度 對(duì)抗體基因進(jìn)行編碼后，將選中的波段與對(duì)應(yīng)的理化值進(jìn)行PLS建模，親和度即以建模得到的校正集均方根誤差（RMSEC）和相關(guān)系數(shù)（R）為變量的函數(shù)值，本文設(shè)定親和度Z為：

1.1.3 計(jì)算抗體濃度 抗體濃度用與抗體i相似的抗體個(gè)數(shù)ci與群體總個(gè)數(shù)N的比值di來(lái)表示，即：

抗體之間相似度達(dá)到一定閾值，即判定抗體相似。相似度是抗體與其他抗體之間的相似程度，其度量方法可基于信息熵、歐氏距離、基因值距離以及海明距離。其中海明距離主要用于計(jì)算離散編碼，符合本文的0－1編碼形式。因此本文采用海明距離作為抗體相似度的度量方法，具體公式如下：

式（5）表示抗體pi和pj的相似度，L為抗體基因編碼長(zhǎng)度；式（4）中分別表示第i個(gè)抗體基因上的第k位基因位和第j個(gè)抗體基因上的第k位基因位。

1.1.4 計(jì)算聚合適應(yīng)度 將親和度和抗體濃度融合成聚合適應(yīng)度。在種群更新的過(guò)程中，抗體被選擇的概率與聚合適應(yīng)度成正比。因此，聚合適應(yīng)度也反映了親和度與抗體濃度在種群更新中的影響力。親和度越高的個(gè)體，聚合適應(yīng)度也越高；抗體濃度越高的個(gè)體應(yīng)被抑制，以保證種群的多樣性，因此聚合適應(yīng)度越低。本文設(shè)定的聚合適應(yīng)度PZ為：

式中λ、μ為某一常數(shù)，分別反映了親和度和抗體濃度期望被選擇到下一代的相對(duì)重要的參數(shù)。

1.1.5 遺傳操作 通過(guò)輪盤(pán)賭選擇法進(jìn)行基因的復(fù)制。在輪盤(pán)賭中，種群中個(gè)體被選擇的概率與對(duì)應(yīng)的聚合適應(yīng)度成正比，即聚合適應(yīng)度越高的個(gè)體被選擇的概率越大，反之，聚合適應(yīng)度越低的個(gè)體被選擇的概率越小。第i個(gè)個(gè)體被選擇的概率Pi可表示為：

式中，PZi和PZj分別為第i個(gè)個(gè)體和第j個(gè)個(gè)體所對(duì)應(yīng)的聚合適應(yīng)度，n為種群大小。通過(guò)單點(diǎn)交叉對(duì)種群個(gè)體進(jìn)行交叉處理。設(shè)定交叉概率為pc，單點(diǎn)交叉操作即交換兩個(gè)個(gè)體的部分基因值。最后設(shè)定變異概率pm，對(duì)于0－1基因鏈碼，隨機(jī)挑選c個(gè)基因位，以變異概率pm對(duì)這些基因位上的值取反，即0變1，1變0。

1.1.6 更新迭代 迭代次數(shù)達(dá)到K時(shí)，選取種群中最優(yōu)個(gè)體作為最終候選個(gè)體，其對(duì)應(yīng)的波段作為最終篩選波段。

1.2 iIGA算法

1.2.1 自適應(yīng)相似度閾值 原算法中固定了抗體相似度閾值，而每次更新迭代，抗體相似度都會(huì)發(fā)生變化，這可能會(huì)導(dǎo)致閾值過(guò)大或者過(guò)小而使抗體濃度不具備差異性，使得引入的抗體濃度失去意義。本文提出自適應(yīng)相似度閾值，將相似度閾值以一定的間隔從群體相似度最小值遍歷至群體相似度最大值，所選取的相似度閾值使得群體抗體濃度的標(biāo)準(zhǔn)差達(dá)到最大，可完全體現(xiàn)出抗體濃度的差異性。

1.2.2 精英保留策略 在種群中另外設(shè)置M個(gè)個(gè)體，作為算法模型中的記憶細(xì)胞角色，并對(duì)其基因進(jìn)行0－1 隨機(jī)編碼，同時(shí)計(jì)算其親和度。由于種群的選擇具有不確定性，可能造成最優(yōu)個(gè)體的丟失。因此將親和度最高的M個(gè)個(gè)體，作為記憶細(xì)胞儲(chǔ)存起來(lái)，直接以概率1 復(fù)制到下一代種群中。在每一次迭代中，將當(dāng)前代數(shù)中親和度最高的M個(gè)個(gè)體，存入到記憶細(xì)胞中以進(jìn)行記憶細(xì)胞更新。

1.2.3 引入貪心算法思想 貪心算法在求最優(yōu)解時(shí)，不必考慮整體最優(yōu)，而是尋找在某種意義上僅是局部最優(yōu)的解，用局部最優(yōu)構(gòu)造全局最優(yōu)，一步一步逼近全局最優(yōu)解。在算法的每一次更新迭代中，對(duì)親和度最優(yōu)個(gè)體基因進(jìn)行局部性探優(yōu)，即隨機(jī)取反一個(gè)基因位并且保證探優(yōu)之后的親和度大于原先個(gè)體的親和度，如果小于，則恢復(fù)被改變的基因位，再隨機(jī)重取一個(gè)基因位進(jìn)行探優(yōu)，直到探優(yōu)之后該個(gè)體的親和度大于探優(yōu)前個(gè)體的親和度。如果重取次數(shù)達(dá)到上限值，則將原先最優(yōu)個(gè)體直接替換最差個(gè)體。

1.3 iIGA算法實(shí)現(xiàn)步驟

Step1：初始化種群。在種群中設(shè)置N個(gè)個(gè)體，作為抗體角色。將光譜波段均分并進(jìn)行0－1 隨機(jī)編碼。

Step2：計(jì)算親和度。對(duì)抗體基因進(jìn)行編碼后，將選中的波段與對(duì)應(yīng)的理化值進(jìn)行PLS 建模，并利用式（1）計(jì)算親和度。

Step3：初始化記憶細(xì)胞。在種群中設(shè)置M個(gè)個(gè)體，作為算法模型的記憶細(xì)胞角色，并對(duì)其基因進(jìn)行0－1隨機(jī)編碼，同時(shí)計(jì)算親和度。

Step4：計(jì)算抗體相似度。采用海明距離作為抗體相似度的度量方法，公式如式（3）～（5）所示。

Step5：自適應(yīng)相似度閾值計(jì)算抗體濃度。采用自適應(yīng)相似度閾值對(duì)抗體進(jìn)行相似判定，將群體抗體濃度標(biāo)準(zhǔn)差達(dá)到最大時(shí)，即抗體濃度差異性最大時(shí)的閾值作為最終自適應(yīng)相似度閾值?？贵w濃度由式（2）計(jì)算。

Step6：計(jì)算聚合適應(yīng)度。如式（6）所示，將親和度和抗體濃度融合成聚合適應(yīng)度。

Step7：引入貪心算法思想。對(duì)每一次迭代中親和度最優(yōu)的個(gè)體基因進(jìn)行局部性探優(yōu)，即隨機(jī)取反一個(gè)基因位并保證探優(yōu)之后的親和度大于原先的個(gè)體。

Step8：更新記憶細(xì)胞。重新計(jì)算群體中個(gè)體的親和度，選取最高的M個(gè)個(gè)體作為記憶細(xì)胞儲(chǔ)存起來(lái)，直接以概率為1復(fù)制到下一代種群中。

Step9：進(jìn)行遺傳操作。通過(guò)輪盤(pán)賭選擇法進(jìn)行基因的復(fù)制。在輪盤(pán)賭中，個(gè)體被選擇的概率如式（7）所示。然后通過(guò)單點(diǎn)交叉、變異對(duì)種群個(gè)體進(jìn)行處理。

Step10：更新迭代。迭代次數(shù)達(dá)K時(shí)，選取種群中最優(yōu)個(gè)體作為最終候選個(gè)體，其對(duì)應(yīng)的波段作為最終篩選波段。

2 實(shí)驗(yàn)部分

2.1 實(shí)驗(yàn)數(shù)據(jù)

2.1.1 數(shù)據(jù)來(lái)源 本研究引用于eigenvector 網(wǎng)站上開(kāi)源的一組標(biāo)準(zhǔn)玉米近紅外光譜數(shù)據(jù)集。該數(shù)據(jù)集由80 個(gè)玉米樣品組成，分別用m5spec、mp5spec、mp6spec 3 種不同的近紅外光譜儀進(jìn)行測(cè)試，包含了其淀粉、蛋白質(zhì)、水分以及油脂含量。以間隔為2 nm在1 100～2 498 nm波長(zhǎng)范圍上收集（700 波長(zhǎng)點(diǎn)）。本文選用mp5spec 儀器掃描得到的光譜數(shù)據(jù)，并選擇其中的淀粉和蛋白質(zhì)含量進(jìn)行研究。

2.1.2 樣本劃分 圖1 為原始光譜圖，數(shù)據(jù)集由80 個(gè)玉米樣品光譜組成。

圖1 原始光譜圖Fig.1 Original spectrogram

采用Kennard－Stone（KS）方法將80個(gè)樣本分成校正集和預(yù)測(cè)集。劃分結(jié)果為校正集樣本52個(gè)，預(yù)測(cè)集樣本28 個(gè)。校正集和預(yù)測(cè)集的淀粉和蛋白質(zhì)含量值統(tǒng)計(jì)見(jiàn)表1。由表1 可知，校正集樣本與預(yù)測(cè)集樣本的平均值（Avg）和標(biāo)準(zhǔn)差（Std）相差不大，通過(guò)KS 方法劃分的數(shù)據(jù)集，保證了校正集樣本均勻分布。

表1 校正集和預(yù)測(cè)集中淀粉和蛋白質(zhì)含量值統(tǒng)計(jì)Table 1 Statistics of starch and protein content in correction set and prediction set（g/100g）

2.2 實(shí)驗(yàn)設(shè)備

本研究使用的是一臺(tái)戴爾計(jì)算機(jī)，處理器是Intel（R）Core i5－9400，CPU 主頻為2.90 GHz，Windows10、64位操作系統(tǒng)，顯卡為GT710，所有計(jì)算均在MATLAB 2016a中進(jìn)行。

2.3 實(shí)驗(yàn)方法

為了考察本文提出的iIGA 算法對(duì)建模預(yù)測(cè)的效果，驗(yàn)證該算法的有效性和優(yōu)越性，分別用全譜波長(zhǎng)、GA、IGA 以及iIGA 算法對(duì)數(shù)據(jù)集中玉米淀粉和蛋白質(zhì)含量的建模效果進(jìn)行了對(duì)比。由于原始光譜波長(zhǎng)點(diǎn)數(shù)較多，如對(duì)波長(zhǎng)點(diǎn)進(jìn)行優(yōu)選組合，運(yùn)算效率將會(huì)非常低，因此將原光譜均分為若干個(gè)子區(qū)間，分別用上述算法進(jìn)化迭代，獲取最大適應(yīng)度值所對(duì)應(yīng)的優(yōu)選子區(qū)間組合。根據(jù)基因優(yōu)選的波段建立PLS 模型，計(jì)算出模型的相關(guān)系數(shù)（R）和預(yù)測(cè)集均方根誤差（RMSEP），適應(yīng)度函數(shù)F為：

將原光譜在1 100～2 498 nm 范圍之間的700 個(gè)波長(zhǎng)點(diǎn)數(shù)劃分為175 個(gè)等距區(qū)間，即遺傳編碼長(zhǎng)度為175，每一個(gè)基因位包含4個(gè)波長(zhǎng)點(diǎn)數(shù)。設(shè)定群體個(gè)數(shù)為50個(gè)，交叉概率為0.85，變異概率為0.05，最大迭代數(shù)為100代。在種群進(jìn)化過(guò)程中尋找最大迭代次數(shù)內(nèi)進(jìn)化過(guò)程的最優(yōu)適應(yīng)度個(gè)體。

對(duì)于iIGA 算法，通過(guò)實(shí)驗(yàn)選取的最佳記憶細(xì)胞數(shù)目為8 個(gè)。相似度閾值從抗體相似度最小值遍歷至抗體相似度最大值，遍歷跨度設(shè)置為0.001，當(dāng)所選取的相似度閾值使得群體抗體濃度的標(biāo)準(zhǔn)差達(dá)到最大時(shí)，即為最終的相似度閾值。經(jīng)過(guò)多次實(shí)驗(yàn)測(cè)試，設(shè)定聚合適應(yīng)度最佳參數(shù)λ為0.7，μ為1.25。設(shè)定貪心尋優(yōu)時(shí)的重取次數(shù)上限為20 次。由于上述算法存在隨機(jī)性，因此將以上3種算法分別運(yùn)行50次求平均。程序地址為https：//github. com/thmguai/Journal_of_analytical_testing_thm.git。

3 結(jié)果與討論

3.1 算法優(yōu)化結(jié)果

3.1.1 自適應(yīng)相似度閾值 相似度閾值從群體相似度最小值遍歷至群體相似度最大值。如圖2 所示，以玉米淀粉含量數(shù)據(jù)集為例，每取一個(gè)閾值，都會(huì)有與之一一對(duì)應(yīng)的群體抗體濃度標(biāo)準(zhǔn)差（STD），當(dāng)所選取的閾值使得標(biāo)準(zhǔn)差達(dá)到最大值時(shí)，表示群體抗體濃度差異性達(dá)到最大，可使IGA算法中提出的保留物種多樣性思想（引入抗體濃度）最大化。

圖2 群體抗體濃度標(biāo)準(zhǔn)差隨抗體相似度閾值的變化Fig.2 Variation of antibody concentration STD with the threshold of antibody similarity

3.1.2 精英保留策略 將親和度最高的8 個(gè)個(gè)體作為記憶細(xì)胞儲(chǔ)存起來(lái)，直接以概率為1 復(fù)制到下一代種群中。如圖3 所示，以玉米淀粉含量數(shù)據(jù)集為例，記憶細(xì)胞中儲(chǔ)存的精英群體的親和度均值遠(yuǎn)高于總?cè)后w，并且隨著迭代次數(shù)逐漸提升，可防止某些最優(yōu)個(gè)體的丟失，使種群朝著正確的方向進(jìn)行更新，有效地避免了陷入局部最優(yōu)。

圖3 精英群體和總?cè)后w的抗體親和度均值隨迭代次數(shù)的變化圖Fig.3 Variation of the mean antibody affinity of elite group and general population with the number of iterations

3.1.3 貪心算法思想 在算法的每一次迭代中，對(duì)抗體親和度最優(yōu)個(gè)體進(jìn)行局部性探優(yōu)，以局部最優(yōu)逐漸構(gòu)建出全局最優(yōu)。如圖4 所示，以玉米淀粉含量數(shù)據(jù)集為例，在對(duì)親和度最優(yōu)抗體進(jìn)行局部尋優(yōu)后，群體中最優(yōu)個(gè)體較之前幾乎都有了提升，突破了局部探優(yōu)之前的最優(yōu)領(lǐng)域，避免了過(guò)早陷入局部領(lǐng)域最優(yōu)，使算法朝著通往全局最優(yōu)的較合理的路徑進(jìn)行更新迭代。

圖4 局部探優(yōu)前后群體中最優(yōu)抗體的親和度隨迭代次數(shù)的變化圖Fig.4 Variation of the affinity of the optimal antibody with the number of iterations before and after local optimization

3.2 結(jié)果比較

iIGA 算法合理地優(yōu)化了算法的性能，同時(shí)增強(qiáng)了算法全局尋優(yōu)能力。圖5 為用GA、IGA 以及iIGA算法篩選之后的變量分別與淀粉和蛋白質(zhì)含量測(cè)試50次進(jìn)行建模分析的預(yù)測(cè)均方根誤差值。

由圖5 可知，用iIGA 算法進(jìn)行建模，預(yù)測(cè)玉米淀粉和蛋白質(zhì)含量的精度較IGA 算法有著明顯的提升，證明了該算法較原算法提升了性能，增強(qiáng)了全局尋優(yōu)能力。表2和表3分別為各算法測(cè)試淀粉含量和蛋白質(zhì)含量的各項(xiàng)性能數(shù)據(jù)。

圖5 GA、IGA以及iIGA算法運(yùn)行50次的淀粉預(yù)測(cè)RMSEP值（A）和蛋白質(zhì)預(yù)測(cè)RMSEP值（B）Fig.5 Starch predicted RMSEP value（A）and protein predicted RMSEP value（B）after GA，IGA and iIGA run 50 times

表2 不同變量篩選算法對(duì)玉米淀粉含量的預(yù)測(cè)結(jié)果Table 2 Prediction results of corn starch content by different variable screening methods

表3 不同變量篩選算法對(duì)玉米蛋白質(zhì)含量的預(yù)測(cè)結(jié)果Table 3 Prediction results of corn protein content by different variable screening methods

從表2可得，在玉米淀粉含量的預(yù)測(cè)上，相比于全譜PLS（F－PLS），GA 算法模型的校正集均方根誤差RMSEC從0.335 7降到0.272 9，校正集預(yù)測(cè)精度提升18.7%。模型的預(yù)測(cè)集均方根誤差RMSEP從0.391 4 下降到0.319 3，預(yù)測(cè)集預(yù)測(cè)精度提升18.4%。而IGA 算法作為GA 的改進(jìn)算法，預(yù)測(cè)精度有所提升。IGA 算法的RMSEC 為0.267 3，RMSEP 為0.312 0，預(yù)測(cè)精度分別較GA 算法提升了2.0%和2.3%。本文提出的iIGA 算法相較于IGA 算法，預(yù)測(cè)速度有所降低，但校正集相關(guān)系數(shù)Rc從0.941 3 提升至0.950 9，預(yù)測(cè)集相關(guān)系數(shù)Rp從0.920 4提升至0.927 7，RMSEC從0.267 3降到0.244 9，校正集預(yù)測(cè)精度提升8.4%；RMSEP 從0.312 0降到0.298 0，預(yù)測(cè)集預(yù)測(cè)精度提升4.5%。對(duì)IGA－PLS和iIGAPLS兩組分別測(cè)試50次的RMSEP值進(jìn)行顯著性檢驗(yàn)，F(xiàn)值為165.22，P值為9.5×10－23，P值小于0.05，因此兩組數(shù)據(jù)存在顯著性差別。由表2 數(shù)據(jù)分析可得，在淀粉含量光譜數(shù)據(jù)集中，iIGA 算法相較于原算法預(yù)測(cè)精度有了明顯的提升。

從表3 可得，在玉米蛋白質(zhì)含量的預(yù)測(cè)上，相比于全譜PLS，GA 算法模型的校正集均方根誤差RMSEC 從0.147 4 降到0.116 0，校正集預(yù)測(cè)精度提升21.3%。模型的預(yù)測(cè)集均方根誤差RMSEP 從0.178 9 下降到0.136 8，預(yù)測(cè)集預(yù)測(cè)精度提升23.5%。而IGA 算法的RMSEC 為0.112 4，RMSEP 為0.124 4，預(yù)測(cè)精度分別較GA 算法提升了3.1%和9.1%。本文提出的iIGA 算法相較于IGA 算法，預(yù)測(cè)速度有所降低，但校正集相關(guān)系數(shù)Rc從0.972 2 提升至0.976 4，預(yù)測(cè)集相關(guān)系數(shù)Rp從0.966 9 提升至0.974 2，RMSEC 從0.112 4降至0.103 7，校正集預(yù)測(cè)精度提升7.7%，RMSEP 從0.124 4降至0.110 3，預(yù)測(cè)集預(yù)測(cè)精度提升11.3%。同樣對(duì)IGA－PLS 和iIGA－PLS 兩組分別測(cè)試50 次的RMSEP 值進(jìn)行顯著性檢驗(yàn)，F(xiàn)值為182.05，P值為4.5×10－24，P值小于0.05，因此兩組數(shù)據(jù)存在顯著性差別。由表3數(shù)據(jù)分析可得，iIGA算法應(yīng)用在蛋白質(zhì)含量光譜數(shù)據(jù)集中相較于原算法預(yù)測(cè)精度有了明顯的提升。

4 結(jié) 論

提出了一種近紅外光譜波長(zhǎng)變量選擇算法，即改進(jìn)的免疫遺傳算法，對(duì)原算法進(jìn)行了3 個(gè)方面的改進(jìn)優(yōu)化，分別是自適應(yīng)相似度閾值、引入精英保留策略和貪心算法思想。采用該算法進(jìn)行波長(zhǎng)選擇后，建立了更加穩(wěn)健、預(yù)測(cè)精度更高的PLS 模型。在對(duì)玉米光譜數(shù)據(jù)集中淀粉和蛋白質(zhì)含量的預(yù)測(cè)建模中，模型的預(yù)測(cè)均方根誤差RMSEP 分別為0.298 0和0.110 3，相比原算法，模型的預(yù)測(cè)精度分別提高了4.5%和11.3%，證實(shí)了iIGA 算法的優(yōu)越性。在要求預(yù)測(cè)精度的場(chǎng)景下，本文提出的改進(jìn)IGA波長(zhǎng)變量選擇算法是一個(gè)更佳的選擇。