Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лекции по математике для студентов энергетических специальностей БНТУ (IV семестр).pdf
Скачиваний:
0
Добавлен:
24.11.2025
Размер:
987.88 Кб
Скачать
ориги-

5

ГЛАВА XIV. ПРЕОБРАЗОВАНИЕ ЛАПЛАСА И ЕГО ПРИМЕНЕНИЕ (ОПЕРАЦИОННОЕ ИСЧИСЛЕНИЕ)

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

Ÿ1. Определение преобразования Лапласа и его основные свойства. Таблица изображений. Теорема Бореля

Формально интегралом Лапласа комплексной функции f(t), определенной на действитель-

ной оси называется функция

+

F (p) = e−ptf(t)dt

(1)

0

комплексного аргумента p.

Естественно интеграл (1) сходится не для любой функции f(t). Для большинства приложе-

ний достаточно рассматривать интеграл Лапласа на множестве функций, называемых

налами.

Определение 1. Оригиналом называется функция f(t) (вообще говоря, комплексная ), определенная на всей действительной оси и обладающая свойствами :

1)f(t) = 0, t ≤ 0;

2)на любом отрезке действительной полуоси [0, +) данная функция кусочно-непрерывна ;

3)существуют действительные числа a è A > 0 такие, что при всех действительных t выполняется неравенство |f(t)| < Aeat.

Число a называется показателем роста, а число α = inf a точным показателем роста

оригинала f(t). Множество всех оригиналов с точным показателем роста α мы обозначим через O(α).

Свойства 2) и 3), как мы увидим ниже, обеспечивают абсолютную сходимость интеграла Лапласа в некоторой полуплоскости, а свойство 1) связано с обратимостью преобразования Лапласа.

Приведем примеры оригиналов, которые мы будем использовать в дальнейшем.

Примеры. 1) Функция Хевисайда

{

1(t) =

0, t ≤ 0;

1, t > 0.

Очевидно, эта функция является оригиналом и точный показатель роста для нее равен

α = 0.

Всюду в дальнейшем мы для удобства будем предполагать, что функция Хевисайда незримо присутствует в записи функции f(t) в качестве множителя, что обеспечивает выпол-

нение условия 1) из определения оригинала. 2) Степенная функция tn, n N.

Для нее по правилу Лопиталя (глава V, Ÿ4) для любого ε > 0

lim

tn

=

lim

(tn)

=

lim

ntn−1

=

lim

n(n − 1)tn−2

= . . . =

lim

n!

= 0,

 

 

εeεt

ε2eεt

 

t→+eεt

 

t→+(eεt)

 

t→+

 

t→+

 

t→+εneεt

 

поэтому найдется такое T > 0, ÷òî ïðè âñåõ t > T выполняется неравенство

tn < 1 tn < eεt.

eεt

Кроме того, tn ≤ T neεt на отрезке [0, T ]. Тогда, полагая A = max{1, T n}, мы получаем tn < Aeεt, t ≥ 0

6

и, значит, произвольное ε > 0 является показателем роста данной функции. Поэтому степенная функция является оригиналом и поскольку, очевидно, отрицательное число не может быть

показателем роста

этой функции, то

t

n

O(0).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λt

, λ C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3) Экспонента e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оригинал

с точным

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и, стало быть, данная функция

показателем роста

α =

Re

λ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eλt

 

=

 

etRe λ+itIm λ

 

= etRe λ

 

 

 

 

 

 

 

 

4) Тригонометрические функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Òàê êàê (глава XIII, Ÿ3)

 

 

 

 

sin λt, cos λt, λ C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin λt =

eiλt − e−iλt

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

òî

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

)

 

 

 

(

 

 

 

 

) |

 

 

 

 

 

 

|

 

 

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

 

 

 

 

 

 

 

 

2i

at

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

ßñíî, ÷òî

 

 

 

=

 

eiλt

 

 

 

 

1

 

 

ïðè

+

 

 

 

 

Im

=

поэтому

 

 

tIm λ

 

 

 

Im

 

 

 

Совершенно

 

sin λt

 

 

− e−iλt

 

 

 

eiλt

 

 

e iλt

 

1 etIm λ + e

 

 

 

 

e

Im λ t, t

0.

 

 

t

 

+

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

sin λt

|

 

 

= +

 

 

 

 

 

 

 

a <

|

 

 

λ

,

 

 

 

sin λt

 

O(

|

 

 

λ

).

 

 

 

 

 

→ ∞ |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

 

 

 

 

 

 

 

 

 

 

 

 

|

 

 

 

 

аналогично мы можем убедиться в том же и для функции cos λt.

5) Гиперболические функции sh λt, ch λt, λ C.

Известно (глава XIII, Ÿ3), что данные гиперболические функции связаны с соответствующими тригонометрическими формулами

sh λt = −i sin iλt, ch λt = cos iλt,

поэтому они также являются оригиналами с точным показателем роста α = |Im()| = |Re λ|. Докажем теперь теорему о сходимости интеграла Лапласа (1).

Теорема 1. Интеграл Лапласа (1) для оригинала f(t) O(α) сходится абсолютно в ком-

плексной полуплоскости Re p > α. Более того, в этой полуплоскости

 

lim F (p) = 0,

(2)

Re p→+

 

интеграл Лапласа является аналитической функцией и допускает дифференцирование по аргументу p под знаком интеграла.

Действительно, выбрав при Re p = x > α показатель a роста оригинала так, чтобы x > a >

α, мы получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и, поскольку несобственный интеграл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e−ptf(t) < Ae(x−a)t, t R

 

 

 

 

 

 

+

(x

 

a)t

 

1

 

(x

 

a)t

 

+

 

1

 

 

 

(x

 

1

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

dt = x a e

 

 

 

0

 

 

 

 

lim

e

a)t +

 

= 0 +

 

= x a

 

 

 

 

 

 

 

= x a t→+

 

x a

x a

сходится, то по признаку сравнения

(глава VII, Ÿ4, пункт 1) интеграл Лапласа (1) сходится

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

абсолютно в полуплоскости Re p = x > α.

Воспользовавшись свойствами определенного интеграла (глава VII, Ÿ1) и предыдущим равенством, мы найдем оценку для интеграла Лапласа:

|

|

 

 

 

 

 

x a

 

F (p)

=

 

+e−ptf(t)dt

 

+e−ptf(t) dt

 

+Ae(x−a)tdt =

A

.

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда, очевидно, следует (2).

Доказательство аналитичности и дифференцируемости интеграла Лапласа ввиду его сложности мы здесь приводить не будем.

Из доказанной теоремы следует, что интеграл Лапласа, как функция комплексного аргумента p, обладает хорошими аналитическими свойствами .

Определение 2. Отображение, по которому каждому оригиналу f(t) формулой (1) ставится в соответствие интеграл Лапласа F (p), называется преобразованием Лапласа. Функцию F (p) мы будем называть изображением оригинала f(t).

7

Обозначать преобразование Лапласа мы будем через f(t) → F (p).

Пример 1. Найти изображения функции Хевисайда и экспоненты.

