Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа справочная книга

.pdf
Скачиваний:
13
Добавлен:
19.10.2023
Размер:
7.61 Mб
Скачать

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 )

104

ОБРАЩЕНИЕ ПРИ ПОМОЩИ КВАДРАТУРНЫХ ФОРМУЛ

[ГЛ. 5

Оно

переведет

точку

z = р в точку

t = 0. В плоскости t

сделаем

разрез

вдоль

положительной

действительной

оси

от точки

t = 1

до бесконечности и будем рассматривать ту

ветвь

У 1 t, для

которой arg (1 — /) = 0 для действитель­

ных t <

1. Контур

интегрирования I перейдет в контур %,

охватывающий

точку

t = 0. Интеграл (5.2.11) после пре­

образования (5.2.12) примет вид

 

 

dn (

еР

\

 

 

 

 

 

dpn п+5"1) ~~

 

 

 

dt

 

 

 

 

~ ( Y \ ~ i + \ )

 

 

п\

 

4 У 1 t

 

 

 

 

 

 

 

 

м f P o r a + O - " " ' "

 

 

 

 

п\ ■

2n+s 1

2ni ( I ) X

(— 1)” n\

ш (I)2n+S 1 |X

|( /i - H - i )

dt

i V \—t {—t)™1

±(Y— t +i)

dt

(| / ^ i— ^ 4 - l ) S 2 | ^ 1 — ^ /"+1 ’

(5.2.13)

Подставляя теперь выражение (5.2.13) вместо п-й произ­ водной в явное выражение (5.2 . 1) для многочленов Рjf} j,

получим

£(УГ=7+0

t = { — \ ) n e- P p n ^ - i {-\)пп\ ______________ dt_

4Jli{ L jn+s~14 { (|ЛП^7 + i)i_2 V\Y~i tn+1

 

22га+5-2л

(/1=7- i)

dt

 

я!_ f

____ e_

 

 

pn2m

({\f\УТ=- t-^ \)s~2 \f\—t *пП’

или

1

25-2g2p_WT=i-\)

 

 

2?nn\

dt

( 5 . 2 . 1 4 )

(KHI7+ 1)S-2 yr= l tnn

x

§ 5.2]

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

105

Рассмотрим функцию

-- _ I (Y— t-1)

F (t) = __?S е____________ •

( Y ^ t + \ ) s-2V T ^ i’

она аналитична в точке t = 0; значит, в окрестности этой точки она может быть представлена своим рядом Тейлора

F (t)= J ] cntn,

«=о

где

 

 

 

 

 

 

 

с,2ш' 5

 

 

 

 

 

 

 

 

 

 

R(fT-~t-i)

dt

 

= J - C

2s~2e‘

 

 

P * [ j V

 

J { V \-t+ \)s-2V ^ - t tnvi

22«re!

Следовательно,

P (/!_ /_ !)

 

 

 

 

 

 

=

 

 

 

 

2s-*e2

 

 

 

 

 

( |/l - ^

+ l)s“2| / l - /

 

 

2гпп I

 

Введя новую

переменную

h — pt/4

и

положив 1—х,

последнюю формулу

можно записать в виде

 

 

J- ( V l — 4 <j c — l )

 

СО . . .

 

 

 

=

=

у

 

п\

(5.2.15)

i V \ - t t x + \У~г V 1—4Мх/jc

bd

 

 

 

 

 

/2 =

0

 

 

где переменная снова обозначена /.

Функция, стоящая слева в выражении (5.2.15), явля­ ется производящей функцией для многочленов Р& (х).

5.2.6. Распределение корней многочленов Р<>> (л;).

В конце предыдущего параграфа мы указали, что для завершения исследования возможности построения квад­ ратурной формулы (5.1.5), точной для многочленов сте­

пени 2п — 1 от

1/р, необходимо показать,

что корни мно­

гочленов аф0

или корни многочленов

(х), отлича­

ющихся от со^ (х) только постоянным множителем, лежат в правой полуплоскости при всех s > 0.

В этом пункте рассмотрим данный вопрос для некото­ рых частных значений s. Докажем следующую теорему.

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

Т е о р е м а 4. Все корни многочленов

= ( - D

при всех целых s ^ 2 лежат в правой полуплоскости,

т. е,

вещественные части всех корней положительны.

При

Д о к а з а т е л ь с т в о . Сначала возьмем s = 2.

доказательстве будем пользоваться некоторыми теоремами, известными из алгебры. Напомним *) их:

а) Для того чтобы вещественные части всех корней многочлена

Qn (х) = Ьпхп+ б*-!*"'1 -1- Ьп^ х п-2 + ... + Ь0

с вещественными коэффициентами были бы одного знака,

необходимо

и достаточно,

чтобы

корни многочленов

/ (х) =

Ьпхп-

Ь ^х"-2+

Ьп-4хл- 4- . . . ,

Ф (х) =

Ь^х»-1-

3xn~s+ V s *"-5 - . . .

были бы все вещественные и разделялись.

б) Если

все

корни

многочленов F (х) = А,/ {х) 4* рф (х)

вещественны

при

любых

вещественных А и р, то корни

многочленов f(x) и ф(х) вещественны и разделяются. Запишем многочлен P'i'(x) в виде

P i' (х) = апхп+ ап-ххп-х+ ап_2хп-2+ . . . + ахх + а0. (5.2.16)

Нам необходимо доказать, что вещественные части всех корней многочлена (5.2.16) одного знака. Для этого, на основании только что сформулированных теорем а) и б), достаточно установить справедливость следующих фактов:

