function?s=?lpf（ADO）
s=fftshift（fft2（A））;
[ab]=size（s）;
a0=round（a/2）;
b0=round（b/2）;
d=DO;
for?i=1:a
????for?j=1:b
distance=sqrt（（i-a0）^2+（j-b0）^2）;
????????if?distance<=d?
??