得出滤波器的传递函数，如何画出其幅频相应

jifang29

滤波器的传递函数如下,下一步我想绘出其幅频响应图，不知道该如何实现，高手指教。
%product.m,阶乘函数%
function y=product(a0)
a=a0;
if (a<21)
    y=factorial(a);
else
   y=1;
   for i=1:a
   y=y*i
   end
end   
%sun.m 求值表达式%
function y=sun(k0,t0,w10,w20,x)
k=k0;t=t0;w1=w10;w2=w20;
y=0;
for p=0:k   
    sum1=0;n=p;
    while( n<=2*t)  
        if (n~=t)
             s=(sin(w2*(n-t))-sin(w1*(n-t)))/(pi*(n-t));
             sum1=sum1+s*product(k+n-p)/(product(n-p)*product(k))*x.^(k+n-2*p);        
        else
             s=(w2-w1)/pi;
             sum1=sum1+s*product(k+n-p)/(product(n-p)*product(k))*x.^(k+n-2*p);
        end   
        n=n+1;
    end
    z=(-1)^(p+k)*product(k)/(product(k-p)*product(p));
    y=y+z*sum1;
end
%myfun1.m求滤波器系数%
function y=myfun(n0)
n=n0;
y=sqrt(1-0.40651^2)*sun(n,16,0,pi/4,0.40651);
end
%transfer1.m传递函数％
function y=transfer1(k0)
k=k0;y=0;syms z;
for j=0:k;
L=sqrt(1-0.40651^2)*(z.^(-1)-0.40651).^j/(1-0.40651*z.^(-1)).^(j+1);
y=y+myfun1(j)*L;
end

以上程序得出我所需要的系统传递函数，现在想要画出其幅频响应图，不知该如何画。