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!