Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лекции №1

.pdf
Скачиваний:
7
Добавлен:
17.03.2016
Размер:
972.7 Кб
Скачать

б)

C1 C2

> C0 C3 . Наличие только положительных коэффициентов в уравнении еще не является

достаточным условием, хотя это условие и необходимо.

 

 

 

 

 

 

4.

Система

 

четвертого

порядка.

Характеристическое

 

 

уравнение

C

0

λ .4 + C λ .3 + C

2

λ .2 + C

3

λ + C

4

= 0 .

Полагая, что

знаки выбраны так что

C

0

> 0 , в этом

 

 

1

 

 

 

 

 

 

 

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

2 > 0 , 3 > 0 и 4 > 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1 C3

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Главный определитель 4

 

=

 

C0 C2 C4

 

0

 

 

> 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

C1 C3

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

C0 C2 C4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Диагональные

 

 

 

 

 

 

 

 

 

 

миноры:

 

 

 

 

 

 

 

1

=

 

C

 

 

> 0 ;

 

 

 

 

 

 

 

 

2

=

 

C1 C3

 

C

 

C

2

C

0

C

3

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C0 C2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

C1 C3

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

C

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

C 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

2

C

 

= = C

 

 

 

 

 

C 2C

 

 

> 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

=

C

0

C

2

C

4

= C

 

2

3

C

C

4

 

C

0

= C

3

1

C

2

C

0

 

4

3

 

2

4

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

0

 

C1 C3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для раскрытия определителя 4

разложим его по адъюнктам первой строки:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1 C3

0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C2

C4

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C0 C4 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C0 C2 C4 0

 

 

 

 

 

 

 

 

(1)1+1

 

 

 

 

 

 

 

(1)2+1

 

 

 

 

 

 

(C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C 2 )

 

 

 

 

 

 

4

=

 

= C

 

 

 

C C

3

0

+ C

2

 

0 C

3

0

 

 

= C

 

 

2

C

3

C

4

 

C

 

 

 

 

 

 

 

 

0

C1 C3 0

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C0

C2 C4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 C2 C4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

C0 C2 C4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[C

 

(C C

 

 

 

 

 

 

 

 

 

 

 

)C 2

 

 

 

 

]=

 

 

 

 

 

C

3

C

0

C

3

C

4

= C C

2

C

3

C

4

 

C 2

 

C

2

C

0

C 2

C

4

= C

4

3

2

C

0

C

 

C

4

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

4

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= C

4

 

[C

3

 

2

C

2

C

4

]= C

4

3

 

> 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как

3

> 0 ,

то

 

 

неравенство

возможно

 

только при условии, что

 

и C4 > 0 .

Из

диагонального

минора

 

 

 

 

3

 

= C

3

2

C 2

C

4

> 0

 

 

выразим

 

C

3

 

> C

2

C

4

/

2

,

которое

 

 

 

будет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

положительно, если C1 > 0 , C4

> 0 и 2

> 0 . При

2

 

> 0 C2

> 0 , так как C2

> C0 C3 / C1 .

 

 

 

 

 

 

 

 

 

 

 

Вывод: для получения устойчивой системы необходимо соблюдение следующих условий: а)

C

0

> 0 , C > 0 , C

2

> 0 , C

3

 

> 0 , C

4

> 0 ; б) C C

2

 

> C

0

C

; в) C C

2

C

3

> C 2

C

4

 

+ C

0

 

C

 

2 .

 

 

 

 

 

 

 

 

 

 

5.

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Система

 

 

 

 

 

 

 

 

пятого

 

 

 

 

 

 

 

порядка.

 

 

 

 

 

 

 

 

Характеристическое

 

 

 

 

 

 

 

уравнение

C

0

λ .5

+ C λ .4 + C

2

λ .3

+ C

3

λ

.2 + C

4

λ + C

5

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

зона устойчивости

а)

б)

Рис. 30. Влияние изменяемого параметра δ на значения определителей: а - система устойчивая; б - система неустойчивая

