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

книги / Методы математической физики и задачи гидроаэродинамики

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
14.42 Mб
Скачать

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

У2(/0=* J v~aty r'dt—[е~"а/у' (О]о + а [e“ tt/# (0 1 * + а2 J &-aty(t)dt,

•V

откуда при условии, что производная */'(7) удовлетворяет огра­ ничению (7.16), т. е. также растет не быстрее экспоненты, имеем

У2(Р) = o?Yip ) - a y { 0 ) - g ' (0).

(7.21 )

Поступая аналогично, методом полной индукции легко устано­ вить, что

Г*(/»=а*Г (р) - а*-1у (0) - а

( 0) - . . . - ау(*-2>(0) -

у/*-1) (0),

£ = 1 ,2 ,..., я,

(7.22)

и при нулевых начальных данных

 

Уь(р)—ь*У(р),

£ = 1 ,2 ,...,л.

 

Таким образом, операции дифференцирования оригинала y(t) соответствует при преобразовании (7.12) алгебраическая опера­ ция умножения изображения на параметр а. Если подставить выражения (7.22) в левую часть операторного уравнения (7.4), то получим

[ а

п 4

- л 1а Л—1+ а .п^. а. -+\ - а п\ Y +

Ф

= F ,

 

где ^ - 1 = ^ - 1

(а) — многочлен степени, не выше п— 1

относи­

тельно а (р)

и, следовательно,

 

 

 

 

Г =

__________ Р -& П -1__________

#

(7.23)

 

 

 

а^аР + аха 7*” 1 + . . . + Дя- 1 а

+ а п

*

 

т.е. изображение решения задачи Коши (7.1), (7.10) найдено. Теперь нужно осуществить обратный переход, т. е. по изве­

стному изображению Y(p) найти соответствующий ему ориги­ нал y(t), который будет являться решением задачи Коши (7.1), (7.10). Иными словами, нужно решить задачу обращения пре­ образования (7.12). Выбор функции а (р), которая пока была произвольной аналитической функцией, как раз и определяется

условиями

обращения преобразования (7.12). Используя обо­

значение

(7.14), преобразование (7.12)

можно переписать

так:

 

 

 

Г ( 5 + /о ) = |

(7.24)

Преобразование (7.24) напоминает известное преобразование Фурье, которое ставит в соответствие некоторой функции <р(t)

181

другую функцию Ф(а) по правилу

Ф (з)=

j

(7.25)

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

(7.25), как известно (см.

(3), (4)],

имеет вид

 

 

!

+ °*

(7.26)

? ( 0 = - ^

] Ф(«)е“ * .

Преобразование (7.25) существует, если <р(t) на любом конеч­

ном интервале непрерывна, или кусочно-непрерывна,

и абсо-

+°®

сходится.

лютно интегрируема на (— оо, +оо), т. е. J \y{t)\dt

Формула (7.26) определяет обращение функции Ф(сг) в каж­ дой точке непрерывности tp(t), С целью получения обращения преобразования (7.24) с помощью обращения преобразования Фурье необходимо:

1) заменить нижний предел /= 0 в интеграле (7.24) на t = —оо, что возможно, если потребовать

y ( t) = 0 при t < 0 .

(7.27)

Требование (7.27) допустимое, поскольку решение задачи Коши однозначно определяется при tf>0 начальными условиями (7.10) н, следовательно, при /< 0 можно доопределить решение произвольно;

2) положить «2 (s, а )= а , и так как аь аг — соответственно

действительная и мнимая части аналитической функции а(р), то тем самым определяются à\(s, a )= s и а(р) = s - f io=p. При этом преобразование (7.12) принимает вид

 

К (/> )= Г е-" 0(/)< #,

 

 

(7.28)

а условие

(7.17) его существования трансформируется

в

нера­

венство —s+So<0 или s>5o» т. е. интеграл

(7.28) существует

в правой

полуплоскости s>So комплексной

плоскости

р

(рис.

7.1);

3) зафиксировать параметр s, положив s = s * > s 0-

При выполнении указанных трех условий

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

(7.24) принимает вид

 

+»•

 

K ( s « + /9 ) = f ÿ(/)e-**'e-'«W

(7.29)

й

является,

очевидно,

преобразованием Фурье функции

<

р

Н

о

тогда

по .формуле обращения (7.26)

+«•

j r(s* + /0)e'»'rf«

----00

y ( t ) = J - С K(s*-l-w)e('+«V«.

di <

УА

5*+ ioo

 

 

 

 

t(t),

 

 

\

 

1-

 

 

I

 

 

 

0

/.|S0

5*

"S

 

 

S■

5*- too

О

t

 

А

Рис, 7.1 Рис. 7.2

Переходя к первоначальной переменной p = s+ io и учитывая, что под интегралом s=s* фиксировано, в силу чего dp=ido, окончательно получаем

со

? (* )= -Ц-

f y(P)epldp, s*>s„,

(7.30)

2я/

J

 

s * — I со

где интегрирование ведется в комплексной плоскости р вдоль любой вертикальной прямой, проходящей правее точки s0 (рис. 7.1).

Преобразование (7.28) называется преобразованием Лапла­ са, интеграл (7.28) — интегралом Лапласа, формула (7.30) дает обращение преобразования Лапласа (в точках непрерывности оригинала y(t)).

При построении преобразования Лапласа и его обращения на оригинал y(i) были наложены ограничения (7.16), (7.27). Если к этим ограничениям добавить упоминавшееся выше требо­ вание непрерывности или кусочной непрерывности функций y(t), Ф(/), то окончательно определение оригинала, не обязательно являющегося решением уравнения (7.1), выглядит следующим

образом.

Функцию y(t) назовем оригиналом по Лапласу, если

1) y(t) непрерывна на (0, оо) или кусочно-непрерывна на любом конечном интервале, принадлежащем (0, оо) ;

2)при /> 0 выполняется условие (7.16);

3)при отрицательных t выполняется условие (7.27). Первые для условия обеспечивают существование преобра­

зования Лапласа (изображения), третье связано с возможно­ стью его обращения.

Функцию Y(p), определяемую равенством (7.28), называют

изображением по Лапласу оригинала y(t). В дальнейшем бу­ дем обозначать оригиналы строчными буквами y(t), f(t) и т. д., а соответствующие им изображения — теми же пропис­ ными буквами Y(p), F(p) и т. д. Принята символическая запись Kt)^F(p), означающая, что «оригинал f(t) имеет изоб­ ражение F(p)»j а также F(p)^f(t), означающая, что «изобра­

жение F(p) имеет оригинал

является (рис. 7.2) еди­

Примером простейшего оригинала

ничная функция

 

1 при / >

О,

О при

0.

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

/ ( / ) при

* > 0 ,

0 при

0

будет удовлетворять всем трем условиям и, следовательно, яв­

ляется оригиналом. В дальнейшем,

рассматривая в качестве

оригиналов функции const, tk, sinotf,

cos ой, ew и т. д.,

мы всег­

да будем предполагать, что они умножены на i\(t).

 

Найдем теперь непосредственно изображение функции еи,

где A,=XJ2 — комплексное число.

Предварительно

устано­

вим показатель роста этой функции. Имеем |ew | = е х*11е/х** |= е х,/,

откуда видно, что функция е** имеет показатель роста sQ—%j. Тогда изображение оригинала е*/ существует в полуплоскости

R ep=s>X j. Обозначая преобразование Лапласа с

помощью

оператора L и учитывая, что R ep>s>À i, получаем

 

 

 

о

1

L

ex/e-'/d t = j srip-Wdt

t

 

 

/ 7

- Х

поскольку при S>Ai И /-»-оо

e - ( s - x t)t

 

 

е-(р -Х )/

 

 

Р — Ь

—ч

 

Итак,

 

 

 

 

 

 

 

 

 

 

 

 

 

ех/_л_

*

_

 

 

 

(7.31)

 

 

 

 

р - Х

 

 

 

 

 

 

 

 

 

 

 

 

 

и, в частности,

 

 

 

 

 

 

 

 

 

 

 

 

-

7

-

 

 

 

 

(7.32)

 

 

 

 

 

 

 

 

Установим некоторые фундаментальные свойства интеграла

Лапласа

(7.28), сходящегося, как

было установлено

выше, в

правой полуплоскости s>so.

интеграла

Лапласа.

Для

любого

Равномерная

сходимость

оригинала

f(t)

при

s > s 0 имеем

(см.

(7.18) и (7.19),

где сле­

дует положить ai = s

и заменить у

на f)

 

 

 

 

 

00

 

 

09

 

 

 

 

 

 

 

j1

 

< JW„ ^e~,s~So>w= -JZ*r-

(7.33)

и если s ^ s i > s 0, то

 

 

 

 

 

 

 

 

 

 

 

0

j e-<*i~~So)*df=

Mo

- •

(7.34)

 

 

 

^ М

 

 

 

Sl — So

Как видим, интеграл Лапласа при s>So мажорируется сходя­ щимся интегралом, зависящим от параметра р (неравенство (7.33), где s = Rep), а при s > s i> s 0 мажорируется сходящимся интегралом, не зависящим от параметра р [неравенство (7.34)]. Следовательно, интеграл Лапласа не только сходится при s>s<j»

но и равномерно сходится при s^Si>So-

Аналйтичность изображения. Из курса теории функций ком­ плексного переменного известно, что если функция ф(р, i) аналитична по р и кусочно-непрерывна по t для всех р в одно­

связной области D и всех (0, оо)

и интеграл

00

 

Ф ( /0 = |

О*#

о

 

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

в полуплоскости Rep^Si>So, кусочной непрерывности оригина­ ла f(t) и аналитичности функции е~р* следует, что F(p) аналитична в указанной полуплоскости и

т

F'(p) = j / V)

(7.35)

О

или

т

(7.36)

О

Асимптотическое поведение изображения на бесконечности. Переходя в неравенстве (7.33) к пределу при R ep = s-> оо, по­ лучаем

Шп F (p )= 0.

(7.37)

Пере­

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

$ 7.2. ПРОСТЕЙШИЕ СВОЙСТВА ПРЕОБРАЗОВАНИЯ ЛАПЛАСА

Iе. Л и н е й н о с т ь .

Если f(t), cp( t ) оригиналы,

причем

f(l)=F(ph ^ ) ф Ф ( р ) ,

то

 

c\ f V)+

с2?(О Ф CiF{p)+ с2Ф(/7),

(7.38)

где ch с2произвольные постоянные.

Свойство линейности интегрального преобразования общего вида (7.2) отмечалось выше и с его помощью было получено изображение (7.4) дифференциального уравнения (7.1). Ис­ пользуем это свойство для получения изображений тригоно­ метрических и гиперболических функций.

С помощью изображения

(7.31)

показательной

функции

и

свойства (7.38)

получаем

 

 

- L

-

U

.

1

sin t = —— (ef< —е- //)

 

1

2i (p — i

21

 

1

P + i )

 

P2 + 1

 

Аналогично получаем

изображения

оригиналов

cos t,

sh t, ch t.

Окончательно имеем

 

 

 

 

 

 

 

sin/ = _ J ___

c o s /=

P

 

 

 

 

 

P

 

P2 "H 1

 

P2 + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.39)

2°. П о д о б и е . Если f(t) оригинал и f(t) ==F(p),

то f(at),

где а любое положительное число, также является оригина­ лом и

/ ( a ^ î r F ( v ) -

<7'40>

Действительно,

I l f (a t ) ] = ^ f { a t ) e r * 4 l t =

1

а

что соответствует формуле (7.40).

Пример. На основании формул (7.39), (7.40) находим изоб­ ражение sin at:

sinctf=—-------- !------- .

*a (p /а)2 + 1

Аналогично находим изображение cos a t, sh at, chat. В резуль­ тате имеем

sin at-

a

 

 

 

 

ch at= — - — .

(7.41)

 

pi — a2

 

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

Если f'(t)

оригинал (тогда f(l) также оригинал) и если f(t)=F(p), то

 

- f'( t) = p F ( p ) - f ( 0 )

(7.42)

и, вообще,

 

 

/ <ft) (/)== P*F (P) —/J*-1 / (0) — ph~2f ’(0) —... —/? /(ft_2) (0)—

 

(0).

(7.43)

Соответствия (7.42), (7.43) получены выше в более общем виде (7.20), (7.22). Если в этих формулах положить а (р )= р и за­ менить обозначения у, Y соответственно на /, F, то придем к соотношениям (7.42), (7.43), которые, как и (7.20), (7.22), по­ казывают, что (в случае /(0) = ... = /< Л_1>(0) = 0 ) операции дифференцирования оригинала соответствует при преобразова­ нии Лапласа умножение изображения на соответствующую сте­ пень р.

Найдем изображение по Лапласу решения задачи Коши (7.1), (7.10). Полагая в формуле (7.23) a *=р, имеем

Г(Р)=*

_________Р (Р ) -

& п -1 (Р )

(7.44)

айр п + щ р п- 1 +

. . . + а я_ ,р + аа

где ^ л -1 (р) — многочлен степени не выше п— 1; F{p) — изобра­

жение правой части / уравнения (7.1).

4°. Умножение оригинала на минус аргумент (дифференци­ рование изображения). Если }{t) — оригинал а если f(t)= F (p )t то

- t f ( t ) = F ' ( p ).

(7.45)

Действительно, записав изображение для —tf(t) и исполь­ зуя формулу (7.36), получаем

£ ! - * / ( * ) ] =

