摘要：為更加快捷且不破壞樣品地對(duì)選后煙葉的化學(xué)成分（煙堿、總糖、還原糖、總氮、鉀、氯）值進(jìn)行檢測(cè)，有效地組織打葉復(fù)烤均值化加工，應(yīng)用選后煙葉中煙堿含量變化與其對(duì)應(yīng)的近紅外光譜變化之間存在相關(guān)關(guān)系，采用漫反射方式獲取校正樣品集近紅外漫反射光譜，選擇微分、平滑、PLS等方法進(jìn)行數(shù)據(jù)預(yù)處理，把校正樣品集近紅外光譜與其相應(yīng)的主要化學(xué)成分含量或指標(biāo)值進(jìn)行關(guān)聯(lián)，建立校正模型，使用驗(yàn)證樣品集對(duì)校正模型有效性進(jìn)行驗(yàn)證。通過(guò)對(duì)檢測(cè)結(jié)果進(jìn)行分析，該建模流程是可行的，為選后煙葉煙堿值的檢測(cè)提供了一種更加快捷且無(wú)損的方法。

關(guān)鍵詞：選后煙葉、煙堿、手持式、近紅外、建模

煙草是我國(guó)的重要經(jīng)濟(jì)作物之一，由于受品種、栽培技術(shù)、土壤、氣候、收購(gòu)眼光等因素影響，使同一產(chǎn)地、同一等級(jí)的煙葉在不同收購(gòu)時(shí)段內(nèi)在品質(zhì)上存在差異。打葉復(fù)烤均值化加工就是通過(guò)一定形式的參配，使打葉復(fù)烤后的煙葉形成一個(gè)內(nèi)在質(zhì)量穩(wěn)定、理化指標(biāo)均勻的成品片煙或配方模塊。在煙草的眾多質(zhì)量評(píng)價(jià)指標(biāo)中煙堿的含量，直接影響煙草制品評(píng)吸質(zhì)量，是配方設(shè)計(jì)的基礎(chǔ)控制因素。早在2015 年中國(guó)煙草總公司就通過(guò)文件明確要求要將成品片煙的煙堿變異系數(shù)控制到 5%以內(nèi)。卷煙加工企業(yè)調(diào)入打葉復(fù)烤企業(yè)加工的煙葉一般要經(jīng)過(guò)原煙入庫(kù)、煙葉分選、配方打葉等流程，而對(duì)選后煙葉煙堿含量的情況掌握是均質(zhì)化加工配方設(shè)計(jì)的關(guān)鍵環(huán)節(jié)。打葉復(fù)烤企業(yè)對(duì)煙葉中煙堿的檢測(cè)一般是使用連續(xù)流動(dòng)分析法，近年來(lái)隨著化學(xué)計(jì)量學(xué)與計(jì)算機(jī)技術(shù)的飛速發(fā)展，實(shí)驗(yàn)室近紅外檢測(cè)煙葉中煙堿的含量得到了廣泛應(yīng)用。但是不論是連續(xù)流動(dòng)分析法還是實(shí)驗(yàn)室近紅外檢測(cè)均需對(duì)樣品進(jìn)行粉碎，且數(shù)據(jù)獲得存在一定的滯后性。

本文通過(guò)對(duì)一種手持式近紅外光譜分析儀建模的研究，提供了一種不破壞樣品且更加快捷的煙葉中煙堿的檢測(cè)方法，有助于打葉復(fù)烤企業(yè)靈活快捷地組織均質(zhì)化加工工作。

1 材料與方法

1.1材料和儀器

1.1.1 材料

云南初烤煙葉（選后煙葉）總計(jì)350個(gè)作為常規(guī)化學(xué)指標(biāo)模型建模樣本，所有樣本產(chǎn)地包括云南昆明、曲靖、文山、麗江、臨滄、玉溪產(chǎn)區(qū)，主要等級(jí)包括B1F、B2F、B3F、C1F、C2F、C3F、C4F、X2F、X4F

