FFT basic question

Log in or register to post comments
Go To Last Post
7 posts / 0 new
Author
Message
npat_avr
npat_avr's picture
Level: Posting Freak
Joined: Mon. Feb 2, 2009
Posts: 1003 View posts
#1
Posted by npat_avr: Fri. Mar 19, 2010 - 01:20 AM
  • 1
  • 2
  • 3
  • 4
  • 5
Total votes: 0

Hi Guys,

It's been a while since I looked at the FFT in detail so I have a basic question.

What information do you need to plot the FFT of a signal?
Sampling frequency
How many samples (N)

What else?

I am capturing some real time ADC data through an AVR and want to plot the FFT of this signal. I know (somewhat) my sampling frequency and number of samples in one cycle of the wave. I am getting about 28 samples in one cycle of a 1Khz sine wave which means my sampling frequency is 28 ksps.

I am starting with a known sine wave input of known frequency to test my FFT and then proceed to a real life signal.
Thanks.

Tags:

AVR Microcontrollers, megaAVR and tinyAVR
Top
Gold Markclawson
clawson's picture
Level: Moderator
Joined: Mon. Jul 18, 2005
Posts: 101062 View posts
Location: (using avr-gcc in) Finchingfield, Essex, England
Posted by clawson: Fri. Mar 19, 2010 - 10:48 AM
  • 1
  • 2
  • 3
  • 4
  • 5
Total votes: 0

http://elm-chan.org/works/akilcd...

Top
npat_avr
npat_avr's picture
Level: Posting Freak
Joined: Mon. Feb 2, 2009
Posts: 1003 View posts
Posted by npat_avr: Fri. Mar 19, 2010 - 01:36 PM
  • 1
  • 2
  • 3
  • 4
  • 5
Total votes: 0

Looks like all FFT roads lead to elm-chan. :)

I have seen that site before and will look at it once more.

So looks like all I need is sampling frequency and number of FFT points = number of samples in one cycle of a periodic wave?

Top
Bronze Markbobgardner
bobgardner's picture
Level: 10k+ Postman
Joined: Wed. Sep 4, 2002
Posts: 26609 View posts
Location: Orlando Florida
Posted by bobgardner: Fri. Mar 19, 2010 - 02:40 PM
  • 1
  • 2
  • 3
  • 4
  • 5
Total votes: 0

You need to sample for one period of the lowest frequency of interest, at a rate twice the highest frequency of interest. Lets say you are interested in voice band maybe 125 to 4000... so you need to sample for 8ms at 8khz... 64 samples.... 125hz spacing. You get real and imaginary ampl and phase every 125 hz. 64 lines gets you from 0 to 8000, but the top half is the same as the bottom half, so you get 32 spectra 125 hz apart from 0 to 4000. Get the magnitude of the re and im parts at ea line and convert to db.

Imagecraft compiler user

Top
npat_avr
npat_avr's picture
Level: Posting Freak
Joined: Mon. Feb 2, 2009
Posts: 1003 View posts
Posted by npat_avr: Sat. Mar 20, 2010 - 05:25 AM
  • 1
  • 2
  • 3
  • 4
  • 5
Total votes: 0

So for example if I had 256 samples of a signal with sampling frequency of say 10 ksps, my lowest frequency would be about 39 Hz? (10000/256)

Top
ka7ehk
ka7ehk's picture
Level: 10k+ Postman
Joined: Fri. Nov 22, 2002
Posts: 21793 View posts
Location: Tangent, OR, USA
Posted by ka7ehk: Sat. Mar 20, 2010 - 05:59 AM
  • 1
  • 2
  • 3
  • 4
  • 5
Total votes: 0

You have to talk about a "sample window" of N samples. You REALLY want N to be 2**J, where J is a positive integer. That is, 16, 32, 64, etc.

If you just take the samples over the raw window, you will get some strange effects that decrease slowly as the number of samples in the window increases. These effects can be reduced by "weighting" the samples near the edged of the window. There are a number of styles of weighting. You MAY be able to get by without doing this, but be prepared for some things that you might not expect.

Jim

 

Until Black Lives Matter, we do not have "All Lives Matter"!

 

 

Top
sjackson
sjackson's picture
Level: Rookie
Joined: Wed. Apr 11, 2007
Posts: 38 View posts
Posted by sjackson: Tue. Jul 20, 2010 - 08:55 PM
  • 1
  • 2
  • 3
  • 4
  • 5
Total votes: 0

Has anyone actually used the elm-chan library and verified its correctness? I tried making some minor tweaks to use a different windowing function and cut out the square root calculation (for speed purposes) and the whole thing seems to break down. Actually comparing the windowed data after fft_input executes looks wrong when comparing its results to another FFT implementation using the same window.

Best I can tell is that the FMULS16 routine is wrong, but I'd be interested in hearing about the experiences of others.

Top