赫尔维兹判据.docVIP

function [ H_list,H ] = Hurwitz( X ) %==============说明============== %Author:zmf;Time:2014.11.29 %使用例子： %输入： %den=[1 12 69 198 200]; %[r2,H2]=Hurwitz(den) %输出： %r2 =12 1 0 0 % 198 69 12 1 % 0 200 198 69 % 0 0 0 200 %H2 = 1.9188e+007 %============end============== n=length(X); H_list=zeros(n-1,n-1); for i=1:n-1 for j=1:n-1 if i==j H_list(i,j)=X(i+1); else if ij2*i-j+1n+1 H_list(i,j)=X(2*i-j+1); else if ji2*i-j+10 H_list(i,j)=X(2*i-j+1); else