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

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

Для остальных начальных значений выполняются соотношения

(n-m-l) (n-m+1) , t*! 1 Г (и m) (n-m)

(H-ffl+2) (n-m+2) /2 1 2

(H-m) ..(п-т)

лп-т) (п-т)

а, г

( -2) x< *

(7.12)

Эти формулы показывают; что только при т = О, т. е. для дифференциального уравнения D(p)x{t) = bf{t) при скачке f{t), начальные значения при /, = + О соответствуют начальным значениям при t= -О.В ([юрмулах (7.12) множитель 1 и.меет ра;-.мерност ь величины f(t). Если воздействие прикладывается в виде скачка, не равно1Ч) единице, то вместо 1 следует поставить величину скачка.

Г1 р и м е р. Найдем реакцию системы на единичную ступенчатую функцию при нулевых начальных значениях, т. е. переходную функцию, если дифференциальное уравнение имеет вид

(0,05/ + ОЛр + i)x(t) = (0,5р + 1) /(0.

Для простоты примем, что пере.менная .г-является безраз.мериой величиной. Ре-1пая характеристическое уравнение 0,05/; + 0,4р +1 = 0, находим корни:

Р1.2=-7±Д = -4±;2.

Согласно заданны.м условня.м .г о = О и x q =0. Так как в данном случае и = 2 и т = 1, то начальные значения для f = +0 в соответствии с (7.11) и (7.12), будут

х,о = х.о = 0, х:о=х:о+.1 = 0 + -; = 10 с-. Определяем установившееся значение искомой координаты:

v -l = l.

Введем новую переменную z(t) = x(t) - 1. Начальные значения для новой нере.менной:

2,о = л-,о-Ху =0 - 1 = -1; 2,0 = xq = W с~\

На основании табл. 7.1 для п = 2 и случая ко.мплексиых корней имеем

2 = {В cos Xt+Csm Xt)e \

-4 + 10 2



Х\о -л-о

=X Q,

.(п-т-2) (п-т-2)

(7.13)

и вместо (7.12) для всех остальных начальных значений

4-1, о

(ff-m-l) (n-m-i)

-*+0 -0

(n-m+i) (п-т[) h . 2

(н-m-i) ( -m-l)

(п-га-1) ,.( -m-l)

(7.14)

В формулах (7.14) единица имеет размерность и.миульса ве;и1чипы/(0. т. е. раз.мерпость /(О, умноженную на вре.мя. Если воздействие тюступает в виде иеединич-иого и.мпульса, то в эти фор.мулы в.место единицы необходн.мо подставить заданную величину импульса.

Как видно из (7.14), при воздействии в виде и.мпульса, в отличие от скачка, даже для диф()ерепциального уравнения вида D(p) x(t) = bj(t), не будет равенства начальных значений для t == +0 и /, = -О, так как будет скачок в значении (п - 1 )-й производной. Скачок же первой производной х, т. е. перелом кривой, будет уже при те = п - 2, а скачок са.мой величины х - при т = п - 1.

Таким образом

2=(- cos 2t. + 3suV2c)e~.

Возврап1аясь к исходной координате, получаем переходную функцию h(t) = х(0 = 1 + z(t.) = 1 - (cos - 3 shi 2t)e-.

AHajioni4HbiM образом можно осуществить переход от неоднородного дифференциального уравнения (7.4) к уравнению без правой части при воздействии типа импульсной функции, в этом случае установившееся значение Ху = О, так как в случае /(t) = 5(0 при i -> °° будет/( =) = 0. Поэтому пет нужды вводить новую смещенную величину 2(0 и задача.заключается только в отыскании начальных значений при- = +0.

Так как едгшичная имну.тьсиая 4уикция является производной от единичного скачка 8(0 = 1() фор.мулы пересчета начальных зпачетш можно получить из (7.11) и (7.12), если заменить в них т на m + 1 и положить й , = 0. Тогда вместо (7.11) для первых п ~ т ~ 2 начальных значений получим



П р и м е р. Найдем реакцию системы иа единичный имиульс при пулевых начальных З1гачеииях, т. е. функцию веса для дифференциального у1)аниепия, приведенного в иредыдуп1ем примере.

Так как в рассматриваемом примере m = н - 1, то в соответствии с (7.14) получим

х,о=х,о+.1 = 0 + - = 10с-, о 0,05

xln = х\ +--1 -Гг -х п1 = 0+- 1 -10 = -60cl

t4.() - .l 0

1 О -

0,05 0,05

В соответствии с табл. 7.1 для и = 2 и комплексных корней

4-10-60

Окончательно получаем функцию веса

дх,о = 10с--, с- -- =±Bz=-mc-.

w(t) = x{t) = 10 (cos 2t - sin 2t)e . Этот результат можно было получить также непосредственным путем для h(l), полученного в предыдущем примере, так как w(t) = h{l).

§ 7.4. Использование преобразований Фурье, Лапласа и Карсона-Хевисайда

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

f{t) = A + Y{A sin k<S)t + coskiiit),

где k - порядок гармоники, a ш = 2л/Г - основная круговая частота. Этот ряд может быть представлен также в комплексной форме:

где комплексный коэффициент Cf, определяется выражением



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