%用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