Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
2ой семестр / Математика / Теория / Магазинников Л.И. - Функции комплексного переменного.pdf
Скачиваний:
92
Добавлен:
18.07.2023
Размер:
1.21 Mб
Скачать

=

(−1)k(2k)! 1

 

. По формуле (8.9) получаем

 

X

 

·

 

 

 

 

 

 

 

 

 

k=0

22k(k!)2

p2k+1

 

t2k

(t/2)2k

 

 

 

 

 

 

 

 

 

 

X

(−1)k

 

 

X

 

 

 

 

f(t) =

22k(k!)2

=

(−1)k

(k!)2

.

 

 

 

 

k=0

 

 

 

 

k=0

 

 

Справа стоит ряд для функции Бесселя нулевого порядка J0(t), сле-

1

довательно, L[J0(t)] = p .

1 + p2

8.4. Некоторые приложения операционного исчисления

8.4.1. Интегрирование линейных дифференциальных уравнений с постоянными коэффициентами

Ограничимся случаем уравнения второго порядка. Дано дифференциальное уравнение

x′′(t) + ax(t) + bx(t) = f(t),

(8.10)

где f(t) оригинал. Найти его решение x(t), удовлетворяющее начальным условиям x(0) = α, x(0) = β.

Пусть L[x(t)] = X(p), L[f(t)] = F (p). Переходим от уравнения (8.10) с неизвестным оригиналом x(t) к уравнению с неизвестным изображением X(p). Пользуясь свойством линейности и правилом дифференцирования оригинала, получаем (p2 + ap + b)X(p) − (pα+ +aα + β) = F (p). Следовательно,

X(p) =

F (p) + pα + aα + β

.

(8.11)

 

 

p2 + ap + b

 

Теперь нужно каким-либо образом найти по известному изображению искомый оригинал x(t).

Пример 8.18. Найти решение x(t) дифференциального уравнения

x′′ + x= cos t, x(0) = 2, x(0) = 0.

 

 

 

 

 

 

 

Решение. По формуле (8.11) находим

 

 

 

X(p) =

p

 

 

1

 

 

 

p3 + 2p2 + 3p + 2

 

+ 2p + 2

 

 

=

 

2

 

 

 

.

p2 + 1

p2 + p

p(p + 1)(p2 + 1)

Разлагая эту дробь на элементарные дроби, получаем

X(p) = p +

2

 

p + 1 + p2

+ 1 p2

+ 1 ,

 

 

2

 

1

 

 

1

 

 

 

1

 

 

p

 

следовательно, x(t) = 2 +

1

(et + sin t − cos t).

 

 

 

 

 

 

 

 

 

2

 

 

 

При решении задачи Коши для уравнения (8.10) с нулевыми начальными данными применяется формула Дюамеля, особенно часто

175

в тех случаях, когда изображение функции f(t) найти затруднитель-

но.

Пусть дано уравнение

 

x′′(t) + ax(t) + bx(t) = f(t)

(8.12)

и поставлены нулевые начальные условия: x(0) = x(0) = 0. С помо-

щью формулы Дюамеля решение уравнения (8.12) с нулевыми начальными данными можно выразить через решение x1(t) уравнения

x1′′(t) + ax1(t) + bx1(t) = 1, x1(0) = x1(0) = 0.

(8.13)

Покажем, как это сделать. От уравнений (8.12) и (8.13) переходим к операторным уравнениям

(

A(p)

·

X(p) = F (p),

 

A(p)

X1(p) =

1

,

(8.14)

·

p

 

 

 

 

 

 

где A(p) = p2 + ap + b, X1(p) = L[x1(t)]. Из системы (8.14) находим

X(p) = pX1(p)F (p). Учитывая, что x1(0) = 0, по формуле Дюамеля

получаем

 

t

 

 

 

x(t) = Z0

 

 

 

f(τ)x1(t − τ)dτ,

(8.15)

 

t

 

 

 

или x(t) = x1

(t)f(0) + x1(τ)f(t − τ)dτ.

1

 

 

R

 

 

 

0

 

 

 

Пример 8.19. Найти решение x(t) уравнения x′′ + x =

 

,

1 + sin2 t

x(0) = x(0) = 0.

 

 

 

Решение. Решаем уравнение x′′1 + x1 = 1, x1(0) = x1(0) = 0 опера-

 

 

1

 

 

 

 

1

 

 

 

p

 

 

 

 

 

ционным способом: X1(p) =

 

 

=

 

 

 

, следовательно,

p(1 + p2)

 

p

p2 + 1

1

cos t. По формуле (8.15) находим

 

 

 

 

 

x1(t) = t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z0

1 + sin2 τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) =

sin(t − τ)dτ

. Вычисляя определённый интеграл, получаем

 

x(t) = sin t arctg(sin t) + 4

 

