Научная статья на тему 'Полиномы, ортогональные по Соболеву, порожденные многочленами Мейкснера'

Полиномы, ортогональные по Соболеву, порожденные многочленами Мейкснера Текст научной статьи по специальности «Математика»

CC BY
172
53
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПОЛИНОМЫ / ОРТОГОНАЛЬНЫЕ ПО СОБОЛЕВУ / ПОЛИНОМЫ МЕЙКСНЕРА / ОРТОГОНАЛЬНЫЕ НА СЕТКЕ / ПРИБЛИЖЕНИЕ ДИСКРЕТНЫХ ФУНКЦИЙ / СМЕШАННЫЕ РЯДЫ ПО ПОЛИНОМАМ МЕЙКСНЕРА / ОРТОГОНАЛЬНЫМ НА РАВНОМЕРНОЙ СЕТКЕ / ORTHOGONAL SOBOLEV POLYNOMIAL / MEIXNER POLYNOMIALS ORTHOGONAL ON THE GRID / APPROXIMATION OF DISCRETE FUNCTIONS / MIXED SERIES IN MEIXNER POLINOMIALS ORTHOGONAL ON A UNIFORM GRID

Аннотация научной статьи по математике, автор научной работы — Шарапудинов И.И., Гаджиева З.Д.

