中药色谱指纹图谱相似度算法的改进(2)
(1)
其中xaj表示样品指纹图谱中非共有峰j的峰面积，nt表示共有峰数。
尽管d值大小可表示样品指纹图谱向量X与对照指纹图谱向量Y的相似程度，但取值范围在(0~∞)之间，通过负指数函数能把取值在(0~∞)变化的函数转换为取值(0~1)函数的特点[10]，将公式（1）转化为公式(2)，其中S代表样品的相似度:
S=e-d (2)
公式(2)中S取值域在(0~1)之间，但如果直接使用公式(2)进行相似度计算，那么不同中药的判断阈值必定各不相同，不能满足“不同中药评价判断阈值统一为0.9，相似度大于0.9就可以判定样品质量合格”这个条件，必须对公式(2)进一步优化。若引入了一个特征系数ω，得到公式(3)：
S=e-ωd (3)
令ω=ln0.9/-d阈（d阈代表样品判断阈值对应的d值），则不同中药的判断阈值都能统一为0.9。具体公式如下：
(4)
对于每一种中药，都有一个对应的ω值，实际应用时ω大小由中药药学实验确定，目前水平下ω的具体数值由d阈[6]（本文指最差合格样品对应的d值）计算。ω可以作为企业或者实验室的保密值，当他人无法得到ω的具体值并且得不到足够多的标准品时，就很难对特定中药进行相似度评价。
2．实验部分
仪器：岛津LC-10ADvp高压泵，岛津SPD-M10Avp二极管阵列检测器，岛津CLASS-vp工作站；柱温箱:TC-100(天津市金洲科学仪器有限公司)，循环微波萃取装置。
药材：香丹注射液（上海通用药业股份有限公司，批号：031082、0310119、0310116、0310127、031062、0310117、0310130、031096、0310112、0310106，共计10批次），香丹注射液（上海中西药业股份有限公司，批号：0305011、010701，共计2批次）；按《中国药典》规定中药取样法于河南省封丘县采集的金银花药材[11]，共11批次，编号：hf1~hf11；同法取自湖南省淑浦和湖北省恩施的金银花药材，编号分别为nga1和nga2。
香丹注射液色谱条件：色谱柱：1. Alltima C18(5 μm，4.6 mm×250 mm), 2. Purospher Start RP-18e(5 μm, 4.6 mm×250 mm)；柱温：30℃；流动相: 1%冰醋酸水溶液(A)-1%冰醋酸甲醇溶液(B)采用梯度洗脱：10%B100%B(055 min), 100%B10%B(5560 min), 10%B(6090 min)；流速： 1.0 mL/min；检测波长：280 nm；进样量：20 L。
金银花色谱条件：A(乙腈)-B(酸水，水-磷酸=100: 0.4)进行梯度洗脱：0~10 min (10%A~10%A), 10~60 min (10%A~30%A), 60~75 min (30%A~90%A), 75~80 min(90%A~90%A), 80~90 min(90%A~10%A),90~120 min(10%A~10%A)；流速1.0 mL/min，进样体积20 μL，检测波长237 nm。
香丹注射液稀释溶液的配置：精确量取0.50 mL香丹注射液于10 mL容量瓶，定容至刻度，之后用微孔滤膜过滤，过滤后即可用于色谱实验。
金银花待测溶液的配制：将金银花粉碎至20目，精确称取已粉碎的金银花1.000 g，置于微循环萃取装置进行萃取，萃取剂为15.00 mL55%乙醇溶液，并将得到的萃取液用55%乙醇溶液定容至25.00 mL。取上述溶液2.50 mL继续用55%乙醇溶液稀释至10.00 mL，过C18固相萃取柱，用微孔滤膜过滤，得到的滤液即可用于色谱实验。
3.结果与讨论
本文所有d值和S值的计算是通过运行于WindowsTM 8.0系统的Microsoft® Excel 2010进行的。
3.1香丹注射液相似度计算
图2-1~2-12为12个批次香丹注射液的指纹图谱，峰面积数据经逐一匹配后得到表1，将上海通用药业的八个批次（批号：031082、0310119、0310116、0310127、031062、0310117、0310130、031096）作为标准系列，得到一个标准，然后将同一厂家的0310106、0310112批次样品看作相近样品与之对比，将上海中西药业的0305011、010701批次样品看作差别样品进行对比。
上一篇： 指示剂置换法用于生物巯基化合物的检测