Решение. Как мы показали в предыдущем примере, эти функции являются оригиналами. Для функции Хевисайда при Re p > 0

+

 

pt 1dt =

 

1

 

 

pt

 

+

 

 

 

 

 

1

 

 

 

 

pt +

1

 

 

 

 

 

1

 

 

1

 

 

e

 

 

e

 

 

 

 

 

 

 

 

lim e

 

=

0 +

=

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

 

 

p

 

 

 

 

 

 

= p t→+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

p

 

 

 

 

 

p p

и, таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1(t)

 

 

, Re p > 0.

 

 

 

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

Для экспоненты eλt, λ C ïðè Re p > Re λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

e−pteλtdt =

e(p−λ)tdt =

1

 

e(p−λ)t

0

=

 

p

λ

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

=

 

 

 

 

lim e(p−λ)t

+

 

 

= 0 +

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

и, значит,

 

 

p − λ t→+

 

 

 

 

 

 

 

 

 

 

− λ

 

 

p − λ p − λ

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eλt

 

, Re p > Re λ.

 

 

 

 

(4)

 

 

 

 

 

 

 

 

p

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Займемся теперь свойствами преобразования Лапласа , связанными с различного рода операциями над оригиналами и изображениями.

1. Линейность преобразования Лапласа.

Åñëè f1(t) O(α1), f2(t) O(α2), è f1(t) → F1(p), f2(t) → F2(p), то для любых постоянных c1, c2 C

c1f1(t) + c2f2(t) → c1F1(p) + c2F2(p), Re p > max1, α2}.

(5)

Доказательство этого свойства немедленно следует из формулы (1).

Во всех следующих свойствах мы будем считать, что f(t) O(α) è f(t) → F (p), Re p > α.

2. Смещение изображения.

Для любого комплексного числа λ

eλtf(t) → F (p − λ), Re p > α + Re λ.

(6)

В самом деле,

 

+

 

 

+

 

 

eλtf(t)

e−pt

(

eλtf(t) dt =

e(p−λ)tf(t)dt = F (p − λ).

 

 

0

 

)

0

 

 

 

3. Запаздывание оригинала.

 

 

 

 

 

 

При любом фиксированном t0 > 0

 

 

 

 

 

 

f(t − t0) → e−t0pF (p), Re p > α.

(7)

Для доказательства заметим, что

 

 

+

 

 

 

+

 

 

 

F1(p) =

e−ptf(t − t0)dt = e−ptf(t − t0)dt

 

 

 

0

 

 

 

 

t0

 

и проведем в последнем интеграле подстановку z = t − t0. В результате мы получим

 

+

 

 

 

 

+

 

F1(p) = 0

e−p(z+t0)f(z)dz = e−t0p

0

e−pzf(z)dz = e−t0pF (p),

 

в чем и требовалось убедиться.

8

4. Дифференцирование оригинала.

Если существует производная f(t) è f(t) O(α), òî

f(t)

pF (p)

f(+0), Re p > α.

(8)

 

 

 

 

Действительно, интегрируя по частям в интеграле Лапласа

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F1(p) = 0

 

e−ptf(t)dt

 

 

 

 

ïðè Re p > α, мы получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

0

 

 

 

+

 

 

 

 

 

 

 

 

(p) =

e

pt

f(t) 0

f(t)de

pt

 

 

lim e

ptf(t)

lim e

ptf(t)+

F1

 

 

 

 

 

= t→+

 

t→+0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 0

f(t)e−ptpdt = 0 − f(+0) + p 0

 

e−ptf(t)dt = pF (p) − f(+0).

5. Интегрирование оригинала.

 

 

 

 

 

 

 

 

 

 

 

Интеграл 0 t

f(s)ds O(β), β = max{α, 0} è

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

(p)

 

 

 

 

 

 

 

 

 

 

 

 

f(s)ds →

F

 

, Re p > β.

 

 

(9)

 

 

 

 

 

 

 

p

 

 

 

t

Покажем сначала, что интеграл f1(t) = f(s)ds является оригиналом. Для этого восполь-

0

зуемся свойствами определенного интеграла (глава VII, Ÿ1) и определением оригинала:

 

|f1(t)| ≤ 0 t |f(s)|ds ≤ A 0 t easds.

β = max α, 0 .

точным показателем роста интеграла

(

 

)

 

 

 

 

 

 

Отсюда при a > 0 |f1(t)| ≤ Aa eas

0t = Aa

eat 1

 

< Aa eat, à ïðè a ≤ 0 |f1

(t)| ≤ At и, стало быть,

Пусть f1(t) → F1(p). Тогда по

предыдущему свойству

 

 

 

 

{ }

 

 

 

 

 

 

оригинала является число

 

 

 

 

f(t)

pF

(p)

f (+0) = pF

(p)

0 = pF

(p).

 

 

1

 

1

 

1

 

 

 

 

1

 

 

1

 

 

 

С другой стороны, f(t) = f(t) и, значит, f

(t)

F (p). Поэтому, pF1(p) = F (p), откуда,

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (p)

 

 

 

 

 

 

 

 

 

 

 

 

F1(p) =

 

 

 

 

, Re p > β.

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

6. Дифференцирование изображения.

 

 

 

 

 

 

 

 

 

 

 

Функция tf(t) O(α) è

 

 

tf(t) → −F (p), Re p > α.

 

 

 

(10)

 

 

 

 

 

 

Поскольку t O(0), f(t) O(α), то при любом ε > 0 найдется Aε > 0 такое, что

|tf(t)| < Aεe(α+ε)t, t ≥ 0

9

и, таким образом, tf(t) O(α). Ïî теореме 1 интеграл Лапласа мы можем дифференцировать

под знаком интеграла и, значит,

 

 

 

=

+

e−ptf(t)

pdt =

+e−pt( tf(t))dt,

 

 

F (p) =

+e−ptf(t)dt

 

 

 

 

 

 

 

 

ò. å.

 

 

0

 

p

 

0 (

 

)

0

 

 

 

 

 

tf(t) → −F (p), Re p > α.

 

 

7. Интегрирование изображения.

 

 

 

 

 

Åñëè

f(t)

O(α), òî

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t)

p

 

 

 

 

 

 

 

 

 

f

 

F (q)dq, Re p > α,

 

(11)

 

 

 

 

t

 

 

причем интегрирование в правой части проводится по любой непрерывной кривой, ведущей из точки p комплексной плоскости в бесконечно удаленную точку, вдоль которой Re q → +∞.

Для доказательства обозначим изображение оригинала

f(t)

через F1(p)

и воспользуемся

t

 

предыдущим свойством. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t ·

 

→ −F1(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

(p), ò. å. F

 

t

 

 

 

 

 

 

 

 

 

 

 

 

и, стало быть, f(t)

→ −

(p) =

F (p). Отсюда, после интегрирования, учитывая

(2), мы получим:

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

F1(q)dq = F1

 

 

 

= Re q→+

1

 

 

1

 

1

 

 

 

1

 

(q) p

(q)

(p) = 0

(p) =

(p) =

 

 

 

 

 

 

lim

F

 

 

F

 

F

 

F

 

F (q)dq.

Значит,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F1(p) = p

F (q)dq, Re p > α.

 

 

 

 

 

 

Найдем теперь, пользуясь доказанными свойствами, изображения оригиналов, рассмотренных в примерах 2), 4), 5).

