Файл: Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа справочная книга.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 22.10.2024

Просмотров: 140

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

100 ОБРАЩЕНИЕ ПРИ ПОМОЩИ КВАДРАТУРНЫХ ФОРМУЛ [ГЛ. 5

коэффициенты при хп+1 в

правой и левой частях

соотно­

шения

(5.2.6),

получим

(ц-fs) (n -f s + 1)... (2n-fs) =

= ап (п +

s — 1) (n -f s)... (2п -f s 2). Отсюда

 

_

(« +

*) ( n f s f l ) ... (2n +

s — 2) ( 2 n fs — 1) (2n + s)

 

n

 

(n +

s — l ) ( n f s ) ... ( 2 n - f s - 2 )

 

 

 

 

 

_ (2n - f s) (2n + s — 1)

 

 

 

 

n f s — 1

Для

определения

bn приравниваем коэффициенты

при хп

в обеих

частях равенства

(5.2.6):

 

(n -f 1) (n -f s) (n -fs — 1) .. . (2n - f s — 1) =

=апп (n -|- s — 1) {ti “I- s)... (2n -4s — 3) -f

Отсюда

следует

-\-bn(n + s -

1) (n + s).. .(2n + s - 2).

 

 

,

(rt+ l) ( n f s )... (2n-f s —2) (2n +

s — 1) .

° n ~

( n f s - 1 )

( n f s ) ... (2 n f s - 2 )

, (2/i f s — 1) (2n + s)

n ( n f s —1) (n + s ) ... (2 n + s —3)

( n f s — 1) ( n f s — 1) (n-f s ) ... (2n-f s — 3) ( 2 n fs — 2) ~~

___ (n + 1) (2n + s — 1) .

(2n+ s — 1) (2n + s) n

_

n - f s — 1

' (2n + s —2) (n + s — 1)

(2n + s — 1) (s —2)

( 2 n f s — 2) (n -f s — 1) •

Наконец, определим cn, сравнивая, например, свободные члены в правой и левой частях равенства (5.2.6):

( _ 1)/и-1 = £я ( _ 1)я + Ся (_ !)« -!,

Сп

_

1 4 и _ 1

(2n + s — 1) (s —2)

 

1 - r

(2n + s - 2 ) ( n - f s - 1 ) ~

( 2 n f s —2) (n -fs — l) — (2 n -fs —1)(s —2) (2n -f s —2) (n -f s — 1)

____ nJ2n + s)_ _ _

(2n-fs — 2) (n -fs — 1)‘

Таким образом, можно сформулировать следующую теорему.

Т е о р е м а 2. Любые три последовательных многочлена Рп](х) связаны рекуррентным соотношением!

(2n + S- 2) (n + s - 1) Р '% 1 (х) =

 

=

[(2п -f- s) (2п -(- s — 1) (2п -}- s — 2) х —■

 

- (s - 2) (2л

+ s - 1)] Рis) (х) + п [2п + s) P{nsl 1(х).

(5.2.7)


§ 5.21

ОРТОГОНАЛЬНЫЕ МНОГОЧЛЕНЫ

101

5.2.3.Дифференциальное уравнение, решением кото­

рого являются многочлены P(ns) (jc). Для многочленов

P„s) (х) можно указать линейное дифференциальное урав­ нение второго порядка с переменными коэффициентами, которому они удовлетворяют.

В самом деле, для P„s) ^ j , как мы уже установили,

имеет место условие ортогональности Е+ /00

ыS ^ - ' П = ’( ~ ) < з ( ф ) ^ = о .

е— /со

где Q

— произвольный многочлен степени не выше п— 1.

Если

перейти

к

переменной х = 1/р,

то это условие

при­

мет вид

 

 

_L

J еих^-2P(s) (х) q {х) dx =

о,

(5.2.8)

 

 

 

 

 

 

 

 

 

с

 

 

 

 

 

 

 

 

 

где

С —окружность

радиуса

1/(2е)

с

центром в точке

х = 1 /(2 е ).

 

 

 

 

 

 

 

 

 

 

 

 

Для

вывода дифференциального уравнения рассмотрим

следующий

интеграл:

 

 

 

 

 

 

 

 

 

 

 

 

 

