%FFT
clear;
clear all;
r_total=[25:25:200];
%定义半径范围
path_fig=['C:\Documents and Settings\oee\桌面\陈立贞编程作业\'];
%定义图像存储路径
for i=1:1
figure;
I=zeros(512,512);
%定义一副512x512的图像数组
[m,n]=size(I);
r=r_total(i);
t=(0:0.001:2)*pi;
x=m/2+r.*cos(t);
y=n/2+r.*sin(t);
%画圆形范围
iptsetpref('ImshowAxesVisible','on') % 显示坐标值
x1=round(x);
y1=round(y);
for j=1:2001
I(x1(j),y1(j))=1;
j;
end
%将圆形边界值设为1
for m=1:512
index=find(I(m,:)>0);
I(m,min(index):max(index))=1;
end
%将圆面值设为1
imshow(I,'notruesize');
%显示定义的图像
hold on
%在当前坐标轴图形里添加绘制图形
name=[path_fig num2str(r) '.png'];
print(gcf,'-dpng',name);
%保存定义的图像
F=fft2(I);
%对自定义图像进行快速傅里叶变换
imshow(fftshift(log(abs(F))),[-15]);
%显示频谱图像,频谱的零频率系数被移到频谱中间
colormap(jet);colorbar
%为频谱图像添加colorbar
hold on
name=[path_fig num2str(r) 'fft2.png'];
print(gcf,'-dpng',name);
%保存fft频谱图像
close;
end