【正文】 close all。for i=1:r1 array(i,:)=transform_fft(array(i,:))。 array(k*t1*2+j+t1+1)=(data1data2)*w(k*t1+1)。 array1(k*t1*2+j+t1+1)=(data1data2)*w(k*t1+1)。 t2=2^(p1)。array1=zeros(1,N,39。double39。endfor j=1:c1 array(:,j)=transform_fft((array(:,j)))。)。subplot(2,3,6)。內(nèi)置FFT2結(jié)果39。FC=fftshift(F)。title(39。array=abs(array)。title(39。Ft=fftshift(array)。double39。)。end[r1,c1]=size(array)。t=log2(c1)。t1=floor(t)。)。subplot(2,3,4)。imshow(image)。,39。ScreenSize39。)。*.bmp。)。subplot(2,3,6)。FC=fftshift(F)。)。*.bmp。)。S=log(1+abs(FC))。image=imread(strcat(pathname,filename))。*.gif39。function pushbutton5_Callback(hObject, eventdata, handles)[filename, pathname]=uigetfile({39。imshow(image)。,39。ScreenSize39。)。*.bmp。)。subplot(2,3,3)。S1=log(1+abs(Ft))。)。 %按列方向計(jì)算時(shí)用到的for m=1:n data_col(m)=exp(1i*2*pi*(m1)/r1)。n=r1/2。t3=floor(t)。t2=ceil(t)。[r,c]=size(image)。Name39。endscrsz=get(0,39。File Selector39。*.tif。自建FFT2函數(shù)結(jié)果39。S1=log(1+abs(Ft))。)。 %按列方向計(jì)算時(shí)用到的for m=1:n data_col(m)=exp(1i*2*pi*(m1)/r1)。n=r1/2。t3=floor(t)。t2=ceil(t)。[r,c]=size(image)。Name39。endscrsz=get(0,39。File Selector39。*.tif。原始圖像39。)。set(gcf,39。%讀取圖片if ndims(image)==3 %ndims(a)表示矩陣a的維數(shù)，即length(size(a)) image=rgb2gray(image)。 },39。*.jpg。endfunction untitled_OpeningF(hObject, eventdata, handles, varargin) = hObject。amp。, [] , ... 39。, untitled_OpeningF, ... 39。, mfilename, ... 39。endarray=array/(r1*c1)。 returnendfunction array=transform_ifft2(array)array=conj(array)。 data2=array1(k*t1+j+t2+1)。 data2=array(k*t1+j+t2+1)。)。 for m=1:n w(m)=exp(1i*2*pi*(m1)/N)。n=N/2。[r1 c1]=size(array)。title(39。array=ifft2(F)。subplot(2,3,5)imshow(S,[])。)。subplot(2,3,3)。)。subplot(2,3,2)。 %按行方向計(jì)算時(shí)用到for m=1:n data_row(m)=exp(1i*2*pi*(m1)/r1)。endn=c1/2。data_col=zeros(1,n,39。t4=ceil(t)。 % ceil(x)：上取整，即加入正小數(shù)至最近整數(shù) if t1~=t2 array(2^t2,c)=0。array=image。title(39。原始圖像39。)。set(gcf,39。)。%圖片路徑、文件名稱image=imread(strcat(pathname,filename))。*.gif39。參考文獻(xiàn)[1]范壽康，DSP技術(shù)與DSP芯片，北京，電子工業(yè)出版社，2008[2]程佩青，數(shù)字信號處理教程，北京，清華大學(xué)出版社出版，2001[3]高西全，丁玉美等，數(shù)字信號處理，北京，電子工業(yè)出版社，2009[4]奧本海姆，2000[5]宗孔德，1997 [6]余成波，陶紅艷，數(shù)字信號處理及MATLAB實(shí)現(xiàn)，北京，清華大學(xué)出版社，2008 [7]曹弋，趙陽，MATLAB實(shí)用教程，北京，電子工業(yè)出版社，2007 附錄I MATLAB設(shè)計(jì)程序 function image_process_FFT()[filename, pathname]=uigetfile({39。程序運(yùn)行調(diào)試初期，曾經(jīng)多次出現(xiàn)錯(cuò)誤、不能產(chǎn)生圖形等問題，但在我翻閱資料認(rèn)真改正及老師同學(xué)的幫助下基本功能還是完成了，經(jīng)過1個(gè)星期的上機(jī)實(shí)習(xí)，程序已得到一些完善，能完成基本的要求的功能。同時(shí)我相信，進(jìn)一步加強(qiáng)對MATLAB的學(xué)習(xí)與研究對我今后的學(xué)習(xí)將會(huì)起到很大的幫助。程序運(yùn)行調(diào)試時(shí)，自己選擇輸入圖片，可以實(shí)現(xiàn)FFT及IFFT的圖片。 如圖所示為點(diǎn)擊“自編FFT”推按鈕得到的只顯示自編IFFT結(jié)果的圖片。系統(tǒng)人機(jī)對話界面如圖52所示：圖52人機(jī)對話界面運(yùn)行GUI，已經(jīng)彈出圖52的系統(tǒng)人機(jī)對話界面，然后點(diǎn)擊“原始圖片”推按鈕，根據(jù)提示選擇圖片，就得到自編程序的原始圖像。控件設(shè)置完成后保存，然后運(yùn)行GUI（操作為ctrl+T），就會(huì)進(jìn)入editor窗口，加入各個(gè)控件功能的函數(shù)代碼。完成人機(jī)設(shè)計(jì)界面如圖51所示：圖51人機(jī)設(shè)計(jì)界面第二步，設(shè)置控件屬性。5系統(tǒng)人機(jī)對話界面圖形用戶界面（GUI），是一種提供人機(jī)交互的工具和方法。同時(shí)利用MATLAB軟件系統(tǒng)內(nèi)置的FFT及IFFT程序進(jìn)行相應(yīng)的運(yùn)算。整個(gè)蝶形運(yùn)算流程圖如圖34所示：圖34整個(gè)蝶形運(yùn)算流程圖4MATLAB程序?qū)崿F(xiàn)使用MATLAB軟件運(yùn)行程序使得圖片實(shí)現(xiàn)采集灰度，并對灰度進(jìn)行FFT及IFFT運(yùn)算。在進(jìn)行第L級運(yùn)算時(shí)，先算出該級不同旋轉(zhuǎn)因子的個(gè)數(shù)(也是該級中各個(gè)蝶形運(yùn)算兩輸入數(shù)據(jù)的間距)，再從R=0開始按序計(jì)算，直到R=B1結(jié)束。應(yīng)用原位計(jì)算，蝶形運(yùn)算可表示成如下形式： 總結(jié)上述運(yùn)算規(guī)律，可采用如下運(yùn)算方法進(jìn)行DITFFT運(yùn)算。蝶形運(yùn)算是分級進(jìn)行的；每級的蝶形運(yùn)算可以按旋轉(zhuǎn)因子的指數(shù)大小排序進(jìn)行；如果指數(shù)大小一樣則可從上往下依次蝶算。為了編寫DITFFT算法的運(yùn)算程序，首先要分析其運(yùn)算規(guī)律，總結(jié)編程思想并繪出程序框圖。FFT對DFT的數(shù)學(xué)運(yùn)算量改進(jìn)：直接采用DFT進(jìn)行計(jì)算，運(yùn)算量為次復(fù)數(shù)乘法和次復(fù)數(shù)乘法。若N不能滿足上述關(guān)系，可以將序列x(n)補(bǔ)零實(shí)現(xiàn)。為了能實(shí)時(shí)處理DFT,要想減少D