Некоторые специальные двумерные ряды по системе {sinx sinkx} и их аппроксимативные свойства Текст научной статьи по специальности «Математика»

This work was supported by the Russian Foundation for Basic Research (project no. 14-01-00417) and by the Grant of the President of the Russian Federation for state support of leading scientific schools (project no. H^-1096.2014.1).

References

1. Saks S. Teoriia integrala [Theory of the integral]. Moscow, Faktorial Press, 2004, 496 p. (in Russian).

2. Stein I. Singuliarnye integraly i differentsial'nye svoistva funktsii [Singular integrals and differential properties of functions]. Moscow, Mir, 1973, 342 p. (in Russian).

3. Lukashenko T. P., Skvortsov V. A., Solodov A. P. Obobshchennye integraly [Generalized integrals]. Moscow, Knizhnyi dom «LIBROKOM», 2010, 280 p. (in Russian).

4. Natanson I. P. Teoriia funktsii veshchestvennoi pere-mennoi [Theory of functions of a real variable]. S.-Peter-burg, Lan', 2008. 560 p. (in Russian).

5. Kashin B. S., Saakyan A. A. Orthogonal series.

Y^K 517.538

Translations of Math. Monographs, vol. 75, Providence, RI, Amer. Math. Soc., 1989. (Rus. ed.: Kashin B. S., Saakian A. A. Ortogonal'nye riady. Moscow, AFC, 1999, 560 p.).

6. Belousov K. V., Lukashenko T. P. O nekotorykh svoistvakh ortorekursivnykh razlozhenii funktsii mnogikh peremennykh po sisteme kharakteristicheskikh funktsii brusov [On some properties autorecording expansions of functions of many variables by system characteristic functions beams]. Sovremennye problemy matematiki i mekhaniki [Modern problems of mathematics and mechanics], 2011, vol. 6, no. 1, pp. 52-60 (in Russian).

НЕКОТОРЫЕ СПЕЦИАЛЬНЫЕ ДВУМЕРНЫЕ РЯДЫ ПО СИСТЕМЕ {sin x sin kx}

И ИХ АППРОКСИМАТИВНЫЕ СВОЙСТВА

И. И. Шарапудинов

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

В настоящей статье вводятся двумерные специальные ряды по системе {sin x sin kx}. Показано, что эти ряды выгодно отличаются от двумерных косинус-рядов Фурье тем, что их частичные суммы вблизи границы квадрата [0, п]2 обладают значительно лучшими аппроксимативными свойствами, чем суммы Фурье. Приводится оценка скорости сходимости частичных сумм специального ряда к функциям f (x,y) из пространства четных 2п-периодических по каждой переменной непрерывных функций.

Ключевые слова: специальные ряды по системе {sin x sin kx}, двумерные ряды, покусочная аппроксимация.

ВВЕДЕНИЕ

Представление функций в виде рядов по тем или иным ортонормированным системам с целью последующего их приближения частичными суммами выбранного ортогонального ряда является, пожалуй, одним из самых часто применяемых подходов в теории приближений и ее приложениях. Наряду с задачами математической физики, для решения которых указанный подход является традиционным, появились и продолжают появляться все новые важные задачи, для решения которых также все чаще применяются методы, основанные на представлении функций(сигналов) в виде рядов по подходящим ортонормированным системам (см., например, [1-9]). При этом часто возникает такая ситуация, когда функция (сигнал, временной ряд, изображение и т.д.) / = /(Ь) задана на достаточно длинном промежутке [0,Т] и нам требуется разбить этот промежуток на части [aj,0^+1 ] (^ = 0,1,..., т), рассмотреть отдельные фрагменты функции, определенные на этих частичных отрезках, представить их в виде рядов по выбранной ортонормированной системе и аппроксимировать каждый такой фрагмент частичными суммами соответствующего ряда. Такая ситуация является типичной для задач, связанных с решением нелинейных дифференциальных уравнений численно-аналитическими методами [4,6], обработкой временных рядов и изображений и других [5-7], в которых возникает необходимость разбить заданный ряд данных на части, аппроксимировать каждую часть и заменить приближенно

исходный временный ряд (изображение) функцией, полученной в результате «пристыковки» функций, аппроксимирующих отдельные части. Но тогда в местах «стыка» возникают нежелательные разрывы (артефакты) (см. [8]), которые искажают общий вид временного ряда (изображения). Такая картина непременно возникает при использовании для приближения «кусков» исходной функции сумм Фурье по классическим ортонормированным системам. В работах автора [1,2] введены некоторые специальные ряды по ультрасферическим полиномам Якоби, частичные суммы о^К/, х) которых на концах отрезка [—1,1] совпадают с исходной функцией / (х), т.е. ±1) = / (±1). В качестве одного из

