[x,y]=meshgrid(0:31);  n=2; D0=200;
D=sqrt((x-16).^2+(y-16).^2);  % 求距離
z=1./(1+D.^(2*n)/D0); mesh(x,y,z), % 計算并繪制濾波器
axis([0,31,0,31,0,1])  % 重新設置坐標系，增大可讀性

surf(x,y,z)   % 繪制三維表面圖