Материал: бархоткин системы автоматического управления

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам
где si , i =1...N

В дальнейшем будем рассматривать вопросы устойчивости стационарных линейных динамических систем, поведение которых представляется линейным дифференциальным уравнением с постоянными коэффициентами:

an x(n) (t) + an−1x(n−1) (t) + K + a2 x¢¢(t) + a1x¢(t) + a0 x(t) =

 

= bm z(m) (t) + bm−1x(m−1) (t)K + b1z¢(t) + b0 z(t), n ³ m.

(4.1)

Будем полагать, что система является устойчивой, если она, будучи выведенной из состояния равновесия, с течением времени возвратится в исходное состояние. Как применить это определение для выяснения вопроса устойчивости САУ?

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

определении устойчивости указывается необходимость прекращения действия возмущения ( z(t) = 0 ).

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

an x(n) (t) + an−1x(n−1) (t) + K + a2 x¢¢(t) + a1x¢(t) + a0 x(t) = 0.

(4.2)

Система будет устойчивой, если с течением времени (при t → ∞ ) собственное движение, определяемое решением уравнения (4.1), будет затухающим:

lim x(t) ® 0.

(4.3)

t →∞

Таким образом, устойчивость является асимптотической, так как характер

движения системы к устойчивому положению не имеет значения, а важен предел x(t) при t → ∞ .

Известно, что общее решение данного дифференциального уравнения имеет

вид:

 

 

 

xc (t) = Cest .

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.4)

Продифференцируем n раз (4.4) и подставим результаты дифференцирования в

(4.2). После сокращения на общий множитель Cest

получим следующее алгебраическое

уравнение:

 

sn + a

 

sn−1 +K+ a

 

 

+ a s + a

 

 

 

a

n

n−1

2

s2

0

= 0.

(4.5)

 

 

 

 

1

 

 

Уравнение (4.5) называется характеристическим уравнением, корни которого

si = ai + jbi определяют характер переходного процесса в системе.

Так как характеристическое уравнение имеет n корней, то его общее решение

при отсутствии кратных корней может быть записано в виде

n

 

 

x(t) = å A esit ,

(4.6)

i=1

i

Ai - константы, значения которых определяются начальными условиями;

- корни характеристического уравнения.

Чтобы узнать, устойчива или неустойчива система, нет необходимости решать характеристическое уравнение и находить его корни. Оказывается, по свойствам корней можно определить необходимые и достаточные условия устойчивости системы.

Рассмотрим влияние корней характеристического уравнения (4.5) на устойчивость системы. В общем случае корни характеристического уравнения являются комплексными. При этом они образуют пары сопряженных корней:

101

si = ai + jbi ;

si+1 = ai jbi .

Каждая такая пара корней дает составляющую переходного процесса, равную

A e1 + jβi )t + A e1

jβi )t = eα1t (A e jβit + A ejβit ) =

i

 

α t

i +1

 

i

i+1

=

*

t + j

),

 

 

A e

