paderm.m

控制系统计算机辅助设计——MATLAB语言与应用(源代码）

function [n,d]=paderm(tau,r,k)
c(1)=1; for i=2:r+k+1, c(i)=-c(i-1)*tau/(i-1); end
Gr=pade_app(c,r,k); n=Gr.num{1}(k-r+1:end); d=Gr.den{1};