Главная >  Дифференцирование и интегрирование по аргументу 

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


создаст в контуре т. ток Ri{t). По примеру п. 8.2.1 можно показать, что достаточно знать функцию R, иначе говоря, импульсную реакцию контура т на контур Z, чтобы определить ток l{t), текущий в m в результате приложения к I любой электродвижущей силы, равной нулю при < 0. £С()

Вместо того чтобы разлагать функцию Efnt) E{t) на ступеньки, как на рис. 8.4, разложим ее на вертикальные полосы (рис. 8.11). £(jut) Каждая часть Е (t) будет при этом заклю- E(Zut) чена между двумя импульсами щириной в Lt. 6CAtj Обе определенные таким образом ступен- чатые функции будут стремиться к E{t) при Lt, стремящемся к нулю, при условии, если Eijt) не имеет резких скачков.

Каждая элементарная реакция будет равна

E(0)Rit), E{M)Rit - M)., E(2M)Rit - 2M), ...

Ток i(t) будет представлять собой сумму этих элементарных реакций.. Поэтому

i(t) = j E(z)R(t - z)dz= j E(t - z)Riz)dz.

(12) показывает, что импульсная реакция контура т на кон-с обобщенным взаимным сопротивлением соотношением

At М М

Рис. 8.11.

Формула тур / связана

Применение формулы обращения

8.3.22. Теорема Меллина - Фурье. Мы видели (п. 1.3.15). что единичная ступень может быть выражена интегралом

+joo

c+jco

- со

Контур интегрирования представляет собой

f-oo


Рис. 8.12.

-joo c-jca

Рис. 8.13.

вещественную ось с маленькой полуокружностью (рис. 8.12), Положим; /ш = р. Тогда

+]со

-уоо



При этом контур интегрирования представляет собой мнимую ось с маленькой полуокружностью, как показано на рис. 8.13. Он может быть заменен любым другим контуром Г, соединяющим точки - уоо и -1- joo и расположенным справа от мнимой оси. Действительно, согласно теореме Коши, так

как между мнимой осью и Г функция - не имеет особых точек, мы можем написать

+ /00

-/с (Г)

Контур г может быть, в частности, прямой, параллельной мнимой оси с положительной абсциссой с (контур Бромвича). Следовательно,

C+JCO


О т:7 г 3 -/7-7 Рис. 8.14.

dp.

(82)

C-JCO

Установив это, рассмотрим функцию h{t), изображенную на рис. 8.14. Повторив рассуждения, сделанные при выводе формул (1) - (6), мы можем разложить эту функцию на ступеньки:

/г (0) Г (О + 1 [/г (т ) - h (t: i)] Г (г - т ).

Если устремить это выражение к пределу, бесконечно уменьшая промежутки т - Cn-i. ступенчатая кривая будет стремиться к* кривой hIJ:), и тогда

+ 00

Заменим в этой формуле - -с) выражением (82). Получим

hit):

27су dt f f

pP 4-t:i

-dph {z)dz.

0 с- joo C+Joo

> = W-i- / ~f.-hi.)d.dp.

c-joo 0

Отсюда, пользуясь формулой преобразования Лапласа, найдем

C+Ja>

г yj. Id

ePtPip)

T. e.

c- joo c+joo

hit)

=w / Pip)dp-

c-joo

(83)

Мы получили формулу обращения. Ее называют формулой Меллина - Фурье.



с- JCO

Таким образом, h{t), неявно выраженная интегральным уравнением

F(/,) = J h{t)e-Pdt, . (84)

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

Формулы (83) и (84) эквивалентны. Можно применять любую из них. Нужно, однако, заметить, что формула (83) более общая, чем формула (84), и может дать результат, даже если формула Лапласа приводит к расходящемуся интегралу.

В качестве примера найдем с помощью формулы (83) некоторые выражения, уже полученные раньше с помощью формулы Лапласа,

Пример 1. Найдем оригинал -рр (п - целое положительное число).

Формула обращения (83) дает

c+jco . .

с-joo

Из первого примера п. 1.3.19 известно, что значение этого интеграла t

Мы снова находим уже известное выражение

1 t

Пример 2. Требуется найти оригинал р2д2 Применим теорему обращения. Искомый оригинал /г(0 имеет вид

c+jco c-jco

Вычисление (см. п. 1.3.12) дает й(0=-~- . а это приводит к уже

1 1

известному выражению -sin □

а р + а

8.3.23. Замечания о применении формулы обращения. Мы доказали с помощью преобразования Лапласа формулу (56):

1. .

p.+i ~ r(v + l)

Эта формула, где v-вещественное число, была доказана для случая, когда v > - 1. Если v < - 1, интеграл Лапласа не сходится. Докажем теперь формулу (56), применяя теорему обращения.

Рассмотрим функцию /(/?). Это изображение h{t), определенное формулой преобразования Карсона. Имеем

с+700



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