Применив к функции Хевисайда свойство дифференцирования изображения, мы благодаря

(3) найдем

 

 

 

 

 

 

1

 

 

t = t · 1(t) → − (

 

)

=

1

,

ò. å.

p

p2

1!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t →

 

.

 

 

 

 

p2

 

 

 

Аналогично,

t2

t3

·

tn

 

 

 

 

 

1!

 

 

 

 

 

 

 

 

 

 

 

 

= t · t → − (

 

 

)

 

= 1! · (2) · p3

p2

 

 

 

 

 

 

2!

 

 

 

3!

 

 

 

 

 

 

= t · t2

→ − (

 

)

 

=

 

 

,

 

 

 

 

p3

 

p4

 

 

 

 

·

 

 

 

· ·

 

 

 

 

 

·

 

 

·

 

 

·

 

= t

·

tn

1

→ − (

(n − 1)!

)

=

n!

.

 

 

 

 

 

 

 

 

 

 

pn

 

 

pn+1

Таким образом, при любом натуральном n

2! = p3 ,

··

tn

n!

, Re p > 0.

(12)

pn+1

10

Для функции sin λt, λ C воспользуемся ее представлением через экспоненту:

 

 

 

 

sin λt =

eiλt − e−iλt

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2i

 

 

 

 

 

 

 

 

 

 

 

Применив свойство линейности преобразования Лапласа и формулу (4), получим:

 

 

 

1

(

 

 

 

 

1

 

 

 

 

 

1

 

)

 

 

 

λ

 

 

 

 

sin λt →

 

 

 

 

 

 

 

 

 

=

 

 

 

.

 

 

 

 

2i

p − λi

p + λi

p2 + λ2

 

 

 

 

Таким образом, для любого комплексного λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin λt

 

 

 

 

 

λ

, Re p > |Im λ|.

 

 

(13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 + λ2

 

 

Аналогично,

 

eiλt + e

iλt

 

 

(

 

 

 

 

 

 

 

 

 

 

 

 

) =

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

1

 

 

 

 

p

 

 

cos λt =

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

,

 

2

 

 

 

 

 

2

p − λi

p + λi

p2 + λ2

ò. å.

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos λt

 

 

 

, Re p > |Im λ|.

 

 

(14)

 

 

 

 

 

 

 

 

 

p2 + λ2

 

 

Для гиперболических функций

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh λt = −i sin iλt, ch λt = cos iλt, λ C.

 

 

 

 

Поэтому, применяя формулы (13) и (14), мы найдем:

 

 

 

 

 

 

 

 

 

 

 

sh λt → −i ·

 

 

 

 

 

 

 

 

 

 

 

λ

, ch λt →

 

 

p

 

 

p

 

 

=

 

 

 

 

=

 

.

p2 + (λi)2

 

p2 − λ2

p2 + (λi)2

p2 − λ2

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

è

sh λt →

λ

, Re p > |Re λ|

(15)

p2 − λ2

ch λt →

 

p

, Re p > |Re λ|.

(16)

p2 − λ2

 

Сведем все найденные изображения (3), (4), (12) (16) в таблицу.

Таблица изображений

 

1

 

 

 

 

 

 

 

 

 

 

 

1)

1(t)

 

 

 

, Re p > 0;

p

2)

eλt

 

1

 

, Re p > Re λ;

p

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

3)

tn

n!

, n N, Re p > 0;

pn+1

4)

sin λt →

 

 

 

λ

 

 

, Re p > |Im λ|;

 

 

 

p2 + λ2

5)

cos λt →

 

 

p

 

 

 

, Re p > |Im λ|;

p2 + λ2

6)

sh λt →

 

 

 

 

λ

 

, Re p > |Re λ|;

 

 

 

p2 − λ2

7)

ch λt →

 

 

 

 

p

 

, Re p > |Re λ|.

p2 − λ2

 

Пользуясь этой таблицей и свойствами 1 7 (соответственно, формулами (5) (11)), мы можем найти изображения и других, более сложных оригиналов.

Пример 2. Найти изображения функций :

 

 

 

t

{

 

 

a) f1

(t) = 2tt sin2 3t; b) f2

(t) =

 

e2s cos 5s

ds; c) f3(t) =

sin t,

t [2, 3];

 

0,

 

 

 

s

t / [2, 3].

0

11

Решение. a) Поскольку

sin2 3t = 12(1 cos 6t),

то по таблице изображений

 

1

1

 

p

sin2 3t →

 

(

 

 

).

2

p

p2 + 36

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

 

 

 

 

 

 

 

 

 

1

1

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

t sin2

3t → −

 

(

 

 

)

=

 

 

 

 

 

 

 

 

2

p

p2 + 36

 

 

=

 

1

p

 

2

1

 

(p2 + 36) − p · 2p

=

1

 

 

 

1

+

36 − p2

.

2

 

·

 

2

(p2

 

 

(

 

 

(p2 + 36)2

 

)

 

(p2 + 36)2 )

Осталось применить свойство смещения изображения:

f (t) = eln 2

 

tt sin2

3t

 

1

 

 

 

1

+

36

(p − ln 2)2

·

2

((p

 

ln 2)2

((p

 

 

1

 

 

 

 

ln 2)2

+ 36)2

 

 

 

 

 

 

 

 

 

 

 

 

)

, Re p > ln 2.

b) Ïî правилу Лопиталя (глава V, Ÿ4)

lim

e2t cos 5t

= lim

(e2t cos 5t)

= lim(

2e2t + 5 sin 5t) =

2,

t 0

t

t

0

t

t

0

 

 

 

 

 

 

 

 

 

 

поэтому функция

e2t cos 5t

t

является оригиналом с точным показателем роста α = 0. Найдем изображение этой функции, применив таблицу и свойство интегрирования изображения1:

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

e

2t

cos 5t

 

 

 

1

 

 

 

 

 

 

 

q

 

 

e2t cos 5t

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

p (

 

 

 

 

 

 

 

 

)dq =

p + 2

p2 + 25

 

 

 

 

t

 

 

 

q + 2

q2 + 25

 

 

d(q + 2)

 

 

 

1

 

 

 

d(q2 + 25)

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= p (

 

 

 

 

 

·

 

 

 

 

 

 

) = p (d ln(q + 2)

 

d ln(q2 + 25)) =

 

 

 

q + 2

2

 

 

 

q2 + 25

 

2

 

 

 

 

q + 2

 

 

 

 

 

 

 

 

 

q + 2

 

 

 

 

 

 

 

 

 

 

 

 

q + 2

 

 

 

 

 

p + 2

 

 

=

d ln

 

 

= ln

 

 

p

=

lim

 

ln

 

 

