def ButworthHighPasshwd0n=2 H=npemptyshape=hwdtype=float mid_x=intw2 mid_y=inth2 for y in range0h for x in range0w d=npsqrtx-mid_x2+y-mid_y2 Hyx=11+d0d+12n

function H = ButworthHighPass(h,w,d0,n) H = zeros(h,w); mid_x = floor(w/2); mid_y = floor(h/2); for y = 1:h for x = 1:w d = sqrt((x-mid_x)^2+(y-mid_y)^2); H(y,x) = 1/(1+(d0/(d+1))^(2*n)); end end end