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

關(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//ffttw_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ì)算大碼為：voidtw_gen_ifft（float*w，intn）{inti，j，k;constdoublePI=3.141592654;for（j=1，k=0;j《=n》》2;j=j《《2）{for（i=0;i《n》》2;i+=j）{#ifdef_LITTLE_ENDIANw［k］=（float）-sin（2*PI*i/n）;//為負(fù)數(shù)的時(shí)候是進(jìn)行IFFT，為正數(shù)的時(shí)候是進(jìn)行FFTw［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）;#elsew［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）;#endifk+=6;}}}注釋：此函數(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;closeall;clearall;方法一：%生成一個(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），：）;%每一行分別做一維FFTfori=1:heightb（i，：）=fft（a（i，：））;end%行列翻轉(zhuǎn)c=transpose（b）;%每一行分別做一維FFT（由于翻轉(zhuǎn)了，此時(shí)相當(dāng)于每一列分別做一維FFT）fori=1:widthd（i，：）=fft（c（i，：））;end%行列翻轉(zhuǎn)e=transpose（d）;%每一行分別做一維IFFTfori=1:heightib（i，：）=ifft（e（i，：））;end%行列翻轉(zhuǎn)ic=transpose（ib）;%每一行分別做一維IFFT（由于翻轉(zhuǎn)了，此時(shí)相當(dāng)于每一列分別做一維FFT）fori=1:widthid（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)如下：//ffttw_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ì)每一行做一次FFTfor（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)行FFTfor（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，1