1.本发明涉及医疗技术领域,具体涉及一种特定蛋白检测方法,以及一种特定蛋白检测装置。
背景技术:2.c反应蛋白(crp)是肝脏合成的一种急性炎症正时相反应蛋白,也是一种特定蛋白。正常人血中的crp浓度很低,在机体突遇紧张、组织创伤和各种炎症刺激时合成快速增加,并从肝细胞中分泌入血液,在感染发生后12
‑
18小时即可检测到高水平的crp。在感染发生后的12
‑
14天升高的crp可降至基线水平。因此多年来—直是评价炎症性疾病的指标之一,并且升高幅度与感染的程度呈相关。crp作为诊断细菌感染的重要标志物之一,已广泛应用于临床。crp还是一个评估心脏病发生率、复发率、死亡率的临床重要指标。近年来的研究发现,炎症在动脉粥样硬化及肿瘤的发生、发展过程中起重要作用。鉴于血清crp的重要作用,其测量的准确性受到了广泛的关注。
3.检测crp的常用方法多种多样,其中包括散射比浊法、透射比浊法,放射免疫测定、化学发光法、elisa法及床旁crp检测(poct)等。目前临床实验室测定血清中crp的方法主要是免疫浊度法,包括乳胶增强透射比浊法和速率散射比浊法,这两种方法主要用于自动化分析系统,速率散射比浊法多用于免疫检测领域的封闭检测系统,乳胶增强透射比浊法多用于生化检测领域的开放检测系统。
4.目前基于透射比浊法对样本的特定蛋白浓度进行检测时,首先得到该样本的特定蛋白反应曲线,然后求得特定蛋白反应曲线的一个曲线特征,接着将该曲线特征代入定标曲线,则可求出样本的特定蛋白浓度。上述方式如果在样本采集(也就是检测样本的特定蛋白浓度)时,如果某一段时间信号存在干扰,特定蛋白反应曲线会出现抖动或鼓包等异常情况,如果选取计算的曲线特征的点正好位于特定蛋白反应曲线抖动或鼓包的位置,就会导致计算结果不准,造成假阴性或假阳性的结果。
技术实现要素:5.根据第一方面,一种实施例中提供一种特定蛋白浓度的检测方法,包括:
6.在预设的采样总时间内基于透射比浊法对待测样本进行特定蛋白检测,以获取待测样本的特定蛋白反应曲线,所述待测样本的特定蛋白反应曲线用于表示采样总时间内由透射比浊法得到的电压的变化;
7.获取所述待测样本的特定蛋白反应曲线的一阶导数的最大值;
8.根据所述一阶导数的最大值与预设阈值的关系,确定所述待测样本的特定蛋白反应曲线是否存在波动,如果不存在波动,在所述采样总时间内选取提取时间段,如果存在波动,则获取所述待测样本的特定蛋白反应曲线的二阶导数,并根据所述二阶导数确定采样总时间内中的采样有效时间,在所述采样有效时间内选取提取时间段;
9.提取所述待测样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征;
10.根据预先基于透射比浊法得到的多个已知浓度样本的特定蛋白反应曲线,提取各已知浓度样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征;
11.建立各已知浓度样本的特定蛋白反应曲线的曲线特征与其对应的特定蛋白浓度间的定标函数关系式;
12.将所述待测样本的曲线特征输入所述定标函数关系式中,得到所述待测样本的特定蛋白浓度。
13.根据第二方面,一种实施例中提供一种计算机可读存储介质,包括程序,所述程序能够被处理器执行以实现第一方面中的方法。
14.根据第三方面,一种实施例中提供一种特定蛋白浓度的检测装置,包括:
15.反应容器,用于容纳待测样本;
16.光源,用于向所述反应容器内的待测样本提供激光;
17.光信号接收器,用于采集激光经待测样本形成的透射光,并将透射光的光信号转换为电信号;
18.存储装置,用于存储各已知浓度样本的特定蛋白浓度和预先基于透射比浊法得到的对应的特定蛋白反应曲线;
19.数据处理装置,与所述光信号接收器和存储装置分别信号连接,用于:
20.根据所述电信号获取待测样本的特定蛋白反应曲线,所述特定蛋白反应曲线用于表示采样总时间内得到的电压的变化;
21.获取所述待测样本的特定蛋白反应曲线的一阶导数的最大值;
22.根据所述一阶导数的最大值与预设阈值的关系,确定所述待测样本的特定蛋白反应曲线是否存在波动,如果不存在波动,在所述采样总时间内选取提取时间段,如果存在波动,则获取所述待测样本的特定蛋白反应曲线的二阶导数,并根据所述二阶导数确定采样总时间内中的采样有效时间,在所述采样有效时间内选取提取时间段;
23.提取所述待测样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征;
24.根据预先基于透射比浊法得到的多个已知浓度样本的特定蛋白反应曲线,提取各已知浓度样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征;
25.建立各已知浓度样本的特定蛋白反应曲线的曲线特征与其对应的特定蛋白浓度间的定标函数关系式;
26.将所述待测样本的曲线特征输入所述定标函数关系式中,得到所述待测样本的特定蛋白浓度。
27.上述实施例中,首先判断特定蛋白反应曲线是否存在波动(例如抖动或鼓包的),如果不存在波动,直接获取曲线特征,如果存在波动,又识别出特定蛋白反应曲线不存在波动的部分,并由此得到待测样本的曲线特征。该方式降低了样本采集时对特定蛋白反应曲线完整性的依赖,在特定蛋白反应曲线发生抖动或鼓包等异常时也能计算出正确的结果,从而提高了检测的准确性和可靠性。
附图说明
28.图1为一种实施例的特定蛋白浓度检测装置的结构示意图;
29.图2为一种实施例的特定蛋白反应曲线的示意图;
30.图3为另一种实施例的特定蛋白反应曲线的示意图;
31.图4为又一种实施例的特定蛋白反应曲线的示意图;
32.图5为一种实施例的特定蛋白浓度检测方法的流程图。
33.10、反应容器;
34.20、光源;
35.30、光信号接收器;
36.40、存储装置;
37.50、数据处理装置。
具体实施方式
38.下面通过具体实施方式结合附图对本发明作进一步详细说明。其中不同实施方式中类似元件采用了相关联的类似的元件标号。在以下的实施方式中,很多细节描述是为了使得本申请能被更好的理解。然而,本领域技术人员可以毫不费力的认识到,其中部分特征在不同情况下是可以省略的,或者可以由其他元件、材料、方法所替代。在某些情况下,本申请相关的一些操作并没有在说明书中显示或者描述,这是为了避免本申请的核心部分被过多的描述所淹没,而对于本领域技术人员而言,详细描述这些相关操作并不是必要的,他们根据说明书中的描述以及本领域的一般技术知识即可完整了解相关操作。
39.另外,说明书中所描述的特点、操作或者特征可以以任意适当的方式结合形成各种实施方式。同时,方法描述中的各步骤或者动作也可以按照本领域技术人员所能显而易见的方式进行顺序调换或调整。因此,说明书和附图中的各种顺序只是为了清楚描述某一个实施例,并不意味着是必须的顺序,除非另有说明其中某个顺序是必须遵循的。
40.本文中为部件所编序号本身,例如“第一”、“第二”等,仅用于区分所描述的对象,不具有任何顺序或技术含义。而本申请所说“连接”、“联接”,如无特别说明,均包括直接和间接连接(联接)。
41.请参考图1,为透射比浊法的原理示意图,包括光源20、反应容器10和光信号接收器30,其中,反应容器10中装有反应液,光源20垂直射入反应容器10,穿过反应容器10壁照射到反应液中的微球上,发生透射,再通过反应容器10壁进入空气,最后到光信号接收器30。
42.特定蛋白反应检测原理是依据抗原、抗体在特定的电解质溶液中反应,快速形成免疫复合物微粒,使反应液出现浊度,随着时间的推移,复合物微粒聚合得就会越大,浊度也逐渐变大。在透射比浊法中,光信号接收器30接收到光信号的强度,随时间的推移越来越小。c反应蛋白是特定蛋白检测项目的一种。
43.请参考图2,为一种特定蛋白反应曲线示意图,坐标横轴是采样时间,坐标纵轴是电压,整个曲线是光信号接收器30接收的光信号强度随时间的变化曲线,通常用数学函数公式表示特定蛋白反应曲线,其包括:
44.v=f(t)。
45.其中,0≤t,v是电压,t是采样时间,整个采样过程的持续时间在本发明中定义为采样总时间,特定蛋白反应曲线为一减函数。
46.请参照图1,本发明提供了一种特定蛋白浓度的检测装置,包括反应容器10、光源
20、光信号接收器30、存储装置40和数据处理装置50。
47.反应容器10用于容纳待测样本,例如可以是透明的反应容器10。
48.光源20作为透射比浊法使用的光源20,其用于向反应容器10内的待测样本提供激光。
49.光信号接收器30设置在激光经待测样本形成的透射光的光路上,其用于采集透射光,并将透射光的光信号转换为电信号。
50.存储装置40用于存储各已知浓度样本的特定蛋白浓度和预先基于透射比浊法得到的对应的特定蛋白反应曲线。
51.一些实施例中,已知浓度样本的数目为n个,即具有n个按浓度梯度排列的已知浓度样本b1,b2,
…
,b
i
,
…
,b
n
;其中,第i个已知浓度样本的特定蛋白浓度为c
i
,且0<c1<c2<
…
<c
i
<
…
<c
n
,1≤i≤n。
52.第i个已知浓度样本的定蛋白反应曲线的公式为:
53.v
bi
=f
i
(t
b
)。
54.其中,0≤t
b
≤t
b
,t
b
∈实数,t
b
为已知浓度样本的采样时间,v
bi
为第i个已知浓度样本采样时获取的电压,t
b
为已知浓度样本的采样总时间的结束时间点。
55.一些实施例中,对于任一已知浓度样本的特定蛋白反应曲线,还判断已知浓度样本的特定蛋白反应曲线是否存在波动,如果存在波动,则对存在波动的特定蛋白反应曲线的已知浓度样本重新检测,以再次获取相应的特定蛋白反应曲线,直到每个已知浓度样本的特定蛋白反应曲线均不存在波动。对于第i个已知浓度样本bi,判断其特定蛋白反应曲线是否存在波动的方式,可以是:获取第i个已知浓度样本的特定蛋白反应曲线的一阶导数f
i’(t
b
),然后获取一阶导数f
i’(t
b
)的最大值max(f
i
'(t
b
)),如果max(f
i
'(t
b
))不大于零,则第i个已知浓度样本的特定蛋白反应曲线不存在波动,否则,第i个已知浓度样本的特定蛋白反应曲线存在波动,那么就再次对该已知浓度样本进行检测,直到该已知浓度样本的特定蛋白反应曲线的一阶导数的最大值不大于0。
56.数据处理装置50与光信号接收器30和存储装置40分别信号连接,其用于根据电信号获取待测样本的特定蛋白反应曲线,特定蛋白反应曲线用于表示采样总时间内得到的电压的变化。
57.获取到的待测样本的特定蛋白反应曲线的公式,包括:
58.v
a
=f(t
a
)。
59.其中,0≤t
a
≤t
a
,t
a
∈实数,t
a
为待测样本的采样时间,v
a
为采样时获取的电压,t
a
为采样总时间的结束时间点。
60.而后获取待测样本的特定蛋白反应曲线的一阶导数f’(t
a
)的最大值max(f'(t
a
)),根据一阶导数的最大值与预设阈值的关系,确定待测样本的特定蛋白反应曲线是否存在波动。一些实施例中,该预设阈值为零,如果max(f'(t
a
))小于等于零,代表待测样本的特定蛋白反应曲线为一减函数,待测样本的特定蛋白反应曲线不存在波动,如果max(f'(t
a
))大于零,待测样本的特定蛋白反应曲线存在波动。
61.如果待测样本的特定蛋白反应曲线不存在波动,在采样总时间内选取提取时间段,如果存在波动,则获取待测样本的特定蛋白反应曲线的二阶导数,并根据二阶导数确定采样总时间内中的采样有效时间,在采样有效时间内选取提取时间段。采样有效时间就是
特定蛋白反应曲线不存在波动的时间段,具体的,数据处理装置50获取待测样本的特定蛋白反应曲线的二阶导数f”(t
a
),根据二阶导数确定采样有效时间。采样有效时间内待测样本的特定蛋白反应曲线需要满足f”(t
a
)大于零。并且,0≤t1≤t2≤t
a
,t1为采样有效时间的起始时间点,t2为采样有效时间的结束时间点。如图3(箭头指向特定蛋白反应曲线处为波动处)所示,通常,满足上述条件的时间段只有一段。如果该时间段有两段或两段以上,选取时间长度最大的时间段作为采样有效时间,例如如图4所示,一满足上述条件的时间段的长度(t'2‑
t'1)较小,就不作为采样有效时间。
62.获取到提取时间段后,数据处理装置50提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
63.一些实施例中,如图2所示,待测样本的特定蛋白反应曲线不存在波动,提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征,包括:
64.计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|d
a1
|,电压差的绝对值|d
a1
|的计算公式包括:
65.|d
a1
|=|f(t
a2
)
‑
f(t
a1
)|。
66.其中,t
a1
为提取时间段的起始时间点,t
a2
为提取时间段的结束时间点,0≤t
a1
<t
a2
≤t
a
,t
a
为采样总时间的结束时间点,将电压差的绝对值|d
a1
|作为待测样本的特定蛋白反应曲线的曲线特征。
67.在另一些实施例中,如图2所示不存在波动,计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积s
a1
,面积s
a1
的计算公式包括:
[0068][0069]
其中,t
a1
为提取时间段的起始时间点,t
a2
为提取时间段的结束时间点,0≤t
a1
<t
a2
≤t
a
,t
a
为采样总时间的结束时间点,将面积s
a1
作为待测样本的特定蛋白反应曲线的曲线特征。
[0070]
一些实施例中,如图3或图4(箭头指向特定蛋白反应曲线处为波动处)所示,待测样本的特定蛋白反应曲线存在波动,提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征,包括:
[0071]
计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|d
a2
|,电压差的绝对值|d
a2
|的计算公式包括:
[0072]
|d
a2
|=|f(t
a2
)
‑
f(t
a1
)|。
[0073]
其中,t
a1
为提取时间段的起始时间点,t
a2
为提取时间段的结束时间点,t1≤t
a1
<t
a2
≤t2,t1为采样有效时间的起始时间点,t2为采样有效时间的结束时间点,将电压差的绝对值|d
a2
|作为待测样本的特定蛋白反应曲线的曲线特征。
[0074]
在另一些实施例中,如图3或图4(箭头指向特定蛋白反应曲线处为波动处)所示存在波动,计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积s
a2
,面积s
a2
的计算公式包括:
[0075][0076]
其中,t
a1
为提取时间段的起始时间点,t
a2
为提取时间段的结束时间点,t1≤t
a1
<
t
a2
≤t2,t1为采样有效时间的起始时间点,t2为采样有效时间的结束时间点,将面积s
a2
作为待测样本的特定蛋白反应曲线的曲线特征。
[0077]
而后,根据存储装置40中存储的n个已知浓度样本的特定蛋白反应曲线,获取各已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征r1,r2,
…
,r
i
,
…
,r
n
,1≤i≤n,r
i
为第i个已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
[0078]
一些实施例中,对于第i个已知浓度样本bi,计算已知浓度样本的特定蛋白反应曲线上不同两点间的电压差的绝对值,以获取电压差的绝对值|d
bi
|,电压差的绝对值|d
bi
|的计算公式包括:
[0079]
|d
bi
|=|f
i
(t
b2
)
‑
f
i
(t
b1
)|。
[0080]
其中,t
b1
和t
b2
是第i个已知浓度样本的特定蛋白反应曲线上两点的采样时间,0≤t
b1
<t
b2
≤t
b
,|d
bi
|为电压差的绝对值,t
b
为已知浓度样本的采样总时间的结束时间点,如果待测样本的曲线特征为|d
a1
|或|d
a2
|,则取t
b2
=t
a2
,t
b1
=t
a1
,将电压差的绝对值|d
bi
|作为第i个已知浓度样本的特定蛋白反应曲线的曲线特征r
i
。
[0081]
一些实施例中,对于第i个已知浓度样本bi,计算已知浓度样本的特定蛋白反应曲线上不同两点间的面积,以获取面积s
bi
,面积s
bi
的计算公式包括:
[0082][0083]
其中,s
bi
为已知浓度样本的特定蛋白反应曲线上不同两点间的面积,t
b1
和t
b2
是已知浓度样本的特定蛋白反应曲线上两点的采样时间,0≤t
b1
<t
b2
≤t
b
,t
b
为已知浓度样本采样总时间的结束时间点,如果待测样本的曲线特征为s
a1
或s
a2
,则取t
b2
=t
a2
,t
b1
=t
a1
,将面积s
bi
作为第i个已知浓度样本的特定蛋白反应曲线的曲线特征r
i
。
[0084]
提取到每个已知浓度样本的特定蛋白反应曲线的特征值后,将每个已知浓度样本的特定蛋白反应曲线的特征值与其对应的特定蛋白浓度建立数学函数对应关系,其中,将第i个已知浓度样本的特定蛋白反应曲线的特征值r
i
与其对应的特定蛋白浓度c
i
建立数学函数对应关系,1≤i≤n,进而获取定标函数关系式:r=f(c),其中,c为已知浓度样本的特定蛋白浓度,r为已知浓度样本的特定蛋白反应曲线的特征值。
[0085]
得到待测样本的特定蛋白反应曲线的曲线特征后,数据处理装置50将待测样本的曲线特征输入定标函数关系式中,从而得到一个特定蛋白浓度,将该特定蛋白浓度作为待测样本的特定蛋白浓度。
[0086]
本发明还提供了一种特定蛋白浓度的检测方法,如图5所示,包括步骤:
[0087]
步骤100、在预设的采样总时间内基于透射比浊法对待测样本进行特定蛋白检测,以获取待测样本的特定蛋白反应曲线,待测样本的特定蛋白反应曲线用于表示采样总时间内由透射比浊法得到的电压的变化。
[0088]
获取到的待测样本的特定蛋白反应曲线的公式,包括:
[0089]
v
a
=f(t
a
)。
[0090]
其中,0≤t
a
≤t
a
,t
a
∈实数,t
a
为待测样本的采样时间,v
a
为采样时获取的电压,t
a
为采样总时间的结束时间点。
[0091]
步骤200、获取待测样本的特定蛋白反应曲线的一阶导数的最大值。即获取待测样本的特定蛋白反应曲线的一阶导数f’(t
a
)的最大值max(f'(t
a
))。
[0092]
步骤300、根据一阶导数的最大值与预设阈值的关系,确定待测样本的特定蛋白反应曲线是否存在波动,如果不存在波动,执行步骤400,如果存在波动,执行步骤500。
[0093]
一些实施例中,该预设阈值为零,如果max(f'(t
a
))小于等于零,代表待测样本的特定蛋白反应曲线为一减函数,待测样本的特定蛋白反应曲线不存在波动,如果max(f'(t
a
))大于零,待测样本的特定蛋白反应曲线存在波动。
[0094]
步骤400、在采样总时间内选取提取时间段后,执行步骤700。
[0095]
步骤500、获取待测样本的特定蛋白反应曲线的二阶导数,并根据二阶导数确定采样总时间内中的采样有效时间后,执行步骤600。
[0096]
采样有效时间内待测样本的特定蛋白反应曲线需要满足f”(t
a
)大于零。并且,0≤t1≤t2≤t
a
,t1为采样有效时间的起始时间点,t2为采样有效时间的结束时间点。通常,满足上述条件的时间段只有一段,如果该时间段有两段或两段以上,则选取时间长度最大的时间段作为采样有效时间,例如如图4所示,一满足上述条件的时间段的长度(t'2‑
t'1)较小,就不作为采样有效时间。
[0097]
步骤600、在采样有效时间内选取提取时间段后,执行步骤700。
[0098]
步骤700、提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
[0099]
一些实施例中,如图2所示,待测样本的特定蛋白反应曲线不存在波动,提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征,包括:
[0100]
计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|d
a1
|,电压差的绝对值|d
a1
|的计算公式包括:
[0101]
|d
a1
|=|f(t
a2
)
‑
f(t
a1
)|。
[0102]
其中,t
a1
为提取时间段的起始时间点,t
a2
为提取时间段的结束时间点,0≤t
a1
<t
a2
≤t
a
,t
a
为采样总时间的结束时间点,将电压差的绝对值|d
a1
|作为待测样本的特定蛋白反应曲线的曲线特征。
[0103]
在另一些实施例中,如图2所示不存在波动,计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积s
a1
,面积s
a1
的计算公式包括:
[0104][0105]
其中,t
a1
为提取时间段的起始时间点,t
a2
为提取时间段的结束时间点,0≤t
a1
<t
a2
≤t
a
,t
a
为采样总时间的结束时间点,将面积s
a1
作为待测样本的特定蛋白反应曲线的曲线特征。
[0106]
一些实施例中,如图3或图4(箭头指向特定蛋白反应曲线处为波动处)所示,待测样本的特定蛋白反应曲线存在波动,提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征,包括:
[0107]
计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|d
a2
|,电压差的绝对值|d
a2
|的计算公式包括:
[0108]
|d
a2
|=|f(t
a2
)
‑
f(t
a1
)|。
[0109]
其中,t
a1
为提取时间段的起始时间点,t
a2
为提取时间段的结束时间点,t1≤t
a1
<t
a2
≤t2,t1为采样有效时间的起始时间点,t2为采样有效时间的结束时间点,将电压差的绝对值|d
a2
|作为待测样本的特定蛋白反应曲线的曲线特征。
[0110]
在另一些实施例中,如图3或图4(箭头指向特定蛋白反应曲线处为波动处)所示存在波动,计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积s
a2
,面积s
a2
的计算公式包括:
[0111][0112]
其中,t
a1
为提取时间段的起始时间点,t
a2
为提取时间段的结束时间点,t1≤t
a1
<t
a2
≤t2,t1为采样有效时间的起始时间点,t2为采样有效时间的结束时间点,将面积s
a2
作为待测样本的特定蛋白反应曲线的曲线特征。
[0113]
步骤800、根据预先基于透射比浊法得到的多个已知浓度样本的特定蛋白反应曲线,提取各已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
[0114]
一些实施例中,已知浓度样本的数目为n个,即具有n个按浓度梯度排列的已知浓度样本b1,b2,
…
,b
i
,
…
,b
n
;其中,第i个已知浓度样本的特定蛋白浓度为c
i
,且0<c1<c2<
…
<c
i
<
…
<c
n
,1≤i≤n。
[0115]
第i个已知浓度样本的定蛋白反应曲线的公式为:
[0116]
v
bi
=f
i
(t
b
)。
[0117]
其中,0≤t
b
≤t
b
,t
b
∈实数,t
b
为已知浓度样本的采样时间,v
bi
为第i个已知浓度样本采样时获取的电压,t
b
为已知浓度样本的采样总时间的结束时间点。
[0118]
一些实施例中,对于任一已知浓度样本的特定蛋白反应曲线,还判断已知浓度样本的特定蛋白反应曲线是否存在波动,如果存在波动,则对存在波动的特定蛋白反应曲线的已知浓度样本重新检测,以再次获取相应的特定蛋白反应曲线,直到每个已知浓度样本的特定蛋白反应曲线均不存在波动。对于第i个已知浓度样本bi,判断其特定蛋白反应曲线是否存在波动的方式,可以是:获取第i个已知浓度样本的特定蛋白反应曲线的一阶导数f
i’(t
b
),然后获取一阶导数f
i’(t
b
)的最大值max(f
i
'(t
b
)),如果max(f
i
'(t
b
))不大于零,则第i个已知浓度样本的特定蛋白反应曲线不存在波动,否则,第i个已知浓度样本的特定蛋白反应曲线存在波动,那么就再次对该已知浓度样本进行检测,直到该已知浓度样本的特定蛋白反应曲线的一阶导数的最大值不大于0。
[0119]
基于上述预先获取的各已知浓度样本的特定蛋白曲线,步骤800中获取各已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征r1,r2,
…
,r
i
,
…
,r
n
,1≤i≤n,r
i
为第i个已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
[0120]
一些实施例中,对于第i个已知浓度样本bi,计算已知浓度样本的特定蛋白反应曲线上不同两点间的电压差的绝对值,以获取电压差的绝对值|d
bi
|,电压差的绝对值|d
bi
|的计算公式包括:
[0121]
|d
bi
|=|f
i
(t
b2
)
‑
f
i
(t
b1
)|。
[0122]
其中,t
b1
和t
b2
是第i个已知浓度样本的特定蛋白反应曲线上两点的采样时间,0≤t
b1
<t
b2
≤t
b
,|d
bi
|为电压差的绝对值,t
b
为已知浓度样本的采样总时间的结束时间点,如果待测样本的曲线特征为d
a1
或d
a2
,则取t
b2
=t
a2
,t
b1
=t
a1
,将电压差的绝对值|d
bi
|作为第i个已知浓度样本的特定蛋白反应曲线的曲线特征r
i
。
[0123]
一些实施例中,对于第i个已知浓度样本bi,计算已知浓度样本的特定蛋白反应曲线上不同两点间的面积,以获取面积s
bi
,面积s
bi
的计算公式包括:
[0124][0125]
其中,s
bi
为已知浓度样本的特定蛋白反应曲线上不同两点间的面积,t
b1
和t
b2
是已知浓度样本的特定蛋白反应曲线上两点的采样时间,0≤t
b1
<t
b2
≤t
b
,t
b
为已知浓度样本采样总时间的结束时间点,如果待测样本的曲线特征为s
a1
或s
a2
,则取t
b2
=t
a2
,t
b1
=t
a1
,将面积s
bi
作为第i个已知浓度样本的特定蛋白反应曲线的曲线特征r
i
。
[0126]
步骤900、建立各已知浓度样本的特定蛋白反应曲线的曲线特征与其对应的特定蛋白浓度间的定标函数关系式。
[0127]
将第i个已知浓度样本的特定蛋白反应曲线的特征值r
i
与其对应的特定蛋白浓度c
i
建立数学函数对应关系,1≤i≤n,进而获取定标函数关系式:r=f(c),其中,c为已知浓度样本的特定蛋白浓度,r为已知浓度样本的特定蛋白反应曲线的特征值。
[0128]
步骤1000、将待测样本的曲线特征输入定标函数关系式中,得到待测样本的特定蛋白浓度。
[0129]
上述实施例中,首先判断特定蛋白反应曲线是否存在波动(例如抖动或鼓包的),如果不存在波动,直接获取曲线特征,如果存在波动,又识别出特定蛋白反应曲线不存在波动的部分,并由此得到待测样本的曲线特征。该方式降低了样本采集时对特定蛋白反应曲线完整性的依赖,在特定蛋白反应曲线发生抖动或鼓包等异常时也能计算出正确的结果,从而提高了检测的准确性和可靠性。
[0130]
本领域技术人员可以理解,上述实施方式中各种方法的全部或部分功能可以通过硬件的方式实现,也可以通过计算机程序的方式实现。当上述实施方式中全部或部分功能通过计算机程序的方式实现时,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:只读存储器、随机存储器、磁盘、光盘、硬盘等,通过计算机执行该程序以实现上述功能。例如,将程序存储在设备的存储器中,当通过处理器执行存储器中程序,即可实现上述全部或部分功能。另外,当上述实施方式中全部或部分功能通过计算机程序的方式实现时,该程序也可以存储在服务器、另一计算机、磁盘、光盘、闪存盘或移动硬盘等存储介质中,通过下载或复制保存到本地设备的存储器中,或对本地设备的系统进行版本更新,当通过处理器执行存储器中的程序时,即可实现上述实施方式中全部或部分功能。
[0131]
以上应用了具体个例对本发明进行阐述,只是用于帮助理解本发明,并不用以限制本发明。对于本发明所属技术领域的技术人员,依据本发明的思想,还可以做出若干简单推演、变形或替换。