摘要：為研究煙葉全面質(zhì)量評價的新方法，2014年于四川省寧南縣采集煙葉樣品，測定分析了烤煙的外觀質(zhì)量、物理特性、化學(xué)成分、中性致香物質(zhì)和感官質(zhì)量，并應(yīng)用主成分分析和聚類分析對煙葉38個質(zhì)量指標(biāo)進(jìn)行了分析和聚類。結(jié)果表明，38個指標(biāo)共提取出了7個主成分，方差累計貢獻(xiàn)率為87.37%，基本上反映了原來變量的全部信息；煙葉樣品經(jīng)過聚類分析共分成3類，經(jīng)過判別分析驗證，正確率達(dá)92.86%；最適宜種植烤煙的海拔高度為1 450 m，其次為海拔1 750 m，海拔2 050 m種植效果最差，黃壤種植效果優(yōu)于紅壤，云煙87和云煙85種植效果優(yōu)于云煙97。

關(guān)鍵詞：主成分分析；聚類分析；煙葉；全面質(zhì)量評價

中圖分類號：S572 文獻(xiàn)標(biāo)識碼：A 文章編號：0439-8114（2016）19-5163-03

DOI：10.14088/j.cnki.issn0439-8114.2016.19.070

The Application of the Principal Component Analysis and Cluster Analysis in Overall Quality Evaluation of Tobacco

L？譈 Peng-hui1，YE Xie-feng1，YUE Dong-lin2，F(xiàn)AN Fang-bing2，TIAN Sheng-guo2，LI Yan-jie2，

LU Jian3，ZOU Hui-yu3，ZHOU Han-hong4，LI Yan4，YU Jian-jun1

（1. College of Tobacco Science， Henan Agricultural University，Zhengzhou 450002， China； 2. Gansu Tobacco Industrial Co.， Ltd.， Lanzhou 730000， China； 3. Ningnan Branch Company， Liangshan Prefecture Company of Sichuan Provincial Tobacco Company， Ningnan 615400， Sichuan， China； 4. Wuhan Cigarette Factory， Wuhan 430000， China）

Abstracts： In order to study the new method of overall quality evaluation of tobacco， tobacco leaves were collected from Ningnan county of Sichuan province in 2014， and their appearance quality， physical properties，chemical components， neutral aroma components and sensory quality were determined and analyzed. 38 quality indexes of tobacco leaf were analyzed and clustered by principal component analysis and cluster analysis. The results showed that，seven principal components were extracted from 38 indexes， and cumulative contribution rate of variance was 87.37%，which basically reflected all information of original variable； tobacco leaves were divided into three categories by cluster analysis，and correct rate reached 92.86% by discriminant analysis. The optimum altitute for planting flue-cured tobacco was 1 450 m，followed by the altitute 1 750 m， and 2 050 m in turn， the planting effects of yellow soil was better than that of red soil，and the planting effects of Yunyan 87 and Yunyan 85 were better than that of Yunyan 97.

Key words： principal component analysis； cluster analysis； tobacco leaf； overall quality evaluation

煙草是中國重要的經(jīng)濟(jì)作物之一，種植面積和總產(chǎn)量居世界第一位，在中國國民經(jīng)濟(jì)中占有重要地位[1]。然而煙草對環(huán)境條件的變化十分敏感，環(huán)境條件的差異不僅影響煙草的形態(tài)特征和農(nóng)藝性狀，而且還直接影響煙葉的化學(xué)成分和質(zhì)量[2]?？緹熎贩N、土壤類型和海拔因素是導(dǎo)致烤煙內(nèi)含物存在地域差異和產(chǎn)量高低的重要生態(tài)因素，對煙葉的產(chǎn)量和產(chǎn)值、化學(xué)成分、致香物質(zhì)含量和感官質(zhì)量有很大的影響[3-5]。當(dāng)前對煙葉進(jìn)行質(zhì)量評價的方法主要有外觀質(zhì)量評價法、常規(guī)化學(xué)成分評價法、感官質(zhì)量評價法等[6-8]。而外觀評價和感官質(zhì)量評價方法受主觀影響較大，化學(xué)成分分析法不能完整如實地反映煙葉的質(zhì)量，因此構(gòu)建一種評價煙葉全部質(zhì)量指標(biāo)的模型顯得尤為重要。