J = i i \ [e'IXxSp{nY (*)]' xk dx■

 

 

 

 

 

 

 

 

С

 

 

 

 

 

 

 

 

Интегрированием

по

частям

получаем

 

 

 

 

 

J

=

[xkel/xxsP ^ ’ (х)

-

2^

j xh- 4 ^ xx sP f

(л:) dx.

 

Первый

член

правой

части

обращается

 

в нуль, так

как

если

перейти

снова

к

переменной

р,

то

выражение

^ ш еРРпУ{^)

ПРИ стремлении р к бесконечности по

пря­

мой

Re р = е будет стремиться

к нулю.

Если

k —0,

то и

интеграл

/

равен

нулю.

 

 

 

еще раз

по частям:

Если

же

& > 0, то интегрируем

+ Д

e1/xxs~2PW (х) [{k -\-s— \) x k —

dx.


102 ОБРАЩЕНИЕ ПРИ ПОМОЩИ КВАДРАТУРНЫХ ФОРМУЛ [ГЛ. 5

Первое слагаемое правой части снова обращается в нуль. Рассмотрим второе слагаемое. В квадратных скобках под знаком интеграла стоит некоторый многочлен степени k, следовательно, в силу условия (5.2.8) интеграл равен нулю

при

k = \,

2,

п — 1.

Таким

образом, мы

доказали,

что

при k — 0,

1, ... , п —

1 интеграл J равен

нулю:

J =

2^- ^

 

(х)]' xk dx =

 

 

 

= 2j5 ^ [— ellxxs~2P ^' {x)-\-seil!Cxs~1P^y (x)-f-

 

 

 

-f- el/xxsP ^" (x)] xk dx =

 

 

= 2^7

e1/A:xs_2

[x2P ^" (x) +

(sx 1) РУУ (x)| xk dx = 0.

 

 

c

 

 

 

 

степени п,

Последнее равенство

означает,

что многочлен

стоящий - в квадратных скобках, ортогонален с весом eHxxs-i к xk для ^ — о, 1, ..., п — 1. Отсюда можно заклю­

чить, что этот многочлен отличается от Рns)( (х) лишь постоянным множителем уп:

x2P(nS)" (X) + (sx - 1) P f (х) = УпР ^ (х).

Для определения множителя уп достаточно сравнить в пос­ ледней формуле коэффициенты при хп:

n ( n — l)(n + s — l)...(2n-fs — 2)-f

4-sn (n -)-s— l ) .. . ( 2n 4-s — 2) = y„(« + s — 1). ..(2n + s 2).

Отсюда n(n — l)-fsn = v„, или yn — n(n + s — 1).

Таким образом, доказана

Т е о р е м а 3. Многочлен Рп'1(х), определяемый форму­ лой (5.2.1), является решением линейного дифференциального уравнения с переменными коэффициентами

x2P{nS)" (х) + (sx - 1) P nsy( ( х )- п (п + s - 1 ) P ns>( (х) = 0. (5.2.9)

5.2.4. Интегральное представление многочленов Р „*(1 х).

Покажем, что многочлены Р{п (х) имеют следующее интег­ ральное представление:

00

Р п [X) =

jj е * * -* (1 - x t f e r l d t. ( 5 . 2 . 1 0 )


§ 5.21 ОРТОГОНАЛЬНЫЕ МНОГОЧЛЕНЫ 163

Проверим это, вычислив интеграл, стоящий в правой

части

(5.2.10):

 

 

 

 

 

 

СО

 

 

 

 

Г (1 + 1 - 1)

\ t™ ~ 4 l-x ty e -‘dt =

 

 

 

о

п

 

 

. .

 

 

со

 

 

-

f (l + ‘. - i ) S (”" " а 2

< -

^

'

 

 

0

/г = 0

/

\

 

 

 

со

 

 

2

( - г -

'

" + - *

 

 

&= 0

2 ( -

'

>

о

 

=

Г И Д !

•)-*

(*) Г(2» + *- * - 1) *”-*•

Полученное выражение совпадает с формулой (5.2.1) для многочленов Рп( (х), что и доказывает утверждение (5 .2 . 10).

5.2.5.Производящая функция для многочленов Р(^ ( х ) .

Многочлены Рп\х) можно рассматривать как коэффициенты разложения в ряд Тейлора некоторой аналитической функ­ ции, которая называется производящей функцией этих многочленов.

Для

нахождения

ее рассмотрим функцию ~^=i-

Она

аналитична во

всей

плоскости

z,

кроме точек

г =

0 и

г = оо;

поэтому

в

любой точке

р

плоскости г,

кроме

указанных двух точек, ее можно представить интегралом Коши:

еР _ 1

С

ег

dz

pn+s-i ~~ 2ni

j

2 re+s_1 г р

 

I

 

 

где / — замкнутый контур,

охватывающий точку р и лежа-

 

 

qZ

 

щий в области аналитичности функции

 

Производная л-го порядка от этой функции будет

представляться формулой

т

 

 

дР_( еР

dz

(5.2.11)

dpn \

2ni

1' P)n+l'

 

Выполним преобразование

г = | 0 / Т 2 Г 7 + 1 ) .

( 5 . 2 . 1 2 )