частных случаев таких рядов возникает ряд вида

Ф(0) = аф(0) + sinO^pk sin kO,

k=1

где

Ф(0) + Ф(п) Ф(0) - Ф(п) аф(0) = „—— + „—— cos O,

2

2

p(O) = ф(0) - аф(O),

(1)

2 Г sin кт

Pk = - р(т) —-ат.

П sin т

В работе [2] исследованы, в частности, аппроксимативные свойства ряда (1) в пространстве С|п, состоящем из четных непрерывных 2п — периодических функций.

В настоящей статье мы введем двумерные ряды вида (1) и рассмотрим некоторые свойства частичных сумм этих рядов. Пусть /(х, у) — четная 2п — периодическая по каждой из переменных х и у и интегрируемая на квадрате [0,п]2 функция, которая в точках (гп^'п), г,^ £ Ж принимает конечные значения. Положим

/(0, V) + /(п^) /(0, V) - /(п^)

S(f) = S(f)(u,v) = f (u,v) -H(f) = H(f)(u,v) = f(u,v) -O(f)= O(f)(u,v) = S(u,v) -

2 2 c

f (u, 0) + f (u,n) f (u, 0) - f (u,n)

2 2

S(u, 0) + S(u, n) S(u, 0) - S (u, n)

соэ v, соэ v,

\ f (0, 0) - f (0, n) - f (п, 0) + f (п,п) e(f) = 6(f)(u,v) = v 7— v 4— v 7 cosucos v+

+ f (0, 0) + f (0, n) - f (П, 0) - f (п,п) + f (0, 0) - f (0, n) + f (П, 0) - f (п,п) + + cos u i cos v I

+

f (0, 0) + f (0, n) + f (П, 0) + f (п,п) 4

и заметим, что

f (x, y) = 6(f)(x, y) + O(f)(x, y) + A(f)(x) + B(f)(x) cos y + c(f)(y) + D(f)(y) cosx, A(f = S(f)(u, 0) + S(f)(u,n), B(f)(u) = S(f)(u, 0) - S(f)(u, n)

2

2

с(/)(и) = Н(/)(0,«) -2 Н(/)(п,у), )(и) = Н(/)(0, у) - Н(/)(п,у)

Определим следующие коэффициенты:

(2)

(3)

(4)

(5)

Ok,i(f) = I I O(f)(u,v)sinkusin dudv,

П2

Э1П И Э1П V

0 0

n n

2 f t, „s, s sinku , , , „, 2 f „.. . sinku ,

ak(f) = - A(f)(ub-du, bk(f) = - B(f)(u)--du,

П J sin u П J sin u

00 П П

Ci(f) = 2 i C(f )(u)sin1v dv, di(f) = 2 i D(f )(u)sin1v dv, П J sin v П J sin v

00

П

0

2

2

П П

и рассмотрим ряды:

те те

O(f) (x, y) ~ sin x sin y ^^ ^^ Ok;l (f) sin kx sin ly, k=i1=0

те те

A(f )(x) ~ sin (f) sin kx, B(f )(x) ~ sin (f) sin kx,

k=i k=i

тете

C (f )(y) ~sin (f )sin ly' D(f )(y) ~sin (f )sin ly-

i=i i=i Тогда мы можем сопоставить функции f(x, y) следующий специальный ряд:

тете

f (x, y) ~ B(f )(x, y) + sin x sin y ^^ ^^ ok;l (f) sin kx sin ly+

k=ii=i

тете

+ sin x (afc (f) + bk (f) cos y) sin kx + sin y (q (f) + d (f) cos x) sin ly. (6)

k=i i=i

Ряд (6) будем называть двумерным специальным рядом по системе {sin x sin nx}^=i- В настоящей работе исследованы аппроксимативные свойства таких рядов. Показано, что специальные ряды вида (6) выгодно отличаются от двумерных косинус-рядов Фурье тем, что их частичные суммы вида

n— i m— i

(f) = Sn,m (f)(x,y) = @(f )(x,y) + sin x sin yS (f )sin kx sin ly+

k=i i=i

V—i /U—i

+ sin xj^ (ak (f) + bk (f)cos y) sin kx + sin yj^ (c (f) + di(f) cos x) sin ly (7)

k=i i=i

вблизи границы квадрата [0,п]2 обладают значительно лучшими аппроксимативными свойствами, чем суммы Фурье вида

n m

Sn,m(f)(x,y) = J^S ak,i (f)cos kx cos ly. (8)

k=0 l=0

Чтобы убедиться в этом, достаточно, например, сопоставить оценку для функции Лебега L^m(x, у) частичных сумм S^mf)(x,y) вида (7) , полученную в настоящей работе, с хорошо известной оценкой функции Лебега сумм Фурье Sn,m(f)(x,y) вида (8).

АППРОКСИМАТИВНЫЕ СВОЙСТВА СПЕЦИАЛЬНЫХ РЯДОВ ПО СИСТЕМЕ {sin x sin kx}

Прежде всего покажем, что частичная сумма Sn,m(f) специального ряда (6) вида (7) является проектором на подпространство четных тригонометрических полиномов вида

n m v

T^m(x,y) = ^ X^Pk.i cos kx cos ly + (pk,0 + Pk,i cos y) cos kx+

k=0 l=0 k=n+i

+ S (P0,i + cos x)cos ly, (9)

l=m+i

т.е. если f(x,y) = T^m (x,y), то

Sn,m(f) = Tn;m(x,y)- (10)

Поскольку Sn,m(f) - линейный оператор, то тождество (10) будет доказано, если мы покажем, что оно верно для f (x,y) = cos kx cos ly с (k, l) <G Gn;m, где

Gnu, = {(k, l) : 0 ^ k ^ n, 0 ^ l ^ m} U {(k, l) : 0 ^ k ^ 1, д - m + 1 ^ l ^ д}и U{(k, l) : 0 ^ v - n + 1 ^ l ^ v, 0 ^ l ^ 1}.

Пусть, например, к, / — четные, 0 ^ к ^ n, 0 ^ / ^ m. Положим к = 2i, / = 2j. Тогда из (2)-(5) для f (x,y) = cos kx cos /y имеем:

S(f)(x, y) = cos /y(cos kx - 1), H(f)(x, y) = cos kx(cos /y - 1), (11)

A(f )(x) = S(f )(x, 0) + S(f )(x,n) = c0s kx - 1, B(f)(x) = 0, (12)

B(f)(y)= H(f)(0,y) + H(f)(n,y) =co,- 1, D(f)(y) = 0, (13)

O(f)(x,y) = (cos kx - 1)(cos /y - 1), B(x, y) = 1. (14)

Из (11)-(14) находим

f (x, y) = cos kx cos /y = 1 + O(f )(x, y) + A(f )(x) + C(f )(y), (15)

причем

)(x) k— • • C(f)(y) 1-1 e . . ,1fi.

