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

книги / Операционное исчисление и обобщенные ряды Лагерра

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
9.58 Mб
Скачать
ckn *

1

2

 

 

3

 

 

 

 

Узлы:

п

 

 

 

у(*)= 2 >

созя*х

 

 

2 k

 

 

 

 

 

tT -c

Xk = 2n+ 1

 

в

_

 

Tl^ k ^ H

C = 2À T

2 - ^ COS2::2ÏÏTÎ

 

f ( x ) — ве­

 

та—о

 

 

щественная

£**■“ <* + С*. 2Л+3

 

 

четная

 

функция

*r2s

tr*

t r zc

t r 3s

Узлы:

Xk ~ -

2л+1

/(JC)~ ве­

щественная

нечетная

функция

Узлы:

—n + l< fe< »

/(х)—ве­ щественная функция

Узлы:

/ ^ - в е ­ щественная четная функция

Узлы:

/ (^ -в е ­ щественная нечетная функция

п

У(* )= 2 dkna{n~kx

 

k ~ i

П

 

ь

 

 

4

^ -о

 

 

V 1

 

 

d i n * *

2 в + Г ^

^"*8Ш * 2в+1

 

 

 

/Я—1

 

 

 

< г* л = < * А + Ч 2 л + 1

 

 

у ( * ) - 2 "

C A n C O S B A X + ^ A B S i n - A *

 

 

а Г о

 

ft= 1

 

 

с». -

-J 2 ' (У'го+/_,п) C0S,tm^ + ^

/n

 

771-0

 

 

 

 

1

 

 

Hîfe

 

^ =

F 2

(/r,n" /- m)siI1~

 

« в = c * + ® M » 5 ‘W ^ f e + ' ^ a n

 

Ся»=с»+(с»в+С*п+ * *

*

 

 

п

 

 

 

 

у ( х ) - Х

CknC03T'hX

 

 

 

fe~0

 

 

 

 

 

2 жк" „

» *

. 4-

 

 

m^Ofm С09Л

Cftn~ Ck+G*-

jjpji

/п= 0

задача *r3s совпадает с задачей t r i$ ..

При /д^О нечетность функции /(*) He

обеспе­

чивает обращения в нуль коэффициентов

условиях задачи tr$ положить у(ха) = Y"(f(xn)-\-f(x—„)), то за­

дача tr3 совпадает с задачей tr'\, т. е. задача tr\ оказывает­ ся алгоритмически эквивалентной задаче trs.

В задачах аппроксимации функций алгебраическими многочленами большую популярность завоевали многочле­

ны П. Л. Чебышева.

созпагссозж ( n = 0,

1, 2 ,...)

образуют

Многочлены Тп (х)=

полную ортогональную

систему с весом

Vl-X-

В

+ 1 ]. Любая функция ограниченной вариации на

[— 1, + 1]

разлагается в

сходящийся ряд по многочленам Чебышева

f(x) =

^ c sT s(x),

dx.

(5.1.8)

 

S=0

У L-X-

 

Для практических целей большую ценность приобрела задача интерполирования по многочленам Чебышева. Су­ ществует два варианта задач указанного типа.

Задача Thu [33]. Построить многочлен

*)

2 ^Ол “Ь с М

х ) -f- с^пТ^х) “Ь •■•Ч” сппТ п(дг),

 

 

 

 

 

 

(5.1.9)

такой, что L Tn+ [f; xk) = f

(xk),

где в качестве узлов интер­

полирования избраны нули многочлена

Тр+1 (д:):

 

Хк =

cos

 

0 ^ k ^

п'

Опираясь на отношение ортогональности

 

 

 

 

'О, если

s

Л+1 2

^ в

T m ( Х к)

1,

если

s = m Ф О (5.1.10)

ft=0

 

 

 

2,

если

s — m ~ 0 ,

получаем ce„ =

_ 5 _ V

 

(0 < s < r a ) .

n+ i j j fkT s(xk)

fe=0

Продолжим по периоду соотношение (10). С этой целью за­

метим, что Та{хъ) = cos . Пусть е = 2г(п+ 1 )+ ц, где

