Реальная физика   ТОЭЭ   ТЭЦ   РЭКТ   ТПОИ   ЭИ  

Теория автоматического управления

Алгебраические критерии устойчивости

1. Понятие устойчивости системы

Под устойчивостью системы понимается способность ее возвращаться к состоянию установившегося равновесия после снятия возмущения, нарушившего это равновесие. Неустойчивая система непрерывно удаляется от равновесного состояния или совершает вокруг него колебания с возрастающей амплитудой.

Устойчивость линейной системы определяется не характером возмущения, а структурой самой системы (рис.61). Говорят, что система устойчива "в малом", если определен факт наличия устойчивости, но не определены ее границы. Система устойчива "в большом", когда определены границы устойчивости и то, что реальные отклонения не выходят за эти границы.

В соответствии с классическим методом решение дифференциального уравнения ищется в виде:

 

y(t) = yвын(t) + yсв(t).

 

Здесь yсв(t) - общее решение однородного дифференциального уравнения, то есть уравнения с нулевой правой частью:

 

aoy(n) + a1y(n-1) + ... + a(n-1)y’ + a(n)y = 0.

 

Физически это означает, что все внешние воздействия сняты и система абсолютно свободна, ее движения определяются лишь собственной структурой. Поэтому решение данного уравнения называется свободной составляющей общего решения. yвын(t) - частное решение неоднородного дифференциального уравнения, под которым понимается уравнение с ненулевой правой частью. Физически это означает, что к системе приложено внешнее воздействие u(t). Поэтому вторая составляющая общего решения называется вынужденный. Она определяет вынужденный установившийся режим работы системы после окончания переходного процесса.

Можно провести аналогию между САУ и пружиной, колебания которой описываются аналогичным дифференциальным уравнением (рис.62). Оттянем пружину, а затем отпустим, предоставив ее самой себе. Пружина будет колебаться в соответствии со свободной составляющей решения уравнения, то есть характер колебаний будет определяться только структурой самой пружины. Если в момент времени t = 0 подвесить к пружине груз, то на свободные колебания наложится внешняя сила Р. После затухания колебаний, описываемых только свободной составляющей общего решения, система перейдет в новый установившийся режим, характеризуемый вынужденной составляющей yвын = y(t ). Если внешнее воздействие само будет изменяться по синусоидальному закону P = Posin(t + ), то после затухания переходного процесса система будет совершать вынужденные колебания с той же частотой, что и вынуждающая сила, то есть yвын = ymaxsin(t + y).

Каждая составляющая общего решения уравнения динамики ищется отдельно. Вынужденная составляющая ищется на основе решения уравнения статики для данной системы для времени t . Свободная составляющая представляет собой сумму из n отдельных составляющих: , где pi корни характеристического уравнения D(p) = a0pn + a1pn-1 + a2pn-2 + ... + an = 0. Корни могут быть либо вещественными pi = ai, либо попарно комплексно сопряженными pi = ai ± ji. Постоянные интегрирования Аi определяются исходя из начальных и конечных условий, подставляя в общее решение значения u, y и их производные в моменты времени t = 0 и t .

Каждому отрицательному вещественному корню соответствует экспоненциально затухающая во времени составляющая yсв(t)i, каждому положительному - экспоненциально расходящаяся, каждому нулевому корню соответствует yсв(t)i = const (рис.63). Пара комплексно сопряженных корней с отрицательной вещественной частью определяет затухающие колебания с частотой i, при положительной вещественной части - расходящиеся колебания, при нулевой - незатухающие (рис.64).

 

 

Так как после снятия возмущения yвын(t) = 0, то устойчивость системы определяется только характером свободной составляющей yсв(t). zПоэтому условие устойчивости систем по Ляпунову формулируется так: в устойчивой системе свободная составляющая решения уравнения динамики, записанному в отклонениях, должна стремиться к нулю, то есть затухать.

Исходя из расположения на комплексной плоскости корни с отрицательными вещественными частями называются левыми, с положительными - правыми (рис.65).

Поэтому условие устойчивости линейной САУ можно сформулировать следующим образом: для того, чтобы система была устойчива, необходимо и достаточно, чтобы все корни ее характеристического уравнения были левыми. Если хотя бы один корень правый, то система неустойчива. Если один из корней равен нулю (в системах, где an = 0), а остальные левые, то система находится на границе апериодической устойчивости. Если равны нулю вещественные части одной или нескольких пар комплексно сопряженных корней, то система находится на границе колебательной устойчивости.

Правила, позволяющие судить о знаках корней характеристического уравнения без его решения, называются критериями устойчивости. Их можно разделить на алгебраические (основаны на составлении по данному характеристическому уравнению по определенным правилам алгебраических выражений, по которым можно судить об устойчивости САУ) и частотные (основаны на исследовании частотных характеристик).

 

8.2. Алгебраические критерии устойчивости

 

8.2.1. Необходимое условие устойчивости

 