ln

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q2 + 25

 

 

 

 

 

 

 

 

 

q2 + 25

 

 

 

Re q→+

 

 

 

 

q2 + 25

 

 

 

 

 

p2 + 25

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 +

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

lim

 

 

ln

 

 

 

 

 

q

 

 

 

 

 

 

ln

p + 2

 

 

= ln 1

 

ln

 

 

p + 2

 

=

 

 

ln

 

p + 2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re q→+

 

1 + 25q2

 

 

p2 + 25

 

 

 

 

 

p2 + 25

 

 

 

p2 + 25

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

1

 

p + 2

 

f2(t) → −

 

ln

 

 

, Re p > 0.

p

 

 

 

p2 + 25

c) С помощью функции Хевисайда функция f3(t) может быть представлена в виде

f3(t) = 1(t − 2) sin t − 1(t − 3) sin t = 1(t − 2)(sin(t − 2) cos 2 + sin 2 cos(t − 2))− −1(t − 3)(sin(t − 3) cos 3 + sin 3 cos(t − 3)).

По таблице и свойству линейности изображения

sin t cos 2 + sin 2 cos t →

cos 2

 

p sin 2

 

cos 2 + p sin 2

 

+

 

=

 

.

p2 + 1

p2 + 1

p2 + 1

Аналогично,

sin t cos 3 + sin 3 cos t → cos 3 + p sin 3. p2 + 1

1Здесь используется главное значение логарифмической функции комплексного аргумента (глава XIII, Ÿ3).

12

Применив свойство запаздывания оригинала, мы найдем:

1(t − 2)(sin(t − 2) cos 2 + sin 2 cos(t − 2)) → e2p cos 2 + p sin 2, p2 + 1

1(t − 3)(sin(t − 3) cos 3 + sin 3 cos(t − 3)) → e3p cos 3 + p sin 3. p2 + 1

Окончательно,

f (t)

e2p(cos 2 + p sin 2) − e3p(cos 3 + p sin 3)

, Re p > 0.

p2 + 1

3

 

Введем теперь одну любопытную интегральную операцию над оригиналами, которая называется сверткой.

Определение 3. Сверткой двух оригиналов f1(t) è f2(t) называется функция, которая обозначается через f1(t) f2(t) и определяется при всех действительных t равенством

t

f1(t) f2(t) = f1(s)f2(t − s)ds.

0

Укажем простейшие свойства свертки . 1) f1(t) f2(t) = f2(t) f1(t).

Для проверки этого свойства выполним в интеграле в левой части данного равенства подстановку z = t − s. Учитывая при этом, что s = t − z, ds = dz и отрезок [0, t] преобразуется

в отрезок [t, 0], мы получим:

t 0 t

f1(t) f2(t) = f1(s)f2(t − s)ds = − f1(t − z)f2(z)dz = f2(z)f1(t − z)dz = f2(t) f1(t).

0 t 0

2) (c1f1(t) + c2f2(t)) f3(t) = c1(f1(t) f3(t)) + c2(f2(t) f3(t)), c1, c2 C.

Это свойство немедленно следует из линейности интеграла.

Покажем, что свертка является оригиналом. Действительно, f1(t) f2(t) = 0 ïðè t < 0, à ïðè t ≥ 0 свертка, как интеграл с переменным верхним пределом кусочно-непрерывной функции, непрерывна (глава VII, Ÿ1, свойство 7)). Далее, если f1(t) O(α1), f2(t) O(α2), то для любого ε > 0 найдутся константы A1 > 0, A2 > 0 такие, что при всех действительных t

|f1(t)| < A1e(α1+ε)t, |f2(t)| < A2e(α2+ε)t

и, значит,

|fi(t)| < Ae(α+ε)t, i = 1, 2; t R, A = max{A1, A2}, α = max1, α2}.

Найдем оценку свертки, использовав свойства 3), 4) определенного интеграла (глава VII, Ÿ1) и предыдущие неравенства:

t t t

|f1(t) f2(t)| ≤ |f1(s)||f2(t − s)|ds ≤ A2 e(α+ε)se(α+ε)(t−s)ds = A2e(α+ε)t ds = A2e(α+ε)tt.

0 0 0

Отсюда, учитывая, что t O(0), à ε > 0 можно выбирать сколь угодно малым, мы заключаем,

÷òî f1(t) f2(t) O(α).

Найдем изображение свертки.

Теорема Бореля (умножения изображений). Åñëè f1(t) O(α1), f2(t) O(α2) è f1(t) → F1(p), f2(t) → F2(p), òî

f1(t) f2(t) → F1(p)F2(p), Re p > α, α = max1, α2}.

13

Д о к а з а т е л ь с т в о. Обозначим изображение свертки через F (p). Тогда по определению изображения

F (p) =

 

f2(t))dt =

 

t

f1(s)f2(t

s)ds dt =

t

e−ptf1(s)f2(t

s)ds.

