padefcn.m

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

   function [nP,dP]=padefcn(c,r,m)
   w=-c(r+2:m+r+1)'; vv=[c(r+1:-1:1)'; zeros(m-1-r,1)];
   W=rot90(hankel(c(m+r:-1:r+1),vv)); V=rot90(hankel(c(r:-1:1)));
   x=[1 (W\w)']; y=[1 x(2:r+1)*V'+c(2:r+1)];
   dP=x(m+1:-1:1)/x(m+1); nP=y(r+1:-1:1)/x(m+1);