ChaN FFT auf dem XMega

[Basti (Gast)]

Hallo,

bin gerade dabei eine FFT für Audiosignale zum laufen zu bekommen.
Davon bin ich aber noch weit entfernt... erstmal hab ich mit 
unterschiedlichen Eingangsparametern gespielt...

Eine C Implementierung hat mir vom Quellcode und der Ausgabe zugesagt. 
Leider war diese deutlich langsamer als der von ChaN...

Okay, nun bin ich bei seiner gelandet und kann mir das Ausgabespektrum 
in folgendem Source nicht ganz erklären:

void test_fft2(void)
{
  
  int16_t capture[FFT_N];    /* Wave captureing buffer */
  complex_t bfly_buff[FFT_N];    /* FFT buffer */
  uint16_t spektrum[FFT_N/2];    /* Spectrum output buffer */
  
  memset(capture,0x00,sizeof(capture));
  memset(bfly_buff,0x00,sizeof(bfly_buff));
  memset(spektrum,0x00,sizeof(spektrum));
  
  for (uint16_t i=0; i<FFT_N; i++)
  {
//    capture[i] = 1000*sin(i*((2*3.1415926535)/32));
    if(i < FFT_N/2) {
      capture[i] = 1000;
      } else {
      capture[i] = -1000;
    }
  }
  
  gpio_toggle_pin(STATE_LED_IR);
  
  fft_input(capture, bfly_buff);
  fft_execute(bfly_buff);
  fft_output(bfly_buff, spektrum);
  
  gpio_toggle_pin(STATE_LED_IR);
  
  for (uint16_t i = 0; i < FFT_N/2; i++)
  {
    udi_cdc_write_buf(&capture[i],2);
    udi_cdc_putc(0x00);
    udi_cdc_write_buf(&spektrum[i],2);
  }

  delay_ms(1000);
  
}

Eingangswert - Nullbyte - Spektrumausgabe

E8 03 00 0A 00 
E8 03 00 AF 02 
E8 03 00 86 01 
E8 03 00 E4 00 
E8 03 00 9C 00 
E8 03 00 89 00 
E8 03 00 64 00 
E8 03 00 62 00 
E8 03 00 4A 00 
E8 03 00 4D 00 
E8 03 00 3A 00 
E8 03 00 3E 00 
E8 03 00 32 00 
E8 03 00 36 00 
E8 03 00 2A 00

c0 ist hier fast 0... aber auch nicht wirklich...
bei dem Rechtecksignal, sollte es doch nur gradzahlige Harmonische 
geben. Wieso sind hier immer welche da?

Bei der C FFT hat das selbe Rechtecksignal wunderbar funktioniert...

Beim Sinus passt soweit alles, außer das der Peak etwas breiter und 
gleichmäßiger in die Ränder geht, was der C Code vorher auch nicht getan 
hat.

Muss ich noch irgendwas einstellen, oder was mach ich hier falsch?

Ein paar digitale Rundungsfehler okay, aber das Rechteck ist nicht so 
wie vorgesehen...

Danke,

Basti

[Basti (Gast)]

Hm, weiß keiner was darüber? Ich habe irgendwie die Vermutung, dass die 
"Seitenbänder" einer Art Rundung/Mittelwertbildung zum Opfer fallen...

Man sieht, dass die gerade Harmonischen kleiner sind als die 
benachbarten ungeraden...

Evtl. ist das sogar für meine Anwendung nicht weiter schlimm... 
verstehen würde ich es trotzdem gern... ;)