CC BY
9ISSN 1814-1196 Научный вестник НГТУ том 56, № 3, 2014, с. 23-36
http://journals. nstu. ru/vestnik Scientific Bulletin of NSTU Vol. 56, No. 3, 2014, pp. 23-36
АВТОМАТИЧЕСКОЕ УПРАВЛЕНИЕ AUTOMATIC CONTROL
И РЕГУЛИРОВАНИЕ AND REGULATION
УДК 511.82
Разложение на простые множители невырожденных
*
полиномиальных матриц
А.Н. КОРЮКИН 1
1 630090, РФ, г. Новосибирск, проспект академика Коптюга, 4, Институт математики им. Соболева,
с. н. с., к. ф.-м. н., е-таЛ: koryukin@sibmail.ru
Предмет исследования. В алгебре полиномиальных матриц (над алгеброй К полиномов от одной буквы) размера п на п изучается вопрос единственности разложения на простые множители. Актуальность. Теория делимости является классической частью математики, вошедшей в процесс обучения в высших учебных заведениях. Автору неизвестно, изучался ли вопрос о единственности разложения на простые множители. Между тем в линейной теории автоматического управления (многоканальные системы) давно и успешно используются понятия наибольшего левого и правого делителей полиномиальных матриц, что может служить основой теории делимости для полиномиальных матриц. Важнейший вопрос теории делимости - существование и единственность разложения на простые множители. Вопрос единственности разложения на простые множители затруднён тем, что в кольце матриц нет коммутативности. Результаты. Показано, что в кольце полиномиальных матриц левые и правые идеалы главные (матрицы над алгеброй полиномов от одной буквы). Отсюда следует, что для любого числа невырожденных матриц существуют наибольший левый делитель, наибольший правый делитель, наименьшее левое кратное, наименьшее правое кратное. Показано, что для любой невырожденной матрицы существует разложение на простые множители. Охарактеризованы неразложимые матрицы: это в точности матрицы с простым (неразложимым) определителем. Проблема единственности разложения на простые множители решается так: для произвольной невырожденной матрицы А описаны все простые правые делители (простые матрицы В, для которых существует матрица С такая, что А = ВС). Это описание сделано в терминах решётки К-подмодулей п-ок (модуля последовательности полиномов из К длины п. Разложение же А = ВС позволяет найти все разложения невырожденной матрицы А на простые множители. Эти поиски записаны в виде двух алгоритмов (для математиков и для инженеров). Описаны некоторые инварианты разложения матрицы на простые множители. Описаны некоторые частные случаи разложения. Приведён один пример разложения. Методика и инструментарий. Теория колец. Теория делимости.
Ключевые слова: полиномиальные матрицы, простые множители, неразложимые элементы, разложение на простые множители, единственность разложения на простые множители, теорема об единственности разложения на простые множители, наибольший общий левый делитель, наибольший общий правый делитель
ВВЕДЕНИЕ
Каждое натуральное число может быть представлено в виде произведения простых одним и только одним способом. Эту теорему называют основной теоремой арифметики [1].
В трудах Евклида [2] она ещё не встречается, но в предложении 31 книги VII было установлено, что любое целое число может быть представлено в виде произведения простых множителей, где некоторые множители могут быть одинаковыми; в предложении 14 книги IX доказано, что если все сомножители различны, то это представление однозначно [3].
Существование наибольшего общего делителя доказывается в самом начале VII книги Евклида [2] с помощью процесса, известного под названием «алгоритма Евклида». Этот алгоритм является основой всех последующих рассмотрений (свойства простых чисел,
* Статья получена 01 июля 2013 г.
существование и нахождение наименьшего общего кратного и т. д.). Однако существование и единственность разложения на простые множители не было доказано в общем виде. Но, во всяком случае, Евклид явно доказывает, что каждое целое делится на простое число (предложение 31 книги VII) [4].
По-видимому, Евклид потому не формулирует общую теорему, что у него не хватает адекватной терминологии и символики для произвольных степеней целого числа. К тому же все свидетельства согласуются с тем, что начиная с этой эпохи разложение заданного числа на простые множители было хорошо известно и им свободно пользовались [4].
Основная теорема арифметики не формулируется точно даже во «Введении в теорию чисел», написанном в 1798 году Лежандром. Вероятно, из-за отсутствия этой теоремы у Евклида она принимается без доказательства во многих школьных учебниках [1].
Первую точную формулировку этой теоремы и её доказательство дал Гаусс в 1801 году в своих знаменитых «Арифметических исследованиях» [5] (во втором разделе этих исследований в п. 16 доказывается однозначность разложения натурального числа на простые множители). По этому поводу Гаусс говорит: «То, что каждое составное число может быть разложено на простые множители, известно из основ; однако то, что этого нельзя сделать несколькими различными способами, совершенно необоснованно по большей части считается само собой разумеещимся».
Вопрос о делении многочленов как об операции обратной умножению, естественно, должен был быть поставлен в самом начале развития алгебры. Ясно, однако, что едва ли возможно решить эту проблему до тех пор, пока не создана последовательная система обозначений для различных степеней переменного. Действительно, до середины XVI века не появляется ни одного примера «евклидова» деления многочленов, и, видимо, Стевин [6] (который, в сущности, использовал обозначение показателя) был первым, кто пришёл к мысли применить «алгоритм Евклида» для отыскания наибольшего общего делителя двух многочленов [7, стр. 113]. В остальном понятие делимости до середины XVIII века считалось присущим только целым числам [4].
В XIX веке в теории делимости был достигнут значительный прогресс. Она была развита для кольца целых чисел, колец полиномов, областей целостности (ассоциативных коммутативных колец без делителей нуля). Были введены евклидовы кольца (кольца с алгоритмом Евклида), кольца главных идеалов, факториальные кольца (области целостности с единственным разложением на простые множители). Показано, что евклидовы кольца являются кольцами главных идеалов, а кольца главных идеалов являются факториальными; что кольцо полиномов над факториальным кольцом факториально [7, 8]. Значительный вклад в развитие теории делимости внёс Гаусс.
Теория делимости стала небольшой и классической частью основ математики, и её стали преподавать студентам высших учебных заведений [9-12].
Но в теории делимости математики ограничились изучением коммутативных колец, причём не имеющих делителей нуля. Между тем в линейной теории автоматического управления (многоканальные системы) по необходимости давно и успешно используются наибольший общий левый, правый делитель полиномиальных матриц [13-18] (полиномы от одной буквы). Понятия наибольшего левого и правого делителей полиномиальных матриц являются стандартными и классическими элементами инструментария теории многоканальных систем, и это преподают студентам высших учебных заведений [19], [20].
Впрочем, оказалось, что понятия наибольшего общего левого и правого делителя полиномиальных матриц были математикам известны [21], [22].
Очевидно, что алгебра матриц над алгеброй полиномов от одной буквы нётерова (как слева, так и справа). Поэтому в ней можно говорить о простых левых и правых делителях, и существует разложение на простые множители (по крайней мере, для невырожденных матриц). Но как быть с единственностью разложения на простые множители? Этот вопрос и исследуется в данной работе.
1. ЛЕВЫЙ/ПРАВЫЙ НАИБОЛЬШИЙ ОБЩИЙ ДЕЛИТЕЛЬ И ЛЕВОЕ/ПРАВОЕ НАИМЕНЬШЕЕ ОБЩЕЕ КРАТНОЕ
Далее k - произвольное поле; 5 - символ; К = ^я] - алгебра (векторное пространство с умножением) полиномов над полем к, п > 1 - произвольное целое положительное число; М = Мп (К) - алгебра матриц размера п х п над алгеброй К; и (М) - множество невырожденных матриц, у которых определитель лежит в k.
и (М) - это в точности множество матриц из М, обратимых в М.
Элементы алгебры М называем просто матрицами. Множество {1,..., п} обозначим через
п . Для i, / е п обозначим через ег / матрицу, у которой в г-ой строке и /-ом столбце единица,
а остальные нули. Такие матрицы будем называть, как это обычно принято, матричными единичками.
Введём обозначения: К п - К-модуль столбцов высоты п полиномов из К; Кп - К-мо-дуль строк длины п полиномов из К.
Напомним, что левый идеал, порождённый одним элементом, называется главным.
Лемма 1. Для п > 1, произвольный подмодуль любого п-порождённого К-модуля порождён не более чем п элементами.
Действительно, алгебра К коммутативна и нётерова (т. е. любая возрастающая цепочка идеалов алгебры К обрывается). Значит, любой конечнопорождённый К-модуль N нётеров. Поэтому все его подмодули конечнопорождены.
Достаточно рассматривать случай, когда К-модуль N свободно порождён своими п элементами. В этом случае элементы К-модуля N можно рассматривать как п-ки полиномов из К.
Рассмотрим произвольный подмодуль Р К-модуля N. Он конечнопорождён. Рассмотрим любое конечное множество п-ок, порождающее К-модуль Р, и следующую «элементарную операцию» над такими множествами п-ок: к одной п-ке можем прибавить другую, помноженную на произвольный полином из К.
Понятно, что в результате этих элементарных операций К-модуль, порождённый множеством п-ок, не изменится.
Будем обращать внимание только на первый элемент п-ок. Вспомним об алгоритме деления с остатком в алгебре К. Теперь элементарными операциями множество п-ок может быть преобразовано так, что только один из элементов этого множества может иметь ненулевую первую координату.
После этого исключим из множества п-ок имеющую ненулевую первую координату. В оставшихся п-ках первая координата нулевая. Исключим первую координату и продолжим процесс.
Теперь индукцией по числу п нетрудно показать, что произвольный подмодуль любого п-порождённого К-модуля порождён не более чем п элементами.
Лемма 2. Любой левый идеал алгебры матриц М главный.
Обозначим через ц произвольный левый идеал алгебры М. Так как К • Е с М (Е -единичная матрица), то ц является К-модулем.
Пространство Кп является К-модулем. Этот модуль порождается своими п элементами. Пространство Кпц является его подмодулем. Значит, согласно лемме podmodul-n-porojd, К-модуль Кпц порождён не более чем п элементами.
Возьмём эти п элементов и составим из них матрицу. Если элементов меньше, чем п, то допишем пустые строки. Эта матрица порождает исходный левый идеал ц. Таким образом, левый идеал ц главный.
Следствие 3. Правый идеал алгебры М главный.
(Из леммы 2 - транспонированием).
Следствие 4. Для любого множества матриц из М существует их правый (левый) наибольший общий делитель и правое (левое) наименьшее общее кратное.
Для А, В е М , правые идеалы АМ п ВМ , АМ + ВМ главные, т. е. СМ = АМ п ВМ DM = АМ + ВМ для некоторых С, Б из М. Тогда С будет наименьшим общим левым кратным матриц А, В; Б будет наибольшим общим левым делителем матриц А, В.
Для правых делителей и кратных Утверждение получается аналогично, или транспонированием (из Утверждения для левых делителей и кратных).
2. СУЩЕСТВОВАНИЕ РАЗЛОЖЕНИЯ НА ПРОСТЫЕ МНОЖИТЕЛИ НЕВЫРОЖДЕННЫХ ПОЛИНОМИАЛЬНЫХ МАТРИЦ
Лемма 5. Любая возрастающая цепочка левых (правых) идеалов алгебры М обрывается.
Алгебра К нётерова; К-модуль М конечно порождён. Поэтому М является нётеровым К-модулем.
Рассмотрим произвольную возрастающую цепочку Р1 с... срп с... левых (правых) идеалов алгебры М. Элементы этой цепочки является К-модулями. Это возрастающая цепочка подмодулей К-модуля М.
Из изложенного выше следует, что любая возрастающая цепочка левых (правых) идеалов алгебры М обрывается.
Напомним, что квадратная матрица над ассоциативным коммутативным кольцом называется невырожденной, если её определитель не ноль.
Обратимыми элементами алгебры М являются матрицы, определитель которых лежит в поле k и отличен от нуля.
Определение 6. В алгебре М простыми (неразложимыми) назовём невырожденные необратимые элементы, у которых при любом разложении в произведение двух сомножителей один из сомножителей будет обратим.
Лемма 7. В алгебре М для любой невырожденной необратимой матрицы любой процесс разложения в произведение необратимых сомножителей можно довести до конца.
Рассмотрим элемент А е М из условия Леммы. Если он простой, то он раскладывается в произведение простых сомножителей.
Далее элемент А не простой. Тогда он разлагается в произведение двух сомножителей А = ВС.
Для произвольной матрицы Q через det(Q) будем, как это общепринято, обозначать определитель матрицы Q.
Матрицы В, С не вырождены (ведь А = ВС и матрица А невырождена). Матрицы В, С необратимы (ведь А = ВС и невырожденная матрица А неразложимая). Поэтому полиномы det(B), det(C) не являются константами. Значит, эти полиномы являются собственными делителями полинома det(А).
Доказательство завершается индукцией по степени полинома det(A).
Следствие 8. В алгебре Млюбая невырожденная необратимая матрица раскладывается в произведение простых сомножителей.
3. К ЕДИНСТВЕННОСТИ РАЗЛОЖЕНИЯ НА ПРОСТЫЕ МНОЖИТЕЛИ -ДЛЯ НЕВЫРОЖДЕННЫХ МАТРИЦ. ДИАГОНАЛЬНОЕ РАЗЛОЖЕНИЕ НЕВЫРОЖДЕННЫХ МАТРИЦ НА ПРОСТЫЕ МНОЖИТЕЛИ
В алгебре К = ^ х]:
1) обратимые элементы - это отличные от нуля константы;
2) выполнена теорема о единственности разложения на простые множители.
В алгебре М любую матрицу А можно разложить так: А = UDV, где и, V - обратимые матрицы из М; Б - диагональная матрица из М [19, § 2, пункты 16, 17]. Её называют формой Смита.
Для невырожденной матрицы А по диагонали матрицы D расположены отличные от нуля полиномы. Разложим их на множители. После этого разложим матрицу D (в произведение наибольшего числа необратимых сомножителей): Б = Б •... • Бт , где т > 1.
Понятно, что сомножители , где г е т , имеют простой определитель.
Теперь понятно, что А разлагается в произведение т сомножителей, имеющих простые определители: А = иБ •... • А •... • БтУ .
Следствие 9. В алгебре М любая невырожденная простая матрица имеет простой определитель.
Так как матрица А простая, то из диагонального разложения следует т = 1, т. е. А = UDV, где матрицы и, V обратимы, а определитель матрицы Б простой. Значит, определитель матрицы А простой.
При этом очевидно, что матрицы из М с простым (неразложимым в К) определителем простые (неразложимые в М). Поэтому:
Лемма 10. Пусть А - невырожденная матрица из М. Тогда матрица А простая (неразложимая) » её определитель простой.
4. ОТ ЕДИНСТВЕННОСТИ РАЗЛОЖЕНИЯ К ЖОРДАНОВОЙ ФОРМЕ
Рассмотрим произвольную невырожденную матрицу А из М и произвольное разложение А = р •...• Рт матрицы А на простые множители.
1. Разложения устойчивы при умножении на матрицы из и(М) и справа, и слева: AQ = р •...• есть разложение матрицы AQ в последовательность т простых множителей р, ..., Р^ (Q еи(М)); ТА = Тр •...• Рт есть разложение матрицы Q в последовательность т простых множителей ТР[, ..., Рт (Т еи (М)).
2. Для любой матрицы разложений на простые множители много. Например, равенство Р • Р2 = Р^ • Q"1P2 даёт бесконечно много разложений матрицы р • Р> в произведение. Поэтому хотелось бы найти в этих разложениях что-то более устойчивое и инвариантное.
3. Рассмотрим цепочку правых идеалов рМ з РР2М з... з р...РтМ:
1) эта цепочка конечна (из т элементов) и заканчивается правым идеалом АМ = р •...• РтМ ;
2) так как матрицы р, ..., Рт невырожденные и не обратимы в М, то во всех звеньях цепочки нет равенства: рМ з РР2М з... з АМ ;
3) так как матрицы р, ..., Рт простые, то эта цепочка правых идеалов плотная.
Так как матрица р простая, то фактор-модуль М / рМ простой. Поэтому цепочка правых идеалов М з> рМ з> р^р2М з> ... з> АМ будет плотной.
4. Теперь понятно, что для единственности разложения матрицы А полезно изучить решётку правых идеалов р алгебры М таких, что р з АМ. Как устроена эта решётка? Она конечна? Дистрибутивна?
5. Обозначим через К" множество всех столбцов из К высоты п .
Понятно, что К п является левым М-модулем (при умножении матрицы на столбец
получаем столбец). В частности, К п также является К-модулем (ведь поле К вкладывается в М с помощью скалярных матриц: К • Е с М).
6. Для любой матрицы В рассмотрим её образ ВКп (при действии матриц слева на столбцы из Кп). Понятно, что ВКп есть подмодуль К-модуля К"; если V - некоторый подмодуль К-модуля М , то VKn является подмодулем К-модуля К".
7. Каждому правому идеалу р сопоставим рКп:.
1) рКп является подмодулем К-модуля Кп;
2) рассмотрим произвольный подмодуль V К-модуля К". Он порождён не более чем п элементами (лемма 1). Составим из этих элементов-столбцов матрицу. Если столбцов меньше, чем п, то в матрицу добавим пустые столбы и сделаем матрицу квадратной. Обозначим её В;
3) ВМКп = ВКп = V , т. е. образ рКп правого идеала р = ВМ алгебры М, порождённого матрицей В, совпадает с произвольно выбранным подмодулем К-модуля К".
8.1. Решётка правых идеалов алгебры М изоморфна решётке подмодулей К-модуля К".
8.2. Отображение f: р ^ рКп, сопоставляющее произвольному правому идеалу р
алгебры М подмодуль рКп К-модуля К", является изоморфизмом решёток.
8.3. Рассмотрим отображение g : V ^ ц , сопоставляющее любому подмодулю V К-модуля Кп множество ц матриц В из М таких, что ВКп с V. Так как ВМКп с ВКп с V, то ц является правым идеалом. Поэтому отображение g можно рассматривать как отображение
решётки подмодулей К-модуля К" в решётку правых идеалов алгебры М.
Отображение g будет обратно отображению /. Поэтому g будет изоморфизмом решёток.
8.4. Пусть р - произвольный правый идеал алгебры М ; Ш = рК" .
Обозначим через R(M, р) решётку правых идеалов алгебры М, содержащих р ; через R(K" ,Ш) - решётку подмодулей К-модуля К", содержащих Ш.
При ограничении отображения / на правые идеалы, содержащие р, и при ограничении отображения g на подмодули К-модуля К", содержащие Ш, получим: решётки R(M,р), R(K",Ш) изоморфны; отображения f \щМ ): R(M,р) ^ R(K",Ш), g \д(Х" Ш):R(K", Ш) ^R(M,р) являются изоморфизмами решёток.
9. Итак, изучать надо подмодули К-модуля К". Точнее: как устроена решётка
R(K" ,Ш)?
Эта решётка изоморфна решётке подмодулей К-модуля К" / Ш; Ш = АК" для некоторой невырожденной матрицы А . Обозначим через А матрицу, дополнительную к А . Имеем: АА = А-Е , где Д = det(A) - определитель матрицы А ; Е - единичная матрица.
Поэтому АЕ с АА , ДЕК" с ААК" , А-К" с АК" = Ш . При этом А ф 0.
Значит, К-модуль К" / Ш сингулярный, т. е. любой его элемент имеет ненулевой
аннулятор. Так как при этом К-модуль К" / Ш конечнопорождён, то К" / Ш является
конечномерным ^пространством. На этом конечномерном ^пространстве К" / Ш символ 5 действует как линейное отображение. И нас интересует решётка подмодулей К-модуля
К" / Ш . То есть нас интересует решётка подмодулей одного линейного преобразования 5 на конечномерном пространстве. Понятно, что нужно вспомнить жорданову форму линейного преобразования.
10. Проще всего линейное преобразование выглядит для алгебраически замкнутого поля. Поэтому интересен случай, когда поле k алгебраически замкнуто. В частности, когда k -поле комплексных чисел X.
11. Вопрос. Пусть k - произвольное алгебраически замкнутое поле; V - произвольное конечномерное k - пространство; £ - произвольное линейное преобразование пространства V. Как устроена решётка инвариантных подпространств пространства V?
5. АЛГОРИТМ ОТЩЕПЛЕНИЯ: ОБОСНОВАНИЕ
Далее А - невырожденная матрица из М. Приведём теперь алгоритм отщепления справа от А. Он позволяет находить все разложения вида ВС = А , где С - простая матрица из М.
1.1. В алгебре М рассмотрим правый идеал АМ . Рассмотрим также подмодуль
АМКп = АКп К-модуля Кп. К-модуль АКп порождается столбцами матрицы А. Рассмотрим все п столбцов А^, ... , Ап матрицы А .
1.2. Рассмотрим фактор-модуль Кп / (АКп). Введём обозначения: Д = det(A) -определитель матрицы А ; Е - единичная матрица (размера п х п).
1.3. Существует матрица А такая, что АА = АА = Д • Е . Так как матрица А невырожденная, то матрица А единственная. Она называется дополнительной к А .
2.1. К-модуль Кп / (АКп) аннулируется полиномом Д = det(A). Действительно, АА = Д • Е , где А - матрица, дополнительная к А . Значит, ААКп = Д • Кп . Поэтому АКп з ААК" з Д • Кп.
2.2. К-модуль Кп / (АКп) является конечномерном ^пространством (ведь АКп з Д • Кп).
2.3. На конечномерном пространстве Кп / (АКп) линейно действует элемент 5.
2.4. Д ^ 0 (ведь матрица А невырожденная).
3.1. Если полином Д степени 0 (константа), то матрица А обратима, и раскладывать на множители нечего. Алгоритм заканчивается.
3.2. Далее: полином Д не ноль, но не обратим в К .
3.3. Разложим полином Д на простые множители.
3.4. Выберем простой сомножитель d полинома Д.
4.1. Рассмотрим фактор-алгебру К = К /^ = 0) алгебры К по идеалу, порождённому элементом d . Так как элемент d простой, то алгебра К является полем.
4.2. Рассмотрим фактор-алгебру матриц М = М / ^ • Е = 0) алгебры М по идеалу, порождённому элементом d • Е .
4.3. Обозначим через А образ матрицы А при каноническом гомоморфизме М ^ М .
5.1. Так как в алгебре К элемент d является делителем элемента Д, то Д = 0, т. е. определитель матрицы А нулевой, т. е. матрица А вырожденная. Значит, её аннулятор
= (А, К") в К -модуле К ненулевой.
5.2. Возьмём произвольный ненулевой элемент X из этого аннулятора.
5.3. Рассмотрим произвольный прообраз X (при каноническом гомоморфизме
Кп ^ Кп) элемента X из Кп.
5.4. В k(5) -модуле k(5)п рассмотрим элемент Y = АХ • d_1. Здесь k(5) - поле дробей
полиномов из ^5]; k(5)п - множество столбцов из k(5) высоты п . По построению, w е Кп .
6.1. Рассмотрим подмодуль V К-модуля К", порождённый столбцами А^, ... , Ап, Y. Этот модуль порождён не более чем п элементами (лемма 1).
6.2. С помощью элементарных преобразований столбцов приведём множество А1, ... , А", У столбцов к ступенчатому виду (нули будут над ступеньками).
6.3. Исключим нулевые столбцы.
6.4. Получим множество N из не более чем п столбцов, порождающее К-модуль V .
7.1. Так как АК" с V и ^пространство К" / АК" конечномерно, то ^пространство
К" / V также конечномерно. Значит, множество N не может содержать менее чем п элементов.
7.2. Поэтому множество N содержит ровно п элементов.
7.3. Составим из п элементов множества N матрицу. Обозначим её через В .
7.4. В силу выбора матрицы В , АМ с ВМ , АМ ф ВМ , и фактор-модуль над М ВМ / АМ простой.
8.1. Из АМ с ВМ следует А е ВМ , т. е. А = ВС для некоторой матрицы С из М .
8.2. Вспомним, что матрица А невырожденная. Значит, матрицы В , С также являются невырожденными.
8.3. По построению, det(B) = det(A) / d . Значит, det(C) = d . Поэтому определитель матрицы С простой. Следовательно, по лемме 10 матрица С простая.
9. Найдём матрицу С по формуле: С = В-1 А .
10. Этим заканчивается алгоритм отщепления. От невырожденной матрицы А отщепляем справа любую возможную простую матрицу С . Произвольность матрицы С следует невырожденности матрицы А , из равенства А = ВС, из произвольности матрицы В .
6. АЛГОРИТМ ОТЩЕПЛЕНИЯ: ВЫЧИСЛЕНИЕ
Далее А - невырожденная матрица из М .
1.1. Выпишем все п столбцов А1, ... , Ап матрицы А .
1.2. Вычислим определитель А = det(A) матрицы А .
2.1. Если полином А степени 0 (константа), то алгоритм заканчиваем (матрица А обратима и раскладывать на множители нечего).
2.2. Далее: полином А не ноль, но не обратим в К .
2.3. Разложим полином А на простые множители.
2.4. Выберем простой сомножитель d полинома А.
3.1. Для всех i, ] е п с помощью алгоритма деления полиномов из К с остатком вычислим остаток а ] при делении полинома а ] (элемент матрицы А в ячейке (/', ])) на полином d . Из этих остатков составим матрицу А (в ячейку (/', ]) этой матрицы помещаем полином а ] ).
— —п
3.2. Методом Гаусса вычислим подпространство всех столбцов X высотой п из К таких, что АХ = 0 (К - это поле К / ^ = 0)).
3.3. Выберем любой ненулевой элемент X из Апп.
3.4. Рассмотрим произвольный столбец X из X е X + d - К".
3.5. Вычислим столбец У = АX - d_1.
По построению У е К" .
4.1. Выпишем столбы А1, ... , Ап, У .
4.2. Элементарными преобразованиями столбцов приведём множество А1, ... , Ап, У столбцов к ступенчатому виду (нули будут над ступеньками).
4.3. Исключим нулевые столбцы.
4.4. Получим множество N из п столбцов.
5. Составим из этих п столбцов матрицу. Обозначим её В .
6. Вычислим матрицу С по формуле: С = В-1 А .
7. Этим заканчивается алгоритм отщепления (от невырожденной матрицы А отщепляем справа любую возможную простую матрицу).
7. ИНВАРИАНТЫ МНОГОКРАТНОГО ОТЩЕПЛЕНИЯ
Применяя теперь алгоритм отщепления к матрице В и т. д., получим возрастающую цепочку АМ с ВМ с... .
В алгебре М возрастающие цепочки идеалов обрываются (лемма 5). Поэтому процесс многократного отщепления рано или поздно закончится. И даже известно число шагов этого процесса: это число делителей полинома det(A) (не обязательно различных).
1. В результате процесса получим плотную цепочку р! с... срт с М правых идеалов алгебры М. Здесь т - число делителей полинома det(M) при его разложении в произведение простых сомножителей (не обязательно различных). Каждому звену цепочки можно приписать простой делитель полинома det(A). Вместо самой цепочки правых идеалов
проще рассматривать цепочку р1К" с... сртКп с К" подмодулей К-модуля К". Эта цепочка также плотная.
2. Для невырожденной матрицы А К-модуль Кп / ( АКп ) является конечномерным ^пространством. На конечномерном ^пространстве действует элемент 5 - как линейное преобразование. Всё сразу сводится к линейной алгебре.
3. Для произвольного элемента f из К обозначим через Vf множество всех элементов
К-модуля К" / (АКп), каждый из которых аннулируется некоторой степенью элемента f.
Множество Vf является К-модулем.
Если полином det(A) является делителем некоторой степени полинома аЬ , то К-модуль
К" / (АКп) есть прямая сумма своих подмодулей Vа, Vb . Поэтому К-модуль К" / (АКп)
разлагается в прямую сумму своих подмодулей Vd , по всем простым делителям d полинома det(А) (здесь делитель d можно брать только один раз).
4. Для разложения на простые множители невырожденной матрицы А нужны следующие данные:
1) плотная цепочка р0 с р1 с... срт сМ правых идеалов алгебры М такая, что
р0 = АМ ;
2) для каждого правого идеала из этой цепочки (не считая р0 и М) нужно выбрать порождающую матрицу.
Этими данными разложение задано однозначно.
Более простые инварианты разложения. Для разложения на простые множители невырожденной матрицы А можно использовать следующие данные:
1) плотная цепочка р0К" с р^К" с... с ртКп с К" правых подмодулей К-модуля К" такая, что р0К" = АК" ;
2) для каждого К-модуля из этой цепочки (не считая р0К", К") нужно выбрать элемент
из Кп , т. е. выбрать порождающее множество не более чем из п элементов (существование этого множества гарантирует лемма 1). Если в этом множестве меньше чем п элементов -
дополнить его до «-элементного любыми элементами из К". Пронумеровать эти п элементов (чтобы из множества получилась п-ка).
Этими данными разложение задано однозначно.
8. ЧАСТНЫЕ СЛУЧАИ
1. Если при разложении полинома на простые сомножители нет одинаковых, то существует совсем простой набор данных разложения матрицы А : разложение на простые множители матрицы А соответствует разложению на простые множители полинома det(A).
2. Своей простотой выделяется более общий случай: когда элемент 5 действует на конечномерное ^-пространство V = К" / (АКп) полупросто, т. е. его жорданова форма диагональна. Эквивалентные инвариантные формулировки этого свойства:
1) любой подмодуль К-модуля V дополняется;
2) для любого простого делителя d полинома det(A), пространство Vd аннулируется полиномом d.
3) конечномерное ^-пространство V есть прямая сумма своих подпространств Vd , где d е D ; D - множество всех простых делителей полинома det(A); Vd - множество элементов из V , аннулируемых элементом d .
3. Если поле k алгебраически замкнуто (например, это поле комплексных чисел), то всё сводится к жордановой форме линейного преобразования, полученного при действии элемента
5 на конечномерное ^пространство К" / АК".
9. АЛГОРИТМ РЕШЕНИЯ уравнения ах + Ьу = d В ПОЛИНОМАХ
Для произвольных полиномов а, Ь из К = ^5] обозначим через d их наибольший общий делитель. Здесь для произвольных а, Ь е К изложен алгоритм нахождения их наибольшего общего делителя d и полиномов х, у из К таких, что ах + Ьу = d , вместе с обоснованием, на базе алгоритма деления полиномов с остатком.
1. Рассмотрим матрицу
при конкретных а, Ь . Например,
Ь = 52 +1:
( 53 52 +1
1 0 0 1
2. Будем рассматривать следующее элементарное преобразование строк: из одной строки можно вычитать другую, предварительно помноженную на полином из К .
С помощью этого элементарного преобразования строк будем уменьшать степени по 5 элементов первой строки до тех пор, пока один из элементов первой строки не станет делителем другого (т. е. если на следующем шаге в первой строке получим 0, то этот шаг делать не надо).
( 53 52 +1
Л
(
-5 52 +1
1 -5^ 01
(к первой строке прибавили вторую, предварительно
помноженную на -5)
(-5 1 1
5 1 - 52
(ко второй строке прибавили первую, предварительно помноженную на 5).
3. В первом столбце получена единица. На следующем шаге в первом столбце должны получить 0. Поэтому элементарные преобразования строк матрицы закончены.
3
а = 5
5
Рассмотрим вектор - строку полученной матрицы, в которой в первом месте стоит
2
делитель первой ячейки другой строки: (1,5,1 - 5 ).
Первый элемент вектора - это наибольший общий делитель d элементов а, Ь ; второй и третий элементы - это полиномы х, у .
Итак, х = 5, у = 1 - 52.
3 2 2
4. Несложно проверить, что ах + Ьу = d : 5 • 5 + (1 + 5 ) • (1 - 5 ) = 1.
5. При желании полином d можно нормировать - вектор можно поделить на старший коэффициент полинома d.
Обоснование алгоритма. Введём обозначения: К3 - множество строчек-троек
2
полиномов из К; V - множество троек (г, х, у) из К3 таких, что г = ах + Ьу; V -множество столбцов высоты 2, элементы которых лежат в V , или множество матриц 2 х 3, у которых строки - элементы из V .
Равенство г = ах + Ьу равносильно тому, что - г + ах + Ьу = 0. Поэтому V есть аннулятор
столбца
V Ь /
Элементарные преобразования строк произвольной матрицы 2 х 3 равносильны умножению этой матрицы справа на матрицу 2 х 2 с единичным определителем (точнее, даже треугольную матрицу с единицами на главной диагонали).
Теперь понятно, что V есть левый М2 (К) -модуль. Значит, при элементарных
2
преобразованиях матриц из V получаем матрицу из ^2 со строками из V.
Поэтому в конце, когда получим тройку (ё, х, у), автоматически будет выполнено равенство d = ха + уЬ .
10. ПРИМЕР
Г5 0^
Рассмотрим n = 2, A = 5 • E = I I. Можно сразу привести два разложения матри-
Г 5 о ^ Г 5 о Y1 о ^ Г 5 о ^ Г1 о Y 5 о ^
цы A : I I = I II |, | | = | || |. Кажется, будто бы все разложения
^ о 51 ^ о 1X о 5 J ^ о 51 ^ о 5 X о 1J
матрицы A этими двумя и ограничиваются, но это не так. Покажем это с помощью алгоритма отщепления.
Рассмотрим произвольное m е n . Тогда AKm = 5Km , Km / (AKm) = Km / (5Km) = km (= - значок изоморфизма). Здесь имеется в виду изоморфизм K-модулей; элемент 5 зануляет (аннулирует) весь K-модуль km .
2 2 2 2 2 Специализируем это для n = 2: AK" = 5K* , K" /(AK ) = k* .
1. Ann = k2 .
— 2
2. Можно (и нужно) выбрать любой ненулевой элемент X из Ann = k . Если выберем
- Г1 ^ - Г о ^
X = II или X = II, то получим оба приведённых выше диагональных разложения
- Г^ -
матрицы A . Выберем X = I I, X = X .
(1 5 0
3. Составим из X и из столбцов матрицы А матрицу: I
4. Элементарными преобразованиями столбцов приведём эту матрицу к ступенчатому
«у: I1 5 0]: I1 0 0]: Р 0 0'
^10 5 ) -5 5) 5 0^
(1 0 ^
5. Исключим нулевой столбец. Полученную матрицу I I обозначим через В .
I1
6. Вычислим определитель матрицы В : det(B) = 5 Ф 0 .
--( 5 0 ^
7. Выпишем для В дополнительную матрицу В : В = I I (для этого в матрице 2 х 2 элементы по диагоняли меняют местами, а у элементов вне диагонали меняют знак).
-1 -1 -1 _ -1 ( 5 0
8. Выпишем матрицу В , обратную к В : В = (det В) • В = 5 -I
V-1 1
1 1 ^ - ( 5 0
9. Вычислим матрицу С = В" 1А : С = В 1А = 5 1В • 5 • Е = В = I ^
Таким образом, получено разложение матрицы А (в произведение простых матриц):
'5 0 ^ (1 0 V 5 0 ^
А = ВС = ВВ, или
V 0 5) V! 5 Д-1 1
Замечание. Полученный пример является частным случаем разложений
ВВ = det(B) • Е, ВВ = det(B) • Е .
Здесь В - любая квадратная матрица; В - дополнительная к ней. Если матрица В размера 2 х 2 и определитель det(B) простой, то матрицы В, В простые, и приведённые выше разложения есть разложения в произведение простых матриц.
ЗАКЛЮЧЕНИЕ
В работе исследован вопрос об единственности разложения невырожденных полиномиальных матриц на простые множители (матрицы п на п над алгеброй полиномов от одной буквы).
СПИСОК ЛИТЕРАТУРЫ
1. Дэвенпорт Г. Высшая арифметика. Введение в теорию чисел. - М.: Наука, Гл. ред. физ.-мат. лит., 1965. -
176 с.
2. Евклид. Начала: пер. с греч. В 3 т. - М.-Л.: ГИТТЛ, 1948-1950. - 507 с. - Т. 1: кн. I-VI; т. 2: кн. VII-IX; т. 3: кн. X-XII.
3. Хрестоматия по истории математики. Алгебра. Теория чисел. Геометрия: учеб. пособие для физ.-мат. фак. пед. ин-в / И.Г. Башмакова, Ю.А. Белый, С.С. Демидов, Б.А. Розенфельд, А.П. Юшкевич. - М.: Просвещение, 1976. -316 с.
4. Бурбаки Н. Очерки по истории математики: пер. с фр. - М.: ИЛ, 1963. - 292 с.
5. Гаусс К.Ф. Труды по теории чисел. - М.: Изд-во АН СССР, 1959. - 979 с.
6. Начала гидростатики. Архимед, Стевин, Галилей, Паскаль / под общ. ред И.И. Агола и др. - М.-Л.: Гос. техн.-теорет. изд-во, 1932. - 413 с.
7. Бурбаки Н. Алгебра. Т. 3. Модули, кольца, формы. - М.: Наука, гл. ред. физ.-мат. лит., 1966. - 554 с.
8. Бурбаки Н. Коммутативная алгебра. - М.: Мир, 1971. - 706 с.
9. Кострикин А.И. Введение в алгебру. Основы алгебры. - М.: Наука, 1994. - 319 с.
10. Алгебра и теория чисел: учеб. пособие для студентов-заочников II курса физ.-мат. факультетов пед. ин-тов / Н.А. Казачек, Г.Н. Перлатов, Н.Я. Виленкин, А.И. Бородин. - М.: Просвещение, 1984. - 193 с.
11. Винберг Е.Б. Алгебра многочленов: учеб. пособие для студентов-заочников III-IV курсов физ.-мат. факультетов пед. ин-тов. - М.: Просвещение, 198о. - 175 с.
12. Куликов Л.Я. Алгебра и теория чисел: учеб. пособие для пед. ин-тов по спец. «Математика», «Математика и физика», «Физика и математика». - М.: Высш. шк., 1979. - 562 с.
13. Anthaklis P.J., Michel A.N. Linaer Systems. - Boston: Birkhauser, 2ооб. - 67о p.
14. Wang Qing-Guo. Decoupling Control. - Berlin Heidelberg: Springer-Vedrlag, 2оо3. - 356 p.
15. Vidya5agar M. Control system sythesis: a factorization approach. Pt. 1. - San Rafael, CA: Morgan & Claypool, 2о 11. - 184 p.
16. Vidya5agar M. Control system sythesis: a factorization approach. Pt. 2. - San Rafael, CA: Morgan & Claypool, 2о 11. - 227 p.
17. Kailath T. Linear Systems. - Englewood Cliffs, New Jork: Prentice-Hall, 198о. - 7о4 p.
18. Wolovich W. Linear Multivariable Control. - Berlin; Heidelberg: Springer-Verlag, 1974. - 368 p.
19. Воевода А.А. Матричные передаточные функции (основные понятия): учеб. пособие. - Новосибирск: изд-во НГТУ, 1994. - 95 c.
20. Воевода А.А. Матричные передаточные функции (синтез): конспект лекций. - Новосибирск: изд-во НГТУ, 1995. - 95 c.
21. MacDuffee C.C. The Theory of Matrices. - Berlin: Springer, 1933. - 119 p.
22. Cahen E. Theory des nombres. Vol. 1. - Paris: Librairie scintifique A. Hermann & Fils, 1914. - 425 p.
Корюкин Анатолий Николаевич, кандидат физико-математических наук, старший научный сотрудник Института математики им. Соболева. Основные направления научных исследований - теория ассоциативных некоммутативных колец, математические основы теории автоматического управления. Имеет более 4о публикаций. E-mail: koryukin@sibmail.ru
The decomposition on simple multipliers of nondegeneratepolinomialny matrixes
A.N. KORYUKIN1
1 Sobolev'5 In5titute of mathematic5, 4 pro5pect Koptyuga, Novo5ibir5k, 630090, Russian Federation, candidate of phy5ical and mathematical 5cience5, 5enior re5earcher, e-mail: koryukin@5ibmail.ru
Object of research. In algebra of polinomialny matrixes (over algebra of K of polynoms from one letter) the size n on n the question of uniqueness of decomposition on simple multipliers is studied. Urgency. The theory of divisibility is a classical part of the mathematics which has entered into process of training in higher educational institutions. But still the question of uniqueness of decomposition on simple multipliers was studied only for commutative algebras, and only in the absence of zero dividers. Meanwhile in the linear theory of automatic control (multichannel systems) long ago concepts of the greatest left and right dividers of polinomialny matrixes that can form a basis of the theory of divisibility for polinomialny matrixes also are successfully used. The most important question of the theory of divisibility - existence and uniqueness of decomposition on simple multipliers. The question of uniqueness of decomposition on simple multipliers is complicated by that in a ring of matrixes there is no commutativity.. Results. It is shown that in a ring of polinomialny matrixes the left and right ideals main (matrixes over algebra of polynoms from one letter). From this it follows that for any number of nondegenerate matrixes there is the greatest left divider, the greatest right divider, the smallest left multiple, smallest right multiple. It is shown that for any nondegenerate matrix there is a decomposition on simple multipliers. Indecomposable matrixes are characterized: it in accuracy of a matrix with simple (indecomposable) determinant. The problem of uniqueness of decomposition is solved on simple multipliers so: for any nondegenerate matrix of A all simple right dividers (simple matrixes of B for which there is a matrix the C such that A = BC) are described. This description is made in terms of a lattice of K-submodules n-ok polynoms of K. The decomposition of A = BC allows to find all decomposition of a nondegenerate matrix of A on simple multipliers. These searches are written down in the form of two algorithms (for mathematicians and for engineers). Some invariants of decomposition of a matrix on simple multipliers are described. Some special cases of decomposition are described. One example of decomposition is given. Technique and tools. Theory of rings. Divisibility theory.
Keywords: Polinomialny matrixes, simple multipliers, indecomposable elements, decomposition on simple multipliers, uniqueness of decomposition on simple multipliers, the theorem of uniqueness of decomposition on the simple multipliers, the greatest general left divider, the greatest general right divider
* Received 01 July 2013.
REFERENCES
1. Davenport H. The higher arithmetic. An Introduction to the theory of numbers. New Jork, Harper & Brothers, 1960. (Russ. ed.: Devenport G. Vysshaya arifmetika. Vvedenie v teoriyu chisel. Moskow, Nauka Publ., Glavnaya redaktsiya fiziko-matematicheskoi literatury, 1965. 176 p.).
2. Evklid. Nachala. V 3 t. Perevod s grecheskogo [Beginnings. Translation from Greek. In 3 vol.]. Moscow, Leningrad, GITTL Publ., 1948-1950. 507 p. Vol. I: books of I-VI. Vol. II: books of VII-IX. Vol. III: books X-XII.
3. Bashmakova I.G., Belyi Yu.A., Demidov S.S., Rozenfeld B.A., Yushkevich A.P. Khrestomatiya po istorii matematiki. Algebra. Teoriya chisel. Geometriya. [Anthology of mathematics history. Algebra. Theory of numbers. Geometry]. Moscow, Prosveshchenie Publ., 1976. 316 p.
4. Bourbaki N. Elements of the History of Mathematics. Translation from French. 1960 (Russ. ed.: Burbaki N. Ocherki po istorii matematiki. Moscow, IL Publ., 1963. 292 p.).
5. Gauss K.F. Untersuchungen über höhere Arithmetik. Translated from Latin. Springer, 1889. (Russ. ed.: Gauss K.F. Trudy po teorii chisel. Moscow, Academy of Sciences ofthe USSR Publ., 1959. 979 p.).
6. Agola I.I., ed., et al. Nachala gidrostatiki. Arkhimed, Stevin, Galilei, PaskaV [Began hydrostatics. Archimedes, Galilei, Pascal, Stevin]. Moscow, Leningrad, Gosudarstvennoe tekhniko-teoreticheskoe izdateltsvo, 1932. 413 p.
7. Bourbaki N. Algebre. Vol. 3. Translation from French], Paris, Heimann, 1962. (Russ. ed.: Burbaki N. Algebra. T. 3. Moduli, kol'tsa, formy. Moscow, Nauka Publ., Glavnaya redaktsiya fiziko-matematicheskoi literatury, 1966. 554 p.).
8. Bourbaki N. Algebre commutative (Russ. ed.: Burbaki N. Kommutativnaya algebra. Moscow, Mir Publ., 1971. 706 p.).
9. Kostrikin A.I. Vvedenie v algebru. Osnovy algebry [Introduction in algebra. Algebra bases]. Moscow, Nauka Publ., 1994. 319 p.
10. Kazachek N.A., Perlatov G.N., Vilenkin N.Ya., Borodin A.I. Algebra i teoriya chisel [Algebra and theory of numbers]. Moscow, Prosveschenie Publ., 1984. 193 p.
11. Vinberg E.B. Algebra mnogochlenov [Algebra of multinomials]. Moscow, Prosveschenie Publ., 1980. 175 p.
12. Kulikov L.Ya. Algebra i teoriya chisel [Algebra and theory of numbers]. Moscow, Vysshaya shkola Publ., 1979.
562 p.
13. Anthaklis Panos J., Michel Antjony N. Linaer Systems. Boston, Birkhauser, 2006. 670 p.
14. Wang Qing-Guo. Decoupling Control. Berlin, Heidelberg, Springer-Vedrlag, 2003. 356 p.
15. Vidyasagar M. Control system sythesis: a factorization approach. Pt. I. San Rafael, CA, Morgan & Claypool, 2011.
184 p.
16. Vidyasagar M. Control system sythesis: a factorization approach. Pt. II. San Rafael, CA, Morgan & Claypool, 2011.
227 p.
17. Kailath T. Linear Systems. Englewood Cliffs, New Jork, Prentice-Hall, 1980. 704 p.
18. Wolovich W. Linear Multivariable Control. Berlin, Heidelberg, Springer-Verlag, 1974. 368 p.
19. Voevoda A.A. Matrichnye peredatochnye funktsii (osnovnye ponyatiya) [Matrix transfer functions (the main concepts)]. Novosibirsk, NSTU Publ., 1994. 95 p.
20. Voevoda A.A. Matrichnye peredatochnye funktsii (sintez) [Matrix transfer functions (synthesis)]. Novosibirsk, NSTU Publ., 1995. 95 p.
21. MacDuffee C.C. The Theory of Matrices. Berlin, Springer, 1933. 119 p.
22. Cahen E. Theory des nombres. Vol. 1. Paris, Librairie scintifique A. Hermann & Fils, 1914. 425 p.
ISSN 1814-1196, http://journals.nstu.ru/vestnik Scientific Bulletin of NSTU Vol. 56, No. 3, 2014, pp. 23-36
