【導(dǎo)讀】高級(jí)技術(shù)計(jì)算語(yǔ)言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。測(cè)試和測(cè)量、財(cái)務(wù)建模和分析以及計(jì)算生物學(xué)等眾多應(yīng)用領(lǐng)域。附加的工具箱（單。內(nèi)特定類(lèi)型的問(wèn)題。它以矩陣運(yùn)算為基礎(chǔ)，把計(jì)算、可視化、程序設(shè)計(jì)融合在一。它可以將聲音文件變換為離散的數(shù)據(jù)文件，然后利用其強(qiáng)大的矩陣運(yùn)。數(shù)字信號(hào)處理是MATLAB重要應(yīng)用的領(lǐng)域之一。NN次復(fù)數(shù)加法運(yùn)算。隨著N的增加，運(yùn)算量將急劇增加，量的時(shí)間和機(jī)器內(nèi)存，不能滿(mǎn)足實(shí)時(shí)的要求。因此，DFT的這種運(yùn)算只能進(jìn)行。理論上的計(jì)算，不適合對(duì)實(shí)時(shí)處理要求高的場(chǎng)合。

　　

【正文】 ft2(array)。 array=abs(array)。 array=array(1:r,1:c)。 subplot(2,3,3)。 imshow(array,[])。 title(39。自建 IFFT2結(jié)果 39。)。 function pushbutton4_Callback(hObject, eventdata, handles) [filename, pathname]=uigetfile({39。*.jpg。*.tif。*.bmp。*.gif39。 },39。File Selector39。)。 image=imread(strcat(pathname,filename))。%讀取圖片 if ndims(image)==3 %ndims(a)表示矩陣 a的維數(shù)，即 length(size(a)) image=rgb2gray(image)。 end scrsz=get(0,39。ScreenSize39。)。 set(gcf,39。Name39。,39?？焖俑道锶~變換 39。)。 subplot(2,3,4)。 imshow(image)。 title(39。原始圖像 39。)。 function pushbutton5_Callback(hObject, eventdata, handles) [filename, pathname]=uigetfile({39。*.jpg。*.tif。*.bmp。*.gif39。 },39。File Selector39。)。 image=imread(strcat(pathname,filename))。%讀取圖片 if ndims(image)==3 %ndims(a)表示矩陣 a的維數(shù)，即 length(size(a)) image=rgb2gray(image)。 end F=fft2(image)。 FC=fftshift(F)。 S=log(1+abs(FC))。 課程設(shè)計(jì)說(shuō)明書(shū) 24 subplot(2,3,5) imshow(S,[])。 title(39。內(nèi)置 FFT2結(jié)果 39。)。 function pushbutton6_Callback(hObject, eventdata, handles) [filename, pathname]=uigetfile({39。*.jpg。*.tif。*.bmp。*.gif39。 },39。File Selector39。)。 image=imread(strcat(pathname,filename))。%讀取圖片 if ndims(image)==3 %ndims(a)表示矩陣 a的維數(shù)，即 length(size(a)) image=rgb2gray(image)。 end F=fft2(image)。 FC=fftshift(F)。 S=log(1+abs(FC))。 array=ifft2(F)。 array=round(abs(array))。 subplot(2,3,6)。 imshow(array,[])。 title(39。內(nèi)置 IFFT2結(jié)果 39。)。 return function figure1_CreateF(hObject, eventdata, handles) function pushbutton7_Callback(hObject, eventdata, handles) [filename, pathname]=uigetfile({39。*.jpg。*.tif。*.bmp。*.gif39。 },39。File Selector39。)。 image=imread(strcat(pathname,filename))。%讀取圖片 if ndims(image)==3 %ndims(a)表示矩陣 a的維數(shù)，即 length(size(a)) image=rgb2gray(image)。 end scrsz=get(0,39。ScreenSize39。)。 set(gcf,39。Name39。,39?？焖俑道锶~變換 39。)。 subplot(2,3,1)。 imshow(image)。 title(39。原始圖像 39。)。 subplot(2,3,4)。 imshow(image)。 title(39。原始圖像 39。)。 [r,c]=size(image)。 課程設(shè)計(jì)說(shuō)明書(shū) 25 array=image。 t=log2(r)。 t1=floor(t)。 % floor(x)：下取整，即舍去正小數(shù)至最近整數(shù) t2=ceil(t)。 % ceil(x)：上取整，即加入正小數(shù)至最近整數(shù) if t1~=t2 array(2^t2,c)=0。 end [r1,c1]=size(array)。 t=log2(c1)。 t3=floor(t)。 t4=ceil(t)。 if t3~=t4 array(r1,2^t4)=0。 end [r1,c1]=size(array)。 n=r1/2。 data_col=zeros(1,n,39。double39。)。 %按列方向計(jì)算時(shí)用到的 for m=1:n data_col(m)=exp(1i*2*pi*(m1)/r1)。 end n=c1/2。 data_row=zeros(1,n,39。double39。)。 %按行方向計(jì)算時(shí)用到 for m=1:n data_row(m)=exp(1i*2*pi*(m1)/r1)。 end array=transform_fft2(array)。 Ft=fftshift(array)。 S1=log(1+abs(Ft))。 subplot(2,3,2)。 imshow(S1,[])。 title(39。自建 FFT2函數(shù)結(jié)果 39。)。 array=transform_ifft2(array)。 array=abs(array)。 array=array(1:r,1:c)。 課程設(shè)計(jì)說(shuō)明書(shū) 26 subplot(2,3,3)。 imshow(array,[])。 title(39。自建 IFFT2結(jié)果 39。)。 F=fft2(image)。 FC=fftshift(F)。 S=log(1+abs(FC))。 subplot(2,3,5) imshow(S,[])。 title(39。內(nèi)置 FFT2結(jié)果 39。)。 array=ifft2(F)。 array=round(abs(array))。 subplot(2,3,6)。 imshow(array,[])。 title(39。內(nèi)置 IFFT2結(jié)果 39。)。 return function array=transform_fft2(array) array=double(array)。 [r1 c1]=size(array)。 for j=1:r1 array(j,:)=transform_fft(array(j,:))。 end for j=1:c1 array(:,j)=transform_fft((array(:,j)))。 end function array1=transform_fft(array) N=length(array)。 n=N/2。 w=zeros(1,n,39。double39。)。 for m=1:n w(m)=exp(1i*2*pi*(m1)/N)。 end p=log2(N)。 array1=zeros(1,N,39。double39。)。 for q=1:p 課程設(shè)計(jì)說(shuō)明書(shū) 27 t1=2^(q1)。 t2=2^(p1)。 for k=0:(t2/t11) for j=0:(t11) if mod(q,2)==1 data1=array(k*t1+j+1)。 data2=array(k*t1+j+t2+1)。 array1(k*t1*2+j+1)=data1+data2。 array1(k*t1*2+j+t1+1)=(data1data2)*w(k*t1+1)。 else data1=array1(k*t1+j+1)。 data2=array1(k*t1+j+t2+1)。 array(k*t1*2+j+1)=data1+data2。 array(k*t1*2+j+t1+1)=(data1data2)*w(k*t1+1)。 end end end end if mod(p,2)==1 return else array1=array。 return end function array=transform_ifft2(array) array=conj(array)。 [r1,c1]=size(array)。 for i=1:r1 array(i,:)=transform_fft(array(i,:))。 end for i=1:c1 array(:,i)=transform_fft(array(:,i))。 end array=array/(r1*c1)。 課程設(shè)計(jì)說(shuō)明書(shū) 28 function pushbutton8_Callback(hObject, eventdata, handles) clc。 close all。 close(gcf)。