王文秀 彭彥昆 鄭曉春 孫宏偉 田 芳 白 京

(1.中國農(nóng)業(yè)大學(xué)工學(xué)院， 北京 100083； 2.國家農(nóng)產(chǎn)品加工技術(shù)裝備研發(fā)分中心， 北京 100083)

便攜式豬肉營養(yǎng)組分無損實時檢測裝置研究

王文秀1,2彭彥昆1,2鄭曉春1,2孫宏偉1,2田 芳1,2白 京1,2

(1.中國農(nóng)業(yè)大學(xué)工學(xué)院， 北京 100083； 2.國家農(nóng)產(chǎn)品加工技術(shù)裝備研發(fā)分中心， 北京 100083)

為了實現(xiàn)豬肉營養(yǎng)組分(脂肪和蛋白質(zhì))的快速、無損、實時檢測，基于近紅外反射光譜設(shè)計了便攜式豬肉營養(yǎng)組分無損檢測裝置。硬件部分包括光譜采集單元、光源單元和控制單元，并開發(fā)了相應(yīng)的檢測軟件，實現(xiàn)樣品光譜信息的有效獲取和實時分析。為了建立穩(wěn)定可靠的預(yù)測模型，考察了波段選擇、樣本分組方式和篩選變量方法對模型的影響。分別基于可見/短波近紅外(Vis/SWNIR)、長波近紅外(LWNIR)及Vis/SWNIR-LWNIR，利用隨機選擇法(RS)、Kennard-Stone法(KS)和基于聯(lián)合X-Y距離的樣本劃分法(SPXY)對樣本進行劃分，建立了脂肪和蛋白質(zhì)質(zhì)量分?jǐn)?shù)的偏最小二乘預(yù)測模型。結(jié)果發(fā)現(xiàn)，基于Vis/SWNIR-LWNIR波段，利用SPXY算法進行樣本分組，取得了最佳的預(yù)測模型。在此基礎(chǔ)上，比較分析競爭性自適應(yīng)加權(quán)算法、隨機蛙跳算法和蒙特卡羅無信息變量消除-連續(xù)投影算法3種算法篩選變量建立的模型效果?；诟偁幮宰赃m應(yīng)加權(quán)算法篩選變量的模型結(jié)果最佳，對脂肪和蛋白質(zhì)建立的模型驗證集相關(guān)系數(shù)分別為0.950 5和0.951 0。結(jié)果表明：基于近紅外反射光譜設(shè)計的便攜式豬肉組分檢測裝置可以對脂肪和蛋白質(zhì)含量進行快速、無損、實時檢測。

豬肉； 營養(yǎng)組分； 無損檢測； 便攜式裝置； 近紅外反射光譜

引言

豬肉營養(yǎng)價值豐富，是優(yōu)質(zhì)蛋白質(zhì)的重要來源，含有人體必需的各種氨基酸，同時也是脂肪能量的主要供給源和脂溶性維生素的重要載體[1]。隨著生活水平的提高和健康意識的不斷加強，人們對生鮮肉的質(zhì)量和品質(zhì)提出了更高的要求，低脂肪、高蛋白的豬肉日益受到人們的青睞[2-3]。依據(jù)營養(yǎng)及品質(zhì)對肉分級，可提高肉品行業(yè)的競爭力，也能確保消費者的滿意度。因此，快速、無損、實時獲取豬肉營養(yǎng)組分信息具有重要的意義。

基于近紅外光譜的檢測方法具有快速、無損、實時、原位的優(yōu)點，國內(nèi)外許多學(xué)者利用該方法對營養(yǎng)成分開展了相關(guān)的研究[4-7]。BARBIN等[8]對完整豬肉和絞碎后豬肉的脂肪和蛋白質(zhì)作了對比研究，利用近紅外反射光譜對絞碎后豬肉建立的偏最小二乘模型驗證集決定系數(shù)R2分別為0.95和0.86。劉魁武等[9]對不同溫度下貯藏的豬肉脂肪和蛋白質(zhì)進行了研究，發(fā)現(xiàn)0～4℃下建立的模型優(yōu)于20℃下模型，相關(guān)系數(shù)分別為0.950和0.713。LIAO等[10-11]對脂肪和蛋白質(zhì)含量進行了在線預(yù)測，經(jīng)多元散射校正和一階導(dǎo)數(shù)處理后建立了最佳的預(yù)測模型，驗證集決定系數(shù)R2分別為0.767和0.757。

綜合現(xiàn)有研究，豬肉營養(yǎng)組分的檢測以靜態(tài)實驗室研究居多，在線檢測的很少，便攜式檢測的幾乎沒有。并且就目前的檢測結(jié)果來看，將樣品絞碎后建立的模型多優(yōu)于完整肉樣的模型結(jié)果，在不損傷樣本的情況下建立的模型效果尚待進一步提高。為了滿足豬肉行業(yè)生產(chǎn)、加工和消費者的需求，實現(xiàn)豬肉脂肪和蛋白質(zhì)含量等營養(yǎng)成分的無損、快速、實時、原位、高通量檢測，本文設(shè)計便攜式豬肉營養(yǎng)成分無損實時檢測裝置，并結(jié)合化學(xué)計量學(xué)，建立脂肪和蛋白質(zhì)含量的無損定量預(yù)測模型。