—в + 1 < ц < в . тогда Т в(хк)= (—1)ГТ [v.\(xk).

Поэтому

 

2,

если s =

т = О

 

 

л-1-1 2

m(Xk) (—1)г, если s =

2(га +

1)г ±

т

*-о

О,

если s Ф 2(ге +

1)г ±

и».

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

GO

Задача Thi2. Построить многочлен

-£<Т„+1 (Л * ) =

с 0л +

х(х ) -j- СгпУзС*) + • • • + “! " СппТп (*)»

 

/(*а), где в качестве узлов интер-

такой, что i r n+; (Л * а) =

Заменой A* = CO S© задача T h u сводится к задаче tr\c (мо­ дификация К. Ланцоша).

Замечания 1. В отличие от задачи Thn в задаче Thn в состав узлов интерполирования включаются концевые точ­

ки интервала интерполирования.

2. Сопоставим зависимость коэффициентов с sn от с г для задач Th,„ Thl2.

Задача ТЬц

Спп. — Сп 2(C2fü+n — CdfT,4-11Н~ Csfn+1) • • • )

(5.1.11)

Сп—1,п = * Св—1 — (Сп+3 4 " Сзп+1 — СЗ»+5 — Сбв+З Н“ • • • )

С п п = С „ — (Сл+2 Hr СзЫ-2 С Зп+4 “ С S n + i + • • • )

Задача Th12

(5.1.12)

Из сопоставления приведенных таблиц следует: достоинство Th, г-метода состоит в том, что величина -g-Cnn прибли-

жает коэффициент с п с точностью до порядка малости коэф­ фициента с зл ; однако в случае, когда коэффициенты cs ряда (8) убывают достаточно быстро, то предпочтение сле­ дует отдать Гйц-методу, так как в целом порядок прибли­ жения коэффициентов к с* в этом случае будет лучше.

Многочлены Чебышева (второго рода)

sin (tt+ l)a rc cosy

(n = 0, 1 , 2, . . . )

удобны для полиномиальной аппроксимации функций, при­ нимающих большие значения на концах интервала

(—1 » +'1)*

Аналогично задачам Thu, Th\z в зависимости от узлов интерполирования могут быть сформулированы два вари­ анта задач. Из них отметим аналог второго варианта.

Задача Th22Построить многочлен

L TT ( f , х) = dlnu x(x) + d2aU2(х) + . . . + d„ -1, nUn-i(x), (5.1.13)

ЛЛ

такой, что Lrr (f , Хъ) = f(x i), где в качестве узлов интерпо-

U П

k

лирования избраны точки xk = cosit — (O ^ k ^ n ).

Замена х = cos0 сводит задачу Тйгг к задаче tns отно­ сительно функции g(©)=sin@/(cos0 ).

Таким образом, искомые коэффициенты d sn исчисляют­ ся по формуле

п— 1

s

2

^

dsn = ~

2л s(xji)$innk — ( l < s < / i - l ) ,

 

fe=i

 

где g(xk) = sin

f(xh).

 

Несколько замечаний о порядке вычисления коэффици­ ентов интерполяционных многочленов.

Если в памяти машины задать матрицу синусов и коси­ нусов, отвечающих интерполяционной задаче с iV-узлами, а вычисление коэффициентов интерпретировать как умно­ жение матрицы на вектор, то потребуется порядка N2 умно­ жений и N2 алгебраических сложений, т. е. суммарно потре­ буется порядка 2N2 арифметических действий. Поэтому при больших значениях N объем вычислений существенно воз­ растает и вопрос об экономии машинного времени становит­ ся актуальным.

Некоторую экономию машинного времени порядка

6 j^-y-j арифметических действий дает рекуррентная схема

вычисления искомых коэффициентов, описанная в работе [108].

Дадим полезную для вычислений интерпретацию этой схемы. С этой целью нормируем интервал интерполирова­ ния [— 1, + 1] к отрезку [0, 1]. Тогда искомые коэффици­ енты принимают вид

Введем в рассмотрение многочлен QN—i(x )= f0 - { - fix -\-fzX2-\-

+•. . + f N-iX1** 1,тогда, очевидно, N -ак}г = QN-i(e " w ). Таким образом, вычисление дискретного коэффициента

Фурье a kN эквивалентно вычислению значения многочлена

 

k

 

k

Этот

Q N - I (X ) в комплексной точке x=cos2jT-^p--Msin2n

 

процесс удобно реализовать по схеме Горнера,

если /г —

комплексные величины или посредством

схемы деления

многочлена QN- г(*) на х2—2cos22n* k +■1

(см.

например,

[53]), если fi — действительные числа.

Описанный метод удобен и тем, что он использует лишь значения cos2it k sin2n k которые в свою очередь могут быть рассчитаны по известным тригонометрическим рекур-

рентным формулам, зная лишь два значения cos

2“

. 2^

 

sin-^.

В силу рекуррентности указанный метод может быть представлен короткой программой. Накопление погрешно­ стей метода незначительно.

Известно несколько приемов минимизации числа умно­ жений [70, 94, 108]. Существенно сократить число арифме­ тических действий (до порядка ZNlgzN) позволяет метод работы [134]. Схему деления многочлена на двучлен удоб­ но использовать и для вычисления коэффициентов с в за­ дачах Thu, Th\z. Это следует из того, что

T r(xk) = cos

= Гг*^ Т(пТ1)~) (задача Tkn)’

Тт(хк) = cos — — Th(xr) (задача Thn).

§ 2. Схема восстановления оригиналов посредством разложения их в ряды но функциям Лагерра

Описываемая ниже схема была использована в предыду­ щих главах книги для построения общей теории операцион­ ного исчисления. В силу большой практической ценности этой схемы изложим ее несколько подробнее на р-языке.

Пусть задано операционное соотношение

 

f(t)+ F (p) (Rep>v>Vo),

 

 

где vo — абцисса абсолютной

сходимости

интеграла Лап­

ласа.

 

 

 

 

 

 

Будем искать разложение f(t) в форме ряда

 

 

 

СО

 

 

 

 

 

т ~ е\ * 2

a ^ t l h ) ,

 

(5.2.1)

 

 

Й=0

 

 

 

где q>n(t) — функции Лагерра,

определенные посредством

многочленов Лагерра £„(£) по формуле

 

 

 

 

=

 

 

 

(5.2.2)

Воспользуемся

тем, что

изображение Лапласа

функции

?а У) имеет вид

 

