[发明专利]熔解曲线重叠峰的分离方法、装置和电子设备

权利要求书

1.一种熔解曲线重叠峰的分离方法，其特征在于，包括：

基于预先获取的PCR反应板各孔熔解曲线的熔点峰位置，利用Savitzky-Golay求导方法确定各孔熔解曲线的二阶导数曲线；

根据所述二阶导数曲线和所述熔点峰位置，确定原始熔解曲线的初始重叠峰位置；

对所述初始重叠峰位置进行叠加拟合，生成两个熔点峰曲线；

根据两个所述熔点峰曲线确定满足误差范围的目标熔点峰，并对所述目标熔点峰进行分离。

2.根据权利要求1所述的熔解曲线重叠峰的分离方法，其特征在于，所述PCR反应板各孔熔解曲线的熔点峰位置，对应所述熔点峰的温度(Tm)值；

利用Savitzky-Golay求导方法确定各孔熔解曲线的二阶导数曲线的步骤，包括：

根据Savitzky-Golay求导方法，利用预先设定的第一求导参数计算各孔熔解曲线二阶导数曲线；

其中，所述第一求导参数包括：具有最高决定系数的窗口大小和多项式阶次。

3.根据权利要求1所述的熔解曲线重叠峰的分离方法，其特征在于，根据所述二阶导数曲线和所述熔点峰位置，确定目标熔解曲线的初始重叠峰位置的步骤，包括：

根据预先给定的第一范围，确定各熔点峰对应的第一区域，并进行归一化处理；所述第一区域包括在所述第一范围内选取的所述熔解曲线的二阶导数区域；

确定所述第一区域内的所述二阶导数曲线的极小值点；

根据所述极小值点的相对幅值，确定满足预先给定的第一阈值的所述极小值点为初始重叠峰位置；预先给定第一阈值的条件包括：在所述第一区域内，所述二阶导数曲线的极大值点幅值相差不超过两倍。

4.根据权利要求3所述的熔解曲线重叠峰的分离方法，其特征在于，预先给定第一范围的步骤包括：

以所述熔点峰位置对应的所述熔点峰的温度(Tm)值为基点，以标准熔点峰分辨率为间隔，在所述基点的左右两端分别确定所述熔点峰对应的波谷；所述标准熔点峰分辨率为最小熔点峰分辨率的四倍；其中，所述最小熔点峰分辨率为2℃；

如果在所述基点的左右两端均不存在所述波谷，则确定所述标准熔点峰分辨率为第一中间值；所述第一范围的最小值为所述熔点峰的温度(Tm)值减去所述第一中间值；所述第一范围的最大值为所述熔点峰的温度(Tm)值加上所述第一中间值；

如果在所述基点的左右两端均存在所述波谷，则确定所述基点与左右两端所述波谷的距离分别为T_左和T_右；比较T_左和所述标准熔点峰分辨率的大小，确定最小值为第二中间值；比较T_右和所述分辨率的大小，确定最小值为第三中间值；

则所述第一范围的最小值为所述熔点峰的温度(Tm)值减去所述第二中间值；所述第一范围的最大值为所述熔点峰的温度(Tm)值加上所述第三中间值。

5.根据权利要求4所述的熔解曲线重叠峰的分离方法，其特征在于，对所述初始重叠峰位置进行叠加拟合，生成两个熔点峰曲线的步骤，包括：

基于初始重叠峰位置设置初始函数的初始参数；

利用Levenberg-Marquardt算法，基于所述初始函数对各熔点峰对应的所述第一区域内的所述二阶导数曲线进行叠加拟合，生成两个待优化函数；

根据所述初始参数对所述两个待优化函数进行优化，生成两个目标函数；其中，所述目标函数包括目标参数；两个所述目标函数分别对应两个所述熔点峰曲线。

6.根据权利要求5所述的熔解曲线重叠峰的分离方法，其特征在于，

所述目标函数为三参数洛伦兹函数：g(v)＝a·c/[c+(v-b)²]；

其中：v用于表示温度，a用于表示系数，b用于表示优化之后的重叠峰的位置，c用于表示方差。