pade_app.m

来自「用MATLAB的脚本语言M语言编写的」· M 代码 · 共 10 行

function [nred,dred]=pade_app(c,nn,nd)
w=-c(nn+2:nn+nd+1)'; 
vv=[c(nn+1:-1:1)'; zeros(nd-1-nn,1)];
W=rot90(hankel(c(nn+nd:-1:nn+1),vv)); 
V=rot90(hankel(c(nn:-1:1)));
x=[1 (W\w)']; dred=x(nd+1:-1:1)/x(nd+1); 
y=[c(1) x(2:nn+1)*V'+c(2:nn+1)]; 
nred=y(nn+1:-1:1)/x(nd+1);
if nargout==1, nred=tf(nred,dred); end