e - ^ = F ' ( p ) ,

что соответствует формуле (7.45). Формула (7.45) показывает, что умножению оригинала на минус аргумент соответствует дифференцирование изображения. Поэтому

( - t)* f(t)= F W {p ) .

(7.46)

Отметим, что в результате умножения оригинала на степенную функцию ik снова получаем оригинал.

Пример. Известно, что ех/==— !—

,

где Я — комплексное

 

 

 

 

 

р — а

 

 

 

число [см. (7.31)].

Умножая

 

экспоненту на (—t)k, на основа­

нии соотношения

(7.46)

получаем

 

 

 

 

( —/)* ех* =

 

 

 

 

( — ! ) ( —2) . . . ( — *)

(7.47)

d p *

Р — ^

 

 

-

Х)*+1

откуда

 

 

 

 

 

 

 

 

 

 

 

 

и, в частности,

7 ,

 

 

* =

0' 1' 2......

 

 

 

 

 

 

 

 

 

 

 

/ft —

^

1

 

_ri 1

о

•••

(7.48)

 

 

* ^—vj Ij

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

 

о р и г и н а л а . Если f ( t ) непре­

рывный оригинал и /(^)==F (р), го

 

 

 

 

 

 

f ü ) d x

 

F(p)

 

(7.49)

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

