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

TAU-Lektsia_2_3

.pdf
Скачиваний:
47
Добавлен:
01.06.2015
Размер:
1 Мб
Скачать
☆

Для n = 1, 2, 3 имеем:

 

 

n = 1: с0 = а0 – a1, с1 = а0 + a1;

 

(8.10а)

n = 2: с0 = а0 – a1 + a2, с1 = 2( а0 – a2), с2

= а0 + a1 + a2;

(8.10б)

n = 3: с0 = а0 – a1 + a2 – a3, с1 = 3( а0 + a3) – a1 – a2,

 

с2 = 3( а0 – a3) + a1 – a2,

с3 = а0 + a1 + a2 + a3.

(8.10в)

В силу утверждения 8.1, если корни исходного характеристического уравнения располагаются внутри единичного круга, то корни преобразованного характеристического уравнения (8.9) располагаются в левой полуплоскости.

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

11

Пример 8.3. Характеристический полином дискретной системы управления имеет вид

Q*(z) = z3 – 0,1 z2 – 0,46 z – 0,08.

Определить ее устойчивость.

Решение. В данном случае

а0 = 1, а1 = – 0,1, а2 = – 0,46, а3 = – 0,08

и в соответствии с (8.10в) коэффициенты преобразованного уравнения равны

с0 = а0 – a1 + a2 – a3 =1 + 0,1 – 0,46 + 0,08 = 0,72,

с1 = 3( а0 + a3) – a1 – a2 = 3(1 – 0,08) + 0,1 + 0,46 = 3,32, с2 = 3( а0 – a3) + a1 – a2 = 3(1 + 0,08) – 0,1 + 0,46 = 3,6, с3 = а0 + a1 + a2 + a3 = 1 – 0,1 – 0,46 – 0,08 = 0,36.

Необходимое условие устойчивости выполняется: все коэффициенты преобразованного характеристического уравнения больше нуля. Определитель Гурвица 2-го порядка есть

2 c1c2 c0c3 3,32 3,6 0,72 0,36 11,69 0

12

Следовательно, система устойчива.

8.2.3. Критерий устойчивости Джури

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

d00 = a0

d01 = a1

…

d0, n – 1 = a n – 1

d0n = a n

d10

d11

…

d1, n – 1

 

………………………………………………………………………………

dn – 1,0

dn – 1,1

dn0

13

Клетки нулевой строки заполняются коэффициентами характеристического уравнения в порядке возрастания нижних индексов: d0k = ak (k = 0,1,..., n).

Элементы первой строки d1k (k = 0,1,..., n – 1) вычисляются следующим образом:

Выписываются элементы нулевой строки и под ними те же элементы в обратном порядке. Из элементов верхней строки вычитаются соответствующие элементы нижней строки, умноженные на отношение последних элементов двух выписанных строк α1 = d0n / d00 = an / a0:

 

a0

a1

…

an – 1

an

 

an

 

 

 

 

 

1

 

–

 

 

 

 

a0

an

an – 1

…

a1

a0

 

d10 = a0 – α1an d11 = a1 – α1an – 1 … d1, n – 1 = a n – 1 – α1a1

Последняя разность обращается в нуль, и она отбрасывается. Поэтому первая строка содержит n элементов: на один элемент меньше, чем нулевая строка.

14

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

(k – 1)-й строки и под ними те же элементы в обратном порядке. Из элементов верхней строки вычитаются соответствующие элементы нижней строки, умноженные на отношение последних элементов выписанных двух строк. Последняя разность, обращающаяся в нуль, отбрасывается. Формула для вычисления i-го элемента k-й строки (k = 1, 2, ... , n) имеет вид

cki ck 1,i k ck 1,n k i 1 , k = 1, 2, … , n, i = 0,1, …, n – k. (8.11)

Критерий Джури. Для того чтобы все нули (корни) характеристического полинома

Q*(z) = a0zn + a1zn – 1 + ... + an

находились внутри единичного круга, необходимо и достаточно, чтобы при a0 > 0 все элементы нулевого столбца таблицы Джури были положительны: di0 > 0,

i = 1, 2,..., n.

Если все элементы нулевого столбца, кроме последнего, положительны: di0 > 0, i = 1, 2,..., n – 1, то положительность последнего элемента, т.е. условие dn0 > 0, эквивалентно необходимому условию устойчивости (8.5). Поэтому если необходимое условие выполняется, то последний элемент dn0 можно не

вычислять.

15

Пример 8.4. Характеристический полином дискретной системы управления имеет вид

Q* (z) z4 0,7z3 0,4z2 0,05z 0,1

Исследовать устойчивость данной системы.

Решение. Сначала проверим необходимое условие устойчивости:

Q* (1) 1 0,7 0,4 0,05 0,1 0,05 0 ( 1)4 Q* ( 1) 1 0,7 0,4 0,05 0,1 1,35 0

Необходимое условие устойчивости выполняется. Вычислим элементы таблицы Джури.

Для нулевой строки имеем

d00 = 1, d01 = – 0,7, d02 = – 0,4, d03 = 0,05, d04 = 0,1.

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

16

1

– 0,7

– 0,4

0,05

0,1

 

 

0,1

0,05

– 0,4

 

– 0,7

1

 

α1 = 0,1

 

 

0,99

– 0,705

– 0,36

0,12

 

 

 

 

0,12

– 0,36

– 0,705

0,99

 

α2 ≈ 0,121

 

 

 

0,975

– 0,69

– 0,275

 

 

 

 

 

 

– 0,275

– 0,69

0,975

 

α3 ≈ –0,282

 

 

 

 

 

 

 

 

 

0,897

–0,883

 

 

 

 

 

 

 

Элементы нулевого столбца (кроме последнего) равны соответственно

 

 

d00 = а0 = 1,

d10 = 0,99,

d20 = 0,975,

d30 = 0,897,

 

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

17

Разделы изучаемые студентами самостоятельно:

• Частотный критерий устойчивости

Ким Д.П. Теория автоматического управления. Т.1. Линейные системы. М.: ФИЗМАТЛИТ, 2003. (стр. 242-247)

18

8.4. Влияние квантования по времени на устойчивость

Рассмотрим влияние дискретизации по времени на устойчивость на примере

АИМ-системы, состоящей из фиксатора нулевого порядка и непрерывной части

сначала с передаточной функцией

 

k

 

WH (s)

 

 

 

 

 

T1s 1

 

 

 

а затем с

 

 

k

WH

(s)

 

 

 

 

 

(T1s 1)(T2 s 1)

 

 

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

В первом случае передаточная функция приведенной непрерывной части есть

W (s)

k(1 e T0s )

s(T1s 1)

П

 

19

Дискретная передаточная функция разомкнутой системы имеет вид

W * (z) WП* (z) ZT {WП

(s)}

k(z 1)

 

1

 

 

ZT

 

 

 

 

 

z

s(T1s 1)

 

 

 

k(z 1)

 

z

 

z

 

k(1 eT0 /T1

)

 

 

 

 

 

 

 

 

 

z

 

 

 

z eT0

 

 

 

z 1

 

z eT0 / T1

 

/ T1

 

Отсюда для характеристического уравнения замкнутой системы имеем z e T0 /T1 k(1 e T0 /T1 ) 0

Коэффициенты преобразованного характеристического уравнения (см. (8.10а))

равны

c0 a0 a1 1 e T0 /T1 k(1 e T0 /T1 ) c1 a0 a1 1 e T0 /T1 k(1 e T0 /T1 )

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

c0 1 e T0 /T1 k(1 e T0 /T1 c1 1 e T0 /T1 k(1 e T0 /T1

) 0

) 0

20

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]