Характеристическое уравнение системы с помощью теоремы Виета может быть записано в виде

 

D(p) = aopn + a1pn-1 + a2pn-2 + ... + an = ao(p-p1)(p-p2)...(p-pn) = 0,

 

где p1, p2, ..., pn - корни этого уравнения. Если система устойчива, значит все корни левые, то есть вещественные части всех корней

отрицательны, что можно записать как ai = -|ai| < 0. Подставим их в уравнение:

 

a0(p + |a1|)(p + |a2| - j2)(p + |a2| + j2)... = 0.

 

Перемножая комплексно сопряженные выражения, получим:

 

a0(p + |a1|)((p + |a2|)2 + (2)2)... = 0.

После раскрытия скобок должно получиться выражение

 

a0pn + a1pn-1 + a2pn-2 + ... + an = 0.

 

Так как в скобках нет ни одного отрицательного числа, то ни один из коэффициентов a0,a1,...,an не будет отрицательным. Поэтому необходимым условием устойчивости САУ является положительность всех коэффициентов характеристического уравнения: a0 > 0, a1 > 0, ... , an > 0. В дальнейшем будем рассматривать только уравнения, где a0 > 0. В противном случае уравнение домножается на -1.

Рассмотренное условие является необходиным, но не достаточным условием. Необходимые и достаточные условия дают алгебраические  критерии Рауса и Гурвица.

 

8.2.1. Критерий Рауса

 

Раус предложил критерий устойчивости САУ в виде алгоритма, по которому заполняется специальная таблица с использованием коэффициентов характеристического уравнения:

1) в первой строке записываются коэффициенты уравнения с четными индексами в порядке их возрастания;

2) во второй строке - с нечетными;

3) остальные элементы таблицы определяется по формуле: ck,i = ck+ 1,i - 2 - rick + 1,i - 1, где ri = c1,i - 2/c1,i - 1, i 3 - номер строки, k - номер столбца.

4) Число строк таблицы Рауса на единицу больше порядка характеристического уравнения.

 

Ri

i\k

1

2

3

4

-

1

c11 = a0

c21 = a2

c31 = a4

...

-

2

c12 = a1

c22 = a3

c32 = a5

...

r3 = c11/cc12

3

c13 = c21-r3c22

c23 = c31-r3c32

c33 = c41-r3c42

...

r3 = c11/c12

4

c14 = c22-r3c23

c24 = c32-r4c33

c34 = c42-r4c43

...

...

...

...

...

...

...

 

Критерий Рауса: для того, чтобы САУ была устойчива, необходимо и достаточно, чтобы коэффициенты первого столбца таблицы Рауса c11, c12, c13,... были положительными. Если это не выполняется, то система неустойчива, а количество правых корней равно числу перемен знака в первом столбце.

Достоинство - критерий прост в использовании независимо от порядка характеристического уравнения. Он удобен для использования на ЭВМ. Его недостаток - малая наглядность, трудно судить о степени устойчивости системы, на сколько далеко отстоит она от границы устойчивости.

 

8.2.2. Критерий Гурвица

 

Гурвиц предложил другой критерий устойчивости. Из коэффициентов характеристического уравнения строится определитель Гурвица по алгоритму:

1) по главной диагонали слева направо выставляются все коэффициенты характеристического уравнения от a1 до an;

2) от каждого элемента диагонали вверх и вниз достраиваются столбцы определителя так, чтобы индексы убывали сверху вниз;

3) на место коэффициентов с индексами меньше нуля или больше n ставятся нули.

Критерий Гурвица: для того, чтобы САУ была устойчива, необходимо и достаточно, чтобы все n диагональных миноров определителя Гурвица были положительны. Эти миноры называются определителями Гурвица.

Рассмотрим примеры применения критерия Гурвица:

1) n = 1 => уравнение динамики: a0p + a1 = 0. Определитель Гурвица: = 1 = a1 > 0 при a0 > 0, то есть условиие устойчивости: a0 > 0, a1 > 0;

2) n = 2 => уравнение динамики: a0p2 + a1p + a2 = 0. Определители Гурвица: 1 = a1 > 0, D2 = a1a2 - a0a3 = a1a2 > 0, так как a3 = 0, то есть условие устойчивости: a0 > 0, a1 > 0, a2 > 0;

3) n = 3 => уравнение динамики: a0p3 + a1p2 + a2p + a3 = 0. Определители Гурвица: 1 = a1 > 0, 2 = a1a2 - a0a3 > 0, 3 = a32 > 0, условие устойчивости: a0 > 0, a1 > 0, a2 > 0, a3 > 0, a1a2 - a0a3 > 0;

Таким образом при n 2 положительность коэффициентов характеристического  уравнения является необходимым и достаточным условием устойчивости САУ. При n > 2 появляются дополнительные условия.

Критерий Гурвица применяют при n 4. При больших порядках возрастает число определителей и процесс становится трудоемким. Имеется ряд модификаций данного критерия, расширяющие его возможности.

