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

1490

.pdf
Скачиваний:
5
Добавлен:
13.02.2021
Размер:
3.68 Mб
Скачать

Сравнивая выражение (3.14) с выражением (3.9) для преобразования Лапласа функции единичного скачка, сформулируем теорему задержки (смещения):

Если функция x(t) обладает преобразованием Лапласа X (s) , то преоб-

разование Лапласа задержанной (смещенной функции) x(t −Τ) определя-

ется как

L[x(t −Τ)] = X (s)exp(st ),

(3.15)

где X (s) есть преобразование Лапласа несмещенной функции, а величина

Τ в показателе экспоненты есть смещение.

Символ L[...] в соотношении (3.15) означает оператор Лапласа, т.е.

L[...]= x(t )exp(st )dt

0

и это символическое обозначение будет использоваться в дальнейшем.

Пример 2. Найдем преобразование Лапласа суммы (разности) двух функций времени

x3 (t) = x1(t) x2 (t) , где

x1 (t) =

1 t >0,

x2 (t) = x1 (t −Τ) =

1 t ,

0 t <0,

0 t .

Функции x1 (t), x2 (t), x3 (t) изображены на рис. 3.4.

Тогда

X3 (s) = [x1(t) x2 (t)]exp(st )dt = x1 (t)exp(st )dt

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

(

 

)

dt = 1

1 exp

(

 

)

 

1 1exp

(

 

 

x

(t −Τ)exp

st

sΤ

=

st

. (3.16а)

 

1

 

 

s

s

 

 

s

 

 

)

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Соотношение (3.16) может быть найдено и прямым интегрированием:

102

X3 (s) = x3 (t)exp(st )dt =

0

= T exp(st )dt = −

1 exp(st )

 

T

=

1 1exp(sT ) .

(3.16б)

 

s

 

0

 

s

 

 

 

0

 

 

 

 

 

 

 

 

 

 

x1 (t )

1

0

x2 (t )

1

0 T

x3 (t )

1

0 T

t

t

t

Рис 3.4. Теорема линейности для преобразования Лапласа: вычитание функций единичного скачка

Нетрудно видеть, что выражение X3 (s) есть разность преобразований Лапласа функций X1(s) и X 2 (s) . Этот пример позволяет сформировать

теорему линейности (сложения).

 

 

Если L[x1(t)] = X1 (s) а L[x2 (t)] = X 2 (s) ,

то преобразование Лапласа взве-

шенной суммы (разности) функции x1 (t)

и x2 (t) равно взвешенной сумме

(разности) преобразований Лапласа этих функций:

 

L[ax1(t) ±bx2 (t)] = aX1 (s) ±bX2 (s) .

(3.17)

Отметим, что в рассмотренном выше примере были использованы как теорема сложения, так и теорема смещения.

103

Рассмотрим теперь область сходимости интеграла (3.16б). Спектральная плотность X3 (s) = 1s 1 exp(st ) не имеет особенностей для всех значе-

ний s , включая точку полюса s =0 , поскольку предел этой функции при s 0 равен (после раскрытия неопределенности)

 

 

1 exp(sT )

 

lim X3

(s) = lim

 

 

= Τ.

s

S 0

S 0

 

 

Таким образом, область сходимости интеграла (3.16б) представляет собой всю s - плоскость, поскольку, как было указано выше, функция x3 (t )

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

3.2.3. Теорема умножения на t

Теорема умножения на t формулируется следующим образом:

Если преобразование Лапласа функции x(t) есть

X (s) , т.е. L[x(t)] = X (s) ,

то преобразование Лапласа функции x1(t) =tx(t) определяется как

L[tx(t)] = −

dX (s)

.

(3.18)

 

 

ds

 

Рассмотрим теперь примеры, демонстрирующие справедливость равен-

ства (3.18).

Пример 1. Зададим функцию времени x(t )= exp(αt ), преобразование Лапласа которой определяется выражением (3.10)

X (s)= s +1α .

Найдем теперь преобразование Лапласа функции x1(t) =tx(t )=t exp(αt ):

X1(s) = t exp(αt )exp(st )dt =

t exp[(s +α)t]dt =

0

0

104

 

= −

 

t

exp[

(s +α)t]

 

+

1

 

exp[(s +α)t]dt =

 

 

 

 

 

 

 

 

 

 

 

 

 

s

+α

 

 

 

 

 

 

0

 

s +α 0

 

 

 

 

 

 

 

 

 

 

 

t

 

exp[(s +α)t]

 

 

 

 

1

exp[(s +α)t]

 

(3.19а)

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

 

.

s +α

(s +α)2

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая тот факт, что даже

при t →∞ любая степень t возрастает мед-

леннее, чем убывает функция exp(βt ), приведем выражение (3.19) к виду

 

 

 

 

 

X

(s)

=

 

1

 

.

 

 

 

(3.19б)

 

 

 

 

(s +α)2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

В общем случае

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L[tn exp(αt )] =

 

1!

 

.

(3.19в)

 

 

 

(s +α)n+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем теперь производную от функции (3.19б)

 

 

 

 