1.1.2儀器

德國(guó)Carl Zeiss公司生產(chǎn)的AURA手持式近紅外光譜儀（以下簡(jiǎn)稱AURA）（如圖1）。

圖1 AURA手持式近紅外光譜儀

美國(guó)ThermoFisher公司生產(chǎn)的 Antaris II 傅里葉變換近紅外光譜儀（以下簡(jiǎn)稱Thermo）（如圖2），此儀器已經(jīng)完成建模，且使用穩(wěn)定。

圖2 Antaris II 傅里葉變換近紅外光譜儀

2 樣品準(zhǔn)備

根據(jù)石林復(fù)烤廠烤季生產(chǎn)原煙情況，選取不同產(chǎn)地、等級(jí)煙葉樣品350個(gè)，具體信息如下圖表1：

圖表1 取樣表：

每個(gè)樣品取6-8片煙葉，放入自封袋中，后續(xù)進(jìn)行手持式近紅外光譜儀建模掃樣，對(duì)樣品進(jìn)行統(tǒng)一編號(hào)。編號(hào)由取樣樣品編號(hào)、產(chǎn)地、等級(jí)組成。如：“ 001#-昆明-C4F”，代表001號(hào)樣品，產(chǎn)地為“昆明”，等級(jí)為“C4F”。

3 光譜采集

建模光譜掃描采用手持式(AURA)近紅外光譜儀，序列號(hào)為131517。手持式近紅外光譜儀有效光斑直徑為18mm，采用與煙葉表面直接接觸的掃描方式。

圖表2 接觸式掃描

如下圖表3所示，將6片煙葉疊放在一起，用手持式近紅外光譜儀（六點(diǎn)掃描）工作流，先掃描3次，分別是葉尖，葉腰，葉基（對(duì)應(yīng)圖表的點(diǎn)1、點(diǎn)2、點(diǎn)3）各掃描一次，得到3條光譜；之后將最下面的三片煙葉抽出放置最上面，再次進(jìn)行葉尖，葉腰，葉基三次掃描。掃描完畢后，每個(gè)煙葉樣本對(duì)應(yīng)6條光譜。保持樣本名稱與取樣的名稱一一對(duì)應(yīng)。6次掃描近紅外的外界環(huán)境保持相對(duì)穩(wěn)定。

圖表3常規(guī)化學(xué)模型檢測(cè)點(diǎn)

由于外界環(huán)境的變化會(huì)引起光譜發(fā)生變化，為了消除這些影響，在每個(gè)樣品測(cè)量之前，儀器將進(jìn)行背景校正，防止光譜發(fā)生偏移，且當(dāng)外部環(huán)境溫度變化超過(guò)2℃時(shí)，儀器會(huì)自動(dòng)進(jìn)行背景的校正。

4 模型建立

4.1 光譜預(yù)處理

手持式近紅外的采集光譜中除含有待測(cè)樣品的原始的化學(xué)信息外，還包含其他外在的干擾信息和噪聲，這些信息噪聲將導(dǎo)致測(cè)得的化學(xué)信息數(shù)值與真實(shí)值之間存在一定的差異。為了提高信噪比盡可能消除誤差，應(yīng)該保持實(shí)驗(yàn)的環(huán)境因素盡量穩(wěn)定一致，除了環(huán)境因素之外還需要各種數(shù)據(jù)處理方法來(lái)減少甚至去除各種干擾光譜信息的影響因素。為下一步的數(shù)據(jù)處理奠定基礎(chǔ)。較常用的光譜數(shù)據(jù)處理方法有均值中心化（mean centering），歸一化處理（normalization）,平滑，求導(dǎo)，標(biāo)準(zhǔn)正態(tài)變量變換（SNV），多元散射校正（MSC），傅立葉變化和其他一些新的方法，化學(xué)指標(biāo)模型建立的過(guò)程中，會(huì)用到一階導(dǎo)數(shù)的預(yù)處理方法。

?? 導(dǎo)數(shù)

