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