得出滤波器的传递函数，如何画出其幅频相应
滤波器的传递函数如下,下一步我想绘出其幅频响应图，不知道该如何实现，高手指教。%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
以上程序得出我所需要的系统传递函数，现在想要画出其幅频响应图，不知该如何画。
页:
[1]