1 材料與方法

1.1 試驗裝置

針對豬肉化學(xué)成分的檢測，本文設(shè)計了波長范圍覆蓋可見/短波近紅外(Visble/short-wavelength near infrared, Vis/SWNIR)和長波近紅外(Long-wavelength near infrared, LWNIR)的便攜式豬肉營養(yǎng)成分無損實時檢測裝置。該裝置包括光譜采集單元、光源單元和控制單元3部分。其中，光譜采集單元包括：輕巧便攜、操作簡單的USB4000型微型光纖光譜儀(美國Ocean Optics公司)，響應(yīng)波長范圍345～1 033 nm；高性能近紅外FLAMENIR型微型光纖光譜儀(美國Ocean Optics公司)，響應(yīng)波長范圍是945～1 665 nm；定制石英采集光纖。針對豬肉樣品不均勻，完整樣品建立模型效果不佳的問題，本研究采用環(huán)形光源設(shè)計，保證采集光纖收集更多位置處的反射光，進而獲取更多反映樣品營養(yǎng)組分的光譜信息。同時，定制石英光導(dǎo)，優(yōu)化光導(dǎo)末端與采集光纖的距離和角度，以避免光不進入樣品而直接反射進入采集光纖。整個裝置工作原理圖如圖1所示。

圖1 檢測裝置工作原理圖Fig.1 Working principle diagram of detection device1.計算機 2.光譜儀 3.采集光纖 4.環(huán)形光導(dǎo) 5.樣品 6.載物臺 7.光源室 8.鹵鎢燈 9.遮光板

豬肉營養(yǎng)成分無損實時檢測裝置軟件基于Visual Studio 2010平臺開發(fā)，采用C語言編程。該軟件的主要功能模塊如圖2所示，主要用于雙波段光譜信息的有效實時采集，對光譜信息進行處理(截去兩端噪聲，平滑及標(biāo)準(zhǔn)化處理)后，代入到預(yù)先植入的預(yù)測模型中，對豬肉營養(yǎng)組分進行實時預(yù)測，同時對樣品的光譜曲線和預(yù)測結(jié)果實時顯示，實現(xiàn)了一鍵式操作、實時人機界面、檢測結(jié)果統(tǒng)計和同步保存等功能。

圖2 檢測軟件功能模塊圖Fig.2 Function modules diagram of detection software

1.2 試驗材料

試驗樣品為經(jīng)過24 h冷卻排酸的新鮮豬肉，分別購買于北京3家不同的超市，修整為10 cm×5 cm×2.5 cm(長×寬×高)的肉樣后，利用蓄冷裝置運送至中國農(nóng)業(yè)大學(xué)無損檢測實驗室，放置于4℃冰箱中待用。為了建立穩(wěn)健可靠的預(yù)測模型，應(yīng)盡可能使脂肪和蛋白質(zhì)覆蓋大的含量范圍，因此本研究中分別選取背最長肌樣本30個、后臀尖樣本45個、純瘦肉樣本25個用于脂肪和蛋白質(zhì)定量預(yù)測模型的建立。試驗時，將樣品從冰箱取出，置于室溫(20℃)20～30 min，至樣品恢復(fù)至室溫并使表面水分蒸發(fā)，來減少對光譜采集的影響。

1.3 光譜信息采集

首先，打開便攜式豬肉營養(yǎng)組分無損實時檢測裝置預(yù)熱30 min，并對檢測軟件初始化。然后，設(shè)置采集參數(shù)(積分時間、平均次數(shù)和平滑度)，積分時間的確定應(yīng)保證參比光強達到總體的80%以上。設(shè)置USB4000型光譜儀的采集參數(shù)為：積分時間400 ms，平均次數(shù)為1，平滑度為4；FLAMENIR型光譜儀的采集參數(shù)為：積分時間1 000 ms，平均次數(shù)為1，平滑度為4。然后，將探頭對準(zhǔn)聚四氟乙烯標(biāo)準(zhǔn)白板采集參比光譜，將黑色遮光板插入光源室，阻止鹵鎢燈光線進入導(dǎo)光纖，采集背景參比。然后，將探頭放置在樣品上，按下按鈕，對樣品的光譜信息進行采集。每個樣品采集5個不同位置處光譜信息，求取平均值后作為該樣品的最終光譜。

1.4 脂肪和蛋白質(zhì)測定

