[发明专利]心脏电兴奋传播的仿真方法

摘要

本发明公开了一种心脏电兴奋传播的仿真方法。按定义仿真步骤如下1.T=To；2.若H(i，j，k)∈w，则对Aijk(l，m，n)执行规则I，II，III，IV；3.记4.令T＝TE；5。如果T＜TL转2；6.结束。本发明的优点是仿真时间短，误差小，时间分辨率高。

主权项

1.一种心脏电兴奋传播的仿真方法：定义1在26邻点立方体网格中，心肌单元记为H(i，j，k)，其邻近单元记为Aijk(l，m，n)＝H(i±l，j±m，k±n)(l，m，n可为1，0，-1)H(i，j，k)的起始兴奋时刻记为T(i，j，k)，其邻近单元的起始兴奋时刻记为Tijk(l，m，n)，当前时刻记为T；定义2H(i，j，k)的兴奋传播速度记为v(i，j，k)，H(i，j，k)与Aijk(l，m，n)之间的距离记为dijk(l，m，n)，则兴奋从H(i，j，k)传播到Aijk(l，m，n)所需的时间为tijk(l，m，n)＝dijk(l，m，n)/v(i，j，k)定义3仿真步长Δt定义为Δt＝min(tijk(l，m，n)i，j，k，l，m，n)定义4H(i，j，k)的有效不应期(或绝对不应期)、相对不应期和超常期分别记为HERP(i，j，k)、HRRP(i，j，k)和HSNP(i，j，k)，Aijk(l，m，n)的有效不应期记为AERP(l，m，n)，H(i，j，k)的自律周期记为CL(i，j，k)；定义5H(i，j，k)在时刻T的动作电位幅度定义为ΦH(i，j，k)，其正常兴奋阈值定义为QH，在相对不应期和超常期内的兴奋阈值分别记为QHRRP和QHSNP，Aijk(l，m，n)在时刻T的动作电位幅度定义为ΦA(l，m，n)，其正常兴奋阈值定义为QA；定义6H(i，j，k)的有效兴奋期(即在此期间有能力使相邻心肌产生兴奋)定义为HEAP(i，j，k)＝T－T(i，j，k)，andΦH(i，j，k)＞QA定义7仿真中遇到的判断条件：C1：T－Δt＜T(i，j，k)＜TC2：T－T(i，j，k)≤HERP(i，j，k)C3：HERP(i，j，k)＜T－T(i，j，k)≤HERP(i，j，k)＋HRRP(i，j，k)C4：HERP(i，j，k)＋HRRP(i，j，k)＜T－T(i，j，k)≤HERP(i，j，k)＋HRRP(i，j，k)＋HSNP(i，j，k)C5：HERP(i，j，k)＋HRRP(i，j，k)＋HSNP(i，j，k)＜T－T(i，j，k)≤CL(i，j，k)C6：T－T(i，j，k)≥CL(i，j，k)C7：ΦA(l，m，n)＜QHC8：ΦA(l，m，n)≥QHSNPC9：ΦA(l，m，n)≥QHRRPC10：H(i，j，k)为传入阻滞单元；定义8在任意时刻，H(i，j，k)可为以下四种状态之一：(1)不应态R：C2∨(C3∧C9)∨(C4∧C8)∨(C5∧C7)(2)波面态W：C1∨C6(3)可兴奋态E：C10∧((C3∧C9)∨(C4∧C8)∨(C5∧C7))(4)待兴奋态P：为一中间跟踪状态仿真规则1：如果满足C6，则令T(i，j，k)＝T(i，j，k)＋CL(i，j，k)仿真规则II：如果H(i，j，k)∈W且Aijk(l，m，n)∈E，则令Tijk(l，m，n)＝T(i，j，k)＋tijk(l，m，n)仿真规则III：如果H(i，j，k)∈W且Aijk(l，m，n)∈R，若T(i，j，k)＋HEAP(i，j，k)＞Tijk(l，m，n)＋AERP(l，m，n)则令Aijk(l，m，n)∈P仿真规则IV：如果H(i，j，k)∈P且H(i，j，k)∈E，则令T(i，j，k)＝Tijk(l，m，n)＋tijk(l，m，n)特征在于其仿真步骤如下：Step1T＝To(To为起始仿真时刻)；Step2若H(i，j，k)∈W，则对Aijk(l，m，n)执行规则I，II，III，IV；Step3记TE=MIN{T(i,j,k)+CL(i,j,k)ifT(i,j,k)<TT(i,j,k),ifT(i,j,k)>TStep4令T＝TE；Step5如果T＜TL转Step2(TL为给定的仿真时间长度)Step6结束。