Практикум по теории управления в среде MATLAB. Никульчев Е.В. - 58 стр.

%получение коэффициентов характеристического полинома
lm= poly(X);
%определение размерности
[L, N] =size(lm);
%создание матрицы с нулевыми значениями
g=zeros(N, N);
%заполнение нечетных строк матрицы Гурвица
s=0;
for i=1:2:N
     j=1;
     j=j+s;
     r=0;
     for r=2:2:N
          g(i, j)=lm(r);
         j=j+1;
     end
     s=s+1;
end
%заполнение четных строк матрицы Гурвица
s=0;
for i=2:2:N
     j=1;
     j=j+s;
     r=0;
     for r=1:2:N
          g(i, j)=lm(r);
         j=j+1;
     end
     s=s+1;
end
g=g(1:N-1, 1:N-1);
%вычисление главных миноров
minor=1;
for i=1:N-1
    dd = det(g(1:i, 1:i));
    if dd<0
      minor=0;
    end
end
%вывод результатов
if minor==0
     disp('СИСТЕМА НЕ УСТОЙЧИВА');
else
     disp('СИСТЕМА УСТОЙЧИВА');
end




                       – 58 –