Галкин С.В. Теория функций комплексного переменного и операционное исчисление- учебное пособие для вузов (2011) (1)
.pdf
|
1 |
|
|
||
f (t) |
|
|
f (t)e (s i )t dt e(s i )t d |
||
|
|||||
|
2 |
|
|
||
|
|
0 |
|
||
|
1 |
|
|
1 |
|
s i |
|
|
|
|
F ( p)e pt d (s i ) |
|
F ( p)e pt dp, |
|
|||||
2 i |
2 i |
|
|||||||
|
|
|
s i |
|
|
||||
|
|
|
|
|
|
|
|||
где p s i для любого s s0 . |
|
|
|
|
|||||
|
|
|
6.4. Лемма Жордана |
|
|
||||
|
|
|
Приведем |
без |
доказательства |
лемму |
|||
|
s0 i |
Жордана (здесь она используется в дока- |
|||||||
|
|
|
зательстве общей теоремы разложения, |
||||||
|
|
|
несколько иная ее форма применена в |
||||||
R |
|
s |
разд. 3.4 для вычисления несобственных |
||||||
|
интегралов). |
|
|
|
|
||||
CR |
|
|
|
|
|
||||
|
|
Построим дугу СR – часть окружно- |
|||||||
s0 i |
|||||||||
|
сти радиусом R с центром в начале коор- |
||||||||
|
|
|
динат (рис. 6.1), |
лежащую в области |
|||||
Рис. 6.1 |
|
s s0 , отметим на ней точки с абсциссой |
|||||||
|
|
|
s0 ( s0 i , s0 i ). |
|
|||||
Лемма Жордана. Пусть функция F ( p) |
– аналитическая в по- |
||||||||
луплоскости |
Re p s s0 . Пусть F ( p) 0 |
при p , |
p CR |
||||||
равномерно по аргументу p, т. е. при p p CR (см. рис. 6.1) выполнено условие max F ( p) 0. Тогда
limR F ( p)e pt dp 0.
CR
6.5. Теоремы разложения
Общая (третья) теорема разложения. Пусть F ( p) – функ-
ция аналитическая за исключением конечного числа особых точек
p1, ..., pn (Re pk s s0 , |
k 1, 2, ..., n), удовлетворяющая доста- |
140 |
|
точным |
условиям изображения. |
Пусть F ( p) 0 при |
p , |
p CR |
равномерно по аргументу |
p. Тогда оригинал можно вы- |
|
числить по формуле |
|
|
|
n
f (t) Res pk F p e pt .
k 1
Доказательство. Пусть в область, ограниченную СR и отрезком, соединяющим точки s i , s i (см. рис. 6.1), попало m
из n особых точек. По общей теореме о вычетах (см. разд. 3.3) запишем соотношение
|
F( p)e pt dp |
s i |
|
|
|
СR |
|
s i |
m
F ( p)e pt dp 2 i Res pk (F( p)e pt ).
k 1
Устремим R к . Тогда внутрь рассматриваемой области войдут все n особых точек. К первому слагаемому может быть применена лемма Жордана (см. разд. 6.4). Его предел при R будет равен нулю. Ко второму слагаемому может быть применена теорема обращения (см. разд. 6.3). Его предел при R будет равен 2 if (t). Следовательно, в результате пре-
дельного перехода, сокращая обе части на 2 i, получаем
|
|
|
|
n |
F p e pt . |
|
|
|
|
|
|
|
f (t) Res pk |
|
|
||
|
|
|
|
k 1 |
|
|
|
|
Теорема доказана. |
|
|
|
|||||
|
|
|
|
|
cm |
|
||
Первая теорема разложения. |
Пусть F ( p) |
, тогда |
||||||
pm 1 |
||||||||
|
|
n |
|
m 1 |
|
|||
t |
|
|
|
|
||||
f (t) cn |
|
|
. |
|
|
|
||
n! |
|
|
|
|||||
n 1 |
|
|
|
|
||||
|
|
|
|
|
|
|
||
Доказательство. Найдем оригинал m (t) от изображения
cm |
. По общей теореме разложения этот оригинал равен |
|
pm 1 |
||
|
||
|
141 |
|
|
|
|
c |
e pt |
1 |
|
|
|
d m |
c e pt |
pm 1 |
|
|
c tm |
|||||
|
m |
(t) Res |
p 0 |
|
m |
|
|
lim |
p 0 |
|
|
|
m |
|
|
m |
. |
|||
|
m 1 |
|
|
|
|
|
|
|||||||||||||
|
|
p |
|
m! |
|
dp |
m |
p |
m 1 |
|
|
|
m! |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
По равномерной сходимости ряда Лорана допустимы почленное интегрирование ряда (при вычислении вычета) и почленный переход к пределу. По общей теореме разложения
n |
|
F p e pt Res p 0 |
|
|
c |
|
|
|
f (t) Res p |
|
|
|
m |
e pt |
|
||
k |
pm 1 |
|||||||
k 1 |
|
m 1 |
|
|
||||
|
|
c |
|
|
c tm |
|
|
Res p 0 |
|
m |
e pt |
|
m |
. |
|
pm 1 |
m! |
||||||
m 1 |
|
|
m 1 |
|
Вторая теорема разложения. Пусть формула |
F ( p) |
( p) |
||||||||||||
( p) |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
имеет в качестве особых точек только полюсы pk , |
|
k 1, ..., n |
||||||||||||
кратности rk . Тогда оригинал можно вычислить по формуле |
|
|||||||||||||
n |
1 |
|
|
|
d rk 1 |
|
( p) |
|
r |
|
|
|
||
f (t) |
|
|
lim p p |
|
|
|
|
p pk |
k |
e pt . |
|
|||
(r 1)! |
k dprk 1 |
( p) |
|
|||||||||||
k 1 |
|
|
|
|
|
|
|
|||||||
k |
|
|
|
|
|
|
|
|
||||||
Доказательство. Применим общую теорему разложения, вычисляя вычеты по формуле вычисления вычета в полюсе rk порядка (см. разд. 3.2). Получим искомую формулу.
142
7. РЕШЕНИЕ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ И СИСТЕМ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ
7.1.Решение дифференциальных уравнений
спостоянными коэффициентами методом операционного исчисления
При решении дифференциальных уравнений и систем используется теорема о дифференцировании оригинала и ее следствие – теорема об изображении n-й производной.
Метод решения основан на том, что преобразование Лапласа сводит дифференцирование в пространстве оригиналов к умножению на комплексное переменное p в пространстве изображений. Поэтому дифференциальное уравнение с постоянными коэффициентами в пространстве оригиналов переходит в алгебраическое уравнение в пространстве изображений. При этом учитываются и начальные условия, что удобно при решении задачи Коши.
Получив решение алгебраического уравнения в пространстве изображений, находим решение в виде некоторого изображения – функции комплексного переменного p. Остается найти соответствующий ему оригинал по свойствам преобразования Лапласа (теоремам подобия, смещения, запаздывания, дифференцирования, интегрирования) или по теоремам разложения.
Пусть задано дифференциальное уравнение с постоянными коэффициентами a0 , a1, ..., an относительно неизвестной функции
x(t) и ее производных, правая часть которого – функция f (t) – является оригиналом:
a0 x n a1x n 1 ... an 1x an x f (t).
Требуется решить задачу Коши для этого уравнения при начальных условиях
143
x(0) x0 , |
|
x10 , |
x |
|
0 |
x20 |
, ..., x |
n 1 |
(0) |
xn 10 . |
x (0) |
|
|
Применим преобразование Лапласа к обеим частям равенства: a0 pn x p pn 1x0 pn 2 x10 ... xn 10
a1 pn 1x( p) pn 2 x0 pn 3 x10 ... xn 20 ...
... an 1( px( p) x0 ) an x( p) F ( p).
Приведем коэффициенты при x( p) в левой части и перенесем члены, зависящие от начальных условий, в правую часть:
M ( p)x( p) F( p) M0 ( p) ,
где
M ( p) a0 pn a1 pn 1 ... an 1 p an
– характеристический многочлен;
M0 ( p) a 0 pn 1 a1 pn 2 ... an 1 x0
a0 pn 2 a1 pn 3 |
... an 2 x10 |
... a0 |
||||
Найдем изображение решения: |
|
|
|
|||
|
|
F ( p) M |
( p) |
|
F ( p) |
|
x( p) |
|
0 |
|
|
|
|
|
M p |
|
M ( p) |
|
||
p a1 xn 20 an xn 10 .
M0 ( p) .
M ( p)
Здесь первое слагаемое – вклад правой части дифференциального уравнения в решение, второе слагаемое – вклад начальных
условий. Если начальные условия нулевые, то M0 ( p) 0 и второе
слагаемое пропадает.
Пример. Решить дифференциальное уравнение x x 1 при начальных условиях x 0 c0 , x 0 c1.
144
Имеем
p2 x( p) c0 p c1 x( p) 1p ;
|
c p c |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
0 |
1 |
|
p |
|
p |
|
|
1 |
|
|
1 |
|
|
||
x( p) |
|
|
|
c |
c |
|
|
|
|
. |
|||||
|
p2 1 |
|
|
|
|
|
p( p2 1) |
||||||||
|
|
0 p2 1 |
1 p2 1 |
|
|
||||||||||
Первые два слагаемых соответствуют |
с0сh t c1sh t, |
оригинал |
|||||||||||||
для третьего слагаемого находим по теореме об интегрировании оригинала:
t
sh t dt ch t 1.
0
Искомое решение имеет вид
x(t) ch t 1 с0сh t c1sh t.
Пример. Решить уравнение x x et |
при начальных усло- |
|||||
виях x(0) 0. |
|
|
|
|
|
|
Имеем |
|
|
|
|
|
|
( p 1)x( p) |
1 |
; |
x( p) |
1 |
; x(t) tet |
|
p 1 |
( p 1)2 |
|||||
|
|
|
|
|||
по теореме о дифференцировании изображения.
Пример. Решить дифференциальное уравнение x x 0 при начальных условиях x(0) 2, x (0) 1.
Имеем
p2 x( p) 2 p 1 x( p) 0; x( p) |
2 p |
|
|
1 |
; |
|
p2 1 |
p2 1 |
|||||
|
|
|
||||
x(t) 2cost sin t.
145
Пример. Решить дифференциальное уравнение x 4x 3et при начальных условиях x(0) 2, x (0) 1.
Имеем
p2 x( p) 2 p 1 4x( p) |
3 |
; |
|
||||||
p 1 |
|||||||||
|
|
|
|
|
|
|
|
||
x( p) |
2 p |
1 |
|
|
|
3 |
|
||
|
|
|
|
|
. |
||||
p2 4 |
p2 4 |
( p 1)( p2 4) |
|||||||
Искомое решение равно
x(t) 2ch2t 12 sh2t 32 et sh2t.
Если свертку в записи решения вычислить трудно, то можно найти оригинал t для последнего слагаемого в изображении решения по теореме разложения:
t Res |
3e pt |
et |
3 e2t |
1 e 2t . |
|
( p 1)( p 2)( p 2) |
|||||
|
|
4 |
4 |
7.2.Решение дифференциальных уравнений
спомощью интеграла Дюамеля
Пусть задано дифференциальное уравнение
a x n a x n 1 ... a |
n 1 |
x |
a |
n |
x f (t) |
|
0 |
1 |
|
|
|
||
с нулевыми начальными условиями. |
|
|
|
|
|
|
Пусть известно решение уравнения x1 |
(t) |
при f (t) 1. Надо, ис- |
||||
пользуяэто решение, найти решение для произвольной правойчасти. Переходя к изображениям, получаем и решаем алгебраическое
уравнение:
M ( p)x( p) F ( p); M ( p)x1( p) 1p ;
146
M ( p) |
1 |
; |
x( p) |
F ( p). |
|
|
px ( p) |
px ( p) |
|
||||
|
|
|
|
|
||
|
1 |
|
1 |
|
|
|
Следовательно, |
|
|
|
|
|
|
x( p) px1( p)F ( p). |
|
|
||||
Отсюда по формуле интеграла Дюамеля |
|
|
||||
x(t) f ( 0)x (t) f (t) x (t) , |
или x(t) x |
(0) f (t) x |
(t) f (t). |
|||
1 |
1 |
|
|
1 |
1 |
|
Для вычисления выбирается одна из этих формул.
7.3. Решение систем дифференциальных уравнений методом операционного исчисления
Задана система дифференциальных уравнений. Надо решить задачу Коши:
x Ax f (t), x(0) x0 , A(n n).
Рассмотрим два способа решения: матричный и координатный. В матричном способе решения к обеим частям системы дифференциальных уравнений (к векторам) применяем преобразование Лапласа и ищем алгебраически вектор изображений
решений:
px( p) x0 Ax( p) F ( p); ( pE A)x( p) F ( p) x0 ; x( p) pE A 1 F ( p) x0 .
Теперь надо найти оригинал x(t) для вектора x( p), определяя
оригиналы для координат вектора по теоремам операционного исчисления или по теоремам разложения.
В координатном способе решения применяем преобразование Лапласа к каждому уравнению системы, получаем систему линейных алгебраических уравнений относительно изображений коор-
147
динат вектора x(t) и решаем ее. Затем находим оригиналы коор-
динат вектора решений по таблице оригиналов или по теоремам разложения.
Пример. Решить систему уравнений
x y 2e t ;
y x t2
при начальных условиях
x(0) c ;
1
y(0) c2 .
Получим решение матричным способом. Запишем систему уравнений в матричном виде
x |
|
0 |
1 x |
|
2et |
|||
|
|
|
|
|
|
|
2 |
|
y |
1 |
0 y |
|
t |
|
|||
|
|
|
||||||
и перейдем к изображениям:
p |
1 |
|
pE A |
1 |
|
|
1 |
|
p |
1 |
|
||||
pE A |
|
|
|
; |
|
|
|
|
|
|
|
|
|
; |
|
1 |
p |
|
p |
2 |
1 |
1 |
|||||||||
|
|
|
|
|
|
|
|
p |
|
||||||
|
|
2 |
|
|
|
|
|
|
|
|
|
||
F |
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
( p) p 1 |
|
; |
|
c1 |
|
; |
|||||||
|
|
2 |
|
|
|
|
0 |
|
c |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
p |
3 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
c |
|
|
||
x( p) |
|
1 |
|
p |
1 |
p 1 |
|
|
|||||||||
|
|
1 |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
p |
2 |
1 |
1 |
2 |
|
|
|
||||||||||
y( p) |
|
|
p |
|
c2 |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
p |
3 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c1 p |
|
|
c2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2 p |
|
|
|
|||||||||
|
|
|
p |
|
1 p |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
3 |
( p |
1) |
2 |
( p 1) |
||||||||||||||
|
|
2 |
1 |
p |
2 |
1 |
||||||||||||||||||||||||
|
p |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
c1 |
|
c2 p |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
2 |
|
. |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
2 |
|
||||||||
|
p 1 |
|
p 1 |
|
|
( p 1) |
( p |
1) |
|
|
p |
( p |
1) |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
148
Отсюда найдем изображения искомых функций:
x( p) |
c1 p |
|
|
c2 |
|
2 |
|
|
|
2 p |
; |
||||||
p2 1 |
|
|
|
p2 1 |
|
|
( p2 1) p3 |
( p2 1)( p 1) |
|||||||||
y( p) |
c1 |
|
|
c2 p |
|
2 |
|
2 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
. |
|||||||
p2 1 |
p2 1 |
( p 1)2 ( p 1) |
p2 ( p2 1) |
||||||||||||||
Запишем оригиналы для слагаемых в изображении функции x( p):
1 |
|
~ ch t 1 |
1 |
t |
2 |
, |
p3 ( p2 |
1) |
2 |
|
|||
|
|
|
|
здесьтрираза примененатеоремаобинтегрированииоригинала;
p |
1 |
1 |
|
|
1 |
1 |
|
|
1 |
1 |
~ |
1 |
et e t 2te t |
|
( p 1)( p 1)2 |
p 1 |
4 |
p 1 |
2 |
( p 1)2 |
4 |
||||||||
4 |
|
|
|
|
||||||||||
Найдем оригинал x(t): |
|
|
|
|
|
|
|
|
|
|||||
x(t) c ch t c |
sh t 2ch t 2 t2 te t sht. |
|||||||||||||
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
Найдеморигиналы для слагаемыхвизображениифункции y(
2 |
~ sht te t ; |
1 |
|
~ sht t. |
|
( p 1)( p 1)2 |
2 p2 ( p2 |
1) |
|||
|
|
Запишеморигинал y(t):
y(t) c1sh t c2ch t 3sh t te t 2t.
Решение системы уравнений получено. Пример. Решить систему уравнений
|
t |
; |
|
x(0) |
1; |
x |
y e |
|
|
||
y |
y 0; |
y(0) |
1. |
||
|
|
|
|
|
|
|
|
|
|
|
|
p):
149
