如何画一阶二维全通滤波器幅频响应

[繁星的光芒]

%z1=exp(jw1) z2=exp(jw2) 程序里的i=j
%画幅频响应曲线
a1=-0.2612
a2=0.2612
a3=0.2612
[w1,w2]=meshgrid(-2:.2:2,-2:.2:2);
h=abs((a1+a2*1./exp(i*w2)+a3*1./exp(i*w1)+1./exp(i*w1).*1./exp(i*w2))/(1+a3*1./exp(i*w2)+a2*1./exp(i*w1)+a1*1./exp(i*w1).*1./exp(i*w2)));
surf(w1,w2,h)

公式在附件
这样画对吗？