mATlab自编理想低通滤波器

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