41

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1 C3

C5

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C0 C2 C4

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Главный определитель Гурвица 5 =

0

C1 C3

C5

0

> 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

C0 C2 C4

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

C1 C3 C5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при условии, что C0

> 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если раскрыть определитель 5

и исследовать его диагональные миноры, то для устойчивой

системы пятого порядка необходимо соблюдение четырех условий: 1)

C0

> 0 ,

C1

> 0 ,

C2

> 0 ,

C

3

> 0 , C

4

> 0 , C

5

> 0 ; 2) C C

2

> C

0

C

; 3) C C

2

C

3

+ C

0

C C

5

> C 2 C

4

+ C

0

C 2

; 4)

 

 

 

 

 

 

 

 

 

 

1

 

 

 

3

 

 

 

1

 

 

 

1

 

1

 

 

3

 

C1 C2 C3 C4 + 2 C0 C1 C4 C5 + C0 C2 ×

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

× C

3

C

5

> C

0

C 2 C

4

+ C C

2

C

5

+ C 2

C 2

+ C

2 C 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

1

2

 

 

1

4

 

0

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Критерий Гурвица не дает возможности оценить запас устойчивости и быстроту затухания

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

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

При графическом изображении зависимостей i от исследуемого параметра δ можно определить области таких значений δ , когда все i оказываются положительными при C0 > 0

(рис. 30). При этом должно сохраняться постоянство значений других параметров, входящих в структуру определителей. В результате построения получим кривые для определителей Гурвица, число которых равно порядку характеристического уравнения системы. Расположение этих кривых относительно друг друга показывает допустимые границы изменения исследуемого параметра δ без нарушения устойчивости системы. В устойчивой системе, при всех значениях δ 1 < δ < δ 2 , всеопределители и C0 больше нуля.

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

 

 

 

 

 

 

Таблица Рауса

 

 

Таблица 3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффи-

Номер

 

Номер n-го столбца

 

 

 

 

циент r

строки

1

2

 

3

 

 

 

 

-

 

1

C1, 1 = C0

C2, 1 = C2

C3, 1 = C4

 

 

 

-

 

2

C1, 2 = C1

C2, 2 = C3

C3, 2 = C5

 

r3

=

 

C1, 1

 

3

C1, 3 = C2, 1 r3 C2, 2

C2, 3 = C3, 1 r3 C3, 2

C3, 3 = C4, 1 r3 C4, 2

 

 

C1, 2

 

 

 

 

 

 

 

 

 

 

 

r4

=

 

C1, 2

 

4

C1, 4 = C2, 2 r4 C2, 3

C2, 4 = C3, 2 r4 C3, 3

C3, 4 = C4, 2 r4 C4, 3

 

 

C1, 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ri =

C1, i2

i

C1, i = C2, i2 ri C2, i1

C2, i = C3, i2 ri C3, i1

C = C

 

r C

 

 

 

 

3, i 4, i2

i 4, i1

 

 

 

C1, i1

 

 

 

 

 

 

 

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

Как и критерий Гурвица, этот критерий представляет собой систему неравенств, составленных по особым правилам из коэффициентов характеристического уравнения замкнутой системы. Он

42

представляет собой некоторое правило (алгоритм), которое наиболее просто пояснено в табл. 3.

В первой строке таблицы записывают коэффициенты C характеристического уравнения, имеющие четный индекс (C0 , C2 , C4 , ), а во второй

строке

-

коэффициенты

характеристического

 

уравнения с нечетными индексами (C1 , C3 , C5 , )

.

В

 

последующие строки

вписывают

коэффициенты

 

Ck ,i

= Ck +1, i2

ri Ck +1, i1 ,

где

ri

= C1, i2

/ C1, i1 ;

i

 

-

 

индекс,

означающий

номер

строки

таблицы;

k

-

 

индекс, обозначающий номер столбца таблицы. Число

 

 

строк

таблицы

 

Рауса

 

равно

степени

Рис. 31. Определение

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

(n +1). После