Рассмотрена задача о конструировании полиномов mα(x, q), (n = 0, 1,...), ортогональных по Соболе-ву и порожденных классическими полиномами Мейкснера. Эти полиномы могут быть определены с помо-[k]r,kщью следующих равенств mαk!(x, q) = x, x[k] = x(x 1) · · · (x k + 1), k = 0, 1,..., r 1,x-rmα 1 ),[r-1] α αr,k+r (x, q) =(r-1)!t=0(x 1 t)mk (t, q), где через mk (t, q) обозначены полиномы Мейкснера, орто-нормированные на сетке Ω = {0, 1,...} с весом ρ(x) = qx Γ(x+α+1) (1 q)α+1. Полиномы mα(x, q) (n = 0, 1,...)Γ(x+1)r,nобразуют ортонормированную систему на Ω = {0, 1,...} относительно скалярного произведения типа Соболева следую-щего вида:(mα r,n, mα r,mr-1r,n) = ∆kmαk=0r,m(0, q)∆kmα∞r,n(0, q) + ∆rmαj=0r,m(j, q)∆rmα(j, q)ρ(j).r,nДля mαn(x, q) мы получили явную формулу, содержащую полиномы Мейкснера Mα-r (x, q):I r-1[ν] lmαr,k+r (x, q) = ( q )rq 1k{hα(q)}-1/2k+rMα-r (x, q) Ar,k,ν x ν!ν=0, k = 0, 1,...,где A =( q 1r,k,ν ν Γ(k + α + 1), Mα(x, q) =Γ(n + α + 1)n xn [k] [k] ( 1 \k1 -,hαn (q) =q n + α -nqn(k + r ν)!Γ(ν r + α + 1) nΓ(α + 1).n!k=0Γ(k + α + 1)k! q

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Sobolev Orthogonal Polynomials Generated by Meixner Polynomials

The problem of constructing Sobolev orthogonal polynomials mα(x, q) (n = 0, 1,...), generated by classical Meixner’s[k]r,kpolynomials is considered. They can by defined using the following equalities mαk!(x, q) = x, x[k] = x(x-1) · · · (x-k +1),k = 0, 1,..., r 1, mα(x, q) = 1x-r), (x 1 t)[r-1]mα(t, q), where mα(t, q) denote Meixner’s polynomialr,k+r(r-1)!k kt=0of degree k, orthonormal on Ω = {0, 1,...} with weight ρ(x) = qx Γ(x+α+1) (1 q)α+1. Polynomials mα(x, q),Γ(x+1)(n = 0, 1,...) are orthonormal on Ω = {0, 1,...} with respect to the inner productr,n(mα r,n, mα r,mr-1r,n) = ∆kmαk=0r,m(0, q)∆kmα∞r,n(0, q) + ∆rmαj=0r,m(j, q)∆rmα(j, q)ρ(j).r,nFor mαn(x, q) we obtain the explicit formula that contains the Мeixner polynomial Mα-r (x, q):I r-1[ν] lmαr,k+r (x, q) = ( q )rq 1k{hα(q)}-1/2k+rMα-r (x, q) Ar,k,ν x ν!ν=0, k = 0, 1,...,r,k,νwhere A = ( q 1 ν Γ(k + α + 1), Mα(x, q) =Γ(n + α + 1)n xn [k] [k] ( 1 \k1 -,q n + αhα -n(k + r ν)!Γ(ν r + α + 1) nn!k=0Γ(k + α + 1)k! qn (q) =q Γ(α + 1).n

Текст научной работы на тему «Полиномы, ортогональные по Соболеву, порожденные многочленами Мейкснера»

ry and their applications : Proc. 18th Intern. Sarat. Winter School]. Saratov, OOO Izd-vo "Nauchnaia kniga", 2016, pp. 322-326. 12. Zygmund A. Trigonometric series: Vol. I, II. Second edition, reprinted with corrections and some additions. London ; New York, Cambridge Univ. Press, 1968, vol. I : 383 p.; vol. II : 364 p. (two

volumes bound as one) (Russ. ed. : Zygmund A. Trigonometricheskie riady. Tom 1. Moscow, Mir, 1965, 616 p.).

13. Korneichuk N. P. Ekstremal'nye zadachi teorii pri-blizheniya [Extremal problems of approximation theory]. Moscow, Nauka, 1976, 320 p. (in Russian).

Please cite this article in press as:

Shakirov I. A. On a Limit Value of a Remainder of the Lagrange Constant Corresponding to the Lagrange Trigonometrical Polynomial. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2016, vol. 16, iss. 3, pp. 302-310 (in Russian). DOI: 10.18500/1816-9791-2016-16-3-302-310.

УДК 517.587

ПОЛИНОМЫ, ОРТОГОНАЛЬНЫЕ ПО СОБОЛЕВУ, ПОРОЖДЕННЫЕ МНОГОЧЛЕНАМИ МЕЙКСНЕРА

И. И. Шарапудинов1, З. Д. Гаджиева2

1 Шарапудинов Идрис Идрисович, заведующий отделом математики и информатики, Дагестанский научный центр РАН; Владикавказский научный центр РАН; Дагестанский государственный педагогический университет, Махачкала, sharapud@mail.ru

2Гаджиева Зульфия Джамалдиновна, научный сотрудник отдела математики и информатики, Дагестанский научный центр РАН, Дагестанский государственный педагогический университет, Махачкала, zulfiya.gadzhieva.1976@mail.ru

Рассмотрена задача о конструировании полиномов ш"п(ж,д), (п = 0,1,...), ортогональных по Соболеву и порожденных классическими полиномами Мейкснера. Эти полиномы могут быть определены с помо-

щью следующих равенств m'k(x,q)

„И

x(x — 1) ••• (x — k + 1), k

0,1,...,r — 1,

х — г

ш"к+г (ж,д) = (г-1)! I] (ж - 1 - ¿)[г-1]где через (£,<?) обозначены полиномы Мейкснера, орто-

нормированные на сетке П = {0,1,...} с весом р(ж) = дх г(гх(+"+)1) (1 - д)а+1. Полиномы ш"п (ж, д) (п = 0,1,...) образуют ортонормированную систему на П = {0,1,...} относительно скалярного произведения типа Соболева следующего вида:

(m'n, mam} = Y^ Akman(0,q)Akm'm(0,q) + Arm'n(j,q)Arm'm(j,q)p(j). k=0 j=0

Для ma, n(x, q) мы получили явную формулу, содержащую полиномы Мейкснера M'-r (x, q):

где Ar

m'k+r (x, q) ' q — 1 \ v

^Г {h'(q)}-1/2

q—1 Г(к + a + 1)

q / (k + r — v )!r(v — r + a + 1)

M'+rr (x,q) — X]

, M' (x, q)

r-1 A i. x[v]

v !

v=0

r(n + a + 1)

n!

, k = 0,1

1 „jfci„ E

r(k + a + 1)k! V q

hn (q) = Г + )q-n r(a + 1).

1—

Ключевые слова: полиномы, ортогональные по Соболеву, полиномы Мейкснера, ортогональные на сетке, приближение дискретных функций, смешанные ряды по полиномам Мейкснера, ортогональным на равномерной сетке.

DOI: 10.18500/1816-9791 -2016-16-3-310-321

k

1

ВВЕДЕНИЕ

Теория полиномов, ортогональных относительно скалярных произведений типа Соболева, получила в последние три десятилетия интенсивное развитие и нашла ряд важных приложений (см. [1-6] и цитированную там литературу). Характерной особенностью скалярных произведений типа Соболева

является, в частности, то что они, как правило, содержат слагаемые, которые «контролируют» поведение соответствующих ортогональных полиномов в одной или нескольких точках числовой оси. Например, часто рассматривают скалярное произведение вида

г 1 р Ь

</,д> = Еf(а)+ / /(г)С0д(г)(1)

V=0 ^а

в котором / и д — функции, заданные на [а, Ь] и непрерывно дифференцируемые там г — 1-раз, для которых /(г-1)(х) и д(г-1)(х) абсолютно непрерывны и /(г)(х), д(г)(х) е Ьр(а, Ь), где Ьр(а,Ь) — пространство Лебега с весом р(х). Следует отметить, что полиномы, ортогональные по Соболеву, по своим свойствам могут весьма существенно отличаться от обычных ортогональных на интервале полиномов. Например, в некоторых случаях оказывается так, что полиномы, ортогональные по Соболеву на интервале (а, Ь), могут иметь нули, совпадающие с одним или с обоими концами этого интервала. Это обстоятельство имеет важное значение для некоторых приложений, в которых требуется, чтобы значения частичных сумм ряда Фурье функции /(х) по рассматриваемой системе ортогональных полиномов совпали в концах интервала (а, Ь) со значениями /(а) и /(Ь). Заметим, что обычные ортогональные с положительным на (а, Ь) весом полиномы этим важным свойством не обладают. Скалярное произведение (1) имеет одну особую точку, а именно точку а, в окрестности которой «контролируется» поведение соответствующих полиномов, ортогональных по Соболеву. Это

г —1

достигается за счет наличия в скалярном произведении (1) слагаемого вида ^ /^На)д(^(а).

V=0

С другой стороны, в работах [7-18] независимо были рассмотрены так называемые смешанные ряды по классическим ортонормированным системам {^к(х)}, которые представляют собой ряды по системам функций {^г,к(х)}, порожденным соответствующей системой {^к(х)} посредством применения равенств вида

^(х) = (х —к,а) , к = 0,1,...,г — 1 (2)

1 Г

Рг,к(х) = ( _ / (х — ^)г—гк = г,г + 1,... . (3)

(г 1) ^ а

Обозначим через ^£2(аЬ) весовое пространство Соболева, состоящее из функций /(х), заданных

на [а, Ь] и непрерывно дифференцируемых там (г—1)-раз, для которых /(г—1) (х) абсолютно непрерывна и /(г)(х) е Ьр(а, Ь). Если / е ^£2(а Ь), то смешанный ряд по системе {^к(х)} имеет вид

г—1 I \к ^

/(х) - £ /(к)(а)(^^ + £ /г,к^г,к(х), (4)

где

к!

к=0 к=г

/г,к = у /(г) (¿)^ГГк= у /(г) (^)^к—г(5)

аа

Пользуясь определением функций ^г,к(х) (см. (2) и (3)), нетрудно увидеть [18], что система {^Г)к(х)}£=0 является ортонормированной в пространстве ^£2(аЬ) относительно скалярного произведения (1), величины /(к)(а)(к = 0,1,..., г — 1) и /г,к(к = г, г +1,...) являются коэффициентами Фурье функции /, а смешанный ряд (4) представляет собой ряд Фурье функции / е ^£2(а Ь) по этой системе. Мы будем называть систему {^Г)к(х)}£=0 ортонормированной по Соболеву относительно скалярного произведения (1) и порожденной ортонормированной системой {^к(х)}£=0. Таким образом, понятие смешанного ряда по ортонормированной системе {^к(х)}£=0, введенное в работах [7-18], совпадает с понятием ряда Фурье по системе {^Г)к(х)}£=0, ортонормированной по Соболеву относительно скалярного произведения (1) и порожденной ортонормированной системой {^к(х)}£=0.

Следует отметить, что смешанные ряды по классическим ортогональным полиномам были введены в работах [7-18] как альтернативный рядам Фурье по тем же полиномам аппарат решения задач, в которых требуется одновременно приближать дифференцируемую функцию и несколько ее производных. В частности, такая задача часто возникает при решении дифференциальных уравнений

Ь

Ь

численно-аналитическими (спектральными) методами [19,20]. Смешанные ряды по ортонормирован-ным системам функций оказались естественным и весьма эффективным средством решения краевых задач для дифференциальных уравнений спектральными методами. По этому поводу мы можем отослать, например, к работе [21].

В настоящей статье мы рассмотрим дискретный аналог скалярного произведения (1) следующего вида:

г —1 те

/,д> = ^ Дк / (о)дк 0(0) + ^ дг / . )дг д. )р.), (6)

к=0 7=0

где функции / и д заданы на множестве О = {0,1,..., }, р = р(.) — дискретная весовая функция, за-

г— 1

данная на множестве О. В случае, когда г = 0, мы будем считать, что Е Дк/(0)Дкд(0) = 0. При г ^ 1

к=0

особой точкой в скалярном произведении (6) является х = 0, в которой контролируется поведение соответствующих ортогональных по Соболеву полиномов дискретной переменной благодаря присут-

г— 1

. к х/гЛ Л к

ствию в (6) выражения Е Дк/(0)Дкд(0). В настоящей статье рассматриваются системы дискретных

к=0

функций, ортонормированных по Соболеву относительно скалярного произведения (6), порожденных заданной ортонормированной системой функций дискретной переменной. Основное внимание будет уделено изучению свойств полиномов, ортогональных по Соболеву, порожденных классическими ортогональными полиномами Мейкснера дискретной переменной.

1. СИСТЕМЫ ДИСКРЕТНЫХ ФУНКЦИЙ, ОРТОНОРМИРОВАННЫХ ПО СОБОЛЕВУ, ПОРОЖДЕННЫЕ ОРТОНОРМИРОВАННЫМИ ФУНКЦИЯМИ

Перейдем к конструированию дискретных функций, ортонормированных по Соболеву относительно скалярного произведения (6), порожденных заданной системой {^к(х)}£=0, ортонормированной на дискретном множестве О = {0,1,...} с весом р(х). Для этого нам понадобятся некоторые обозначения и понятия. Если целое к ^ 0, то положим а[к] = а(а — 1) ■ ■ ■ (а — к + 1), а[0] = 1 и рассмотрим следующие функции:

х[к]

^г,к(х) = к = 0,1,... ,г — 1, (7)

к!

х — г

Е (х — 1 — *)[г—1]^к—г(*), Г ^ х,

^г,к(х) = { (г—1)! ^ 7 г У " (8)

0, х = 0,1,..., г — 1,

которые определены на сетке О = {0,1,...}. Рассмотрим некоторые важные разностные свойства системы функций ^г>к(х), определенных равенствами (7) и (8). Введем оператор конечной разности Д/: Д/(х) = /(х + 1) — /(х) и положим Д^+1 /(х) = ДД^/(х). Имеет место следующая

Лемма 1. Имеют место равенства

Д"^г,к (х) = <

^)к—^(х), если 0 ^ V ^ г — 1, г ^ к,

г (х), если V = г ^ к, (9)

^)к—и(х), если V ^ к < г,

