hurwitz.m

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

function H=hurwitz(den)
n=length(den)-1;
for i=1:n, i1=floor(i/2);
   if i==i1*2, 
      hsub1=den(1:2:end); i1=i1-1;
   else, hsub1=den(2:2:end); end
   l1=length(hsub1); l2=n-i1-l1;
   H(i,:)=[zeros(1,i1),hsub1,zeros(1,l2)];
end