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

t4=ceil(t); if t3~=t4

array(r1,2^t4)=0; end

[r1,c1]=size(array); 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)); subplot(2,4,7); imshow(S1,[]);

title('自建FFT2函数结果'); array=transform_ifft2(array); array=abs(array); array=array(1:r,1:c); subplot(2,4,8); imshow(array,[]);

title('自建IFFT2函数结果'); F=fft2(image); FC=fftshift(F); S=log(1+abs(FC)); subplot(2,4,3) imshow(S,[]);

title('内置FFT2函数结果'); array=ifft2(F);