0, если к < V ^ г.

Доказательство. Справедливость утверждения леммы 1 при г = 1 почти очевидна, поэтому мы будем считать, что г ^ 2. Прежде всего заметим, что если /(х) = (х — 1 — £)[г—1], то Д/(х) = (г — 1)(х — 1 — £)[г—2], поэтому для г ^ х, к в силу (8) имеем:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

.. /х — г+1 ^ .

Д^г,к (х) = —-— (х — ^)[г —1] ^к —г (*) — ^ (х — 1 — ^)[г —1] ^к —г (*)

(г )! V + —п + —п

(г — 1) ¿=0

г=0 ¿=0

х —г+1

х

.. х —г+1

Е ((х — ^)[г —1] — (х — 1 — *)[г—111) ^к —г(*) =

1

x — r+1

X

(r - 2)! ^

v ' t=0

(x - 1 - í)[r —2]^k— 1 — (r — i)(t) = ^ — i)fc — i(x). (10)

Отсюда убеждаемся в справедливости первого из равенств (9) для г ^ х. Справедливость равенства Д^г,к(х) = ^г-1;к-1(х) для х = 0,1,... ,г — 2 очевидна. Остается проверить первое из равенств (9) для х = г — 1. Но в этом случае мы имеем:

1 x—r+1 (r_ i)[r—i]

A^,k(x) = ^r,fc(x + 1) = --T V (x - í)[r—1]^fc—r(í) = -^k—r(0) = ^fc—r(0),

(r - 1)! (r - 1)!

