基于matlab数字图像处理之低通滤波器6页.doc

实践一:理想低通滤波器、Butterworth低通滤波器、高斯低通滤波器1.1.1理想低通滤波器实践代码： I=imread(couple.bmp);%I=rgb2gray(I); subplot(221),imshow(I);title(原图像);s=fftshift(fft2(I);subplot(223),imshow(abs(s),);title(图像傅里叶变换所得频谱);subplot(224),imshow(log(abs(s),); title(图像傅里叶变换取对数所得频谱); a,b=size(s);a0=round(a/2);b0=round(b/2);d=10;for i=1:afor j=1:b distance=sqrt(i-a0)2+(j-b0)2); if distance=d h=1; else h=0; end; s(i,j)=h*s(i,j);