基于gui的fft運(yùn)算-資料下載頁(yè)
2025-06-27 17:31本頁(yè)面
【正文】 ndims(image)==3 %ndims(a)表示矩陣a的維數(shù)，即length(size(a)) image=rgb2gray(image)。endF=fft2(image)。FC=fftshift(F)。S=log(1+abs(FC))。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)。endF=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。)。returnfunction 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)。endscrsz=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)。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)。endn=c1/2。data_row=zeros(1,n,39。double39。)。 %按行方向計(jì)算時(shí)用到for m=1:n data_row(m)=exp(1i*2*pi*(m1)/r1)。endarray=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)。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。)。returnfunction array=transform_fft2(array)array=double(array)。[r1 c1]=size(array)。for j=1:r1 array(j,:)=transform_fft(array(j,:))。endfor j=1:c1 array(:,j)=transform_fft((array(:,j)))。endfunction 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)。endp=log2(N)。array1=zeros(1,N,39。double39。)。for q=1:p 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 endendif mod(p,2)==1 return else array1=array。 returnendfunction array=transform_ifft2(array)array=conj(array)。[r1,c1]=size(array)。for i=1:r1 array(i,:)=transform_fft(array(i,:))。endfor i=1:c1 array(:,i)=transform_fft(array(:,i))。endarray=array/(r1*c1)。function pushbutton8_Callback(hObject, eventdata, handles)clc。close all。close(gcf)。28
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1