This paper proposes a higher order cumulant-based method for harmonic retrieval in non-Gaussian ABMA noise.Hilbert transform is used to transform the real measruements into their complex counter-part

then a kind of elaborately defined fourth-order cumulant of this complex measurements is employed to iden-tify the All parameters of the non-Gaussian noise. After prefiltering the noisy measurements with the identifiedAR polynomial

SVD-TLS method can be applied to retrieve harmonics.This method releases the asymmetrically-distributed assumption of the non-Gaussian noise and is also effcient when there is quadratic phase coupling inharmonics.Simulation examples are presented to support the result of this paper.