( ^ (X - ir J^k —r(í) = , 1

(r - 1)! (r - 1

1 x—-T+1 (r - 2)[r —2]

^r— 1,k — 1(x) = (r - 2)! t=0 (x - 1 - í)[r —2]^(k —1) —(r—1) (í) = (r - 2)! ^k —r (0) = ^k —r (0),

поэтому мы снова находим A^r>k(x) = 1;k—1(x). Таким образом, полностью доказано первое из равенств (9) для 0 ^ v ^ 1. Его справедливость для 2 ^ v ^ r - 1 выводим по индукции.

Рассмотрим второе из равенств (9). В силу уже доказанного первого из равенств (9) и того, что второе из них для r = 1 легко проверяется, имеем

Ar ^r,k (x) = AAr—1Vr,k (x) = A^1,k—r+1 (x) = ^k—r (x).

Тем самым мы доказали второе из равенств (9). Третье и четвертое равенства из (9) непосредственно вытекают из (7). Лемма 1 полностью доказана. □

Пусть 0 ^ r — целое. Обозначим через lp пространство дискретных функций f = f(x), заданных на сетке О = {0,1,...}, в которых скалярное произведение < f,g > определено с помощью равенства (6). Рассмотрим задачу об ортогональности, нормированности и полноте в lp системы {^r,k(x)}£=0, состоящей из функций, определенных равенствами (7) и (8).

Теорема 1. Предположим, что функции (x) (k = 0,1,...) образуют полную в lp ортонорми-рованную систему c весом p(x). Тогда система {^r,k (x)}£=0, порожденная системой {^k (x)}£=0 посредством равенств (7) и (8), полна в lp и ортонормирована относительно скалярного произведения (6).

Доказательство. Из равенства (8) следует, что если r ^ k и 0 ^ v ^ r - 1, то Av^r>k(x))x=0 = 0, поэтому в силу (6) и (9) имеем:

те те

<^r,k > = Ar ^r,k (x)Ar (x)p(x) = ^k —r (x)^¿ — r (x)p(x) = , k, l ^ r, (11)

x=0 x=0

r—1

<^r,k> = Av^r,k(0)AV^r,í(0) = , k, l < r. (12)

v=0

Очевидно также, что

(^г,к, ) = 0, если к<г ^ I, или 1<г ^ к. (13)

Это означает, что функции ^Г)к(£) (к = 0,1,...) образуют в 1р ортонормированную систему относительно скалярного произведения (6). Чтобы проверить полноту этой системы в 1р предположим, что для функции / = /(х) £ 1р имеют место равенства

(^,к,/) = 0, к = 0,1,....

Тогда, во-первых, в силу того что 0 = (^Г;к, /) = Дк/(0) при к = 0,..., г — 1, имеем /(^) = 0 для всех ^ = 0,..., г — 1. Во-вторых, из равенств (^Г;к, /) = 0, к = г, г + 1,... и полноты в 1р исходной системы (£) (к = 0,1,...) следует, что Дг/(х) = 0 (х £ О) и поэтому / совпадает с алгебраическим полиномом степени не выше г — 1. Из этих двух фактов вытекает, что /(х) = 0 (х £ О). Теорема 1 доказана. □

Систему функций ^Г)к(£) (к = 0,1,...) мы будем называть системой, ортонормированной по Соболеву относительно скалярного произведения (6).

1

Из теоремы 1 следует, что система дискретных функций ^г>к(£) (к = 0,1,...) является ортонорми-рованным базисом (ОНБ) в пространстве , поэтому для произвольной функции /(х) £ мы можем записать равенство

те

/(х) = £</,^г,к>^г,к (х), (14)

к=0

которое представляет собой ряд Фурье функции /(х) £ по системе {^г,к(¿)}те=0, ортонормированной по Соболеву. Поскольку коэффициенты Фурье </, ^г,к > имеют вид

г—1

/г,к = </, ^г,к> = Е Д"/(0)Д"^г,к(0) = Дк/(0), к = 0,..., г — 1, (15)

V=0

те те

/г,к = </, ^г,к > = £Дг / (. )Дг ^г,к (.)р(.) ^Дг / (.>к —г (. )р(. ), к = г, г +1,..., (16)

7=0 7=0

то равенство (14) можно переписать в следующем смешанном виде

г—1

/(х) = £ Дк/+ Е /г,к^г,к(х), х £ О. (17)

к=0 к=г

Поэтому ряд Фурье по системе {^г>к(£)}£= мы будем, следуя [7-18], называть смешанным рядом по исходной ортонормированной {^к (£)}те=0. Отметим некоторые важные свойства смешанных рядов (17) и их частичных сумм вида

г—1

^г>п(/,х) = 5>к/+Е /г,к^г,к(х). (18)

к=0 к=г

Из (17) и (18) с учетом равенств (9) мы можем записать (0 ^ V ^ г — 1, х £ О)

г—^1 х[к] те Д"/(х)= Е Дк+:/(0)^Т+ Е /г,к+:^г —:,к (х), (19)

к=0 к=г—V

г —: — 1 х[к] п—: Д" (/,х)= ^ Дк + : / (0)^ + Е /г,к + :^г — :,к (х), (20)

к=0 к=г—V

Д Yг,те(/, х) ^Т —:,п—V(Д х). (21)

2. НЕКОТОРЫЕ СВЕДЕНИЯ О ПОЛИНОМАХ МЕЙКСНЕРА

При конструировании полиномов, ортогональных по Соболеву и порожденных классическими многочленами Мейкснера, нам понадобится ряд свойств этих многочленов, которые мы приведем в настоящем параграфе. Для произвольного а и 0 < д < 1 положим

р(х) = р(х; а) = дх Г(^+ + +)1) (1 — д)а+1, (22)

Ма (х,д) = прх)ДП{ р(х)х[п]}, (23)

где Дп/(х) — конечная разность п-го порядка функции /(х) в точке х, т.е. Д0/(х) = /(х), Д1 /(х) = Д/(х) = /(х + 1) — /(х), Дп/(х) = ДДП—1 /(х) (п ^ 1), а[0] = 1, а[к] = а(а — 1) ■ ■ ■ (а — к + 1) при к ^ 1. Для каждого 0 ^ п равенство (23) определяет [22], алгебраический полином степени п, для которого (0,д) = (п+а).

Полные доказательства приведенных после свойств полиномов Мейкснера Ма(х, д) можно найти, например, в [22,23]. Прежде всего отметим, что полиномы Ма(х,д) допускают следующее явное представление:

М? (*,,) = Г(п + а +1' Е /"Л,, (1 — ^ ■ (24)

^ п! к=0 Г(к + а + 1)к!\ д) К '

Если а > —1, то полиномы М^(х) = ма(х,д) (п = 0,1,...) образуют ортогональную с весом р(х) (см. (22)) систему на множестве О = {0,1,...}:

ма(х)ма(х)р(х) = ¿п,кй£(д) 0 < д < 1, а > —1, (25)

где

^(д) = Е р(х) {ма(х)}2 = (п + а)д_пГ(а + 1). (26)

ж=0 V п /

Нам понадобятся также следующие свойства:

дма(х) = ма (х+1) — ма (х) = 1 ма+11 (х), (27)

д

м:(х, д) = Е ма(х, д), (28)

к=0

' . г [Л

ма+-гг (х,д) = 1)' --Г" ма+г-л (х,д), (29)

; =0 ^'

м-1 (х, д) = (п —1)' — Л (—х) м\ _ (х — 1,д), I -целое, 1 ^ I ^ п, (30)

п' \д у

где (а)о = 1, (а) = а(а + 1)... (а + I — 1). Лемма 2. Пусть 0 ^ г. Тогда

(х + г)М Л N 1 ^ г

