work21.m

通过迭代算法来实现快速付立叶变换

clear all
clc 
close all
A =  imread('Lena.bmp');     %读入图像
figure
imshow(A);                   %显示原图像
A = double(A);
B = Myfft2(A);               %快速付立叶变换

B = abs(B);
B = B/max(B(:));
B = imadjust(B);                      
hh  = log10(B(:));           %求系数的幅度分布(对数分布)
figure
hist(hh)                     %显示系数分布
figure 
imshow(uint8((B)));      %显示付立叶变换后的图