cos t "ln

2

 

 

 

 

 

 

1

!#.

 

cos t − ln

 

2

 

 

 

 

2

cos t +

2

 

2

+ 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что требование нулевых начальных данных не является существенным. Путём замены искомой функции можно перейти от произвольных начальных данных к нулевым, не нарушая линейности уравнения.

По такой же схеме можно решать задачу Коши и для систем дифференциальных уравнений.

176

Пример 8.20. Найти решение системы

x(t) + y(t) − y(t) = et,

2x(t) + y(t) + 2y(t) = cos t, x(0) = y(0) = 0.

Решение. Пусть L[x(t)] = X(p), L[y(t)] = Y (p). Тогда

 

 

 

pX(p) + pY (p) − Y (p) =

p

1 1p,

 

 

 

 

 

 

 

 

 

 

 

2pX(p) + pY (p) + 2Y (p) =

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4p2 + 2

 

 

 

p2 + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из этой

системы находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X(p) = −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

p(p

4)(p

1)(p2 + 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y (p) =

 

 

 

 

 

 

 

 

p2 + p + 2

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(p − 4)(p − 1)(p2 + 1)

 

 

 

 

 

 

 

 

Функции X(p) и Y (p) разлагаем на элементарные дроби:

 

 

 

X(p) =

 

 

 

1

 

 

 

+

 

 

 

 

 

11

 

 

 

 

 

 

 

3p/17 − 5/17

 

 

 

1

,

 

 

 

 

p

 

 

1

 

34(p

 

 

4)

 

 

p2 + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2p

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

22

 

 

 

 

 

 

 

4p/17 − 1/17

 

 

 

 

 

Y (p) =

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

+

 

.

 

 

 

 

3(p

 

 

1)

 

51(p

 

 

 

4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 + 1

 

 

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) = et

 

 

11

e4t

3

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

cos t +

 

sin t −

 

 

,

 

 

 

 

 

34

17

17

2

 

 

 

 

 

 

 

 

 

 

 

 

2

et

 

 

 

 

 

22

e4t +

 

 

4

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

y(t) = −

 

+

 

 

 

 

cos t −

 

sin t.

 

 

 

3

51

17

17

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

= 0,

 

Пример 8.21. Найти решение системы: x′′y′′= 2 sin t,

если

x(0) =

1, x(0) = y(0) = y(0) = 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Обозначим L[x(t)] = X(p), L[y(t)] = Y (p). Применяя

теорему о дифференцировании оригинала, находим:

L[x′′(t)] = p2X(p)+p−1, L[y(t)] = pY (p)−1, L[y′′(t)] = p2Y (p)−p−1.

Получаем систему на неизвестные изображения X(p) и Y (p):

 

X(p) p2Y (p) = 1 − p − p

 

 

− p

,

 

 

pX(p) − Y (p) = −1,

 

2

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + p2

 

 

 

 

 

 

находим:

 

 

 

 

 

 

 

 

 

 

 

решая которую 1

 

 

 

 

p

1

 

 

 

 

 

p

X(p) =

 

 

 

, Y (p) =

 

 

+

 

.

1 + p2

1 + p2

1 + p2

1 + p2

Следовательно,

x(t) = sin t − cos t,

 

 

 

 

 

 

 

 

 

 

y(t) = sin t + cos t.

 

 

 

 

 

 

 

 

177

8.4.2. Решение интегральных уравнений типа свёртки

Уравнение вида

x

 

 

ϕ(x) = f(x) + Z0

k(x − t)ϕ(t)dt,

(8.16)

где f(x) и k(x) известные функции, ϕ(x) неизвестная функция, называется интегральным уравнением типа свёртки. Будем предполагать, что функции f(x) и k(x) являются оригиналами. При этом предположении можно довольно просто найти решение ϕ(x)

этого уравнения. Применим к обеим частям в (8.16) преобразование Лапласа и воспользуемся теоремой об умножении изображений. Обозначим L[ϕ(x)] = Φ(p), L[f(t)] = F (p), L[k(x)] = K(p). Тогда

Φ(p) = F (p) + K(p)Φ(p). Отсюда находим изображение неизвестной

функции Φ(p) =

F (p)

.

 

 

 

 

1 − K(p)

 

Пример 8.22. Решить интегральное уравнение

 

 

x

1

 

 

R

 

ϕ(x) = sin x +

(x − t)ϕ(t)dt.

 

0

 

Решение. Переходим к операторному уравнению Φ(p) = p2 + 1 +

1

 

p2

 

1

1

1

 

+

 

Φ(p). Отсюда Φ(p) =

 

=

 

 

 

+

 

,

p2

(p2 − 1)(p2 + 1)

2

p2 + 1

p2 − 1

следовательно, ϕ(x) = 12 (sin x + sh x).

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

178