(хтгй м (х,д) = Е (—д^м* (х,д).

Доказательство. Полагая в (29) а = г, имеем:

м°(х,д) = ¿(—1)^ ^м;-^ (х,д).

;0(х,д) =

V=0

В то же время мы можем записать

( + ) [г] к+г ^^ м (х,д) = 3=0 7; м3,(х,д)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

х,д) = > 7; м0(х,д), (31)

где

1 ^ (+ + г)[г] ^ (4 + г)[г]

= хщ Е o, мГ ^ д)м°в д) = д; (1 — д) Е д4 (^й мг ^ д)м; ^ д)

Е д4 Г(Г (+++)1)(1 — д)г+1мг (*, д)м?(^, д),

(к + г)[г](1 — д)г 4=0у г(* + 1) в частности, 7; = 0 при ^ < к. Из (31) и (32) имеем:

= (к+г)и (1—д)г Е р(',г,д)мг (^,д) Е(—1Г( (х,д) = = ,—к)Е р(^,г,д)[мк (',д)]2 =

= д;(—1);-к Г г У(д) = (—д);-к ( г ук + Лг(г + 1)

(к + г) [г] (1 — — кДк(д) (к + г) [г] (1 — — кД к Дг + ^

(32)

Отсюда и из (31) находим

[г] к+г

(1 — ^ Мкг(х,д) = 7=(—— ^М°(х,д) = ¿0(—д)'(г)м0+г(х,д).

Лемма 2 доказана. □

Из (25) следует, что полиномы

та(х) = та(х, д) = (Ла(д))—1/2ма(х, 5) (33)

образуют ортонормированную систему на множестве О с весом р(х) = р(х, а, д), т.е.

У^ та(х)та(х)р(х) = 5П)к. 0 < д < 1, а > —1. (34)

хе^

3. ОРТОГОНАЛЬНЫЕ, ПО СОБОЛЕВУ, ПОЛИНОМЫ, ПОРОЖДЕННЫЕ МНОГОЧЛЕНАМИ МЕЙКСНЕРА

Из равенства (34) следует, что если а > —1, то полиномы та(х, д) (п = 0,1,...) образуют ортонормированную на О = {0,1,...} с весом р(х) систему. Эта система порождает на О систему полиномов так+г(х, д) (к = 0,1,...), определенных равенством

1 х—г

так+г(х,д) = ( _ 1)! Е(х — 1 — ¿)[г—1]т^, д). (35)

(г )! ¿=0

Кроме того, определим полиномы

х

та,к (х,д) = , к = 0,1,...,г — 1. (36)

Покажем, что полином т0,к+г(х,д) обращается в нуль, если х £ {0,1, ...,г — 1}. С этой целью рассмотрим следующий дискретный аналог формулы Тейлора (х £ {г, г + 1,...}):

1 х—г

¥ (х) = 1(¥, х) + --— У> — 1 — *)[г—1] Дг ¥ (*), (37)

(г — 1)! ^

х

¿=0

где

^^х) = ¥ (0) + х + « х[2] + ... + х[г—1]. (38)

Так как для функции ¥(х) = х[1+г], где целое I ^ 0, имеем Дг¥(х) = (I + г)[г]х[1] и 1(¥,х) = 0 , то из (37) следует, что

х—г

[1+г]

--- V(х — 1 — £)[г—1]г[1] = --^-- V(х — 1 — £)[г—1]Дг¥т = ,, (39)

(г — 1)! ¿=0 (I + г)[г](г — 1)! ¿0 (I + г)[г]

С другой стороны, функция х[1+г] обращается в нуль в узлах х £ {0,1,...,г — 1}, поэтому наше утверждение вытекает из того, что полином та+г(х,д) в силу (24) и (33) можно представить в виде линейной комбинации функций вида х[1]. Поэтому из теоремы 1 и свойства (29) непосредственно выводим следующий результат.

Теорема 2. Если а > —1, то система полиномов та к(х,д) (к = 0,1,...), порожденная мно-

а г

гочленами Мейкснера та(х, д) (п = 0,1,...) посредством равенств (35) и (36), полна в и ортонормирована относительно скалярного произведения (6).

4. ДАЛЬНЕЙШИЕ СВОЙСТВА ПОЛИНОМОВ так(х, я)

Перейдем к исследованию дальнейших свойств полиномов та к(х, д) (к = 0,1,...). Речь, в первую очередь, идет о том, чтобы получить представление полиномов та к(х, д), которое не содержит знаков

суммирования с переменным верхним пределом типа (35). С этой целью применим формулу (37) к полиному Г(х) = мк+7(х,д) и запишем

г (х) = дг_1 (г, х) + — — 1 — ;)[г_1] Дг ма+г (*, д).

( )! 4=0

(40)

Вместо Дгмк+гг(х, д) подставим его значение, которое согласно формуле (27) равно ()гма(х,д), тогда из (37) получим:

Г (х) — ^г_1(Г,х) =

д— 1 г 1

д ^ (г — 1)! 4=0

— — 1 — ^)[г-1]ма(^,д).

Сопоставляя (35) и (41) с (33), находим

(^т) {ха(д)}1/2 та,к+г(х, д) = Г(х) — ^г_1 (Г, х). Из (42) имеем (Г(х) = м^+Г(х,д))

та, к+г(х,д)= (д^) {Ха(д)}_1/2 [Г(х) — (Г,х

Далее, в силу (27)

Д' ма+_гг (х, д) =

д ^ ма+_гг_+^^(х, д),

поэтому из (24) находим

Д' ма+_гг (0, д) =

д — 1

Г(к + а + 1)

д у (к + г — V )'Г(^ — г + а + 1)

= А

Равенства (38) и (44), взятые вместе, дают

^ А к хИ Г(х) — ^г_1 (Г, х) = ма+_гг(х, д) — Е

v=0

Подставив это выражение в (43), находим

тг,к+г (х,д) = ^ д — 1 ^

К (д)}

-1/2

ма+7 (х, д) —

ЕАг,к^х[ ]

V!

V=0

, к = 0,1,.

(41)

(42)

(43)

(44)

(45)

(46)

Еще одно важное представление для полиномов так+г(х,д) можно получить, если мы обратимся к равенствам (24) и (33) и запишем

т (х,д) =

Е

к[1] х[1]

Г(к + а + 1) к' К(д)}1/2 ¿0 Г(1 + а + д

1

1 - -

(47)

Подставим это выражение в (35) и воспользуемся равенством (39). Это приводит к следующему

явному виду для полиномов так+г(х, д) (к = 0,1,...)

Г(к + а + 1)

та , к+г (x, д) =

к к[1] х[г+1]

Е

к' К(д)}1/2 ¿0 Г(1 + а + 1)1'(г + 1)[г] V д

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1

(48)

V

г

к

к

1

5. ПОЛИНОМЫ {т°к(х, я)}^

Рассмотрим частный случай, который соответствует выбору а = 0. Заметим, что если а = 0, то из (44) имеем = 0 при всех V = 0,1,..., г — 1. Поэтому из (46) имеем:

т0,к+г(х,д)= (д^) К(д)} —1/2 М—+гг(х, д), к = 0,1,.... (49)