устойчивости по критерию

заполнения

таблицы

по

ней

можно

судить

 

 

об

 

 

Найквиста: 1 - астатическая

устойчивости системы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

устойчивая САУ четвертого

Условия устойчивости Рауса: чтобы САУ была

устойчивой, необходимо и достаточно, чтобы

порядка: 2 - астатическая

коэффициенты первого столбца таблицы Рауса имели

неустовчивая САУ третьего

один и тот же знак, т. е. были бы положительными, так

порядка; 3 - статическая

как

всегда

можно

сделать

C0

> 0 :

 

C1, 1 = C0

> 0,

устойчивая САУ третьего

C1, 2

= C1

> 0,

C1, 3

> 0, , C1, n+1 > 0 .

Если не

 

все

порядка; 4 - статическая

коэффициенты первого столбца положительны, т. е.

неустовчивая САУ четвертого

если система неустойчива, то число правых корней

порядка

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

 

 

знака в первом столбце таблицы Рауса. Критерий Рауса удобен, когда заданы численные значения коэффициентов характеристического уравнения. В этом случае определение устойчивости можно выполнить быстро даже при характеристических уравнениях высокого порядка. Так как форма алгоритма, с помощью которого составляют таблицу Рауса, очень удобна для программирования на ЭВМ, то критерий Рауса широко применяют при исследовании с помощью ЭВМ влияния на устойчивость либо коэффициентов характеристического уравнения, либо отдельных параметров системы.

3.Критерий Найквиста

САУ устойчива в замкнутом состоянии, если годограф разомкнутой системы не охватывает точкус координатами (1, j 0) на комплексной плоскости (рис. 31).

Физическое толкование критерия Найквиста. Представим себе некоторую САУ (рис. 32).

При x(t) = 0 и отрицательной обратной связи x = x(t)y(t) = −y(t), т. е. обратная связь

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

при условии, что на

частоте

ω CP

сигнала

 

WP ( j ω CP ) = −1 = −e j p ,

входной

и

выходной

 

сигналы имеют одинаковые амплитуды, но

 

сдвинуты по фазе на 180° (т.е. на π

радиан). Таким

 

образом, раз возникшее колебание будет

 

существовать без изменения амплитуды. В самом

 

деле, сигнал как бы лишь дважды смещается по

 

фазе, каждый раз по 180°; результирующий сдвиг на

 

входе системы равен нулю, ослабления амплитуды

 

нет.

 

 

 

 

Очевидно, годограф разомкнутой системы на

Рис. 32. Физическое представленне

частоте ω CP пересекает ось действительных величин

критерия Найквиста

43

в точке (1, j 0). Когда модуль передаточной функции на частоте, где фазовый сдвиг равен 180°,

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

если модуль частотной передаточной функции A(ω ) < 1. Это следует из того, что сигнал,

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

Рассмотрим несколько годографов, характер переходного процесса и фрагменты ЛАЧХ и

ЛФЧХ (рис. 33).

 

 

 

1. Годограф не охватывает точку (1, j 0) (рис. 39,а), запас по фазе ϕ З = ϕ + π

> 0 , так как

π < ϕ < 0 . Система устойчива,

колебательный процесс затухающий. На частоте,

при которой

модуль комплексного коэффициента усиления равен единице, фазовый угол ϕ

< π ,

т. е. имеется

некоторый запас по фазе. Там,

где ЛФЧХ проходит через значение ϕ = −π

, она находится в

области отрицательных значений; ординатуЛАЧХ, где ϕ = −π , называют запасом по амплитуде. 2. Годограф проходит через точку (1, j 0) (рис. 33,б), ϕ З = ϕ + π = 0 , так как ϕ = −π .

Возникший колебательный процесс может существовать сколь угодно долго, амплитуда колебаний сохраняет свое значение, определяемое высотой начального импульса неизменным. Такую систему называют консервативной. ЛАЧХ

пересекает

ось

частот

там, где

ϕ = −π

 

, т.

е.

