[發(fā)明專利]熔解曲線重疊峰的分離方法、裝置和電子設備

在審

• 申請?zhí)枺?/td>: 202210002067.3
• 申請日: 2022-01-04
• 公開（公告）號: CN114386461A
• 公開（公告）日: 2022-04-22
• 發(fā)明（設計）人: 楊智;李冬;余海;賀賢漢
• 申請（專利權）人: 杭州博日科技股份有限公司
• 主分類號: G06K9/00
• 分類號: G06K9/00;G06F17/15
• 代理公司: 北京超凡宏宇專利代理事務所(特殊普通合伙) 11463
• 代理人: 董艷芳

1.一種熔解曲線重疊峰的分離方法，其特征在于，包括：

基于預先獲取的PCR反應板各孔熔解曲線的熔點峰位置，利用Savitzky-Golay求導方法確定各孔熔解曲線的二階導數(shù)曲線；

根據(jù)所述二階導數(shù)曲線和所述熔點峰位置，確定原始熔解曲線的初始重疊峰位置；

對所述初始重疊峰位置進行疊加擬合，生成兩個熔點峰曲線；

根據(jù)兩個所述熔點峰曲線確定滿足誤差范圍的目標熔點峰，并對所述目標熔點峰進行分離。

2.根據(jù)權利要求1所述的熔解曲線重疊峰的分離方法，其特征在于，所述PCR反應板各孔熔解曲線的熔點峰位置，對應所述熔點峰的溫度(Tm)值；

利用Savitzky-Golay求導方法確定各孔熔解曲線的二階導數(shù)曲線的步驟，包括：

根據(jù)Savitzky-Golay求導方法，利用預先設定的第一求導參數(shù)計算各孔熔解曲線二階導數(shù)曲線；

其中，所述第一求導參數(shù)包括：具有最高決定系數(shù)的窗口大小和多項式階次。

3.根據(jù)權利要求1所述的熔解曲線重疊峰的分離方法，其特征在于，根據(jù)所述二階導數(shù)曲線和所述熔點峰位置，確定目標熔解曲線的初始重疊峰位置的步驟，包括：

根據(jù)預先給定的第一范圍，確定各熔點峰對應的第一區(qū)域，并進行歸一化處理；所述第一區(qū)域包括在所述第一范圍內選取的所述熔解曲線的二階導數(shù)區(qū)域；

確定所述第一區(qū)域內的所述二階導數(shù)曲線的極小值點；

根據(jù)所述極小值點的相對幅值，確定滿足預先給定的第一閾值的所述極小值點為初始重疊峰位置；預先給定第一閾值的條件包括：在所述第一區(qū)域內，所述二階導數(shù)曲線的極大值點幅值相差不超過兩倍。

4.根據(jù)權利要求3所述的熔解曲線重疊峰的分離方法，其特征在于，預先給定第一范圍的步驟包括：

以所述熔點峰位置對應的所述熔點峰的溫度(Tm)值為基點，以標準熔點峰分辨率為間隔，在所述基點的左右兩端分別確定所述熔點峰對應的波谷；所述標準熔點峰分辨率為最小熔點峰分辨率的四倍；其中，所述最小熔點峰分辨率為2℃；

如果在所述基點的左右兩端均不存在所述波谷，則確定所述標準熔點峰分辨率為第一中間值；所述第一范圍的最小值為所述熔點峰的溫度(Tm)值減去所述第一中間值；所述第一范圍的最大值為所述熔點峰的溫度(Tm)值加上所述第一中間值；

如果在所述基點的左右兩端均存在所述波谷，則確定所述基點與左右兩端所述波谷的距離分別為T_左和T_右；比較T_左和所述標準熔點峰分辨率的大小，確定最小值為第二中間值；比較T_右和所述分辨率的大小，確定最小值為第三中間值；

則所述第一范圍的最小值為所述熔點峰的溫度(Tm)值減去所述第二中間值；所述第一范圍的最大值為所述熔點峰的溫度(Tm)值加上所述第三中間值。

5.根據(jù)權利要求4所述的熔解曲線重疊峰的分離方法，其特征在于，對所述初始重疊峰位置進行疊加擬合，生成兩個熔點峰曲線的步驟，包括：

基于初始重疊峰位置設置初始函數(shù)的初始參數(shù)；

利用Levenberg-Marquardt算法，基于所述初始函數(shù)對各熔點峰對應的所述第一區(qū)域內的所述二階導數(shù)曲線進行疊加擬合，生成兩個待優(yōu)化函數(shù)；

根據(jù)所述初始參數(shù)對所述兩個待優(yōu)化函數(shù)進行優(yōu)化，生成兩個目標函數(shù)；其中，所述目標函數(shù)包括目標參數(shù)；兩個所述目標函數(shù)分別對應兩個所述熔點峰曲線。

6.根據(jù)權利要求5所述的熔解曲線重疊峰的分離方法，其特征在于，

所述目標函數(shù)為三參數(shù)洛倫茲函數(shù)：g(v)＝a·c/[c+(v-b)²]；

其中：v用于表示溫度，a用于表示系數(shù)，b用于表示優(yōu)化之后的重疊峰的位置，c用于表示方差。