求理想高通滤波算法  哪位给写一下啊！  谢谢了~~

可以直接用工具箱操作,那个就很简单了.!

截止频率d0选成50后 图片没什么大的变化，只是变成绿色 不知怎么回事啊！ 帮我看一下啊！ 谢谢～～

I=imread('4.jpg');
imshow(I);
f=double(I);
g=fft2(f);
g=fftshift(g);
[M,N]=size(g);
nn=2; % 理想高通通滤波器
d0=50;
m=fix(M/2); n=fix(N/2);
for i=1:M
for j=1:N
d=sqrt((i-m)^2+(j-n)^2);
if (d<=d0)
h=0;
else

h=1; % 计算低通滤波器传递函数
end
g(i,j)=h*g(i,j);
end
end
g=ifftshift(g);
J2=ifft2(g);
J3=uint8(real(J2));
Y=imresize(J3,[200 200] )
figure,imshow(Y);