iir_pade.m

高等引用数学的matlab求解的配套程序

   function H=iir_pade(r0,a,n,T)
   syms x r b; aas=((1-x^2)/(1+b*x)^2)^r; b0=(3+a-2*sqrt(3*a))/(3-a);
   aas=subs(aas,{r,b},{r0,b0}); c=taylor(aas,x,2*n);
   c=sym2poly(c); c=c(end:-1:1); [n,d]=padefcn(c,n-1,n);
   H=tf(n(end:-1:1)*(6*b0/T/(3-a))^r0,d(end:-1:1),'Variable','z^-1','Ts',T);