CN105527361B - 一种高效微流电色谱相似度分析以及计算影响度的方法 - Google Patents
一种高效微流电色谱相似度分析以及计算影响度的方法 Download PDFInfo
- Publication number
- CN105527361B CN105527361B CN201410514223.XA CN201410514223A CN105527361B CN 105527361 B CN105527361 B CN 105527361B CN 201410514223 A CN201410514223 A CN 201410514223A CN 105527361 B CN105527361 B CN 105527361B
- Authority
- CN
- China
- Prior art keywords
- peak
- similarity
- sample
- response time
- alignment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明为一种高效微流电色谱相似度分析以及计算影响度的方法,该方法包括:通过对齐步骤将不同图谱上的波峰进行对齐,之后再通过相似度公式的计算,得到相关色谱图中的各检测物质的相似度数值,根据得到的数值来判断被检测物质与标准样品是否属于同一种物质。最后,通过影响程度公式来计算出各个峰在相似度计算时的影响程度。本发明可以有效的比对被检测物质与标准样品是否属于同一种物质,同时还能获取同一个被检测物质中的各个峰对相似度计算的影响程度,从而提高了高效微流电色谱的分析、辨别样品的能力。
Description
技术领域
本发明涉及高效微流电色谱样品分析领域,具体说涉及一种高效微流电色谱相似度分析以及计算影响度的方法。
背景技术
高效微流电色谱是一种非常重要的分离分析仪器,它的最大作用就是确定待测样品的名称,以及分析待测样品的成分组成。
为了能够让高效微流电色谱准确地分析出待测样品的名称,目前色谱分析领域常常在计算模块引入相似度的概念,如中国专利CN 104007212。由于仪器、待测样品、色谱柱以及溶剂的不确定性,在进行相似度计算之前,往往先要把所有待测样品得到的色谱数据进行一个对齐操作,保证所有样品在同一个条件下进行相似度计算,如中国专利CN101929989 A。
同时,为了能够了解到同一个样品中各个峰在相似度计算时的影响程度,需要引入影响度计算概念。影响度计算方法的引入,可以直观地让用户了解待测物质与标样之间存在差异的主要原因,帮助用户提升物质检测的效率和速度。
然而,目前现有技术还没有将对齐步骤,相似度计算,影响度计算等结合起来组成完整的技术方案,而影响度的计算法方式也比较简单,缺少高信度的计算方法。
发明内容
本发明的目的是为了提高高效微流电色谱在比对样品时的准确率以及帮助用户了解各个峰对相似度计算的影响程度,提高高效微流电色谱物质检测的准确率,提升高效微流电色谱物质检测的效率。
为了实现本发明的目的而采用的技术方案是:一种高效微流电色谱相似度分析以及计算影响度的方法,包括以下步骤:
S1,数据获取步骤,采集标准样品与未知样品的色谱信号数据;
S2,数据对齐步骤,通过对齐步骤将标准样品与未知样品的数据进行对齐;
S3,相似度计算步骤,根据设定的相似度计算公式,得到未知样品的相似度数值;
S4,影响度计算步骤,通过设定影响度计算公式,得到各个色谱峰对相似度计算的影响程度。
优选地,所述数据对齐步骤包括:
S21,从数据源选择数据列表;
S22,选择要对齐的峰,初始值为最后的峰,每一行只能选一个峰;
S23,以第一个数据为基准,进行对齐运算:
Ty=Tx*T1/T2;
T1:第一行对齐峰的响应时间;
T2:当前行对齐峰的响应时间;
Tx:对齐前峰的响应时间;
Ty:对齐后峰的响应时间;
S24,显示对齐后的数据列表。
优选地,所述相似度计算步骤包括:
S31.导入需要分析的所有样品数据,需要的项目是响应时间、峰名称、峰高/峰面积;
S32.根据时间窗/时间带的设置,遍历所有样品,计算出所有的响应时间,若某一样品在某个响应时间内没有峰的话,则峰高/峰面积设为0;
S33.以标准样品为基础,计算所有样品的相似度,预设相似度公式如下:
其中xi是第一个样品的峰高/峰面积,yi是另一个样品的峰高/峰面积,n是响应时间的个数。当其中某一项为0的时候xi/yi设为0,当两项都为0的时候xi/yi设为1。当xi>2yi时xi/yi设为0,反之当yi>2xi时xi/yi也设为0。
优选地,所述影响度计算步骤包括:
S41.取所有样品的相似度的平均值,大于该平均值的分为第一类,小于该值的分为另一类;
S42.分别计算出这两类的各个响应时间的峰高/峰面积的平均值;
S43.计算出各个峰在相似度计算时的影响程度,预设影响度公式如下:
其中xi是第一类样品的峰高/峰面积的平均值,yi是另一类样品的峰高/峰面积的平均值。
附图说明
图1是标准样品与被测样品进行相似度计算和影响度计算的步骤流程图。
具体实施方式
实施例:
本实施例的实验环境为:流动相:0.02mol/L磷酸二氢钠二乙胺缓冲溶液+甲醇(70:30);流速:1mL/min;检测波长:220nm;进样量:50μL;柱温:30℃;色谱柱:C18(l50mm×4.6mm,粒径5μm)。
本实施例具体步骤为:
1.从数据源选择数据列表;给出数据源和具体列表。同理,下同的每一步都要给出具体的数值。
2.选择要对齐的峰,初始值为最后的峰,每一行只能选一个峰;
3.以第一个数据为基准,进行对齐运算:
T1:第一行对齐峰的响应时间;
T2:当前行对齐峰的响应时间;
Tx:对齐前峰的响应时间;
Ty:对齐后峰的响应时间;
Ty=Tx*T1/T2;
4.显示对齐后的数据列表。
完成了标准样品与待测样品相关数据的对齐之后,即进入到相似度计算的步骤。相似度计算的步骤为:
1.导入需要分析的所有样品数据,需要的项目是响应时间、峰名称、峰高/峰面积;
2.根据时间窗/时间带的设置,遍历所有样品,计算出所有的响应时间,若某一样品在某个响应时间内没有峰的话,则峰高/峰面积设为0;
3.以标准样品为基础,计算所有样品的相似度,公式如下:
其中xi是第一个样品的峰高/峰面积,yi是另一个样品的峰高/峰面积,n是响应时间的个数。当其中某一项为0的时候xi/yi设为0,当两项都为0的时候xi/yi设为1。
当xi>2yi时xi/yi设为0,反之当yi>2xi时xi/yi也设为0。
完成相似度计算之后,如果相似度的数值越接近于1,说明被测样品与标准样品越有可能是同一种物质。反之,如果相似度的数值与1的差距越大,说明被测样品与标准样品越不可能是同一种物质。
完成以上这些步骤之后,即可开始进行影响度的计算,方法如下:
1.取所有样品的相似度的平均值,大于该平均值的分为第一类,小于该值的分为另一类;
2.分别计算出这两类的各个响应时间的峰高/峰面积的平均值;
3.计算出各个峰在相似度计算时的影响程度,公式如下:
其中xi是第一类样品的峰高/峰面积的平均值,yi是另一类样品的峰高/峰面积的平均值。
得到影响度计算的结果之后,如果结果越大,说明这个色谱峰对相似度计算的结果影响越大。反之,如果结果越小,说明这个色谱峰对相似度计算的结果影响越小。通过影响度的计算,用户可以很清楚地知道被测样品中哪一种物质对样品的测定起到了最大的影响。
通过对齐步骤、相似度计算步骤以及影响度计算步骤,用户能够很直观地了解标准样品与被测样品之间所存在的联系,以及是何种物质对样品的测定造成了最大的影响。从而提高了高效微流电色谱在比对样品时的准确率以及帮助用户了解了各个峰对相似度计算的影响程度。
以上所述仅为本发明的较佳实施方式,本发明的保护范围并不以上述实施方式为限,但凡本领域普通技术人员根据本发明所揭示内容所作的等效修饰或变化,皆应纳入权利要求书中记载的保护范围内。
Claims (3)
1.一种高效微流电色谱相似度分析以及计算影响度的方法,其特征在于包括以下步骤:
S1,数据获取步骤,采集标准样品与未知样品的色谱信号数据;
S2,数据对齐步骤,通过对齐步骤将标准样品与未知样品的数据进行对齐;
S3,相似度计算步骤,根据设定的相似度计算公式,得到未知样品的相似度数值;相似度公式如下:
其中xi是标准样品的峰高与峰面积的比值,yi是未知样品的峰高与峰面积的比值,n是响应时间的个数,当其中一项为0的时候xi/yi设为0,当两项都为0的时候xi/yi设为1,当xi>2yi时xi/yi设为0,反之当yi>2xi时xi/yi也设为0;
S4,影响度计算步骤,通过设定影响度计算公式,得到各个色谱峰对相似度计算的影响程度:
S41.取所有样品的相似度的平均值,大于该平均值的分为第一类,小于该值的分为另一类;
S42.分别计算出这两类的各个响应时间的峰高与峰面积的比值的平均值;
S43.计算出各个峰在相似度计算时的影响程度,公式如下:
其中xi是第一类样品的峰高与峰面积的比值的平均值,yi是另一类样品的峰高与峰面积的比值的平均值。
2.根据权利要求1所述的高效微流电色谱相似度分析以及计算影响度的方法,其特征在于:所述数据对齐步骤包括:
S21,从数据源选择数据列表;
S22,选择要对齐的峰,初始值为最后的峰,每一行只能选一个峰;
S23,以第一个数据为基准,进行对齐运算:
Ty=Tx*T1/T2;
T1:第一行对齐峰的响应时间;
T2:当前行对齐峰的响应时间;
Tx:对齐前峰的响应时间;
Ty:对齐后峰的响应时间;
S24,显示对齐后的数据列表。
3.根据权利要求1所述的高效微流电色谱相似度分析以及计算影响度的方法,其特征在于:所述相似度计算步骤包括:
S31.导入需要分析的所有样品数据,需要的项目包括响应时间、峰名称、峰高与峰面积的比值;
S32.根据时间窗与时间带的比值的设置,遍历所有样品,计算出所有的响应时间,若某一样品在某个响应时间内没有峰的话,则峰高与峰面积的比值设为0;
S33.以标准样品为基础,计算所有样品的相似度,公式如下:
其中xi是标准样品的峰高与峰面积的比值,yi是未知样品的峰高与峰面积的比值,n是响应时间的个数,当其中一项为0的时候xi/yi设为0,当两项都为0的时候xi/yi设为1,当xi>2yi时xi/yi设为0,反之当yi>2xi时xi/yi也设为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410514223.XA CN105527361B (zh) | 2014-09-29 | 2014-09-29 | 一种高效微流电色谱相似度分析以及计算影响度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410514223.XA CN105527361B (zh) | 2014-09-29 | 2014-09-29 | 一种高效微流电色谱相似度分析以及计算影响度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105527361A CN105527361A (zh) | 2016-04-27 |
CN105527361B true CN105527361B (zh) | 2017-06-20 |
Family
ID=55769715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410514223.XA Active CN105527361B (zh) | 2014-09-29 | 2014-09-29 | 一种高效微流电色谱相似度分析以及计算影响度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105527361B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110161161A (zh) * | 2019-07-01 | 2019-08-23 | 汕头出入境检验检疫局检验检疫技术中心 | 基于高效液相指纹图谱和聚类分析的红木鉴定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101603959A (zh) * | 2009-06-29 | 2009-12-16 | 上海应用技术学院 | 一种快速测定猪肉香精相似度的方法 |
CN101929989A (zh) * | 2009-06-24 | 2010-12-29 | 陈喆 | 一种代谢谱峰位置分辨与对齐的方法 |
CN104007212A (zh) * | 2014-06-23 | 2014-08-27 | 华中科技大学 | 一种基于局部相似性匹配的色谱特征峰波段提取与分析方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07218491A (ja) * | 1994-01-31 | 1995-08-18 | Shimadzu Corp | クロマトグラフ用検出装置 |
SE0302559D0 (sv) * | 2003-09-25 | 2003-09-25 | Astrazeneca Ab | Method |
US7983874B2 (en) * | 2008-06-10 | 2011-07-19 | National University Of Ireland, Galway | Similarity index: a rapid classification method for multivariate data arrays |
-
2014
- 2014-09-29 CN CN201410514223.XA patent/CN105527361B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101929989A (zh) * | 2009-06-24 | 2010-12-29 | 陈喆 | 一种代谢谱峰位置分辨与对齐的方法 |
CN101603959A (zh) * | 2009-06-29 | 2009-12-16 | 上海应用技术学院 | 一种快速测定猪肉香精相似度的方法 |
CN104007212A (zh) * | 2014-06-23 | 2014-08-27 | 华中科技大学 | 一种基于局部相似性匹配的色谱特征峰波段提取与分析方法 |
Non-Patent Citations (3)
Title |
---|
中药指纹图谱相似度计算的规范化研究;乔善磊等;《第二军医大学学报》;20041031;第25卷(第10期);全文 * |
化学指纹图谱的相似性测度及其评价方法;程翼宇;《化学学报》;20021231;第60卷(第11期);全文 * |
相似系统理论用于中药色谱指纹图谱的相似度评价;刘永锁;《色谱》;20050331;第23卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105527361A (zh) | 2016-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2013128285A3 (en) | Improved analytical methods for analyzing and determining impurities in dianhydrogalactitol | |
MX2018001134A (es) | Aparato y metodo de deteccion de gas. | |
CN103743844B (zh) | 甲醇中乙醇含量的测定方法 | |
CN105518456A (zh) | 色谱数据处理装置 | |
CN103308634B (zh) | 含气量测量装置及方法 | |
MX2016006457A (es) | Analisis de cromatografia liquida de alto rendimiento (hplc) de impurezas en dianhidrogalactitol. | |
CN102680499A (zh) | 基于mems技术的盐度检测传感器及其检测方法 | |
CN105527361B (zh) | 一种高效微流电色谱相似度分析以及计算影响度的方法 | |
CN103105362A (zh) | 一种多量程水质在线检测仪 | |
CN203164406U (zh) | 基于gtem室的超高频传感器灵敏度测试装置 | |
CN204008578U (zh) | 一种气相色谱柱箱的温度保护系统 | |
CN204594947U (zh) | 一种制备离子色谱仪 | |
CN103852642A (zh) | 一种检测微量固体导电性的方法 | |
CN202676589U (zh) | 一种便携式矿井水源快速检测识别设备 | |
CN104237427A (zh) | 一种液相色谱流动相预警装置 | |
CN103884789A (zh) | 一种快速测定灵芝产品中多糖肽的方法 | |
CN105044263B (zh) | 一种4-乙氧基-2,3-二氟苯酚的气相色谱检验方法 | |
CN104391055B (zh) | 一种检测、分析粮食中苯并[a]芘的方法 | |
Tudyka et al. | Very low cost multichannel analyzer with some additional features | |
CN203629601U (zh) | 一种火灾探测器智能测试装置 | |
MX2018006425A (es) | Metodos para analisis y resolucion de preparaciones de dianhidrogalactitol y derivados o analogos del mismo. | |
CN104597159A (zh) | 基于色谱法的污染预警方法及装置 | |
CN105628797B (zh) | 艾考糊精原料药分子量及其分布测定方法 | |
CN105865973A (zh) | 一种流体密度测量装置 | |
CN105203691A (zh) | 共聚酯热熔胶网膜中特种含磷阻燃剂的分析测定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |