paderm.m

来自「包含大量遗传算法程序」· M 代码 · 共 4 行

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