脂肪的測定參照GB/T 9695.7—2008，取3～5 g樣品，加入50 mL 2 mol/L鹽酸，煮至無明顯肉塊，過濾，濾渣105℃干燥1 h后用石油醚為抽提劑抽提6 h，將接收器103℃干燥1 h，冷卻稱量，重復(fù)干燥冷卻步驟直至質(zhì)量變化小于樣品質(zhì)量的0.1%[12]。蛋白質(zhì)的測定參照GB 5009.5—2010，利用全自動凱氏定氮儀對絞碎后的肉樣進行處理。

圖3 模型建立過程Fig.3 Model building process

1.5 數(shù)據(jù)處理與模型評價

近紅外光譜常受其他與待測組分無關(guān)的因素干擾，如儀器響應(yīng)、隨機噪聲、雜散光等，導(dǎo)致光譜曲線不重復(fù)或產(chǎn)生基線漂移，因此首先要對原始光譜信息進行預(yù)處理。本文采用懲罰偏最小二乘平滑方法濾除原始光譜中的噪聲信號，該方法的優(yōu)化目標(biāo)函數(shù)既包括原始光譜和去噪光譜的最小二乘項(用來控制擬合誤差)，又包括了最小二乘的懲罰項(用來限定去噪光譜的平滑程度)，可較好地保持信號真實度[13]。同時采用標(biāo)準(zhǔn)正態(tài)變量變換(Standard normal variate transformation, SNVT)對濾波去噪后的光譜數(shù)據(jù)做中心化和標(biāo)度化處理，消除固體顆粒大小、表面散射以及光程變化對漫反射光譜的影響。

為建立穩(wěn)健的預(yù)測分析模型，本文對波段選擇、校正集和驗證集的劃分方法以及變量篩選方法進行比較研究。采用隨機選擇(Random selection, RS)法、Kennard-Stone(KS)法和基于聯(lián)合X-Y距離的樣本劃分(Sample set partitioning based on jointX-Ydistances, SPXY)法將預(yù)處理后的全部樣本按照3∶1的比例劃分為校正集和驗證集，則校正集中含75個樣本，驗證集中包括25個樣本。然后，分別基于Vis/SWNIR、LWNIR及Vis/SWNIR-LWNIR建立脂肪和蛋白質(zhì)含量的偏最小二乘(Partial least square, PLS)預(yù)測模型。在最佳模型的基礎(chǔ)上，采用競爭性自適應(yīng)加權(quán)算法(Competitive adaptive reweighted sampling, CARS)、隨機蛙跳算法(Random Frog)和蒙特卡羅-無信息變量消除-連續(xù)投影算法(Monte Carlo-uninformative variable elimination-successive projections algorithm method, MC-UVE-SPA)對變量進行多次篩選，消除無關(guān)信息對模型的影響。光譜分析過程如圖3所示。利用相關(guān)系數(shù)R、校正集標(biāo)準(zhǔn)分析誤差(Standard error of calibration, SEC)和驗證集標(biāo)準(zhǔn)分析誤差(Standard error of prediction, SEP)對預(yù)測模型進行評價和驗證。上述分析過程在Matlab R2012a和復(fù)雜多變量數(shù)據(jù)智慧處理軟件系統(tǒng)(大連達碩信息技術(shù)有限公司)中完成。

2 結(jié)果與討論

2.1 樣品光譜與理化值分析

從圖3所有樣品原始光譜曲線很難看出不同部位樣品光學(xué)信息的差異，因此對Vis/SWNIR和LWNIR區(qū)域光譜數(shù)據(jù)作主成分分析，前6個主成分的累積貢獻率為99.39%，然后利用SPSS軟件中Kruskal-Wallis檢驗對前6個主成分作顯著性差異分析，其漸進顯著性為p<0.05，表明不同部位樣品的光譜信息存在差異，為建立穩(wěn)健模型奠定了基礎(chǔ)[14]。由于Vis/SWNIR區(qū)域光譜兩端噪聲信號較大，因此只截取494.29～951.79 nm之間2 438個波長下光譜信息進行后續(xù)分析。背最長肌、后臀尖和純瘦肉3個部位樣本在Vis/SWNIR和LWNIR區(qū)域的原始光譜曲線如圖4a所示，可見3類樣本總體變化趨勢一致，僅在540～580 nm之間有所不同，該波段與氧合肌紅蛋白和氧合血紅蛋白的吸收峰密切相關(guān)[15]。經(jīng)過平滑和SNVT預(yù)處理后的Vis/SWNIR區(qū)域和LWNIR區(qū)域光譜曲線如圖4b所示。脂肪和蛋白質(zhì)的理化值統(tǒng)計結(jié)果如表1所示。

圖4 樣品光譜曲線Fig.4 Sample spectrum curves

%

2.2 不同波段、不同分組方式建模結(jié)果

