数字信号处理基于MATLAB的FFT算法设计课设说明书(5)

课程设计说明书

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)=(data1-data2)*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);

20

课程设计说明书

附录Ⅱ

以下是用GUI界面实现的程序编码：

function pushbutton1_Callback(hObject, eventdata, handles) axes(handles.axes1); x=imread('sxy.jpg'); imshow(x); title('原图');

function pushbutton2_Callback(hObject, eventdata, handles) axes(handles.axes2); x=imread('sxy.jpg'); f=rgb2gray(x); imshow(f); title('灰度图');

function pushbutton3_Callback(hObject, eventdata, handles) axes(handles.axes3); img=imread('sxy.jpg'); f=rgb2gray(img); F=fft2(f); FS=fftshift(F); S=log(1+abs(FS)); imshow(S,[]); title('FFT图');

function pushbutton4_Callback(hObject, eventdata, handles) axes(handles.axes4); img=imread('sxy.jpg'); f=rgb2gray(img); F=fft2(f); FS=fftshift(F); S=log(1+abs(FS)); fr=real(ifft2(ifftshift(FS))); ret=im2uint8(mat2gray(fr)); imshow(ret,[]); title('IFFT图');

function pushbutton5_Callback(hObject, eventdata, handles)

21

课程设计说明书

axes(handles.axes5); x=imread('sxy.jpg'); imshow(x); title('原图');

function pushbutton6_Callback(hObject, eventdata, handles) axes(handles.axes6); image=imread('sxy.jpg'); if ndims(image)==3 image=rgb2gray(image); end

imshow(image); title('灰度图');

function pushbutton7_Callback(hObject, eventdata, handles) axes(handles.axes7); image=imread('sxy.jpg'); if ndims(image)==3 image=rgb2gray(image); end

[r,c]=size(image); array=image; t=log2(r); t1=floor(t); t2=ceil(t); 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);

22

课程设计说明书

n=r1/2;

data_col=zeros(1,n,'double'); for m=1:n

data_col(m)=exp(-1i*2*pi*(m-1)/r1); end n=c1/2;

data_row=zeros(1,n,'double'); for m=1:n

data_row(m)=exp(-1i*2*pi*(m-1)/r1); end

array=transform_fft2(array); Ft=fftshift(array); S1=log(1+abs(Ft)); imshow(S1,[]); title('FFT图');

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,'double'); for m=1:n

w(m)=exp(-1i*2*pi*(m-1)/N); end p=log2(N);

array1=zeros(1,N,'double');

23

课程设计说明书

for q=1:p t1=2^(q-1); t2=2^(p-1); for k=0:(t2/t1-1) for j=0:(t1-1) 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)=(data1-data2)*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)=(data1-data2)*w(k*t1+1); end end end end

if mod(p,2)==1 return else

array1=array; return end

function pushbutton8_Callback(hObject, eventdata, handles) axes(handles.axes8); image=imread('sxy.jpg'); if ndims(image)==3 image=rgb2gray(image); end

[r,c]=size(image); array=image; t=log2(r);

24

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库数字信号处理基于MATLAB的FFT算法设计课设说明书(5)在线全文阅读。