The invention relates to a distributed harmonic source identification method based on singe-frequency current transmission characteristics. The method comprises the following steps that (1) a current harmonic source identification model is established; (2) constraint conditions are exported; (3) a current harmonic observation matrix is established and preprocessed; (4) the harmonic source maximum identification frequency T and a non-Gaussianity threshold value F in an electricity network are set, and k is made to be equal to 0; (5) the k is made to be equal to k+1, and a FastICA method is used for solving a mixing matrix Wh; (6) the inverse matrix Wh-1 of the Wh is calculated, the phase position and the row order of the Wh-1 are adjusted according to the constraint conditions, the Wh-1 is denoted as Ph, if the Ph meets the constraint conditions, the step (7) is executed, and otherwise the step (5) is executed again; (7) the reverse matrix of the Ph is calculated, a demixing matrix Wh=Ph-1 is updated, the non-Gaussianity sum phi of identification results is calculated, if phi is no smaller than the threshold value F, the demixing matrix Wh is output, the harmonic source identification result is calculated, the positions of main harmonic sources are analyzed, iteration is ended, and otherwise the step (8) is executed; (8) when k is smaller than or equal to T, the step (5) is executed again, when k is larger than T, the demixing matrix Wh and the non-Gaussianity threshold value are output, and the positions of the main harmonic sources in the network are further analyzed.