= > a sm jx, —:- = > рг sm iy, (lb)

Z—/ J cm 1/ z—/

sm x z—' sm y

j=i y г=1

O(f )(х,У) 1-1 о • • • • Л7\

_ /./.Oj sm jx sm iy, (17)

sm x sm y

y j=1i=1

где

2 sin jt 2 sin it

aj = -/ A(f)(t)s-j-tdt = Aj (f), в< = -/ C (fXt^dt = Ci(f), (18)

П J о si— t П J о si— t

4 fn fn si— ju si— iv

a7- вг = ^ O(f)(u,v)-dudv = oi?-. (19)

n2 J 0 J0 si— u si— v

Кроме того,

bj(f) = 0 (1 ^ j ^ k - 1), di = 0 (1 ^ i ^ l - 1). (20)

Из (15)-(20) получаем:

f (x,y) = cos kx cos ¿y = B(x, y) + (f )(x,y) + (f )(x) + (f )(y) = (f )(x,y).

Тем самым мы доказали, что тождество (10) справедливо для T^^y) = cos kx cos ly в том случае, когда v = k, д = l, k,1 — четные. Его справедливость в остальных случаях (k,1) <G G^m проверяется аналогично.

Заметим, что если функция f (x, y) обращается в нуль в точках множества P = {(0,0), (0, п), (п, 0), (п,п)}, то B(x, y) = 0 и, как следствие, равенство (7) принимает следующий вид:

