Butterworth滤波器Matlab代码

%
% File_name:巴特沃思低通滤波器的实现
% Author:
% Date:2009年10月15日
%
Blurred=imread（‘home3.bmp‘）;
Blurred=rgb2gray（Blurred）;%I = RGB2GRAY（RGB） converts the truecolor image RGB 
                          %to the grayscale intensity image I.
figure（1）;
imshow（ Blurred）;
title（‘原图像‘）;
%设计巴特沃思低通滤波器
k=5;
[mn]=size（Blurred）;
H=zeros（mn）;
D0=30;
for u=1:m
      for v=1:n
       Duv=（（u-round（m/2））^2+（v-round（n/2））^2）^（1/2）;%ROUND（X） rounds the elements of X to the nearest integers.
       H（uv）=1/（1+（Duv/D0）^（2