10
Tôi là người mới bắt đầu trong lập trình và hiện đang cố gắng làm việc trên một dự án yêu cầu triển khai Fast Fourier Transform.Cải thiện tốc độ thực hiện FFT
Tôi đã cho đến nay được quản lý để thực hiện những điều sau:
Có ai có bất kỳ lựa chọn thay thế và gợi ý để cải thiện tốc độ của chương trình mà không mất đi độ chính xác.
short FFTMethod::FFTcalc(short int dir,long m,double *x,double *y)
{
long n,i,i1,j,k,i2,l,l1,l2;
double c1,c2,tx,ty,t1,t2,u1,u2,z;
/* Calculate the number of points */
n = 1;
for (i=0;i<m;i++)
n *= 2;
/* Do the bit reversal */
i2 = n >> 1;
j = 0;
for (i=0;i<n-1;i++) {
if (i < j) {
tx = x[i];
ty = y[i];
x[i] = x[j];
y[i] = y[j];
x[j] = tx;
y[j] = ty;
}
k = i2;
while (k <= j) {
j -= k;
k >>= 1;
}
j += k;
}
/* Compute the FFT */
c1 = -1.0;
c2 = 0.0;
l2 = 1;
for (l=0;l<m;l++) {
l1 = l2;
l2 <<= 1;
u1 = 1.0;
u2 = 0.0;
for (j=0;j<l1;j++) {
for (i=j;i<n;i+=l2) {
i1 = i + l1;
t1 = u1 * x[i1] - u2 * y[i1];
t2 = u1 * y[i1] + u2 * x[i1];
x[i1] = x[i] - t1;
y[i1] = y[i] - t2;
x[i] += t1;
y[i] += t2;
}
z = u1 * c1 - u2 * c2;
u2 = u1 * c2 + u2 * c1;
u1 = z;
}
c2 = sqrt((1.0 - c1)/2.0);
if (dir == 1)
c2 = -c2;
c1 = sqrt((1.0 + c1)/2.0);
}
/* Scaling for forward transform */
if (dir == 1) {
for (i=0;i<n;i++) {
x[i] /= n;
y[i] /= n;
}
}
return(1);
}
Trừ khi bạn cần tự viết nó cho mục đích hiểu, FFTW (http://www.fftw.org/) là một thư viện tuyệt vời. Đó là một việc thực hiện tự điều chỉnh, siêu nhanh và đáng tin cậy và bạn có thể gọi nó từ C++ tốt (xem http://www.fftw.org/faq/section2.html#cplusplus) –
Tôi thích FFTReal rất nhiều. http://ldesoras.free.fr/prod.html –
Tại sao bạn viết thực hiện của riêng bạn thay vì sử dụng một trong vô số thư viện trên mạng, có khả năng tất cả nhanh hơn, được kiểm tra tốt hơn, chính xác hơn và có nhiều tính năng hơn? – PlasmaHH