function result = myfilter (f , H , D0);

%根據輸入圖象f求輸出圖象result
F = fft2(f);
F1 = movetocenter(F);
[x,y] = size(F1);
%G = F.*H;

G =zeros(x,y);
for i = 1:x
    for j = 1:y
        G(i,j) = F1(i,j) * H(i,j);
    end
end

result = ifft2(G);

result = uint8( abs(result) );