книги / Операционное исчисление и обобщенные ряды Лагерра
..pdf1 |
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 |
|
|
|