Действительно, пусть

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

| / ( т ) Л

=

ф(/>),

 

(7.50)

 

о

 

 

 

 

 

 

 

 

гле Ф(р) подлежит определению.

Дифференцируя

оригинал

(7.50), получаем

 

 

 

 

 

 

 

 

 

/ ( T ) d x

рф{р)—^/(х)ах,

0

/ У ) ф р ф ( р ) .

Но по условию f(t) ~=F(p), поэтому рф(р) = F(p) и Ф(р) = *=F(p)/p, что и требовалось доказать. Из соотношений (7.42), (7.49) видно, что если дифференцированию оригинала соответ­

ствует умножение изображения на параметр р, то обратному действию — интегрированию оригинала — также соответствует

обратное действие — деление изображения на р.

 

6°. Д е л е н и е

о р и г и н а л а на

а р г у м е н т . Если f(t)/t

оригинал и если f ( t = F ( p ) , j o

 

 

 

 

J M . =

^F{g)dq,

(7.51)

 

 

Р

 

 

где (рис. 7.3)

 

 

 

М

 

Ит

(7.52)

f F {q)dq=

Г F(q)dq.

J

R e ^ + «• i

 

P

 

 

p

 

Действительно, пусть

 

 

 

 

М О . ~ Ф ( р ) ,

(7.53)

где ф (р) подлежит определению. Умножая оригинал f(t)/t на минус аргумент, на основании соотношения (7.45). получаем —/(*)==:ф'(р). Но по условию f{t)=F{p), в силу чего Ф'(р) = е=з—F(p). Интегрируя это тождество, находим

&

<НЗ>)-<ь{р)= — J F(g)dq.