Далее, если мы обратимся к равенству (30), то можем записать

к / 1 \г

Мк"+гг(x, д) = (^1 — д) х[г]Мкг(х — г д). (50)

Из (49) и (50) находим

т0,к+г(x, д) = (к +' г)! ВД 2 х[г]мк(х — г, д),

С учетом (26) и (33) этому равенству можно придать также следующий вид:

т0,к+г(х,д) = ((к + г)[г]) —1/2х[г]тк(х — г,д), к = 0,1,.... (52)

к!

+ (д)} 2 ^Мк(х — г,д), к = 0,1,.... (51)

Воспользовавшись леммой 2 и равенством (51), мы также можем получить для т0к+г(х,д) представления в виде линейных комбинаций полиномов Мейкснера т^+Дх, д) (г = 0,1,..., г). Наконец, если 0 ^ к ^ г — 1, то в силу определения (36)

х[к]

т°,к (х,д) = ~кг. (53)

г,к

Из теоремы 1 следует, что система полиномов т0 к(х, д) (к = 0,1,...) является ортонормированным базисом (ОНБ) в пространстве , поэтому для произвольной функции /(х) £ мы можем записать равенство

те

/(х) = Е</, т0,к>т0,к(x, g), (54)

к=0

которое представляет собой ряд Фурье функции /(х) £ по системе {т0,к(х, д)}те=0, ортонорми-рованной по Сободлеву относительно скалярного призведения (6). Поскольку коэффициенты Фурье </, т°к> имеют вид

г-1

/г ,к = </, т0,к> = ^ /(0^т0,к(0, д) = Дк/(0), к = 0,..., г — 1, (55)

V=0

те

/г,к = </, т0,к> = ^ Дг/(.)т°—г(., д)р(.), к = г, г + 1,..., (56)

7=0

то равенство (54) можно переписать в следующем смешанном виде:

г —1 х[к] те /(х) = £ Дк/(0)— + ^ /г,кт0,к(х, д), х £ О. (57)

к=0 к=г

6. РАЗНОСТНЫЕ СВОЙСТВА ЧАСТИЧНЫХ СУММ ФУРЬЕ ПО СИСТЕМЕ {т^ (х, д)}°=0

Основные разностные свойства сумм Фурье по полиномам так(х,д), которые согласно (18) имеют

вид

^ , х[к] ^ а к а

^(/,х) = £ Дк/(0) — + ^/г,кта,к(х, д), (58)

к=0 к=г

где

0

/г,к = </, т0,к> = ^ Дг/(.)т°—г(., д)р(.), к = г, г + 1,....

7=0

выражены равенствами (19)-(21). Для системы {так(х,д)}£=0 они принимают вид (0 ^ V ^ г — 1)

r — V — 1

Avf(x) = E A'+Vf(°)^T + E fr,k+vma—V,k(x,q),

k=0 r —v —1

k=r —v n—v

AvYran(f,x) = E Ak+Vf(°)^r + E fr,k+vma—v,k(x,q),

k=0

k=r —v

A Yr,n (f, x) Yr — v,n—v (A f,x)-

Из (19) и (20) мы также можем записать для n ^ r > v ^ 0

(59)

(60) (61)

Avf (x) - Av Y" (f, x) =

k=n—v+1

fr,k+v mr — v,

(62)

Равенство (62) дает выражение для погрешности, проистекающей в результате замены конечной разности ДV/(х) ее приближенным значением ДV(/,х). При решении задачи об оценке этой погрешности возникает вопрос об асимптотических свойствах полиномов т0^к(х,д). Этот вопрос, в свою очередь, сводится, как это было показано в пп. 4 и 5 настоящей статьи, к задаче об асимптлотических свойствах полиномов т^ (х, д), которые весьма подробно исследованы в [22].

Библиографический список

1. Iserles A., Koch P. E., Norsett S. P., Sanz-Serna J. M. On polynomials orthogonal with respect to certain Sobolev inner products // J. Ap-prox. Theory. 1991. Vol. 65, iss. 2. P. 151-175. DOI: 10.1016/0021-9045(91)90100-0.

2. Marcellan F., Alfaro M, Rezola M. L. Orthogonal polynomials on Sobolev spaces: old and new directions // J. Comput. Appl. Math. 1993. Vol. 48, iss. 1-2. P. 113-131. D0I:10.1016/0377-0427(93)90318-6.

3. Meijer H. G. Laguerre polynimials generalized to a certain discrete Sobolev inner product space // J. Approx. Theory. 1993. Vol. 73, iss. 1. P. 1-16. DOI: 10.1006/jath.1993.1029.

4. Kwon K. H., Littlejohn L. L. The orthogonality of the Laguerre polynomials (x)} for positive integers k // Ann. Numer. Anal. 1995. Vol. 2. P. 289303.

5. Kwon K. H., Littlejohn L. L. Sobolev orthogonal polynomials and second-order differential equations // Ann. Numer. Anal. 1998. Vol. 28. P. 547594.

6. Marcellan F., Yuan Xu On sobolev orthogonal polynomials. arXiv: 6249v1 [math.C.A]. 25 Mar 2014. P. 1-40.

7. Шарапудинов И. И. Приближение дискретных функций и многочлены Чебышева, ортогональные на равномерной сетке // Матем. заметки. 2000. Т. 67, № 3. С. 460-470. DOI: 10.4213/mzm858.