(4

 

 

 

<?n(f) 4" T 7------уГЛ>и тем, чторазложение (1)

 

 

( т + » )

 

 

 

эквивалентно разложению

 

 

 

 

 

е -2Т* fQ h t) ~

2

ah%(2t).

 

(5.2.3)

 

 

 

ft=0

 

 

 

Тогда в пространстве изображений формальному разло­

жению (3) отвечает разложение

 

 

 

 

 

 

 

00

(1- р ) к

 

 

 

( 2йР

 

2

 

(5.2.4)

 

=

а* (1+р)*+1’

 

 

 

W-0

 

 

 

Осуществим конформное отображение

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

Вер^О в круг 121^ 1 с помощью функции

 

 

(5.2.5)

Тогда ряд (4) примет вид

со

 

 

S(z) = 2

«***.

(5-2.6)

^(г) — Л(1+г)

+ 2А 1+г)’

(5.2.7)

‘Следовательно, задача вычисления коэффициентов а% раз­ ложения (1) сводится к вычислению коэффициентов сте­ пенного ряда (6) функции g(z). В общем случае разложение

(1) — это обобщенный ряд Лагерра в смысле определения главы 2. Для практики численного обращения представля­ ет интерес, когда обобщенный ряд Лагерра аппроксимиру­ ет реальную (а не обобщенную) функцию.

Достаточно широкий класс таких функций f(t) описыва­ ется ограничением

e - i h* f(h t)e L 2(0, оо).

(5.2.8)

Поскольку функции Лагерра (2) образуют

в пространстве

1 /2(0, оо) полную ортонормированную систему, то коэффи­ циенты разложения (1) в этом случае должны удовлетворять условию

2 | a * l 2= f e - W \ f( h t ) \ 2d t < o o .

 

s-о

g

 

 

 

Говорят [106], что функция g(z)= 2

a kzk ( |z I <1)

 

 

 

ft=o

 

 

 

 

00

 

 

принадлежит классу Я 2, если

2 I

I 2< оо.

 

 

 

О

 

 

Следовательно, если оригинал f(t) таков, что выполнено

условие (8), то

 

 

 

 

*<3>= д

а

+

г й ) е я 2,

(5.2.9)

где F(p) изображение Лапласа функции f(t). Отметим не­ которые свойства функций класса Н2 [106].

Если g{z) 6 # 2, то