導(dǎo)數(shù)（微分）可以消除樣品光譜的基線漂移、增強(qiáng)譜帶的特征信息并克服譜帶重疊，是最為常用的光譜數(shù)據(jù)預(yù)處理方法。一階導(dǎo)數(shù)可以去除同波長(zhǎng)無(wú)關(guān)的基線漂移；二階導(dǎo)數(shù)可以去除同波長(zhǎng)線性相關(guān)的基線漂移。在原煙的數(shù)據(jù)采集中，導(dǎo)數(shù)計(jì)算可以減少基線偏移、漂移和背景干擾造成的數(shù)據(jù)偏差。

4.2異常樣本剔除

異常樣本在近紅外整個(gè)模型的構(gòu)建過(guò)程中會(huì)對(duì)模型存在很大的影響，不僅會(huì)會(huì)誤導(dǎo)近紅外變量的選擇，而且這些奇異值還會(huì)給模型的參數(shù)估計(jì)帶來(lái)偏離，降低模型的精度與穩(wěn)健性；在實(shí)際近紅外構(gòu)建模型的過(guò)程中，有很多情況會(huì)引起近紅外樣本的異常，大體上分為近紅外光譜的異常與近紅外化學(xué)值的異常；按照引起異常樣本的原因又可具體的分為如下幾種情形：

?①環(huán)境引起的異常光譜，比如近紅外儀器背影變化產(chǎn)生的影響或者掃到了非信號(hào)物質(zhì)，檢測(cè)條件的變化，如樣本的溫度，濕度；

?②儀器自身的不穩(wěn)定引起的異常光譜。比如儀器的檢測(cè)器兩端做的比較粗糙所儀器的問(wèn)題；近紅外隨著長(zhǎng)時(shí)間的使用儀器的老化，以及更換近紅外儀器的零部件帶來(lái)的光譜差異；
??? ③被檢測(cè)物質(zhì)自身所引起的異常光譜；這類異常光譜比較特殊，這是近紅外應(yīng)用模型維護(hù)所必須需要考慮到的問(wèn)題，由于地域以及年度之間的差異導(dǎo)致了未來(lái)的光譜與模型里面產(chǎn)生的光譜不一致；
??? ④由于工藝參數(shù)的改變儀器的異常光譜，近紅外光譜受外界環(huán)境影響比較大，如改變某些工藝參數(shù)，導(dǎo)致近紅外儀器檢測(cè)條件發(fā)生變化，比如真空回潮等參數(shù)會(huì)改變煙葉物質(zhì)中所含有的水分；這個(gè)時(shí)候外界的濕度因素就會(huì)發(fā)生變化，產(chǎn)生異常光譜；
??? ⑤基礎(chǔ)數(shù)據(jù)產(chǎn)生的異常樣本，當(dāng)基礎(chǔ)數(shù)據(jù)操作失誤，或者樣本混淆的時(shí)候，化學(xué)值與樣本所對(duì)應(yīng)的近紅外光譜不一致所產(chǎn)生的異常樣本；

手持式近紅外異常光譜的常用剔除方法是半數(shù)重采樣法和蒙特卡洛偏最小二乘交叉檢驗(yàn)法。

?? 半數(shù)重采樣法

基于對(duì)原始光譜的隨機(jī)半數(shù)重采樣統(tǒng)計(jì)出現(xiàn)奇異長(zhǎng)度的樣本。從原始光譜矩陣中隨機(jī)選擇部分（一般選擇總樣本數(shù)的一半）樣本作為采樣子集，計(jì)算每個(gè)采樣子集矩陣的均值和方差，再根據(jù)均值和方差計(jì)算采樣子集中每個(gè)樣本的向量長(zhǎng)度（向量長(zhǎng)度計(jì)算公式與數(shù)據(jù)標(biāo)準(zhǔn)化公式相同)。對(duì)光譜數(shù)據(jù)進(jìn)行多次隨機(jī)采樣，并記錄每次采樣后計(jì)算的向量長(zhǎng)度。對(duì)樣本的向量長(zhǎng)度進(jìn)行排序，距離最大的一定概率（如 5%或 10%)的樣本得分為 1，其余為 0。最后對(duì)各樣本的總得分進(jìn)行統(tǒng)計(jì)，得分最高的部分樣本就為奇異樣本。

