Power from a Fourier transform

[neil.thompson]

So I have been away from education for a little while now and I'm going through some refresher stuff - in particular I have been playing around with FFTs.

If i take (with MATLAB notation):

time = 0:0.01:10
y = fft(sin(2*pi*f*time))

with f = 5
then the maximum amplitude of the fft output is about 498.

with f = 10
the maximum amplitude of fft output is 492.

I understand the amplitude is 'halved' in both cases because this fft is ambiguous so the energy is spread over two peaks. But why is the energy less when the frequency increases? I have more cycles in the case with more frequency, but I suppose this means I have less samples. Also, is it usual to normalise this in some way? It seems like this is something you wouldn't want if you were dealing were plotting energy return from doppler shifts.

 

[DrClaude]

But why is the energy less when the frequency increases?

How is the energy related to the FT? Have you checked the width of the peaks?

Also, is it usual to normalise this in some way? It seems like this is something you wouldn't want if you were dealing were plotting energy return from doppler shifts.

There are different ways of normalizing the FFT, but it is customary that doing FFT^-1(FFT(f)) will return N × f, although MATLAB returns the original result.

 

[neil.thompson]

How is the energy related to the FT? Have you checked the width of the peaks?
.

Good point - my language was imprecise. What I meant was peak power, or peak amplitude -- I guess what I'm asking is: if I have a frequency resolution on the FFT such that all of the signal should be confined to a signal frequency bin (e.g., each bin spans is 2kHz, my input signal is a single tone at 4.5kHz so all of the signal should fold into the third bin) so why does the height of the amplitude peak depend on the frequency of the signal?

edit: and into the bargain, related to your point, if I have enough data to get a very fine resolution FT then I guess this will reveal that there will be some spread across bins so the peak isn't the same because the curve is wider/narrower.

There are different ways of normalizing the FFT, but it is customary that doing FFT^-1(FFT(f)) will return N × f, although MATLAB returns the original result.

Right, thanks.

 

[DrClaude]

Have you tried to take a finer grid an see what you get then?