選擇合適的波段對建立穩(wěn)健模型至關(guān)重要，同時利用合理有效的方法選擇具有代表性的樣本構(gòu)成校正集也會影響到模型的適應(yīng)性和預(yù)測性能以及模型的精簡程度。基于Vis/SWNIR、LWNIR和Vis/SWNIR-LWNIR，本文利用RS法、KS算法和SPXY法3種方法將樣本劃分為校正集和驗證集。RS法遵循隨機采樣的原則，操作簡單但隨機性強。KS算法通過計算兩兩樣本之間的歐氏距離，保證校正集樣本在空間距離上分布均勻。SPXY法是在KS算法的基礎(chǔ)上發(fā)展起來的，且同時考慮y空間的權(quán)重，以此改善模型的預(yù)測能力[16]。不同光譜范圍、不同分組方式下脂肪和蛋白質(zhì)的建模結(jié)果如表2和表3所示。

分析表中脂肪和蛋白質(zhì)的建模結(jié)果可以發(fā)現(xiàn)，不管采用何種校正集驗證集劃分方式，基于Vis/SWNIR波段的建模結(jié)果均優(yōu)于LWNIR波段建模結(jié)果，而基于Vis/SWNIR-LWNIR建立的模型結(jié)果又優(yōu)于前兩者。這可能是由于脂肪和蛋白質(zhì)屬于含C—H鍵及N—H鍵的物質(zhì)，Vis/SWNIR-LWNIR波段相比單一Vis/SWNIR波段或LWNIR波段涵蓋了更多C—H鍵和N—H鍵的伸縮振動區(qū)域，包含了更多與之相關(guān)的重要信息。而Vis/SWNIR的有效光譜范圍是500～950 nm，LWNIR的有效波段是900～1 650 nm，前者包含的與C—H鍵和N—H鍵相關(guān)的信息多于后者，因此基于Vis/SWNIR波段的建模結(jié)果優(yōu)于基于LWNIR波段的結(jié)果。此外，SPXY算法分組后建立的模型結(jié)果優(yōu)于RS算法和KS算法分組后建模結(jié)果，對脂肪和蛋白質(zhì)的驗證集相關(guān)系數(shù)分別達到0.934 8和0.934 2。這可能是由于SPXY算法在分組過程中，同時將X變量和y變量考慮在內(nèi)，確保樣本在X和y空間具有相同的權(quán)重，因此選出了更加具有代表性的校正集樣本來訓(xùn)練預(yù)測模型。

表2 不同分組方式下脂肪的建模結(jié)果Tab.2 Modeling results for fat using different grouping methods

表3 不同分組方式下蛋白質(zhì)的建模結(jié)果Tab.3 Modeling results for protein using different grouping methods

2.3 不同篩選變量方法建模

從上述討論可知，基于Vis/SWNIR-LWNIR光譜范圍對脂肪和蛋白質(zhì)含量取得了最佳的模型效果。但是，脂肪和蛋白質(zhì)屬于2類不同的物質(zhì)，在近紅外光譜區(qū)域具有不同的吸收譜帶，因此繼續(xù)利用3種不同的變量篩選方法在已建立的最佳模型基礎(chǔ)上提取各個物質(zhì)的特征變量，提高模型性能。

圖5 對脂肪含量數(shù)據(jù)進行CARS篩選變量過程圖及特征變量分布圖Fig.5 Process diagram of variables selection using CARS for fat and variable distribution graph

2.3.1CARS算法篩選變量

CARS算法基于模型集群思想，通過自適應(yīng)重加權(quán)采樣篩選PLS模型中回歸系數(shù)絕對值大的變量，利用指數(shù)衰減函數(shù)去掉絕對值小的變量，采樣交互檢驗選出交互驗證均方根誤差最小的變量組合[17-20]。本研究中設(shè)置采樣次數(shù)為100次，最大潛變量數(shù)為10，預(yù)處理方法選擇“center”，并選擇十折交叉檢驗法。圖5a所示為對脂肪進行100次CARS篩選變量過程圖，可以看出隨著運行次數(shù)的增加，篩選變量的數(shù)量呈現(xiàn)由快到慢的變化趨勢，十折交互驗證均方根誤差開始時隨著無關(guān)變量的剔除誤差逐漸減小，然后誤差變大，表明與待測組分相關(guān)的變量被剔除。由各個變量回歸系數(shù)的變化情況可看出，絕對值較大的變量為最終篩選的關(guān)鍵變量。圖中垂直星點對應(yīng)的位置為交互驗證均方根誤差最小的一次，對應(yīng)的變量為最終入選的特征變量。利用該方法篩選脂肪和蛋白質(zhì)的特征變量，分別由最初的2 574個減少到56個和49個，特征變量分布圖如圖5b所示。

