關(guān)于DSP中fft函數(shù)調(diào)用方法

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

1、通過(guò)dsplib庫(kù)調(diào)用DSPF_sp_ifftSPxSP函數(shù)來(lái)實(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)因子，計(jì)算大碼為：

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

}

}

}

注釋?zhuān)捍撕瘮?shù)中的N為此行的點(diǎn)數(shù)，而不是行列的總數(shù)，在在fft2d中出現(xiàn)了錯(cuò)誤，就是由這個(gè)N導(dǎo)致的。

2、通過(guò)matlab代碼來(lái)驗(yàn)證fft和fft2d的關(guān)系。

clc;

close all;

clear all;

方法一：

%生成一個(gè)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)了，此時(shí)相當(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)了，此時(shí)相當(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、通過(guò)上述matlab的驗(yàn)證可以確定fft和fft2d的關(guān)系，針對(duì)16*16的數(shù)據(jù)在DSP中進(jìn)行實(shí)驗(yàn)如下：

//fft

tw_gen_fft （w_array， col）;

//如果圖像的行列不同的話在第一次fft之前n應(yīng)為對(duì)應(yīng)行的點(diǎn)數(shù)，則應(yīng)在行列翻轉(zhuǎn)后再次調(diào)用tw_gen_fft函數(shù)，并將n改為對(duì)應(yīng)的列的點(diǎn)數(shù) //對(duì)每一行做一次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）; //對(duì)每一列（翻轉(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平臺(tái)下的fft2d的內(nèi)聯(lián)函數(shù)調(diào)用問(wèn)題