研究人員對于煙葉質(zhì)量評價方法進(jìn)行了很多研究[9-12]，已有研究雖然利用了主成分分析、聚類分析和判別分析開辟了煙葉評價的新途徑，但是評價的指標(biāo)不夠全面。為此，本研究通過主成分分析、聚類分析和判別分析對煙葉質(zhì)量進(jìn)行全面評價，以期探討海拔、烤煙品種和土壤類型對煙葉質(zhì)量的影響，為煙葉生產(chǎn)的栽培管理和品種規(guī)劃提供參考。

1 材料與方法

1.1 樣品采集

試驗于2014年在四川省寧南縣取樣。取樣時記錄每個取樣點的烤煙品種、海拔高度和土壤類型，共計14個樣品。每個樣品設(shè)3次重復(fù)，均各取2 kg C3F等級烤后煙葉備用。試驗所取樣品信息具體見表1。

1.2 指標(biāo)測定及方法

烤煙外觀質(zhì)量的評價按照GB 2635-1992[13]進(jìn)行，由3名烤煙分級專家對外觀質(zhì)量進(jìn)行評價。

各物理特性指標(biāo)測定方法參考吉書文等[14]的方法進(jìn)行。

樣品中常規(guī)化學(xué)成分總糖、還原糖、煙堿、總氮、氯含量測定參考王瑞新[15]的方法進(jìn)行，鉀含量測定按照YC/T 173-2003[16]進(jìn)行。

中性致香物質(zhì)由國家煙草栽培生理生化研究基地利用GC/MC采用內(nèi)標(biāo)法測定分析。

煙葉感官質(zhì)量由河南中煙工業(yè)有限責(zé)任公司技術(shù)中心和國家煙草栽培生理生化研究基地組織評吸委員會進(jìn)行評定。評分后由有利于吸食品質(zhì)的指標(biāo)得分減去不利于吸食品質(zhì)的指標(biāo)得分，即得到感官質(zhì)量的總得分。

1.3 主成分分析

把不同樣品的各個質(zhì)量指標(biāo)構(gòu)成主成分分析相關(guān)矩陣，根據(jù)主成分分析后所確定的不同主分量的線性組合計算不同樣品每個主成分的得分，進(jìn)而計算每個煙葉樣品各個主成分的綜合得分。主成分分析采用SPSS 19.0軟件進(jìn)行。具體步驟為：將n個樣本的p個指標(biāo)變量構(gòu)建成數(shù)據(jù)集，通過變換將原變量xi 轉(zhuǎn)換成因子變量，因子變量是原變量的線性組合，記原變量指標(biāo)為X1，X2，…，Xp，設(shè)它們降維處理后的綜合指標(biāo)，即新變量為f1，f2，f3，…，fm（m≤p），這m個變量之間互不相關(guān)，因子變量分別稱作原始變量的第1，第2，第3，…，第m個主成分，表示為：

f1=l11X1+l12X2+…+l1pXpf2=l21X1+l22X2+…+l2pXp …………fm=lm1X1+lm2X2+…+lmpXp

1.4 聚類分析

對用主成分分析提取出的主因子（fm）根據(jù)其組成重新命名，采用Q型聚類中類間評價鏈鎖法和歐氏距離平方對各個樣品的各主成分得分（Zm）進(jìn)行聚類。聚類分析采用SPSS 19.0軟件進(jìn)行。

1.5 判別分析

以聚類分析所確定的類別為基礎(chǔ)，建立以各主成分得分（Zm）為基礎(chǔ)的判別函數(shù)，將每一個主成分得分值（Zm）分別代入Fisher判別函數(shù)dik（為第k個判別函數(shù)所求得的第i個樣品的值，k為判別函數(shù)組中判別函數(shù)的個數(shù)），以函數(shù)dik的大小作比較，函數(shù)值最大值表明該主成分得分值（Zm）屬于該組，也即該樣品屬于該組[10]。判別分析采用SPSS 19.0軟件進(jìn)行。