d

 

X (s) =

d 1

 

 

 

 

 

1

 

 

.

(3.20)

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

 

 

ds

 

 

 

 

 

(s

+α)

2

 

 

 

ds s +α

 

 

 

 

 

 

Сравнивая выражения (3.19б) и (3.20) нетрудно убедится, что

 

 

 

 

 

 

 

 

L[tx(t)] = −dX (s) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ds

 

 

 

 

Пример 2. Зададим единичный скачок Хевисайда

 

 

 

 

x(t) =

1 при t >0,

 

 

 

 

 

 

 

 

 

 

 

 

0 при t <0,

 

 

 

 

 

 

 

 

 

обладающий преобразованием Лапласа (3.9)

 

 

 

 

 

 

 

 

L[x(t)] = exp(

st )dt =

1.

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

s

 

Умножим функцию Хевисайда x(t)

 

на t

 

и найдем преобразование Лапласа

функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t при t > 0,

 

 

 

 

 

 

 

 

x1(t) =tx(t) =

 

 

 

< 0.

 

 

 

 

 

 

 

 

 

 

 

0 при t

 

 

 

 

 

105

Тогда

X1(s) = t exp(st )dt = −

t

exp(st )

 

+

 

1 exp(st )dt =

 

 

 

 

 

0

 

 

 

 

 

 

s

 

 

 

 

 

 

0

 

s 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

t

exp(st )

 

1

 

exp(st )

 

=

1

.

(3.21)

 

 

 

 

 

 

s

 

 

 

0

 

 

s2

 

 

 

 

 

0

s2

 

 

 

 

 

 

 

 

 

 

 

Сравнивая выражения (3.21) и (3.9), видим что

 

 

 

 

 

 

 

X1(s) =

1

 

 

 

d 1

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

.

 

 

 

 

 

 

 

s

2

 

 

 

 

 

 

 

 

 

 

 

 

 

ds s

 

 

 

 

 

3.2.4. Теоремы о преобразовании Лапласа производной и интеграла

Зададим функцию x(t) , обладающую преобразованием Лапласа X (s) и

найдем преобразование Лапласа производной dx(t) / dt . Поскольку, по оп-

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

x (t) = dx(t)

= lim

x(t) x(t δt)

 

1

dt

δt0

δt

 

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

X1

(s) = x1 (t )exp(st )dt =

δlimt0

x(t) x(t δt)

exp(st )dt.

(3.22)

 

 

0

 

 

0

δt

 

Используя теорему смещения, преобразуем выражение (3.22) к виду

 

 

1 exp(

sδt )

 

 

 

 

X1

(s) = lim

 

 

X (s)

(3.23)

 

δt

 

 

 

δt0

 

 

 

 

 

и раскроем неопределенность типа 0/0, имеющую место при δt 0 . Тогда

X1

(s) = lim[s exp(sδt)]X (s) = sX (s).

(3.24а)

 

δt0

 

106

Преобразование Лапласа производной порядка n имеет вид

 

L[ d n x(t)

] = sn X (s),

(3.24б)

dtn

 

 

где X (s) есть преобразование Лапласа функции x(t) .

Поскольку исходная функция x(t) есть первообразная (неопределенный интеграл)

x(t) = x1 (t)dt,

 

 

то из выражения (3.24а) следует, что

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

интеграла

может быть определено в виде

 

 

 

 

t

 

X (s)

 

 

L x(τ )dτ =

 

.

(3.25)

s

0

 

 

 

3.3.Обратное преобразование Лапласа

Рассмотрим теперь возможность перехода от спектральной плотности

X (s) к функции времени x(t) , определенной на интервале t (0,)

с ис-

пользованием обратного преобразования Лапласа (3.8) [8, 16, 17]

 

(1)

 

 

1

σ + j

 

x(t) = L

 

[ X (s)] =

 

X (s)exp(st )ds.

 

 

2π j

 

 

 

 

 

σ j

 

Это выражение можно переписать в виде:

 

x(t) =

 

1

X (s)exp(st)ds,

(3.26)

2π j

 

C

 

 

 

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

107

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

относится именно к этому классу.

Запишем правильную рациональную функцию

X (s) =

a

N

sn + a

n1

sn1

+... + a

, (n m)

(3.27)

 

 

 

0

 

 

sn +b

 

 

+...b

 

 

 

 

sm1

 

 

 

 

 

m1

 

0

 

 

и, полагая, что все корни полинома знаменателя различны (так называе-

мый случай простых полюсов), представим функцию X (s)

в виде суммы

простых дробей:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X (s) =

a

sn + a

n1

sn1 +... + a

=

k

 

+

k

2

 

+...

+

k

m

 

.

(3.28)

n

 

0

 

 

1

 

 

 

 

 

(s sP1 )(s sP 2 )...(s sPm )

 

 

 

 

 

 

 

 

 

 

 

 

s sP1

s sP 2

 

s sPm

 

В выражении (3.28) величины sPi

(i =1,

2, ..., m) есть корни полинома зна-

менателя (т.е. полюсы), а константы ki

 

(i =1, 2,

..., m)