分析圖5中各個參數(shù)特征變量的分布，與脂肪相關(guān)的變量中730～760 nm與脂肪族烴的4級倍頻相關(guān)，850～950 nm為C—H鍵伸縮的第3級倍頻區(qū)域，1 300～1 400 nm為亞甲基C—H對稱伸縮振動的第2級組合頻區(qū)。與蛋白質(zhì)相關(guān)的變量多集中在800～850 nm和1 400～1 500 nm之間，前者與N—H鍵伸縮相關(guān)，后者為伯醇的1級倍頻區(qū)域[21]。從圖5中看出，與C—H鍵和N—H鍵相關(guān)的變量，Vis/SWNIR波段多于LWNIR波段，這也驗證了表2和表3中的建模結(jié)果，同時也說明CARS算法可以客觀地篩選出各個組分的特征變量。以特征變量為因變量，以脂肪和蛋白質(zhì)質(zhì)量分?jǐn)?shù)作為自變量，建立新的PLS模型，各參數(shù)驗證集模型結(jié)果如圖6所示。與使用全部變量相比，模型結(jié)果有所提高。

圖6 基于CARS算法篩選變量建模的驗證集結(jié)果Fig.6 Model results for prediction set based on selected variables using CARS

2.3.2Random Frog算法篩選變量

Random Frog算法采用改進的逆跳馬爾科夫鏈蒙特卡洛進行模型采樣，得到多個不同的子模型，統(tǒng)計分析各個變量在模型中的選擇頻率，以此來評價該變量的重要程度，作為選擇變量的重要依據(jù)?；谡龖B(tài)分布隨機維數(shù)轉(zhuǎn)換機制的逆跳馬爾科夫改進方法，計算速度快，對模型的預(yù)測能力強[22]。本研究中設(shè)置運算次數(shù)為10 000次，交互驗證最大潛變量個數(shù)為10，預(yù)處理方法選擇“center”，評價指標(biāo)為回歸系數(shù)。圖7所示為利用Random Frog算法對蛋白質(zhì)含量數(shù)據(jù)提取特征變量的過程，設(shè)定0.1為閾值，閾值分割線以上變量為最終選擇的特征變量。

圖7 利用Random Frog算法對蛋白質(zhì)含量數(shù)據(jù)篩選特征變量Fig.7 Characteristic variables selection for protein using Random Frog algorithm

利用上述方法對脂肪和蛋白質(zhì)篩選特征變量，分別選出57和49個關(guān)鍵變量，占全部變量的2.22%和1.91%。以篩選的特征變量作為因變量，以脂肪和蛋白質(zhì)質(zhì)量分?jǐn)?shù)作為自變量，建立PLS模型，其驗證集結(jié)果如圖8所示。與基于全部變量建立的PLS模型相比，蛋白質(zhì)質(zhì)量分?jǐn)?shù)的預(yù)測模型性能有所提高，驗證集相關(guān)系數(shù)上升至0.948 4，脂肪質(zhì)量分?jǐn)?shù)的模型結(jié)果并無明顯改善，但是使用的變量數(shù)目僅是全部變量的2.22%，模型有了很大的簡化，表明該方法的可行性。

2.3.3MC-UVE-SPA算法篩選變量

UVE算法在原始光譜矩陣中加入一個相同大小的噪聲矩陣，采用刀切法基于交互驗證構(gòu)建多個PLS模型，基于各個變量的回歸系數(shù)bj，根據(jù)公式sj=mean(bj)/std(bj)，其中mean(bj)為系數(shù)向量的平均值，std(bj)為系數(shù)向量的標(biāo)準(zhǔn)偏差，以sj的絕對值作為評價指標(biāo)，顯著大于噪聲的變量保留下來作為特征變量[23]。本研究中設(shè)置最大潛變量數(shù)為10，預(yù)處理方法選擇“centter”，閾值為0.99。圖9a所示為利用UVE算法對蛋白質(zhì)質(zhì)量分?jǐn)?shù)篩選變量圖，圖中左側(cè)為真實變量，藍線兩側(cè)為入選的變量。由于該方法篩選出來的變量數(shù)較多，不能很好地簡化模型，因此在其基礎(chǔ)上繼續(xù)使用SPA算法篩選特征變量。圖9b為利用SPA算法篩選變量過程均方根誤差隨變量數(shù)變化圖。

圖8 基于Random Frog算法篩選變量建模的驗證集結(jié)果Fig.8 Model results for prediction set based on selected variables using Random Frog

圖9 利用MC-UVE-SPA算法對蛋白質(zhì)含量數(shù)據(jù)篩選特征變量Fig.9 Characteristic variables selection for protein using MC-UVE-SPA

利用SPA算法再次篩選變量后，大大降低了模型的復(fù)雜程度，建立的PLS模型驗證集結(jié)果如圖10所示。針對脂肪參數(shù)，共篩選出23個特征變量，占全部變量的0.89%，驗證集相關(guān)系數(shù)與使用全部變量相當(dāng)，但是誤差略有降低。針對蛋白質(zhì)參數(shù)，篩選出30個特征變量，占全部變量的1.17%，但是模型相關(guān)系數(shù)提高，預(yù)測誤差降低。

圖10 基于MC-UVE-SPA篩選變量建模的驗證集結(jié)果Fig.10 Model results for prediction set based on selected variables using MC-UVE-SPA

2.4 討論

