以下主要是通過代碼調(diào)用fft函數(shù)

1、通過dsplib庫調(diào)用DSPF_sp_ifftSPxSP函數(shù)來實現(xiàn)fft

//fft

tw_gen_fft （w_array， N）;

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

其中函數(shù)tw_gen_fft（）為旋轉(zhuǎn)因子，計算大碼為：

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）;//為負(fù)數(shù)的時候是進(jìn)行IFFT，為正數(shù)的時候是進(jìn)行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;

}

}

}

注釋：此函數(shù)中的N為此行的點數(shù)，而不是行列的總數(shù)，在在fft2d中出現(xiàn)了錯誤，就是由這個N導(dǎo)致的。

2、通過matlab代碼來驗證fft和fft2d的關(guān)系。

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

%行列翻轉(zhuǎn)

c = transpose（b）;

%每一行分別做一維FFT（由于翻轉(zhuǎn)了，此時相當(dāng)于每一列分別做一維FFT）

for i = 1:width

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

end

%行列翻轉(zhuǎn)

e = transpose（d）;

%每一行分別做一維IFFT

for i = 1:height

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

end

%行列翻轉(zhuǎn)

ic = transpose（ib）;

%每一行分別做一維IFFT（由于翻轉(zhuǎn)了，此時相當(dāng)于每一列分別做一維FFT）

for i = 1:width

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

end

%行列翻轉(zhuǎn)

ie = transpose（id）;

方法二：

afft = fft2（a）;

aifft = ifft2（afft）;

方法三：

f = fft（fft（a）?！??！?

ift = ifft（ifft（f）?！?。’;

3、通過上述matlab的驗證可以確定fft和fft2d的關(guān)系，針對16*16的數(shù)據(jù)在DSP中進(jìn)行實驗如下：

//fft

tw_gen_fft （w_array， col）;

//如果圖像的行列不同的話在第一次fft之前n應(yīng)為對應(yīng)行的點數(shù)，則應(yīng)在行列翻轉(zhuǎn)后再次調(diào)用tw_gen_fft函數(shù)，并將n改為對應(yīng)的列的點數(shù) //對每一行做一次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）; } //行列翻轉(zhuǎn) DSPF_dp_mat_trans（（double *）y_array， height， width， （double *）x_array）;

tw_gen_fft （w_array， row）; //對每一列（翻轉(zhuǎn)后變?yōu)樾校┻M(jìn)行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）; } //行列翻轉(zhuǎn) DSPF_dp_mat_trans（（double *）y_array， height， width， （double *）x_array）;

至此，完美解決dsp平臺下的fft2d的內(nèi)聯(lián)函數(shù)調(diào)用問題