Р

Отсюда

при

Re

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

равенства

(7.37)

lim

ф ( ^ ) = 0 ,

получаем

 

 

 

 

00

 

 

 

 

<Hp)=^F{q)dq,

(7.54)

где интеграл понимается в смысле (7.52). Из соотношений (7.54), (7.53) вытекает соотношение (7.51) и требуемое утверж­ дение доказано.

Пример. Известно, что sin/

--------- . Отношение

sin t

t

 

/>2 + 1

1 Ё

оригинал, так как lim - sln— = 1. Найдем изображение функции /-*.о t

sin t . Согласно свойству (7.51), t

sin t

dq

sin t

= - | - —arctg p.

q2 + \

t

 

 

sin i

 

(7.55)

 

= arct gp.

 

~T~

Интегрируя оригинал

in /

согласно свойству (7.49), полу-

— •

С

sin т . . 1

.

I

------- dx ==—

arctg p.

J *

P

о

 

Последний интеграл называют интегральным синусом и обозна­ чают

 

sU = |- 2 ^ - r f t .

(7.56)

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

 

 

s I / = — arctg/>.

(7.57)

 

P

 

и

7°. З а п а з д ы в а н и е о р и г и н а л а . Если H t ) — оригинал

го

 

 

/ ( / _ т ) = е - ^ ( р ) ,

(7.58)

где х любое положительное число. Оригинал /( / х)

называ­

ют запаздывающим (рис. 7.4).