называются вычета-

ми. Для случая простых полюсов вычеты определяются в виде

 

 

 

 

 

ki =[X (s)(s sPi )]

 

s=sPi

 

 

 

 

 

 

 

 

 

(3.29а)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или, в некоторых случаях, как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ki = lim [X (s)(s sPi )].

 

 

 

 

 

(3.29б)

 

 

 

 

SSPi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отметим, что нормировка функции x(t) (см. выражение (3.27)) выпол-

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

Подставим теперь разложение (3.28) функции x(t) в выражение (3.26) для обратного преобразования Лапласа. Тогда

108

x(t) =

k1

exp(st )

ds +

 

k2

 

exp(st )

ds +... +

km

exp(st )

ds. (3.30)

 

 

 

 

 

 

 

 

 

2π j

C s sP1

2π j

 

C

 

s sP2

 

2π j

C s sPm

Каждый из интегралов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

C

 

exp(st )

ds,

 

 

 

 

 

 

 

 

 

 

2π j

 

s sPi

 

 

 

 

входящий в выражение (3.30), может быть вычислен с использованием интегральной формулы Коши [18]

1 f (z)

2π j Cz z0 dz = f (z0 ) .

Таким образом, выражение (3.30) даёт представление функции времени x(t) в виде

x(t) = k1 exp(sP1t )+ k2 exp(sP2t )+... + km exp(sPmt ), t > 0.

(3.31)

Выражение (3.31), а также его обобщение для случая кратных полюсов, называют теоремой разложения Хевисайда (или просто разложением Хевисайда).

Рассмотрим несколько примеров использования разложения Хевисайда для определения функции x(t) по её спектральной плотности X (s) .

Пример 1. Зададим спектральную плотность в виде

X (s) = ss2++2s = s(ss++21) = ks1 + sk+2 1.

Эта функция имеет два полюса sP1 = 0 и sP2 = −1. Найдём вычеты:

k1 =[x(s)s]

 

=

s + 2

 

 

 

 

 

= 2,

 

S =0

s +1

 

 

 

 

 

 

 

 

 

S =0

 

 

k

2

= x(s)

(

s +

1

 

 

 

=

s + 2

 

= −1.

 

 

 

 

 

 

 

 

 

)

 

 

S =−1

 

 

s +1

 

 

 

 

 

 

 

 

 

 

 

 

S =−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

109

Таким образом

x(t) = 2 exp(t ) при t > 0 .

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

X (s)= 2 exp(t ) exp(st )dt =

0

= −2s exp(st )0 + s 1+1exp (s +1)t 0 = s(ss++11) ,

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

Пример 2. Пусть спектральная плотность X (s) имеет вид

X (s) =

 

ω0

=

ω0

.

s2

+ω02

(s + jω0 )(s jω0 )

 

 

 

Эта функция имеет комплексно-сопряженные полюсы sP1 = − jω0 и

sP2 = jω0 .

 

 

 

 

 

 

 

 

 

 

 

 

Соответствующие вычеты определим как

 

 

 

 

 

 

 

k1 =[x(s)(s + jω)]

 

 

 

=

 

 

ω0

|S =− jω = −

1

,

 

 

 

 

 

 

S =− jω0

 

s jω0

 

2 j

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

k2 =[x(s)(s jω)]

 

 

=

 

 

ω0

 

|S = jω =

1

.

 

 

 

 

 

 

 

S = jω0

s

+ jω0

2 j

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку полиномы числителя и знаменателя правильной рациональной функции X (s) имеют только действительные коэффициенты, то полю-

сы, равно как и вычеты, имеют вид комплексно-сопряженных пар.

Запишем теперь функцию времени

x(t) = 21j exp( jω0t )exp(jω0t ) =sinω0 , t > 0 .

110

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

3.2.1 и имеет вид X (s) = ω0 .

s2 +ω2

Пример 3. Рассмотрим теперь спектральные плотности, обладающие кратными полюсами. Пусть спектральная плотность имеет вид

X (s) =

1

,

(s +α)2

 

 

откуда следует, что функция X (s) обладает кратными полюсами с кратностью K = 2 : sP1 = sP2 = 0 .

Для оценки формы функций времени, соответствующих преобразованиям Лапласа с кратными полюсами, необходимо найти функцию, преобра-

зование Лапласа которой имеет вид x(s) =

1

где n – кратность по-

(s +α)n

 

 

люса.

 

 

Напомним теорему об умножении на t (см. п. 3.2.3), в соответствии с которой преобразование Лапласа функции вида t exp(αt ) определяется в

виде

 

 

 

 

 

 

 

 

 

 

L t exp(αt ) =

 

1

 

 

.

 

 

(s +α)2

 

 

 

 

 

 

 

 

В общем случае преобразование Лапласа произведения

tn exp(αt ) мож-

но записать как

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

n!

 

 

 

 

 

 

 

 

 

 

 

L t

 

exp(αt ) = (s +α)n+1 .

 

 

Таким образом, формула

обратного

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

Лапласа функции

x(s) , обладающей полюсом кратности n +1, может быть записана как

111

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