n —1 m —1 V—1

sn,m(f )(x,y) = si— x si— y^ ofc)í(f) si— kx si— ly + si— xj^ (afc (f) + (f)cos y) si— kx+

k=1 1 = 1 k=1 д—1

+ si— y (cí(f) + d¿ (f) cos x) si— ly. (21)

1=1

n— 1

Полагая Rn(t,s) = ^ si— kt si— ks, мы получаем из (21) следующее интегральное представление:

k=1

п п

OVW Л 4si— x si— y f Rn(x,u)Rm (y,v) , , ,

Sn,m(f)(x,y) = -2- O(f)(u,v)-:-:-dudv+

n2 J J si— u si— v

о 0

+ Г[A(f)(u) + B(f)(u) cos y]Rs—-u)du+

П Jo si— u

+

2 sin y

[C(/)(v) + D(/)(v) cosx] (y,v)

dv.

sin v

(22)

Перейдем к вопросу об аппроксимативных свойствах частичных сумм (/)(х,у). Прежде всего заметим, что если Т^т(х, у) — тригонометрический полином вида (9), то в силу (10) имеет место равенство

/ (x, y) - sv,m (/)(x, y) = / (x, y) - Tvm (x, y)+Tvm (x, y) - sv:m(/)(x, y) =

= f (x, y) — Tn,m (x, y) + Sn',m (Tn,m — f )(x, y).

(23)

Если, кроме того, полином Т^т (х, у) вида (9) совпадает с /(х, у) при (х, у) £ Р, то из (22) следует, что

S(T- /)(x y) =

n,m\ n,m J/V 'У/

2 sin x ^

4 sin x sin y

O(Tnvm - /)(u, v)-n(x,u)-m(y,v) du dv+

sin u sin v

0 0

+-

— (x u)

[A(Tn:m - /)(u) + - /)(u) cos y] du+

sin u

2 sin y

—«(y,v)

+-- I [C(T^m - /)(v) + D(Tnv;m - /)(v)cosx]^^^-dv.

sin v

(24)

Обозначим через пространство четных 2п-периодических по каждой из переменных х и у и непрерывных функций / = /(х,у), для которых определена норма ||/1| = шахх у)€К2, ТП^^г(/) — пространство тригонометрических полиномов вида (9), совпадающих с /(х,у) при (х,у) £ Р. Тогда в силу (22)-(24) имеем:

|S(T^-m - /)(u,v)| ^ 3^ - /1|, |H(T^ - /)(u, v)| ^ 3^ - /

и, как следствие,

Аналогично

io(Tn:m - / )(u,v)i ^ 9|Tn:m - /

