理想低通滤波器（频率域滤波）

%频率域滤波
clc;
close all;
img=imread(‘Fig3.10(b).jpg’);
img=im2double(img);
Fimg=fft2(double(img));%傅里叶变换
Fimg=fftshift(Fimg);%将变换的原点移到频率矩形的中心
[M,N]=size(img);
%半径为5的理想低通滤波处理
dist1=5;
z1=zeros(M,N);
for i=1:M
for j=i:N
if(sqrt(((i-M/2)^2+(j-N/2)2))<dist1)
z1(i,j)=1;
end
end
end
g1=Fimg.*z1;
g1=ifftshift(g1);
img1=real(ifft2(g1));
%半径为15的理想低通滤波处理
dist2=15;
z2=zeros(M,N);
for i=1:M
for j=i:N
if(sqrt(((i-M/2)^2+(j-N/2)2))<dist2)
z2(i,j)=1;
end
end
end
g2=Fimg.*z2;
g2=ifftshift(g2);
img2=real(ifft2(g2));
%半径为30的理想低通滤波处理
dist3=20;
z3=zeros(M,N);
for i=1:M
for j=i:N
if(sqrt(((i-M/2)^2+(j-N/2)2))<dist3)

z3(