1 ) ее тейлоровское разложение 2

в силу ортогонально­

г о

 

сти системы функции 1, в, г2, . . . на окружности С(|г|=1 ) сходится в среднем к некоторой функции gi(z) класса Ь% на С ;

2) g](z) является граничным значением функции g(z) в следующем смысле: функция g(z) при стремлении г к С п о любому пути в углах с вершинами на С и со сторонами,, являющимися хордами, принимает почти всюду значения gi(z);

3) ак = 2Гг I

g(z) ~ î k ’>

J

g{z)zndz = 0

(га = 0, 1, 2, . . . ) ,

c

z

с

 

 

 

 

I g ( ^ ) | 2d9 = 2 ^ 2 И * I 2,

0

 

 

ft-0

 

 

 

 

g(z), если I

г I

< 1

c

 

0,

если I 2 I

>

1 .

Вывод: если изображение Лапласа F(p) удовлетворяет требованию (9), то отвечающий ему оригинал f(t) однознач­ но восстанавливается, в общем случае в метрике Lz(0, оо), рядом

оо

/(At) = etht 2 а*?* (0 ,

(5.2.10)

й=0

 

где коэффициенты ак могут быть вычислены по формуле

IT S * (*) I*-» * = 0, 1 , 2 , . . .

(5.2.11).

либо по формуле

аъ=

е-Ч9

(5.2.12)

Это сравнительно простой метод определения лагерровского спектра функций, в частности, может быть с успехом ис­ пользован для вывода большого количества новых разложе­ ний специальных функций по многочленам Лагерра.

П р и м е р

1. Разложить функцию

„(ЗУв ряд по

обобщенным многочленам Лагерра LW (?)

[21, 93].

Воспользуемся операционным соответствием

 

2fi—v

 

i ' - ' W

» J ® (8 V m w - j à z g * ' { ? ? )

(Rev > — 1, Re|A > — 1).

Применяя теорему подобия это соотношение приведем le виду:

00

Образуем функцию

 

«Г(г) =

= (Х г)~ ^ (2]/ Б ‘),

J {n\z) I tml =

( - 1) ^ J„ +,(2l/Xi) |г=1.

Так как

 

(5.2.13)

то

M*V

 

^ 'n Z rT 1

( Х » Г ^ , ( З ^ Т ) = Х - ^

 

п=-0

 

П р и м е р

2. Вычислить свертку (f —

 

 

О

где «î»„ (*) = 2

M O ,

 

А»0

 

 

(здесь Р„ (ft) — произвольный многочленстепени п от цело­ численного аргумента ft).

Найдем изображение заданной функции

СО

t

е~ир 2

Pn(k)

J е~Р* J

(t — х) t'P j,(2К t ) d t d t =

_»+ 2

„ А!

 

 

ОО

 

 

=

F (P ) .

 

 

 

Так как

 

 

 

 

 

 

 

2- , - ?_р(_1_)=е- г2

is ®

(2 -

1)\

а л = - Ç

Д*Р„(0),

 

 

ft=0

 

 

 

 

 

 

 

Д*Р„ (0) = 0

при k

>ге,

 

 

то в силу (13):

 

 

 

 

 

 

t

 

 

 

 

 

 

 

f ♦ , <f -

ф -12J .(2^

) *

=

4 l

гц‘У + з)

щ * 11(*>•

-

 

 

 

А=0

 

 

П р и м ер 3. Найти сумму ряда

2 *r(*+v+i)

(0 =

Фп>V(t, х),

 

 

 

 

4=0

 

 

где P n (ft)— произвольный многочлен степени

п целочисленного аргумента й. Имеем

 

 

Г e-pt г ф п>, (f, х) d t =

- é p

2

( т

—1)

=

О

 

^

4=о

4

 

 

== ~^+i exp .К1 * )] 2

(“ 1)‘ -р^ ) 1 г ( т '_ 1)

ехр[К^ ~

- 1 ) ] = ^ e x p [ K

l - i ) ] 2 ( - i ™ ^

( - ± - i f =

= 1 - i2-0

 

 

D*+/+r

 

 

 

Следовательно, Фл, v(t, #) =

(x ty -^ e*^ (—1)*

qj

(я)

(я*У/2Х

 

 

 

i-o

 

 

 

Соседние файлы в папке книги