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

2294

.pdf
Скачиваний:
1
Добавлен:
15.11.2022
Размер:
1.37 Mб
Скачать

Для функции a справедливо равенство: a 1 a a , которое позволяет свести вычисление функции a при a 1

к вычислению функции a при

a 0,1 . В частности, при

 

 

 

n 1 n!,

1 2

 

 

 

натуральных значениях a ( a n ):

 

,

3 2

 

2.

 

 

 

 

 

 

 

 

 

 

 

Найдем тепеть изображение

функции

f t ta ,

где

a 1 – любое действительное число.

 

 

 

 

 

 

По определению F p f t e pt dt . Считая

p 0 дей-

0

 

ствительной переменной, сделаем в интеграле замену пере-

менной: t1 pt , t t1 p , dt dt1

 

p . В результате получим:

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a 1

 

 

 

F p

 

t1ae t1 dt1

 

 

 

 

 

 

 

 

 

.

 

 

pa 1

 

 

pa 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, если a 1, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t a

=

a 1

.

 

 

 

 

 

 

 

 

 

 

 

 

(6.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

a 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В частности, если a n – целое положительное число, то

t n =

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, что совпадает с результатом (6.20).

 

pn 1

 

 

При a 1 2 и a 1 2 получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

,

 

 

t =

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 p

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

p

 

Эти результаты принято записывать в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

=

 

1

,

 

2

 

t

 

=

 

 

 

 

1

 

 

.

(6.25)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

p

 

 

 

 

p

p

81

7. ОБРАТНОЕ ПРЕОБРАЗОВАНИЕ ЛАПЛАСА

7.1. Формула Римана–Меллина

Рассмотрим теперь обратную задачу. Дана функция F p комплексного переменного p. Требуется выяснить, при

каких условиях функция F p является преобразованием Ла-

пласа для некоторой функции f t и как найти оригинал f t .

Для решения этой задачи сравним преобразование Лапласа с преобразованием Фурье. Если функция f t абсолютно интегрируема на всей числовой оси t и на каждом отрезке

T,T

конечной длины кусочно-непрерывна и

кусочно-

монотонна, то для нее существует преобразование Фурье:

 

 

 

1

 

 

 

 

 

 

F

 

 

f t e i t dt .

(7.1)

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция f t по функции F находится с помощью

обратного преобразования Фурье:

 

 

 

 

 

1

 

 

 

 

 

 

f t

 

 

 

F ei t d .

(7.2)

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть f t – оригинал,

 

 

F p – преобразование Лапласа

функции

f t . Т.к. f t 0

при t 0 и

p i ,

то форму-

лу для

определения

F p можно

записать

в виде:

F p f t e pt dt ;

82

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

F i f t e

i t

 

 

 

 

 

 

 

 

 

t e

 

 

 

 

 

 

 

 

 

t

i t

 

 

dt

 

 

 

 

2 f

e

 

dt .

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда следует, что при

 

каждом постоянном

0

функция

F i являетсмя

преобразованием

Фурье для

 

 

 

f t e t . Поэтому, в силу (7.2):

 

 

 

 

функции

 

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f t e t

 

F i ei t d .

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда находим, что

f t

1

2

 

 

 

F i e i t d .

(7.3)

 

 

 

 

Преобразуем интеграл в правой части (7.3) в интеграл по

комплексного

переменного

p i . Т.к.

const , то

 

1

p ; d

1

dp . Если меняется от –∞ до ∞, то p меня-

 

 

 

 

 

 

i

 

 

i

 

 

 

 

ется по прямой const

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

 

 

 

 

 

 

 

 

1

i

 

 

 

 

 

 

 

 

f t

 

 

F p e pt dp .

(7.4)

 

 

 

 

 

 

 

2 i

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

Функция

f t , определяемая по формуле (7.4), называ-

ется обратным преобразованием Лапласа для функции F p .

Формула (7.4) называется формулой обращения Римана– Меллина.

Вышеизложенные рассуждения являются нестрогим доказательством следующей теоремы.

Теорема. Если функция F p комплексного переменно-

го p i

аналитична в

некоторой

полуплоскости

Re p 0

и удовлетворяет

условию

lim F p 0

, то

 

 

 

p

 

83

функция

F p является преобразованием Лапласа для функ-

ции f t , причем

 

 

1

i

 

f t

 

 

F p e pt dp при t 0 и

f t 0 при t 0 .

2 i

 

 

 

i

 

Непосредсвтенное вычисление интеграла (7.4) часто сопряжено с большими трудностями. Рассмотрим случаи, когда обратное преобразование Лапласа можно найти более простыми способами.

7.2. Применение вычетов для отыскания обратного преобразования Лапласа

Теорема. Если функция F p аналитична в плоскости комплексного переменного p, кроме конечного числа особых

точек p1 , p2

,..., pn , и если lim F p 0 , то функция F p вля-

 

p

 

 

 

ется преобразованием Лапласа для функции

f t , причем

 

n

 

 

 

 

 

 

 

 

 

f t Re s F p e

pt

, pk

при t 0 и

f t 0 при t 0 .

 

k 1

(7.5)

Доказательство. Построим на плоскости комплексного переменного p замкнутый контур C, состоящий из отрезка AB прямой Re p 0 и дуги CR радиусом R с центром в начале

координат.

Отрезок AB и радиус R дуги CR выберем так, чтобы все особые точки p1 , p2 ,..., pn функции F p лежали внутри об-

ласти, ограниченной контуром C. Тогда, в силу основной теоремы о вычетах:

 

 

n

 

 

 

 

 

 

 

 

F p e

pt

dp 2 i Re s F p e

pt

, pk

, или

 

 

C

 

k 1

 

 

 

84

F p e pt dp

AB

F p e

 

n

 

 

 

 

 

 

pt

dp 2 i Re s F p e

pt

, pk .

 

 

 

CR

 

k 1

 

 

 

 

 

 

В силу леммы Жордана при t 0 : lim F p e pt dp 0 .

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

Т.к. lim

 

F p ept dp

 

 

F p e pt dp , то

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

AB

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

i

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

f t

 

 

 

 

 

F p e pt dp

 

lim

 

F p e pt dp

 

F p e pt dp

 

2 i

 

 

 

 

 

 

 

 

 

 

 

 

2 i R

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

AB

 

 

CR

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re s F p e

pt

, pk .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f t 0.

 

 

 

Также можно доказать, что при t 0

 

 

 

Пример.

F p

 

 

 

p

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

pe pt

 

 

 

 

 

 

pe pt

 

f t Res

 

 

 

 

 

 

 

,1

Res

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

p2 1

 

p2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pe

pt

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

 

 

 

t sh t.

 

 

 

p 1 2

 

2

 

 

 

p 1

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

pe

pt

 

 

 

, 1

lim

 

 

 

 

 

 

 

p 1 2

 

 

p 1

 

 

 

 

 

 

 

 

 

p

 

7.3. Оригиналы рациональных изображений

Рассмотрим случай, когда изображение F p является правильной дробно-рациональной функцией переменнного p:

F p

A p

, где

A p и

B p – многочлены, причем сте-

B p

 

 

 

 

85

пень числителя меньше степени знаменателя. Функцию F p

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

1.

 

 

M

 

, 2.

 

 

 

M

 

 

 

,

3.

 

 

Mp N

,

4.

 

 

 

 

 

Mp N

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p a

p a k

 

p2 ap b

 

 

p2 ap b 2

где

 

a2 4b 0 .

Оригиналы этих функций находятся с помо-

щью свойств преобразовнаия Лапласа:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.

 

 

M

 

=Meat ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

Mt k 1eat

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

 

 

 

 

=

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p a k

k 1 !

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4b a2

 

 

 

 

 

 

Mp N

 

 

 

 

M

p

 

 

 

 

 

 

2N

Ma

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 ap b

 

 

 

 

 

2

 

 

4b a2

 

 

 

 

 

 

 

 

 

2

 

4b a2

 

 

 

 

a

 

 

 

4b a2

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

4

 

 

 

 

 

 

 

 

2

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=e

 

 

 

 

 

2

N Ma

 

 

 

 

 

 

 

 

 

 

 

 

4b a

2

 

 

 

 

 

 

 

2

 

M cos t

sin t , где

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4b a2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

4.

Оригинал дроби вида

 

 

 

 

 

 

Mp N

 

проще всего находится

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 ap b 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с помощью вычетов.

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

 

 

Пример. F p

 

 

 

 

 

1

 

 

 

 

 

A

 

B

 

Cp D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

p

 

p2

4

 

p

p 1

 

p2 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

1

 

1

 

1

 

 

1

 

 

p

 

 

1 1

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20 p2 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 p 5 p 1

 

 

5 p2 4

 

 

 

 

 

 

 

 

 

86

ft 14 15 et 201 cos 2t 101 sin 2t .

7.4.Свертка функций

Пусть

функции

f t и

g t являются

оригиналами.

 

Функция

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h t t

f g t d

 

 

 

(7.6)

 

 

 

 

0

 

 

 

 

 

 

 

 

 

называется сверткой оригиналов и обозначается

 

 

 

 

 

 

h t f t g t .

 

 

 

 

 

Операция свертывания оригиналов обладает свойствами

 

коммутативности и ассоциативности:

 

 

 

 

 

1. g t f t f t g t ,

 

 

 

 

 

 

 

 

 

2. h t f t g t h t f

t g t .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проверим свойство коммутативности.

 

 

 

 

 

 

t

 

 

 

 

1 t

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

g t f t

 

g f

t d

 

 

 

g t f d

 

 

 

 

 

d d 1

 

 

1

1

1

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t f 1 g t 1 d 1 f t g t .

 

0

 

 

 

Свертка оригиналов также является оригиналом.

 

Теорема

об

умножении

изображений.

Если

f t =F p ,

g t =G p , то

 

 

 

 

 

 

f t g t =F p G p .

(7.7)

 

 

 

 

При свертывании оригиналов их изображения перемножаются. Доказательство. Пусть

h t f t g t t f g t d .

0

87

t

Тогда H p e pt dt f g t d .

0 0

Если в повторном интеграле изменить порядок интегрирования, то получим:

H p

0

 

 

t t1

 

 

 

 

f d g t e pt dt

 

 

 

dt dt1

0

 

 

 

 

 

 

 

 

 

 

 

f d g t1 e p t1 dt1 f e p d g t1 e pt1 dt1

 

0

0

0

0

 

 

F p G p , что и требовалось доказать.

 

 

Теорема об умножении оригиналов. Если f t

и g t

оригиналы, порядок роста которых соответственно равен 1

и

2, и если

f t =F p ,

g t =G p , то

 

 

 

 

 

 

 

 

 

1

a i

 

 

 

f t g t =

 

F q G p q dq ,

(7.8)

 

2 i

 

 

 

a i

 

 

где a σ1 , Re p a σ2 .

Этот интеграл называется сверткой изображений. Таким образом, умножение оригиналов соответствует свертыванию изображений и деление результата на 2 i .

Теорема (7.8) доказывается аналогично теореме (7.7).

7.5. Интеграл Дюамеля

Пользуясь теоремой об умножении изображений, найжем

оригинал

для изображения

вида

p F p G p , где

f t =F p , g t =G p . Для этого заметим, что

 

 

p F p G p p F p f 0 G

p f 0 G p .

 

 

 

 

Отсюда следует, что

88

p F p G p =f t g t f 0 g t

t

 

 

 

 

 

 

f

 

g t d f 0 g t .

(7.9)

 

0

 

 

 

 

 

 

Если воспользоваться коммутативностью свертки, то со-

отношение (7.9) принимает вид:

 

 

 

 

t

 

t d f 0 g t .

(7.10)

p F p G p =

g f

 

 

 

 

0

 

 

 

Интеграл в правых частях соотношений (7.9) и (7.10) называется интегралом Дюамеля. В частности, если f 0 0, то

 

 

t

 

 

 

t g t f

 

g t d

p F p G p = f

 

 

 

 

0

 

 

t g f t d .

0

Выражение p F p G p можно переписать другим спо-

собом:

 

 

 

 

 

p F p G p pG p g 0 F p g 0 F p .

 

 

 

 

 

 

Отсюда следует, что

 

 

 

 

 

p F p G p =g t f t g 0 f t =

 

t

 

 

 

 

 

 

f t d g 0 f t .

(7.11)

g

0

 

 

 

 

 

Если воспользоваться коммутативностью свертки, то со-

отношение (7.11) принимает вид:

 

 

t

 

t d g 0 f t .

(7.12)

p F p G p =

f g

 

 

0

 

 

 

Этот интеграл также нразывается интегралом Дюамеля. В частности, если g 0 0 , то

89

p F p G p =g t f t t g f t d

 

0

 

t

 

 

 

t d

.

f g

 

0

 

 

90

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]