Недостаток критерия Гурвица - малая наглядность. Достоинство - удобен для реализации на ЭВМ. Его часто используют  для  определения  влияния одного из параметров САУ на ее устойчивость. Так равенство нулю главного определителя n = ann-1 = 0 говорит о том, что система находится на границе устойчивости. При этом либо an = 0 - при выполнении остальных условий система находится на границе апериодической устойчивости, либо предпоследний минор n-1 = 0 - при положительности всех остальных миноров система находится на границе колебательной устойчивости. Параметры САУ определяют значения коэффициентов уравнения динамики, следовательно изменение любого параметра Ki влияет на значение определителя n-1. Исследуя это влияние можно найти, при каком значении Ki определитель n-1 станет равен нулю, а потом - отрицательным (рис.67). Это и будет предельное значение исследуемого параметра, после которого система становится неустойчивой.

Контрольные вопросы

  1. Что понимают под устойчивостью САУ в малом и в большом?
  2. Какой вид имеет решение уравнения динамики САУ?
  3. Как найти вынужденную составляющую решения уравнения динамики САУ?
  4. Какой вид имеет свободная составляющая решения уравнения динамики САУ?
  5. Что такое характеристическое уравнение?
  6. Какой вид имеют корни характеристического уравнения?
  7. Чем отличаются правые и левые корни характеристического уравнения?
  8. Сформулируйте условие устойчивости систем по Ляпунову.
  9. Что такое граница устойчивости?
  10.  Что такое критерии устойчивости?
  11.  Сформулируйте необходимое условие устойчивости САУ.
  12.  Сформулируйте критерий Рауса.
  13.  Сформулируйте критерий Гурвица.
  14.  В чем достоинства и недостатки алгебраических критериев устойчивости?
CAD   ТПОИ   ЭИ   Эконометрика   к библиотеке  

Знаете ли Вы, как разрешается парадокс Ольберса?
(Фотометрический парадокс, парадокс Ольберса - это один из парадоксов космологии, заключающийся в том, что во Вселенной, равномерно заполненной звёздами, яркость неба (в том числе ночного) должна быть примерно равна яркости солнечного диска. Это должно иметь место потому, что по любому направлению неба луч зрения рано или поздно упрется в поверхность звезды.
Иными словами парадос Ольберса заключается в том, что если Вселенная бесконечна, то черного неба мы не увидим, так как излучение дальних звезд будет суммироваться с излучением ближних, и небо должно иметь среднюю температуру фотосфер звезд. При поглощении света межзвездным веществом, оно будет разогреваться до температуры звездных фотосфер и излучать также ярко, как звезды. Однако в дело вступает явление "усталости света", открытое Эдвином Хабблом, который показал, что чем дальше от нас расположена галактика, тем больше становится красным свет ее излучения, то есть фотоны как бы "устают", отдают свою энергию межзвездной среде. На очень больших расстояниях галактики видны только в радиодиапазоне, так как их свет вовсе потерял энергию идя через бескрайние просторы Вселенной. Подробнее читайте в FAQ по эфирной физике.

НОВОСТИ ФОРУМА

Форум Рыцари теории эфира


Рыцари теории эфира
 13.08.2020 - 13:54: ЭКОЛОГИЯ - Ecology -> Биологическая безопасность населения - Карим_Хайдаров.
13.08.2020 - 13:47: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Проблема народного образования - Карим_Хайдаров.
13.08.2020 - 13:45: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> РАСЧЕЛОВЕЧИВАНИЕ ЧЕЛОВЕКА. КОМУ ЭТО НАДО? - Карим_Хайдаров.
13.08.2020 - 13:44: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Пламена Паскова - Карим_Хайдаров.
13.08.2020 - 06:34: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> КОМПЬЮТЕРНО-СЕТЕВАЯ БЕЗОПАСНОСТЬ ДЛЯ ВСЕХ - Карим_Хайдаров.
13.08.2020 - 06:33: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> ЗА НАМИ БЛЮДЯТ - Карим_Хайдаров.
13.08.2020 - 06:29: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Андрея Фурсова - Карим_Хайдаров.
13.08.2020 - 05:43: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> Проблема государственного терроризма - Карим_Хайдаров.
13.08.2020 - 05:42: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Вячеслава Осиевского - Карим_Хайдаров.
12.08.2020 - 21:59: СОВЕСТЬ - Conscience -> РУССКИЙ МИР - Карим_Хайдаров.
12.08.2020 - 21:56: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от схиигумена Сергия (Николая Романова) - Карим_Хайдаров.
11.08.2020 - 13:54: ЭКОЛОГИЯ - Ecology -> ЭКОЛОГИЯ ДЛЯ ВСЕХ - Карим_Хайдаров.
Bourabai Research Institution home page

Bourabai Research - Технологии XXI века Bourabai Research Institution