?? 蒙特卡洛偏最小二乘交叉檢驗(yàn)

基于蒙特卡羅交叉驗(yàn)證（MCCV)的一類奇異樣本識(shí)別方法。利用 MCCV 隨機(jī)劃分校校正集與預(yù)測(cè)集，如果奇異樣本在校正集中，整個(gè)模型的質(zhì)量將受到影響；相反，如果奇異樣本在預(yù)測(cè)集中，僅此樣本的預(yù)測(cè)結(jié)果受到影響。盡管這種情況對(duì)預(yù)測(cè)結(jié)果都有影響，但效果明顯不同。本文就利用奇異樣本出現(xiàn)在校正集或預(yù)測(cè)集時(shí)模型預(yù)測(cè)誤差的差異，通過(guò) MCCV 及統(tǒng)計(jì)分析來(lái)進(jìn)行奇異樣本的識(shí)別。根據(jù)預(yù)測(cè)集中奇異樣本的預(yù)測(cè)殘差會(huì)明顯大于正常樣本的預(yù)測(cè)殘差也提出了一種基于MCCV的奇異樣本識(shí)別方法?；?MCCV的奇異樣本識(shí)別方法充分利用統(tǒng)計(jì)學(xué)的性質(zhì)，能夠在一定程度上降低由掩蔽效應(yīng)帶來(lái)的風(fēng)險(xiǎn)，檢出光譜陣和性質(zhì)陣方向的奇異點(diǎn)，有望在奇異樣本檢測(cè)中得到更廣泛的應(yīng)用。

4.3 波長(zhǎng)變量選擇

在校正模型的建立過(guò)程中，選取參與校正的樣本和光譜信息變量對(duì)建立穩(wěn)定的模型時(shí)十分必要的。光譜信息變量選擇是從原始變量中挑選出一些有代表性的特征變量，代替原始變量進(jìn)行數(shù)據(jù)分析和處理。在煙草無(wú)損檢測(cè)試驗(yàn)中，近紅外光譜儀每次可獲取大量的光譜數(shù)據(jù)，應(yīng)著不同原始光譜數(shù)據(jù)對(duì)待測(cè)煙葉樣品的品質(zhì)信息的貢獻(xiàn)率不能完全相同，有些光譜信息變量反映的信息較為豐富，有些光譜信息變量反映的信息量較少，甚至與待測(cè)煙葉樣品品質(zhì)成分含量無(wú)關(guān)。如果將近紅外光譜儀所獲取的光譜數(shù)據(jù)都用于建立模型，則建模計(jì)算時(shí)間很長(zhǎng)，計(jì)算量也很大，建立的近紅外光譜預(yù)測(cè)模型復(fù)雜，模型的穩(wěn)定性差。研究發(fā)現(xiàn)，通過(guò)特定的光譜信息變量篩選方法對(duì)原始光譜信息變量進(jìn)行優(yōu)選，其作用是可減少模型的建立時(shí)間，簡(jiǎn)化建模過(guò)程，最重要的是可以剔除無(wú)信息變量或非線性變量，最終可以建立跟隨性強(qiáng)、預(yù)測(cè)能力好的近紅外定量校正模型。

?手持式近紅外常用的波長(zhǎng)選擇方法是CARS（競(jìng)爭(zhēng)性自適應(yīng)加權(quán)取樣法），具體的計(jì)算過(guò)程如下：