從上述分析可知，基于3種方法篩選的特征變量建立的PLS模型效果均優(yōu)于使用全部變量建模結(jié)果，表明建立模型過程中篩選變量的重要性，這也符合“奧克姆剃刀”原則[24]，即以模型結(jié)果為導(dǎo)向，利用高效簡潔的特征變量來建立最為穩(wěn)健可靠的預(yù)測模型。從所選波長點數(shù)目來看，利用MC-UVE-SPA篩選到的變量數(shù)目最少，分別占全部變量的0.89%和1.17%，但其模型結(jié)果與Random Frog算法篩選到的2.22%和1.91%變量的建模結(jié)果相當(dāng)。就模型結(jié)果而言，CARS算法利用與Random Frog算法相當(dāng)?shù)淖兞繑?shù)目，給出了最優(yōu)的模型結(jié)果，具有最低的誤差。綜合比較3種方法，CARS算法能找到最佳的變量組合方式，取得最佳的模型效果，是一種比較穩(wěn)定的算法。

3 裝置試驗驗證

通過上述討論，找到了脂肪和蛋白質(zhì)質(zhì)量分?jǐn)?shù)的特征波長，并建立了PLS預(yù)測模型。為了檢驗?zāi)Ｐ偷姆€(wěn)定性和可靠程度，另取20塊豬肉樣本用于對脂肪和蛋白質(zhì)模型的驗證。將樣本從冰箱取出后，置于室溫30 min，利用便攜式豬肉營養(yǎng)組分無損檢測裝置實時采集樣本的反射光譜曲線。編寫的軟件自動實時進行光譜預(yù)處理，并通過內(nèi)置模型得到脂肪和蛋白質(zhì)質(zhì)量分?jǐn)?shù)的預(yù)測值。然后參照國家標(biāo)準(zhǔn)方法，對脂肪和蛋白質(zhì)質(zhì)量分?jǐn)?shù)進行測定，得到其標(biāo)準(zhǔn)值。試驗驗證結(jié)果如表4所示。

表4 試驗驗證結(jié)果Tab.4 Results of prediction tests

4 結(jié)論

(1)基于近紅外反射光譜設(shè)計了便攜式豬肉營養(yǎng)組分無損快速檢測裝置。該裝置由光譜采集單元、光源單元和控制單元等硬件部分和編寫的控制程序組成，可以實現(xiàn)對豬肉脂肪、蛋白質(zhì)等營養(yǎng)組分的無損、快速、實時檢測。

(2)利用該檢測裝置，分別基于Vis/SWNIR、LWNIR和Vis/SWNIR-LWNIR建立了脂肪和蛋白質(zhì)的偏最小二乘預(yù)測模型。結(jié)果表明，利用Vis/SWNIR波段建模結(jié)果優(yōu)于LWNIR，而Vis/SWNIR-LWNIR波段的建模結(jié)果最佳。

(3)為建立穩(wěn)健的模型，比較了RS算法、KS算法和SPXY算法3種分組方式以及CARS算法、Random Frog算法和MC-UVE-SPA 3種變量篩選方法對模型效果的影響。結(jié)果表明，SPXY算法優(yōu)于RS算法和KS算法，可以得到更加具有代表性的樣本作為校正集。在脂肪和蛋白質(zhì)質(zhì)量分?jǐn)?shù)的最佳模型基礎(chǔ)上，利用CARS篩選的變量建立更為簡化的偏最小預(yù)測模型取得了最佳的預(yù)測效果，具有最小的誤差，對脂肪和蛋白質(zhì)建立的模型驗證集相關(guān)系數(shù)分別為0.950 5和0.951 0。結(jié)果表明，該裝置可以實現(xiàn)對豬肉脂肪和蛋白質(zhì)的無損、快速、實時預(yù)測及分析，具有較高的應(yīng)用前景。

1 呂自治．豬肉的營養(yǎng)價值及其科學(xué)食用[J]．肉類研究，2003(2)：49-50．

2 王重龍，陶立，張東紅，等．安全優(yōu)質(zhì)豬肉的生產(chǎn)[J]．安徽農(nóng)學(xué)通報，2005，11(6)：121-123．

3 王笑丹．吉林省優(yōu)質(zhì)豬肉品質(zhì)評定方法[D]．長春：吉林大學(xué)，2004． WANG Xiaodan．The evaluation methods of high quality pork in Jilin [D]．Changchun：Jilin University，2004．(in Chinese)

4 成芳，樊玉霞，廖宜濤．應(yīng)用近紅外漫反射光譜對豬肉肉糜進行定性定量檢測研究[J]．光譜學(xué)與光譜分析，2012，32(2)：354-359． CHENG Fang，F(xiàn)AN Yuxia，LIAO Yitao．Qualitative and quantitative detection of minced pork quality by near infrared reflectance spectroscopy [J]．Spectroscopy and Spectral Analysis，2012，32(2)：354-359．(in Chinese)