запас

по

фазе и

амплитуде отсутствует.

 

 

 

3. Годограф охватывает точку

(1, j 0)

 

(рис.

 

33,в),

ϕ З

= ϕ

+ π < 0 ,

 

 

так

 

как

π

< ϕ

< −2 π .

 

 

 

Система

неустойчива,

процесс

 

носит

расходящийся

 

 

характер,

теоретически амплитуда колебаний

способна расти до бесконечности.

На частоте среза, где модуль равен

единице и ЛАЧХ пересекает ось

частот

 

(L = 0)

,

фазовый

угол

ϕ З

= −π

ϕ

, т.

е.

запас по фазе

отрицательный.

 

 

Нет

и

положительного

 

запаса

по

амплитуде,

так

как

при

ϕ

= −π

LЗ

> 0 .

 

 

 

 

 

 

 

Примечания: 1. Для нахождения

частоты среза, на которой модуль

равен 1, проводят окружность с

центром в начале координат и

радиусом, равным 1. Пересечение

этой

окружности

с

годографом

определяет частоту среза - частоту,

Рис. 33. Характеристики амплитудно-

где модуль равен единице. 2. Если

фазовые, переходные, ЛАЧХ и ЛФЧХ САУ: а

имеется запас по фазе, но он мал,

то возникщие колебания затухают

– устойчивая; б - консервативная; в –

медленнее, чем в случае, когда

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

 

44

запас по фазе достаточен (рис. 34). Принято считать, что для удовлетворительно работающих систем автоматического регулировании запас по

фазе ϕ З

= 35 45O , а

запас

по

амплитуде

LЗ ≥ 10 дБ.

И,

как

дополнительное условие, ЛАЧХ должна пересекать ось частот с наклоном 20 дБ/дек.

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

Рис. 34. Влияние угла запаса по фазе на переходной процесс

45

y(t)

12. КАЧЕСТВО ПРОЦЕССА АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ

Показателями качества функционирования САУ называют количественные величины, характеризующие поведение системы в переходном процессе при поступлении на её вход единичного ступенчатого воздействия.

Пусть задачей САУ является обеспечение равенства управляемой величины заданной величине yУСТ при действии возмущения z(t). Любая материальная система

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

ε (t) = yУСТ y(t).

10.1.Методы анализа качества переходного процесса

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

Из прямых показателей качества переходного процесса наиболее часто используют следующие величины (рис. 35).

1. Время регулирования tP , в течение которого, начиная с момента приложения воздействия на систему, отклонения управляемой величины y от её установившегося значения yУСТ будут больше заданной величины ε (оценка быстродействия системы). Обычно принимают, что по истечении времени tP отклонение управляемой величины от установившегося значения должно быть не более

ε= 5% .

2.Перерегулирование σ - максимальное отклонение ymax управляемой

величины от установившегося значения, выраженное в процентах от yУСТ

(характеризует

колебательность

переходного

процесса).

Абсолютное

значение

ymax

определяют из кривой переходного процесса ymax = ymax yУСТ . Соответственно перерегулирование

σ= [( ymax yУСТ ) / yУСТ ] 100% .

3.Установившаяся ошибка -

отклонение установившегося значения выходной величины y(t)

от заданного значения.

4. Время достижения первого максимума tmax .

Рис. 35. Качественные показатели переходного процесса

46

OB = WP ( j ω C ), т.е. он представляет собой комплексный

5. Время нарастания переходного процесса tH - минимальное время, за которое переходная характеристика системы пересекает уровень установившегося значения.

6.Частота колебаний ω = 2 π / T где T - период колебаний.

7.Коэффициент колебательности М - отношение модуля передаточной функции замкнутой системы при ω C к модулю передаточной функции при ω = 0 , т.

е. M = WЗ ( j ω C ) /WЗ ( j 0) .

На рис. 36 представлен фрагмент годографа

некоторой

разомкнутой устойчивой системы.

Очевидно, что отрезки OA , OB , AB - векторные

величины,

