butterworth低通滤波 matlab 程序

clear all;
I=imread('D:\1.jpg');
I=rgb2gray(I);
subplot(221),imshow(I);xlabel('(a)原始图像');

fftI1=fft2(I);                       %二维离散傅立叶变换
sfftI1=fftshift(fftI1);              %直流分量移到频谱中心
RR1=real(sfftI1);                    %取傅立叶变换的实部
II1=imag(sfftI1);                   %取傅立叶变换的虚部
A1=sqrt(RR1.^2+II1.^2);          %计算频谱幅值
A1=(A1-min(min(A1)))/(max(max(A1))-min(min(A1)))*255;%归一化
subplot(222),imshow(A1);xlabel('(b)原始图像的傅里叶变换图像');
f=double(I);%数据类型转换
g=fftshift(fft2(I));%图像傅里叶转换
[N1,N2]=size(I);%傅里叶变换图像尺寸
n=2;%参数赋初始值
d0=50;
n1=fix(N1/2); %数据圆整
n2=fix(N2/2); %数据圆整
for i=1:N1 %遍历图像像素
    for j=1:N2;
    d=sqrt((i-n1)^2+(j-n2)^2);
    if d==0
        h=0;
    else h=1/(1+(d/d0)^(2*n));
    end
    g(i,j)=h*g(i,j); %图像矩阵计算处理
    end
end