+e−pt(f1(t)

 

+e−pt

 

 

+dt

 

 

 

0

 

 

0

0

 

 

 

 

0

0

 

 

 

 

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

+

+

F (p) = ds

e−ptf1(s)f2(t − s)dt.

0s

Проведем в интеграле по переменной t замену переменной z

= t − s и поскольку при этом

t = z + s, dt = dz и полуось [s, +) отображается в полуось [0, +), òî

 

 

 

 

 

 

 

 

+

+

 

 

 

 

 

 

+

 

 

 

 

+

 

 

 

 

 

 

F (p) = 0

ds

0

e−p(z+s)f1(s)f2(z)dz = 0

e−psf1(s)ds 0

 

e−pzf2(z)dz = F1(p)F2(p),

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 3. Найти изображение интеграла f(t) = 0 t

5s sin(t + s)ds.

 

 

 

Решение. Переписав подынтегральную функцию в виде

 

 

 

 

 

 

 

 

 

 

 

 

5s sin(t + s) = 5s sin((t − s) + 2s) = 5s sin 2s cos(t − s) + 5s cos 2s sin(t − s),

мы, воспользовавшись определением свертки, замечаем, что

 

 

 

 

 

 

 

 

 

f(t) = 0 t

5s sin 2s cos(t − s)ds + 0 t

5s cos 2s sin(t − s)ds = (5t sin 2t) cos t + (5t cos 2t) sin t.

По таблице и свойству смещения изображения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5t sin 2t = eln 5·t sin 2t

 

 

 

 

2

 

 

 

, 5t cos 2t

 

 

 

p − ln 5

.

 

 

 

 

 

 

 

(p

 

ln 5)2

+ 4

(p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln 5)2

+ 4

 

 

Осталось применить теорему Бореля:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(t)

 

 

 

 

2

 

 

p

+

 

p − ln 5

 

 

1

 

=

 

 

 

 

3p − ln 5

 

 

, Re p > ln 5.

(p

 

 

ln 5)2 + 4 · p2 + 1

(p

 

· p2 + 1

((p

 

 

 

 

 

 

 

 

ln 5)2

+ 4

 

 

 

ln 5)2 + 4)(p2 + 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Докажем полезную в некоторых приложениях формулу Дюамеля. Пусть в предположениях теоремы Бореля существует также производная f2(t) è f2(t) O(α2). Воспользовавшись теоремой Бореля и свойством дифференцирования оригинала, получим:

(f1(t) f2(t))→ p(F1(p)F2(p)) (f1(t) f2(t))|t=+0 = pF1(p)F2(p) 0 = pF1(p)F2(p).

Свертка представляет собой интеграл по переменной s с переменным верхним пределом t от функции, которая зависит от s è t. Применив теорему о дифференцировании под знаком

интеграла (глава VIII, Ÿ2), мы найдем:

 

 

(f1(t) f2(t))=

t

f1(s)f2(t − s)ds

= f1(t)f2(0)+t

f1(s)f2(t−s)ds = f1(t)f2(0)+f1(t) f2(t).

Значит,

0

 

0

 

 

f1(t)f2(0) + f1(t) f2(t) → pF1(p)F2(p) формула Дюамеля.

 

 

(17)

Найдем в завершение этого параграфа изображение периодического оригинала.

Теорема 1.
В этом параграфе мы рассмотрим

14

 

 

 

 

 

 

 

 

 

 

Теорема 2. Изображение F (p) T -периодического оригинала f(t) находится по формуле

 

 

T

 

 

 

 

 

 

 

 

 

 

F (p) = 0

e−ptf(t)dt

, Re p > 0.

 

 

 

(18)

 

 

1 − e−pT

 

 

 

 

 

 

Для доказательства проведем в интеграле

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

F (p) = 0

e−ptf(t)dt

 

 

 

 

 

подстановку s = t + T. Тогда t = s − T, dt = ds и полуось [0, +] преобразуется в полуось

[T, +]. Следовательно,

 

T )ds = epT

+e−psf(s)ds = epT

F (p)

 

T e−ptf(t)dt ,

F (p) = +e−p(s−T )f(s

 

 

 

 

 

 

 

 

T

 

 

T

 

 

 

 

0

откуда и следует формула (18).

 

 

 

 

 

 

 

 

 

Пример 4. Найти изображение функции f(t) = {t}, ãäå {t} дробная часть числа t, ò.å.

разность между числом и его целой частью.

 

 

 

 

 

 

 

Решение. Ясно, что эта функция является 1-периодической и на промежутке [0, 1) {t} = t.

Y

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

2

3

 

4

5

 

6

 

X

O

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

T e−ptf(t)dt = 1 e−pttdt.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

Применяя метод интегрирования по частям, получим:

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

tde−pt

 

 

 

 

 

 

 

 

 

 

 

1

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e−pttdt = p

= −p e−ptt 0

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

1

 

 

 

1

 

 

 

e

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

=

p

(e−p +

p

e−pt

0)

=

p

(e−p +

p

p

) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e−ptdt =

1 − e−p(p + 1). p2

Воспользовавшись формулой (18), находим:

f(t)

1 − e−p(p + 1)

.

 

p2(1 − e−p)

Ÿ2. Нахождение оригинала по изображению. Обратное преобразование Лапласа

методы, которые позволяют по известному изображению F (p) найти соответствующий оригинал , т. е. функцию f(t), для которой f(t) → F (p).

Обратное соответствие мы будем обозначать точно также, как и прямое, т. е. F (p) → f(t).

Сформулируем сначала важную для приложений теорему, гарантирующую единственность оригинала, найденного по изображению.

Åñëè f1(t) è f2(t) два оригинала, имеющие одно и тоже изображение, то f1(t) = f2(t) во всех точках непрерывности этих функций.

15

Доказательство этого утверждения можно найти, например, в учебнике Свешникова А.Г. è Тихонова А.Н. по теории функций комплексной переменной, приведенному в списке литературы.

Простейшим из методов нахождения оригинала по изображению является табличный, который основан на таблице изображений, перечисленных в предыдущем параграфе свойствах преобразования Лапласа и теореме Бореля. Для удобства перепишем справа налево приведенную в предыдущем параграфе таблицу изображений.

Таблица оригиналов

1)p1 1(t), Re p > 0;

2)

 

1

 

→ eλt, Re p > Re λ;

p

λ

 

 

 

 

1tn−1

3)pn (n − 1)!, n N, Re p > 0; 1 1

4)p2 + λ2 λ sin λt, Re p > |Im λ|;

5)

p

cos λt, Re p > |Im λ|;

 

p2 + λ2

 

1

1

 

6)

 

 

sh λt, Re p > |Re λ|;

p2 − λ2

λ

7)

p

ch λt, Re p > |Re λ|.

p2 − λ2

Табличным методом мы можем эффективно найти оригинал для изображения, являющегося

рациональной функцией. В самом деле, такая функция, как следует из теоремы 1 предыдущего

параграфа, непременно является правильной и поэтому (глава VI, Ÿ3) ее можно разложить в

сумму простейших дробей вида

a

 

 

 

 

(p − λ)r ,

(1)

ãäå a, λ комплексные числа, 1 ≤ r N. Оригинал же для дроби (1) легко находится по таблице оригиналов и свойству смещения изображения. В самом деле, при r = 1

 

a

 

→ aeλt.

p

λ

 

 

 

Åñëè æå r > 1, òî

a

→ ae

λt tr−1

 

 

 

.

(p − λ)r

(r − 1)!

Если все коэффициенты рациональной функции действительные числа, то мы можем использовать разложение этой функции на сумму простейших дробей вида (1) с действительны-

ìè a, λ и дробей вида

bp + c

, b, c, k, l R, k2 4l < 0, 1 ≤ s N.

(2)

(p2 + kp + l)s

Найдем оригинал для дроби (2). При s = 1 после выделения в знаменателе полного квадрата, мы получим

bp + c

= b

·

p − µ

+ c

1 ·

ν

 

,

p2 + kp + l

(p − µ)2 + ν2

(p − µ)2

+ ν2

 

 

 

ãäå µ, ν > 0, c1 некоторые действительные числа. Тогда по таблице оригиналов и свойству смещения изображения

 

bp + c

→ eµt(b cos νt + c1 sin νt).

 

 

 

 

p2 + kp + l

 

 

 

Åñëè æå s > 1, то найдя оригиналы f1(t) è f2(t) для дробей

bp + c

è

1

, ñîîò-

 

 

p2 + kp + l

p2 + kp + l

ветственно, достаточно по теореме Бореля (s − 1) раз последовательно использовать свертку

16

оригиналов, а именно,

 

 

 

bp + c

 

 

=

 

bp + c

 

 

·

 

 

 

 

1

 

 

 

 

→ f3(t) = f1(t) f2(t),

 

 

(p2 + kp + l)2

p2 + kp + l

p2 + kp + l

 

 

 

bp + c

 

=

 

 

 

bp + c

 

 

 

 

·

 

 

1

 

 

 

 

→ f4

(t) = f3(t) f2(t),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(p2 + kp + l)3

(p2 + kp + l)2

p2 + kp + l

· ·

 

·

·

·

 

·

 

 

 

 