причем OA + AB = OB и AB = OB OA .

Вектор OA = −1, а вектор

Рис. 36. Фрагмент годографа разом разомкнутой устойчивой системы

коэффициент

усиления

разомкнутой системы на

частоте

среза.

Тогда вектор

AB = 1 +WP ( j ω C ). Если

рассматриваемую САУ замкнуть,

тогда

комплексный

коэффициент

замкнутой системы на частоте среза

WЗ ( j ω ) = WP ( j ω )/[1+ WP ( j ω C )].

Выражение в знаменателе - вектор АВ, который оказывается тем меньше, чем меньше запас по фазе. Очевидно и то обстоятельство, что чем меньше знаменатель, тем больше WЗ ( j ω ).

Если построить амплитудные характеристики для различных 1 +WP ( j ω C ), то

они будут иметь вид, представленный на рис. 37. Пик характеристики тем выше, чем меньше 1 +WP ( j ω C ). Чем выше пик, тем сильнее выражены колебательные свойства

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

Колебательные свойства системы оценивают по отношению модуля частотной передаточной функции замкнутой системы на частоте среза к модулю комплексного

коэффициента

усиления

на частоте ω

= 0 .

Для

удовлетворительного

протекания

переходного процесса система должна иметь M = 1,2 1,4 (рис. 37). При больших

значениях M

колебательный процесс

затухает медленно. При

значении

коэффициента

M < 1,2

колебательный

процесс

становится апериодическим и

быстродействие системы, как правило, низкое.

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

10.2. Косвенный метод определения показателей качества переходного процесса по распределению корней характеристического уравнения

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

47

Ряс. 37. Влияние угла запаса на коэффициент колебательности

Рассмотрим

характеристическое

уравнение

C0 λ .n + C1 λ .n1 + … …+ Cn1 λ

+ Cn

= 0 .

Если λ .1 ,

 

λ .2 ,

…, λ ..n

-

корни этого

уравнения, а

переменная

y

 

 

 

 

 

 

 

n

 

 

характеризует процесс

управления,

то

y(t) = Ai

eλ .i t .

 

 

 

 

 

 

 

i=1

 

 

Требуется написать условия, при которых величина

y

за

время

регулирования

tP

стала

бы

равной 1/ m

уста-

новившегося значения. Таким образом, косвенно задается время стабилизации переходного процесса.

В этом случае все корни λ .1 , λ .2 , …, λ ..n характеристического уравнения должны удовлетворять не только условиям устойчивости, но и иметь отрицательную вещественную часть по абсолютному значению не меньше α . Величину α , определяемую на комплексной плоскости

корней как расстояние от мнимой оси до ближайшего к ней корня, находят из соотношения

1/ m = eα t p , откуда, логарифмируя, получаем α

= ln (m)/ tP . Следовательно,

чтобы отклонение

управляемого параметра уменьшилось за время

tP в m раз, необходимо,

чтобы все корни

характеристического уравнения находились в левой полуплоскости на расстоянии не меньше чем ln (m) / tP от мнимой оси.

Введем в

характеристическое уравнение

новую переменную z = λ

+ ln (m) / tP .

Для

переменной z мнимой является, очевидно, ось,

сдвинутая влево на величину

ln (m) / tP .

Тогда

преобразованное

характеристическое

уравнение

будет

иметь

вид

C0 (z ln (m) / tP )n + C1 (z ln (m) / tP )n1 + …+ Cn1 (z ln (m) / tP )+ Cn = 0 .

Каждая степень разности в данном уравнении может быть раскрыта в следующем виде:

z ln m n

= zn n zn1 ln m +

n (n 1)

zn2 ln m 2

2!

 

 

tP

 

 

tP

 

tP

 

 

n (n

1) (n 2)

zn3 ×

ln m + …+

(1)

ln m

n .

 

 

 

 

 

 

 

3!

 

 

 

 

 

 

 

 

 

tP

 

 

tP

 

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

