Файл: Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа справочная книга.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 22.10.2024
Просмотров: 110
Скачиваний: 0
Г Л А В А 2
НЕКОТОРЫЕ АНАЛИТИЧЕСКИЕ МЕТОДЫ ОБРАЩЕНИЯ ПРЕОБРАЗОВАНИЯ ЛАПЛАСА
§ 2.1. Нахождение оригинала с помощью формулы обращения
В гл. 1 было показано, что если f(t) является оригиналом, a F (р) его изображением по Лапласу, то для вычисления оригинала по изображению можно пользоваться комплекс ным интегралом
с-\- iоо |
|
|
f(t) = 2*7 |
§ np)eptdp, |
(2Л.1) |
С — |
£ СО |
|
где с есть абсцисса в полуплоскости абсолютной сходимости интеграла Лапласа. Правда, для непосредственного вычис ления функции f(t) использовать формулу (2.1.1) затруд нительно, так как она требует знания функции F (р) для
комплексных |
значений p = c-\-iy (— оо < .у < о о ) и инте |
||
грал является |
несобственным, с колеблющимся ядром. Но |
||
поскольку |
(2.1.1) |
является интегралом от аналитической |
|
функции, |
взятым |
по контуру в комплексной плоскости, |
его можно преобразовать, применив для этого методы, известные из теории функций комплексного переменного, например изменение пути интегрирования, вычисление вычетов и т. п. Такого рода преобразования в некоторых случаях позволяют получить практически удобное выра жение для оригинала, из которого можно получить важные свойства функции, определяемой комплексным интегралом.
Методы вычисления оригинала при помощи таких пре образований комплексного интеграла (2.1.1) рассмотрим
вследующих пунктах настоящего параграфа.
2.1.1.Разложение оригинала в ряды по показательным функциям. Для одного важного класса изображений F (р) можно получить разложение оригинала в ряд, члены кото-
§ 2.1] |
н а х о ж д е н и е |
о р и г и н а л а |
п о ф о рм у л е о б р а щ е н и я |
25 |
|||||
рого соответствуют особым точкам изображений. |
А именно, |
||||||||
справедлива следующая |
|
|
F (р) мероморфнсс, |
||||||
Т е о р е м а 1. |
Пусть 1) функция |
||||||||
2) функция F (р) |
аналитична в некоторой полуплоскости |
||||||||
Re р > |
а; |
3) |
существует система окружностей |
|
|
||||
|
Сп : |
| р \= |
Rn, |
Ri < R 2< ... |
{Rn-+ со), |
|
|||
на которой F (р) |
стремится к нулю равномерно относи |
||||||||
тельно argр; |
4) |
для любого с > а абсолютно |
сходится |
||||||
интеграл |
с-\- i оо |
|
|
|
|
|
|
||
$ |
F (р) dp. |
|
|
|
|
|
|||
|
|
с — £ оо |
|
|
|
|
|
|
|
Тогда оригиналом F (р) служит функция |
|
|
|||||||
|
|
|
|
f(t) = ^ resF(p)ept, |
|
(2.1.2) |
|||
|
|
|
|
|
р* |
рк |
|
|
|
где вычеты вычисляются |
по |
всем полюсам функции |
F(p) |
исуммирование выполняется по группам полюсов, лежащих
вкольцах между соседними окружностями Сп.
До к а з а т е л ь с т в о . В условиях теоремы 1 справед лива теорема 7 гл. 1, согласно которой F (р) является
изображением функции
С-j- i ОО
£(0 = й 5 |
$ eptF(P)dp. |
(2.1.3) |
|
с — i оо |
|
Обозначим через С'п часть окружности С„, расположенную
левее |
прямой |
Re р = с, через |
с ± ibn —точки пересечения |
этой |
прямой |
с окружностью |
Сп и через Г„ — замкнутый |
контур, составленный из отрезка [c— ibn, с + ibn\ и дуги С'п и проходимый против часовой стрелки. Так как по лемме Жордана при t > О
lim \ ept F (р) dp = 0,
« - * оо „г
то интеграл в формуле (2.1.3) можно заменить следующим интегралом:
I ( 0 = П т о2й J eP‘ F(P)dp. |
(2.1.4) |
|
Теперь, применяя теорему |
ГЛ |
получим |
Коши о вычетах, |
||
/ ( / ) = Пт |
^ res F(p)ept, |
|
26 |
АНАЛИТИЧЕСКИЕ МЕТОДЫ ОБРАЩЕНИЯ |
[ГЛ. 2 |
где вычеты берутся во всех особых точках функции F (р), лежащих внутри Тп. Полученное равенство доказывает теорему.
2.1.2. Частные случаи разложения оригинала в ряды по показательным функциям. Рассмотрим случай, когда функция F (р) является дробно-рациональной функцией. Тогда имеет место
Т е о р е м а 2. Если функция F(p) = |
дробно-рацио |
нальна, причем степень многочлена А (р) меньше степени многочлена В (р), то оригиналом ее служит функция
пУ1*-1
/ ( 0 = 2 — — Ига \ - r { F { p ) ( p - P k ) nkep% (2.1.5)
где pk — полюсы F (р), a nk— ux кратности, и сумма берется по всем полюсам.
Д о к а з а т е л ь с т в о . Прежде всего заметим, что функ ция F (р) будет изображением. Это следует из теоремы о разложении дробно-рациональной функции на простейшие дроби, линейности преобразования Лапласа и справедли вости формулы
t nep«t |
я! |
|
~ { Р — Ро)п+ 1 |
||
|
(где справа |
стоит |
изображение, |
слева —его оригинал). |
||||
Следовательно, |
|
c - \ - i со |
|
|
|||
|
|
|
|
|
|
||
|
|
|
f(t) = 2ni 5 |
ept Р (Р) dp, |
(2.1.6) |
||
где c > m a x R e p b |
pk —полюсы функции F (р). |
||||||
Так |
k |
как |
и в предыдущей теореме, |
интеграл (2.1.6) |
|||
же, |
|||||||
можно |
заменить |
интегралом |
(2.1.4), ибо в силу того, что |
||||
F(p)->- 0 при p -v o o , применима лемма Жордана. |
|||||||
Применяя к интегралу (2.1.4) теорему о вычетах и |
|||||||
формулу (1.3.4) |
для вычисления |
вычетов |
в полюсах, при |
дем к формуле (2.1.5).
В частности, если все полюсы простые, то формула (2.1.5) упрощается:
/(о= k=2 |
A (Pk) Pkt |
(2.1.7) |
\ B'(Pk) |
|
§ 2.2] |
р а з л о ж е н и е о р и г и н а л а в с т е п е н н ы е р я д ы |
27 |
(мы воспользовались формулой (1.3.6) для вычисления вычетов в простых полюсах).
З а м е ч а н и е . Если многочлен В (р) имеет действительные коэф фициенты, то каждому его комплексному корню р отвечает комп лексно сопряженный корень р. Если и многочлен А (р) имеет дей ствительные коэффициенты, то тогда
|
|
(Р) |
,pt |
_ |
Л (р) |
|
|
|
|
' (Р) |
|
|
В' (р) |
|
|
|
|
|
|
„ |
А (р) |
, |
|
и, следовательно, сумма выражении —■ |
- |
ер1, вычисленная для комп |
|||||
лексно сопряженных корней |
pk |
и р к, |
будет равна 2 Re ^ j f k\ -ePkt. |
||||
o |
r |
. |
|
|
|
|
^ \Pk) |
Значит, |
формулу (2.1.7) можно в этом |
случае представить в виде |
|||||
|
|
|
|
+ 2 Re V А (р*) / Р |
|||
|
|
|
|
|
|
L |
b ' (Pk) e |
где в первом слагаемом суммирование ведется по всем действительным корням В (р), а во втором слагаемом —по всем комплексным корням с положительными мнимыми частями.
§ 2.2. Разложение оригинала в степенные ряды
Предположим, что изображение F (р) аналитично в бес конечно удаленной точке. Тогда, как известно из опера ционного исчисления, F (оо) = 0. Разложим функцию F (р) в ряд Лорана в окрестности бесконечно удаленной точки и покажем, что ее оригинал можно получить, взяв сумму оригиналов членов этого разложения. Зная, что оригиналом
функции \/рпявляется функция |
1)!, сформулируем |
следующую теорему. |
|
Т е о р е м а 3. Если F (р) аналитична в бесконечно уда ленной точке и имеет в ее окрестности разложение Лорана
СО
( 2. 2. 1)
ft —1
то оригиналом F (р) служит функция
со
|
|
(2.2.2) |
k = \ |
|
|
которая является целой функцией. |
|
|
Д о к а з а т е л ь с т в о . |
По условию теоремы |
функция |
F [р) аналитична в круге |
\ p \ ^ R . Положим |
p = \ / q и |
28 АНАЛИТИЧЕСКИЕ МЕТОДЫ ОБРАЩЕНИЯ 1ГЛ. 2
F (p) = f (— ) = Ф (?). Функция Ф (?) = ^ будет анали-
тична в круге j ? ! С 1IR, и для ее коэффициентов на осно вании неравенства Коши *) верны неравенства
( £ = 1 , 2 , . . . ) .
Из полученных неравенств для любого t находим
СО с о
k= \ |
Iс‘ 1 Й ^ а д 2 |
т |
= а д е 5 1 ' 1 (2-2-3) |
k=0 |
|
|
Отсюда видно, что ряд (2.2.2) сходится во всей пло скости /, т. е. f(t) является целой функцией переменной t.
Из неравенства (2.2.3) непосредственно следует, что для ^ > 0
|/ ( 0 1 < С е « И
Таким образом, для функции f(t) выполняется нера венство вида (1.1.3) и она является оригиналом.
Умножив ряд (2.2.2) на e~pt, получим ряд, равномерно сходящийся для всех значений t. Значит, его можно по членно проинтегрировать по 1 в пределах от нуля до бесконечности. Тогда
со со
2С*(Ёту.е~Р>dt=
оk = l
со |
оо |
|
“ 2 |
1 -pt |
cJl |
° к § (ft-l)!' |
< 4 - 1 nk > |
|
k = \ |
о |
i = l |
функция / (г) границе этого