1) все корни многочленов

Qn (х) = апхп - ап^хп-2+ a„_4x"-4- . . . ,

Rn-i (х) = — an_3xn~3+ ап^ьхп~ъ—...

вещественны и разделяются; 2) все корни многочленов

Р$’|г) (х) = Аапхп4- р а ^ л :”- 1 — Aa„-2*"-2

-рал_3хл- 34- karl-ixn-i 4- 1Шп_ъхп~5

будут вещественны при любых вещественных А и р.

*) См.

Д. К. Ф а д д е е в , И.

С. С о м и н с к и й, Сборник задач

по высшей

алгебре, М., «Наука»,

1968, стр, 106, 108.

§ 5.2]

 

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

107

Из (5.2.7) видно, что

многочлен Р'п (лг) удовлетворяет

рекуррентному соотношению

 

 

 

 

РТ (х) = 2 (2л - 1) х Р

(х) +

Рп-ъ (х),

(5.2.17)

причем Р0— 1,

Р1(х) = 2х— 1.

 

 

 

 

Теперь

найдем

 

рекуррентное соотношение для много­

членов Р<£’ (х). Запишем многочлен Р'п (х) в виде

где

Рп (х) — А„ (х) -j- Вп (х) + Сп {x)-\-Dn (х),

 

 

 

 

 

 

 

 

 

 

 

Ап(х) =

апхп+ are_4x"-4+

a ^ 8xn~s+ . . . ,

 

Вп (х) =

a,,-!*"-1 +

ап-ъхп~ъ+

а„_9хи- 9+ . . . ,

(5.2.18)

Сп {х) = ап-2хп~2+

ап-ьхп~ь+

а„_wxn~10+ . . . ,

 

Dn (х) = а„_sxn-s+ ап-7хп~7+ а ^ х '1-11 +■■■ .

 

Тогда многочлен

Р^> и) (*) можно

записать следующим

образом:

 

 

 

 

 

 

 

 

 

 

р (к Ю { х ) = Я [4 (х) - Сп ( х ) ] + р [Вп ( х ) - Dn (* )] .

Рекуррентное соотношение (5.2.17) примет вид

 

А п ( х ) + Вп ( х ) + Сп ( х ) 4- Dn(х) =

 

 

 

= 2 (2л — 1) х

[ А п-1 ( х ) -j- Вп- 1 (х) -f- Cn-i ( х ) 4- Dn~i (х)] -j-

 

 

 

 

4“ Ап-%(х)+

Вп-2(х) + Сп- 2 (х) +

Dп—2 (X).

Отсюда

и

из

равенств (5.2.18)

и подобных равенств для

P'n-i (х)

и Рп’ -ъ (х)

получаем

 

 

 

 

 

 

Ап (х) =

2 (2л -

1) x4„_j (х) 4 - Сп- 2 (лг),

 

 

 

(х) = 2 ( 2 л - \)х Вп-х (х) 4- Dn (х),

 

 

 

С„ (х) =

2 (2л -

1) хСя_х (х) 4- 4 „_2 (х),

 

 

 

Dn(х) =

 

2 (2л - 1) xD„_x (х) 4 я_2(х).

 

Умножим

первое

и третье равенства соответственно на Я

и — Я, а

второе и четвертое на

р и — р и сложим их:

Я [Л„ (х) -

С„ (х)] 4-р [Я„ (х) - Dn(х)] =

= 2 (2л -

1) хЯ [Ап-! (х) - Сп-! М ] + Я [С„_2 (х) - А„-г (х)] 4

4-2 (2л - 1) хр [Вп-! (X) - Dn-! (*)] + Р [Dn- 2 (х) - Вп - 2 (х)].

108

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

Полученному

равенству придадим вид

 

 

 

Я [Ап(х) - Сп(а)] + р [Вп {х) - Dn(а)] =

 

 

 

=

2(2я — \)х{%[Ап-1 (х) - Сл_а(х)] + р

 

(х) -

(л:)]}—

 

 

-

{Я [Ля_2 (х) Сл—^(а)] +

р [В„_2(л:) - ZV, (а)]}.

