книги / Методы математической физики и задачи гидроаэродинамики
..pdfи интегрируя дважды по частям, получим
У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), то получим
[ а |
0а |
п 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, оо) ;
Итак, |
|
|
|
|
|
|
|
|
|
|
|
|
|
ех/_л_ |
* |
_ |
|
|
|
(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 и всех tŒ (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), получаем |
|
|
|
|
|
|
|
|
|
Отсюда |
при |
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). |
|