Главная ->  Повышение запаса устойчивости 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 [ 70 ] 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248

в:,-2=ь1

Ь а

k % )

b -2 а ~1 bo йо

{К-2 а Л

1*0 0 J

I о j

*о йо j

й-4

1*0 йо j

B\ = bl

bo 0;

b a

1*0 0

b ~i = *o

I bo 0

(8.66)

При иостуилении на вход системы единичного и.миульса 5(Г) = Y {t), изображение которого но Лапласу равно 1, изображение управляемой величины можно также представить в виде дробно-рациональной функции (8.60). Разница будет заключаться только в том, что степень числителя т возрастает на единицу, а последний коэффициент числителя = 0. Это обусловлено те.м, что получение реакции систе.мы па единичный импульс (весовой функщи) эквивалентно дифференцированию переходной функции, получающейся при действии единичного скачка. В области изображений это эквивалентно умножению на комплексную величину р.

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

(8.67)

где ж(г) - весовая функция системы но задающему или возмущающему во.здействию, x{t) - отклонение управляемой величины от нового установивпгегося состояния в переходпо.м процессе при действии едигшчпой ступеньки задающего или возмуща-юпюго воздействия.

Таким образо.м, техника вычисления оце1и<и / полностью совпадает с вычислением оценки / но формуле (8.61) или (8.65). Совпадает при этом и значение определителя Д (8.62). Отличаться в вычислениях будут определители До, А, и козф-



фшшстъ1Во,...,В или В, В , что обусловлено повьннением степени тввыра-жении (8.60) на единицу при вычислении / по сравнению со случае.м вычисления /.

Интегральная оценка / также может использоваться в безраз.мериом виде аналогично формуле (8.57):

4=-- (8.68)

Интегральные оценки In Г (или выражения квадратичных динамических ошибок) применяются для выбора структуры и параметров систем автоматического управления. При это.м наилучшими гшра.метрами считаются такие, при которых величина I или / имеет минимальное значение. Вычисление квадратичных тггеграль-

ных оценок / и / можно также производить иа основании так называемой формулы Релея, которая будет дока,запа ниже, в главе 11, Здесь она будет приведена без доказательства.

Если X(j(}}) есть изображение Фурье функции времени x(t), то существует зависимость, определяемая теоремой Парсеваля

]x\t)clt =j-]\X(Mf rfw=i ]X(7C0)f clay,

т. е. интегрирование квадрата функции но времени в пределах от пуля до бесконечности можно за.мепить интегрированием квадрата модуля изображения Фурье этой функции по всем частотам. При нахождении интегральной оценки /, соответствую-Н1ей peaKiuiH системы на входное задающее воздействие типа 1(f), изображение Фурье исследуемого отклонения х(С) = г/(°°)У(0 будет

где Ф(/со) - частотная передаточная функция замкнутой системы. Тогда

, ГгФ(0)-Ф(7Со)- ,

/ = -[ со. (8.69)

Л (if

В астатических системах и статических системах с пеедипичпой обратной связью или с масштабированием (см, § 9.3) установившееся значение у (°°) = 1 и Ф (0) = 1. Тогда формула (8.69) будет иметь вид

1=-]-Щ- dbx (8.70)

где Ф.(;(й) = 1 - ФОсо) - частотная передаточная функция замкнутой системы по ошибке.

Аналогичным обра;юм для входного задающего воздействия типа единичного импульса 5(f), изображение которого равно 1, изображение Фурье исследуемого откло-




нения х(0 = -y{t) равно частотной передаточной функции за.мкнутой систе.мы: X(j(a) = = Ф(/со) 1. В резу;нугате получаем

r = i Ф(;со)(о. (8.71)

Подобные выражения могут быть нолучетш! и для входного возмущающего воздействия, если вместо частотной передаточной функции Ф(/со) использовать передаточную функцию по возмущающему во.эдействию Ф (;(0), Недостатком интегральных оценок является то, что здесь ничем не ограничивается форма кривой переходного процесса. Оказывается, например, что три совершенно различных но форме процесса, изображешгых па рис. 8.15, имеют одно и то же значение квадратичной интегральной оценки (8.56). Часто оказывается, что выбранные по минимуму этой оценки пара.метры сисгемы соответствуют сильно колебательному процессу, ибо отмечавшееся уже при это.м стремление приблизить процесс к идеальному скачку вызывает большую скорость процесса при подходе к установивше.муся значению х = 0.

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

Если, например, подать на вход системы единичный скачок, то оишбка в переходно.м процессе определится заштрихованной частью на рис. 8.16, а. Очевидно, что величина интегральной оценки (8.56) будет те.м меньше, че.м ближе будет кривая переходного процесса к ломаной линии АОВС. Но приближение процесса к этой линии требует увеличения угла наклона кривой в начальной стадии процесса (приближение части кривой 0D к отрезку ОБ).

Увеличение же пача-аьной скорости .может вызвать значительное перерегулирование и, следовате.тьпо, малый запас устойчивости.

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

1 = 1{х-+ТЧ)Л,

(8.72)

где Г - некоторая постоянная времени.

Выясним, какой вид переходного процесса будет получаться при выборе пара.метров системы по минимуму улучшенной интегральной оценки (8.72). Для этого проделаем следуюпню преобразования:

/, = (х Txf dt - 2Гх xdt= ](х Txf dt - Тх = Txl + ](х -t- Tifdt,

ООО о ,

гдехц - начальное значение отклонения в переходном процессе.



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 [ 70 ] 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248