|A(Tn,'m /)| ^ 3|Tn,:m /|B(Tn,:m /)| ^ 3||Tn,:m |C(Tn,:m - /)| ^ 3|Tn,:m - /|D(Tn,:m - /)| ^ 3||Tn,:m -

(25)

(26)

(27)

(28)

Пусть -Е^т(/) = ^ II/ — Т^т|| — наилучшее приближение функции / £ С|П тригономет-

Тп, т ^ ^гг , т ( ^ )

рическими полиномами Т^т = Т^т(х,у) £ ТП'т(/). Из (25) и (26) с учетом (25)-(28) получаем:

где

|/(x, y) - sn:m(/)(x, y)| ^ :m(/)[1 + 9An(x)Am(y) + 6Av(x) + 6AM(y)],

|-s(t,u)|

(29)

2

As(t) = | sint|-

2

/1 j 0 sin u ni j 0

Далее, пусть E^m(/) _ наилучшее приближение функции / e С|П тригонометрическими полиномами TVm вида (9), не обязательно совпадающими с /(x, y) в точках множества P. Тогда не трудно заметить, что

EV-m (/) ^ K,m(/) ^ 2Ev,m. (зо)

Сопоставляя (30) с (29), мы можем записать

-du = | sint| —

s-1

Е

k=1

sin kt sin ku

sin u

du.

|/(x, y) - SV:m(/)(x, y)| ^ 2EV:m(/)[1 + 9An(x)Am(y) + 6Av(x) + 6AM(y)].

(31)

В связи с оценками (29) и (31) возникает задача об исследовании поведения величины An(t) при n ^ га и -1 ^ t ^ 1. Эта задача была решена в работе автора [2], в которой получена следующая оценка:

An(t) ^ c(1 + ln(1 + n| sint|)). (32)

0

п п

2

0

п

0

С учетом (32) из (31) мы выводим следующий результат: Теорема. Имеет место оценка

|f (x, y) - (f)(x, y)| ^ cE^m(f)[1 + ln(! + n| sinx|)ln(1 + m| siny|)+ + ln(1 + v| sinx|) + ln(1 + д| siny|)].

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

1. Шарапудинов И. И. Предельные ультрасферические ряды и их аппроксимативные свойства // Ма-тем. заметки. 2013. Т. 94, вып. 2. С. 295-309. 001: 10.4213/шгшЮ292.

2. Шарапудинов И. И. Некоторые специальные ряды по ультрасферическим полиномам и их аппроксимативные свойства // Изв. РАН. Сер. матем. 2014. Т. 78, № 5. С. 201-224. 001: 10.4213/!Ш8117.

3. Дедус Ф. Ф., Махортых С. А., Устинин М. Н., Дедус А. Ф. Обобщенный спектрально-аналитический метод обработки информационных массивов. Задачи анализа изображений и распознавания образов. М. : Машиностроение, 1999.

4. Пашковский С. Вычислительные применения многочленов и рядов Чебышева. М. : Наука, 1983. 384 с.

5. Арушанян О. Б., Волченскова Н. И., Залеткин С. Ф. О вычислении коэффициентов рядов Чебышева для решений обыкновенных дифференциальных уравнений // Сиб. электрон. матем. изв. 2011. Т. 8. С. 273-283.

6. Trefethen L. N. Spectral methods in Matlab. Philadelphia : SIAM, 2000.

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

8. Mukundan R., Ramakrishnan K. R. Moment functions in image analysis. Theory and Applications. Singapore : World Scientific, 1998.

9. Malvar H. S. Signal processing with lapped transforms. Boston ; London : Artech House, 1992.

Some Special Two-dimensional Series of {sin x sin kx} System and Their Approximation Properties

1.1. Sharapudinov

Daghestan Scientific Centre of Russian Academy of Sciences, 45, Gadgieva str., Makhachkala, Republic of Dagestan, 367000, Russia, sharapudinov@mail.ru

In present paper there were introduced two-dimensional special series of the system {sin x sin kx}. It's shown that these series have the advantage over two-dimensional cosine Fourier series, because they have better approximation properties near the bounds of the square [0,1]2. It's given convergence speed estimate of special series partial sums to functions f (x, y) from the space of even 2n-periodic continuous functions.

Key words: special series of system {sin x sin kx}, two-dimensional series, piecewise approximation.

References

1. Sharapudinov I. I. Limit Ultraspherical Series and Their Approximation Properties. Math. Notes, 2013, vol. 94, iss. 2, pp. 281-293. DOI: 10.1134/S0001434613070274.

2. Sharapudinov I. I. Some special series in ultraspherical polynomials and their approximation properties. Izv. Math., 2014, vol. 78, no. 5, pp. 1036-1059. DOI: 10.1070/IM2014v078n05ABEH002718.

3. Dedus F. F., Mahortyh S. A., Ustinin M. N., Dedus A. F. Obobshchennyi spektral'no-analiticheskii metod obrabotki informatsionnykh massivov. Zadachi analiza izobrazhenii i raspoznavaniia obrazov [Generalized spectral and analytic method of data arrays processing. Problems of image analysis and pattern recognition]. Moscow, Mashinostroenie, 1999 (in Russian).

4. Pashkovskiy S. Numerical applications of polynomials and Tchebychev series [Vychislitel'nye primeneniia

mnogochlenov i riadov Chebysheva]. Moscow, Nauka, 1983 (in Russian).

5. Arushanyan O. B., Volchenskova N. I., Zaletkin S. F. On calculation of Chebyshev series coefficients for the solutions to ordinary differential equations. Sib. Elektron. Mat. Izv., 2011, vol. 8, pp. 273-283 (in Russian).

6. Trefethen L. N. Spectral methods in Matlab. Fhiladelphia, SIAM, 2000.

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

8. Mukundan R., Ramakrishnan K. R. Moment functions in image analysis. Theory and Applications. Singapore, World Scientific, 1998.

9. Malvar H. S. Signal processing with lapped transforms. Boston ; London, Artech House, 1992.