Inverse DFT in C++, was mache ich falsch?

[Piti V. (piti)]

Hallo zusammen,
ich möchte die Daten eines Sensors auswerten und habe dazu Daten 
eingelesen, DFT transformiert (mit evtl. Befensterung) und möchte die 
Daten nun per inverser Fourier rücktransformieren. Code für diese 
Funktion befindet sich unten.
Funktioniert auch alles ganz gut, allerdings stimmt das Ergebnis der 
iDFT nicht. Bekomme dabei das Bild im Anhang (1. Input, 2. Platzhalter, 
3. DFT, 4. iDFT. dh Bild 4 stimmt nicht). Hier müsste ja ohne 
Befensterung die selbe Funktion erhalten werden wie auch der Input war, 
stattdessen habe ich einen merkwürdigen Phasenverschiebung und eine 
niedrigere Amplitude. Kann jemand einen Fehler erkennen?

Habe in einem Matlab-Skript das Ganze nachgebaut, da funktioniert alles 
wie es soll.

Vielen Dank im Voraus.

// Durchführung der inversen DFT
void DFT_C::do_inv_dft(void)
{
int i,k;
double temp_r_inv,temp_i_inv; 

if (!data_in) return;

data_out_r_inv=new double[d_length];
data_out_i_inv=new double[d_length];
data_out_b_inv=new double[d_length];

if (data_out_r_inv&&data_out_i_inv&&data_out_b_inv) data_out_flag=1;
else return;

for (i=0;i<d_length;i++)
  {
  temp_r_inv=0;
  temp_i_inv=0;
  for (k=0;k<d_length;k++)
    {
    temp_r_inv+=data_out_r[k]*cos(2*M_PI*k*i/d_length); //data_out kommt aus der DFT
    temp_i_inv+=data_out_i[k]*sin(2*M_PI*k*i/d_length);
    }
  data_out_r_inv[i]=temp_r_inv;
  data_out_i_inv[i]=temp_i_inv;
  data_out_b_inv[i]=sqrt(pow(temp_r_inv,2)+pow(temp_i_inv,2));
  }
}

[Quarkltasche (Gast)]

Hi,
Dein Codeschnipsel ist ohne weiteren Zusammenhang etwas
schwer zu durchschauen.
Ist "d_length" die Gesamtlaenge des DFT Vektors, der
von 0 ... N/2 -1 die positiven und von N/2 ... N-1
die negativen Frequenzen enthaelt (bei einer N-Punkt DFT)?
Oder hast Du ueberhaupt nur (wie Dein DFT Bild zeigt) den
positiven Anteil betrachtet?

Insgesamt sieht es so aus, als haettest Du die negativen
Frequenzen vergessen (dafuer spricht auch die halbierte
Amplitude) oder irgendwie falsch behandelt (Deine Schleife
mit dem Zaehler k erzeugt fuer k >=N/2 falsche Frequenzen).

Gruss & viel Erfolg bei der Fehlersuche
 Quarktasche

[Dieter (Gast)]

Nimm doch einfach die von fftw.org. Gibt kaum eine bessere 
Implementierung.

Es sei denn, es ist eine Übung an der Uni, dann einfach mal einen Sinus 
reingeben und im Debugger nachschauen, wo der fehler liegt.

[Michael W. (Gast)]

>Nimm doch einfach die von fftw.org.
Das ist doch eine FFT, oder?
Der TE will ja DFT.

[2^5 (Gast)]

Eine FFT ist eine mögliche Implementation einer DFT.

[Schnellhilfe (Gast)]

Tipp: Lass dir aus deinem Matlab-Code den C++-Code von Matlab 
autogenerieren und schau ihn dir dann an.