TAU-Lektsia_2_3
.pdf
Для 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