1 sin(b

 

(4.7

 

1

i

i

 

 

где A1* , ji - новые постоянные интегрирования. Можно доказать, что

A* =

 

 

 

 

= arctg

A1

+ A1+1

.

A2

+ A 2

; j

i

 

 

1

1

1+1

 

 

A1

- A1+1

 

 

 

 

 

 

Составляющая (4.7) представляет собой синусоиду с амплитудой, изменяющейся во времени по экспоненте. Исследуем ее поведение при различных

значениях ai . Для первых трех случаев будем полагать, что bi ¹ 0 , т.е. на выходе системы имеет место колебательный процесс.

Первый случай: ai < 0 - колебания с течением времени затухают (рис.4.2); устойчивая система возвращается в состояние покоя.

Второй случай: ai >

неустойчивая.

X(t)

0

t

Рис.4.2. Затухающие колебания

(ai < 0, bi ¹ 0)

0 - колебания расходящиеся (рис.4.3); система

X(t)

 

0

t

 

Рис.4.3. Расходящиеся колебания

 

(ai > 0, bi ¹ 0)

 

Третий случай: ai = 0 (корни чисто мнимые) - колебания незатухающие с постоянной амплитудой (рис.4.4); состояние системы нейтральное, причем незначительное изменение параметров может привести систему к неустойчивости.

102

X(t)

 

0

t

Рис.4.4. Незатухающие колебания

 

(ai = 0, bi ¹ 0)

Четвертый случай: βi = 0 (корни действительные) - соответствующая ему составляющая переходного процесса Aieait представляет собой экспоненту, которая будет затухать или увеличиваться только в зависимости от знака αi (рис.4.5).

X(t)

 

 

1

 

2

0

t

Рис.4.5. Апериодический процесс: 1 - расходящийся (a1 > 0); 2 - сходящийся (a2 < 0)

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

Следовательно, для затухания переходного процесса необходимо, чтобы вещественные части корней характеристического уравнения были отрицательными.

lim x (t) → 0

.

Только тогда каждое слагаемое будет стремиться к нулю и в целом t →∞

Таким образом, необходимым и достаточным условием устойчивости системы

является отрицательность вещественных частей всех корней характеристического уравнения.

Корни характеристического уравнения - комплексные величины, состоящие в общем случае из вещественной и мнимой частей. Их можно представить на плоскости,

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

(рис.4.6).

103

 

jβ

 

s1

s*2

s

 

s*

 

s4

5

s3

 

α

 

 

s5

 

s2

 

 

 

s1*

 

 

 

Рис.4.6. Комплексная плоскость корней характеристического уравнения с нанесенными корнями:

( s1, s1* ) - комплексно-сопряженные с отрицательной

вещественной частью; ( s2, s2* ) - комплексно-

сопряженные с положительной вещественной частью; ( s3 , s4 ) - положительный и отрицательный

вещественные; ( s5, s5* ) - чисто мнимые

Из рисунка видно, что все корни, лежащие слева от мнимой оси, имеют отрицательные вещественные части. С учетом этой особенности расположения корней условие устойчивости системы можно сформулировать следующим образом: для устойчивой системы необходимо и достаточно, чтобы все корни характеристического уравнения находились в левой полуплоскости комплексной плоскости s.

Если хотя бы один из корней характеристического уравнения находится в правой полуплоскости, система неустойчивая.

Если хотя бы два корня расположены на мнимой оси, т.е. являются чисто мнимыми, система находится на границе устойчивости. В этом случае незначительными изменениями можно сделать систему или устойчивой, или неустойчивой.

Следовательно, для суждения об устойчивости линейной системы нет

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

Процессы в реальных САУ описываются, как правило, нелинейными дифференциальными уравнениями. Однако исследование устойчивости движения нелинейных систем "в малом" часто осуществляется по уравнениям первого приближения. Линеаризация проводится путем отбрасывания членов ряда Тейлора порядка выше первого. Возникает вопрос: как влияют принятые при линеаризации

допущения на достоверность суждений относительно устойчивости исходной нелинейной системы?

А.М. Ляпунов впервые доказал допустимость суждения об устойчивости "в малом" нелинейной системы по устойчивости линейной системы, полученной путем линеаризации исходной системы.

Теорема 1. Нелинейная система устойчива "в малом", если все корни

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

104

Теорема 2. Нелинейная система неустойчива "в малом", если хотя бы один

корень характеристического уравнения линеаризованной системы имеет положительную вещественную часть.

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

Определение устойчивости системы по виду корней характеристического уравнения возможно лишь в наиболее простых случаях. Обычно устойчивость системы

определяется косвенными методами с помощью так называемых критериев устойчивости.

4.2. Алгебраический критерий устойчивости Рауса – Гурвица

Алгебраическими называются такие критерии устойчивости, которые позволяют вынести суждение об устойчивости системы путем вычислений, производимых над коэффициентами ее характеристического уравнения. Из алгебраических критериев устойчивости наибольшее распространение получил критерий Рауса - Гурвица. В разной форме он был предложен сначала английским математиком Е. Раусом (Routh), а затем швейцарским математиком А. Гурвицем (Hurwitz) в конце XIX века.

Приведем формулировку критерия Рауса - Гурвица. Для того чтобы динамическая система была устойчива (т.е. чтобы все корни характеристического уравнения лежали в левой полуплоскости комплексной плоскости s), необходимо и

достаточно, чтобы при an > 0 определитель Рауса - Гурвица и все его диагональные миноры были больше нуля.

Таким образом, критерий Рауса - Гурвица позволяет судить о положении корней характеристического уравнения на комплексной плоскости s, а следовательно, об устойчивости системы, не решая характеристического уравнения.

Правило составления главного определителя Рауса - Гурвица, имеющего n строк

иn столбцов, где n - порядок характеристического уравнения, состоит в следующем: 1) сначала заполняется главная диагональ определителя, начиная с коэффициента an−1

изаканчивая коэффициентом a0 ;

2)затем заполняются остальные места в строках матрицы таким образом:

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

все недостающие коэффициенты, т.е. коэффициенты с индексами, большими n или меньшими нуля, заменяются нулями.

Всоответствии с указанным правилом определитель Рауса - Гурвица имеет следующий вид:

105