反卷积 C++

Question

我从http://rosettacode.org/wiki/Deconvolution/1D#C获取了这个 一维去卷积的示例代码，但它似乎不能正常工作。 在尝试构建此项目时，我首先遇到这个错误，'_fft'：非法使用类型“void”。如果有帮助的话，这段代码最初是为C编写的，我不得不进行一些更改，例如将include从complex.h更改为complex，并使用命名空间std。谢谢你的帮助。

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <complex>
using namespace std;

double PI;
complex<double> cplx;




void _fft(cplx buf[], cplx out[], int n, int step)
 {
if (step < n) {
    _fft(out, buf, n, step * 2);
    _fft(out + step, buf + step, n, step * 2);

    for (int i = 0; i < n; i += 2 * step) {
        cplx t = cexp(-I * PI * i / n) * out[i + step];
        buf[i / 2]     = out[i] + t;
        buf[(i + n)/2] = out[i] - t;
    }
}
}

void fft(cplx buf[], int n)
{
cplx out[n];
for (int i = 0; i < n; i++) out[i] = buf[i];
_fft(buf, out, n, 1);
}

cplx *pad_two(double g[], int len, int *ns)
{
int n = 1;
if (*ns) n = *ns;
else while (n < len) n *= 2;

cplx *buf = calloc(sizeof(cplx), n);
for (int i = 0; i < len; i++) buf[i] = g[i];
*ns = n;
return buf;
}

void deconv(double g[], int lg, double f[], int lf, double out[]) {
int ns = 0;
cplx *g2 = pad_two(g, lg, &ns);
cplx *f2 = pad_two(f, lf, &ns);

fft(g2, ns);
fft(f2, ns);

cplx h[ns];
for (int i = 0; i < ns; i++) h[i] = g2[i] / f2[i];
fft(h, ns);

for (int i = 0; i >= lf - lg; i--)
    out[-i] = h[(i + ns) % ns]/32;
free(g2);
free(f2);
}

int main()
{
PI = atan2(1,1) * 4;
double g[] = {24,75,71,-34,3,22,-45,23,245,25,52,25,-67,-96,96,31,55,36,29,-43,-7};
double f[] = { -3,-6,-1,8,-6,3,-1,-9,-9,3,-2,5,2,-2,-7,-1 };
double h[] = { -8,-9,-3,-1,-6,7 };

int lg = sizeof(g)/sizeof(double);
int lf = sizeof(f)/sizeof(double);
int lh = sizeof(h)/sizeof(double);

double h2[lh];
double f2[lf];

printf("f[] data is : ");
for (int i = 0; i < lf; i++) printf(" %g", f[i]);
printf("\n");

printf("deconv(g, h): ");
deconv(g, lg, h, lh, f2);
for (int i = 0; i < lf; i++) printf(" %g", f2[i]);
printf("\n");

printf("h[] data is : ");
for (int i = 0; i < lh; i++) printf(" %g", h[i]);
printf("\n");

printf("deconv(g, f): ");
deconv(g, lg, f, lf, h2);
for (int i = 0; i < lh; i++) printf(" %g", h2[i]);
printf("\n");
}

- Logan Adams

Answer 1

如果没有运行这个代码，cplx被当作类使用，但它实际上是一个变量。你能更改这一行吗：

complex<double> cplx;

like this

typedef complex<double> cplx;

- abellina

---

我也发现了同样的问题。看起来这个有相当多的额外损坏。我不确定这是否曾在任何环境中编译过。 - marko

就像我说的那样，我在网上找到了它，并且无法验证它是否在任何环境中编译过，但是我按照建议更改了typedef，但是cexp函数不起作用。 有什么想法吗？ 感谢您的帮助。 - Logan Adams

我发现exp被重载为复数。所以你应该替换为cexp。第40行也有一个无效的转换，你需要将其转换为cplx*。我还发现了一个未定义的I，我想这实际上是1。最后一个问题是第59行的operator/。我建议你做更多的研究，如果你有更多的问题，请在这里重新发布，我们都愿意帮助你解决问题。 - abellina