Файл: Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа справочная книга.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п

со+ • • • +

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