Файл: Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа справочная книга.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 22.10.2024
Просмотров: 122
Скачиваний: 0
52 МЕТОДЫ ЧИСЛЕННОГО ОБРАЩЕНИЯ [ГЛ. 3
другого удовлетворял перечисленным выше условиям. Если будут справедливы равенства
|
lim |
pF(p) = lim/(*) = /„, |
|
|||
|
р-*-оо |
t -*■ О |
|
|
|
|
|
lim |
pF {p) = lim / (^) = |
/«, |
|
||
p -* 0 |
t -*• со |
|
|
|
||
и пределы конечны, |
то F (р) можно |
представить в виде |
||||
F W = F M + (!f + т т г ) - |
|
|||||
Оригинал для слагаемого в скобках |
вычисляется точно, |
|||||
а оригинал fx(t) |
для функции |
Fl (p) |
равен |
f(t) — f<x, — |
||
— (fo— fooje-t и, |
как |
функция |
0, |
обращается |
в нуль на |
|
концах отрезка |
[0, |
я]. Синус-ряд Фурье для функции |
||||
fx(t) будет сходиться |
быстрее, чем |
ряд для функции f(t). |
§3.2. Обращение преобразования Лапласа
спомощью ряда Фурье по синусам
Вэтом параграфе мы изложим еще один способ на хождения оригинала по значениям изображения в равно отстоящих точках на действительной оси. Этот способ основан на двух допущениях, не ограничивающих, однако, его общности. Во-первых, предполагается, что изображе ние F (р) существует при Re р > 0. Это всегда можно сделать, если рассматривать вместо изображения F (р)
изображение F (р + а) при достаточно больших а. Послед
нее равносильно |
умножению оригинала / (/) на e~at. Во- |
||
вторых, предполагается, что / (0) = 0. |
Этого можно достичь, |
||
положив /i (0 = |
/ (0 — f (О)а^, что |
равносильно |
замене |
изображения F (р) изображением F (р) — |
|
||
Интеграл Лапласа |
|
|
|
|
СО |
|
|
|
F (р) — \ e~ptf (0 di |
(3.2.1) |
|
|
о |
|
|
преобразуем при |
помощи подстановки |
|
|
|
е~ш — cos 6, |
|
(3.2.2) |
§ 3.21 ПРИМЕНЕНИЕ РЯДА ФУРЬЕ ПО СИНУСАМ 53
где а есть произвольное число, большее нуля. Тогда
m = / ( - | l n c o s 6 ) |
= cp(0) |
( о < 0 < - £ ) , |
(3.2.3) |
|
я/2 |
|
|
|
|
oF (or)= ^ |
(cos 0)р/«—»sin 6 ф (О) £/Э. |
(3.2.4) |
||
о |
|
|
|
|
Функцию ср (0) разложим в |
ряд Фурье по синусам |
нечет |
||
ных кратных дуг |
СО |
|
|
|
|
|
|
|
|
ф(0)= |
2] |
cftsin(2£-|- 1) 0, |
(3.2.5) |
|
|
k= 0 |
|
|
|
коэффициенты ск которого определяются обычным спо собом :
|
|
я/2 |
|
|
|
|
|
ск = ~ ^ ф (6) sin (2/fe + 1)0 <Й. |
|
(3.2.6) |
|||
|
|
о |
|
|
|
|
Чтобы выразить |
значения коэффициентов ск через зна |
|||||
чения |
изображения F (р), |
поступим следующим |
образом. |
|||
В интеграле (3.2.4) |
положим р = (2п-\-\)а (п = 0, |
|
1,...), |
|||
тогда |
получим |
|
я/2 |
|
|
|
|
|
|
|
|
|
|
|
o [/'(2 n -f |
1)а] = |
^ cos2" 0 sin0 ф (0) d0. |
|
(3.2.7) |
|
|
|
|
о |
|
|
|
Ядро этого интеграла может быть представлено |
в виде |
|||||
линейной комбинации функций sin (2^ +1)0: |
|
|
|
|||
cos2n 0 sin 0 = |
|
|
2n |
|
|
|
|
|
|
sin [2 (n — k) -f 1] 0, |
|
= 0. |
|
|
|
|
|
— 1 |
|
|
|
|
|
|
|
(3.2.8) |
В интеграл (3.2.7) подставим выражения (3.2.5) и (3.2.8). Так как
я/2 |
( 0 |
ц Ф у |
5 |
sin (2^ + 1) 0 sin (2v + 1) 8 п!0 = | |
^ = ^ |
то при фиксированном п останутся только члены, ДЛЯ
которых y = n — k (k = 0, ! , . . . , « ) , |
т. е. |
a f [ ( 2 « + 1) a] = 2~2п П |
Cn-kI |
*= о |
|
54 |
|
МЕТОДЫ ЧИСЛЕННОГО ОБРАЩЕНИЯ |
[ГЛ. 3 |
||
ИЛИ |
|
|
|
|
|
\12П' |
2п |
со+ • • • + |
2л |
Cn-k+ • • • + |
сл — |
|
k - \ |
||||
1\ п |
|
|
|
|
|
4я+1 |
|
oF[(2n+ 1 )а], |
или |
|
k = 0 |
<з ' 2'9» |
|
Подставив в это равенство последовательно п = О, 1, ... , получим линейную систему уравнений с треугольной ма трицей для определения коэффициентов ск:
со = ^°F (o),
42
Со 4-С! = - огГ (За),
2с0+ 3q -f с2 = ^ а/7 (5а),
Выбор значения а обусловливается величиной промежутка, для которого необходимо вычислить значение оригинала f{t), а следует выбирать малым для больших t и, наобо рот, большим для малых t.
З а м е ч а н и е . Равенством (3.2.3) функция |
ф(0) |
определена на |
|||||||
отрезке |
0 ^ |
0 < л/2. |
Тригонометрический |
ряд, стоящий справа |
|||||
в (3.2.5), |
представляет функции, обладающие двумя |
особенностями: |
|||||||
эти функции |
являются |
нечетными, |
или, |
иначе |
говоря, |
имеющими |
|||
в плоскости с декартовыми осями координат |
(0, |
<р) график, симмет |
|||||||
ричный относительно начала координат, |
и |
такими, |
что их график |
||||||
является |
симметричным относительно прямой |
0 = |
п/2. |
Тригонометри |
|||||
ческий ряд (3.2.5) дает продолжение функции |
ф (0), |
обладающее |
|||||||
указанными |
свойствами |
симметрии, |
кроме |
того, |
2л-периодическое, |
||||
и является для продолженной функции рядом Фурье. |
К. определению |
его сходимости применимы все известные теоремы о сходимости ряда Фурье.
§3.3. Обращение преобразования Лапласа
спомощью рядов по обобщенным многочленам
Чебышева — Лагерра
В предыдущих параграфах этой главы были рассмот рены методы обращения преобразования Лапласа, в кото рых оригинал f(t) находился по значениям изображения F (р) в равноотстоящих точках действительной оси. В на
§ 3.3) ОБОБЩЕННЫЕ МНОГОЧЛЕНЫ ЧЕБЫШЕВА - ЛАГЕРРА |
55 |
стоящем параграфе будет рассмотрен метод восстановления оригинала, использующий значение изображения F (р) и значения его производных в одной точке.
Пусть задано преобразование Лапласа
|
|
СО |
|
|
|
|
|
|
F{p)=\e~ptf{t)dt. |
|
(3.3.1) |
||
|
|
О |
|
|
|
|
Предположим, |
что функция f(t) удовлетворяет условию |
|||||
|
|
с о |
|
|
|
|
|
|
\e~,t l \f(t)fd t< o o , |
|
(3.3.2) |
||
где |
Я > — 1. |
О |
|
|
|
|
|
|
|
|
ряды |
||
по |
Тогда из общей теории разложения функций в |
|||||
ортогональным |
многочленам |
известно, |
что ряд |
Фурье |
||
для |
функции g{t) = t~xf{t) по |
обобщенным |
многочленам |
|||
Чебышева — Лагерра |
|
|
|
|
||
|
|
СО |
|
|
|
|
|
|
|
|
|
|
(3.3.3) |
|
|
4 = 0 |
|
|
|
|
сходится в среднем к этой функции. Это |
означает, что |
|||||
для |
частичных сумм |
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
s * ( f ) = 2 |
|
|
|
|
|
|
4= 0 |
|
|
|
|
имеет место равенство |
|
|
|
|
||
|
|
с о |
|
|
|
|
|
П т |
\ | g (t) — SN (t) |2 tle~* dt = |
0. |
|
||
|
W — со 5 |
|
|
|
|
При некоторых дополнительных ограничениях, накла дываемых на g (t), которые мы здесь не приводим *), имеет место равенство
СО |
|
g ( t ) = 2 7 7 ^ (0- |
(3.3.4) |
Обобщенные многочлены Чебышева — Лагерра ортогональ ны на полуоси 0 < х < о о с весом р (х) = хке~х и могут
*) См. Г. Сеге , Ортогональные многочлены, гл. IX, теорема
9.1.5, М., Физматгиз, 1962,
56 МЕТОДЫ ЧИСЛЕННОГО ОБРАЩЕНИЯ [ГЛ. 3
быть представлены формулами
4 Х) (o = i - r V ^ - ( e- ^ ) )
тР> /Л = |
V r(fe+ X+l) (-0"» |
к КЧ |
jL Г ( т + Ь + 1) т\ (га—т)\ |
|
т = 0 |
Свойство ортогональности многочленов L*'1(/) имеет вид
\ e ,fL'£) {t)L$)(t)dt = 0 (Ифт), |
(3.3.5) |
|
о |
|
|
со |
|
|
rk= § |
[LV (0]а dt = ^ № - +1) . |
(3.3.6) |
о |
|
|
Если функция f(t) удовлетворяет условию (3.3.2), то пре образование Лапласа F (р) будет аналитической функцией в полуплоскости R e p > 1/ 2 . В самом деле, в силу нера венства Шварца — Буняковского можно записать цепочку неравенств:
СО |
0 0 |
|
|
|
§ | e~ptf (/)! dt |
jj ^е~*е~(Rep-1) <| g у) j dt sg |
|
||
о |
о |
|
|
|
( co |
W / 2 |
Гео |
|
|
j $ /W |
2 <Rep - » {dt\ |
U tKe-‘ j g (t) j2 dt |
||
|
oo |
1 1/2 |
(oo |
|
|
5 t*r*<2 Re^ - » dt\ |
{^ |
I f (0 i2 dt |
|
|
о |
> |
lo |
|
Отсюда видно, что интеграл (3.3.1) при условии (3.3.2) сходится абсолютно и равномерно в полуплоскости R e p >
> 1 /2 . |
Следовательно, функция |
F (р) аналитична в |
полу |
плоскости Rep > 1 / 2 . |
|
|
|
Для определения коэффициентов ak разложения (3.3.4) |
|||
сделаем |
следующее. Функцию |
e~{'p~1)t разложим |
в ряд |
по многочленам Чебышева —Лагерра: |
|
||
|
0- ( р- IX. |
l*k гР) //\ |
(3.3.7) |
|
Tk Lk w> |
4=0