本发明涉及一种贵重香料产地类别的鉴别方法，特别涉及一种沉香类别的鉴别方法。
背景技术：
沉香是指白木香属(aquilariaspp.)树木在生长过程中形成的由木质部组织及其分泌物共同组成的混合物质。沉香以其药用价值和芳香气味而著称，广泛应用于中药、香薰、香水、宗教庆典等领域。中药沉香，为瑞香科植物白木香aquilariaspp.含有树脂的木材，分布于广东、海南、广西、福建等地。具有行气止痛，温中止呕，纳气平喘之功效。常用于胸腹胀闷疼痛，胃寒呕吐呃逆，肾虚气逆喘急。奇楠是沉香中最为昂贵和珍贵的一种，主要产于惠安系地区，被认为是市面上最高品质的沉香，具有神秘的东方香气，与其他普通沉香不同，奇楠不需要点燃既能闻到香味。广泛认为奇楠和惠安系沉香可以通过香气来区分。沉香的主要成分由色酮、倍半萜和其他芳香化合物构成，其中倍半萜和其他芳香化合物是主要发香成分。文献记载，目前发现的色酮物质有80多种，倍半萜有70多种。目前，国内对沉香成分检测已有研究，成分检测大多集中于色酮物质，而对倍半萜研究相对较少，对沉香分类研究更是凤毛麟角。在我国，现行的沉香规范为《中华人民共和国药典》。2015年版《中华人民共和国药典》(一部)通过对沉香的乙醇抽提物、显色反应、薄层色谱以及液相特征图谱几项测试对药用沉香的品质加以规范，其中液相特征图谱项目通过对特征峰(沉香四醇以及几种色酮)对沉香进行判别。在我国尚无对沉香分类的规定。目前对沉香类别的鉴定主要是人工定级，即定级人员凭借经验通过烘烤、燃烧等方法对产生的气味进行判别，从而对沉香分类。但是人工辨别容易受主观性和偶然性干扰，且经验不足者，不能进行判定。即使定级人员进行鉴定，其鉴定的准确性由于受多种主观或客观因素的影响，鉴定结果的准确性也不高，而且需要消耗大量的沉香，对沉香样品需求量大，导致鉴定成本高，鉴定时间长。本发明从气味角度出发，模拟人工香味鉴别过程，对沉香进行分级。通过对气味成分的分析，用数据取代定级人员的经验，避免主观性干扰，同时无定级经验者亦可以通过该方法对沉香类别进行鉴定。本发明首先对沉香进行热处理，即模拟沉香加热/燃烧；并收集产生的气体；接着通过气相色谱-质谱(gc-ms)联用的方法对气体成分进行分析，根据或依据沉香热挥发成分中的气味成分及其相对含量，选择沉香中热挥发成分中具有鉴别特征的化学成分(包括倍半萜及其他芳香族化合物)作为观测值并建立fisher判别函数、bayes判别函数；然后将鉴定沉香样品中的热挥发成分中具有鉴别特征的化学成分分别代入建立的fisher判别函数、bayes判别函数，进行计算；最后根据计算值进行判断。利用热处理-气相色谱质谱联用对沉香进行分类的方法，更具体地，是指利用气相色谱质谱联用分析沉香热处理时产生的化学成分，主要针对倍半萜和其他芳香化合物，从而区分沉香类别。技术实现要素：本发明的目的是针对现有沉香类别鉴定中存在的技术缺陷，提供一种奇楠沉香与惠安系普通沉香的鉴别方法，本发明方法利用微量沉香样品即可对沉香的类别进行鉴定和分类。本发明方法从气味角度出发，通过模拟人工气味鉴别过程，对沉香进行分级。本发明方法避免了鉴别过程中人为主观性和偶然性带来的干扰，鉴定结果准确；而且对于气味品鉴无定级经验者也可以通过本发明的方法对沉香进行分类，避免了沉香鉴别必须需要定级人员进行鉴别的缺陷，提高了鉴定效率、鉴定的准确性，降低了沉香鉴定的成本，缩短沉香鉴定时间。为实现本发明的目的，本发明一方面提供一种奇楠沉香与惠安系普通沉香的鉴别方法，包括如下顺序进行的步骤：1)将待鉴别沉香样品置于密封的容器中，进行加热处理；2)对加热处理后的容器内的气体进行气相色谱-质谱分析，并计算特征指标化合物的相对含量；3)fisher判别法将待鉴别沉香的gc-ms分析的特征指标化合物的相对含量代入惠安系普通沉香与奇楠沉香的fisher判别分析方程(即公式fh-q)中，获得相应的fisher判别值yfh-q，然后计算待测样品的判别值yfh-q与奇楠、惠安系普通沉香标准品的fisher判别重心坐标值之间的差距；如果待测样品的判别值yfh-q与奇楠沉香标准品的fisher判别重心坐标值之间的差值的绝对值小于待测样品的判别值yfh-q与惠安系普通沉香标准品的fisher判别重心坐标值之间的差值的绝对值，则判定该样品为奇楠沉香；如果待测样品的判别值yfh-q与奇楠沉香标准品的fisher判别重心坐标值之间的差值的绝对值大于待测样品的判别值yfh-q与惠安系普通沉香标准品的fisher判别重心坐标值之间的差值的绝对值，则判定该样品为惠安系普通沉香。其中，步骤1)中所述加热处理温度为120-180℃，优选为160℃；所述加热处理时间≥50min，优选为60min。特别是，还包括将待鉴定沉香样品研磨成粒度为40-60目的沉香粉末后，再进行所述的加热处理。尤其是，进行加热处理的沉香样品的质量为10-30mg，优选为20mg。特别是，加热处理是时沉香中的热挥发性成分挥发至密闭容器中。其中，步骤2)中所述特征指标化合物共15种化合物，其编号依次为x1-15,15种化合物的casno(即cas登录号，即美国化学会的美国化学文摘服务社(chemicalabstractsservice,cas)化学物质登录号)依次为000104-20-1、000128-37-0、000489-29-2、000637-69-4、003321-66-2、003691-12-1、005951-61-1、005956-09-2、005956-12-7、017066-67-0、028976-67-2、071596-72-0、1000192-43-5、1000374-19-7、117066-77-0。即cas登录号依次为000104-20-1、000128-37-0、000489-29-2、000637-69-4、003321-66-2、003691-12-1、005951-61-1、005956-09-2、005956-12-7、017066-67-0、028976-67-2、071596-72-0、1000192-43-5、1000374-19-7、117066-77所代表的化合物。特别是，步骤3)中所述惠安系普通沉香与奇楠沉香的fisher判别分析方程为：yfh-q＝-1.306×x1+0.159×x2+8.000×x3+2.890×x4+0.172×x5+9.884×x6+16.271×x7-0.845×x8-0.516×x9+0.762×x10-1.102×x11+10.947×x12+38.369×x13+2.749×x14+1.714×x15-6.446(公式fh-q)其中，x1-x15为编号依次为1-15的15种化合物gc-ms分析的相对含量。尤其是，所述奇楠沉香与惠安系普通沉香的fisher判别分析方程按照如下步骤建立：3a)将42个沉香标准品分为2组，其中第1组为奇楠沉香组，记为1.00，共有16个沉香标准品；第2组为惠安系普通沉香组，记为2.00，共有26个沉香标准品；3b)将沉香标准品分别置于密封的容器中，进行加热处理；3c)对加热处理后的容器内的气体分别进行气相色谱-质谱分析，以追踪到的150种化学成分作为观测值，并分别计算各化学成分的相对含量；3d)采用逐步判别法从150种化学成中筛选得到15种特征指标化合物；3e)对两组沉香标准品中的15种特征指标化合物的相对含量进行fisher判别分析的计算，获得所述的惠安系普通沉香与奇楠沉香的fisher判别分析方程(即公式fh-q)。特别是，步骤3b)中将奇楠、惠安系普通沉香标准品分别研磨成粒度为40-60目的沉香粉末后，再进行所述的加热处理。尤其是，进行加热处理的沉香标准品的质量为10-30mg，优选为20mg。其中，步骤3)中所述奇楠沉香标准品的fisher判别重心坐标值为6.546；所述惠安系普通沉香标准品的fisher判别重心坐标值为-4.028。特别是，步骤3)中所述星洲系标准沉香、惠安系标准沉香的fisher判别分析方程(即公式fh-x)按照如下方法获得：对两组42个沉香标准品，及15种特征指标化合物的相对含量进行fisher判别分析的计算，代入fisher判别分析方程(即公式fh-q)公式计算，得到所述奇楠标准沉香、惠安系标准沉香的判别分析方程，其中：a为使δ(a)达到极大值的列向量(a可结合定义及拉格朗日乘子法计算出，具体数学表达为式中a为组内离差阵，b为组间离差阵，λ为a-1b的特征根，a即a-1b的特征向量)，a’为a的行向量(即a的转置)，为42个样品的总均值(42个样品各特征化合物的相对含量平均值组成的列向量)。特别是，步骤3)中所述奇楠标准沉香、惠安系标准沉香的fisher判别重心坐标值按照如下方法获得：对两组42个沉香标准品，及15种特征指标化合物的相对含量进行fisher判别分析的计算，代入fisher判别分析法的重心坐标公式计算，分别得到所述奇楠标准沉香、惠安系标准沉香的重心坐标值，其中：dj为样品重心坐标值(d1表示奇楠沉香重心坐标值，d2表示惠安系沉香重心坐标值)，a为使δ(a)达到极大值的列向量(a可结合定义及拉格朗日乘子法计算出，具体数学表达为式中a为组内离差阵，b为组间离差阵，λ为a-1b的特征根，a即a-1b的特征向量)，a’为a的行向量(即a的转置)，为42个样品的总均值(42个样品各特征化合物的相对含量平均值组成的列向量)，为各组样品的平均值(即本发明中表示16个奇楠标准沉香样品各特征化合物的相对含量平均值组成的列向量、表示26惠安系样品各特征化合物的相对含量平均值组成的列向量)。本发明另一方面提供一种奇楠沉香与惠安系普通沉香的鉴别方法，包括如下顺序进行的步骤：1)将待鉴别沉香样品置于密封的容器中，进行加热处理；2)对加热处理后的容器内的气体进行气相色谱-质谱分析，并计算特征指标化合物的相对含量；3)bayes判别分析法将待鉴别沉香的gc-ms分析的特征指标化合物的相对含量代入奇楠沉香与惠安系普通沉香的bayes判别分析方程(即公式b1h-q和b2h-q)中，计算待鉴定样品的yb1、yb2值，如果计算值yb1>yb2，则判定该样品属于第1类，为奇楠沉香；如果计算值yb2>yb1，则判定该样品属于第2类，为惠安系沉香。其中，步骤1)中所述加热处理温度为120-180℃，优选为160℃；所述加热处理时间≥50min，优选为60min。特别是，还包括将待鉴定沉香样品研磨成粒度为40-60目的沉香粉末后，再进行所述的加热处理。尤其是，进行加热处理的沉香样品的质量为10-30mg，优选为20mg。特别是，加热处理是使沉香中的热挥发性成分挥发至密闭容器中。其中，步骤2)中所述特征指标化合物共15种化合物，其编号依次为x1-15,15种化合物的casno(即cas登录号，即美国化学会的美国化学文摘服务社(chemicalabstractsservice,cas)化学物质登录号)依次为000104-20-1、000128-37-0、000489-29-2、000637-69-4、003321-66-2、003691-12-1、005951-61-1、005956-09-2、005956-12-7、017066-67-0、028976-67-2、071596-72-0、1000192-43-5、1000374-19-7、117066-77-0。其中，步骤3)中所述奇楠沉香与惠安系普通沉香的bayes判别函数(即公式b1h-q和b2h-q)为：yb1h-q＝-15.302×x1+2.720×x2+105.280×x3+39.591×x4+2.261×x5+131.868×x6+202.707×x7-10.335×x8-6.704×x9+9.831×x10-12.109×x11+142.119×x12+496.505×x13+35.467×x14+22.588×x15-88.193(公式b1h-q)yb2h-q＝-1.496×x1+1.036×x2+20.684×x3+9.029×x4+0.447×x5+27.344×x6+30.649×x7-1.404×x8-1.249×x9+1.773×x10-0.460×x11+26.354×x12+90.760×x13+6.393×x14+4.466×x15-6.712(公式b2h-q)特别是，步骤3)中所述奇楠沉香与惠安系普通沉香的bayes判别函数(即公式b1h-q和b2h-q)按照如下步骤建立：3a)将42个沉香标准品分为2组，其中第1组为奇楠沉香组，记为1.00，共有16个沉香标准品；第2组为惠安系沉香组，记为2.00，共有26个沉香标准品；3b)将沉香标准品分别置于密封的容器中，进行加热处理；3c)对加热处理后的容器内的气体分别进行气相色谱-质谱分析，以追踪到的150种化学成分作为观测值，并分别计算各化学成分的相对含量；3d)采用逐步判别法从150种化学成中筛选得到15种特征指标化合物；3e)对两组沉香标准品中的15种特征指标化合物的相对含量进行bayes计算，获得奇楠沉香与惠安系普通沉香的bayes判别分析方程(即公式b1h-q和b2h-q)，其中，公式b1h-q、b2h-q具体如下：yb1h-q＝-15.302×x1+2.720×x2+105.280×x3+39.591×x4+2.261×x5+131.868×x6+202.707×x7-10.335×x8-6.704×x9+9.831×x10-12.109×x11+142.119×x12+496.505×x13+35.467×x14+22.588×x15-88.193(公式b1h-q)yb2h-q＝-1.496×x1+1.036×x2+20.684×x3+9.029×x4+0.447×x5+27.344×x6+30.649×x7-1.404×x8-1.249×x9+1.773×x10-0.460×x11+26.354×x12+90.760×x13+6.393×x14+4.466×x15-6.712(公式b2h-q)特别是，步骤3b)中将奇楠、惠安系普通沉香研磨成粒度为40-60目的沉香粉末后，再进行所述的加热处理。尤其是，进行加热处理的沉香标准品的质量为10-30mg，优选为20mg。与现有技术相比，本发明的优点如下：本发明通过在密闭容器中加热处理沉香标准品，模拟人工气味鉴别，沉香中的挥发性物质挥发，对沉香挥发性物质进行gc-ms分析，再对沉香挥发性物质及其相对含量采用逐步判别法进行筛选，确认奇楠与惠安系沉香的特征指标成分，并对特征指标成分的相对含量进行fisher判别分析的计算和bayes判别分析计算，分别获得相应的fisher判别方程、重心坐标和bayes判别方程，然后将待鉴定沉香样品的特征指标成分的gc-ms分析的相对含量代入相应的fisher判别方程和bayes判别方程进行判断。本发明方法判别结果准确，鉴定效率高，而且降低了沉香鉴定的成本，缩短了鉴定时间；克服了现有人工气味鉴别过程中人为主观性和偶然性的缺陷。具体实施方式下面结合具体实施例来进一步描述本发明，本发明的优点和特点将会随着描述而更为清楚。但这些实施例仅是范例性的，并不对本发明的范围构成任何限制。本领域技术人员应该理解的是，在不偏离本发明的精神和范围下可以对本发明技术方案的细节和形式进行修改或替换，但这些修改和替换均落入本发明的保护范围内。下述实施例中所使用的实验方法如无特殊说明，均为常规方法。下述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。本发明实施例中使用的沉香标准品“星洲系、惠安系、奇楠”均购自印尼、马来西亚、文莱、巴布亚新几内亚、苏门答腊等地；惠安系和奇楠购自中国、越南、印度、泰国、柬埔寨等地；其他待分类鉴别的沉香样品a、b由北京沉香协会提供；沉香样品c、d购自流通领域市场；沉香样品e、f采香于产地，6个样品均具备报关单、申购单等资料证明，且经过北京沉香协会专家团鉴定。奇楠是沉香中的上品，因其独特的味道和稀缺的资源，而价值连城。奇楠主要产自中国、越南、印度、泰国、柬埔寨等地，也称为惠安产区。该产区除了产出奇楠，也产出惠安系普通沉香，然而二者价值悬殊。没有经验的消费者很难将二者区分，本实施例将通过本发明对奇楠和普通沉香进行分类。实施例1“惠安系普通沉香”、“奇楠沉香”样品预处理沉香分为普通沉香和奇楠两种，普通沉香分为星洲系和惠安系普通沉香，奇楠沉香产于惠安产区，是沉香中的名贵品种，也是沉香中价值最高的香品。奇楠的研究是沉香研究的重大课题，将惠安产区的奇楠和普通沉香有效区分能保护沉香市场，防止不良商家以次充好，欺骗消费者。本实施例通过本发明对奇楠和惠安系普通沉香进行区分。一、试验器具的预处理选取密封性能良好且对实验结果无干扰的实验器具(即玻璃瓶：本试验选用的为安捷伦5188-275920ml顶空进样瓶；玻璃针管(即气相色谱分析的采样针管，100ml)；金属针头等)于烘箱内180℃预热2小时以上，以祛除试验器具和空气中的挥发性气体对试验的干扰，以使空白峰面积不能超过样品峰面积的2％，即视为实验器具准备完备，可以用于试验。二、沉香样品的预处理将42个沉香标准样品，其中奇楠沉香标准品16个(编号为第1组，记为1.00)，惠安系普通沉香标准品26个(编号为第2组,记为2.00)和4个样品(a、b、e、f)，分别取少量研磨，制成粒度为40-60目的标准沉香粉末，备用。实施例2“惠安系普通沉香”、“奇楠沉香”样品的fisher判别方程yh-q的建立一、沉香标准样品的加热处理1、分别取实施例1制备的惠安系普通沉香、奇楠沉香标准品粉末各20mg，分别置于不同的、且预热处理后的玻璃顶空进样瓶中，拧紧密封；2、将装有沉香粉末的顶空瓶分别置于烘箱中，加热升温至160℃，并在保持温度为160℃的条件下保持加热1h，即对沉香粉末于160℃的条件下加热1h。本发明加热处理的温度除了160℃之外，其他加热温度120-180℃均适用于本发明，加热时间除了1h之外，通常加热时间≥50min均适用于本发明。本发明预热温度比沉香样品加热处理过程中的加热温度高，在180℃时祛除挥发性物质，确保对试验没有干扰。二、gc-ms分析1、气体收集将加热处理1h后的每个顶空瓶分别通过两个金属针头，同时连接氮气及预处理后的空玻璃针管，其中氮气针头插入玻璃瓶底部，空玻璃针管针头插入玻璃瓶2-3mm。氮气缓慢通入玻璃瓶中，将玻璃瓶中气体推入玻璃针管(采样玻璃针管)，至针管收集满100ml气体；2、gc-ms分析将采样针管内的气体注入气质联用仪的tenax采样管，进行气相色谱-质谱分析(即gc-ms分析)，其中：气质联用仪选用：agilenttechnologies7890a/5975c；气相色谱分析条件如下：色谱柱：dm726641a：325℃:50m×320μm×1μm进样口温度300℃，氦气流速3ml/min，不分流。采用程序升温进行分离：初温40℃，保持1min,以5℃/min的速率升到150℃，再以1℃/min的速率升到170℃，保持3min，再以1℃/min的速率升到180℃，保持3min，再以10℃/min的速率升到200℃，再以20℃/min的速率升到280℃，保持1min。数据处理：输出最小峰面积为最大峰面积的1％；质谱数据库选择为nist14数据库(其他质谱数据库也适用于本发明)。三、fisher判别分析模型1、建立判别函数1a、筛选特征化学成分42个沉香标准品分为两组，奇楠沉香16个(第1组，记为1.00)；惠安系普通沉香26个(第2组，记为2.00)；通过gc-ms分析，共追踪150种(用v001-150表示)化学成分作为观测值，并计算其相对含量。针对追踪到的150种化学成分的相对含量采用逐步回归分析法(又称逐步判定法)筛选出15种化学成分作为鉴定惠安系普通沉香和奇楠沉香的特征化学成分，其中，15种特征化学成分的分子量、casno.如表1所示。其中，逐步判别分析法筛选特征化学成分首先将42个样品的150种化学成分建立矩阵。分别计算其λ统计量、f值，筛选遵循两个原则：a、每次选择λ统计量(广义方差，也称为wilks的lambda统计量)最小者进入模型。首先选择λ最小进入模型，然后重新计算未被选中的其他变量与选中变量的λ统计量，剩余统计量最小者与已经进入模型的变量搭配进入模型。b、当p＜α入时进入，当p＞α出时剔除，其中f值是组间均方与组内均方的比。其中规定α入＝0.05，α出＝0.10。边进入边剔除，一步一步以此类推，逐步筛选出特征化学成分，筛选结果见表1。表1惠安系普通沉香、奇楠沉香fisher判别函数系数、特征化学成分分子量及其casno.其中，casno.又称为cas登录号，即美国化学会的美国化学文摘服务社(chemicalabstractsservice,cas)化学物质登录号。1b、建立fisher判别方程fisher判别法(又称费希尔判别法、费舍判别法)是由fisher在1936年提出，是根据方差分析的思想建立起来的一种能较好区分各种总体的线性判别法，其基本思想是投影(降维)。将k类p维数据投影到某一个方向，使得投影后类与类之间尽可能分开。fisher判别借助了一元方差分析的思想衡量类与类之间是否分开。即不同类之间离散程度最大化，并使每类内部离散程度最小，进而导出判别函数。fisher判别式按照各类总体中方差尽可能小，不同类中的均值之间差距尽可能大的原则，即每个总体内样本离散度越小越好。因此，依据两个原则建立数学模型：1、希望各类总体的平均值相差越大越好；2、希望每个总体的离差平方和越小越好。对42个沉香标准样品及其15种特征化合物的相对含量建立矩阵模型gmn(其中m＝42，n＝15)，42个样品分为2组，即j＝1，2(j表示组别总体；1表示奇楠，2表示惠安系)。因此建立基本表达式：其中：a为使δ(a)达到极大值的向量，是投影法线；ssr为组间平方和；sse为组内平方和；a为组内离差阵(又称交叉乘积阵)；b为组间离差阵。fisher判别函数(即判别分析方程)定义为：y＝a′x+c重心坐标公式：其中a’为a的行向量(即a的转置)，为各组样品的平均值(本实施例中表示16个奇楠样品各特征化合物的相对含量平均值组成的列向量、表示26惠安系样品各特征化合物的相对含量平均值组成的列向量)，为42个样品的总均值(42个样品各特征化合物的相对含量平均值组成的列向量)。然后对fisher判别函数求解，并采用不加权法计算奇楠、惠安系普通沉香两组的重心，得到判别函数值及重心坐标。本实施例中根据42个沉香标准样中的15种特征化学成分的相对含量，采用spss19.0软件对特征成分的相对含量进行fisher判别分析的计算，获得惠安系普通沉香与奇楠沉香的fisher判别分析方程(公式fh-q)，公式中的各样品的特征化学成分(x1-15)以15种特征化学成分的相对含量表示：yfh-q＝-1.306×x1+0.159×x2+8.000×x3+2.890×x4+0.172×x5+9.884×x6+16.271×x7-0.845×x8-0.516×x9+0.762×x10-1.102×x11+10.947×x12+38.369×x13+2.749×x14+1.714×x15-6.446(公式fh-q)2、判别函数的特征描述spss19.0软件同时给出了本发明实施例中fisher判别分析模型函数(即公式fh-q)的特征跟和λ统计量，分别如表2、3所示。表2奇楠系与惠安系沉香的fisher判别函数yfh-q特征值判别函数特征根方差的％累积％相关系数公式fh-q27.690100.0100.00.982表3奇楠系与惠安系沉香的fisher判别的λ统计量函数λ统计量卡方自由度p公式fh-q0.035109.088150.000由表2、3中的参数对fisher判别分析模型的描述可知，本发明方法建立的惠安系普通沉香与奇楠沉香的fisher判别分析模型函数的特征根为27.690，该特征根完全解释了方差，相关系数是0.982；概率p值为0.000，表明该函数判别效果显著，可见本发明中关于判别奇楠、惠安系普通沉香的fisher判别分析模型适用。3、计算奇楠沉香、惠安系普通沉香标准品的重心坐标值采用spss19.0软件计算奇楠沉香标准品、惠安系普通沉香标准品的重心坐标值，计算结果如表4：表4奇楠、惠安系普通沉香组重心处的函数值分类重心坐标值第1组(奇楠)6.546第2组(惠安系普通)-4.028实施例3“惠安系普通沉香”、“奇楠沉香”样品的判别一、待鉴定沉香样品的加热处理除了分别取用实施例1制备的沉香样品a、b、e、f粉末各20mg之外，其余与实施例2的步骤(一)相同。二、gc-ms分析与实施例2的步骤(二)相同。测定计算待鉴定沉香样品中15种特征化学成分的相对含量，测定结果如表5所示。表5奇楠沉香与惠安系普通沉香fisher鉴别的特征化学成分的相对含量、判定结果三、沉香的fisher判别将样品a、b、e、f的经gc-ms分析后的15种化学成分的相对含量分别代入公式fh-q，计算得到yfh-q值，分别为ye：6.4398；yf：6.5593；ya：-3.0500；yb：-6.1083。分别计算每个样品的yfh-q值与两组(即奇楠、惠安系普通沉香标准品组)重心坐标值，相近者即为分类。例如，ye＝-35.9385，分别计算与第1组重心坐标(6.546)和第2组重心坐标(-4.028)的距离，经计算|ye-6.546|＝0.1026＜|ye-(-4.028)|＝10.4678，可以与第1组重心坐标更相近，确定e为奇楠。同理可知，a为惠安系，b为惠安系，f为奇楠沉香。如果待测样品的fisher判别值yfh-q与奇楠标准沉香样品的重心坐标值之间的差值的绝对值等于待测样品的fisher判别值yfh-q与惠安系标准沉香样品的重心坐标值之间的差值的绝对值，则使用bayes判别分析法进行判别。实施例4“惠安系普通沉香”、“奇楠沉香”样品的bayes判别方程bh-q的建立一、沉香标准样品的加热处理与实施例2相同。二、gc-ms分析与实施例2相同。三、bayes判别分析模型1、建立判别函数42个沉香标准品分为两组，奇楠沉香16个(编号为第1组，记为1.00)，惠安系普通沉香26个(编号为第2组,记为2.00)；通过gc-ms分析，共追踪150种(用v001-150表示)化学成分作为观测值，并计算其相对含量。经逐步回归法筛选15种化学成分为鉴定奇楠沉香与惠安系普通沉香的特征化学成分，其分子量、casno.如表1所示。bayes判别分析(又称贝叶斯判别分析)，其思想源于贝叶斯统计。贝叶斯判别方法依据如下两个原则建立数学模型：1、后验概率最大化准则；2、平均损失最小化准则。得到判别函数。在样品分类判别时，把对应项的数据代入判别函数中，得到的数值较大者，即为该样品的分类。bayes数学模型其中，r(d)为平均错判损失若存在d*满足上述数学模型，则判定准则d*为bayes判别准则。bayes判别函数定义为其中称为判别系数，为常数项式中：s为总体样品协方差阵j为样品的分类组别，本发明中42个样品分为2组，即j＝1，2(j表示组别总体；1表示奇楠，2表示惠安系)；为各组样品的平均值，本实施例中表示16个奇楠样品各特征化合物的相对含量平均值组成的列向量、表示26惠安系样品各特征化合物的相对含量平均值组成的列向量；qj为j类样品的先验概率在样品分类判别时，把对应特征化学成分的相对含量代入判别函数中，得到的数值较大者，即为该样品的分类。由于计算量较大，本实施例采用spss19.0软件对样品中化学成分的相对含量进行bayes计算，获得bayes判别函数系数(如表6所示)，进而得到奇楠沉香与惠安系普通沉香的bayes判别函数(即公式b1h-q、b2h-q)，公式中的各样品观测值(x1-15)以15种特征化学成分的相对含量表示。表6奇楠沉香、惠安系普通沉香的bayes判别函数系数yb1h-q＝-15.302×x1+2.720×x2+105.280×x3+39.591×x4+2.261×x5+131.868×x6+202.707×x7-10.335×x8-6.704×x9+9.831×x10-12.109×x11+142.119×x12+496.505×x13+35.467×x14+22.588×x15-88.193(公式b1h-q)yb2h-q＝-1.496×x1+1.036×x2+20.684×x3+9.029×x4+0.447×x5+27.344×x6+30.649×x7-1.404×x8-1.249×x9+1.773×x10-0.460×x11+26.354×x12+90.760×x13+6.393×x14+4.466×x15-6.712(公式b2h-q)2、bayes判别分析模型评价采用初始验证正确率及交叉验证正确率对奇楠沉香与惠安系普通沉香的bayes判别公式进行评价，采用spss19.0软件进行计算，spss19.0软件给出了奇楠沉香与惠安系普通沉香的bayes判别公式的判别结果(如表7)。其中：初始验证正确率，即将42个沉香标准样品各自对应的15种特征化学成分的相对含量分别带入公式b1h-q和b2h-q，经计算比较得到的预测结果与实际的样品分类相比较得出的数据；交叉确认法正确率，即将42个样品中的某个样品单独提出，用其他剩余41个样品建立bayes判别模型和判别函数，将该单独提出的样品对应物质的相对含量带入bayes判别函数，经计算比较得到预测结果与实际的样品分类相比较，42个样品分别重复上述过程，而得到的数据。本实施例中，采用spss19.0软件进行计算，spss19.0软件同时给出了本实施例的bayes判别结果验证。可以看出初始验证正确率100％，交叉验证正确率100％，具体结果如表7所示。可见本发明中关于判别奇楠沉香、惠安系普通沉香的bayes判别分析模型适用于判别奇楠沉香、惠安系普通沉香的判别。表7判别结果验证实施例5“惠安系普通沉香”、“奇楠沉香”样品的判别一、待鉴定沉香样品的加热处理除了分别取用实施例1制备的沉香样品a、b、e、f粉末各20mg之外，其余与实施例2的步骤(一)相同。二、gc-ms分析与实施例2的步骤(二)相同。计算待鉴别沉香样品的特征化学成分的相对含量，测定结果如表8所示。三、沉香的bayes判别将样品a、b、e、f的经gc-ms分析后的15种特征化学成分的相对含量分别代入公式b1h-q、b2h-q，分别计算得到yb1、yb2值。计算结果见表8。表8奇楠沉香与惠安系普通沉香bayes鉴别的特征化学成分的相对含量、判定结果分别对每个待鉴定样品的yb1、yb2进行比较，较大数值作为分类。例如a样品中特征化学成分的相对含量代入公式b1h-q、b2h-q，计算得到yb1＝-34.5221，yb2＝11.0282，yb1＜yb2，则属于第2类，属于惠安系普通沉香；再如e样品，经计算yb1＞yb2，则属于第1类，属于奇楠沉香。通过比较可知a、b为惠安系普通沉香；e、f为奇楠沉香。如果待测样品的bayes判别值yb2＝yb1，则使用fisher判别分析法。当前第1页12