A fast and accurate method for Generalized Harmonic Analysis is presented. Thefrequency of sinusoid is estimated around a peak of the Fourier spectrum using binarysearch. The amplitude and the phase are estimated to minimize the squared sum of theresidual signal. This process is repeated till the prescribed number of sinusoids areextracted. Sinusoidal parameters are recalculated to reduce errors introduced by thepeak detection using the windowed Discrete‐Time Fourier Transform. Audio signals areanalyzed by the proposed method, which confirms the accuracy compared to existingmethods. The proposed algorithm has high degree of concurrency and is suitablyimplemented on Graphical Processing Unit (GPU).