In this paper the use of frequency-matched linear numerical integration techniques for the simulation of systems modelled by the dynamic phasors approach (DPA) is proposed and investigated. Such methods show an increased accuracy and computational efficiency around the matched frequency. Such frequency matched methods are derived by analyzing the local truncation error in the frequency domain and determining the coefficients of the method in such a way so that the error is minimized around the matched frequencies. The aim in this way is to increase the accuracy and the efficiency of the numerical calculation, during the inherent fast oscillations featured in such systems modelled with dynamic phasors approach.