图像处理中的低通滤波器

%理想低通滤波器产生模糊，振铃现象
clear;clc;close all;
I=imread（‘lena.bmp‘）;
J=imread（‘lena1.bmp‘）;
subplot（231）;imshow（I）;title（‘原图像‘）;
subplot（232）;imshow（J）;title（‘含噪声图像‘）
J=double（J）;
f=fft2（J）;       %傅立叶变换
g=fftshift（f）;   %数据矩阵平衡
subplot（233）;imshow（log（abs（g））[]）color（jet）;title（‘噪声图像频谱图‘）
[mn]=size（f）;
n1=floor（m/2）;
n2=floor（n/2）;
D0=35;                    %截断半径
for i=1:m
    for j=1:n
        d=sqrt（（i-n1）^2+（j-n2）^2）;          
        if d        else h=0;
        end
        g1（ij）=h*g（ij）;
    end
end
b=ifftshift（g1）;
c=uint8（real（ifft2（b）））;
subplot（234）;
imshow（c）;title（‘理想低通滤波‘）

% %梯形低通滤波器   截断半径35-50
% D0=35;D1=50;                    %截断半径
% for i=1:m
%     for j=1:n
%         d=sqrt（（i-n1）^2+（j-n2）^2）;          
%         if d