Главная ->  Логарифмическое определение устойчивости 

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 249 250 251 252 253 254

а на участке BD

Решение уравнения (16.12) будет

(16.13)

(16.14)

откуда получаем

е = -ГАе -kict + C.

(16.15)

Условимся для простоты отсчитывать время t от начала участка АВ (рис. 16.5, а). Тогда начальные условия будут

где пока неизвестно. Используя начальные условия, находим произвольные постоянные для уравнения (16.15):

Ci = + кгс с = ъ + TCi.

(16.16)

Аналогично для участка BD согласно (16.13), отсчитывая время t тоже от начала этого участка (рис. 16.5, б), получим решение

dQ dt

= С[е +к,с, t

6= -TiC[e +hct + C, C = %-hc, C,= ~b+T,C[. }

(16.17)

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

значениями величин С, Cg, 9а, С[, С, G. Заметим, что величины G. и б,. необходимые для определения произвольных постоянных, находятся как

значения G в конце предшествующих им участков. Поэтому, если будет

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

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

значения 6 и G, какие были в начале его А. Легко заметить, что при этом оба полупериода {АВ и BD) должны быть одинаковыми вследствие симметрии характеристики (рис. 16.4, с). Поэтому для определения автоколебаний



Рис. 16.5.



достаточно рассмотреть только один участок АВ и потребовать, чтобы

ё=-ёл. (16.18)

Обозначив период искомых автоколебаний через 2Т, а длительность участка АВ, следовательно, через Т, из (16.14) найдем

Подставляя сюда (16.18) и замечая, что из (16.16) = - kiC, получаем выражение

Ci(l+e ri)=2/cic.

(16.19)

в котором содержатся две неизвестные: Ci и Т. Величину Т (длительность, участка АВ) можно выразить из (16.15), так как известно, что в конце участка 6 = -6. Из (16.15) и (16.16) при этом находим

riCi(l -е ri) - kicT~2b.

Подставив сюда значение Ci из (16.19), получим уравнение для определения полупериода автоколебаний:

2Ti 2Ti kicTi

(16.20)

f \ 1

г / T b Щ

Это трансцендентное уравнение для Т легко решается графически (рис. 16.6) пересечением двух кривых:

Г Г ъ

Если найдено вещественное положительное значение для Г, то это свиде- Рис. 16.6. тельствует о наличии периодического решения в даннойсистеме. Чтобы доказать, что это соответствует автоколебаниям, нужно исследовать их устойчивость, т. е. показать, что в переходном процессе система ведет себя, как изображено на рис. 16.3, с, но не так, как на рис. 16.3, б. Это будет показано ниже.

Амплитуда найденных автоколебаний определяется как Эщах на участке АВ (рис. 16.5, с) путем исследования функции (16.15) на максимум обычным путем.

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

r = Oi{xi, Xz, Хп, f,g).

= OziXi, Xz, ...,Xn,f, g),

dx ~di

P = On{0Ci, Xz, Xn, f,g)

(16.21)




<с начальными условиями

1 = *10. 2 = аго. ; - по при t = о,

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

Пусть, например, в уравнениях (16.21) будет п = 3 (система третьего порядка). Переменные ж, х, Xg здесь могут иметь любой физический смысл.

Но условно их можно представить мысленно как прямоугольные координаты некоторой точки М (рис. 16.7, с).

В реальном процессе регулирования в каж-aj Ч J Дый момент времени величины х, х, Хд имеют

-- вполне определенные значения. Это соответствует

тч вполне определенному положению точки М в

f I \Мо(1-о) пространстве (рис. 16.7, с). С течением времени

--! у * В реальном процессе величины Ж], х, х опреде-

ленным образом изменяются. Это соответствует щ определенному перемещению точки М в простран-

стве по определенной траектории. Следовательно, Рис. 16.7. траектория движения точки М может служить

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

Точка М называется изображающей точкой, ее траектория называется фазовой траекторией, а пространство (х, х, х называется фазовым пространством.

Так как производные по времени от координат точки представляют проекции ее скорости v на оси координат, то дифференциальные уравнения системы в форме (16.21) представляют собой выражения для проекций скорости V изображающей точки М (рис. 16.7, с) на оси координат. Следовательно, по значениям правых частей уравнений (16.21) в каждый момент времени моншо судить о направлении движения изображающей точки М, а вместе с тем и о поведении соответствующей реальной системы в процессе регулирования.

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

Если переменных в уравнениях (16.21) будет всего две: х- и х (система второго порядка), то изображающая точка будет двигаться не в пространстве, а на плоскости {фазовая плоскость).

Если переменных будет любое число тг > 3 (система тг-го порядка), то фазовое пространство будет не трехмерным, а тг-мерным.

Итак, фазовое пространство и фазовые траектории представляют собой лишь геометрический образ динамических процессов, протекающих в системе. В этом геометрическом представлении участвуют координаты и исключено время. Фазовая траектория сама по себе дает лишь качественное представление о характере поведения системы. Чтобы определить количественно положение изображающей точки (а значит, и состояние системы) в любой момент времени, нужно найти решение заданных дифференциальных уравнений (16.21) во времени.

Если уравнения (16.21) составлены в отклонениях от установившегося состояния, то последнее характеризуется значениями ж = 2 = . . . =а; =0. Следовательно, изображением установившегося состояния системы является начало координат фазового пространства.

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



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 249 250 251 252 253 254