Forum: PC-Programmierung fftw in c++ mit 2-dim messdaten


von Martin F. (mael15)


Lesenswert?

hallo,

ich habe ein messsignal einer sonde, welches ich auf dem bildschirm in x 
und y koordinaten als punkt sinnvoll darstellen kann.
davon möchte ich eine fft mit hilfe von fftw machen, alle beispiele 
beziehen sich aber auf bilder oder eindimensionale werte, helfen also 
nicht so passend weiter.
muss ich jetzt die amplitude mit fftw_plan_dft_r2c_1d oder doch x und y 
werte mit fftw_plan_dft_r2c_2d nehmen?
und falls 2d, wäre folgendes dann korrekt, besonders das füllen des 
arrays mit x und y werten?
1
double* samples = (double*) malloc(2 * smplCount * sizeof(double));
2
fftw_complex* spectrum = (fftw_complex*) fftw_malloc((smplCount / 2 + 1) * sizeof(fftw_complex));
3
memset(samples, 0, 2 * smplCount * sizeof(double));
4
memset(spectrum, 0, (smplCount / 2 + 1) * sizeof(fftw_complex));
5
6
fftw_plan forward = fftw_plan_dft_r2c_2d(smplCount, 2, samples, spectrum, 0);
7
for (int i = 0; i < smplCount; i++) {
8
  samples[i] = input->get(i).x;
9
  samples[smplCount + i] = input->get(i).y;
10
}
11
12
fftw_execute(forward);
13
// auswertung, in meinem fall darstellung als bild
14
free(samples);

danke!

Bitte melde dich an um einen Beitrag zu schreiben. Anmeldung ist kostenlos und dauert nur eine Minute.
Bestehender Account
Schon ein Account bei Google/GoogleMail? Keine Anmeldung erforderlich!
Mit Google-Account einloggen
Noch kein Account? Hier anmelden.