chy_lingjidian.m

来自「butterworth+chebcv的matlab 程序」· M 代码 · 共 14 行

%求切比雪夫滤波器系统函数的零点、极点和增益
function[z,p,k]=chy_lingjidian(n,epi)
gamma=1/epi+sqrt(1/(epi*epi)+1);
a=(gamma^(1/n)-gamma^(-1/n))/2
b=(gamma^(1/n)+gamma^(-1/n))/2
z = [];
p =-a*sin((pi*(1:2:n-1))/(2*n))+...
    j*b*cos((pi*(1:2:n-1))/(2*n)); 
p = [p; conj(p)];
p = p(:);
if rem(n,2)==1   % n is odd
    p = [p; -1];
end
k = real(prod(-p))/sqrt(1+epi*epi);