Под устойчивостью системы понимается способность ее возвращаться к состоянию установившегося равновесия после снятия возмущения, нарушившего это равновесие. Неустойчивая система непрерывно удаляется от равновесного состояния или совершает вокруг него колебания с возрастающей амплитудой.
Устойчивость линейной системы определяется не характером возмущения, а структурой самой системы (рис.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). Это и будет предельное значение исследуемого параметра, после которого система становится неустойчивой.
Когда тот или иной физик использует понятие "физический вакуум", он либо не понимает абсурдности этого термина, либо лукавит, являясь скрытым или явным приверженцем релятивистской идеологии.
Понять абсурдность этого понятия легче всего обратившись к истокам его возникновения. Рождено оно было Полем Дираком в 1930-х, когда стало ясно, что отрицание эфира в чистом виде, как это делал великий математик, но посредственный физик Анри Пуанкаре, уже нельзя. Слишком много фактов противоречит этому.
Для защиты релятивизма Поль Дирак ввел афизическое и алогичное понятие отрицательной энергии, а затем и существование "моря" двух компенсирующих друг друга энергий в вакууме - положительной и отрицательной, а также "моря" компенсирующих друг друга частиц - виртуальных (то есть кажущихся) электронов и позитронов в вакууме.
Однако такая постановка является внутренне противоречивой (виртуальные частицы ненаблюдаемы и их по произволу можно считать в одном случае отсутствующими, а в другом - присутствующими) и противоречащей релятивизму (то есть отрицанию эфира, так как при наличии таких частиц в вакууме релятивизм уже просто невозможен). Подробнее читайте в FAQ по эфирной физике.