基于虚拟心脏的心电逆问题求解论文.docVIP

基于虚拟心脏的心电逆问题求解 Mission 摘要心电逆问题研究的目标是要从无损测得的体表电位来找出心脏的对应状态。本文提出了一种新的心电逆问题求解方法, 即从体表电位求出虚拟心脏的模型参数, 由所求得的虚拟心脏的模型参数即可确定心脏所处的状态。该法把欠定的逆问题求解问题转化为正问题的参数优化问题, 避免了以往心电逆问题求解方法中的病态问题。 关键字：　虚拟现实; 心电逆问题;W PW 预激综合征 SOLUTION FOR ELECTROCARD IOGRAPHIC INVERSE PROBLEM BASED ON VIRTUAL HEART Yangwenhua 200710407112 ABSTRACTW e in t roduced a new app roach to electrocardiographic inverse problem studies, i.e. , solution s in term s of virtual heart model parameters. The heart state can be determined from examination of the recovered virtual heart model parameters. Thu s the ill-conditioned inverse problem was turned to be a forward problem with parameter optimization, and the ill-posed nature w as avoided. Key words: 　V irtual Reality; Electrocardiographic Inverse Problem; W PW preexcitation syndrome 　 基于虚拟心脏的心电逆问题求解方法 具备了虚拟心脏后, 通过设置其模型参数我们就可以仿真不同心脏状态下的体表电位。例如, 通过去除虚拟心脏模型中某一区域的心肌单元可以仿真心肌梗塞; 通过减慢束支的兴奋传导速度可以仿真束支传导阻滞; 通过设置预激点可以仿真预激综合征; 通过增加心室的心肌单元可以仿真心室肥大等等。反过来, 对于给定的体表电位分布数据, 如果能找出对应的虚拟心脏模型参数, 那么我们就可以定量地确定心脏的对应状态。这就是我们基于虚拟心脏的心电逆问题研究新方法的基本原理。其原理框图如图1 所示。从病人身上测量得到的体表电位数据经一个基于知识的预诊断程度大致定性地确定出属于哪一类心脏病; 然后模型编辑器自动设置虚拟心脏的模型参数并启动仿真运算, 仿真得出的体表电位与测量的体表电位经匹配准则进行比较, 如果符合匹配要求, 则由所设置的虚拟心脏模型参数即可确定对应的心脏状态; 否则, 执行调整模型参数→仿真→比较判断的参数优化循环过程直至符合匹配准则为止。 图1　基于虚拟心脏心电逆问题求解方法的原理框图 　 一．　匹配准则 仿真得出的体表电位与实测的体表电位是不可能完全一致的, 那么仿真得出的体表电位与实测的体表电位到底符合到何种程度, 我们就要用一定的准则来判断它。 本研究采用了三条匹配准则, 它们基于下述三种误差: (1) 总平方误差, squerr, 由下式给出 　　Squerr= [(t,i)- (t,i)]2 (1) 式中（t,i）代表体表第i 点在时刻t 的电位测量值, (t,i)为相应的电位仿真值, n 为测量点个数, T 为心动周期。squerr 的量纲为mV 2。 (2) 总最大值位置误差,maxperr, 定义为 maxp err = C‖Pmmax ( t) - Psmax ( t) ‖ (2) 　　 式中Pmmax ( t) 代表 t 时刻电位测量值的最大值位置, Psmax ( t) 为t 时刻电位仿真值的最大值位置, C 为常数, ‖?‖为“网格距离”, 定义为 ‖a - b‖ =　　 (3) 式中 m = ( in t)- ( in t) n = |a - b|% 20, 　if n 10, then　n = 20 – n (4) 　　 上式中, a、b 代表图2 中200 个电极位置分布的序列编号(0～ 199)。|….|是绝对值操作, % 是余数操作, 也就是说, n 是指a 与b 之差的绝对值除以20 所得的余数, 如果该余数大于10, 则n由20 减去该余数所得。 例如, 图2 中A 与B 之间的“网格距离”为 ‖A - B ‖ = ‖67 - 109‖ = (3 - 5) 2 + 22 = 2 2 图2 中C 与D 之间的“网格距离”为 ‖C - D ‖ =