5 樊玉霞．豬肉肉糜品質(zhì)與安全可見/近紅外光譜快速檢測方法的實驗研究[D]．杭州：浙江大學(xué)，2011． FAN Yuxia．Rapid detection of minced pork quality and safety based on visible and near infrared spectroscopy [D]．Hangzhou：Zhejiang University，2011．(in Chinese)

6 湯修映，牛力釗，徐楊，等．基于可見/近紅外光譜技術(shù)的牛肉含水率無損檢測[J]．農(nóng)業(yè)工程學(xué)報，2013，29(11)：248-254． TANG Xiuying，NIU Lizhao，XU Yang，et al．Nondestructive determination of water content in beef using visible/near-infrared spectroscopy [J]. Transactions of the CSAE，2013，29(11)：248-254．(in Chinese)

7 KAMRUZZAMAN M，MAKINO Y，OSHITA S．Non-invasive analytical technology for the detection of contamination, adulteration, and authenticity of meat, poultry, and fish: a review [J]．Analytica Chimica Acta，2015，853(1)：19-29．

8 BARBIN D F，ELMASRY G，SUN D W，et al．Predicting quality and sensory attributes of pork using near-infrared hyperspectral imaging [J]．Analytica Chimica Acta，2012，719(10)：30-42．

9 劉魁武，成芳，林宏建，等．可見/近紅外光譜檢測冷鮮豬肉中的脂肪、蛋白質(zhì)和水分含量[J]．光譜學(xué)與光譜分析，2009，29(1)：102-105． LIU Kuiwu，CHENG Fang， LIN Hongjian，et al．Visible/NIR analysis of fat, protein and water in chilled pork[J]．Spectroscopy and Spectral Analysis，2009，29(1)：102-105．(in Chinese)

10 LIAO Yitao，F(xiàn)AN Yuxia，CHENG Fang．On-line prediction of fresh pork quality using visible/near-infrared reflectance spectroscopy [J]．Meat Science，2010，86(4)：901-907．

11 廖宜濤，樊玉霞，伍學(xué)千，等．豬肉肌內(nèi)脂肪含量的可見/近紅外光譜在線檢測[J]．農(nóng)業(yè)機械學(xué)報，2010，41(9)：104-107． LIAO Yitao，F(xiàn)AN Yuxia，WU Xueqian，et al．On line prediction of intramuscular fat content in pork muscle with visible/near infrared spectroscopy [J]．Transactions of the Chinese Society for Agricultural Machinery，2010，41(9)：104-107．(in Chinese)

12 張偉力，曾勇慶．豬肉肌內(nèi)脂肪測定方法及其誤差分析[J]．豬業(yè)科學(xué)，2008，25(7)：102-103．

13 楊慧中，陳定三．局部懲罰加權(quán)核偏最小二乘算法及其應(yīng)用[J]．控制工程，2011，18(6)：886-889． YANG Huizhong，CHEN Dingsan．Local penalized weighted kernel partial least squares algorithm and its application[J]．Control Engineering of China，2011，18(6)：886-889．(in Chinese)

14 樊書祥，黃文倩，郭志明，等．蘋果產(chǎn)地差異對可溶性固形物近紅外光譜檢測模型影響的研究[J]．分析化學(xué)，2015，43(2)：239-244． FAN Shuxiang，HUANG Wenqian， GUO Zhiming，et al．Assessment of influence of origin variability on robustness of near infrared models for soluble solid content of apples[J]．Chinese Journal of Analytical Chemistry，2015，43(2)：239-244．(in Chinese)

15 COZZOLINO D，BARLOCCO N，VADELL A，et al． The use of visible and near-infrared reflectance spectroscopy to predict colour on both intact and homogenised pork muscle[J]．LWT—Food Science and Technology，2003，36(2)：195-202．

16 石力安．基于可見/近紅外光譜的生鮮牛肉水分快速無損檢測方法的研究[D]．北京：中國農(nóng)業(yè)大學(xué)，2016． SHI Li’an．Study on rapid and nondestructive detection method of beef moisture content based on VIS/NIR spectroscopic technique [D]．Beijing：China Agricultural University，2016．(in Chinese)

17 李江波，彭彥昆，陳立平，等．近紅外高光譜圖像結(jié)合CARS算法對鴨梨SSC含量定量測定[J]．光譜學(xué)與光譜分析，2014，34(5)：1264-1269． LI Jiangbo，PENG Yankun，CHEN Liping，et al．Near-infrared hyperspectral imaging combined with CARS algorithm to quantitatively determination soluble solids content in “Ya” pear[J]．Spectroscopy and Spectral Analysis，2014，34(5)：1264-1269．(in Chinese)

18 陸婉珍．現(xiàn)代近紅外光譜分析技術(shù)[M]．北京：中國石化出版社，2006．