·

 

 

 

 

· ·

 

· · ·

 

 

 

bp + c

 

 

 

 

 

=

 

bp + c

 

 

 

·

 

 

 

1

 

→ fs(t) f2(t).

 

 

 

(p2 + kp + l)s

(p2 + kp + l)s−1

p2 + kp + l

 

Пример 1. Найти оригинал, соответствующий изображению

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3p2 2p + 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p3 + 1

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Поскольку

 

 

 

 

3p2 2p + 1

 

 

 

 

 

 

3p2 2p + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p3 + 1

 

 

 

 

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

то разложение данной дроби на простейшие имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3p2 2p + 1

=

 

 

a

 

+

 

 

bp + c

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p3 + 1

 

 

 

 

 

p + 1 p2 − p + 1

Найдем неизвестные коэффициенты разложения. Так как

 

 

 

 

 

 

 

a

 

+

 

 

 

 

bp + c

 

=

a(p2 − p + 1) + (bp + c)(p + 1)

,

 

 

 

 

p + 1

 

 

 

p2 − p + 1

 

 

 

 

 

 

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

òî

3p2 2p + 1 = a(p2 − p + 1) + (bp + c)(p + 1)

ïðè âñåõ p C. Отсюда при p = 1 мы находим 6 = 3a и, значит, a = 2. Далее при p = 0 1 = a + c, откуда c = 1. Наконец, при p = 1 мы получаем 2 = a + 2b + 2c и, стало быть, b = −c = 1. Таким образом,

3p2 2p + 1

=

2

+

p − 1

.

p3 + 1

 

 

 

p + 1 p2 − p + 1

Найдем оригиналы для каждой из простейших дробей. Для первой по таблице оригиналов

2

2e−t.

p + 1

Вторую предварительно преобразуем, выделив в знаменателе полный квадрат:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p − 1

 

 

p − 1

 

 

 

p −

1

 

1

 

 

 

 

 

p −

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

3

 

 

 

 

 

 

 

 

=

 

 

=

 

2

 

2

=

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

p + 1

 

p

21

)

+ 43

 

(p

21 )

+

43

 

 

 

 

 

 

+ (

 

 

)

2 3

·

 

 

 

 

 

(

 

)

2

 

 

 

 

 

 

p −

1

3

p −

1

+

3

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

