Файл: Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа справочная книга.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 22.10.2024
Просмотров: 141
Скачиваний: 0
5. ■ТЕОРИЯ КВАДРАТУРНЫХ ФОРМУЛ $5
Умножим |
равенства |
(5.1.11) |
на Г (s + n + /n— 1) и запи |
|||||||||
шем полученную систему в виде |
|
|
|
|
|
|||||||
s-\-n— 1 + Gi + (s + ^ — 2 )a 2- f . . . |
|
|
|
|
|
|||||||
|
|
... |
(s |
|
n — 2) (s |
n — 3)... san= |
0, |
|
|
|||
s_j_rt + ai + (s + n ~ 1) a2+ • • • |
|
|
|
|
(5.1.12) |
|||||||
... + |
(s + n — 1) (s -j- n — 2 )... (s + |
1) an= 0, |
||||||||||
|
|
|||||||||||
's+ 2n- 2 |
+ |
G1 + (s + |
2n — 3) a2 + ... |
|
|
|
|
|
||||
... + (s + 2n — 3) (s + |
2n — 4).. ,(s + « — 1) an = 0. |
|
|
|||||||||
Определитель этой |
системы есть |
|
|
|
|
|
||||||
1 |
s + n —2 |
... [(s+ n —2)(s + n —3) |
|
... |
s] |
|
|
|||||
А = 1 |
s + |
n — 1 |
... [(s + n — 1) (s + n — 2) |
... |
(s+1)] |
|
||||||
1 |
s + |
2n —3 |
... [(s + 2n —3) (s+ 2n —4) |
... (s + n — 1)] |
||||||||
Достаточно |
убедиться в том, |
что А ф 0, |
так как |
тогда |
||||||||
система |
(5.1.12) будет иметь решение alt |
а2, .... |
ап и |
|||||||||
только |
одно. |
Рассмотрим систему |
п функций |
|
xs+n~2, |
хS+,I_1.......... xs+2«-3_ Они линейно независимы на любом отрезке, не приводящемся к точке. Построим линейное дифференциальное уравнение порядка п, для которого эти функции образуют полную систему независимых ре шений:
У |
У' |
... y<n) |
2 |
(xs+n~2y |
(д-5+Л-2)(Ш |
Х^+п-1 |
(x s + n - iy |
(x s+ n~iyn) |
j^J+2/1-З |
(^ s+ a n - з у |
(^s+2ra-3)frt |
Если разложить определитель |
по элементам первой строки |
||
и разделить обе части уравнения на |
то последнее |
||
уравнение можно записать так: |
|
|
|
п) + с?дг"-У |
+ . . . ■+ спу = |
0, |
(5.1.13) |
где с? —const. |
|
|
|
96 ОБРАЩЕНИЕ ПРИ ПОМОЩИ КВАДРАТУРНЫХ ФОРМУЛ [ГЛ. 5
Действительно, рассмотрим алгебраическое дополнение элемента у:
(x s+n~2)' |
( x s+n~ 2 )" |
. . . |
(x S+n-2yn> |
(^■S+n-1 ) ' |
^ S + 4 - l j '( |
_ |
(x s + n - i y m |
( Xs + 2 h - 3 y |
(x S + 2 a - 3 y ' |
. . . |
(х З -г 2 П -З у т |
Из элементов первого, второго, ..., п-го столбцов выне
сем за знак определителя соответственно xs+”~3, xs+"~4, ...
... , |
xs~2. Затем из элементов строк оставшегося определи |
||||||||
теля |
вынесем |
множители |
1, х, ... , |
х ^ 1. R итоге полу |
|||||
чим определитель xn(n+s~3) D, |
где D |
равно |
|
||||||
s + n — 2 |
(s + |
n —2)M |
.. |
[(s + |
n —2) |
... |
( s - i ) ] |
|
|
S - j - t l — 1 |
(s + |
л — l) M |
.. |
[(s + n — 1 ) |
... |
s] |
= С Л |
||
|
|
|
|
|
|
|
|
|
|
s+ 2 n —3 |
(s + 2 n - 3 )M ... |
[(s+ 2n —3) |
. .. |
(s + n —2)] |
|
и буквой М обозначен множитель, на единицу меньший предыдущего.
Совершенно аналогично устанавливается, что алгебраи
ческие дополнения |
элементов у1, у".........у(п) равны соот |
||||||||||
ветственно |
|
|
|
|
|
|
|
|
|
|
|
|
С П ^ д /л (5 + л -3 )+ 1 ) |
с п |
|
3)+2 |
|
£ П х п ( з * п - % ) ^ п |
|||||
Этим доказывается разложение (5.1.13). |
Эйлера и |
||||||||||
Уравнение (5.1.13) |
является |
уравнением |
|||||||||
имеет две |
особые |
точки: х = 0, |
х = |
оо. |
|
|
|||||
Выпишем определитель Вронского для решений xi+n~2, |
|||||||||||
xs+n- \ ... , |
xs+2n~3 этого |
уравнения: |
|
|
|
||||||
Г (х 5+Я-2, |
.... |
xs+2n~3) = |
|
|
|
|
|
|
|||
|
s + n - 2 |
(s + |
n — 2) X s + n ~3 |
. |
. [(s + |
n —2) |
. .. sxs-1] |
|
|||
■ |
S+n-1 |
(s + |
n — l)x s+n_2 |
. |
. [(s + |
n - 1 ) |
. ••(s + l)* 4 |
|
|||
-s+ая-й |
(s + |
2n — 3) X s t 2« - i. |
. [(s + 2n —3). . ( s + n — l)xs4n"2] |
||||||||
|
|
|
|
|
|
|
|
|
|
|
(5.1.14) |
Так |
как |
решения |
х*4”-2, |
х5+я“1......... |
х542"-3 |
уравнения |
(5.1.13) |
линейно |
независимы, то определитель (5.1.14) |
|||
может |
равняться |
нулю только в особых точках |
уравне |
||
ния, т. е. в точках |
х = 0 и х = |
оо. В остальных |
точках |
||
он отличен от нуля, |
в частности, |
и в точке х = 1. |
Но при |
§ 5.2] |
ОРТОГОНАЛЬНЫЕ МНОГОЧЛЕНЫ |
97 |
х= 1 определитель W (xs+n~2, ... , xs+2n~s) совпадает с опре делителем А, и поэтому А Ф 0, следовательно, система (5.1.12) имеет единственное решение.
Таким образом, существование и единственность мно
гочлена со„ |
доказаны. |
Р
Так как весовая функция в (5.1.7) зависит от пара
метра s, то и многочлен со„^—J тоже будет зависеть от s.
Будем обозначать его “ «’ (у )- Чтобы закончить исследо
вание возможности построения |
формулы (5.1.5), точной |
||
для |
многочленов степени 2п— 1 |
от 1/р, |
необходимо пока |
зать, |
что все корни многочленов со® |
при любых s > О |
лежат в правой полуплоскости. Этот вопрос будет рас смотрен в следующем параграфе.
Квадратурная формула (5.1.5) наивысшей степени точ ности для интеграла (5.1.4) в частном случае s = l была построена Г. Солзером (см. [8]).
§5.2. Ортогональные многочлены, связанные
сквадратурной формулой наивысшей степени точности
5.2.1. |
Явное |
выражение |
многочленов |
|
Для |
||
получения явного выражения ю® |
рассмотрим следую |
||||||
щий многочлен степени п: |
|
|
|
|
|
||
р п ( } ) = |
( - 1Ye-Pp^s-' |
(е”р-п- ^ ) = |
|
|
|||
|
(п\ |
(— l)"-ft( tt+ s — 1) |
... (ra+ s + fe-2) |
, |
(5.2.1) |
||
=£2= 0U. |
|
|
|
|
|
|
|
который |
мы записали в виде, |
аналогичном формуле Род- |
|||||
рига для многочленов Лежандра. |
|
|
|
||||
Покажем, что для этого многочлена выполняется усло |
|||||||
вие ортогональности |
|
|
|
|
|
||
8 - j - i 0 0 |
|
|
|
|
|
|
|
§ epp-sP f(~ ^ p -mdp = 0 |
|
(m = 0, 1 ,..., п - |
1), |
(5.2.2) |
|||
е — / с о |
|
|
|
|
|
|
|
что равносильно |
условию (5.1.7). |
|
|
|
4 |
В. И , К р ы л о в , Н . С . С к о б л я |
98 ОБРАЩЕНИЕ ПРИ ПОМОЩИ КВАДРАТУРНЫХ ФОРМУЛ [ГЛ, 8
Для доказательства этого в интеграле (5.2.2) заменим ( у ) его выражением (5.2.1) и выполним интегриро
вание по частям:
е + iсо
S epp-sP f { j ) p - mdp =
е — гео
8 + гоо
Pn- ffl-1 i № |
- " |
- s+V p = |
|
||||
г — гео |
|
dpn |
|
|
|
|
|
|
|
|
|
е + |
г*оо |
|
|
= (— 1)'у г- т -1 |
^ |
№ |
П-*+1) |
|
|||
8 — г*со |
|
||||||
|
dp- |
|
|
|
|
||
|
е 4* г оо |
|
|
|
|
|
|
■(— \)п (п — т — 1) |
^ |
г ,п - т - 2 . |
^ |
( е Рр-П-ЗД)Ф . |
|||
|
|
|
|
dp- |
|
|
|
Легко показать, что первый |
член |
полученного |
выра |
||||
жения равен нулю, так |
как |
каждое его слагаемое, |
полу |
ченное после дифференцирования произведения epp~n~s+1,
будет |
иметь вид epp~k (k ^ s ) |
и ввиду того, что функция |
||||||
ер остается |
ограниченной |
на линии интегрирования, будет |
||||||
стремиться |
к |
нулю |
при |
удалении точки р в бесконечность |
||||
вдоль |
прямой |
интегрирования. |
по частям п — т — 1 |
|
||||
Выполнив |
интегрирование |
раз, |
||||||
для интеграла (5.2.2) получим |
следующее выражение: |
|||||||
|
|
|
8 + гео |
|
|
|
|
|
± { п — т — 1)! |
^ |
^ i ( e pp~n's+1)dp== |
|
|||||
|
|
|
8 — /С О |
|
|
|
т - 1)! ^dm (epp-"-s+1) 84“ /со |
|
|
|
|
|
= |
± |
(п - |
||
|
|
|
|
|
|
|
8 |
— /С О |
А так |
как |
s > О, п ^ |
1, |
то полученное выражение |
тоже |
|||
будет |
равно |
нулю. |
Таким образом, многочлен Рns)( |
^ |
заданный равенством (5.2.1), удовлетворяет условию орто гональности (5.1.7), а ввиду единственности многочлена, удовлетворяющего этому условию, можно заключить, что
многочлен Р(/-1 |
будет отличаться от (oj,s) |
только |
постоянным множителем, равным старшему коэффициенту
n U ,
Pln ( I ) = («■+ s - 1) (п + S)... (2л + s - 2) со<5)! 1
§ 5.2] |
ОРТОГОНАЛЬНЫЕ МНОГОЧЛЕНЫ |
99 |
|||
5.2.2. |
Рекуррентное |
соотношение |
для |
многочленов |
|
Р |
Для многочленов Рhs) |
можно получить рекур |
|||
рентную формулу, связывающую три многочлена P„Li |
|||||
P (ns)^yJ, |
подо®но Т0МУ> |
как эт0 Делается для |
|||
обычных |
ортогональных |
многочленов. |
Возьмем произведе |
ние — PnV -M ; оно является многочленом степени я 4-1
Р\ р Г
от 1/р и может быть представлено как линейная комби
нация многочленов Ро° (j^j, p iS){j)> •••> P"+l{ j ) :
i Pi"(j)=■2 |
c-‘Pi“(?)• |
<5-2-3> |
fc = |
0 |
|
Коэффициенты этого разложения можно определить по формуле
|
|
|
8+ to |
|
|
|
|
|
|
0Л= |
- ' Т + ,„ -------------------------------, |
(5.2.4) |
|||||
|
|
|
55 |
\ |
[ '’I" ( ? ) ] ' * |
|
||
|
|
|
|
е — /со |
|
|
|
|
из которой видно, что если |
£ < n — 1, |
то |
— |
|||||
многочлен степени ниже я и по условию |
(5.1.7) интеграл |
|||||||
в числителе |
(5.2.4) |
равен нулю. |
Таким |
образом, |
в соот |
|||
ношении |
(5.2.3) |
могут быть |
отличными |
от нуля |
только |
|||
Сп,п-1 »Сши Сп, п+1» |
|
|
|
|
|
|||
|
|
+ cmP? [ j) + cn,M P^+i ( 1 ) . |
(5.2.5) |
|||||
Переменную |
1/р |
обозначим х |
и соотношение (5.2.5) пере |
|||||
пишем в |
виде |
|
|
|
|
|
|
|
|
Р ^\Д х) =(anx + bn) P f (x) + cnPn.i (х). |
(5.2.6) |
||||||
Так как нам известно явное |
выражение Pjf (х), |
коэффи |
||||||
циенты |
а„, |
Ьп, |
сп |
определяются |
просто. Приравнивая |
4*