19 呂程序，姜訓(xùn)鵬，張銀橋，等．基于變量選擇的小麥粗蛋白含量近紅外光譜檢測[J/OL]．農(nóng)業(yè)機械學(xué)報，2016，47(增刊)：340-346．http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=2016s052&journal_id=jcsam．DOI：10.6041/j.issn.1000-1298.2016.S0.052. Lü Chengxu，JIANG Xunpeng，ZHANG Yinqiao，et al．Variable selection based near infrared spectroscopic quantitative analysis on wheat crude protein content [J/OL]．Transactions of the Chinese Society for Agricultural Machinery，2016，47(Supp.):340-346．(in Chinese)

20 劉燕德，施宇，蔡麗君，等．基于CARS算法的臍橙可溶性固形物近紅外在線檢測[J/OL]．農(nóng)業(yè)機械學(xué)報，2013，44(9)：138-144．http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20130925&journal_id=jcsam．DOI：10.6041/j.issn.1000-1298.2013.09.025． LIU Yande，SHI Yu，CAI Lijun，et al．On-line NIR detection model optimization of soluble solids content in navel orange based on CARS [J/OL]．Transactions of the Chinese Society for Agricultural Machinery，2013，44(9)：138-144．(in Chinese)

21 梁逸曾，許青松．復(fù)雜體系儀器分析[M]．北京：化學(xué)工業(yè)出版社，2012．

22 陳立旦，趙艷茹．可見-近紅外光譜聯(lián)合隨機蛙跳算法檢測生物柴油含水量[J]．農(nóng)業(yè)工程學(xué)報，2014，30(8)：168-173． CHEN Lidan，ZHAO Yanru．Measurement of water content in biodiesel using visible and near infrared spectroscopy combined with Random-Frog algorithm [J]. Transactions of the CSAE，2014，30(8)：168-173. (in Chinese)

23 李倩倩，田曠達，李祖，等．無信息變量消除法變量篩選優(yōu)化煙草中總氮和總糖的定量模型[J]．分析化學(xué)，2013，41(6)：917-921． LI Qianqian，TIAN Kuangda，LI Zu，et al．Model of total nitrogen and total sugar in tobacco optimizing after uninformative variable elimination [J]．Chinese Journal of Analytical Chemistry，2013，41(6)：917-921．(in Chinese)

24 BLUMER A，EHRENFEUCHT A，HAUSSLER D，et al．Occam’s razor[J]．Information Processing Letters，1987，24(6)：377-380．

PortableNondestructiveDetectionDeviceforNutrientComponentsofPork

WANG Wenxiu1,2PENG Yankun1,2ZHENG Xiaochun1,2SUN Hongwei1,2TIAN Fang1,2BAI Jing1,2

(1.CollegeofEngineering,ChinaAgriculturalUniversity,Beijing100083,China2.NationalR&DCenterforAgro-processingEquipment,Beijing100083,China)

In order to realize fast, nondestructive and real-time detection of nutrition components (fat and protein) for pork, a portable nondestructive detection device based on near infrared reflectance spectra was designed and developed. The hardware part included spectrum acquisition unit, light source unit and control unit. The corresponding detection software was developed to realize the effective acquisition and real-time analysis of the sample spectrum information. In order to establish a stable and reliable forecasting model, the research focused on the effects of band selection, different sample grouping methods and variables selection methods on the models. Based on visible/short wavelength near infrared (Vis/SWNIR), long wavelength near-infrared (LWNIR) and Vis/SWNIR-LWNIR, all the samples were divided by random selection (RS) method, Kennard-Stone (KS) algorithm and sample set partitioning based on jointX-Ydistances (SPXY) algorithm, and then partial least square prediction models for fat and protein content were built, respectively. The results showed that the best prediction models for fat and protein were built based on Vis/SWNIR-LWNIR by using SPXY algorithm. On the basis of the best model for each parameter, comparative analysis of competitive adaptive weighted algorithm, Random Frog algorithm and uninformative variable elimination-successive projection algorithm were employed to screen variables. The results showed that the simplified model based on competitive adaptive weighting algorithm was the best with correlation coefficients in the prediction set of 0.950 5 and 0.951 0 for fat and protein, respectively. The results indicated that the designed portable detection device based on near infrared reflectance spectroscopy was able to realize fast, nondestructive and real-time detection of fat and protein content for fresh meat and had certain application potential and market prospects.

pork; nutrition component; nondestructive detection; portable device; near infrared reflectance spectroscopy

S379.9

A

1000-1298(2017)09-0303-09

10.6041/j.issn.1000-1298.2017.09.038

2017-02-17

2017-03-24

國家重點研發(fā)計劃項目(2016YFD0401205)和公益性行業(yè)(農(nóng)業(yè))科研專項(201003008)

王文秀(1989—)，女，博士生，主要從事農(nóng)畜產(chǎn)品無損檢測技術(shù)與裝置研究，E-mail: Godlovexiu@163.com

彭彥昆(1960—)，男，教授，博士生導(dǎo)師，主要從事農(nóng)畜產(chǎn)品無損檢測技術(shù)與裝置研究，E-mail: ypeng@cau.edu.cn