2 結(jié)果與分析

2.1 煙葉質(zhì)量的主成分分析

2.1.1 煙葉質(zhì)量各指標(biāo)主成分分析相關(guān)矩陣的特征值 對14個煙葉樣品的外觀質(zhì)量、物理特性、化學(xué)成分、中性致香物質(zhì)和感官質(zhì)量共38個指標(biāo)進(jìn)行主成分分析，所得相關(guān)矩陣的特征值見表2。由表2可知，根據(jù)特征值大于1的原則，煙葉質(zhì)量的38個指標(biāo)共提取出了7個主成分，方差累計貢獻(xiàn)率為87.37%，大于85%，說明提取的7個主成分基本上反映了原來變量的全部信息。

2.1.2 煙葉質(zhì)量各指標(biāo)主成分分析因子載荷和因子得分 對14個煙葉樣品的外觀質(zhì)量、物理特性、化學(xué)成分、中性致香物質(zhì)和感官質(zhì)量共38個指標(biāo)進(jìn)行主成分分析，所得因子載荷矩陣見表3。

由表3可知，根據(jù)煙葉質(zhì)量各指標(biāo)在主成分中的載荷量，取載荷絕對值大于0.5可以看出，身份、外觀總得分、單葉重、含梗率、填充值、拉力、總糖、還原糖、總氮、苯丙氨酸類、新植二烯、香氣物質(zhì)總量、香氣質(zhì)、透發(fā)性、細(xì)膩程度、柔和程度、圓潤感、余味、刺激性、干燥感和評吸總得分在主成分f1中有較高的載荷；單葉重、含梗率、拉力、總糖、氯、煙堿、香氣質(zhì)、圓潤感、余味和干燥感在主成分f2中有較高的載荷；顏色、成熟度、身份、油分、色度、外觀總得分、苯丙氨酸類和棕色化產(chǎn)物類在主成分f3中有較高的載荷；葉寬、類西柏烷類、煙氣濃度、勁頭、香氣量和透發(fā)性在主成分f4中有較高的載荷；葉長、氯、鉀、類胡蘿卜素類和新植二烯主成分f5中有較高的載荷；雜氣在主成分f6中有較高的載荷；結(jié)構(gòu)在主成分f7中有較高的載荷。

2.2 基于主成分分析的煙葉質(zhì)量的聚類分析

基于表3的因子載荷矩陣，利用“材料與方法”中的公式求得各主成分的得分，然后對14個煙葉樣品采用Q型聚類中類間評價鏈鎖法和歐氏距離平方對各個樣品的各主成分得分（Zm）進(jìn)行聚類，結(jié)果共分成3類，分別定義為好、中、差，分別用1、2、3表示，結(jié)果見表4。

由表4可知，煙葉樣品共分成3類。第1類的煙葉樣品有4個，分別為T1、T2、T5、T6；第2類的煙葉樣品有9個，分別為T3、T4、T7、T8、T9、T10、T11、T12、T14；第3類的煙葉樣品有1個，為T13。

2.3 基于主成分分析和聚類分析進(jìn)行的不同煙葉樣品的判別分析

根據(jù)各個煙葉樣品的7個主成分得分和表4的聚類結(jié)果，對14個煙葉樣品進(jìn)行判別分析，以確定聚類分析結(jié)果的可靠性，判別分析結(jié)果見表5。由表5可知，經(jīng)過判別分析后，第1類的煙葉樣品有4個，第2類的煙葉樣品有10個，第3類的沒有煙葉樣品，共有13個樣品被正確分類，第2類有一個樣品分類錯誤。同聚類分析的結(jié)果比較發(fā)現(xiàn)，聚類分析的結(jié)果中92.86%的煙葉樣品分類是正確的。

3 小結(jié)