競(jìng)爭(zhēng)性自適應(yīng)權(quán)重取樣（competitive adaptive reweighted sampling, CARS)競(jìng)爭(zhēng)性自適應(yīng)加權(quán)取樣(CARS)是一種基于回歸系數(shù)進(jìn)行波長(zhǎng)點(diǎn)選擇的方法。該方法模仿達(dá)爾文進(jìn)化論中的“適者生存”原則，將每個(gè)波長(zhǎng)看作一個(gè)個(gè)體，對(duì)波長(zhǎng)實(shí)施逐步淘汰。利用回歸系數(shù)絕對(duì)值的大小作為衡量波長(zhǎng)重要性的指標(biāo)，同時(shí)，引入了指數(shù)衰減函數(shù)來(lái)控制波長(zhǎng)的保留率。每次通過(guò)自適應(yīng)重加權(quán)采樣技術(shù)篩選出 PLS 模型中回歸系數(shù)絕對(duì)值大的波長(zhǎng)點(diǎn)，去掉權(quán)重小的波長(zhǎng)點(diǎn)，利用交差驗(yàn)證選出模型交叉驗(yàn)證均方根誤差(root mean square error of crossvalidation, RMSECV)值最低的子集，可有效選擇與目標(biāo)值相關(guān)的最優(yōu)波長(zhǎng)組合。CARS-Monte-Carlo-sampling算法采用蒙特卡羅采樣法采樣N次，每次從樣品集中隨機(jī)抽取 80%的樣品作為校正集，利用抽取的光譜矩陣，和濃度矩陣，分別建立PLS模型，采用指數(shù)衰減函數(shù)EDF強(qiáng)行去除回歸系數(shù)值相對(duì)較小的波長(zhǎng)點(diǎn)，第i次采樣時(shí)，波長(zhǎng)點(diǎn)的保留率r=e-bt，第一采樣時(shí)，所有的m個(gè)變量都被用于建模，故第N次采樣時(shí)僅兩個(gè)波長(zhǎng)被使用，通過(guò)N 次采樣時(shí)篩選出 PLS 模型中回歸系數(shù)絕對(duì)值大的波長(zhǎng)點(diǎn)，用每次產(chǎn)生的新變量子集建立PLS回歸模型，計(jì)算每個(gè)模型的交互驗(yàn)證均方差 RMSECV，選擇 RMSECV 值最小的變量子集，即為最優(yōu)變量子集。

?4.4 建模方法

定量分析是對(duì)被研究對(duì)象所含成分的數(shù)量關(guān)系或所具備性質(zhì)件的數(shù)量關(guān)系進(jìn)行量化的分析過(guò)程。采用無(wú)損檢測(cè)技術(shù)獲取煙葉的化學(xué)成分信息時(shí)，通常只能獲取與待測(cè)樣品化學(xué)成分相關(guān)的間接信息（如光譜信息），如果要進(jìn)一步了解待測(cè)樣品的品質(zhì)信息（如煙堿、總糖、總氮等），則需要將無(wú)損檢測(cè)方法所獲取的待測(cè)樣品信號(hào)特征與常規(guī)方法（如流動(dòng)分析法或者其他化學(xué)方法）獲取的信息建立相應(yīng)的定量校正分析模型。

偏最小二乘（partial least squares, PLS）也是一種基于因子分析的多變量校正方法，在主成分回歸中，只對(duì)光譜矩陣X進(jìn)行分解，消去了自變量光譜陣X中的噪聲信息，而因變量濃度陣Y也會(huì)含有不相關(guān)信息，因此因變量濃度陣Y也應(yīng)同樣處理。與PCR不同的是，在PLS中，自變量光譜陣X和因變量濃度陣Y的分解應(yīng)該同時(shí)進(jìn)行，應(yīng)該是將光譜陣X信息引入到濃度陣Y的分解過(guò)程中，在每計(jì)算光譜陣一個(gè)新的主成分之前，交換光譜陣X與濃度陣Y的得分，從而使自變量主成分直接與被分析組分含量關(guān)聯(lián)。20世紀(jì)80年代開始PLS就應(yīng)用于化學(xué)分析研究，現(xiàn)在化學(xué)計(jì)量學(xué)中多變量校正方法中最受推崇的之一就是PLS，廣泛應(yīng)用在化學(xué)分析測(cè)量和相關(guān)的研究中。

