发明名称 一种色谱指纹图谱相似度计算方法
摘要 一种色谱指纹图谱相似度评价方法,包括以下步骤:采用标准样品建立对照指纹图谱;采用与建立对照指纹图谱相同的色谱条件建立待评价样品的样品指纹图谱;将样品指纹图谱的每一个峰值与对照指纹图谱中对应峰值相比,得到样品的比值指纹图谱向量和对照比值指纹图谱向量;计算这两个比值向量的相似度S。本发明既能反映待评价样品指纹图谱与标准指纹图谱之间共有峰的相似性,也考虑了非共有峰数目及其大小对产品质量的影响,能整体反应各样品指纹图谱之间相似程度,该相似度的大小能够灵敏地定量表达样品指纹图谱与标准指纹图谱定性和定量差异,可靠地鉴别药品、食品的真伪,更好地从质量稳定性和安全性上控制药品、食品的质量。
申请公布号 CN103278591B 申请公布日期 2015.08.26
申请号 CN201310183164.8 申请日期 2013.05.16
申请人 江苏师范大学 发明人 孟庆华;闫艳;李广超;李亮;蒋德林
分类号 G01N30/86(2006.01)I 主分类号 G01N30/86(2006.01)I
代理机构 南京天华专利代理有限责任公司 32218 代理人 徐冬涛
主权项 一种色谱指纹图谱相似度计算方法,其特征是它包括以下步骤:(1)、采用标准样品建立对照指纹图谱,所述对照指纹图谱的向量为<img file="FDA0000735410870000016.GIF" wi="114" he="91" /><img file="FDA0000735410870000017.GIF" wi="446" he="81" />1≤i≤n,n为对照指纹图谱的色谱峰数即标准样品中的成分总数;y<sub>i</sub>代表标准样品中i成分的峰面积或峰高值;(2)、采用与建立对照指纹图谱相同的色谱条件建立待评价样品的样品指纹图谱,样品指纹图谱的向量为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>X</mi><mo>&RightArrow;</mo></mover><mo>=</mo><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>x</mi><mn>3</mn></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mo>...</mo><msub><mi>x</mi><mi>n</mi></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>a</mi><mn>3</mn></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>x</mi><mrow><mi>a</mi><mi>j</mi></mrow></msub><mo>,</mo><mo>...</mo><msub><mi>x</mi><mrow><mi>a</mi><mi>m</mi></mrow></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000735410870000011.GIF" wi="1015" he="103" /></maths>其中,1≤i≤n,n为样品指纹图谱中与对照指纹图谱共有峰的色谱峰数,x<sub>i</sub>代表待评价样品中i成分即样品指纹图谱中与对照指纹图谱共有峰的峰面积或峰高值,其中,对照指纹图谱中存在而在样品指纹图谱中不存在的成分,样品指纹图谱中该共有峰的取值为x<sub>i</sub>=0;1≤j≤m,m为样品指纹图谱中与对照指纹图谱非共有峰的色谱峰数,x<sub>aj</sub>代表待评价样品中j成分即样品指纹图谱中与对照指纹图谱非共有峰的峰面积或峰高值;(3)、将样品指纹图谱的每一个峰值与对照指纹图谱中对应峰值相比,样品的比值指纹图谱表示为n+1维向量:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><msup><mi>X</mi><mo>&prime;</mo></msup><mo>&RightArrow;</mo></mover><mo>=</mo><mrow><mo>(</mo><mrow><mfrac><msub><mi>x</mi><mn>1</mn></msub><msub><mi>y</mi><mn>1</mn></msub></mfrac><mo>,</mo><mfrac><msub><mi>x</mi><mn>2</mn></msub><msub><mi>y</mi><mn>2</mn></msub></mfrac><mo>,</mo><mfrac><msub><mi>x</mi><mn>3</mn></msub><msub><mi>y</mi><mn>3</mn></msub></mfrac><mo>,</mo><mo>...</mo><mo>,</mo><mfrac><msub><mi>x</mi><mi>i</mi></msub><msub><mi>y</mi><mi>i</mi></msub></mfrac><mo>,</mo><mo>...</mo><mo>,</mo><mfrac><msub><mi>x</mi><mi>n</mi></msub><msub><mi>y</mi><mi>n</mi></msub></mfrac><mo>,</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>x</mi><mrow><mi>a</mi><mi>j</mi></mrow></msub></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>y</mi><mi>i</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000735410870000012.GIF" wi="828" he="261" /></maths>对照比值指纹图谱表示为n+1维向量:<img file="FDA0000735410870000013.GIF" wi="574" he="99" />采用下述公式计算这两个比值向量的相似度S:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><mn>1</mn><mo>-</mo><msqrt><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msub><mi>&delta;</mi><mi>i</mi></msub><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><mfrac><msub><mi>x</mi><mi>i</mi></msub><msub><mi>y</mi><mi>i</mi></msub></mfrac></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>&delta;</mi><mi>a</mi></msub><msup><mrow><mo>(</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>x</mi><mrow><mi>a</mi><mi>j</mi></mrow></msub></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>y</mi><mi>i</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></mrow><msub><mi>n</mi><mi>t</mi></msub></mfrac></msqrt><mo>&times;</mo><mrow><mo>(</mo><mrow><mn>2</mn><mo>-</mo><mfrac><mrow><mo>|</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mo>|</mo></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>|</mo></mrow></mrow></mfrac></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000735410870000014.GIF" wi="1255" he="411" /></maths>其中:δ<sub>i</sub>、δ<sub>a</sub>分别是各共有峰、非共有峰总面积的权重系数,权重系数的取值范围δ<sub>i</sub>≥0、δ<sub>a</sub>≥0;<img file="FDA0000735410870000015.GIF" wi="130" he="142" />表示样品指纹图谱中样品指纹图谱与对照指纹图谱非共有峰的峰总面积或峰总高值,<img file="FDA0000735410870000021.GIF" wi="104" he="124" />作为求非共有峰比值的标准,表示对照指纹图谱峰总面积或峰总高值,n<sub>t</sub>表示实际的共有峰峰数,n<sub>t</sub>≤n,当对照指纹图谱存在的成分在样品指纹图谱中不存在时,样品指纹图谱中该共有峰的取值为x<sub>i</sub>=0。
地址 221116 江苏省徐州市铜山区上海路101号