uchb1ap.m

来自「数字信号处理的比较全面的程序代码」· M 代码 · 共 12 行

function [b,a]=u_chb1ap(N,Rp,Omegac);
[z,p,k]=cheblap(N,Rp);
     a=real(poly(p));
   aNn=a(N+1);
     P=p*Omegac;
     a=real(poly(p));
   aNu=a(N+1);
     k=k*aNu/aNn;
     b0=k;
     B=real(poly(z));
     b=k*B;