我发现了一种使用FFTW3的FFT库函数计算离散正弦变换的方法。但我应该说它仍然比MATLAB dst（）函数慢（我仍然不知道为什么）。让我简单地解释一下我所做的：

//x->Your first array with some values
//y->FFT input array
//z->FFT output array
//n->First array(x) size

#define REAL 0
#define IMAG 1

function_Name(){
const int n = 1000;
double x[n][2];

//Fill x with some value
for(i){
x[i][REAL]=i+1;
x[i][IMAG]=0;
}

fftw_complex y[(n+1)*2];
fftw_complex z[(n+1)*2];

y[0][REAL]=0;
y[0][IMAG]=0;

for(){
//Add x to y from 2. element(y[1])
}

y[n+1][REAL]=0;
y[n+1][IMAG]=0;

for(){
//Invert x and multiply by "-" and add to y
}

fftw_plan plan = fftw_plan_dft_1d((n+1)*2, y, z, FFTW_FORWARD, FFTW_ESTIMATE);

fftw_execute(plan);

**//You can use z but first you have to Invert (fliplr in MATLAB) it and take only first n elements and only IMAG parts.**

fftw_destroy_plan(plan);
fftw_cleanup();

}

我知道这个算法并不完美，但我正在努力加快速度，并使一个干净。对于速度差异MATLAB的fft（）函数和FFTW的FFT函数仍然是一个问题。我将打开一个新的问题。