The invention provides a power system frequency calculation method of windowing spectral line interpolation. The method includes the following steps that 1, voltage sampling data of ten continuous sampling cycles of a public connection point are obtained through a voltage measurement circuit; 2, time domain windowing function processing is conducted on the voltage sampling data; 3, local discrete Fourier transformation is conducted on the data obtained after windowing, and voltage amplitude values of specific harmonic wave times are obtained; 4, for the calculated specific voltage amplitude values, interpolation is conducted between the two adjacent maximal amplitude values; 5, calculation is conducted through a polynomial curve fitting method, a frequency coefficient is obtained, and thus fundamental wave frequency of a power grid is determined. By means of the method, a blackman window function is applied to N sampling points, local discrete Fourier transformation (DFT) is only conducted on designated harmonic frequency , it is proved through experiments that by means of the method, the spectrum leakage inhibition effect is good, the calculation speed is high, the calculation process is simple, precision is high, universality is high, and the method can be used for high-precision electric energy metering and electric energy quality analysis.