基于加权模糊识别的中药材质量评估及产地分析模型

1、基于加权模糊识别的中药材质量评估及产地分析模型作者：兰继斌 林健 宋玉风毕业论文2010-11-19 2:06:08【摘要】 研究通过采集药材HPLC指纹图谱数据，运用加权模糊模式识别方法评估中药材的质量及其产地分布情况，得到了很好的分辨率。最后通过一个实例说明了此模型实际应用的可行性与优越性。 【关键词】 指纹图谱； 模糊模式识别； 加权； 阀值Abstract As the quality of the medicines on the market varies greatly, based on the weighted fuzzy pattern recognition of qua

2、lity assessment and origin analysis model of Chinese herbal medicines, a method using HPLC fingerprint of the medicines is presented to assess the quality and origin in this paper . The numerical experiments show that the method has good resolution .By an example, it is verified that this model is m

3、ore feasible and effective than the previous model in practice. Key words fingerprint; fuzzy pattern recognition; weighted; threshold 我国是中草药资源丰富的国家，对中药材的药性、成分等方面的研究也有着优良的传统。中药指纹图谱1，2就是在对中药系统的不断深入研究中发展起来的，是表征中药所含成分与其质量关系的有效手段。目前，在应用中药指纹图谱研究药材成分的过程中，得到了许多不错的结果3，4。特别是与交叉学科的相互渗透方面5，6，其研究水平和应用范围也在不断扩大。由于

4、气候及产地等因素的影响，中药材的成分和质量等级也有着差异性。鉴于此，我们有必要探询一种有效的药材质量评估方式，来识别其等级及产地分布情况。1 仪器及分析条件71.1 仪器与试剂 岛津LC-10ATVP高效液相色谱仪；岛津SPD10AVP紫外检测器；TB328B型电光分析天平（上海天平仪器厂）；HH-2型电子恒温水浴锅（金坛市新航仪器厂）；RE-52AA型旋转蒸发仪（上海亚荣生化仪器厂）。甲醇为色谱纯，水为重蒸水，其余试剂为分析纯。1.2 色谱条件 色谱柱：岛津ULTRONVX-ODS(4.6mm250mm,5l)柱；流动相：甲醇-水-冰乙酸（60：40：1）；流速：0.8ml/min；检测波长

5、：254nm；柱温：室温；进样量：20l。2 模型的建立 Step1. 设所考虑的对象为A药材，取m种已知质量等级并有代表性的该药材A1（产地1），A2（产地2），Am（产地m），它们相应的等级评价为“一等、二等、m等”。测定不同产地此药材的HPLC指纹图谱，采集各等级药材在实验中的指纹图谱的保留时间和峰面积。 Step2. 选择t个共有色谱峰作为指纹图谱的特征峰，其中药材Aj观测到的第i个特征峰的保留时间为aji，峰面积为bji，i=1,2,，t, j=1,2,m。建立表征指纹图谱中保留时间及峰面积的信息矩阵Hj=aj1 aj2 ajtbj1 bj2 bjt，并对矩阵Hj每行分别应用公式rj

