【正文】 的 x(n)的順序數(shù)是用 M 位二進制數(shù)表示，則倒序數(shù)只需將順序數(shù)的二進制位倒置即可，按照這一規(guī)律用硬件電路和匯編語言很容易產(chǎn)生倒序數(shù)。十進制倒序數(shù)用 J 表示，與之對應(yīng)的二進制數(shù)用 JB 表示。依此類推，直到完成最高位加 1，逢 2 向右進位的運算。 課程設(shè)計說明書 8 5 Matlab 程序?qū)崿F(xiàn) Matlab 提供的 FFT 函數(shù)是一個計算 DFT 的智能程序，能自動選擇快速算法進行 DFT 運算。 圖 51 代碼輸入窗口 3. 運行后就會出現(xiàn)如圖 52 所示的窗口。 GUI 是包含圖形對象，如窗口、圖標、菜單和文本等圖文并茂的用戶界面。 以本設(shè)計要求 為例介紹 ：首先運行 Matlab 軟件，輸入 GUIDE，出現(xiàn)一個窗口，如圖 61 所示。 課程設(shè)計說明書 11 圖 62 設(shè)計界面 第二步， 設(shè)置控件屬性 ： 雙擊組件可以設(shè)置文本框，推按鈕的屬性，如顯示大小，名稱 和字體 等?？丶O(shè)置完成后保存，然后運行 GUI(操作為 ctrl+T)，就會進入 editor 窗口， 加入各個控件功能的函數(shù)代碼。 圖 63 系統(tǒng)人機對話界面 課程設(shè)計說明書 12 運行調(diào)試 運行 GUI，已經(jīng)彈出圖 63 所示的系統(tǒng)人機對話界面，然后點擊各推按鈕就會得到本設(shè)計的要求。從通過設(shè)計任務(wù)的分析到功能的劃分再到確定總體的設(shè)計方案，再到最后的每個功能部分的具體實現(xiàn)，整個過程都需要我們組的人充分利用所學(xué)的知識進行思考借鑒。通過本課程設(shè)計使學(xué)生了解數(shù)字圖像的基本概念，掌握數(shù)字圖像處理的基本內(nèi)容， FFT 基本原理和 Matlab 實現(xiàn)方法。最后經(jīng)過努力，又深入學(xué)習(xí)了圖形用戶界面（ GUI），完成了選做要求的人機對話界面。通過這次課設(shè)我學(xué)會了，綜合運用課程中所學(xué)到的理論知識去獨立完成一個設(shè)計課題，通過查閱手冊和文獻資料，培養(yǎng)了我獨立分析 和解決實際問題的能力。 課程設(shè)計說明書 14 參考文獻 [1] 范壽康 .《 DSP 技術(shù)與 DSP 芯片》 .北京電子工業(yè)出版社， 2020 [2] 程佩青 .《數(shù)字信號處理教程》 .北京清華大學(xué)出版社出版， 2020 [3] 高西全，丁玉美等 .《數(shù)字信號處理》 .北京電子工業(yè)出版社， 2020 [4] 余成波，陶紅艷 .《數(shù)字 信號處理及 MATLAB 實現(xiàn)》 .北京清華大學(xué)出版社， 2020 [5] 曹弋，趙陽 .《 MATLAB 實用教程》 .北京電子工業(yè)出版社， 2020 [6] 奧本海姆 .《 離散時間信號處理 》 .科學(xué)出版社， 2020 [7] 宗孔德 ， 胡廣書 .《 數(shù)字信號處理 》 .清華大學(xué)出版社， 1997 課程設(shè)計說明書 15 附錄 1 Matlab 源程序： function image_process_FFT() x=imread(39。ScreenSize39。,[0 0 scrsz(3)1 scrsz(4)])?？焖俑道锶~變換 39。 title(39。 end subplot(2,3,4)。)。 t1=floor(t)。 t=log2(c1)。 end [r1,c1]=size(array)。)。double39。 Ft=fftshift(array)。 title(39。 array=abs(array)。 title(39。 FC=fftshift(F)。內(nèi)置 FFT2 函數(shù)結(jié)果 39。 subplot(2,3,6)。)。 end for j=1:c1 array(:,j)=transform_fft((array(:,j)))。double39。 array1=zeros(1,N,39。 t2=2^(p1)。 array1(k*t1*2+j+t1+1)=(data1data2)*w(k*t1+1)。 array(k*t1*2+j+t1+1)=(data1data2)*w(k*t1+1)。 for i=1:r1 array(i,:)=transform_fft(array(i,:))。 gui_State = struct(39。, gui_Singleton, ... 39。, untitled_OutputF, ... 39。, [])。 end if nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。 function varargout = untitled_OutputF(hObject, eventdata, handles) varargout{1} = 。)。)。)。 S=log(1+abs(FS))。)。)。 S=log(1+abs(FS))。 title(39。 x=imread(39。 imshow(f)。 function pushbutton5_Callback(hObject, eventdata, handles) axes()。 if ndims(image)==3 image=rgb2gray(image)。 t1=floor(t)。 t=log2(c1)。 end [r1,c1]=size(array)。)。double39。 Ft=fftshift(array)。自建 FFT 圖 39。 for j=1:r1 array(j,:)=transform_fft(array(j,:))。 w=zeros(1,n,39。 end p=log2(N)。 for q=1:p t1=2^(q1)。 課程設(shè)計說明書 21 array1(k*t1*2+j+1)=data1+data2。 array(k*t1*2+j+1)=data1+data2。 image=imread(39。 end [r,c]=size(image)。 t2=ceil(t)。 t3=floor(t)。 n=r1/2。 for m=1:n data_col(m)=exp(1i*2*pi*(m1)/r1)。)。 S1=log(1+abs(Ft))。 imshow(array,[])。 function array=transform_fft2(array) array=double(array)。 n=N/2。 for m=1:n w(m)=exp(1i*2*pi*(m1)/N)。)。 data2=array(k*t1+j+t2+1)。 data2=array1(k*t1+j+t2+1)。 return end function array=transform_ifft2(array) array=conj(array)。 end array=array/(r1*c1)。 。 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。 end for j=1:c1 array(:,j)=transform_fft((array(:,j)))。自建 IFFT 圖 39。 array=abs(array)。 end array=transform_fft2(array)。 data_row=zeros(1,n,39。double39。 if t3~=t4 array(r1,2^t4)=0。 end [r1,c1]=size(array)。 t=log2(r)。)。 end end end end if mod(p,2)==1