hurwitz.m

来自「matlab在控制系统中的应用（内附使用方法及每章内容、程序）」· M 代码 · 共 14 行

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