[Improvement on simulation algorithm of excitation propagation in heart modeling]

It is important to simulate the excitation propagation process of cardiac bio-electricity in the research of ECG forward problem. Traditional methods describe them with wave simulation algorithm such as LFX simulation algorithm and vector propagation algorithm etc, these methods have some problems to certain extent, due to the presence of discreteness of space and time and asymmetry of the myocardium. This paper discussed the simulation algorithm in 2-dimension space under the circumstance of layered and non-layered structure of myocardium. By calculating the theoretic values of simulating time based on Huygen's principle, we found that there were errors in LFX algorithm and no errors in vector propagation algorithm under the circumstance of non-layered structure of the myocardium, no mater what myocardium is isotropic or anisotropic. However, there exist errors from both algorithms when the myocardium has the layered structure. An improved algorithm is proposed and the simulations have been performed to examine the efficacy of the new algorithm, and the errors are reduced obviously. By increasing the number of myocardial blocks in the model, we also analyzed its influence on the error of simulation algorithm.