Применение любого критерия устойчивости (например, критерия Гурвица) к видоизмененному характеристическому уравнению дает возможность установить минимальное значение отрицательной вещественной части у наименее удаленного от мнимой оси корня характеристического уравнения, что дает возможность судить о времени затухания процесса или «степени устойчивости» α .

Распределение корней на комплексной плоскости можно характеризовать не одной, а

Рис. 38. Влияние коэффициента колебательности на переходной

48

несколькими величинами: расстоянием α ближайшего корня от мнимой оси и углом ϕ , в

который вписываются наиболее отдаленные от мнимой оси комплексные корни (рис. 45). Величину cosϕ = β называют колебательностью системы (коэффициентом затухания колебаний).

Угол ϕ для уменьшения колебательности следует уменьшить. Его значение характеризует время регулирования: чем больше α , тем меньше время регулирования, а значение β - колебательность системы; чем меньше β , тем более система склонна к колебаниям. Следовательно, для

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

10.3.Метод приближенного аналитического определения корней характеристического уравнения

Если

λ .1 ,

λ .2 ,

…, λ ..n

являются

корнями

характеристического

уравнения

C0 λ .n + C1 λ .n1 + … + Cn1 λ + Cn

= 0 , то для

нахождения

переходного процесса

в системе

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

Имеем

характеристическое

уравнение

вида

C0

λ .n + C1 λ .n1 + … …+ Cn1 λ

+ Cn = 0 . В

уравнении,

образованном из

трех

последних

членов

характеристического

уравнения,

Cn2 λ .2 + Cn1 λ + Cn = 0 определяем

корни

λ .01

и

λ .02 .

Если эти корни

оказываются

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

Вычисление вещественного корня. Процесс вычисления первого вещественного корня

состоит в следующем: задают первое приближенное

значение

искомого

корня

в виде

λ .1.* = −Cn / Cn1 и делят характеристическое уравнение на

разность

(λ λ .1.* ) до

тех пор,

пока в

остатке не окажется двучлен вида Cn*1 λ + Cn , который нельзя разделить без остатка на разность

(λ λ ..*1 ).

 

 

 

 

В

качестве

второго приближения для первого искомого корня берут значение λ ..*1* ,

определяемое как

отношение вида -

Cn / Cn*1 . Затем характеристическое

уравнение делят на

разность

(λ λ ..*1* )

до тех пор, пока

не останется двучлен вида Cn**1 λ

+ Cn . Берут третье

приближение корня λ ..*1** , определяемое выражением λ ..*1** = −C1 / Cn**1 . Чаще всего достаточно двухтрех приближений для того, чтобы остаток от деления характеристического уравнения на соответствующую разность (λ λ .1 ) был бы близок к нулю. Это означает, что первый искомый корень λ .1 определен, после чего степень характеристического уравнения понижается на единицу.

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

Вычисление комплексной пары кор-

ней. В этом случае в качестве первого приближения берут трехчлен вида

λ .2 + (Cn1 / Cn2 ) λ + Cn / Cn2

и характеристическое уравнение делят на этот трехчлен до тех пор, пока в остатке не окажется трехчлен

Рис. 39. Области расположения корней с заданными значениями α и β

49

вида C* n2 λ .2

+ C* n1 λ + Cn , который не делится без остатка на трехчлен первого

приближения.

Затем

берут

второе

приближение

λ .2

+ (C *n1 / C* n2 ) λ + Cn / C* n2

и

на него

снова

делят

характеристическое

уравнение до

тех

пор, пока в остатке не

получится трехчлен

вида

C** n2

λ .2 + C **n1 λ + Cn . Третье принижение будет λ .2 + (C **n1 / C** n2 ) λ

+ Cn / C**n2 .

 

Обычно двух, трех приближений достаточно для получения корректного результата. Получив удовлетворительное приближение вида λ .2 + A λ + B , определяют первые два корня

исходного характеристического уравнения, т. е. λ .1, 2 = −A / 2 ± A2 / 4 B . После этого

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

(Нелинейные системы)

50