(一)技术领域
[0001] 本
发明涉及基于多种检验技术的武夷岩茶产地智能鉴别方法，即联合
近红外光谱、稳定同位素、微量元素、
氨基酸、
儿茶素以及
电子舌数据鉴别武夷岩茶产地的方法，属于地理标志产品真实性识别技术领域。(二)背景技术
[0002] 根据GB/T 17924-2008的定义，地理标志产品是指利用产自特定地域的原材料，按照传统工艺在特定地域内所生产的，
质量、特色或者声誉在本质上取决于其原产地域地理特征，并按法定程序经审核批准以原产地域名称命名的产品。茶叶是典型的地理标志保护产品，有武夷岩茶、安溪
铁观音、大红袍、永春佛手、西湖龙井、安吉白茶、祁
门红茶、普洱茶、碧螺春茶等近50种地理标志产品茶叶。
[0003] 目前，国内外已开展成品茶产地识别鉴定研究，仪器检测结合化学计量学分析方法为最主要产地识别方法，仪器检测方法主要有
近红外光谱、同位素质谱、液相色谱、
传感器等；常用计量学方法包括偏最小二乘、主成分分析、
人工神经网络、
支持向量机等。
[0004] 近红外光谱检测技术广泛应用茶叶原产地检测中，Zhou等利用近红外光谱仪对25个西湖龙井和70个浙江龙井样品进行检测并建立费氏判别函数模型，训练集合、交叉验证集合和测试集合的识别准确率分别96.7％、95.3％和96.7％。周健等对4个龙井茶品种(龙井43、群体种、迎霜和乌
牛早)进行近红外检测，并采用PLS建立模型，其4个品种茶叶准确率分别为89.8％、90.9％、96.1％和99.5％。占茉莉等采用近红外光谱仪扫描10份西湖龙井和18份浙江龙井茶叶样品，
聚类分析表明，西湖龙井茶具有独特特征而自成一类，缙
云和新昌表现出类似光谱特征存在交叉，富阳产的也自成一类。赵杰文等采用近红外光谱仪对龙井、碧螺春、毛峰和铁观音(各20份)进行检测，并采用主成分建立
预测模型，校正集和预测集的鉴别准确率分别为98.75％和95.0％，但龙井样品易误判为碧螺春。陈全胜等采用近红外光谱对龙井、碧螺春、铁观音和祁门红茶进行检测，并建立SIMCA(基于主成分分析的一种分类方法)判别模型，龙井、碧螺春、祁门红茶和铁观音的识别准确率分别为90％、80％、100％和
100％；李晓丽等对西湖龙井、浙江龙井、羊岩勾青、
雪水云绿和庐山云雾茶(各30份)进行近红外检测，并采用主成分分析法建立模型，除西湖龙井和浙江龙井存在部分重叠外，其余均能很好的区分。
[0005] 同位素是所有
生物的一个自然标签，与生物的生长环境密切相关，因此同位素质谱(IRMS)为茶叶原产地识别提供科学的、可靠的鉴别方法。IRMS目前已经广泛应用在各种农产品的原产地检测中， 等利用同位素质谱和
核磁共振光谱，结合主成分分析，能够理想地鉴别斯洛文尼亚三个不同地区的红酒。Brescia等利用IRMS测定牛奶中的δ13C、δ15N和Ba元素含量，区分了不同原产地的牛奶，证明了IRMS适用于乳制品。Martinelli等对来自美国、南美、欧洲和澳大利亚的气泡葡萄酒进行同位素检测，发现具有显著差异。Tamara等测定了43份印度、23份斯里兰卡和12份中国茶叶中稳定同位素，非线性分析表明茶叶原产地的判定易受判别函数的影响，很好的区分不同国家产区的茶叶。
[0006] 王睿等采用电感耦合
等离子体原子发射光谱法，对新疆6个主要产地的36个石榴样品的可食部分(果肉)和籽中12种金属元素的含量进行测定，采用主成分分析PCA和线性判别分析LDA对石榴可食部分和籽中金属元素进行综合评价。结果表明：PCA得出2个三因子模型，分别解释了石榴可食部分和籽中金属元素数据的84.29％和60.33％；通过对石榴可食部分中金属元素组成进行PCA，PCA可将36个石榴样品划分为6类，与实际产地吻合。陈辉等应用电感耦合等离子体质谱测定了北京顺义、河北阜平和河北平山三个地区65个荆条蜜样品中38种元素含量，并采用PCA和反向传输人工神经网络按照不同产地对荆条蜜进行了分析，交叉检验整体准确率为95.4％。
[0007] 不同农产品因所处地理、
气候、环境等的不同，常规和特殊化学成分的种类、含量都有所不同，可通过色谱技术测定其种类和含量的差异将来自不同产地的同一种农产品区分开，达到产地溯源的目的。Longobardi等采用高效液相色谱测定希腊的初榨
橄榄油中的15种氨基酸等含量，最后利用方差分析得到26个样品之间存在极显著差异(p＜0.01)，该方法能有效的将不同来源的橄榄油区分开来。Kücük等研究了来自不同地区的三种土
耳其蜂蜜的12种氨基酸含量，对比发现由于产地的不同，各种不同产地氨基酸的含量有较大差异(p＜0.05)，可进行不同产地来源判断。Collomb等对瑞士不同海拔地区生产的牛奶进行色谱分析，测定其氨基酸组分及含量，发现三个不同海拔地区牛奶的氨基酸组分及含量存在明显差异。
[0008] LF等采用反相高效液相色谱法测定绿茶、红茶和黑茶中儿茶素、咖啡因等，来自5个不同国家28批次的茶叶(产自中国的红茶、绿茶、黑茶；产自日本的绿茶；产自斯里兰卡的红茶；产自肯尼亚的红茶以及产自印度的红茶)可采用PCA进行区分区分。Kodama等采用毛细管
电泳测定了日本静冈(n＝4)、鹿儿岛(n＝4)、三重县(n＝4)产的茶叶中7种儿茶素(+C、-C、EC、CG、ECG、EGC、EGCG)和咖啡因含量，采用PC(主成分分析)，识别准确率100％。
[0009] 电子舌是感官仿生技术在仪器分析中的一种应用，其原理是用多传感器阵列代替生物体的味觉系统，并由计算机代替大脑进行分析。当检测溶液
吸附在传感器上时，传感器的电容发生改变，产生特征响应
信号，由电子舌自带的专家
自学习系统进行
模式识别，做出定性和定量分析。传感器技术是一种快速简便的检测手段，且需要的样品量非常少；缺点是与其他检测方法相比，重复性较差，无法表征单一组分的差异。传感器技术也已有用来鉴定原产地。Guadarrama等利用电子舌对橄榄油进行质量控制，成功辨别出不同品种和不同原产地的西班牙橄榄油；Penza等利用膜传感器阵列结合人工神经网络，能够有效地分辨出9个不同原产地的意大利红酒。
[0010] 从上述例子可以看出，国内外对地理标志产品的鉴别方法已有很多，但是很多研究尚有不足之处，如
采样不充分样品数量少，无法保证样品的准确性和代表性；样品空间选择跨度大，往往选自不同国家、不同地区，本身就具有很大的区别；再加之甚至选择了不同品种样品进行比较，不同品种间本身差异较大，因此这类鉴别方法对小范围的地理标志产品产地判别借鉴意义不大；建模方法采用单种检测数据结合计量学方法进行，单种检测数据无法代表产地溯源的全部信息，导致产地识别率较低，上述这些严重影响了地理标志产品保护技术的创新和突破。针对如上情况，有必要建立基于多种检验技术的武夷岩茶产地智能鉴别模型，即一种联合近红外、稳定同位素、微量元素、氨基酸、儿茶素以及电子舌数据的武夷岩茶产地鉴别方法。(三)发明内容
[0011] 本发明目的在于解决单种检测数据无法代表产地溯源的全部关键信息以及传统计量学方法存在的
缺陷等问题，提供基于多种检验技术的武夷岩茶产地智能鉴别方法，联合近红外光谱、稳定同位素、微量元素、氨基酸、儿茶素及电子舌数据建立的武夷岩茶产地识别模型技术方法，该方法基于偏最小二乘判别模型，将不同产地岩茶(包括地理标志产地内和产地外岩茶)近红外特征光谱数据、稳定同位素数据、微量元素数据、氨基酸、儿茶素及电子舌数据融合在一起，建立分析模型，提取试样后利用模型客观、准确的判定岩茶产地。
[0012] 本发明采用的技术方案是：
[0013] 基于多种检验技术的武夷岩茶产地智能鉴别方法，即融合近红外光谱、稳定同位素、微量元素、氨基酸、儿茶素以及电子舌数据鉴别武夷岩茶产地的方法，所述方法包括：
[0014] (A)采集不同产地岩茶样品：
[0015] 武夷岩茶产区外样品数＞100份，且产区周边50公里范围内样品占比＞50％；武夷岩茶产区内样品数为产区外样品的2～3倍，采样范围
覆盖主要产区各生产企业，且每家企业应不少于3个样品；
[0016] (B)测定不同产地岩茶样品的近红外特征光谱数据：
[0017] 64次扫描，特征光谱条带取其平均值，扫描范围为12000-4000cm-1，数据点的间隔为1.928em-1，室温25℃，湿度保持稳定，
无损检测，无需采用
粉碎等前处理，采用相同加料方法、加料量，加料完毕即可进行检测。(C)测定不同产地岩茶样品的氢、
氧、氮、
碳四种稳定同位素质谱数据：
[0018] δ13C、δ15N、δ18O、δ2H、δ86Sr等稳定同位素含量测定，每个样品至少重复分析3次以上，取平均值作为最终结果。
[0019] 通过SVM-RFE(支持向量机回归特征消去法)对武夷岩茶稳定同位素数据进行训练和预测，随机重复100次，并对各变量的模型特征进行排序，筛选出岩茶原产地的同位素特征变量，其排序顺序为氢、氧、氮、碳、锶。并利用预测集计算模型的灵敏度增维
精度、
分辨率增维精度、识别率增维精度，均按重复计算100次平均结果，氢、氧、氮、碳四种数据组成的模型，其识别率最高，达93.93％，因此建模只需要选择氢、氧、氮、碳四种数据即可，无需对锶等其他稳定同位素含量进行检测。
[0020] (D)测定不同产地岩茶样品的铯、
铜、
钙、铷微量元素数据
[0021] 用原子吸收光谱仪测定Ca、Mg、Mn元素含量，用电感耦合等离子体质谱测Ti、Cr、Co、Ni、Cu、Zn、Rb、Cd、Cs、Ba、Sr元素含量。干茶样品
微波消解，消解完毕，观察消解液是否澄清，若浑浊，则重复压
力消解步骤，若完全澄清，则定容后用上述仪器进行测定。
[0022] 通过SVM-RFE法对微量元素数据进行训练和预测，随机重复100次，并对各变量的模型特征进行排序，筛选出岩茶原产地的微量元素特征变量，并通过预测集计算每一维变量累加之后的模型增维精度，得到铯、铜、钙、铷、锶、钡特征排序顺序。然后对特征变量按自然序逐级进行组合，并利用预测集计算模型的灵敏度增维精度、分辨率增维精度、识别率增维精度，由铯、铜、钙、铷四者微量元素组成的模型，其识别率增维精度最高，说明此四种微量元素间的信息具有互补性，只需要选择建模的铯、铜、钙、铷四种微量元素进行检测即可，无需对其它微量元素进行检测。
[0023] (E)测定不同产地岩茶样品的氨基酸数据：
[0024] 利用HPLC法对不同产地岩茶样品中的27种氨基酸进行检测，平行测定两次，取平均值。
[0025] 通过SVM-RFE法对武夷岩茶氨基酸组分数据进行训练和预测，随机重复100次，并对各变量的模型特征进行排序，筛选出茶叶原产地的特征变量，并通过预测集计算每一维变量累加之后的模型增维精度，确定其排序顺序为天冬酰胺、脯氨酸、
色氨酸、磷
乙醇胺、尿素和缬氨酸。然后对特征变量按自然序逐级进行组合，并利用预测集计算模型的灵敏度增维精度、分辨率增维精度、识别率增维精度，由天冬酰胺、脯氨酸、色氨酸、磷乙醇胺四种氨基酸组成的模型，其识别率最高，说明此四种氨基酸间的信息具有互补性，只需要选择建模的天冬酰胺、脯氨酸、色氨酸、磷乙醇胺四种氨基酸进行检测即可。
[0026] (F)测定不同产地岩茶样品的儿茶素数据：
[0027] 利用HPLC法对不同产地岩茶样品中的6种儿茶素以及咖啡
碱进行检测，平行测定3次，取平均值。
[0028] 采用支持向量机回归特征消去法对儿茶素和咖啡碱共7个特征变量进行逐级组合后，各儿茶素和咖啡碱对于地理特征的贡献率从高到低依次为表没食子儿茶素(EGC)、儿茶素(C)、表没食子儿茶素
没食子酸酯(EGCG)、没食子酸(GA)、表儿茶素(EC)、表儿茶素没食子酸酯(ECG)和咖啡碱。然后对特征变量按自然序逐级进行组合，并利用预测集计算模型的灵敏度、分辨率、识别率，最高模型识别率为0.8596，模型中包含EGC、C、EGCG、GA和EC，该模型的灵敏度为0.9322，分辨率为0.6734。基于儿茶素和咖啡碱数据建立的SVM-RFE模型，其灵敏度增维精度均高于0.9000，说明对于保护区内的岩茶样品判别结果较为可靠。而其分辨率均较低，说明对于地理标志保护区外的假冒岩茶样品容易发生误判。在EGC和C的SVM模型中，当增加了EGCG和GA后其识别率略有下降，说明EGC和C这两个儿茶素与EGC和C之间相关性亦较强，但无法提供增益的产地信息。但当EC变量被纳入模型时，模型识别率达到最高，说明EC是对EGC和C的产地信息表征的有益补充。ECG和咖啡碱再纳入模型中，识别率增维精度不断下降，说明ECG和咖啡碱不能对产地判别提供增益效益，且可能与前面的五种儿茶素之间存在一定的负相关关系，因此建模采用EGC、C、EGCG、GA和EC五种儿茶素含量。(G)测定不同产地岩茶样品的电子舌数据：
[0029] 本研究采用ASTREE Electronic电子舌系统，该套系统具有一个
电极(Ag/AgCl)以及7个独立传感器(ZZ、BA、BB、CA、GA、HA、JB)，电子舌系统使用前用NaCl和C5H8NO4Na标准溶液对传感器进行活化，并依次进行信号校准和仪器
自诊断，自诊断通过后可进行样品滋味的分析检测，检测完毕后，将电子舌数据以文本形式导出，取各传感器最后十秒(第91-100秒)稳定数据的均值作为检测数据进行后续分析。
[0030] (H)联合近红外、稳定同位素、微量元素、氨基酸、儿茶素以及电子舌数据建立不同产地岩茶鉴别
数据库[0031] (1)将每条近红外数据(Y轴数据)拼接在Excel数据表格中，每行所有列数据组成每条近红外数据；
[0032] (2)将每个样品的稳定同位素数据按氢、氧、氮、碳顺序拼接在近红外数据后，将微量元素数据按铯、铜、钙、铷拼接在稳定同位素后，将氨基酸数据按天冬酰胺、脯氨酸、色氨酸、磷乙醇胺顺序拼接在微量元素后，将儿茶素数据按EGC、C、EGCG、GA和EC顺序拼接在氨基酸数据后，最后将电子舌7个独立传感器(按ZZ、BA、BB、CA、GA、HA、JB顺序)最后十秒(第91-100秒)稳定数据的均值拼接在儿茶素数据后，武夷岩茶产区内样品组成的Excel数据表，以data1命名；武夷岩茶产区外样品组成的Excel数据表，以data2命名；
[0033] (3)运行MATLAB
软件中的edit功能，打开data1.xls、data2.xls，以Mat文件格式保存，文件名对应为data1.mat、data2.mat；
[0034] (4)数据分割：参考R.W.Kennard和L.A.Stone在Computer aided design of experiments中的方法，随机选取武夷岩茶产区内总数的65～70％作为原产地内模型数A1，随机取武夷岩茶产区外65～70％作为原产地外模型数A2，建立kenstone分割程序，设计思路：首先计算两两样本之间距离，选择距离最大的两个样品；然后分别计算剩余的样本与已选择的两个样本之间的距离；对于每个剩余样本而言，其与已选样品之间的最短距离被选择，然后选择这些最短距离中相对最长的距离所对应的样本，作为第三个样品；重复上述步骤，直至所选的样品的个数等于事先确定的数目为止。
[0035] (5)K-折交互验证法：将样本数据集随机划分为K个子集(一般是均分)，将一个子集作为验证集，其余的K-1组子集作为训练集；轮流将K个子集作为验证集，交叉重复K次，得到K次的结果，并用K次结果的平均值作为分类器或模型的性能指标。K-折法下，每一个样本数据都被用作
训练数据，也被用作验证数据，避免了过度学习和欠学习状态的发生。
[0036] (6)偏最小二乘法鉴别模型的建立：对步骤(4)和(5)数据分割后的融合近红外、稳定同位素、微量元素、氨基酸、儿茶素以及电子舌数据，采用偏最小二乘法分析并建立PLSDA模型；
[0037] (I)取待测未知产地样品按照上述步骤B、C、D、E、F和G，测定近红外特征光谱数据、稳定同位素质谱数据、微量元素数据、氨基酸数据、儿茶素以及电子舌数据，将测得数据代入上述PLSDA模型，若预测结果小于0，则判断待测样品为武夷岩茶产地外样品；若预测结果大于0，则判断待测样品为武夷岩茶产地内样品。
[0038] 具体的，所述步骤(H)中分割程序分别为：[model1，test1]＝kenstone(data1，A1)和[model2，test2]＝kenstone(data2，A2)，得到model1、test1、model2、test2。
[0039] 偏最小二乘法的建模思路：能够在自变量存在严重多重相关性的条件下进行回归建模，在偏最小二乘回归模型中，设法去掉不太重要的相关性变量，剩余的自变量回归系数将更容易解释，且更易于辨识系统信息与噪声，避免将本应保留的系统信息舍弃。具体的，所述步骤(H)中偏最小二乘法鉴别模型的建立过程如下：
[0040] (a)合并训练集：xxxc＝[data1(model1，：)；data2(model2，：)]；
[0041] (b)合并预测集：xxxp＝[data1(test1，：)；data2(test2，：)]；
[0042] (c)求训练集平均光谱：mx＝mean(xxxc)；
[0043] (d)训练集减去平均光谱：xxxc＝xxxc-ones(A，1)*mx；
[0044] A为：A1+A2；
[0045] (e)预测集减去平均光谱：xxxp＝xxxp-ones(B，1)*mx；
[0046] B为：原产地内测试集数B1与原产地外测试集数B2之和；
[0047] (f)响应变量：yyc＝-ones(A，2)；yyc(1：A1，1)＝1；yyc(A1+1：A，2)＝1；
[0048] A1与B1之和为原产地内样品总数C1；
[0049] A2与B2之和为原产地外样品总数C2；
[0050] (g)用K-折交互验证：
[0051] indices＝crossvalidation(′Kfold′，x，k)；
[0052] (h)建模过程：
[0053] [betattt，www，BETAPLS1]＝plsbasetotal(xxxc，yyc(：，1)，lvp1)；
[0054] [betattt，www，BETAPLS2]＝plsbasetotal(xxxc，yyc(：，2)，lvp2)；
[0055] cy＝[xxxc*BETAPLS1(：，lvp1)，xxxc*BETAPLS2(：，lvp2)]；
[0056] py＝[xxxp*BETAPLS1(：，lvp1)，xxxp*BETAPLS2(：，lvp2)]；
[0057] [rrt，cyy]＝max(cy′)；
[0058] [rwwrt，pyy]＝max(py′)；
[0059] (i)计算训练过程中模型的灵敏度和分辨率：
[0060] err01＝length(find(cyy(1：A1)＝＝1))/A1；
[0061] err02＝length(find(cyy(A1+1：A1+A2)＝＝2))/110；
[0062] (j)计算预测未知样品过程中模型的灵敏度和分辨率：
[0063] err1a＝length(find(pyy(1：B1)＝＝1))/B1；
[0064] err1b＝1-length(find(pyy(B1+1：B1+B2)＝＝1))/B2；
[0065] (k)保存预测结果：save cyy cyy；save pyy pyy；
[0066] (l)py的第一列即为预测结果。
[0067] 可画图给出详细结果：
[0068] bar(cy(：，1))；
[0069] figure
[0070] bar(py(：，1))
[0071] 其中kenstone程序如下所示：
[0072]
[0073]
[0074]
[0075] 其中Kfoldcv程序如下所示：
[0076]
[0077]
[0078]
[0079] 其中plsbasetotal程序如下所示：
[0080]
[0081] 本发明的有益效果主要体现在：本发明基于偏最小二乘判别模型，将不同产地岩茶(包括地理标志产地内和产地外岩茶)近红外特征光谱数据、稳定同位素数据、微量元素数据、氨基酸数据、儿茶素以及电子舌数据融合在一起，建立分析模型，提取试样后利用模型客观、准确的判定岩茶产地，识别率最高，可达100.0％，高于单种数据PLSDA的判别结果。(四)具体实施方式
[0082] 下面结合具体
实施例对本发明进行进一步描述，但本发明的保护范围并不仅限于此：
[0083] 实施例1：
[0084] A、采集不同产地岩茶样品
[0085] 国标(GB/T 18745-2006)中规定了武夷岩茶的地理保护范围，即福建省武夷山市行政区划内，本发明在武夷岩茶地理标志保护区武夷街道、崇安街道、上梅、星村、五夫、岚谷、新丰街道、洋庄、兴田、下梅、吴屯11个行政区域内进行样品采集，每个行政区域内随机选择3个取样点(分别以A、B、C标示)，共33个取样点，取样范围基本涵盖主要产区，每个采样点取样15份(分别以A-1、A-2......A-15进行标示)，获得495份地理标志保护区武夷岩茶样品，另在福建省除武夷山市外其他县市(建阳、建瓯、漳州、泉州、松溪、政和)，以及广西、贵州、江西(婺源、赣州)等11个地点收集保护区外岩茶样品，每个地点取样15份(分别以1、2......15进行标示)，获得165个非地理标志保护区岩茶样品。地理标志产地内样品数与地理标志产地外样品数之比为3∶1。
[0086] B、不同产地岩茶近红外特征光谱数据
[0087] 无损检测、布鲁克TENSOR37，采用相同加料方法、加料量(目测)，加料完毕即进行检测。表1为15个兴田A样品部分近红外数据表，其中X轴为
波长范围，Y轴为吸光度。
[0088] 表1：15个兴田A样品部分近红外数据表
[0089]
[0090]
[0091] C、不同产地岩茶氢、氧、氮、碳四种稳定同位素质谱数据
[0092] δ13C、δ15N、δ18O、δ2H、δ86Sr由Thermo Fisher MAT253稳定同位素质谱仪测定。地理标志产区内、外岩茶样品同位素比率经上述方法检测，表2节选部分岩茶样品同位素比率数据表。
[0093] 表2：部分岩茶样品同位素比率统计表
[0094]
[0095]
[0096] 通过SVM对武夷岩茶稳定同位素数据进行训练和预测，随机重复100次，并对各变量的模型特征进行排序，筛选出岩茶原产地的同位素特征变量顺序为氢、氧、氮、碳、锶；并利用预测集计算模型的灵敏度、分辨率、识别率，均按重复计算100次平均结果，见表3。
[0097] 表3：同位素特征变量组合结果情况
[0098]
[0099]
[0100] 由表3可知，当氢和氧的同位素数据进行组合后，模型识别率下降，说明氧和氢对原产地特征的贡献具有较强的相关性；而加入碳和氮同位素数据后，模型识别率上升，达到93.93％，说明氮和碳具有较好的互补，因此建模只需要选择氢、氧、氮、碳四种数据即可，建模中无需增加锶的数据，在实际检测中，同位素锶的含量无需检测。
[0101] D、测定不同产地岩茶样品的铯、铜、钙、铷四种微量元素数据
[0102] 待测茶叶微波消解后，观察消解液是否澄清，若浑浊，则重复压力消解步骤，若完全澄清，使用日立180-50原子吸收光谱仪测定样液中Ca、Mg、Mn元素含量，使用Thermo Fisher XSeries II电感耦合等离子体质谱测定微波消解液中Ti、Cr、Co、Ni、Cu、Zn、Rb、Cd、Cs、Ba、Sr元素含量。表4为节选部分地理标志产区内、外岩茶样品微量元素数据表。
[0103] 表4：部分岩茶样品微量元素统计表
[0104]
[0105]
[0106] 通过SVM-RFE对微量元素数据进行训练和预测，随机重复100次，并对各变量的模型特征进行排序，筛选出岩茶原产地的微量元素特征变量，并通过预测集计算每一维变量累加之后的模型增维精度，得到铯、铜、钙、铷、锶、钡特征排序顺序。然后对特征变量按自然序逐级进行组合，并利用预测集计算模型的灵敏度增维精度、分辨率增维精度、识别率增维精度，由铯、铜、钙、铷四者微量元素组成的模型，其识别率增维精度最高达0.8121，说明此四种微量元素间的信息具有互补性，只需要选择建模的铯、铜、钙、铷四种微量元素即可。
[0107] E、测定不同产地岩茶样品的氨基酸数据
[0108] 采用高效液相衍生化法对不同产地岩茶样品中的27种氨基酸进行检测，平行测定两次，取平均值，部分岩茶样品7种氨基酸含量数据见表5。
[0109] 表5 不同产地部分岩茶7种氨基酸含量(单位：百分比)
[0110]
[0111]
[0112] 通过SVM对武夷岩茶氨基酸组分数据进行训练和预测，随机重复100次，并对各变量的模型特征进行排序，筛选出茶叶原产地的特征变量，并通过预测集计算每一维变量累加之后的模型增维精度，确定其排序顺序为天冬酰胺、脯氨酸、色氨酸、磷乙醇胺、尿素和缬氨酸。然后对特征变量按自然序逐级进行组合，并利用预测集计算模型的灵敏度增维精度、分辨率增维精度、识别率增维精度，由天冬酰胺、脯氨酸、色氨酸、磷乙醇胺四者氨基酸组成的模型，其识别率增维精度最高达0.78，说明27种氨基酸间的信息互补性较弱，只需要选择建模的四种氨基酸即可。
[0113] F、测定不同产地岩茶样品的儿茶素数据
[0114] 采用高效液相法对不同产地岩茶样品中的6种儿茶素和咖啡碱进行检测，平行测定两次，取平均值，部分岩茶样品儿茶素和咖啡碱含量数据见表6。
[0115] 表6 不同产地岩茶儿茶素和咖啡碱含量
[0116]
[0117]
[0118] 采用支持向量机回归特征消去法对儿茶素和咖啡碱共7个特征变量进行逐级组合后，对于地理特征的贡献率从高到低依次为EGC、C、EGCG、GA、EC、ECG和咖啡碱。然后对特征变量按自然序逐级进行组合，并利用预测集计算模型的灵敏度、分辨率、识别率，模型最高识别率为0.8596，灵敏度为0.9322，分辨率为0.6734，模型中包含EGC、C、EGCG、GA和EC。基于儿茶素和咖啡碱数据建立的模型，灵敏度均高于0.9000，说明对于保护区内的岩茶样品判别结果较为可靠；而分辨率均较低，说明对于地理标志保护区外的假冒岩茶样品易发生误判。在EGC和C的模型中，增加了EGCG和GA后，识别率略有下降，说明EGC和C与EGC和C之间相关性亦较强，但无法提供增益的产地信息。但当EC变量被纳入模型时，模型识别率达到最高，说明EC是对EGC和C的产地信息表征的有益补充。ECG和咖啡碱再纳入模型中，识别率增维精度不断下降，说明ECG和咖啡碱不能对产地判别提供增益效益，且可能与前面的五种儿茶素之间存在一定的负相关关系，因此建模采用EGC、C、EGCG、GA和EC五种儿茶素含量。
[0119] G、测定不同产地岩茶样品的电子舌数据
[0120] ASTREE Electronic电子舌系统，7个独立传感器名称分别为ZZ、BA、BB、CA、GA、HA、JB。各传感器最后十秒(第91-100秒)稳定数据的均值作为检测数据。其数据见表7。
[0121] 表7部分岩茶7个传感器最后10秒均值数据表
[0122]
[0123] H、融合近红外、稳定同位素、微量元素、氨基酸、儿茶素和电子舌数据建立不同产地岩茶鉴别数据库
[0124] (1)将每条近红外数据拼接在Excel数据表格中，每行所有列数据组成每条近红外数据。
[0125] (2)将每个样品的稳定同位素数据按氢、氧、氮、碳顺序，微量元素按铯、铜、钙、铷顺序，氨基酸按天冬酰胺、脯氨酸、色氨酸、磷乙醇胺顺序，儿茶素按EGC、C、EGCG、GA和EC顺序，电子舌按ZZ、BA、BB、CA、GA、HA、JB顺序，依次拼接在近红外数据后，地理标志产区内样品组成495行、4172列(近红外Y轴数据共4148列，依次增加氢、氧、氮、碳、铯、铜、钙、铷、天冬酰胺、脯氨酸、色氨酸、磷乙醇胺、EGC、C、EGCG、GA、EC、ZZ、BA、BB、CA、GA、HA、JB变为4172列)的Excel数据表，以data1命名；地理标志产区外样品组成165行、4172列的Excel数据表，以data2命名。
[0126] (3)运行MATLAB软件中的edit功能，打开data1.xls、data2.xls，以Mat文件格式保存，文件名对应为data1.mat、data2.mat；
[0127] (4)数据分割：参考R.W.Kennard和L.A.Stone在Computer aided design of experiments上所述方法，原产地内随机取330(A1)个样品作为模型，原产地外随机取110(A2)个样品作为模型，建立kenstone分割程序，[model1，test1]＝kenstone(data1，330)和[model2，test2]＝kenstone(data2，110)，得到model1、test1、model2、test2。
[0128] (5)对上述数据建模PLSDA：
[0129] ①合并训练集：xxxc＝[data1(model1，：)；data2(model2，：)]；
[0130] ②合并预测集：xxxp＝[data1(test1，：)；data2(test2，：)]；
[0131] ③求训练集平均光谱：mx＝mean(xxxc)；
[0132] ④训练集减去平均光谱：xxxc＝xxxc-ones(440，1)*mx；
[0133] ⑤预测集减去平均光谱：xxxp＝xxxp-ones(220，1)*mx；
[0134] ⑥响应变量：yyc＝-ones(440，2)；yyc(1：330，1)＝1；yyc(331：440，2)＝1；
[0135] ⑦用K-折交互验证：
[0136] indices＝crossvalidation(′Kfold′，x，k)；
[0137] ⑧建模过程：[betattt，www，BETAPLS1]＝plsbasetotal(xxxc，yyc(：，1)，lvp1)；
[0138] [betattt，www，BETAPLS2]＝plsbasetotal(xxxc，yyc(：，2)，lvp2)；
[0139] cy＝＝[xxxc*BETAPLS1(：，lvp1)，xxxc*BETAPLS2(：，lvp2)]；
[0140] py＝[xxxp*BETAPLS1(：，lvp1)，xxxp*BETAPLS2(：，lvp2)]；
[0141] [rrt，cyy]＝max(cy′)；
[0142] [rwwrt，pyy]＝max(py′)；
[0143] ⑨计算训练过程中模型的灵敏度和分辨率：
[0144] err01＝length(find(cyy(1：330)＝＝1))/330；
[0145] err02＝length(find(cyy(331：440)＝＝2))/110；
[0146] ⑩计算预测未知样品过程中模型的灵敏度和分辨率：
[0147] err1a＝length(find(pyy(1：165)＝＝1))/165；
[0148] err1b＝1-length(find(pyy(166：220)＝＝1))/55；
[0149] 保存预测结果：save cyy cyy；save pyy pyy；
[0150] py的第一列为预测结果，小于0为原产地外样品，大于0为原产地内样品，画图给出详细结果：
[0151] bar(cy(：，1))；
[0152] figure
[0153] bar(py(：，1))；
[0154] I、PLSDA模型识别率
[0155] 按上述建模方法分别对近红外光谱、稳定同位素、微量元素、氨基酸、儿茶素、电子舌以及六者融合数据进行建模分析，其结果见表8。
[0156] 表8：PLSDA分类模型判别结果汇总
[0157]
[0158] 从表8可知，各检测技术所代表的特征指标之间具有互补性，采用本发明PLSDA建模方法分别对近红外光谱、稳定同位素、微量元素、氨基酸、儿茶素、电子舌以及六者融合数据进行分析，以六者融合数据为模型数据时，识别率最高，达100.0％，远高于单种数据的PLSDA判别结果。
[0159] J、盲样检测
[0160] 盲样监督小组从武夷岩茶农户处购买岩茶样品，监控晒青、做青、杀青等步骤，确保岩茶样品的原产地域属性，上述样品作为盲样中的地理标志产区内样品；从建阳、建瓯、婺源等地购买岩茶，作为盲样中的地理标志产区外样品，上述盲样与建模岩茶样品出自不同厂家。分析检测人员事先未得知待测盲样的产地属性，随机选取若干份，检测，再按本发明方法对盲样产地属性进行判断，并与盲样监督小组进行核对，确定盲样识别率。将20、60、100个盲样按照步骤B、C、D、E和F获得的盲样数据代入上述PLSDA模型，判断其地理标志属性，其判别准确率均达到100.0％。
[0161] 实施例2：
[0162] 采用与实施例1相同的建模方法，数据分割用Kenstone分割程序，用K-折交互验证，分别建立PLSDA、神经网络ELM和最小二乘支持向量机LS-SVM模型，近红外数据不变，稳定同位素、微量元素、氨基酸、儿茶素和电子舌分别按照氢、氧、氮、碳、锶、Cs、Cu、Ca、Rb、Sr、Ba、天冬酰胺、脯氨酸、色氨酸、磷乙醇胺、尿素、缬氨酸、EGC、C、EGCG、GA、EC、ECG、咖啡碱、ZZ、BA、BB、CA、GA、HA、JB拼接在近红外数据后，其模型识别率分别为89.5％、83.2％、87.7％。
[0163] 实施例3：
[0164] 采用与实施例1相同的建模方法，数据分割用Kenstone分割程序，用K-折交互验证，分别建立PLSDA、神经网络ELM和最小二乘支持向量机LS-SVM模型，近红外数据不变，稳定同位素、微量元素、氨基酸、儿茶素和电子舌分别按照氢、氧、氮、碳、锶、Cs、Cu、Ca、Rb、Sr、Ba、天冬酰胺、脯氨酸、色氨酸、磷乙醇胺、尿素、缬氨酸、EGC、C、EGCG、GA、EC、ZZ、BA、BB、CA、GA、HA、JB拼接在近红外数据后，其模型识别率分别为90.1％、83.8％、88.9％。
[0165] 实施例4：
[0166] 采用与实施例1相同的建模方法，数据分割用Kenstone分割程序，用K-折交互验证，分别建立PLSDA、神经网络ELM和最小二乘支持向量机LS-SVM模型，近红外数据不变，稳定同位素、微量元素、氨基酸、儿茶素和电子舌分别按照氢、氧、氮、碳、锶、Cs、Cu、Ca、Rb、Sr、Ba、天冬酰胺、脯氨酸、色氨酸、磷乙醇胺、EGC、C、EGCG、GA、EC、ZZ、BA、BB、CA、GA、HA、JB拼接在近红外数据后，其模型识别率分别为96.2％、85.3％、89.6％。
[0167] 实施例5：
[0168] 采用与实施例1相同的建模方法，数据分割用Kenstone分割程序，用K-折交互验证，分别建立PLSDA、神经网络ELM和最小二乘支持向量机LS-SVM模型，近红外数据不变，稳定同位素、微量元素、氨基酸、儿茶素和电子舌分别按照氢、氧、氮、碳、锶、Cs、Cu、Ca、Rb、天冬酰胺、脯氨酸、色氨酸、磷乙醇胺、EGC、C、EGCG、GA、EC、ZZ、BA、BB、CA、GA、HA、JB拼接在近红外数据后，其模型识别率分别为98.7％、86.9％、90.2％。
[0169] 实施例6：
[0170] 采用与实施例1相同的建模方法，数据分割用Kenstone分割程序，用K-折交互验证，分别建立PLSDA、神经网络ELM和最小二乘支持向量机LS-SVM模型，近红外数据不变，稳定同位素、微量元素、氨基酸、儿茶素和电子舌分别按照氢、氧、氮、碳、Cs、Cu、Ca、Rb、天冬酰胺、脯氨酸、色氨酸、磷乙醇胺、EGC、C、EGCG、GA、EC、ZZ、BA、BB、CA、GA、HA、JB拼接在近红外数据后，其模型识别率分别为100.0％、88.0％、89.6％。
[0171] 实施例7～11：
[0172] 岩茶样品、近红外数据、同位素数据、微量元素数据、氨基酸数据、儿茶素数据以及电子舌数据等与实施例1相同，分别采用不同方法进行鉴别。实施例9-11分别采用现有发明
专利CN103630528A(
申请号201210307144.2)、CN102455320A(申请号201010526790.9)、CN103245713A(申请号201310095950.2)所述方法进行鉴别，实施例2～6与实施例1指标区别见表9。
[0173] 表9：实施例1与实施例7～11比较情况表
[0174]
[0175]
[0176] 从上述比较结果可以得出，采用本发明方法，其盲样检测识别率达100.0％以上，可作为武夷岩茶产地溯源识别技术方法。