Последнее

равенство

есть

не что иное,

как

рекуррентное

соотношение для Р&>м (х), а именно:

 

 

 

 

P f. в) (х) =

2 (2п -

1) хР ^ >

(х) -

P*Lri (х),

(5.2.19)

где

 

 

 

 

Р? "|Х>М = 2Ях — р

 

Я ф. О,

 

Ptf- ^ W = Я,

 

при

 

В<°. и) (х) =

0,

 

Р<°'В)(х) =

— р

 

при

Я =

0.

Случай р = О,

Я =

0 тривиален, и его

из

рассмотрения

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

многочленов

 

 

 

 

 

 

 

 

 

 

 

 

Р(К ю (д-);

p<Mj) (д-),

. . . ,

р(К я) (х)

при Я ф 0

(5.2.20)

или

^

 

(

*

) .

 

Р

^

Д

?

w

> 0,

5. 2( . •21)• • >

 

Р фп

 

 

причем во второй последовательности Р®>»*>(х) будет мно­

гочленом

степени

п 1 .

 

 

 

 

 

 

 

 

Эта последовательность многочленов обладает следую­

щими свойствами:

 

 

 

 

 

 

 

 

 

1)Последний многочлен последовательности есть отлич­ ная от нуля постоянная, а именно, в первой последова­ тельности это есть Я, во второй р.

2)Ни при каком значении х два рядом стоящих мно­ гочлена последовательности не обращаются в нуль.

Всамом деле, пусть хг является корнем Р®* -и>(х) и

PP-L-f (*); тогда ввиду (5.2.19) хх будет корнем и много­ члена (х). Рассуждая так дальше, мы пришли бы в конце концов к тому, что этот общий корень был бы корнем и многочлена ^ (х) в последовательности (5.2 .20) или Р<*-■^ (х) в последовательности (5.2.21), что невоз­ можно, так как эти многочлены есть отличные от нуля постоянные.

3) Если какой-нибудь из многочленов последователь­ ности обращается в нуль при вещественном значении х, то два соседних с ним многочлена имеют при этом х зна­ чения разных знаков. Действительно, из формулы (5.2.19)

§ 5.2]

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

109

видно,

что если Р^'_^ (х) = 0, то

 

 

/М. Р) (X) = — Рп'-У^ (X).

 

Таким образом, эта последовательность многочленов

образует обобщенный ряд Штурма. Из алгебры

известно,

что если обозначить через и (х) число перемен знака в ряду

этих многочленов при заданном значении

х, а через г

число вещественных корней многочлена

^ (х), то будем

иметь г ^ и (— со) —и (со). А так как

последовательность

р(К р) (х)

содержит многочлены

всех

степеней с коэффи­

циентами

при старших членах одного

и того же знака, то

и (со) —0,

а

и(— со) = ц в

последовательности (5.2.20) и

и (— оо) =

п — 1 в последовательности (5.2.21). Следователь­

но, все корни многочлена

^ (х) не только вещественны

при любых

вещественных к и ц,

но,

кроме того, просты

и корни

 

(х) взаимно разделяются с корнями *) Р&. р ) (*).

Отсюда

на основании

теорем

а)

и б)

(см. стр. 106)

мы можем сделать вывод, что вещественные части всех

корней

многочлена

Р„’ (х) одного знака. Этот знак может

быть

только

плюсом,

так

как

коэффициент

при

х’1"1

в Рп'(х) отрицательный. Таким

образом, доказано, что

все

корни многочлена

Р ^ (х) при s = 2 лежат в правой

полуплоскости.

 

 

 

 

 

 

 

Перейдем теперь

к

параметру

s = 3. Из явного выра­

жения

(5.2.1)

для

многочленов Р& (х) видно,

что много­

член

Р’п (х) можно

получить

из

многочлена Р'п(х),

если

коэффициенты ak последнего умножить на ft I k I 1 ' Пусть

для

s =

2 многочлен РТ (х)

имеет вид (5.2.16), тогда для

s —3

многочлен Рп'(х) имеет следующий вид:

 

р п' М =

 

апх" + jjqpj- ап^х»-' +

 

 

 

 

2 я - 1

 

n-j-2

п + 1

 

 

+

я +

1 а11-гхп~2-f ...

я + 1 ахх

я -f-1

(5.2.22)

Для дальнейшего нам потребуется следующая

известная

из анализа теорема**).

 

 

 

 

*) См. Ф.

Г а н т м а х е р ,

Теория матриц, М., «Наука», 1967,

стр. 471.

 

П о л н а , Г.

С е г е , Задачи

и теоремы из анализа,

**) См. Г.

т. II,

М.,

Гостехиздат, 1956,

стр.

77.

 

 

Соседние файлы в папке книги из ГПНТБ