[发明专利]一种应用近红外光谱技术快速预测远安黄茶闷黄时间的方法

权利要求书

1.一种应用近红外光谱技术快速预测远安黄茶闷黄时间的方法，其特征在于，扫描获得不同闷黄时间鲜叶样品的近红外光谱，对样品光谱进行预处理剔除噪声信息后，样品光谱转化为成对的数据点；再将全部光谱数据均划分为20个子区间，分别建立每个子区间数据的最小二乘支持向量机方法模型，筛选出建模的最佳子区间数据；再应用主成分分析法对最佳子区间信息数据点进行压缩抽提，以主成分得分为输入数值，通过调节隐含层神经元个数和信息传递函数，建立不同闷黄时间样品的人工神经网络预测模型，用于预测黄茶的闷黄时间，具体包括以下步骤：

步骤一、鲜叶样品采集与分类

采集湖北省远安县安吉白茶品种单芽、一芽一叶和一芽二叶三个不同部位鲜叶样品，经杀青后对样品进行闷黄，同时精准记录闷黄时间；根据闷黄时间不同，将样品分为校正集和验证集2个集合，分别用于建立校正集近红外光谱预测模型和对校正集预测模型稳健性进行检验；

其中，一芽一叶由单芽、第一叶和较长梗构成，一芽二叶由单芽、第一叶、第二叶和长梗构成；

步骤二、光谱扫描

应用近红外光谱仪扫描获得全部闷黄样品的近红外光谱，光谱扫描范围4000-10000cm^-1、分辨率8cm^-1、检测器为InGaAs，每个样品采集3次光谱、每次扫描64次，然后对3次采集的光谱进行平均，以平均光谱作为该样品的最终光谱用于后续建立模型；

步骤三、光谱噪声信息预处理

应用化学计量学软件对步骤二中得到的近红外光谱采用矢量归一化方法进行去噪声预处理，提高光谱的信噪比，有利于建立稳健的预测模型；光谱去噪声后，再将样品光谱转化为成对的数据点；

步骤四、光谱子区间划分

将全部光谱信息数据均等分为20个光谱信息子区间，精准筛选反映样品闷黄时间的光谱信息子区间，用于后续建立最小二乘支持向量机模型；

步骤五、最小二乘支持向量机LS-SVM模型建立

本发明应用LS-SVM方法分别建立每个光谱信息子区间的预测模型，比较模型相关系数Rc和交互验证均方根方差RMSECV大小，初步筛选出建模的最佳光谱信息子区间数据，达到精准筛选反映闷黄时间的光谱信息的目的，该步骤中筛选出建模的最佳光谱信息子区间数据为6406.4-6703.4cm^-1；其中，Rc最大、RMSECV最小，表示建立的最小二乘支持向量机模型结果最佳，

其中，RMSECV计算公式为：

Rc计算公式为：

式中，n表示样本数，yi和y_i’分别为样品集中第i个样品的实测值和预测值，为样品集中第i个样品的实测值的平均值，式中i≤n；

步骤六、unsupervised Kohonen结构人工神经网络预测模型建立

应用非线性的人工神经网络方法进一步精准预测样品的闷黄时间，包括：

1)最佳光谱信息子区间主成分分析

采用主成分分析方法PCA对筛选的最佳光谱信息子区间数据进行主成分分析，得到每个主成分的单独贡献率值、累计贡献率值和主成分得分；在主成分累计贡献率≥85％才可以有效代表样品建模光谱子区间信息；该步骤中在采用主成分分析方法PCA对筛选的最佳光谱信息子区间数据进行主成分分析时，采用前3个主成分代表最佳光谱子区间的光谱信息；

2)人工神经网络预测模型建立

以最佳光谱信息子区间的主成分得分为输入值、以样品不同闷黄时间为输出值，应用Neuro Shell 2软件建立unsupervised Kohonen结构的人工神经网络预测模型，unsupervised Kohonen结构人工神经网络含有1个隐含层、3、4、5个神经元数和logistic、linear[0,1]和tanh三种信息传递函数，比较模型相关系数Rc和交互验证均方根方差RMSECV大小，得到最佳近红外光谱预测模型，其中，Rc越大、RMSECV越小，表示模型预测效果越好；该步骤中选择4个神经元和logistic信息传递函数建立不同闷黄时间样品的unsupervised Kohonen结构人工神经网络模型；

其中RMSECV计算公式为：

Rc计算公式为：

式中，n表示样本数，yi和y_i’分别为样品集中第i个样品的实测值和预测值，为样品集中第i个样品的实测值的平均值，式中i≤n；

其中以相关系数Rc最大和交互验证均方根方差RMSECV最小的模型为最佳模型，经比较后得到最佳校正集模型；

步骤七、模型稳健性检验

应用全部验证集样品对不同闷黄时间的人工神经网络预测模型效果进行检验，所得结果用相关系数Rp和验证均方差RMSEP表示，其中Rp越大、RMSEP越小则表示模型稳健性越好，可以准确的预测样品的闷黄时间；

其中RMSEP计算公式为：

Rp计算公式为：

式中，n表示样本数，yi和y_i’分别为样品集中第i个样品的实测值和预测值，式中i≤n。

2.根据权利要求1所述的一种应用近红外光谱技术快速预测远安黄茶闷黄时间的方法，其特征在于：所述的步骤一中鲜叶样品数量为120份，鲜叶样品按照3:1的比例划分为校正集和验证集，其中校正集样品90个、验证集样品30个。