A method for improving accuracy of estimating natural frequency of a distribution system. The estimation is performed in order to identify a rectifier load that generates harmonics due to commutation oscillation. It includes the steps of measuring a charging current flowing into a power factor improving capacitor to obtain first data, removing fundamental frequency from the first data to obtain second data, determining, from the second data, a first time at which the amplitude of commutating oscillation components takes a maximum value, picking up third data in the interval beginning from the first time and ending at a second time later than the first time by T/N, where T is the period of the fundamental frequency, and N is the number of rectifier phases, performing frequency analysis on the third data to form fourth data by expanding the third data into a Fourier series, and determining a local peak frequency, that is, the natural frequency in the fourth data. This prevents frequency components of particular orders from being emphasized, and reduces interference between the commutation oscillations, thereby making it possible to obtain a local peak value, that is, the natural frequency with ease at high accuracy.