(

 

 

 

 

 

(

 

 

)

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

Отсюда по таблице оригиналов и свойству

 

 

(

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

смещения изображения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p − 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

3

 

1

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e2

 

cos

 

 

t

 

 

sin

 

 

 

t .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 − p + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно, получаем:

 

 

(

 

2

 

3

 

 

 

2

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2p + 1

 

 

 

 

 

 

t

 

 

 

 

 

 

 

1

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2e−t + e2 (cos

3

t

 

sin

 

 

 

 

t).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p3 + 1

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим теперь общий метод нахождения оригинала по изображению.

Теорема 2. Предположим, что функция F (p) комплексного аргумента p = x + yi удовле-

творяет следующим условиям :

a) F (p) аналитическая функция в полуплоскости Re p > α и здесь lim F (p) = 0;

p→∞

и во всей комплексной

17

b) при любом фиксированном Re p = x > α сходится несобственный интеграл

x+i∞

|F (p)|dy,

x−i∞

где интегрирование проводится по прямой оси. Тогда функция F (p) в полуплоскости определяемой равенством

p = x + yi, y (−∞, +), параллельной мнимой Re p > α является изображением функции f(t),

 

 

x+i∞

 

f(t) =

1

eptF (p)dp, Re p > α

(3)

2πi

x−i∞

èэта функция является оригиналом.

Ñдоказательством этой теоремы также можно ознакомиться в учебнике Свешникова А.Г. è Тихонова А.Н. по теории функций комплексной переменной, имеющемся в списке литературы.

Соотношение (3) называется формулой Меллина èëè формулой обратного преобразования Лапласа, которая в общем виде решает задачу нахождения оригинала по известному изобра-

жению. Интеграл в правой части формулы (3) мы будем называть интегралом Меллина.

При определенных условиях мы можем вычислить интеграл Меллина, использовав теорию вычетов (глава XIII, Ÿ7).

Теорема 3. Пусть функция F (p) комплексного аргумента p удовлетворяет условиям теоремы 2 и, сверх того, она допускает аналитическое продолжение в полуплоскость Re p ≤ α,

где она имеет конечное число изолированных особых точек p1, p2, . . . , pn плоскости удовлетворяет условию

lim F (p) = 0.

p→∞

Тогда эта функция является изображением оригинала

n

f(t) =

res(eptF (p)).

(4)

k=1 pk

Для доказательства заметим, что по теореме 2 оригинал для данного изображения находится по формуле (3). Применяя к интегралу Меллина следствие из теоремы о вычетах (глава XIII, Ÿ7), мы получим

x+i∞

 

n

x i

 

 

 

pt

res(eptF (p)),

e

 

F (p)dp = 2πi k=1 pk

− ∞

 

 

 

откуда и следует формула (4).

В частности, формула (4) позволяет находить оригинал для изображения, являющегося рациональной функцией, не разлагая ее на простейшие дроби. Как известно (глава XIII, Ÿ6),

для такой функции все особые точки, которые мы обозначим через p1, p2, . . . , pn, являются

полюсами. Порядки этих полюсов мы обозначим через r1, r2, . . . , rn, соответственно. Эти по-

люсы совпадают с корнями соответствующей кратности знаменателя рациональной функции. Осталось воспользоваться полученной в главе XIII, Ÿ7 формулой для вычисления вычетов в

полюсах:

 

 

n

 

 

 

 

 

 

 

 

drk1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

(p

p )rk eptF (p) .

 

 

 

 

 

 

 

 

(rk

 

 

 

 

dprk1

 

 

 

 

 

 

 

f(t) =

k=1

1)! p→pk

(

 

 

k

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Åñëè все корни p1, p2, . . . , pn знаменателя рациональной функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (p) =

F1(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F2(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

простые, то вычеты в них равны

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pt

 

 

 

 

 

pt F1(p)

 

 

pkt

·

 

 

 

F1

(pk)

 

F1

(pk)

pkt

res(e F (p)) = lim (p

p

k

)e

 

 

= e

 

 

 

 

 

 

 

 

 

 

=

 

 

e

 

 

 

 

 

 

 

F2(p)

 

F2(pk)

 

F

 

 

pk

p

 

pk

 

 

 

F2(p)

 

 

 

 

 

 

 

 

 

 

(pk)

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p→pk

 

 

p−pk

 

 

 

 

 

18

и, значит, в этом случае

f(t) = n F1(pk)epkt.

k=1 F2(pk)

Пример 2. Найти оригинал для изображения

1

ep F (p) = p2 .

Решение. Функция F (p) удовлетворяет всем условиям теоремы 3, так как, во-первых, она является аналитической во всей комплексной плоскости, за исключением особой точки p = 0,

во-вторых,

1

 

 

 

lim F (p) = lim

ep

= 0

 

p→∞

p→∞ p2

 

и, в-третьих, при Re p = x > 0 интеграл

x+i∞

|F (p)|dy

x−i∞

сходится по признаку сравнения (глава VII, Ÿ4, пункт 1), так как на данной прямой p = = x + yi, x > 0, y (−∞, +)

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

x

y

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

+y2

x2+y2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

p

 

 

 

 

ex+yi

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ex

 

è

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

F (p)

|

=

 

 

2

 

=

 

2

 

 

2

 

=

 

 

 

x

2

+ y

2

 

 

 

2

2

 

 

p

 

x + y

 

 

 

 

 

 

 

 

 

x + y

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

x

 

 

 

 

 

y

−∞

 

πex

 

 

 

 

 

 

e

dy =

e

 

arctg

 

 

 

 

 

=

 

 

 

.

 

 

 

 

 

x2 + y2

x

x

 

 

 

 

 

x

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Воспользуемся формулой (4). Точка p = 0 является существенно особой для функции eptF (p).

Найдем вычет в ней, разложив функцию в ряд Лорана. Применив разложение экспоненты в ряд Маклорена, получим:

1

 

p2t2

p3t3

p4t4

1

 

1

 

1

 

1

 

eptF (p) =

 

(1 + pt +

 

+

 

+

 

+ . . .)(1 +

 

+

 

+

 

+

 

+ . . .).

p2

2!

3!

4!

p

2!p2

3!p3

4!p4

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

res(eptF (p)) = t +

t2

 

t3

 

 

t4

 

 

tn

 

+

 

 

+

 

 

+ . . . =

 

 

 

 

 

 

 

(n

1)!n!

0

2!

2!3!

 

3!4!

 

 

n=1

и, значит,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tn

 

 

 

 

 

 

f(t) =

 

 

.

 

 

 

 

 

(n 1)!n!

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, оригинал для данного изображения представляется степенным рядом. Известно, что этот степенной ряд выражается через модифицированную функцию Бесселя первого

ðîäà I1(z) специальную функцию, являющуюся решением линейного однородного дифференциального уравнения второго порядка

z2y′′ + zy(z2 + 1)y = 0

(уравнения Бесселя), по формуле

(2t).

f(t) = tI1

19

Ÿ3. Применение преобразования Лапласа к решению дифференциальных уравнений

В некоторых случаях преобразование Лапласа может быть с успехом использовано для решения дифференциальных и интегральных уравнений, а также в задачах математической физики. Алгоритм применения преобразования Лапласа здесь следующий: во-первых, пользуясь свойствами преобразования Лапласа (Ÿ1), мы переходим от дифференциального или интегрального уравнения к алгебраическому уравнению относительно неизвестного изображения искомого решения, из которого оно, т. е. изображение, и находится; во-вторых, применяя методы нахождения оригинала по изображению (Ÿ2), мы по изображению восстанавливаем соответствующий оригинал, который и будет являться решением дифференциального или интегрального уравнения.

Особенно удобно использовать указанный алгоритм в задаче Коши для линейного дифференциального уравнения n-го порядка с постоянными коэффициентами:

x(n) + a1x(n−1) + . . . + an−1x+ anx = f(t);

 

(1)

x(0) = x0, x(0) = x, . . . , x(n−1)

(0) = x(n−1).

 

(2)

0

 

0

 

 

Здесь x = x(t) неизвестная функция аргумента t

 

R; a1, a2, . . . , an; x0

, x

, . . . , x(n−1)

 

 

0

0

действительные числа, f(t) оригинал.

Решение x(t) задачи Коши (1) (2) существует, единственно (глава IX, Ÿ3) и является оригиналом вместе со всеми своими производными до n-ой включительно. Последнее вытекает из общей теории линейных дифференциальных уравнений. Обозначим изображение искомого решения через X(p), ò. å. x(t) → X(p). Использовав последовательно n раз свойство дифференцирования оригинала, мы получим:

x(t) → pX(p) − x(0) = pX(p) − x0,

x′′(t) → p(pX(p) − x0) − x(0) = p2X(p) − x0p − x0,

· · · · · · · · · ·

x(n)(t) → pnX(p) − x0pn−1 − x0pn−2 − . . . − x(0n−2)p − x(0n−1).

Пусть f(t) → F (p). Применив преобразование Лапласа к обеим частям дифференциального уравнения (1), мы, учитывая предыдущие соотношения, найдем

Pn(p)X(p) = F (p) + Qn−1(p),

ãäå Pn(p) = pn + a1pn−1 + . . . + an−1p + an характеристический полином соответствующе-

го линейного однородного дифференциального уравнения, Qn−1(p) полином степени n − 1.

Отсюда

X(p) = F (p) + Qn−1(p). Pn(p)

По найденному изображению X(p) мы, использовав табличный метод или формулу Меллина, найдем соответствующий оригинал x(t), который и будет служить решением поставленной

задачи Коши (1) (2).

Замечание 1. Если требуется найти общее решение дифференциального уравнения (1), то следует считать начальные условия (2) произвольными постоянными, т. е.

x(0) = c1, x(0) = c2, . . . , x(n−1)(0) = cn.

Замечание 2. Если начальные условия (2) заданы в некоторой точке t0 ̸=,0то достаточно в дифференциальном уравнении (1) провести замену независимой переменной по формуле s = t − t0. В результате мы получим дифференциальное уравнение точно такого же вида, но

ñначальными данными в нуле.

Пример 1. Решить задачу Коши

4x′′ + x = 3 cos 2t ; x(0) = 1, x(0) = 2.

20

Решение. Пусть x(t) → X(p). Тогда

x(t) → pX(p) − x(0) = pX(p) + 1,

x′′(t) → p(pX(p) + 1) − x(0) = p2X(p) + p − 2.

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

4 (p2X(p) + p − 2)

 

3p

+ X(p) =

 

.

p2 + 41

Отсюда

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3p

 

 

 

 

 

 

4

(p2 +

 

)X(p) = 4(2 − p) +

 

 

 

 

 

 

и, значит,

 

4

p2 + 41

 

 

 

 

 

 

 

 

 

 

 

2 − p

 

 

 

 

 

 

 

 

 

 

3p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X(p) =

 

+

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

p2 + 41

 

4 p2 + 41 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

части. Для первой воспользуемся таблицей

Найдем оригиналы для каждой из дробей в правой (

 

 

 

 

 

 

)

 

 

 

 

 

 

 

оригиналов:

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 − p

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

t

 

 

 

t

 

 

= 4

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

4 sin

 

cos

.

 

p2 + 41

· p2 + 41

p2 + 41

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

Для второй дроби заметим, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3p

 

 

 

 

3

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 p2 + 1

2

 

4

p2 +

1

 

 

 

 

 

 

и так как по таблице оригиналов (

 

 

 

 

4 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 + 41

2

 

 

 

 

 

 

 

 

 

 

 

 

 

то, применив свойство дифференцирования изображения, мы получаем

 

 

 

 

 

 

 

 

 

 

 

3p

 

 

 

 

3

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t sin

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

p2 + 1 2

 

4

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, решением

поставленной задачи Коши является функция

 

 

 

(

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) = 4 sin 2t cos 2t + 34t sin 2t .

Åñëè â задаче Коши (1) (2) все начальные условия равны нулю, то может оказаться полезной формула Дюамеля, полученная в Ÿ1.

Найдем сначала в этом случае решение более простой задачи Коши

x(n) + a1x(n−1) + . . . + an−1x+ anx = 1; x(0) = x(0) = . . . = x(n−1) = 0.

(3)

Пусть x1(t) решение этой задачи, а X1(p) его изображение. Применяя преобразование Лапласа к обеим частям данного дифференциального уравнения, мы получим

Pn(p)X1(p) =

1

= X1(p) =

1

.

 

 

p

pPn(p)

Аналогично находится изображение X(p) решения x(t) исходной задачи:

X(p) = F (p) . Pn(p)

Очевидно,

X(p) = pF (p)X1(p).

Применяя формулу Дюамеля, мы находим:

t

X(p) → x(t) = f(t)x1(0) + f(t) x1(t) = f(s)x1(t − s)ds.

0

 

 

 

21

Таким образом, решение задачи Коши (1) (2) с нулевыми начальными условиями

может

быть получено по формуле

0 t

 

 

x(t) =

f(s)x1(t − s)ds,

(4)

ãäå x1(t) решение задачи Коши (3). Преимуществом этой формулы является тот факт, что îíà

избавляет нас от необходимости нахождения изображения правой части дифференциального

уравнения, что само по себе может быть сложной задачей.

Пример 2. Решить задачу Коши

x′′ + 2x+ x =

e−t

; x(0) = 0, x(0) = 0.

t2 + 1

 

 

Решение. Решим сначала задачу

x′′ + 2x+ x = 1; x(0) = 0, x(0) = 0.

Обозначая через x1(t) решение этой задачи, а через X1(p) его изображение, мы находим

1

1

1

 

1

→ x1(t) = 1 − e−t − te−t.

X1(p) =

 

=

 

 

 

p(p + 1)2

p

p + 1

(p + 1)2

Тогда по формуле (4), учитывая, что x1(t) = te−t, мы получим

t

 

t

t

(

 

)

 

 

e−s

 

 

 

 

ds

1

 

 

d(s2 + 1)

 

 

 

1

 

0

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

x(t)=

 

s2 + 1(t

 

s)e−t+sds=e−t t

 

 

s2 + 1

2

 

 

s2 + 1

=e−t

t arctg t

 

2 ln(t2 + 1) .

Аналогично с помощью преобразования Лапласа мы можем решить и задачу Коши для системы n линейных дифференциальных уравнений с постоянными коэффициентами :

 

x1

= a11x1 + a12x2 + . . . + a1nxn + f1(t),

 

 

 

 

x2

= a21x1 + a22x2 + . . . + a2nxn + f2(t),

 

·

· · · · · · · · · · ·

(5)

 

·

 

 

 

 

 

 

 

 

xn= an1x1 + an2x2 + . . . + annxn + fn(t),

x1(0) = x10, x2(0) = x20, . . . , xn(0) = xn0,

(6)

ãäå x1 = x1(t), x2 = x2(t), . . . , xn = xn(t) неизвестные функции; aij, i, j = 1, n; xi0, i = 1, n

действительные числа; fi(t), i = 1, n оригиналы.

Пусть x1(t), x2(t), . . . , xn(t) решение задачи Коши (5) (6). Обозначим

x1(t) → X1(p), x2(t) → X2(p), . . . , xn(t) → Xn(p).

Применяя преобразование Лапласа к обеим частям каждого из уравнений системы дифференциальных уравнений (5), мы для нахождения неизвестных изображений получим систему линейных алгебраических уравнений

 

 

(a11 − p)X1(p) + a12X2(p) + . . . + a1nXn(p) = −x10 − F1(p),

 

 

 

(p) + (a22

 

p)X2(p) + . . . + a2nXn(p) = x20

 

F2(p),

 

a21X1

 

 

 

 

 

 

 

 

ãäå

 

· ·

· · ·

· · · · · ·

· ·

 

 

соответственно. Ре-

 

 

 

изображения оригиналов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an1X1(p) + an2X2(p) + . . . + (ann − p)Xn(p) = −xn0 − Fn(p),

 

F1(p), F2(p), . . . , Fn(p)

 

 

 

f1(t), f2(t), . . . , fn(t),

 

шив эту систему, нам останется найти соответствующие изображениям оригиналы.

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

{

x= 2x + y + t, y= x − 2y − 3,

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

22

Решение. Если для искомого решения

x(t) → X(p), y(t) → Y (p),

òî

x(t) → pX(p) − x(0) = pX(p) 2, y(t) → pY (p) − y(0) = pY (p).

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

дифференциальных уравнений, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pX(p) 2 = 2X(p) + Y (p) + p12 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pY (p) = X(p)

 

 

 

 

 

2Y (p)

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

èëè

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(p − 2)X(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y (p) = 2 +

1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решим эту систему

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−X(p) + (p + 2)Y (p) = p.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

линейных уравнений по формулам Крамера (глава I, Ÿ5, пункт 2). Так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∆ =

 

 

p

 

 

 

2 1

 

 

 

 

= p2

5,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

p + 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 +

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 =

 

 

 

 

 

 

 

3

 

 

p + 2

 

 

 

= 2p + 4 p + p2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

2

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 +

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

òî

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 =

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1 + p + p2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2p + 4

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X(p) =

 

 

= p2 5 p(p2 5) + p2(p2 5),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y (p) =

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 5

 

 

 

 

p(p2 5)

 

 

 

 

p2(p2

5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

=

1 (5 − p2) + p2

=

1

 

 

 

 

 

p

 

 

 

 

 

 

 

1

 

 

 

,

 

 

 

 

 

 

1

=

1

 

 

1

1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

)

 

 

 

 

 

 

òî

 

p(p2 5)

 

 

5

 

p(p2 5)

 

 

 

 

5

 

p2 5

p

 

 

p2(p2

5)

 

5

p2 5

p2

 

 

 

 

 

 

 

 

2p + 4 2

 

p

 

1

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

8

 

 

 

 

 

 

 

 

p

22

1

 

 

 

2 1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X(p) = p2 5

5 (p2 5 p) + 5 (p2 5 p2 ) = 5 · p2 5 + 5 · p2 5 + 5 (p p2 ),

 

 

 

 

1

 

 

 

6

 

 

 

p

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

p

4

1

 

 

 

1 6

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y (p) =

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 (p2

 

5 p2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 5 (p2 5 p)

 

 

 

 

 

 

 

 

 

5 · p2 5 5 · p2 5 5 (p p2 )

и, следовательно, по таблице оригиналов

 

 

 

 

55 sh 5 t +

 

 

5(1 − t),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) = 5 ch 5 t +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

55

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6 + t).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(t) = ch 5 t

 

 

 

 

 

 

 

 

 

 

 

 

 

sh 5 t