4.5 模型評(píng)價(jià)

通常校正集和驗(yàn)證集中都會(huì)有相關(guān)系數(shù)（R）、校正集中有校正標(biāo)準(zhǔn)偏差（SEC/RMSEC）、驗(yàn)證集中有預(yù)測(cè)標(biāo)準(zhǔn)偏差（SEP/RMSEP）等。

（1）? 相關(guān)系數(shù)

兩個(gè)定量校正結(jié)果相關(guān)程度的一個(gè)統(tǒng)計(jì)量是用決定系數(shù)來(lái)描述的，它主要是用來(lái)判斷定量校正模型與待測(cè)組分的線性關(guān)系的好壞

（2）? 校正標(biāo)準(zhǔn)偏差和預(yù)測(cè)標(biāo)準(zhǔn)偏差

建立近紅外光譜定量校正模型時(shí)，通常是要將樣品分為校正集和預(yù)測(cè)集兩個(gè)部分，校正集樣本是用來(lái)建立定量校正模型，預(yù)測(cè)集是用來(lái)對(duì)模型進(jìn)行驗(yàn)證的，通過(guò)比較真實(shí)值與預(yù)測(cè)值的差異，來(lái)判斷模型預(yù)測(cè)能力的好壞。SEC和SEP是計(jì)算模型得出的真實(shí)值與其預(yù)測(cè)值之間的誤差平方和的均方根值。

在校正模型中，采用留一交互驗(yàn)證法來(lái)建立定量校正模型。留一交互（叉）驗(yàn)證法就是：每次從煙葉校正樣本集中取出一個(gè)煙葉樣本，然后用余下的煙葉樣本來(lái)建立校正模型，用建好的校正模型來(lái)預(yù)測(cè)之前取出的這個(gè)煙葉樣本，直到煙葉校正樣本集中每個(gè)煙葉樣本都被取出過(guò)一次。

對(duì)同一批次樣本,SEC和SEP值越小說(shuō)明模型的精度越高,兩者值越接近說(shuō)明模型穩(wěn)定性越好。

4.6?? 手持式模型分析

（1） 光譜預(yù)處理

通過(guò)對(duì)建模光譜的分析，由于手持式近紅外采取了直接緊貼煙葉表面的掃描方式，且外部環(huán)節(jié)保持相對(duì)穩(wěn)定，采集的光譜復(fù)雜程度并不高，只需要利用一介導(dǎo)數(shù)的預(yù)處理方法來(lái)提高模型的精度即可；

?圖表4 原始光譜??????圖表5 光譜一階導(dǎo)數(shù)

（2） 異常樣本剔除

采用蒙特卡洛交叉驗(yàn)證剔除化學(xué)值異常的樣本，對(duì)下圖比較分散的樣本進(jìn)行剔除，由于是縱坐標(biāo)比較分散的樣本點(diǎn)。

圖表6 煙堿異常統(tǒng)計(jì)??? ? 圖表7 氯異常統(tǒng)計(jì)?

圖表8 總糖異常統(tǒng)計(jì)??? ? ?圖表9 還原糖異常統(tǒng)計(jì)

圖表10總氮異常統(tǒng)計(jì)? ??????圖表11 鉀異常統(tǒng)計(jì)?

（3）波長(zhǎng)變量選擇

建模發(fā)現(xiàn)采用競(jìng)爭(zhēng)性自適應(yīng)權(quán)重取樣 (CARS)方法選擇的波長(zhǎng)建立的模型效果最佳，成分?jǐn)?shù)選擇13

（4）模型建立

采用偏最小二乘法建立的回歸模型效果最佳，下圖分別是6個(gè)指標(biāo)的PLS模型預(yù)測(cè)值與真實(shí)值對(duì)比圖和真實(shí)值絕對(duì)誤差圖。從圖中可以看出，所有模型都集中的回歸直線附近，誤差也近似正態(tài)分布，說(shuō)明建立的模型效果較好。

