Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Галкин С.В. Теория функций комплексного переменного и операционное исчисление- учебное пособие для вузов (2011) (1)

.pdf
Скачиваний:
342
Добавлен:
23.03.2016
Размер:
1.65 Mб
Скачать

 

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