Три теоремы о матрицах Вандермонда
Автор: Артисевич Анжела Евгеньевна, Шабат Алексей Борисович
Журнал: Владикавказский математический журнал @vmj-ru
Статья в выпуске: 1 т.22, 2020 года.
Бесплатный доступ
Рассматриваются алгебраические вопросы, связанные с дискретным преобразованием Фурье, определенным при помощи симметричной матрицы Вандермонда Λ. Основное внимание в первых двух теоремах уделяется выработке формулировок, независящих от размера N×N матрицы Λ и явных формул для элементов матрицы Λ через корни уравнения λN=1. В третьей теореме рассматриваются рациональные функции f(λ), λ∈C, удовлетворяющие условию "вещественности" f(λ)=f(1λ) на всей комплексной плоскости и связанные с известной задачей о коммутировании симметричных матриц Вандермонда Λ с (симметричными) трехдиагональными матрицами T. Показано, что уже несколько первых уравнений коммутирования и указанное выше условие вещественности определяют вид рассматриваемых рациональных функций f(λ), а найденные уравнения для элементов трехдиагональных матриц T не зависят от порядка N коммутирующих матриц. Полученные уравнения и приведенные примеры позволяют высказать гипотезу о том, что рассматриваемые рациональные функции являются обобщением многочленов Чебышева. В определенном смысле аналогичная гипотеза была высказана в недавно опубликованной в журнале "Теоретическая и математическая физика" работе В. М. Бухштабера с соавторами, где обсуждаются приложения этих обобщений в современной математической физике.
Матрица вандермонда, дискретное преобразование фурье, условия коммутирования, многочлены лорана
Короткий адрес: https://sciup.org/143170629
IDR: 143170629 | DOI: 10.23671/VNC.2020.1.57532
Текст научной статьи Три теоремы о матрицах Вандермонда
/1 X1x i x2
...
...
...
\ 1 XN
N -1 1
N -1 2
N -1
N
\
/
x j ∈ C ,
симметричной, если A T = A. В этом случае x i = 1, а все остальные x j являются степенями одного и того же комплексного числа λ. Этот частный случай матрицы Вандермонда (1) имеет следующий вид:
/ 1 1 1 ... 1\
-
1 А А 2 ...A
λ ∈ C .
... ... ... ......
\ 1 A N -1 A 2( N -1) ... A (n - 1)(n -1) /
Известное свойство унитарности матриц дискретного преобразования Фурье (см. [1]), приводит нас к следующему уравнению для матриц Вандермонда (1):
AA* = A* A = N • E,(3)
где * обозначает эрмитово сопряжение A * = A T , а E = diag(1,..., 1) — единичную матрицу.
Теорема 1. Матрица Вандермонда (1) удовлетворяет уравнению (3) в том и только в том случае, если x j для j = 1, 2,... , N являются N различными корнями уравнения x N = e iNY , y G R . При этом матрица (1) записывается в виде произведения симметричной матрицы Вандермонда (2) с A N = 1 и диагональной A = Л • diag(1, e iY , e i 2 Y ,...) .
Последнюю формулу в теореме поясним следующим примером.
Пример 1.1. Выбор x i = e iY находится в нашем распоряжении и, выбрав при N = 3
1 1
xi = ^(V3 + i) ^ xi = ^(1+ iV3), X3 = i; A = e3 ,(4)
мы получаем «унитарную» матрицу (1), удовлетворяющую уравнению (3):
/1 Xi x2A /1 1 1A /1 00\
A = I 1 x 2 x 2 = I 1 A A 2 I I 0 x i 0 I , A 3 = 1.
1 x3 x23 1 A2 A 0 0x
Для доказательства теоремы 1 рассмотрим диагональные элементы с номерами
и 33 произведения A * A. Уравнение (3) дает для разностей y j = | x j | 2
- 1
y i + y 2 + • • • + y N = y 2 + y 2 + • • • + y N = 0 ^ y j = 0 (Vj G [N]).
Поэтому справедливо утверждение
Лемма 1. В условиях теоремы 1 выполняется | x j | = 1 для любого j.
Таким образом независимо от размера матрицы Вандермонда доказательство теоремы 1 сводится в силу формул Виета к проверке импликации xi + ... + xn — 0, x2 + ... + xN = 0, ............
тN -i 4- 4- TN-i — 0
x i —+ . . . т x n — v
x NN = a ( V j G [N]).
Остается заметить, что | a | = 1 в силу леммы 1 и что уравнение z N = a сводится к уравнению A N = 1 заменой z = e iY при a = e iNY . Доказательство обратного утверждения, 2 iπ
-
т. е. A = e n ^ (3), можно извлечь из цитированной выше монографии [1].
Пример 1.2. Пусть N = 3. Тогда восемь уравнений матричного равенства A * A = 3E можно представить в виде трех систем
I
x 1 + x 2 + x 3 = 0; x 2 1 + x 2 1 + x 2 3 = 0,
I
y 1 + y 2 + y 3 = 3; y 1 2 + y 2 2 + y 3 2 = 3,
I
У 1 Х 1 + У 2 Х 2 + У 3 Х 3 = 0;
У 1 Х 1 + У 2 Х 2 + У 3 Х 3 = 0,
def где мы сменили обозначения леммы 1 и yj = Xj Xj. В силу леммы yj = 1 О Xj Xj = 1 для любого j = 1, 2, 3. При этом последняя система уравнений совпадает с первой, а вторая превращается в тождества вида 3 = 3. Аналогично в случае N = 4
1 X1 X21
1 x2 x2
1 X3 X23
у 1 x4 x4
и 16 уравнений матричного равенства A * A = 4E, при условии что y j = X j X j = 1,
Замечание 1. Множество корней из единицы при N = 3 состоит из трех элементов, образующих циклическую группу. При этом любой из элементов Ах = e ", А2 = e 4П1 можно использовать в качестве образующей циклической группы. Если решать рассматриваемую полиномиальную систему из 8 уравнений при помощи вычислительной техники, то решения записываются в виде, аналогичном (4), но со свободным параметром. Компьютер выдает два таких решения (с учетом перестановок). Следует заметить, что уже для случая N = 5 вычислительная техника испытывает затруднения и не доводит решение до конца.
Следующая теорема показывает, что для симметричных матриц Вандермонда Л вида (2) условие AA * = N • E из теоремы 1 можно заменить условием Л 2 = N • Q, где Q — матрица перестановок, состоящая из нулей и единиц:
1 0 ... 00
.
0 0 .. 01
.... . . .. ..
. ..
\0 1 . ... 0/
Теорема 2. Симметричная N х N матрица Вандермонда (2) является матрицей дискретного преобразования Фурье и удовлетворяет условию A N = 1 тогда и только тогда, когда Л 2 = N • Q, где Q — симметричная матрица перестановок (5) .
⊳ Доказательство этой теоремы приведено в монографии [1], и мы ограничимся замечанием, что для вывода основного уравнения A N = 1 достаточно приравнять к нулю элемент с номером 12 в матрице Л 2 , так как при А = 0
-
2 . Канонические многочлены Лорана
Следуя работам [2, 3], рассмотрим теперь связанное с матрицами Вандермонда (1) коммутационное уравнение TA = AT , где T — трехдиагональная матрица общего вида:
/ ai bi 0 ... 00\
T def bi a2 b2 ... 00
... ... ... ......
\ 0 0 ...... 6 n -1 gn)
элементы которой можно выразить через элементы матрицы (1). В дополнение к работе [2] мы покажем, что в условиях теоремы 2, т. е. для симметричных матриц Вандермонда вида (2), диагональные элементы a j матрицы (6) записываются в виде «канонических» многочленов Лорана:
P n (X) = П ( X+ Y j + x ) .
Отметим, что коэффициенты этих «многочленов» степени 2n определяются только их номером n (см. ниже) и не зависят от размера N рассматриваемой симметричной матрицы (2).
-
2.1. Необходимые условия коммутирования. Схема вывода формул, выражающих элементы треугольной матрицы через элементы матрицы Вандермонда (1), мало зависит от предположения симметричности A T = A, и мы ограничимся ниже именно этим случаем, предполагая для простоты искомую матрицу (6) также симметричной.
Коммута тор симметрических матриц кососимметричен и его первая строка с элементами 12 , 13 , . . . позволяют выразить разности диагональных элементов a 2 - a 1 , a 3 - a 1 , . . .
через недиагональные. В результате получаем
12 : a i + b i X = b i + a 2 + b 2 , * 13 : a i + b i X 2 = b 2 + a з + Ь з
14 : a i + b i X3 = Ь з + a 4 + b 4
a 2 - a i = Xb i - b i - 6 2 , а з — a i = X 2 b i — b 2 — Ь з , a4 — a i = X3b i — Ьз — b 4 ,...
Вторая строка коммутатора с элементами 23 , 24 , 25 . . . дает
: b i + a 2 X 2 + b 2 X 4 = b 2 X + a з X 2 + b з X 3 ,
: b i + a 2 X 3 + b 2 X 6 = b з X 2 + a 4 X 3 + b 4 X 4 , : b i + a 2 X 4 + b 2 X 8 = b 4 X 3 + a 5 X 4 + b 5 X 5
и, подставив сюда найденные из предыдущих уравнений разности диагональных элементов, находим уравнения для выражения bj, j ^ 2, через bi и b2. В частности, ьз(х) = fx + X + 1) b2(X) — fX + X + Т2) bi(X)’
4 \ (10)
b4(X) = ( X2 + X + 2 + V + V2 ) b2(X) — ( X2 + X + 1 + Т + "\2' + Тз ) bi(X).
λ λ 2 λ λ 2 λ 3
Итак, можно считать доказанной следующую теорему 3.
Теорема 3. Первые две строки уравнений коммутативности (8) и (9) позволяют найти вид коэффициентов трехдиагональной матрицы T в форме многочленов Лорана от независимой переменной А и выразить их через b i (A) , 6 2 (A) .
Сравнив полученные выражения с уравнением, полученным из элемента 34 коммутатора b2A3 + азА6 + ЬзА9 = A% + A6a4 + A9b4,
приходим к выводу (ср. [2]), что условие b i = 0 является необходимым условием для выполнения коммутационных соотношений TA = AT в симметричном случае (2).
-
2.2. Формулы Виета. Рассматривая условия разрешимости коммутационных соотношений в кольце многочленов Лорана c j λ j от формальной переменной λ, будем использовать следующие обозначения:
^ cj Aj = [^ cj Aj ]- + [^ cj Aj ]+; [^ cj Aj ]+d=f jE cj Aj, где квадратные скобки [...]+ и [•••]- обозначают сумму соответственно членов с отрицательными и положительными степенями λ.
Лемма 2. Пусть b i = 0, b 2 = 1. Тогда все многочлены Лорана b m (A) и a m (A), найденные из уравнений (8) и (9) , инвариантны относительно замены λ ↔ λ -1 , и для многочлена a N +2 (A) при любом N > 0 имеем
— [a N +2 ] + = A N + 2A N 1 + 3A N 2 + • • • + NA + N, [ A a N +2]— = [a N +i ] - , N > 0. **
-
<1 Уравнения первой строчки при b i = 0, b 2 = 1 и a i = 1 дают
a 2 = 0, a 3 + b 3 = 0, a m = 1 — b m — b m -i , m > 3. (13)
Учитывая формулы (8), (9), получаем теперь (ср. [2]), что при m > 0
b m +i = Ab m + 1 + A + • • • + A m 1 , b i = 0, b 2 = 1
Найденные по этим формулам многочлены Лорана b m (A) инвариантны относительно замены λ ↔ λ -1 и для доказательства уравнений (12) остается применить индукцию к b m + b m -i — 1 >
Следствием полученных выше формул является
Лемма 3. При m E 3 комплексные корни A m = 1 из единицы являются комплексными нулями многочленов Лорана a m (A) , m E 3 , из леммы 3 .
-
< Заменив при помощи уравнения A m = 1, m = N + 2, отрицательные степени A в формуле (12) на положительные
P N (A) = A N + 2A N -i + 3A N -2 + ••• + NA + N + N + ••• + A N , \N +2 _ i _ \2 1 _ x3 i _ xN +i
A 1 ^ an a ’ an-i A ’•••’ a A ’ после приведения подобных членов получаем
P N (A) = N ( A N +i + A N + ••• + 1 ) = 0, mod A N +2 = 1. ▻ (14)
Очевидно, общее число 2n комплексных корней многочленов Лорана P n (A) из формулы (7) превосходит при n > 2 число n + 2 соответствующих комплексных корней из единицы. Дополнительные корни уравнения P n (A) = 0 можно найти, используя обобщенные формулы Виета , приспособленные к произведениям вида (7):
(А + Y 1 + 1)(
А + Y 2 + 1)
= А 2 + Д 2 + (Y 1 + Y 2 ) (а + 1)
+ Y 1 Y 2 + 2,
P n ( A ) = П (А + Y j + а) = A n + А П + " (A n 1 + А " — 1) + °" 2 (A n 2 + А п - 2) + "‘+ "
Здесь коэффициенты " j выражаются через элементарные симметрические многочлены ° i , i ^ j,
"1 = °1 = У? Yi; "2 = °2 + n = n + У2YiYj, n> 2, i "j (n + 1) = "j (n) + Yn+1°j-1(n) + "j-2(n), "0 = 1. Пример 3. При n = 3 комплексные корни уравнения А5= 1 дают 4 из 6 корней рассматриваемого уравнения: Р3(А) = А3 + Дз + 2^А2 + д2 ^ + 3(А + А) + 3 = П (А + Yj+ А) = 0 Формулы Виета приводят в этом случае к следующей системе уравнений для γj : Y1 + Y2 + Y3 = 2, Y1Y2 + Y1Y3 + Y3Y2 = 0, Y1Y2Y3 = —1- Можно проверить независимо, что рассматриваемый многочлен Рз(А) является приводимым и факторизуется следующим образом: Р3(А)= (А + 1 + 1) (А2 + А2 + А + 1 + 1) .
Список литературы Три теоремы о матрицах Вандермонда
- Бурланков Д. Е., Кузнецов М. И., Чирков А. Ю., Яковлев В. А. Компьютерная алгебра. Нижний Новгород: Нижегородский гос. ун-т им. Н. И. Лобачевского, 2002. 105 с.
- Grunbaum F. A. The eigenvectors of the discrete Fourier transform: a version of the Hermite functions // J. Math. Anal. Appl. 1982. Vol. 88, № 2. P. 355-363. DOI: 10.1016/0022-247X(82)90199-8
- Шабат А. Б. Симметрические многочлены и законы сохранения // Владикавк. мат. журн. 2012. Т. 14, вып. 4. С. 83-94. DOI: 10.23671/VNC.2012.14.11014
- Бухштабер В. М., Тертычный С. И. Семейство явных решений уравнений резистивной модели перехода Джозефсона // Теор. и мат. физика. 2013. Т. 176, № 2. С. 163-188. DOI: 10.4213/tmf8512