圖表12 煙堿模型

圖表13 總糖模型

圖表14 還原糖模型

圖表15 總氮模型

圖表16 鉀模型

圖表17 氯模型

?（5）模型評(píng)價(jià)

圖表18 模型評(píng)價(jià)結(jié)果

5?模型外部驗(yàn)證

煙堿指標(biāo)

圖表19 煙堿外部驗(yàn)證結(jié)果

圖表20 煙堿外部驗(yàn)證比對(duì)圖

總糖指標(biāo)

圖表21 總糖外部驗(yàn)證結(jié)果

圖表22 總糖外部驗(yàn)證比對(duì)圖

還原糖指標(biāo)

圖表23 還原糖外部驗(yàn)證結(jié)果

圖表24 還原糖外部驗(yàn)證比對(duì)圖

總氮指標(biāo)

圖表25 總氮外部驗(yàn)證結(jié)果

圖表26 總氮外部驗(yàn)證比對(duì)圖

鉀指標(biāo)

圖表27 鉀外部驗(yàn)證結(jié)果

圖表28 鉀外部驗(yàn)證比對(duì)圖

氯指標(biāo)

圖表29 氯外部驗(yàn)證結(jié)果

圖表30 氯外部驗(yàn)證比對(duì)圖

驗(yàn)證結(jié)論

圖表31 50個(gè)驗(yàn)證樣本驗(yàn)證結(jié)果

對(duì)上述表格數(shù)據(jù)進(jìn)行分析，手持式近紅外對(duì)原煙相關(guān)化學(xué)指標(biāo)的檢測(cè)結(jié)果基本都達(dá)到了實(shí)際使用的精度需求。模型內(nèi)部建模樣本化學(xué)值梯度基本能夠覆蓋常規(guī)檢測(cè)樣本的化學(xué)值范圍，模型對(duì)不同等級(jí)產(chǎn)地或者時(shí)間梯度的原煙具有有穩(wěn)定的預(yù)測(cè)能力，且對(duì)絕大部分樣本的預(yù)測(cè)誤差落在一定的范圍內(nèi)，整體誤差呈正態(tài)分布；

建模樣本分析

此次模型中包含350個(gè)建模樣本，主要包括云南昆明、曲靖、文山、麗江、臨滄、玉溪主要產(chǎn)區(qū)，建模樣本基本覆蓋了云南煙葉復(fù)烤有限責(zé)任公司石林復(fù)烤廠烤季原煙的所有產(chǎn)區(qū)。其中，上部煙樣本數(shù)量為79，中部煙數(shù)量為184，下部煙數(shù)量為87，建模樣本在部位分布上基本滿足模型建立的需求，模型對(duì)不同產(chǎn)地不同等級(jí)的原煙都具有較好的預(yù)測(cè)能力。

圖表32 產(chǎn)地分布圖

對(duì)建模樣本基礎(chǔ)數(shù)據(jù)分析，煙堿、總糖、還原糖、總氮等指標(biāo)的建模集化學(xué)值分布基本都呈正態(tài)分布，能夠滿足模型建立的要求。

圖表33建模集煙堿分布圖 ????????圖表34 建模集總糖分布圖

圖表35 建模集還原糖分布圖 ????圖表36 建模集總氮分布圖?

6?結(jié)論

綜上所述，基于350個(gè)樣品所建立的6種常規(guī)化學(xué)成分的預(yù)測(cè)模型，其相關(guān)系數(shù)及偏差均符合項(xiàng)目要求。其中，煙堿預(yù)測(cè)模型的決定系數(shù)高達(dá)0.96，平均絕對(duì)偏差值為0.22。建立的預(yù)測(cè)模型可應(yīng)用于煙葉卸車、入庫(kù)、挑選等不同環(huán)節(jié)的原煙化學(xué)信息檢測(cè)，其檢測(cè)數(shù)據(jù)可為原料煙葉的化學(xué)成分波動(dòng)、質(zhì)量分析等提供可靠的數(shù)據(jù)，進(jìn)行一次調(diào)控，為均質(zhì)化配方打葉提供數(shù)據(jù)支撐。

