Matlab图像的巴特沃斯低通滤波

%用MATLAB实现Butterworth低通滤波器的代码所示：
I=imread（‘9.jpg‘）;
figure（1）;
subplot（22[1 3]）
imshow（I）;                          %显示原始图像
title（‘原始图像‘）
I1=imnoise（I‘salt‘）;%添加椒盐噪声
subplot（222）;
imshow（I1）;
title（‘椒盐噪声图像‘）;
f=double（I1）;
g=fft2（f）;               %对图像傅里叶变换
g=fftshift（g）;           %图象频谱中心从矩阵的原点移到矩阵的中心
subplot（224）;
imshow（log（abs（g））[]）;
title（‘快速傅立叶变换（FFT）图像‘）;
for k=1:4
figure（2）;
[N1N2]=size（g）;
n=2;  
d0=1+（k-1）*25;              %截止频率=1265176
n1=fix（N1/2）;
n2=fix（N2/2）;
for i=1:N1
  for j=1:N2
        d=sqrt（（i-n1）^2+（j-n2）^2）;
          h=1