4-27.m

matlab数字信号与应用源代码4

%例程4－27实现Butterworth低通滤波器

I = imread('eight.tif');
figure, imshow(I),
D = imnoise(I,'gaussian');
figure, imshow(D)
D=double(D);
F=fft2(D);                      %傅立叶变换
F=fftshift(F);                   %转换数据矩阵
[N1,N2]=size(F);
n=2;
d0=40;
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/(1+(d/d0)^(2*n));    %计算低通转换函数
       FD(i,j)=h*F(i,j);        %低通滤波
end
end
FD=ifftshift(FD);
FD=ifft2(FD);
FD=uint8(real(FD));
figure,imshow(FD)