CC BY
32
2019 Теоретические основы прикладной дискретной математики №43
УДК 519.113.3
ДИСКРЕТНЫЕ ПЕРИОДИЧЕСКИЕ ФУНКЦИИ БЕРНУЛЛИ
М. С. Беспалов
Владимирский государственный университет имени Александра Григорьевича и Николая Григорьевича Столетовых, г. Владимир, Россия
В пространстве дискретных периодических функций с условием нормировки (сумма значений отсчётов по периоду равна нулю) рассматриваются операции циклической свертки, конечной разности и дискретного преобразования Фурье. Приводится обзор свойств дискретных периодических функций Бернулли, выделенных в качестве основного объекта изучаемой структуры. Дискретные периодические функции Бернулли положительного порядка практически идентичны конструкции специальных чисел и многочленов, введённых М. С. Беспаловым и Н. М. Коробовым.
Ключевые слова: дискретное преобразование Фурье, циклическая свертка, конечная разность, производящая функция, числа и многочлены Коробова.
DOI 10.17223/20710410/43/2
BERNOULLI'S DISCRETE PERIODIC FUNCTIONS
M. S. Bespalov
Vladimir State University named after Alexander and Nikolay Stoletovs, Vladimir, Russia
E-mail: [email protected]
This paper is a survey of known and some new properties of the discrete periodic Bernoulli functions bn(j) of order n introduced by V. N. Malozemov and viewed as elements x = x(0)x(1) ...x(N — 1) £ С CN with the normalization condition
N-1
Y x(k) = 0. It is proved that the operator A : CN ^ CN where A[x] = y = k=0
= y(0)y(1).. .y(N — 1), y(k) = x(k + 1) — x(k), is a bijection and A[bn] = bn-1. Moreover, according to Malozemov's result, the set of the discrete periodic Bernoulli functions is an infinite cyclic group relative to the cyclic convolution x * y(s) =
N-1
= S x(j)y(s — j) with a neutral element bo, and b n * bm — bn+m. It is proved j=0
that either the set of N — 1 cyclic shifts xfc^(j) = x(j — k) of any discrete periodic Bernoulli function or the set {bm, bm+1,..., bm+N-2} yields a basis of the space CN.
те
The generating function ^ bntn of a sequence of discrete periodic Bernoulli functions
n=0 N-1
is calculated. Formulas ^ sin2m(nk/N), m £ Z, for calculating the sums of even de-
k=1
grees of sinuses at equidistant nodes of a circle are found by means of these functions and the discrete Fourier transform. It has been established that a cyclic shift by 1 and the multiplication by —N transform these functions of positive order into special polynomials Pn(k), which were introduced by Bespalov and Korobov and have become popular as the Korobov polynomials of the first kind in the form Kn(x) = n!Pn(x).
We have calculated the Korobov numbers Kn = — n! ■ N ■ bn(1) up to K13 and the Korobov polynomials up to K7(x) for any array size (parameter) N.
Keywords: discrete Fourier transform, cyclic convolution,finite difference, generating function, Korobov numbers and Korobov polynomials.
Введение
Дискретные периодические функции Бернулли (далее именуются как сигналы Бернулли) введены в [1] и изучались в отдельных докладах семинара «Дискретный гармонический анализ и геометрическое моделирование» [2]. Полученные для них результаты собраны в учебном пособии [3]. Представление В-сплайнов на сетке через сигналы Бернулли предложено в [4]. Конструкция сигналов Бернулли положительного порядка практически идентична конструкции специальных чисел и многочленов, введённых в [5, 6] и ставших популярными (после публикации работ А. В. Устинова [7, 8]) в виде чисел и многочленов Коробова. Предложенное в [1] название объясняется тем, что дискретные периодические функции Бернулли служат аналогом на сетке чисел и многочленов Бернулли.
Приведённый в данной работе обзор свойств сигналов Бернулли предваряется напоминанием свойств чисел и многочленов Бернулли, описанием пространства сигналов и операций в нём. В качестве практического применения сигналов Бернулли приведён способ вычисления через них тригонометрических сумм, в п. 6 указана связь сигналов Бернулли с позднее введёнными числами и многочленами Коробова.
1. Числа и многочлены Бернулли
Числа Бернулли Bk вычисляются по рекуррентной формуле
n-1
Е сПвк = 0, (1)
k=0
к n!
где Cn = —-^уу, с начальным условием B0 = 1. В частности, B1 = —1/2, B2 = 1/6,
B3 = 0, B4 = —1/30. Все числа Бернулли с нечётными номерами, кроме B1, равны нулю.
Многочлены Бернулли определяются через их разложение по степеням x, где участвуют числа Бернулли:
n
Bn(x) = £ CnBkxn-k. (2)
к=0
1 1 3 1
В частности, B0 (x) = 1, B1(x) = x — -, B2(x) = x2 — x + -, B3(x) = x3 — -x2 + -x.
2 6 2 2
Из (2) вытекает, что Bn = Bn(0).
Приведём отдельные свойства многочленов Бернулли:
Bn(1 — x) = (—1)nBn(x) — теорема дополнения; (3)
n
Bn(x + y) = £ CnBk(y)xn k — теорема сложения аргументов; (4)
к=0 n
Bn (x) = nBn-1(x);
/ Bn (x) dx = 0 — нормировка; 0
__^ _ _ ,,„n- 1
Jn(x + 1) Bn( '
Bn(x + 1) — Bn(x) = nxn 1 — вычисление конечной разности.
В последнее время более удобным считается [8] определение чисел и многочленов Бернулли через экспоненциальную производящую функцию:
+ те в ^ В (х)
ь _ Вп _ ^ вга(х) ^га
в* - 1 п=0 П! в* - 1 п=0 П!
2. Пространство сигналов
Основной класс дискретных функций — функции, заданные на множестве целых чисел 2, то есть двусторонние последовательности. При цифровой обработке сигналов (ЦОС) рассматривают дискретные периодические (с периодом N: х(к + N _ х(к)) функции, которые для сокращения записи будем называть сигналами. Пространство (периодических) сигналов обозначим , выделяя основную часть сигнала х Е CN вида х _ (х(0) х(1) ... х^ — 1)).
Матрицу прямого дискретного преобразования Фурье (ДПФ) порядка N определим поэлементно (далее N фиксируем)
р _ ^ = )£=о, (6)
где ш _ _ exp(2п¿/N). При ЦОС исходный сигнал х преобразуется линейным оператором (6) в выходной сигнал (спектр Фурье сигнала х): у _ х • Р; перепишем его в развернутом виде
N-1
у(«) _ Е х(?, (7)
¿=о
откуда вытекает периодичность сигнала у.
Оператор обратного дискретного преобразования Фурье определяется комплексно-сопряжённой матрицей :
Р-1 _ 1 Р" _ 1 (к )N-1
Р _ N ^ ^ _ N(Ш ^=о .
Операторная запись для (7): у _ Т[х]; для восстановления: х _ Т 1[у]. Другие основные операции над х, у Е CN:
1) линейные операции — сложение сигналов и умножение на число;
2) умножение сигналов в виде покоординатного умножения
х О у _ (х(0) • у(0) х(1) • у(1) ... х^ — 1) • y(N — 1));
3) циклическая свертка
N-1
х * у (з) _ Е х(?)у(«— , з Е
¿=о
4) циклический сдвиг (вправо): ) _ х(? — к),
в частности, xN^ _ х, а xN_ х^ —операция обратного сдвига;
5) конечная разность А сигнала х Е CN:
Д[х](?) _ х(? + 1) — х(?).
Для второй операции использован известный [9] знак о операции умножения по Адамару (покоординатного умножения) векторов.
Строки матриц ^, ^ называются [10] дискретными экспоненциальными функциями (ДЭФ). Основную ДЭФ обозначим г = (1 ш ш2 ... -1 ). Остальные ДЭФ (в матричной форме) следующие:
г2° = г 0 г г3° = г2° 0 г г4° гЫ-1° гЫ° = г0°
где знак операции в показателе принято опускать. В матрице (6) ДЭФ следуют в порядке г0,гЫ-1,гЫ-2,...,г2,г, а в матрице (8)—по возрастанию степени. Вместо обозначения г0° будем использовать символ 1 для основного постоянного сигнала: 1 = (11... 1).К основным сигналам относят единичный импульс 5 = 5Ы = (1000 ... 0).
Двуместные операции о и * коммутативны, ассоциативны и дистрибутивны, в пространстве СЫ обладают нейтральными элементами 1 и 5 соответственно. Приведённые операции линейны.
Ы-1
В пространстве СЫ введём норму ||х|| = * £ |х(к)|2.
У к=0
Утверждение 1. В обозначениях X = Т [х], У = Т [у] имеем Т [х * у] = X о У, ЖТ [х о у] = X * У; выполняется равенство Парсеваля
N ||х||2 = ||Х ||2.
Это утверждение (см. доказательство в [3]) и следующие две леммы легко проверяются.
Далее рассмотрим подпространство С^ пространства СЫ, состоящее из сигналов х € СЫ, удовлетворяющих условию нормировки (аналог (5)):
N-1
£ х(к) = 0. (9)
к=0
Лемма 1. Приведённые операции, за исключением умножения сигналов о и
N 0
ДПФ, замкнуты в пространстве СЫ.
Лемма 2. В структуре (СЫ; *) единственный нейтральный элемент ^5 — — .
Лемма 3. Если Ь € СЫ, то решение х € СЫ разностного уравнения Д[х] = Ь единственно.
Доказательство. Приведём алгоритм вычисления х € СЫ. Полагаем х(0) = С, не определяя пока С. Уравнение Д[х] = Ь перепишем в виде х(2) = х(2 — 1) + Ь(2 — 1), откуда вытекает, что х(2) = С+Ь(0) + Ь(1) +... + Ь(2 — 1) для 2 от 1 до N — 1. По условию нормировки (9) составляем уравнение
N-2
ЖС + £(Ж — 2 — 1)Ь(7 ) = 0, (10)
¿=0
позволяющее однозначно вычислить С, через которое восстанавливаем все отсчёты сигнала х. ■
Итак, в классе сигналов СЫ определена обратная к операции вычисления конечной разности операция Д-1 : СЫ ^ СЫ .
Выделим возможные специальные виды сигналов класса С^. Сигнал х назовём симметричным с центром в т и обозначим х € 5(т), если х(т — к) = х(т + к) для всех к € Аналогично, если х(т — к) = х(т + к + 1), то х € $(т + 1/2). Сигнал х назовём антисимметричным с центром в т и обозначим х € (т), если х(т — к) = —х(т + к) для всех к € Аналогично, х € (т + 1/2), если имеем х(т — к) = —х(т + к + 1).
За счёт периодичности сигнала его центры повторяются через N/2: х € 5(т) влечёт х € 5(т ± N/2) и т. д.
Утверждение 2. Пусть х € С^, £ = п/2 при п € Если х € 5(£), то Д[х] € (£ —1/2), А-1 [х] € (£ + 1/2). Если х € (£), то д[х] € 5 (* — 1/2), А-1 [х] € 5 (£ + 1/2).
Доказательство, которое можно провести проверкой всех восьми случаев отдельно (хотя явно выписаны четыре случая, помним, что £ может быть целым или дробным), продемонстрируем на примере одного из более сложных случаев.
Пусть х € £(0), у = Д-1[х]. Предложим у(0) = —х(0)/2 как вариант начального значения. Формула у(з + 1) = у(з) + х(з) из леммы 3 допускает обращение в виде у(] — 1) = у(3) — х(3 — 1). По этим формулам: у(к) = х(0)/2 + х(1) + х(2) + ... + х(к — 1), у(1 — к) = —у(к) при к = 2, 3,..., ^ + 1)/2 и у(1) = х(0)/2, что влечёт у € (1/2). Осталось доказать справедливость условия нормировки (9) для построенного у, что означает верность выбора начального значения у(0).
Если N чётное, то условие (9) выполняется автоматически. Пусть N нечётное. Получили, что
у (^ + 1)/2) = х(0)/2 + х(1) + х(2) + ... + х ((N + 1)/2 — 1) (11)
и у(^ + 1)/2) = —у(1 — (N + 1)/2). По условию периодичности у(^ + 1)/2) = = у(1 — (N + 1)/2). Значит, должно быть y((N + 1)/2) = 0. Проверим — так ли это? При х € 5(0) имеем х € 5(N/2), что влечёт у(^ + 1)/2) = х^ — 1) + ... + +х(^ + 1)/2) + х(^/2 из (11). Сложим эти два представления: 2у^ +1)/2) =
N -1
= х(к) = 0 по условию х € CN. к=0
Ввиду периодичности считаем, что доказан случай
х € 5(т) ^ Д-1[х] € (т +1/2).
3. Сигналы Бернулли
Основным представителем класса CN служат дискретные периодические функции Бернулли [3, с. 171] порядка в, координаты которых определяются равенством
1 N-1 /2ПЛ
Мз) = N £ (^ — 1)- ■ , в, 3 € Z, ^ = ех^ ■—\ . (12)
Будем называть их короче — сигналами Бернулли порядка в (периода N, где N считаем фиксированным). В частности, начальный сигнал Бернулли (для наглядности добавили разделительные запятые)
Ьо = 5 — N1 = ^ — 1, —1, —1,..., —1).
По свойствам корней Ж-й степени из единицы
N-1 N-1
£ ш^' = 0 при всех 1 ^ 2 < Ж, £ ш^ = N (13)
к=0 к=0
проверяется указанный вид Ь0, принадлежность всех ДЭФ (кроме 1) пространству С,^ и следующая лемма.
Лемма 4. Для сигналов Бернулли любого целого порядка Ь5 € С,^. Лемма 5. Для сигналов Бернулли Ь5 = Д[Ь5+1] и Ь5+1 = Д-1[Ь5]. Доказательство. Вычислим координаты сигнала
1 N-1 1 N-1
Д[Ь*+1 ](2) = 1 £ (шк — 1)-(5+1)шк(^+1) — - £ (шк — =
^ к=1 ^ к=1 1 N-1
= 1 Е (шк — 1)-(*+1)шк^ (шк — 1) = Ьв(2). ^ к=1
По лемме 3 получаем второе утверждение леммы 5. ■
Замечание 1. В качестве альтернативного определения сигналов Бернулли можно предложить формулу Ь0 = 5 — — 1 для сигнала Бернулли нулевого порядка и соотношения леммы 5.
Следствие 1. Для сдвига сигнала Бернулли верно = Ь5 + Ь5-1. Доказательство. В координатах доказанная формула Ь5-1 = Д[Ь8] выглядит так: Ьв-1 (2) = Ьв(2 + 1) — Ь5(2) или ЬД2 + 1) = Ь5(2) + Ьв_1(2). ■
Следствие 2. Отсчёты всех сигналов Бернулли вещественные. Более того, отсчёты всех сигналов Бернулли отрицательного порядка и координаты векторов Nк+1 • Ьк при неотрицательных к — целые числа.
Доказательство. Имеем Ь-1 = Д[Ь0] = (—1 00... 00 1). Так как отсчёты Ь-1 целые, то по лемме 5 получаем, что отсчёты каждого следующего сигнала Бернулли отрицательного порядка целые.
Координаты вектора N • Ь0 = N5 — 1 суть целые числа. Для сигналов Бернулли следующих (больших) порядков, согласно формуле Ь5+1 = Д-1[Ь5], вычисляем координаты по алгоритму леммы 3, решая при этом уравнение (10) для С и накапливая множитель N в знаменателе. ■
Утверждение 3. При к ^ N имеет место следующее представление для сдви-
к
гов сигнала Бернулли: Ьк+к = £ С^Ь^+т. В частности, при к = N получаем аналог
т=0
формулы (1), который справедлив лишь для N, равного размеру сетки:
N-1
£ стЬ5+т = 0. (14)
т=0
Доказательство. Применив формулу Ь^+1 = Ь5+1 + Ь5 следствия 1 трижды, получаем базу индукции
Ь2+2 = Ь«+2 + Ь^+1 = Ь«+2 + Ьз+1 + Ьз+1 + Ь*.
Аналогично доказываем индуктивный переход от к — 1 к случаю к:
ьк^ _ Ь(к-1)^ , 7 (к-1)^ _ 1 ст ь + 1 ст ь _
Ь«+к _ Ь«+1+к-1 + Ь«+к-1 _ Ск-1Ь«+1+т + Ск-1Ь«+т _
т=0 т=0
к к-1 к _ Е Ск-А+. + Е Ск-А+. _ Е Ск Ь«+7,
.7=1 .=0 .=0
где применили известное свойство биномиальных коэффициентов
с£1 + ск-1 _ С. (15)
Формула (14) вытекает из равенства ЬN^ _ Ь для Ь € . ■
Замечание 2. Формулу (14) при небольших N можно использовать как следующий рекуррентный способ вычисления значений сигналов Бернулли:
1 N -1
Ьз _ — ^ £ ^/^-к. ^ к=1
В частности, Ь5 _ — -Ь8-1 при N _ 2, Ь5 _ —Ь8-1 — -Ь3-2 при N _ 3 и т. д.
2 3
Следствие 3 (теорема сложения аргументов, аналог (4)). При к ^ 0 верно
Ьи(з + к)_ Е Ск Ьп-г(з). 1=0
Замечание 3. В следствии 2 приведено утверждение о целочисленности отсчётов сигналов Nк+1Ьк, которое довольно грубое. Анализ биномиальных коэффициентов в соотношении линейной зависимости (14) позволяет получить следующее утверждение:
Если N _ р^1 РГ .. .-РТ, где р., 2 _ 1,... , к, — простые числа, п ^ 0, то отсчёты сигналов
[ п ] [ п ] [ " ]
^Г-1 ] рР2-1 ] ...ркРк-1 ]Ьп целочисленные. Доказательство не приводим из-за его громоздкости.
Лемма 6. Дискретное преобразование Фурье сигналов Бернулли имеет вид
Т[Ьп](0) _ 0, Т[Ьп](2) _ (ш7 — 1)-п при 2 _ 1, 2,..., N — 1,
что короче записывается при п € N двумя формулами:
Т[Ь-п] _ (г — 1)п и условно Т[Ьп] _ (г — 1)-п
Доказательство. По лемме 4 имеем Т[Ьп](0) _ 0. Для остальных отсчётов
N-1 1 N-1 1 N-1 N-1
Т[Ьп](2) _ Е 1 Е (шк — 1)-п^ш-7 _ 1 Е (шк — 1)-п Е ^(к-7> «=0 ^ к=1 ^ к=1 «=0
Применяя (13), получаем Т[Ьп](2) _ (ш7 — 1)-п. ■
Следующее свойство сигналов Бернулли первым, по-видимому, отметил В. Н. Мало-земов.
Теорема 1. Циклическая свертка сигналов Бернулли равна соответствующему сигналу Бернулли:
Ьп * Ьт Ьп+т.
Доказательство. Первый способ — проверяем непосредственно:
N-1 1 N-1 N-1
(Ьп * Ьт)(з) = Е -ш Е (^ — Е (^ — 1)-т^'-в)г =
3=0 ^ '=1 г=1
1 N-1 N-1 N-1
= ^ Е (^ — 1)-п Е (^ — 1)-т^7 Е
^ '=1 1=1 3=0
Продолжаем по формуле (13):
1 N-1
(Ьп * Ьт)(3) = — Е — 1)-П(^ — = Ьп+тЗ).
Второй способ — через ДПФ. По утверждению 1 и лемме 6 с учётом того, что нулевой отсчёт всех спектров Фурье равен нулю, получаем
Т[Ьп * Ьт] = Т[Ьп] ◦ Т[Ьт] = (г — 1)-П о (г — 1)-т = Т[Ьп+т].
Теорема доказана. ■
Есть третий (более длинный) способ доказательства через итерации операторов А и А-1. Его укажем в виде следствия теоремы.
Следствие 4. Для любого т € Z конечная разность всех порядков и её обращение для х € С,^ представима как
Дт[х] = Ь-т * х.
Эту формулу можно переписать в виде х = Дт[х] * Ьт, полученным в [1] и применённым в [4] для анализа В-сплайнов.
Следствие 5. Сигнал Бернулли любого порядка есть соответствующая степень образующего сигнала Ь1 :
Ьп = Ьп*, где п € Z. Лемма 7 (аналог теоремы дополнения (3)). Для сигналов Бернулли верно
Ьт(т — з ) = (—1)тЬт(з). Доказательство. По определению (12) и условию UJN = 1 имеем
1 N-1 1 N-1
Ьт(т — з) = — Е (^ — 1)-т^(т-л = - Е (1 — ^=
1 N-1 1 N-1
= Тг Е (1 — ^)-т^= I1 Е (1 — ^= (—1)тЬт(з).
Лемма доказана. ■
Из утверждения 2 вытекает
Утверждение 4 (свойство симметрии, антисимметрии). Сигнал Бернулли чётного порядка симметричен, нечётного — антисимметричен:
Ь2т € 5(т), Ь2т+1 € (т + 1/2).
N — 1 N — 1
Следствие 6. Для отдельных отсчётов верно Ь1(0) =--^ и Ь1(1) = ^ .
, ^ п N — 1 — 2к Для остальных отсчётов Ь1(1 + к) =--.
N — 1
Доказательство. Имеем Ь0 € 5(0), Ь1 € (1/2), Ь0(0) = ———; обозначим N — 1
2С = ———. Тогда Ь1(0) = —С и Ь1(1) = С; далее по лемме 3. ■
При 1 ^ к < N следующая формула для даёт явный вид начальных отрицательных сигналов Бернулли в виде сдвинутого набора биномиальных коэффициентов с чередующимися знаками. Её частные случаи:
Ь-1 = (—1 0 0 ... 0 0 1), Ь-2 = (1 0 0 ... 0 0 1 — 2), Ь-3 = (—10 ... 0 1 — 3 3), Ь-4 = (10 ... 0 1 — 4 6 — 4) и т.д.
Утверждение 5. Для сигналов Бернулли отрицательного порядка имеет место
Ь-* = Е (—1)' С' 5(7; (16)
.7=0
= 4' (2к — 1)!! (17) ||Ь-'1 = (2к)!! . (17)
Доказательство. Следствие 1 влечёт базу индукции для (16): Ь-1 = Ь^ — Ь0 = 5^ — 5 — -1(1^ — 1) = 5^ — 5, Ь-2 = ЬГ1 — Ь-1 = 52^ — 25^ + 5.
Индуктивный переход (во второй сумме заменяем в = з +1, переобозначаем символом з и применяем (15)):
Ь-' = ЬГ'+1 — Ь-'+1 = Е(—1). С^7 — Е(—1). С^-7 =
. =0 . =0
= 5^ + Е(—1)7 7 + С'-!)^-7^ + (—1)' 5 = Е (—1)7 С 5('-7)^.
7=1 7=0
Из (16) вытекает приведённый явный вид сигналов и квадрат нормы:
1|Ь-'I2 = Е С)2.
7=0
По тождеству Вандермонда [11, с. 27] 7 — п'
Л ^ = С = М = (2к)!!(2к — 1)!! = 2'(2к — 1)!! = 4'(2к — 1)! 7=0(к!)2 (к!)(к!) к! (2к)!!
7=0
Утверждение доказано
Следствие 7. Набор из N — 1 сигналов Бернулли подряд идущих порядков составляет базис пространства С^.
Доказательство. Подпространство С^ ^мерного пространства См имеет размерность N — 1. Если отсчёты с номером 0 временно удалим, как допускающие восстановление по свойству нормировки (9), то заметим, что векторы 6_1,6_2,..., +1 линейно независимы, как обладающие диагональной структурой. Значит, они составляют базис пространства С^. По формуле (14) линейной зависимости этот базис можно сдвигать на одну позицию в любую сторону. ■
Замечание 4. Очевидно, что для любого х Е С^ полный набор его сдвигов ли-
N _1
нейно зависим: У] = 0. Зато набор из ^ — 1)-го сдвига любого сигнала Бернулли к=0
составляет базис пространства С^^. Для сигнала Бернулли порядка — 1 это легко проверяется, а далее доказывается по лемме 3.
Рассмотрим простейшие частные случаи сигналов Бернулли. Случай N = 2. Пространство СО одномерно и при всех целых в
• 1)1 (0 — 1)= (— 1ьо.
' 2'+^ ' \ 2/ Случай N = 3. Пространство СО двумерно с базисом
Ьо = 1(2 — 1 — 1) и Ь_1 = (—1 0 1)
и соотношением 6к+2 = ——6^.
3
Случай N = 4. Хотя пространство С трёхмерно, сигналы Бернулли группируются не по тройкам, а по парам (чётные и нечётные порядки), а именно (отделили координаты запятыми):
62' = 2'+' (2'+1 + 1, —1, —2'+1 + 1, — 1Г , ( 1)'_1
б2'_1 = (2' + 1, 2' — 1, —2' + 1, —2' — 1)'^ .
Для остальных N общих формул представления всех сигналов Бернулли нет. По теореме 1 в кольце (С^^;+, *) легко вычисляется обратный элемент к сигналам Бернулли. Для произвольного обратимого х Е С^^ есть две стандартные схемы вычисления обратного, то есть такого у Е С^^, что х * у = 60. Один приём «в лоб» методом неопределённых коэффициентов — составление и решение системы ^го порядка. С использованием циркулянтной матрицы С(х) для сигнала х (её определение см. в [9, с. 40]) уравнение записывается в матричном виде у • С(х) = 60.
Второй приём покажем на примере сигнала х = (1, 0, 3, —4). Вычислим ДПФ от него: X = Т[х] = (0, —2 — 4г, 8, —2 + 4г). Найдём У, такой, что X о У = 1 — 8 =
= (0,1,1,1): У = (0, —1 + 2г,1, —1 — 2 ). Через обратное ДПФ от У найдём искомый
10 8 10 1
У =-(—3, —21,13,11).
у 1601 ' ..
Предложим третий приём. Произвольный х Е С^^ легко раскладывается по базису
{6_1, 6_2,... , +1}. В случае х = (1,0, 3, —4) имеем х = 26_1 + 36_2. Искомый у Е С$,
такой, что x * y = bo, будем искать в сдвинутом базисе в виде y = abi + ßbo + 7^-1-По линейности циклической свертки и теореме 1 имеем
x * y = 2abo + (3a + 2ß)b- + (3ß + 27)b-2 + 37b-a.
Здесь на помощь приходит соотношение (14) линейной зависимости, которое добавим с некоторым неопределённым коэффициентом s:
x * y = bo + s(4bo + 6b_i + 4b_2 + Ь_з). Приравняв коэффициенты, получаем систему, которая решается последовательно:
s Л 10 34 27 17 3 9
Y = -, ß = — s, a = — s, s = —— —> a =--, ß = — -, y = — —.
' 3> 9 ' 27 ' 40 20 4' ' 40
При N ^ 5 первые два приёма слишком сложны для использования без привлечения компьютера, а с помощью третьего можно решать и вручную.
4. Производящая функция сигналов Бернулли
Производящей функцией последовательности a = {ak}£=o называется [11]
те
p(x) = Е afcxk.
fc=o
Введём производящие функции для N отсчётов последовательности сигналов Бернулли (12) неотрицательного порядка:
те
№(x) = E bs(k)xs, k = 0,1,..., N - 1.
s=o
те
Теорема 2. Производящая вектор-функция Ф^) = £ bsxs для всех отсчётов сиг-
s=o
налов Бернулли неотрицательного порядка равна
Ф(x) = S + p(x; N) ((1 + x)N-1,1, (1 + x), (1 + x)2,..., (1 + x)N_2) ,
где ^(x; N )=1 - (1x+ x)N.
Доказательство. По формуле bs(k) = bs(k — 1) + bs-1(k — 1) из следствия 1 имеем
тете
(x) = bo(k) + £ bs(k)xs = bo(k) + E(bs(k — 1)xs + bs-i(k — 1)xs).
s=1 s=1
При 2 ^ k ^ N — 1 получим
1 те те
^fc(x) = —— + E bs(k — 1)xs + x £ bs(k — 1)xs = <£fc-1(x) + x^fc-1(x). N s=1 s=o
Небольшое отличие при k = 1 и k = N :
те
<^(x) = bo(1) + E(bs(0)xs + bs-1(0)xs) = bo(1) — bo(0) + (1 + x)^o(x),
s=1
(x) = bo(N ) — bo(N — 1) + (1 + x)^N-1 (x).
Получили формулы
<£1(х) = —1 + (1 + х)^0(х), ^0(х) = ^ (х) = 1 + (1 + х)^-1(х). Выразим все рассмотренные производящие функции через одну:
^(х) = (1+ х)"-1^ (х) при 1 ^ к ^ N — 1,
Ых) = 1 + ^\(х), ^0(х) = 1 + (1+ x)N-1^(х).
х + 1
Приравняв две последние, получим
<^(х; N) = ^(х) х
1 — (1 + x)N' Теорема доказана. ■
Следствие 8. Отсчёт сигналов Бернулли с номером 1 порядков в ^ 2 как функция от размера сетки N есть многочлен нечётного порядка; гарантируется, что среди его корней присутствуют корни N = ±1.
х
1 — (1 + х);
Доказательство. Для ((x; N) = ^i(x) = --¡r¡—;—вычислим
x
((x;-N) = 1 - (1 + x)-n = x - ((x;N)'
что влечёт первое утверждение. Более того, из доказательства следствия 2 вытекает наличие N в знаменателях всех отсчётов сигналов Бернулли положительного порядка.
Рассмотрим величины — Nbs(1) как чётные многочлены степени s от переменной N. Из условия ((x; 1) = —1 вытекает, что при s ^ 2 числа N = ±1 служат корнями этого многочлена. ■
5. Тригонометрические суммы
Через производящую функцию теоремы 2, взятую в виде —N((x; N), в [5, 12] вычисляются суммы чётных отрицательных степеней синусов в равноотстоящих узлах окружности. В. Н. Малоземов в [2] предложил более красивое изложение этих результатов с распространением на положительные степени. Эти суммы имеют практическое применение [13].
Теорема 3. Суммы чётных степеней синусов в равноотстоящих узлах выражаются через нормы сигналов Бернулли (s > 0):
N—1 nk N N-1 1
Е sin2s N = N ИМ2, I— ■ 2s(1 k/N) = 4sNllfesM2. k=i N 4s fc=i sin (nk/N)
В частности, при s < N имеем
n—i (2s - 1)!!
ЕГ N = N w"1- (18)
Сумма E . 2 — /A7N при натуральном s равна чётному многочлену от N степени 2s,
N —1 1
V _-
¿1 sin2s(nk/N)
содержащему множитель ^2 — 1). В качестве формулы для вычисления этих многочленов предлагается
V 1 = ((* — ^ — Ц«**^ — 1) \ (19)
¿1 8ш2*(пк^) (в!)2 V (^ — 1)(^ — !)(** — 1) У . ( 9)
Доказательство. По равенству Парсеваля (утверждение 1) N||68||2 = [&«]||2. Правую часть вычислим согласно результатам леммы 6:
N — 1 N— 1
l|F[Ml2 = ||(r - i)-sII2 = Е l^fc - i|—2s = Е l^fc/2 - k/2|—2s ■ M—.
fc=i fc=i
По формуле Эйлера = cos ^ + i sin ^ получим
N1
IIF [bs]|2 = E fc=1
2sin N
nk
Далее считаем s > 0 и запишем две формулы:
N-1 nk N—1
N ИМ2 = 4s E sin2' Пк, N ||bs||2 = 4—' N- sin-2' Пк.
fc=i N fc=i N
Для первой из них по формуле (17) при s < N получаем (18). Сигналы Бернулли положительного порядка восстанавливаются через производящую вектор-функцию теоремы 2: = ^ (Ф(ж))(а) |x=0, а для отдельных отсчётов
bs(k) = i lim (-— (un = — lim | ---(—)
'v ; s! du'^feV 7 s! dz' 1 - zN
При вычислении квадрата нормы вместо нулевого отсчёта берем отсчёт с номером N:
||b ||2 = (- V2 f lim -^L ( zk-1(z - 1)tk-1(t - -)
И"*|| V«V ¿iz-1'^1 dz' Sí* V (zN - 1)(íN - 1)
что равносильно (19).
Так как при замене N на - N в правой части (19) ничего не меняется, эта функция чётная. К тому же эта функция является многочленом, что нетрудно проверить (в [5] это доказано для очень близкой конструкции). Один из корней этого многочлена N =1, так как при N =1 правая часть (19) обращается в нуль. Из чётности многочлена вытекает, что и N = -1 есть корень. ■
Приведём частные случаи результатов теоремы 3:
N-1 ^ АТ N-1 4 nk 3N
„ АТ \ о __„„„ АТ \ о
fc=1 N 2 fc=1
2 nk N
' N " = У
nk 5N
N = 16
N 8 к
N-1 «nk 5N . N-1 8 nk 35N
i — =-при N ^ 5.
k=1 n 16 fc=1 N 128
Из курса математического анализа известно [14, с. 434], что
Г Г( ^s+i) (2s — 1)!!
X sin xdx = ^^ =П W1' (20)
Сравнивая (20) с формулой (18), замечаем, что для интеграла /0 sin2' x dx обе формулы прямоугольников (при Ax = n/N) точны при достаточном (N ^ s + 1) числе узлов:
N — 1 N Гп
Ax £ sin2'(kAx) = Ax£ sin2'(kAx) = sin2' x dx. fc=0 fc=1 Jo
Е . 2з/—, /А7Ч по произвольным заданным N и в.
Верны следующие частные случаи:
^ 1 _ N2 — 1 ^ 1 _ (N2 — 1)(N2 + 11)
¿1 яп2(пк^) = 3 , ,£1 яп4(пк^) = 45 ,
^ 1 _ (N2 — 1)(2^ + 23N2 + 191) '=1 ¡Ш^к/^) = 945 ,
^ 1 _ ^2 — 1)(3N6 + 43N4 + 337N2 + 2497) ¡Ь^к/^) = 14175 ,
^ 1 _ (N2 — 1)(2N8 + 35N6 + 321N4 + 2125N2 + 14797) ¡¡п0^/^) = 93555 .
Из этого набора формул легко вычисляются сумммы чётных степеней котангенсов в равноотстоящих узлах, приведённые в [15] и востребованные в механике.
По алгоритму леммы 3 и формуле теоремы 3 составлена программа [16] вычисления
N-1 1
¿=1 8Ш2в(7к/^)
6. Связь с числами и многочленами Коробова
При построении производящей функции в теореме 2 более удобный вид получается при рассмотрении обратного сдвига. Поэтому далее рассмотрим следующую последовательность дискретных периодических функций (называемых сигналами):
—N ■ Ь£Г, —N ■ Ь|", —N ■ Ьт,..., —N ■ ЬТ,...
Будем использовать предложенные в [6] названия и обозначения:
— специальные числа Р5 для начальных отсчётов этих сигналов:
Р = —N ■ ЬТ(0) = —N ■ Ьв(1);
— специальные многочлены Р8(к) для остальных отсчётов этих сигналов (за одним исключением: последняя координата начального сигнала выпадает из общей конструкции, Р0(х) = 1):
Рв(к) = —N ■ Ь5(к + 1), к = 0,1,..., N — 1.
Специальные числа Р5 и специальные многочлены Р8(х), а также производящие функции для них
те те
= Е РХ, Ф(и, х) = Е Р.(х)м"
«=0 «=0
введены и вычислены в [5, 6]. В [5] обозначения другие. После естественного видоизменения эти специальные числа и многочлены приобрели вид чисел Коробова и многочленов Коробова первого рода [8]
К = в!Р5, Кв(х) = в!Рв(х)
с теми же экспоненциальными (термин из [11]) производящими функциями
те К те К (х)
^(«) = ЕК < ф(«,х)=ЕКпх)и5. «=0 в! «=0 в!
Из теоремы 2 вытекает
Утверждение 6. Экспоненциальные производящие функции чисел и многочленов Коробова первого рода имеют следующий вид:
^ ^(м + 1)х
^(м) = ^(м; N) = --—^--, Ф(м,х) = Ф(м,х; N) =
(и + 1^ — 1' К У ' У ' ' ' (и + 1)N — 1'
Основной результат работы [5] сформулируем в принятых обозначениях, используя понятие обобщённой степени к(1) = к(к — 1)... (к — в + 1) .
Теорема 4. Решение бесконечной (так как в Е М) системы разностных уравнений
АР ] = Р5-1
N-1
с условием нормировки Е Р?(к) = 0 и начальным условием Р0(к) = 1 имеет следую-
'=0
щий вид:
к(') к('-1) к(1) Р'(к) = Р0 ^ + Р1-7-—тгг + ... + Р'-1 V + Р-. (21)
в! (в — 1)! 1!
Здесь специальные числа Р8 имеют вид многочленов степени в от параметра N со старшим коэффициентом В8/в!, где В — числа Бернулли. Начиная с в = 2, это чётные многочлены (то есть многочлены от N2), для которых числа 1 и —1 служат корнями.
Рекуррентная формула для вычисления чисел Р5 как многочленов от N имеет следующий вид:
N ('+1) N(5) N(2)
Р0 (7+ЩТ + Р1 7Т +... + Р'-1 + = 0. (22)
Доказательство. По условию нормировки из соображений симметрии имеем N — 1 N — 1 Р1(к) = к---—, что можно представить в виде Р1(к) = Р0к(1) + Р1, где Р1 =---—.
Это база индукции для доказательства (21) и (22). Предположим, что
к('-1) к('-2) к(1) Р-1<к) = Р0 (7ГЩ7 + Р1 (Т—^ + ■ ■ ■ + Р'-2 - + Р-1-
Так как А[к(п)] = пк(п-1), решение уравнения А[Р8] = Р5-1, согласно лемме 3, имеет вид (21), где Р5 пока не найдено. Ищем его методом прогонки:
Р5(0) = Р5, Р5(1) = Р5 + Р5-1, Р5(2) = Р, + 2Р5-1 + Р5-2,...
Общий вид формулы для 0 ^ к ^ N — 1:
Р-(к) = Е С'Р«-7.
7=0
N-1 N-1 '
По условию нормировки Е Р?(к) = Е Е С'Р- = 0. Отделим старшее слагаемое и
'=0 '=0 7=0 поменяем пределы суммирования:
N-1 / N-1 \
+ Е- Е- СМ Р'-7 = 0.
7=1 \ г=7 /
N-1
Согласно [11, задача 1.2], получили £ Сд^Р- + N^5 = 0. Запись этого выражения
.7 = 1
в виде многочлена от N приводит к (22).
Теперь рассмотрим специальные числа Р5 как многочлены от N степени 5, где старший коэффициент при N5 может оказаться нулём. Этот коэффициент найдём анализом старших слагаемых.
База индукции: Р0 = -0 = 1, Р1 = + О^1-1), где -0 = 1, — = -1/2-числа Бернулли.
—
Индуктивное предположение: Р- = —-п- + O(N--1) при всех к от 1 до п — 1. Так как
к!
формулу (22) анализируем для всех N, отмечаем, что с ростом N она обрывается на п:
1 п
Р =__1 V —_Р
N 7=1 (г + 1)!Рп-7 •
Следим только за старшим коэффициентом:
Р" = — I ((п-—Т)^ + = Ц5 (Т?—.)
( — 1 п п! \ в
"Г £ г + пй-л?Nп + О(^-1) = —■ + О(^-1)
п! 7=1 (г + 1)!(п — г)! у п!
по свойству (1).
Чётность Рп как многочлена от N и наличие множителя (^ — 1) (при п ^ 2) доказаны в следствии 8. ■
В [12] отмечено, что решение данной системы разностных уравнений уже получено, найдены производящие функции для них и приведена формула (19) для тригонометрических сумм.
В [6] специальные числа Р5 вводятся формулой (22) с начальным условием Р0 = 1, а специальные многочлены Р5(ж), для которых Р5(к) = — N6^ + 1) на сетке, — формулой (21) с начальным условием Р0(ж) = 1, без мотивировки происхождения этих формул. В [8] предложен подход к определению Р5 и Р5(ж) через производящие функции.
Утверждение 7. Три следующих способа определения специальных чисел Р5 и специальных многочленов Р5(ж) (а следовательно, также чисел и многочленов Коробова первого рода) эквивалентны:
1) в виде решения системы разностных уравнений с нормировкой (9);
2) формулами (21), (22);
^ ^(м + 1)х
3) через производящие функции ^(м) = + — 1 и У(м,ж) = + _ 1
соответственно.
Переходы 1 ^ 2, 1 ^ 3 доказаны в [5] (см. также теорему 4 и утверждение 6), переход 2 ^ 3 — в [6], а переход 3 ^ 1 — в [8].
Замечание 5. Существует четвёртый (рекуррентный) способ определения специальных многочленов Рп(к) (от целочисленного аргумента к на сетке) с помощью формулы (14). Этот способ аналогичен формуле (22) для вычисления специальных чисел. Однако для него нужно слегка поправить определение Р0(к): оставляем неизменным Р0(к) = 1 при 0 ^ к ^ N — 2 и изменяем один отсчёт Р0 ^ — 1) = 1 — N .За счёт
этой поправки, которая не влияет на вычисление последующих многочленов Рп(к), добились того, что и для начального многочлена верно условие нормировки
N-1
Е Рп(к) = 0. (23)
'=0
Введение поправки позволяет определить и вычислять, применяя операцию А, «специальные многочлены отрицательного порядка» по аналогии с сигналами Бернулли. В качестве начальных данных предлагаемого четвёртого способа достаточно вычислить N — 2 таких «многочленов», начиная от Р-1 до -2). Формулы для вычисления их отсчётов, удовлетворяющих условию (23): Р-'(3) = 0 при 0 ^ ] ^ N — к — 2, Р-'^ — к — 1 + в) = (—^^С' при 0 ^ в ^ к.
Наличие величин, которые названы «многочленами», позволяет запустить процесс вычисления специальных многочленов по формуле
1 N-1
Ps(з) = — -1 Е С^Ч-'(з) = 0, в = 1, 2, 3,... (24)
^ '=1
Приведём несколько начальных чисел Коробова, вычисленных по формуле (22) и соотношению К, = в!Ps (отметим, что в [6-8] вычисления выполнены до третьего порядка):
N — 1 N2 — 1 N2 — 1 1
К = — --1, К = --1, Кз = — --1, К = — — (N2 — 1)(^ — 19),
1 2 ' 2 6 ' 3 4 ' 4 3Г Л
К = (N2 — ^ — 9), К6 = -Ш(N2 — 1)(2N4 — 145^ + 863),
4 84
5
К7 = — — (N2 — 1)(^ — 25)(2N2 — 11).
Следующие числа Коробова вычислены методом неопределённых коэффициентов с использованием результатов теоремы 4:
К8 = — ^^2 — 1)^6 — 497^ + 9247N2 — 33953),
7
К9 = — (N2 — 1)^2 — 49)(3N4 — 50N2 + 167),
20
K10 =-(N2 - 1)(10N8 - 2993N6 + 114597N4 - 1184767N2 + 3250433),
660
15
Kn = - — (N2 - 1)(N2 - 9)(N2 - 81)(2N4 - 49N2 + 173), 8
N 2 - 1
K12 =--(1382N10 - 653818N8 + 42418211N6 - 845983589N4 + 6117468907N2-
5460
-13695779093),
K13 = ^^(N2 - 1)(N2 - 121)(1382N8 - 77096N6 + 1336965N4 - 8756954N2 + 18382103). 840
Из наблюдения за видом чисел Коробова вытекает Гипотеза 1.
lim K2n+1 = 1 - 4n2 N ^те K2n 2
Если гипотеза верна, то нам известны старшие ненулевые коэффициенты всех чисел Коробова. В качестве упражнения можно попробовать методом неопредёленных коэффициентов продолжить вычисления чисел Коробова, общий вид которых (при верности гипотезы) следующий:
7 455
К14 = -(^ — 1)Т12 К15 = — — (^ — 1)(N2 — 169)Тю(^, 6 4
К16 = — 1)Tl4(N), К17 = -7-^2 — 1)(N2 — 9)^2 — 2-)(N2 — 22-)Т8(^.
510 4
Здесь ) —чётный многочлен степени п со старшим коэффициентом 1, коэффициенты которого и следует вычислить. Для применения метода неопределённых коэффициентов предлагается таблица значений отдельных специальных чисел Рп = Кп/п! в зависимости от размера сетки N.
п N
2 3 4 5 6 7
14 1/214 2/2187 —255/16384 8/125 71698571/35831808 445/49
15 —1/215 —1/2187 255/32768 —21/125 —71698571/71663616 357/49
16 1/216 1/6561 1/65536 76/625 —429909599/429981696 —948/49
17 —1/217 0 —513/131072 0 262143/131072 474/49
Для вычисления К16 добавим недостающее значение Р16 = —6684671/216 при N = 8. Приведем также вид начальных многочленов Коробова:
N — 1 о N2 — 1
^^ 1 ^ ~
К0(ж) = 1, К1(ж) = ж — ^ ^ , К2(ж) = ж2 — ^ +
6
^ / ч 3 3(N +1) 2 N(N + 3) N2 — 1 Кз(ж) = ж3 — 1 2 ; ж2 + -ж--4—,
К4(ж) = ж4 — (2N + 4)ж3 + (N2 + 6N + 4)ж2 — 2N^ + 2)ж — (N2 — 1) (^2 — 19),
30
5 5 5
К5(ж) = ж5 — -(N + 3)ж4 + 2 + 9N + 11)ж3 — -(3N2 + 1Ш + 6)ж2 — 2 3 2
— — --N2 — 90Мж + 2 — 1)(^ — 9),
64 5
Кб (ж) = ж6 — 3^ + 4)ж5 + 2 + 12N + 21)ж4 — -(4N2 + 2Ш + 20)ж3+
+ 1(—N + 10-N2 + 300N + 144)ж2 + (2N4 — -0N2 — 72N )ж + Кб,
7 7 35
К7(ж) = ж7 — - (N + -)ж6 + 2 + 1-N + 34)ж5 — — (-N2 + 34N + 4-)ж4+
77
+ -(—N4 + 170N2 + 67^ + -48)ж3 + -(-^ — 22^2 — -48N — 240)ж2+ 64
+ 6 — 119^ + 1918N2 + 2-20Мж + К7. 6
По свойству многочленов Коробова первого рода ДК5(ж) = зК5-1(ж) и числам Коробова К5 = К5(0) в качестве начальных условий вычисляется следующий общий вид старших коэффициентов многочленов Коробова первого рода:
К (ж) = т« _ n(N + п — 2) - + С2 2N2 + 6(п — 2)N + 3п2 — 13п +12 т„-2_
Кп(ж) — ж ж | V.-'п л гл ж
2 ■ п 12 12
У
3 (2п — 4)N2 + (3п2 — 13п + 12)N + (п3 — 7п2 + 14п — 8) п-3
сСп ж I • • •
При фиксированном N численное значение чисел и многочленов Коробова на сетке удобнее вычислять не по этим формулам, а с помощью леммы 3 или формулы (24). В [8] предложен метод вывода многочленов Бернулли из многочленов Коробова
lim = ад.
Заключение
Дискретные периодические функции Бернулли служат естественным и достаточно удобным аппаратом при более подробном изложении теории дискретного преобразования Фурье. Поэтому не случайно на них практически одновременно и независимо друг от друга вышли проф. В. Н. Малоземов и автор работы. В рамках теории чисел дискретные периодические функции Бернулли представлены в виде многочленов Коробова и активно исследуются. Например, работа [17] в течение 2017г. процитирована 5 раз. Кроме возможных приложений в технике, их востребованность обусловлена ещё и тем, что числа и многочлены Коробова служат одним из основных примеров для изложения теневого исчисления, развитого в работах Рота и его последователей [18]. В рамках теневого исчисления изучаются [8] также числа и многочлены Коробова второго рода, которые в данной работе не затрагиваются. Близким объектом служат многочлены Карлица.
ЛИТЕРАТУРА
1. БерМ.Г., Малоземов В. Н. Наилучшие формулы для приближенного вычисления дискретного преобразования Фурье // Ж. вычислит. матем. и матем. физики. 1992. Т. 32. №11. С. 1709-1719.
2. Избранные главы дискретного гармонического анализа и геометрического моделирования. Ч.2 / под ред. В.Н. Малоземова. СПб.: Изд-во ВВМ, 2014. 605 c.
3. Малоземов В. Н., Машарский С. М. Основы дискретного гармонического анализа. СПб.: Лань, 2012. 304с.
4. Малоземов В. Н., Певный А. Б. Дискретные периодические сплайны и их вычислительные применения // Ж. вычислит. матем. и матем. физики. 1998. Т. 38. №8. С. 1235-1246.
5. Беспалов М. С. Представление для сумм четных отрицательных степеней синусов в равноотстоящих узлах // Изв. вузов. Математика. 1996. Т. 8 (411). С. 6-12.
6. Коробов Н. М. Специальные полиномы и их приложения // Диофантовы приближения. Математические записки. 1996. Т. 2. С. 77-89.
7. Устинов А. В. О формулах суммирования и интерполяции // Чебышевский сб. 2001. Т. 1. №1. С. 52-71.
8. Устинов А. В. Полиномы Коробова и теневой анализ // Чебышевский сб. 2003. Т. 4. №4(8). С. 137-152.
9. ХорнР., Джонсон Ч. Матричный анализ. М.: Мир, 1989. 655 с.
10. Т'рахтман А. М., Трахтман В. А. Основы теории дискретных сигналов на конечных интервалах. М.: Сов. радио, 1975. 208с.
11. Комбинаторный анализ. Задачи и упражнения / под ред. К. А. Рыбникова. М.: Наука, 1982. 368с.
12. Беспалов М. С. О производящих функциях для некоторых тригонометрических сумм // Научные исследования института — техническому и культурному прогрессу. Материалы XXV научн. конф. ВПИ. Ч. 1. Владимир: ВПИ, 1990. С. 40.
13. Гузев М. А., Устинов А. В. Механические характеристики модели молекулярной динамики и полиномы Коробова // Дальневост. матем. журн. 2016. Т. 16. №2. С. 39-43.
14. Зорич В. А. Математический анализ. Ч. II. М.: Наука, 1984. 640с.
15. Беспалов М. С. Тригонометрические суммы для задач молекулярной динамики // Меж-дунар. конф. по матем. теории управления и механике. Тез. докл. Суздаль, 7-11 июля 2017. Владимир: ООО «Аркаим», 2017. С. 36-37.
16. Беспалов М. С., Панина Н. А. Программа вычисления точного значения сумм четных отрицательных степеней синусов в равноотстоящих узлах окружности. Зарегистрированная программа для ЭВМ. Свидетельство о гос. регистрации программ для ЭВМ №2011616549 от 22 августа 2011г.
17. Долгий Д. В., Ким Д. С., Ким Т. О полиномах Коробова первого рода // Матем. сб. 2017. Т. 206. №1. С. 65-79.
18. Roman S. M. and Rota G.-C. The umbral calculus // Adv. Math. 1978. V. 27. P. 95-188.
REFERENCES
1. Ber M. G. and Malozemov V. N. The best formulas for approximate calculastion of a discrete Fourier transform. Comput. Math. and Math. Physics, 1992, vol.32, no. 11. pp. 1533-1544.
2. Izbrannye glavy discretnogo garmonicheskogo analiza i geometricheskogo modelirovaniya. Ch. 2. [Selected Chapters of Discrete Harmonic Analysis and Geometric Modeling. P. 2]. Ed. V.N. Malozemov. St. Petersburg, VVM Publ., 2014. 605p. (in Russian)
3. Malozemov V. N. and Macharskiy S. M. Osnovy discretnogo garmonicheskogo analiza [Basics of Discrete Harmonic Analysis]. St. Petersburg, Lane Publ., 2012. 304 p. (in Russian)
4. Malozemov V. N. and Pevniy A. B. Discrete periodic splines and their computational application. Comput. Math. and Math. Physics, 1998, vol.38, no.8. pp. 1181-1192.
5. Bespalov M. S. Representation for sums of even negative degrees of sines in equivdistant nodes. Russian Mathematics, 1996, vol. 40, no. 8, pp. 4-10.
6. Korobov N. M. Special'nye polinomy i ikh prilojeniya [Special polynomials and its application]. Diofantovy Pribligeniya. Matematicheskie Zapiski, 1996, vol.2, pp. 77-89. (in Russian)
7. Ustinov A. V. O formulakh summirovaniya i interpolyacii [On summation and interpolation formulas]. Chebyshevskii Sbornik, 2001, vol. 1, no. 1, pp. 52-71. (in Russian)
8. Ustinov A. V. Polynomy Korobova i tenevoi analis [Polynomials of Korobov and umbral analysis]. Chebyshevskii Sbornik, 2003, vol.4, no.4(8), pp. 137-152. (in Russian)
9. Horn R. A. and Johnson C. R. Matrix Analysis. Cambridge, Cambridge University Press, 1986. 662 p.
10. Trakhman A. M., TrakhmanV.A. Osnovy teorii signalov na konechnikh intervalakh [Fundamentals of the Theory of Signals of Finite Intervals]. Moscow, Soviet Radio Publ., 1975. 208 p. (in Russian)
11. Kombinatorniy analis. Zadzchi i upragneniya [Combinatorial Analysis. Tasks and Exercises]. Ed. K. A. Rybnikov. Moscow, Nauka Publ., 1982. 368 p. (in Russian)
12. Bespalov M. S. O proizvodyachikh funkciyakh dlia nekotorikh trigonometricheskikh summ [On generating functions for same trigonometric sums]. Proc. XXV Sci. Conf. VPI. P. 1. Vladimir, VPI Publ., 1990, p. 40. (in Russian)
13. Guzev M. A. and Ustinov A. V. Mechanical characteristics of molecular dynamics model and Korobov polynomials. Dal'nevostochnyi Matematicheskii Zhurnal, 2016, vol.16, no. 2, pp. 39-43. (in Russian)
14. ZorichV.A. Matematicheskii analis [Mathematical Analysis]. P.2. Moscow, Nauka Publ., 1984. 640 p. (in Russian)
15. Bespalov M. S. Trigonometricheskie summy dlia zadach molekuliarnoy dinamiki [Trigonometric sums for problems of molecular dynamics]., Intern. Conf. Math. Control Theory and Mechanics. Vladimir, Arcaim Publ., 2017, pp. 36-37. (in Russian)
36
М. С. Becna^OB
16. Bespalov M. S. and Panina N. A. Program calculation exact value for sums of even negative degrees of sines in equivdistant nodes. Program for computation no. 2011616549. 22.08.2011.
17. Dolgy D. V., Kim D. S., and Kim T. Korobov polynomials of the first kind. Sb. Math. 2017, vol.208, no. 1. pp. 60-74.
18. Roman S. M. and Rota G.-C. The umbral calculus. Adv. Math., 1978, vol.27, pp. 95-188.
