本发明涉及一种对清香型原酒品质分级的方法，属于白酒鉴别领域。
背景技术：
：：白酒作为中国传统饮料之一，有悠久的历史。白酒的生产属于固态发酵，发酵过程中很难对各种数据进行实时监测，因此目前仍未完全实现数字化生产。不同批次生产出来的白酒存在或多或少的差异，因此，酒厂需要对不同原酒进行品质品鉴分级。目前，对原酒的分级只能通过专业人员凭借多年的经验，根据感官对原酒的品质进行鉴定，该方法不仅效率较低，费用高，而且准确性受主观因素影响较大。因此，通过科学手段，客观、准确地区分不同批次原酒之间的差异，并进行准确分级，更加快速科学高效。主成分分析(principalcomponentsanalysis，以下简称pca)是最重要的降维方法之一。对数据进行降维的同时损失最小的信息，即用主要的特征来代替整体数据。原酒扫描出来的数据图谱包含的信息非常多，相对弱化了品质相关信息，会对品质分级起到干扰作用。线性判别分析(lineardiscriminantanalysis，以下简称lda)算法一种监督学习的降维方法，lda算法的思想是将数据投影到低维空间之后，使得同一类数据尽可能的紧凑，不同类的数据尽可能分散。这点和pca不同，pca是不考虑样本类别输出的无监督降维技术。也就是说它的数据集的每个样本是有类别输出的。利用这个方法处理后即可得到原酒品质的较好分级。技术实现要素：[技术问题]目前对原酒的分级都是通过专业人员根据经验、感官进行评定，效率低且准确率受主观影响大。[技术方案]为了解决上述问题，本发明采用三维荧光光谱分析不同等级的原酒酒样，获得各个酒样的三维荧光光谱数据；之后采用平行因子分析酒样光谱数据的差异性，提取4条光谱作为新的特征量，并采用主成分分析、线性判别降维，结合支持向量机建立白酒原酒分级模型，实现了原酒的分级。本发明的方法操作简单，无污染，而且准确率高。本发明的第一个目的是提供一种利用三维荧光光谱结合pca-lda-svm对清香型原酒品质分级的方法，包括如下步骤：(1)采用三维荧光光谱仪对待测原酒酒样进行扫描，得到三维荧光光谱光谱数据；(2)将步骤(1)得到的光谱数据进行矩阵缩小处理，得到53×61的光谱矩阵；(3)提取步骤(2)53×61的光谱矩阵的原始光谱中激发波长为290nm、310nm、340nm、360nm处的发射光谱数据，作为新的特征参数；(4)对步骤(3)中新的特征参数进行主成分分析，对数据进行降维，得到65个主成分的数据；(5)将步骤(4)得到的主成分的数据进行线性判别分析后，得到原酒的分级模型输入的特征量；(6)将步骤(5)得到的特征量代入原酒的分级模型中，得到待测原酒酒样的级别。在本发明的一种实施方式中，步骤(1)所述扫描的参数包括：激发与发射狭缝宽度为5nm，积分时间为1s，激发波长为200-500nm，步长5nm，发射波长为210-500nm。在本发明的一种实施方式中，步骤(1)扫描之后得到的三维荧光光谱是一个47641×3的光谱矩阵。在本发明的一种实施方式中，步骤(2)是将步骤(1)得到的47641×3的光谱原始数据中提取第3列的荧光强度值，共47641个数据，每隔781个数据作为新的一列，共61列得到781×61的光谱矩阵；之后每隔15个点取一个发射波长，得到一个53×61的光谱矩阵。在本发明的一种实施方式中，步骤(6)所述的原酒的分级模型的构建方式为：步骤1：使用三维荧光光谱仪采集某一品牌不同等级白酒的原酒酒样的三维荧光光谱光谱数据；步骤2：光谱数据前处理a、对所有原酒酒样进行三维荧光光谱扫描，每一个原酒酒样扫描后得到一个47641×3的光谱矩阵，提取第3列的荧光强度值，共47641个数据，每隔781个数据作为新的一列，共61列得到781×61的光谱矩阵；之后每隔15个点取一个发射波长，得到一个53×61的光谱矩阵；b、建立不同等级的清香型原酒的三维荧光光谱数据库；c、利用平行因子分析所有原酒酒样的光谱数据得到在激发波长为290nm、310nm、340nm、360nm处的光谱数据，作为新的特征参数；d、对新的特征参数进行主成分分析，对数据进行降维，得到65个主成分的数据；步骤3：将得到的65个主成分的数据作为新的特征参数进行线性判别分析后，采用线性支持向量机建模，得到原酒的分级模型。在本发明的一种实施方式中，步骤(6)所述的原酒的分级模型的构建方式中原酒酒样采集具体为：采集清香型大茬特级、优级和一级原酒；采集清香型二茬特级、优级和一级原酒，每个等级的酒收集65个原酒酒样样本，共390个原酒酒样样本。在本发明的一种实施方式中，步骤(6)所述的原酒的分级模型的构建方式中三维荧光光谱仪的扫描条件为：激发与发射狭缝宽度为5nm，积分时间为1s，激发波长为200-500nm，步长5nm；发射波长为210-500nm。在本发明的一种实施方式中，步骤(1)所述的三维荧光光谱仪为smartfluo-pro三维荧光光谱仪。本发明的第二个目的是本发明所述的方法在食品领域中的应用。在本发明的一种实施方式中，所述的应用是在清香型原酒分级中的应用。[有益效果](1)同一厂家生产的原酒的光谱有高度相似性，主成分提取的90％的信息可能是共性较多，无法体现其样本间的差异性，而本发明使用三维荧光光谱技术分析不同等级的原酒酒样，获得各个酒样的三维荧光光谱数据；采用平行因子分析酒样光谱数据的差异性，提取4条光谱作为新的特征量，并采用主成分分析、线性判别降维，结合支持向量机建立白酒原酒分级模型，之后进行未知原酒酒样的检测，大茬单独建精度在92.5％以上，预测未知酒样准确率为92.3％以上；大茬、二茬混合建模的精度为86.9％以上，预测未知酒样准确率在85.9％以上。(2)本发明是一种的新的白酒原酒质量控制方法，操作简单，无污染。附图说明图1为6个不同等级的清香型大、二茬三维荧光光谱图，其中(a)为大茬特级酒；(b)为大茬优级酒；(c)为大茬一级酒；(d)为二茬特级酒；(e)为二茬优级酒；(f)为二茬一级酒。图2为平行因子分析360个光谱数据得到的3个组分图；其中(a)为310nm/335nm荧光峰；(b)为290nm/310nm荧光峰；(c)为340nm/400nm荧光峰。图3为由激发波长290nm、310nm、340nm、360nm形成新的特征参数光谱图。图4为前6个主成分的成分累积贡献图。图5为采用线性支持向量机建模得到的大茬原酒的分级模型。图6为采用线性支持向量机建模得到大茬、二茬原酒的分级模型。图7为采用bp神经网络建模得到的大茬原酒的分级模型；其中(a)为训练集拟合效果图；(b)为验证集拟合效果图；(c)为预测集拟合效果图；(d)为总数据集拟合效果图。图8为采用bp神经网络建模得到大茬、二茬原酒的分级模型；其中(a)为训练集拟合效果图；(b)为验证集拟合效果图；(c)为预测集拟合效果图；(d)为总数据集拟合效果图。图9为一种利用三维荧光光谱结合pca-lda-svm对清香型原酒品质分级的方法的流程图。图10为大茬单独建模准确率的测试结果。图11为大茬、二茬混合建模的准确率的测试结果。具体实施方式以下对本发明的优选实施例进行说明，应当理解实施例是为了更好地解释本发明，不用于限制本发明。实施例中提及的平行因子分析方法可参考文献(c.m.andersen,r.bro.practicalaspectsofparafacmodelingoffluorescenceexcitation-emissiondata[j].journalofchemometrics,2003,17:200-215.)；主成分分析方法可参考文献(allegrettai,marangonib,manzarip,etal.macro-classificationofmeteoritesbyportableenergydispersivex-rayfluorescencespectroscopy(ped-xrf),principalcomponentanalysis(pca)andmachinelearningalgorithms[j].talanta,2020,212:120785.)；支持向量机建模方法可参考文献(guj,lus.aneffectiveintrusiondetectionapproachusingsvmwithnavebayesfeatureembedding[j].computers&security,2020,103(3):102158.)；bp神经网络建模方法可参考文献(zhun,wangk,zhangsl,etal.applicationofartificialneuralnetworkstopredictmultiplequalityofdry-curedhambasedonproteindegradation[j].foodchemistry,2020,344(12):128-586.)；线性判别分析可参考文献(武斌,王大智,嵇港,等.基于正交线性判别分析和电子鼻技术的食醋分类[j].食品与发酵工业,2020,046(006):263-268.)。实施例1酒的分级模型的构建-向量机建模酒的分级模型的构建方法，包括如下步骤：(1)使用smartfluo-pro三维荧光光谱仪采集清香型某一品牌不同等级的白酒原酒的三维荧光光谱光谱数据；a、酒样采集：采集清香型大茬特级、优级和一级原酒；采集二茬特级、优级和一级的原酒，每个等级的酒收集65个样本，共390个样本；b、扫描条件：激发与发射狭缝宽度为5nm，积分时间为1s，激发波长为200-500nm，步长5nm；发射波长为210-500nm；(2)光谱数据前处理a、对390个样本进行三维荧光光谱扫描，每一个样本扫描后得到一个47641×3的光谱矩阵，进行变换后得到781×61的光谱矩阵，由于发射波长步长太密导致数据量过大，可以去除一些发射波长的点进行适当缩放减小数据量仍然保持光谱的有效信息，每隔15个点取一个发射波长，得到一个53×61的光谱矩阵，去除由瑞丽散射和水拉曼形成的峰；b、建立不同等级的清香型原酒的三维荧光光谱数据库，如图1所示；c、利用平行因子分析390个样本的光谱数据的差异性，390个样本的光谱数据主要由3个组分图中的荧光峰组成，每个样本的荧光光谱中这3个荧光峰的含量不同，3个组分图如图2所示，共有3个峰值，分别是在激发波长为290nm、310nm、340nm处，再额外多选取一个360nm处的发射光谱，提取其3个组分的峰值激发波长处的发射光谱数据和激发波长为360nm处的发射光谱数据，共4条发射光谱数据作为新的特征参数，特征参数谱图如图3所示；d、对新的特征参数进行主成分分析，对数据进行降维，得到前65个主成分的数据，65个主成分的累计贡献率为99.92％；图4为前6个主成分的成分累积贡献图；(3)先进行大茬原酒单独分级：线性判别结合支持向量机建模利用大茬3个等级酒共195个酒样，提取4条光谱后的数据，经主成分分析得到65个主成分作为新的特征参数，线性判别分析后，采用支持向量机建模，核函数为高斯函数，最佳c值为987.78，最佳g值为474.74，建模精度为92.5％，得到原酒的分级模型(图5)。(4)大茬、二茬原酒混合建模：线性判别结合支持向量机建模利用大茬3个等级和二茬3个等级酒，共390个酒样，提取4条光谱后的数据，经主成分分析得到65个主成分作为新的特征参数，线性判别分析后，采用支持向量机建模，核函数为线性，最佳惩罚系数c为998.39，建模精度为86.9％，得到酒的分级模型，其结果如图6所示。对照例1bp神经网络建模将实施例1步骤(2)得到的前65个主成分的数据进行建模：(1)先进行大茬原酒单独分级：bp神经网络利用大茬3个等级酒共195个酒样，提取4条光谱后的数据，经主成分分析得到65个主成分作为新的特征参数，线性判别分析后，采用bp神经网络建模，得到原酒的分级模型(图7)；(2)大茬、二茬原酒混合建模：bp神经网络利用大茬3个等级和二茬3个等级酒，共195个酒样，提取4条光谱后的数据，经主成分分析得到65个主成分作为新的特征参数，线性判别分析后，采用bp神经网络建模，得到原酒的分级模型，如图8所示。对照例2调整实施例1步骤(2)c中的4条发射光谱数据为290nm和310nm两条发射光谱，其他和实施例1保持一致，得到大茬原酒单独的分级模型和大茬、二茬原酒的分级模型。实施例2调整实施例1中主成分的个数为6、25、65、75、全部数据，其他和实施例1保持一致，得到检测结果。结果如表1所示，从表1可以看出：6个主成分代表着原始数据的信息95.01％，但其准确率相比较于原始数据有所降低；当主成分个数为65时，代表原始数据的99.92％的信息，其准确率为最高，与原始数据相比提高2.3％；而当主成分个数继续增加时，准确度相较于65个时有所下降。可见，主成分累计贡献率为95％时的信息其相似度较高，差异性并不体现在其中，更多地体现在累计贡献率95.01％到99.92％的贡献率的成分中。因此建模过程中进行主成分降维时选取的个数为65个。表1主成分个数对建模准确率的影响注：表中的准确率是采用五折交叉验证法得到的。实施例3一种利用三维荧光光谱结合pca-lda-svm对清香型原酒品质分级的方法(如图9)，包括如下步骤：(1)取78个新产原酒作为待测酒样，采用三维荧光光谱仪对待测酒样进行扫描，得到三维荧光光谱光谱数据；(2)将步骤(1)得到的三维荧光光谱光谱数据中进行变换，得到781×61的光谱矩阵；之后每隔15个点取一个发射波长，得到一个53×61的光谱矩阵；(3)提取步骤(2)53×61的光谱矩阵的原始光谱中激发波长为290nm、310nm、340nm、360nm处的发射光谱数据，作为新的特征参数；(4)对步骤(3)中新的特征参数进行主成分分析，对数据进行降维，得到65个主成分的数据；(5)将步骤(4)得到的65个主成分的数据进行线性判别分析后，得到原酒的分级模型输入的特征量；(6)将步骤(5)得到的特征量输入实施例1得到的原酒的分级模型中，得到待测酒样的级别。具体如下：大茬单独建模预测准确率为92.3％(图10)；大茬、二茬混合建模的预测准确率为85.9％(图11)。对照例3将实施例3得到的65个主成分的数据输入对照例1的原酒的分级模型，得到的结果如下：大茬单独建模预测准确率为90％(表2)，大茬和二茬混合建模为83.3％(表3)。表2lda-bp神经网络大茬模型部分预测结果表3lda-bp神经网络大茬和二茬混合建模部分预测结果对照例4将实施例3得到的65个主成分的数据输入对照例2的原酒的分级模型，得到的结果如下：大茬单独建模预测准确率为80％，大茬和二茬混合建模为77.95％。综上所述，本发明通过扫描清香型大、二茬各3个等级原酒的三维荧光光谱，利用平行因子算法对光谱数据进行分析得到4个组分，其中共出现峰值点4个，分别在激发波长290nm、310nm、340nm和360nm处。从原有的光谱数据中提取出4个峰值波长处的发射光谱作为新的特征参数，进行主成分降维，发现各光谱样本之间的差异性主要集中在主成分贡献率的后10％中，而前90％的贡献率的主成分中进行建模误差率较大。提取99.92％贡献率的主成分维度作为建模特征产量，利用线性判别结合支持向量机算法建模，建模分为两部分，一是大茬建模，以195个大茬酒样的荧光光谱数据作为训练集建模，交叉验证后模型精度为92.5％，再采集39个未参与建模的大茬酒样的光谱数据作为测试集，验证模型，预测结果为92.3％％；二是大茬和二茬混合建模，以390个大、二茬酒样作为训练集建模，交叉验证后模型精度为86.9％，再采集39个未参与建模的大茬酒样的光谱数据作为测试集，验证模型，预测结果为85.9％，较好的实现了对清香型大、二茬各等级间的分类。虽然本发明已以较佳实施例公开如上，但其并非用以限定本发明，任何熟悉此技术的人，在不脱离本发明的精神和范围内，都可做各种的改动与修饰，因此本发明的保护范围应该以权利要求书所界定的为准。当前第1页12当前第1页12