This project presents the theoretical background in electrophysiology that is used as a basis in the development of numerical methods for the simulation of heart’s electrical activity. We discuss the mathematical models currently used in cardiac electrophysiology research to simulate the propagation of an action potential wave, and then focus on the description of the ionic currents. This study allows us to anticipate the numerical model sensitivity with regard to the discretization, with the aim of recovering the physiological behaviors of the tissue (i.e. conduction velocity, electrical restitution, etc). Thereafter, we introduce the pseudo-ECG signal reconstruction method from the action potential map. This leads us to examine the approximation of such a signal using the finite element basis coming from the discretization of the monodomain equations. We discuss then a particular cardiac rhythm disorder: action potential alternans. In a second step, we perform an extended simulation study in order to validate the previous numerical methods. We carefully check the convergence of the conduction velocity with respect to the mesh size, and we bring out the advantage of phenomenological approach in terms of electrical restitutions and computational cost. We reproduce typical rhythm disorders and show that the pseudo-ECG signal is able to detect them. Finally, we discuss the influence of anisotropy in cardiac wave propagation and alternans development, characterizing the effect of fibers direction and pacing site in the resulting spatio-temporal distribution.