@借鉴https://blog.csdn.net/charlene_bo/article/details/70599105,在此谢谢
f=imread('coulomb.bmp');
h=double(f(:,:,1));
[M,N] = size(h);
P=2*M;Q=2*N;
D=zeros(P,Q);
D(1:M,1:N)=h;
figure,imshow(D,[0,255]);
u = 0:(P-1);
v = 0:(Q-1);
[V,U] = meshgrid(v,u);
j=fft2(D.*(-1).^(U+V));
figure;imshow(abs(j),[0,100000]);
D0 = 10;
n = 6;
v0 = 113;
u0 = 128;
v1 = 146;
u1 = 128;
D1 = sqrt((V-v0).^2+(U-u0).^2);
D2 = sqrt((V-v1).^2+(U-u1).^2);
H = 1./(1+(D0^2./(D1.*D2)).^n);
figure;imshow(H,[])
% 滤波
G = j.*H;
g = real(ifft2(G));
L=g(1:M,1:N);
figure;
imshow(L,[0,255]);