%理想低通濾波器產生模糊，振鈴現象
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