8. Шарапудинов И. И. Приближение функций с переменной гладкостью суммами Фурье -Лежандра // Матем. сб. 2000. Т. 191, № 5. С. 143160. DOI: 10.4213/sm480.

9. Шарапудинов И. И. Аппроксимативные свойства операторов Yn+2r (f) и их дискретных аналогов // Матем. заметки. 2002. Т. 72, вып. 5. С. 765-795. DOI: 10.4213/mzm466.

10. Шарапудинов И. И. Смешанные ряды по ультрасферическим полиномам и их аппроксимативные свойства // Матем. сб. 2003. Т. 194, № 3. С. 115— 148. 001: 10.4213Дт723.

11. Шарапудинов И. И. Смешанные ряды по ортогональным полиномам / Дагестан. науч. центр РАН. Махачкала, 2004. 176 с.

12. Шарапудинов И. И. Смешанные ряды по полиномам Чебышева, ортогональным на равномерной сетке // Матем. заметки. 2005. Т. 78, вып. 3. С. 442-465. 001: 10.4213/тгт2599.

13. Шарапудинов И. И. Аппроксимативные свойства смешанных рядов по полиномам Лежандра на классах Wr // Матем. сб. 2006. Т. 197, № 3. С. 135-154. 001: 10.4213Дт1539.

14. Шарапудинов Т. И. Аппроксимативные свойства смешанных рядов по полиномам Чебышева, ортогональным на равномерной сетке // Вестн. Дагестан. науч. центра РАН. 2007. Т. 29. С. 12-23.

15. Шарапудинов И. И. Аппроксимативные свойства средних типа Валле - Пуссена частичных сумм смешанных рядов по полиномам Лежандра // Мат. заметки. 2008. Т. 84, № 3. С. 452-471. 001: 10.4213/тгт5541.

16. Шарапудинов И. И., Муратова Г. Н. Некоторые свойства г-кратно интегрированных рядов по системе Хаара // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2009. Т. 9, вып. 1. С. 68-76.

17. Шарапудинов И. И., Шарапудинов Т. И. Смешанные ряды по полиномам Якоби и Чебышева и их дискретизация // Матем. заметки. 2010. Т. 88, вып. 1. С. 116-147. 001: 10.4213/тгт6607.

18. Шарапудинов И. И. Системы функций, ортогональных по Соболеву, порожденные ортогональными функциями // Современные проблемы теории функций и их приложения : материалы 18-й

междунар. Сарат. зимней школы. Саратов : ООО Изд-во «Научная книга», 2016. C.329-332.

19. Trefethen L. N. Spectral methods in Matlab. Philadelphia : SIAM, 2000. 160 p.

20. Trefethen L. N. Finite difference and spectral methods for ordinary and partial differential equation. Cornell University, 1996. 325 p.

21. Магомед-Касумов Р. Г. Приближенное решение обыкновенных дифференциальных уравнений с использованием смешанных рядов по системе Ха-ара // Современные проблемы теории функций и их приложения : материалы 18-й междунар. Сарат.

зимней школы. Саратов : ООО Изд-во «Научная книга», 2016. C. 176-178.

22. Шарапудинов И. И. Многочлены, ортогональные на дискретных сетках. Махачкала : Изд-во Даг. гос. пед. ун-та, 1997. 252 с.

23. Gasper G. Positivity and special function // Theory and Application of Special Functions : Proceedings of an Advanced Seminar Sponsored by the Mathematics Research Center, the University of Wisconsin-Madison, March 31-April 2, 1975 / ed. R. Askey. N.Y. ; San Francisco ; L. : Academic Press, Inc., 1975. P. 375-433.

Образец для цитирования:

Шарапудинов И. И., Гаджиева З. Д. Полиномы, ортогональные по Соболеву, порожденные многочленами Мейкс-нера // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2016. Т. 16, вып. 3. С. 310-321. 001: 10.18500/1816-9791-2016-16-3-310-321.

Sobolev Orthogonal Polynomials Generated by Meixner Polynomials 1.1. Sharapudinov1, Z. D. Gadzhieva2

1 Idris I. Sharapudinov, Dagestan Scientific Center RAS, Vladikavkaz Scientific Center RAS, Dagestan State Pedagogical University, 45, M. Gadzhieva st., 367032, Makhachkala, Russia, sharapud@mail.ru

2Zulfiya D. Gadzhieva, Dagestan Scientific Center RAS, Dagestan State Pedagogical University, 45, M.Gadzhieva st., 367032, Makhachkala, Russia, zulfiya.gadzhieva.1976@mail.ru

The problem of constructing Sobolev orthogonal polynomials mr,n(x,q) (n = 0,1,...), generated by classical Meixner's

[k]

polynomials is considered. They can by defined using the following equalities m"k (x,q) = ^r, x = x(x -1) ■ ■ ■ (x-k+1),

x — r

k = 0,1,... ,r - 1, m"k+r(x,q) = (r-n! ^ (x - 1 - i)[r-1]m£(i,q), wherem^(t,q) denote Meixner's polynomial

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(r )! t=0

