Курсовая работа: Анализ системы автоматического регулирования угловой скорости вращения турбины

0 0 0 0.5000 2.5400

d5 = 1.7531e+041

A4 = 1.0e+010 *

0.0000 0.0000 0 0

0.2531 0.0075 0.0000 0.0000

2.5400 0.1753 0.2531 0.0075

0 0.5000 2.5400 0.1753

d4 = 1.3753e+031

A3 = 1.0e+010 *

0.0000 0.0000 0

0.2531 0.0075 0.0000

2.5400 0.1753 0.2531

d3 = 1.3757e+022

A2 = 1.0e+009 *

0.0001 0.0000

2.5308 0.0750

d2 = 5.4398e+012

A1 = 75030

d1 = 75030

Так как все определители Гурвица больше нуля, то система УСТОЙЧИВА.

Критерий Михайлова

Построим годограф Михайлова – кривую, которая описывается характеристическим вектором на комплексной плоскости. Характеристический вектор получим, подставив в характеристический полином (знаменатель передаточной функции замкнутой системы):

Программа анализа устойчивости САУ:

disp (' *** Анализ устойчивости по критерию Михайлова ***');

% знаменатель характеристического уравнения замкнутой системы

% a(7)*p^6+a(6)*p^5+a(5)*p^4+a(4)*p^3+a(3)*p^2+a(2)*p^1+a(1), гдевектор a найденранее

for i=1:1101

w(i)=i-1;% вектор значений частот