# Extracting amplitudes of components from FFT

I have a 4096-point FFT of a continuous pipe organ sound (no windowing). I wish to extract the relative amplitudes of the various harmonically-related sinusoidal components for voicing analysis. I can easily find the bins with the local maximum power for each component, but since the frequencies are not exactly in the middle of a bin, some of these peaks are spread between two adjacent bins, and the amplitude in each bin is less than it would be if the frequency were shifted to exactly hit that bin. For frequency analysis I have used various interpolation formulas with great success. But now I am doing amplitude analysis, and I want to know if there is a similar interpolation formula to estimate the amplitude of a component in a way that is independent of whether that component is exactly in the center of a frequency bin or is between two bins. My current best guess is the take the sum of the powers (magnitude squared) of the bin with the peak power together with the bins immediately adjacent to it. The sum of those three bin powers would then represent the power of that one component of the sound. Is there a better way?

• Perhaps a window function would be good, but it won't solve the problem I have stated. It is still possible to get an actual frequency that is in-between two bins, and in that case the power in those bins will be less than if the frequency was exactly on a bin - windowing or not. – Robert Scott Feb 22 '15 at 17:37