摘要:运用MATLAB 5.0分析了云南地区3个烤烟品种的主要化学成分。结果表明,在0.05显著性水平下,云系列与K326品种间有显著差异的化学指标为果胶、木质素、水溶性灰分碱度、石油醚提取物;云系列和红大品种间有显著差异的化学指标为挥发酸、细胞壁物质含量、氨基氮、多酚;红大和K326品种间有显著差异的化学指标为氨基氮、多酚。
关键词:烤烟;化学成分;MATLAB 5.0;拟合分布;统计分析
中图分类号:S572;TS41+1文献标识码:A文章编号:0439-8114(2011)01-0165-04
Application of MATLAB in Discrepancy Analysis of chemical components of Different Flue-cured Cultivars
WANG Hao-ya,ZHANG Qiang,DONG Gao-feng,WANG Li-min,HE Zhi-jun,XIANG Ming
(Yunnan Academy of Tobacco Science, Kunming 650106,China)
Abstract: The principal chemical components of 3 flue-cured tobacco cultivars in Yunnan were analyzed by MATLAB 5.0. The results indicated that at the 0.05 significance level, the chemical indicators of significant differences between Yun series and K326 were volatile acids, fruit glue and polyphenol; between Yun series and Hongda were volatile acids, cell wall substances, amino nitrogen and polyphenol; and between Hongda and K326 are volatile acids,amino nitrogen and polyphenol.
Key words: flue-cured tobacco; chemical components; MATLAB 5.0; fitting distribution; statistical analysis
烟叶的化学成分是烟叶质量形成的物质基础,是决定烟叶品质的内在要素,也是卷烟配方应该考虑的重要指标。烟草种植的生态环境、栽培措施、采收调制、发酵加工以及陈化都会对作为卷烟原料的烟叶化学成分产生影响。研究烤烟的化学成分,对卷烟原料的品质具有非常重要的意义[1]。而在以往的研究中,大多数的试验数据处理用的是DPS、SAS、SPSS等统计软件中的菜单式统计分析方法,没有针对每个成分的数据做拟合分布,没有看各组成分适合哪种统计方法,直接进行方差分析等分析方法。而MATLAB是开放式的软件,可以根据数据的不同特点,自行编写程序,可以先进行拟合分布,然后根据分布特点选择适合的统计方法。
MATLAB(Matrix Laboratory)是由美国Math Works公司推出的一款集数值分析、矩阵计算、信号处理和图形显示于一体的工具软件。
自1984年推向市场以来,现已成为国际认可(EEE)的最优化的科技应用软件之一[2],广泛应用于矩阵代数、数值计算、数字信号处理、振动理论、神经网络控制、动态仿真等领域。由于目前的多数试验中,测定的物质主要是烤烟的总糖、还原糖、总氮、烟碱、氯、淀粉、钾、pH值、蛋白质,所以分析测定中避开这些常测定的物质,而测定烤烟的挥发酸、挥发碱、细胞壁物质含量、果胶、全纤维素、木质素、总灰分、水溶性灰分碱度、乙醇提取物、氨基氮、二氯甲烷提取物、石油醚提取物、多酚,并将MATLAB引入烟草领域,初步将其应用在不同烤烟品种部分化学成分和烟气成分差异性分析中,为后续进一步挖掘MATLAB在烟草领域中更强大的功用打下一定的基础。
1材料与方法
1.1试验材料
试验于2006年进行,地点设在云南烟草科学研究院。供试品种为K326、云烟85和云烟87(简称云系列)、红花大金元(简称红大)。这3个系列品种来自云南红河、玉溪地区的主要烟叶收购站。每个品种取C2F等级,每一样品取样质量至少为5kg。K326样品数139份,云系列样品数361份,红大样品数71份。
1.2试验方法
烟叶化学成分含量测定方法:挥发酸、挥发碱按照中华人民共和国烟草行业标准YC/T 35-1996[3]进行;细胞壁物质含量、果胶、全纤维素、木质素采用文献[4,5]方法测定;总灰分、水溶性灰分碱度、乙醇提取物、氨基氮、二氯甲烷提取物采用文献[6]方法测定;石油醚提取物、多酚分别按照中华人民共和国烟草行业标准YC/T 176-2003[7]、YC/T 202-2006[8]方法测定。
1.3统计分析方法
所有数据采用Excel 2003和MATLAB 5.0统计软件进行统计分析。
2结果与分析
用MATLAB 5.0对所有数据先进行拟合分布的分析,对服从正态分布和非正态分布的烟叶的品质特性,采用的方差分析方法是不同的,因此,数据拟合分布是方差分析的基础。
2.1不同烤烟品种挥发酸含量的差异显著性分析
用MATLAB 5.0编程,做关于3个烤烟品种挥发酸的方差分析,可以得到两个图,一个是盒子图(图1)、一个是图形(图2)。图1为云系列、K326、红大的烟叶挥发酸含量的盒子图,可以直观地看出3个品种挥发酸含量的取值范围和均值,K326挥发酸均值最大,云系列次之,红大最小。方差分析还输出一个图形窗口(图2),3条线比较直观地表示出3个品种挥发酸的平均值。根据拟合分布的结果可知挥发酸含量不服从正态分布,所以,根据图2可以看出,云系列与K326挥发酸含量无显著差异,红大与云系列挥发酸含量有显著差异。
2.2不同烤烟品种挥发碱含量的差异显著性分析
图3显示了云系列、K326、红大的烟叶挥发碱的盒子图,可以直观地看出3个品种挥发碱含量的取值范围和均值,红大挥发碱均值最大,K326次之,云系列最小。图4中3条线比较直观地表示出3个品种挥发碱的平均值。根据拟合分布的结果可知挥发碱含量不服从正态分布,所以,根据图4可以看出,3个品种之间挥发碱含量无显著差异。
2.3不同烤烟品种细胞壁物质含量的差异显著性分析
图5显示了云系列、K326、红大的烟叶细胞壁物质含量的盒子图,可以直观地看出3个品种细胞壁物质含量的取值范围和均值,云系列细胞壁物质含量均值最大,K326次之,红大最小。图6中3条线比较直观地表示出3个品种细胞壁物质含量的平均值。根据拟合分布的结果可知细胞壁物质含量服从正态分布,所以,根据图6可以看出,云系列与红大细胞壁物质含量有显著差异,K326与云系列及红大的无显著差异。
2.4不同烤烟品种果胶含量的差异显著性分析
图7显示了云系列、K326、红大的烟叶果胶含量的盒子图,可以直观地看出3个品种果胶的取值范围和均值,K326果胶均值最大,红大次之,云系列最小。方差分析还输出一个图形窗口(图8),3条线比较直观地表示出3个品种果胶含量的平均值。根据拟合分布的结果可知,果胶含量不服从正态分布,所以,根据图8可以看出,云系列与K326果胶含量有显著差异,红大与其余两个品种果胶含量无显著差异。
2.5不同烤烟品种全纤维素含量的差异显著性分析
图9显示了云系列、K326、红大的烟叶全纤维素含量的盒子图,可以直观地看出3个品种全纤维素的取值范围和均值,云系列全纤维素均值最大,K326次之,红大最小。方差分析还输出一个图形窗口(图10),3条线比较直观地表示出3个品种全纤维素的平均值。根据拟合分布的结果可知,全纤维素含量不服从正态分布,所以,根据图10可以看出,3个品种全纤维素含量之间无显著差异。
2.6不同烤烟品种木质素含量的差异显著性分析
图11显示了云系列、K326、红大的烟叶木质素含量的盒子图,可以直观地看出3个品种木质素的取值范围和均值,红大木质素均值最大,K326次之,云系列最小。方差分析还输出一个图形窗口(图12),3条线比较直观地表示出3个品种木质素含量的平均值。根据拟合分布的结果可知木质素含量不服从正态分布,所以,根据图12可以看出,红大与其余两个品种木质素含量无显著差异,云系列与其K326木质素含量有显著差异。
2.7不同烤烟品种其他化学成分的差异显著性分析
根据上述的方法,以此类推,总灰分含量不服从正态分布,3个品种总灰分含量之间无显著差异;水溶性灰分碱度不服从正态分布,云系列与K326水溶性灰分碱度有显著差异,红大与其余两个品种水溶性灰分碱度无显著差异;乙醇提取物服从正态分布,3个品种之间乙醇提取物无显著差异;氨基氮含量服从正态分布,云系列与K326氨基氮含量无显著差异,红大与其余两个品种氨基氮含量有显著差异;二氯甲烷提取物不服从正态分布,3个品种之间二氯甲烷提取物无显著差异;石油醚提取物不服从正态分布,云系列与K326石油醚提取物有显著差异,红大与其余两个品种石油醚提取物无显著差异;多酚含量不服从正态分布,云系列与K326多酚含量无显著差异,红大与其余两个品种多酚含量有显著差异。
3讨论
在品种间差异的研究中,在0.05显著性水平下得到以下结论:云系列与K326在品种间有显著差异的化学指标为果胶、木质素、水溶性灰分碱度、石油醚提取物;云系列和红大品种间有显著差异的化学指标为挥发酸、细胞壁物质含量、氨基氮、多酚;红大和K326品种间有显著差异的化学指标为氨基氮、多酚。
一般的关于烟叶主要化学成分研究中,都默认烟叶的主要化学数据符合正态分布,而没有真正的去验证试验所涉及的数据是否符合正态分布,有些情况下,检测指标的数据需要进行一些常规的转换(如正弦、余弦等),从而使其符合正态分布,然后再进行以后的统计分析。试验从云南烟叶资源的化学成分入手,试图通过MATLAB 5.0的分布拟合工具(Distribution Fitting Tool),对云南省不同烤烟品种的烟叶主要化学成分进行差异性分析。更加准确地揭示不同烤烟品种在主要化学成分方面的差异,为以后MATLAB在烟草领域的应用打下一定的基础。
参考文献:
[1] 周翔,梁洪波,董建新, 等.山东烟区烤烟化学成分含量变化及聚类分析[J].中国烟草科学,2009,30(6):13-17.
[2] 夏文芳,吴鸿修,李享元.基于USB的实时数据采集系统及其在MATLAB中的应用[J].电子技术应用, 2003, 29(2):14-17.
[3] YC/T35-1996, 烟草及烟草制品总挥发酸、碱的测定[S].
[4] 张槐苓,葛翠英,穆怀静,等. 烟草分析与检验[M]. 郑州:河南科学技术出版社,1994.
[5] 李兴波,闫克玉,丁海燕,等. 河南烤烟(40级)细胞壁物质含量及其规律性研究[J]. 郑州轻工业学院学报,1999,14(3):27-30.
[6] 肖协忠.烟草化学[M].北京:中国农业出版社,1997.
[7] YC/T 176-2003,烟草及烟草制品石油醚提取物的测定[S].
[8] YC/T 202-2006,烟草及烟草制品多酚的测定[S].
推荐访问: 烤烟 差异性 化学成分 品种 分析