致謝：本研究得到了云南銘帆科技有限公司的技術(shù)支持，特此表示感謝。

參考文獻(xiàn)

[1]樂(lè)俊明,陳鷹,丁映,等.近紅外光譜分析法測(cè)定煙草化學(xué)成分[J].貴州農(nóng)業(yè)科學(xué)，2005, 33(3):62-63.

YUE Jun-ming, CHEN Ying, DING Ying，et al. Determination of Tobacco Chemical Components by Near Infrared Spectrum Method[J]. Guizhou Agricultural Sciences , 2005, 33(3):62-63.

[2]邱軍,張懷寶,宋巖,王允白,許武,付中會(huì),李乃會(huì),等. 近紅外光譜分析技術(shù)在煙草行業(yè)的應(yīng)用[J].中國(guó)煙草科學(xué),2008,29(1):55-59.

QIU Jun，ZHANG Huaibao1，SONG Yan，WANG Yunbai，XU Wu，FU Zhonghui，LI Naihui，et al. Application of Near Infrared Spectroscopic Analytical Techniques in Tobacco

Industry[J].China Tobacco Science,2008,29(1):55-59.

[3]張辭海，胡蕓，劉娜，彭黔榮，邵學(xué)廣，等. 烤煙煙堿近紅外定量模型的適用性[J].煙草科技,2019,52(1):53-58.

ZHANG Cihai1，HU Yun1，LIU Na1，PENG Qianrong，SHAO Xueguang，et al. Adaptability of near-infrared spectroscopy based quantitative model for nicotine determination in flue-cured tobacco[J].Tobacco Science & Technology,2019,52(1):53-58.

[4]孫陽(yáng),劉翠玲,孫曉榮,聞世震，等. 基于便攜式近紅外光譜儀的櫻桃番茄糖分快速分析模型研究[J]. 食品與發(fā)酵工業(yè)，2021.

SUN Yang，LIU Cuiling，SUN Xiaorong，WEN Shizhen，et al. Rapid analysis model of sugar content in cherry tomatoes based on portable near infrared spectrometer[J]. Food and Fermentation Industries，2021.

[5]王東丹，李天飛，吳玉萍，張承聰，葉蘭欽，陳潤(rùn)瓊，楊金輝，等. 近紅外光譜分析技術(shù)在煙草化學(xué)分析上的應(yīng)用研究[J].云南大學(xué)學(xué)報(bào)（自然科學(xué)版），2001，23（2）：135—137.

WANG Dongdan，LI Tianfei，WU Yuping，ZHANG Chengcong，YE Lanqin，CHEN Runqiong，YANG Jinhui，et al. Applied Study of Near-infrared Spectroscopy Analytical Technique

on Chemical Analysis of Tobacco[J] Journal of Yunnan University，2001，23（2）：135—137.

[6]閔順耕，李寧，張明祥，等.近紅外光譜分析中異常值的判別與定量模型優(yōu)化[J].光譜學(xué)與光譜分析，2004,24（10）:1205—1209.

MIN Shungeng，LI Ning，ZHANG Mingxiang，et al. Outlier Diagnosis and Calibration Model Optimization for Near InfraredSpectroscopy Analysis，2004,24（10）:1205—1209.

[7]江蘇，馬翔，陳永福，繆明明，李天飛，王亞明，等. 近紅外光譜分析技術(shù)及其在煙草行業(yè)中的應(yīng)用[J]. Chinese Journal of Spectroscopy Laboratory，2006,23（3）：633—637.

JIANG Su，MA Xiang，CHEN G YongFu，MIAO MingMing，LI TianFei，WANG YaMing[J]. Near Infrared Spectroscopic Technology and Its Applications in Tobacco Industry，2006,23（3）：633—637.