根據(jù)特征值大于1的原則，煙葉質(zhì)量的38個指標(biāo)共提取出了7個主成分，方差累計貢獻(xiàn)率為87.37%，大于85%，說明提取出的7個主成分基本反映了原來變量的全部信息，有效剔除了重復(fù)信息。

煙葉樣品經(jīng)過聚類分析共分成3類，第1類的煙葉樣品有4個，第2類的煙葉樣品有9個，第3類的煙葉樣品有1個。判別分析結(jié)果顯示，聚類分析結(jié)果共有13個樣品被正確分類，正確率達(dá)到92.86%。表明聚類分析和判別函數(shù)分析方法有較好的一致性，所建立的煙葉評價模型是可行的。

聚類分析結(jié)果表明，最適宜種植烤煙的海拔高度為1 450 m，其次為海拔1 750 m，海拔2 050 m的種植效果最差，黃壤種植效果優(yōu)于紅壤，云煙87和云煙85種植效果優(yōu)于云煙97。

參考文獻(xiàn)：

[1] 劉國順，王彥亭，陳海如，等.煙草栽培學(xué)[M].北京：中國農(nóng)業(yè)出版社，2003.

[2] 劉國順.國內(nèi)外煙葉質(zhì)量差距分析和提高煙葉質(zhì)量技術(shù)途徑探討[J].中國煙草學(xué)報，2003（增刊）：54-58.

[3] 程呂新，盧秀萍，許白成，等.基因型和生態(tài)因素對煙草香氣物質(zhì)含量的影響[J].中國農(nóng)學(xué)通報，2005，21（11）：137-139.

[4] 許自成，劉國順，劉金海，等.銅山煙區(qū)生態(tài)因素和煙葉質(zhì)量特點[J].生態(tài)學(xué)報，2005，25（7）：1748-1753.

[5] 邵 麗，晉 艷，楊宇虹，等.生態(tài)條件對不同烤煙品種煙葉產(chǎn)質(zhì)量的影響[J].煙草科技，2002（10）：40-45.

[6] MOSELEY J M. The relationshlip of maturity of the leaf at harvest and certain properties of the cured leaf of flue-cured tobacco[J].Tobacco Science，1963（7）：67-75.

[7] 彭黔榮，蔡元青，王東山，等.根據(jù)常規(guī)化學(xué)指標(biāo)識別煙葉品質(zhì)的BP神經(jīng)網(wǎng)絡(luò)模型[J].中國煙草學(xué)報，2005，11（5）：19-24.

[8] 杜詠梅，郭承芳，張懷寶.水溶性糖、煙堿、總氮含量與烤煙吃味品質(zhì)的關(guān)系研究[J].中國煙草科學(xué)，2000，21（1）：7-10.

[9] 常壽榮，吳 濤，羅華元，等.烤煙品種、部位及生態(tài)環(huán)境對煙葉致香物質(zhì)的影響[J].云南農(nóng)業(yè)大學(xué)學(xué)報，2010，25（1）：58-62.

[10] 葉協(xié)峰，魏躍偉，楊宇熙，等.基于主成分分析和聚類分析的烤煙質(zhì)量評價模型構(gòu)建[J].農(nóng)業(yè)科學(xué)與綜合研究，2009，25（3）：268-271.

[11] 鄒淵淵，伍岳慶，姚 宇.基于主成分和和回歸分析的煙葉質(zhì)量分析模型[J].計算機應(yīng)用，2011，31（S2）：81-89.

[12] 郭建華，宋紀(jì)真，王廣山，等.基于主成分分析和聚類分析的煙葉物理特性區(qū)域歸類[J].煙草科技，2014（8）：14-17.

[13] GB 2635-1992，烤煙國家標(biāo)準(zhǔn)[S].

[14] 吉書文，騰兆波.煙草物理檢測[M].鄭州：河南科學(xué)出版社，1997.

[15] 王瑞新.煙草化學(xué)[M].北京：中國農(nóng)業(yè)出版社，2003.

[16] YC/T 173-2003，煙草及煙草制品鉀的測定火焰光度法[S].