The invention provides an electromagnetic transient simulation algorithm containing multiple switching elements, a system, a medium and equipment, and relates to the technical field of electromagnetictransient simulation of a power system. The method comprises the following steps: if switching action is carried out in (t-delta t, t], the switching action time is td1, and performing interpolationcalculation on the node voltage at the moment t and the voltage and current of each element by utilizing the node voltage at the moment td1 + delta t/2 and the voltage and current of each element at the moment td1 + delta t; when entering the next simulation step length (t, t + delta t), if a switching action exists in the (t-delta t, t), calculating the node voltage and the voltage and current ofeach element at the t + delta t moment by adopting two backward Euler methods with the step length of delta t/2 at the t moment; otherwise, calculating by adopting a trapezoidal integral method withthe step length of delta t; and repeating the steps in sequence to obtain the voltage of the whole time node and the voltage and current of each element. The electromagnetic transient simulation algorithm can process multiple switching actions with the same step length, and is high in simulation precision and high in universality.