iir_pade.m

来自「薛定宇高等应用数学问题的MATLAB求解例子源码」· M 代码 · 共 6 行

   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);