This article presents the implementation of digital filters used in digital relay protection current measuring elements. Mathematical descriptions of the considered digital filters as well as the computer programs for their coefficients calculation are described. It has been shown that in order to reliable estimate the digital filter performance its input signals waveforms must be close to the actual secondary current waveform of the current transformer to which the digital protection with the estimated digital filter is connected. For these purposes in MatLab–Simulink dynamic simulation environment the power system and the current measuring element models were developed. Performed calculations allowed to reveal that the exponentially decaying DC component which in some cases contains in primary fault current drives the current transformer core into saturation even when its nominal parameters are not exceeded. This results in distortion of the current transformer secondary current which in this case contains higher and inter-harmonics. Moreover, such harmonic content is not completely taking into account during coefficients calculation of the considered digital filters what results in signal magnitude estimation inaccuracy. Comparison of the digital filters response to the above-mentioned input signals allowed to find out such digital filter implementations which enable signal magnitude estimation with a minimum error. Ways of filtering quality improvement concerned with the window functions are proposed. Thus, the joint usage of digital filter and Hamming window allows to achieve the zero value of the signal magnitude gain factor in high-frequency range and substantially suppress all spectral components above 100 Hz. The increasing of the signal magnitude settling time in this case can be reduced by choosing the most optimal parameters of the all components of the current measuring element.