

28. Интегральные оценки качества
Интегральные оценки качества являются интегралами по времени от некоторых функций координат системы (выходной координаты, сигнала ошибки) и оценивают одним числом как величину отклонения, так и время регулирования. В качестве исследуемого процесса обычно выбирается разность между установившимся процессом в системе и самой координатой. Рассмотрим замкнутую систему управления стандартной структуры, на вход которой поступает единичный ступенчатый сигнал v(t) = 1[t]. Тогда реакция системы
будет представлять собой переходную функцию
h(t) , которая в соответствии с (4.10)
определяется выражением h(t) = hy + hn (t) , где hy = KN (0) = F(0) – устано-
D(0)
.
вившаяся составляющая; hn (t) – переходная составляющая, характеризующая
переходной процесс.
Введем отклонение Dh(t) = hy - h(t) процесса h(t) от его установившего-
ся значения. Очевидно, что Dh(t) = -hn (t) .
Простейшими интегральными оценками качества являются следующие:
¥ |
|
|
|
|
|
I0 = ò Dh(t)dt , |
(7.8) |
||||
0 |
|
|
|
|
|
¥ |
|
|
|
|
|
I1 = ò |
|
Dh(t) |
|
dt , |
(7.9) |
|
|
||||
0 |
|
|
|
|
|
¥ |
|
|
|
|
|
I 2 = ò Dh2 (t)dt . |
(7.10) |
||||
0 |
|
|
|
|
|
Оценка I0 носит название линейной интегральной оценки, |
I1 – абсолют- |
ной интегральной и I 2 – квадратичной интегральной оценки.
Значение интегралов будет конечной величиной только в том случае, если
lim Dh(t) = lim[-hn (t)] = 0 , т. е. только для асимптотически устойчивых систем.
t ®¥ t ®¥
Поясним физический смысл оценок (7.8) – (7.10), для чего обратимся к рис. 7.3. Для h(t) , соответствующих кривым 1, 2, 3 (см. рис. 4.2), построены графики отклонения Dh(t) (на рис. 7.3 соответственно кривые 1, 2, 3).
Величина I0 для кривой 1 есть величина площади, ограниченной этой

кривой и координатными осями.
Очевидно, чем меньше I0 , тем меньше текущие отклонения h(t) от установившегося значения и тем меньше будет время регулирования в системе. В идеальном случае, если I0 = 0 , то время регулирования бу-
дет равно нулю. Для кривых 2, 3 в силу того, что они меняют свой знак, оценка I0 неприменима, так
Рис. 7.3
как величина интеграла может оказаться очень малой (даже равной нулю), но процессы будут затухать
медленно. Поэтому линейные интегральные оценки можно применять, если заведомо известно, что переходная составляющая имеет монотонный характер.
Для колебательных процессов обычно применяются оценки I1 , I 2 , кото-
рые имеют аналогичный смысл: чем меньше величина I, тем меньше время регулирования и меньше отклонения координаты системы от установившегося процесса.
Любые интегральные оценки носят качественный и сравнительный характер, т. е. по величине I нельзя определить, например, время регулирования или перерегулирование в системе. Но если для двух вариантов проектируемой системы окажется, что I ¢ < I ¢¢, то считается, что качественные показатели первой системы лучше, чем второй.
Наиболее просто вычисляются интегральные оценки I0 и I 2 . Пусть передаточная функция замкнутой системы имеет вид
|
|
|
F(s) = |
KN (s) |
= |
b0 sm + b1s m -1 + ... + bm |
. |
|
|
|
|
|
|
|
(7.11) |
|||||||||||||||
|
|
|
D(s) |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
a |
sn + a sn -1 |
+ ... + a |
n |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Найдем |
изображение |
|
отклонения |
Dh(t) = hy - h(t) |
с |
учетом |
того, что |
|||||||||||||||||||||||
L{h(t)} = F(s) / s , |
hy = F(0) = const и {hy } = F(0) s : |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
DH (s) = L{Dh(t)}= |
F(0) - F(s) |
. |
|
|
|
|
|
|
|
(7.12) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¥ |
|
|
|
|
¥ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Так как |
I 0 = ò Dh(t)dt = lim |
ò Dh(t)e- st dt = lim DH (s) , |
то с учетом (7.11) |
|||||||||||||||||||||||||||
имеем |
|
|
0 |
|
|
|
s ®0 0 |
|
|
|
|
|
|
|
s ®0 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I0 |
= lim |
|
b |
(a sn + ... + a |
) - a |
(b sm + ... + b |
) |
|
= |
|
|
|
|
|
||||||||||||||
|
|
|
m |
0 |
|
|
|
|
n |
|
n |
0 |
|
|
|
|
m |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
san (a0 s |
n |
+ a1s |
n -1 |
+ ... |
+ an ) |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
s®0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s ébm (a0 sn -1 + ... + an -1 ) - an (b0sm -1 + ... + bm-1)ù |
|
|
b a |
|
- a b |
|
|
|
|||||||||||||||||||||
= |
ë |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
û |
|
= |
|
m |
|
n -1 |
|
n m -1 |
. |
(7.13) |
|||
|
|
sa (a sn + a sn -1 + ... + a ) |
|
|
|
|
|
|
|
|
|
a |
2 |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
n |
0 |
1 |
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|

Квадратичная интегральная оценка I 2 может быть определена на основе
формулы Парсеваля (или Релея).
В частности, для астатических систем
|
|
I 2 |
= |
1 |
¥ |
|
Fe ( jw) |
|
2 |
dw. |
|
|
(7.14) |
||
|
|
|
|
|
|
||||||||||
|
|
|
ò |
|
|
|
|
|
|
|
|
||||
p |
|
|
|
w2 |
|
|
|
|
|||||||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Изображение отклонения DH (s) |
всегда можно представить как отношение |
||||||||||||||
двух полиномов |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
DH (s) = |
1 [F(0) - F(s)]= |
b0 s m + b1sm -1 |
+ ... + bm |
. |
(7.15) |
||||||||||
a0 s n + a1sn -1 |
|
||||||||||||||
|
s |
|
|
|
|
|
|
+ ... + an |
|
При этом оценка I 2 может быть аналитически вычислена при m < n через коэффициенты ai , bi (7.15). Выражение для вычисления I 2 имеет достаточно
сложный вид и здесь не приводится. Для наиболее распространенного случая m = n -1 приведем несколько конечных выражений для вычисления I 2 :
|
b2 |
|
|
|
b2a |
0 |
+ b2a |
2 |
|
n = 1, I2 = |
0 |
; n = 2, |
I2 |
= |
1 |
0 |
; |
||
2a0a1 |
2a0a1a2 |
|
|||||||
|
|
|
|
(7.16) |
|||||
|
|
|
|
|
|
|
|
n = 3, I2 = |
b2a |
a |
3 |
+ (b2 |
- 2b |
b |
2 |
)a |
a |
3 |
+ b2a |
a |
1 |
|
0 2 |
|
1 |
0 |
|
0 |
|
2 0 |
|
. |
|||||
|
|
|
2a0a3 (a1a2 - a0a3 ) |
|
|
|||||||||
|
|
|
|
|
|
|
Наряду с оценками I0 , I 2 употребляются и более сложные интегральные оценки, учитывавшие не только само отклонение Dh(t) , но и его производные.
При использовании интегральных оценок можно выделить два направления: анализ системы – получение оценки для заданной системы и синтез системы – минимизация оценки по каким-либо параметрам.