6、i=aji t i=1aji，Rji=bji t i=1bji归一化得：hj=rj1 rj2 rjtRj1 Rj2 Rjt。 Step3. 将药材Aj归一化后的信息矩阵hj转化成信息向量：Yj=(rj1 rj2 rjt Rj1 Rj2 Rjt）, j=1，2,m。 Step4. 在同样实验条件下，对一未知质量等级的待测A药材A进行观测，得到t个共有色谱峰的保留时间ai和峰面积bi，i=1,2,t。建立表征指纹图谱中保留时间及峰面积的信息矩阵H=a1 a2 atb1 b2 bt,对矩阵H的每一行分别应用公式ri=ai t i=1ai，Ri=bi t i=1bi归一化得：h=r1 r2 rtR1

7、R2 Rt。 Step5. 将药材A归一化后的信息矩阵h转化成信息向量：X=(r1 r2 rt R1 R2 Rt）。 Step6. 根据药材中各成分所对应色谱峰面积的大小，以及该成分对药效的重要程度，综合考虑给出1t号峰所对应的一个权重向量w=(w1,w2,wt)。并根据实验中保留时间的稳定性优于峰面积的稳定性这一事实，考察两者稳定性的相对程度，给出指纹图谱中保留时间和峰面积的权重指标=(1,2)。 Step7. 利用综合加权最小平均贴近度公式 N(X,Yj)=1t i=1wi(rirji)+2t i=1wi(RiRji) 1 21t i=1wi(ri+rji)+2t i=1wi(Ri+Rji

8、)计算N(X，Y1)，N(X，Y2)，,N(X，Ym)的值，其中“”为取小算子。 Step8. 计算N(X，Yj0)=maxN(X，Y1)，N(X，Y2)，N(X，Ym)，考察实验中不同产地该药材的差异性程度，设置阀值。当N(X，Yj0)时，则认为X与Yj0最贴近。由择近原则8，最后判断X与Yj0为一类。当N(X，Yj0)时，则判断X的质量不合格。3 实例分析 白芷在我国生长区域分布十分广泛，其主要产于浙江A1（一等品）、河南A2（二等品）、四川绵阳A3（三等品）、湖北A4（四等品）、湖南A5（五等品），建立HPLC法测定上述各产地何首乌的HPLC指纹图谱，如图1所示，并采集各等级白芷在实验中

9、指纹图谱的保留时间和峰面积7。3.1 选择10个共有色谱峰作为指纹图谱的特征峰，其中药材Aj观测到的第i个特征峰的保留时间为aji，峰面积为bji，i=1,2,10，j=1,2,5。3.2 建立表征指纹图谱中保留时间及峰面积的信息矩阵，并对该矩阵每行分别归一化并转化成信息向量：Yj=(rj1 rj2 rj10 Rj1 Rj2 Rj10)，rj1 、Rj1具体数值分别如表1、表2所示，i=1,2,10，j=1,2,5。表1 110号峰保留时间归一化后数值（略）表2 110号峰峰面积归一化后数值（略）3.3 现有一个等待定级的白芷样本A，在同样实验条件下，得到10个共有色谱峰的保留时间ai和峰面积

10、bi，i=1,2,10。建立表征指纹图谱中保留时间及峰面积的信息矩阵并对每行分别归一化得：按上述公式转化成信息向量X=(r R)，其中： r=(0.0177 0.0303 0.0429 0.0591 0.0735 0.0857 0.0998 0.1372 0.1776 0.2762 ) R=( 0.1936 0.0235 0.1168 0.0286 0.0029 0.0325 0.0449 0.3183 0.1225 0.1165)3.4 根据白芷各成分所对应色谱峰面积的大小，以及该成分对药效的重要程度，综合考虑给出110号峰所对应的一个权重向量w=(0.16,0.04,0.12,0.04,0

11、.04,0.08,0.08,0.20,0.12,0.12)。并根据实验中保留时间的稳定性优于峰面积的稳定性这一事实，考察两者稳定性的相对程度，给出白芷指纹图谱中保留时间和峰面积的权重指标=(0.6,0.4)。3.5 利用综合加权最小平均贴近度公式计算得：N(X，Y1)=96.782%，N(X，Y2)=94.804%，N(X，Y3)=97.645%，N(X，Y4)=91.634%，N(X，Y5)=96.790%。 根据白芷的成分组成以及评估人员对药材质量的要求，给出阀值=0.85。 综上可知该待测白芷样品为三等品，其产于四川一带。4 结束语 本研究通过采集中药材HPLC指纹图谱数据，建立符合药材

12、特征的数学模型，分析待测药材的质量等级及产地分布，得到了很好的辨别效果，该模型可以广泛应用于各种中药材的分析和评估。【参考文献】 1 王瑾,胡璟华. 中药指纹图谱的研究概述. 解放军药学学报 ,2004,3 : 214216.2 杨桂芳,王鹏. 中药指纹图谱的建立方法和作用. 时珍国医国药, 2004,1 :21.3 郑云枫,彭国平. 栀子药材HPLC指纹图谱的研究. 南京中医药大学学报(自然科学版) ,2005,1: 4547.4 易中宏,胥秀英,郑一敏,王龙. 茯苓药材指纹图谱模糊模式识别研究. 现代中药研究与实践 ,2004,4:62.5 袁鹏,张铭光,袁敏等. 模糊聚类分析法在郁金裂解色谱测定指纹图谱中的应用研究. 计算机与应用化学 ,2005,