of degree k, orthonormal on Q = {0,1,...} with weight p(x) = qxr(rx("X+|)1)(1 - q)a+1. Polynomials mr,n(x,q), (n = 0,1,...) are orthonormal on Q = {0,1,...} with respect to the inner product

r-1 oo

(m?,n!m?,m} = Akmr,n(0,q)Akm?,m(0,q) + ^ Arm?in(j', q)Arm?,m (j,q)p(j). k=0 j=0

For mr,n(x, q) we obtain the explicit formula that contains the Мeixner polynomial M^-r (x, q):

(z,q) = ( —^TY (hfc (q)} 1/2

q — I

Г— A ,

-r

v=0

k = 0, I,...,

k

(q - 1V r(k + a + 1) „,a, , r(n + a + 1) ^A n[klx[kl / 1

where= - -r—--; ,M^(x,q) = -:-V , 1 - -

V q / (k + r — v)!r(v — r + a + 1) n! 0 r(k + a + 1)k! \ q

ha(q)= ^ + ^q-n r(a + 1).

Key words: orthogonal Sobolev polynomial, Meixner polynomials orthogonal on the grid, approximation of discrete functions, mixed series in Meixner polinomials orthogonal on a uniform grid.

References

1. Iserles A., Koch P. E., Norsett S. P., Sanz- directions. J. Comput. Appl. Math., 1993, vol. 48, Serna J. M. On polynomials orthogonal with re- iss. 1-2, pp. 113-131. DOI: 10.1016/0377-0427(93) spect to certain Sobolev inner products. J. Approx. 90318-6.

Theory, 1991, vol. 65, iss. 2, pp. 151-175. DOI: 3. Meijer H. G. Laguerre polynimials generalized to a

10.1016/0021-9045(91)90100-0. certain discrete Sobolev inner product space. J. Ap-

2. Marcellan F., Alfaro M., Rezola M. L. Orthogo- prox. Theory, 1993, vol. 73, iss. 1, pp. 1-16. DOI: nal polynomials on Sobolev spaces: old and new 10.1006/jath.1993.1029.

4. Kwon K. H., Littlejohn L. L. The orthogonality of the Laguerre polynomials {Ll-k)(x}} for positive integers k. Ann. Numer. Anal., 1995, vol. 2, pp. 289-303.

5. Kwon K. H., Littlejohn L. L. Sobolev orthogonal polynomials and second-order differential equations. Ann. Numer. Anal., 1998, vol. 28, pp. 547594.

6. Marcellan F., Xu Y. On sobolev orthogonal polynomials. arXiv:6249v1 [math.C.A]. 25 Mar 2014, pp. 1-40.

7. Sharapudinov I. I. Approximation of discrete functions and Chebyshev polynomials orthogonal on the uniform grid. Math. Notes, 2000, 67, iss. 3, pp. 389-397. DOI: 10.1007/BF02676675.

8. Sharapudinov I. I. Approximation of functions of variable smoothness by Fourier - Legendre sums. Sb. Math., 2000, vol. 191, no. 5, pp. 759-777. DOI: 10.1070/SM2000v191n05ABEH000480.

9. Sharapudinov I. I. Approximation properties of the operators Yn+2r (f} and of their discrete analogs. Math. Notes, 2002, vol. 72, iss. 5, pp. 705-732. DOI: 10.1023/A:1021421425474.

10. Sharapudinov I. I. Mixed series in ultraspherical polynomials and their approximation properties. Sb. Math., 2003, vol. 194, no. 3, pp. 423-456. DOI: 10.1070/SM2003v194n03ABEH000723.

11. Sharapudinov I. I. Smeshannye rjady po ortog-onal'nym polinomam [Mixed series of orthogonal polynomials]. Makhachkala, Dagestan Scientific Center RAS, 2004, 176 p. (in Russian).

12. Sharapudinov I. I. Mixed series of Chebyshev polynomials orthogonal on a uniform grid. Math. Notes, 2005, vol. 78, iss. 3, pp. 403-423. DOI: 10.1007/s11006-005-0139-3.

13. Sharapudinov I. I. Approximation properties of mixed series in terms of Legendre polynomials on the classes Wr. Sb. Math, 2006, vol. 197, no. 3, pp. 433-452. DOI: 10.1070/SM2006v197n03A BEH003765.

14. Sharapudinov T. I. Approximative properties of mixed series by Chebyshev polynomials, orthogonal on an uniform net. Vestnik Dagestan. nauch. centra RAN, 2007, vol. 29, pp. 12-23 (in Russian).

15. Sharapudinov I. I. Approximation properties of the Vallee - Poussin means of partial sums of a mixed series of Legendre polynomials. Math.

Notes, 2008, vol. 84, iss 3-4, pp. 417-434. DOI: 10.1134/S0001434608090125.

16. Sharapudinov I. I., Muratova G. N. Same Properties r-fold Integration Series on Fourier - Haar System. Izv. Saratov Univ. (N.S), Ser. Math. Mech. Inform., 2009, vol. 9, iss. 1, pp. 68-76 (in Russian).

17. Sharapudinov I. I., Sharapudinov T. I. Mixed series of Jacobi and Chebyshev polynomials and their discretization. Math. Notes, vol. 88, iss. 1-2, pp. 112139. DOI: 10.1134/S0001434610070114.

18. Sharapudinov I. I. Sistemy funktsii, ortogonal'nykh po Sobolevu, porozhdennye ortogonal'nymi funktsi-iami [System functions orthogonal on Sobolev generated orthogonal functions]. Sovremennye prob-lemy teorii funktsii i ikh prilozheniia : materi-aly 18-i mezhdunar. Sarat. zimnei shkoly [Modern problems of function theory and their applications : Proc. 18th Intern. Sarat. Winter School]. Saratov, OOO Izd-vo "Nauchnaia kniga", 2016, pp. 329-332.

19. Trefethen L. N. Spectral methods in Matlab. Fhiladelphia, SIAM, 2000, 160 p.

20. Trefethen L. N. Finite difference and spectral methods for ordinary and partial differential equation. Cornell University, 1996, 325 p.

21. Magomed-Gasimov M. G. Priblizhennoe reshe-nie obyknovennykh differentsial'nykh uravnenii s ispol'zovaniem smeshannykh riadov po sisteme Khaara [Approximate solution of ordinary differential equations with the use of mixed series on the Haar system]. Sovremennye problemy teorii funktsii i ikh prilozheniia : materialy 18-i mezhdunar. Sarat. zimnei shkoly [Modern problems of function theory and their applications : Proc. 18th Intern. Sarat. Winter School]. Saratov, OOO Izd-vo "Nauchnaia kniga", 2016, pp. 176-178.

22. Sharapudinov I. I. Mnogochleny, ortogonal'nye na dickretnyh setkah [Polynomials orthogonal on discrete grids]. Mahachkala, Izd-vo Dag. Gos. Ped. Un-ta, 1997, 252 p. (in Russian).

23. Gasper G, ppositivity and special function. Theory and Application of Special Functions : Proc. of an Advanced Seminar Sponsored by the Mathematics Research Center, the University of Wisconsin-Madison, March 31-April 2, 1975, ed. R. Askey, New York, San Francisco, London, Academic Press, Inc., 1975, pp. 375-433.

Please cite this article in press as:

Sharapudinov I. I., Gadzhieva Z. D. Sobolev Orthogonal Polynomials, Generated by Meixner Polynomials. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2016, vol. 16, iss. 3, pp. 310-321 (in Russian). DOI: 10.18500/1816-97912016-16-3-310-321.

i Надоели баннеры? Вы всегда можете отключить рекламу.