以下主要是通過代碼調用fft函數

1、通過dsplib庫調用DSPF_sp_ifftSPxSP函數來實現fft

//fft

tw_gen_fft （w_array， N）;

DSPF_sp_fftSPxSP （N， x_array， w_array， y_array， brev， rad， 0， N）;

其中函數tw_gen_fft（）為旋轉因子，計算大碼為：

void tw_gen_ifft （float *w， int n）

{

int i， j， k;

const double PI = 3.141592654;

for （j = 1， k = 0; j 《= n 》》 2; j = j 《《 2）

{

for （i = 0; i 《 n 》》 2; i += j）

{

#ifdef _LITTLE_ENDIAN

w［k］ = （float） -sin （2 * PI * i / n）;//為負數的時候是進行IFFT，為正數的時候是進行FFT

w［k + 1］ = （float） cos （2 * PI * i / n）;

w［k + 2］ = （float） -sin （4 * PI * i / n）;

w［k + 3］ = （float） cos （4 * PI * i / n）;

w［k + 4］ = （float） -sin （6 * PI * i / n）;

w［k + 5］ = （float） cos （6 * PI * i / n）;

#else

w［k］ = （float） cos （2 * PI * i / n）;

w［k + 1］ = （float） -sin （2 * PI * i / n）;

w［k + 2］ = （float） cos （4 * PI * i / n）;

w［k + 3］ = （float） -sin （4 * PI * i / n）;

w［k + 4］ = （float） cos （6 * PI * i / n）;

w［k + 5］ = （float） -sin （6 * PI * i / n）;

#endif

k += 6;

}

}

}

注釋：此函數中的N為此行的點數，而不是行列的總數，在在fft2d中出現了錯誤，就是由這個N導致的。

2、通過matlab代碼來驗證fft和fft2d的關系。

clc;

close all;

clear all;

方法一：

%生成一個16x16的單位矩陣

%a = eye（16）;

%［m， n］ = size（a）;

a1 = ［0:15］;

height = 16;

width = 16;

%［m， n］ = size（a）;

a = a1（ones（1，16）， ：）;

%每一行分別做一維FFT

for i = 1:height

b（i，：） = fft（a（i，：））;

end

%行列翻轉

c = transpose（b）;

%每一行分別做一維FFT（由于翻轉了，此時相當于每一列分別做一維FFT）

for i = 1:width

d（i，：） = fft（c（i，：））;

end

%行列翻轉

e = transpose（d）;

%每一行分別做一維IFFT

for i = 1:height

ib（i，：） = ifft（e（i，：））;

end

%行列翻轉

ic = transpose（ib）;

%每一行分別做一維IFFT（由于翻轉了，此時相當于每一列分別做一維FFT）

for i = 1:width

id（i，：） = ifft（ic（i，：））;

end

%行列翻轉

ie = transpose（id）;

方法二：

afft = fft2（a）;

aifft = ifft2（afft）;

方法三：

f = fft（fft（a）。‘）。’;

ift = ifft（ifft（f）。‘）。’;

3、通過上述matlab的驗證可以確定fft和fft2d的關系，針對16*16的數據在DSP中進行實驗如下：

//fft

tw_gen_fft （w_array， col）;

//如果圖像的行列不同的話在第一次fft之前n應為對應行的點數，則應在行列翻轉后再次調用tw_gen_fft函數，并將n改為對應的列的點數 //對每一行做一次FFT for（i = 0; i 《 height; i++） { DSPF_sp_fftSPxSP（16， （float *）x_array + 2*i*width ， w_array， （float *）y_array + 2*i*width， brev， 4， 0， 16）; } //行列翻轉 DSPF_dp_mat_trans（（double *）y_array， height， width， （double *）x_array）;

tw_gen_fft （w_array， row）; //對每一列（翻轉后變為行）進行FFT for（j = 0; j 《 width; j++） { DSPF_sp_fftSPxSP（16， （float *）x_array + 2*j*width， w_array， （float *）y_array + 2*j*width， brev， 4， 0， 16）; } //行列翻轉 DSPF_dp_mat_trans（（double *）y_array， height， width， （double *）x_array）;

至此，完美解決dsp平臺下的fft2d的內聯函數調用問題