РУБРИКИ |
Диплом: Многочленные матрицы |
РЕКОМЕНДУЕМ |
|
Диплом: Многочленные матрицыДиплом: Многочленные матрицыСодержание Введение .......................... 3 Глава I. Многочленные матрицы. §1. Элементарные преобразования многочленной матрицы... 4 §2. Канонический вид λ-матрицы.............. 5 §3. Наибольшие общие делители миноров.......... 10 §4. Условия эквивалентности λ-матриц........... 15 §5. Элементарные делители многочленной матрицы.... 19 Глава II. Матричные многочлены. §1. Деление матричных многочленов. Обобщенная теорема Безу. 23 §2. Скалярная эквивалентность.............. 28 §3. Характеристический многочлен матрицы........ 29 §4. Минимальный многочлен матрицы.......... 31 §5. Критерий подобия матриц............... 33 §6. Нормальная форма Жордана.............. 36 Глава III. Функции от матриц. §1. Многочлен от жордановой матрицы........... 40 §2. Скалярные функции................. 42 §3. Представление значений функций многочленами.... 46 §4. Элементарные делители функций........... 48 §5. Степенные ряды.................. 49 Литература ........................ 53 Введение Дипломная работа состоит из трех глав. Первые два параграфа I главы посвящены изучению свойств многочленных матриц. В §3, этой же главы, вводятся понятия наибольших общих делителей миноров и инвариантных множителей многочленной матрицы. На основе этого в последующих двух параграфах рассматриваются условия эквивалентности λ- матриц и строится аналитическая теория элементарных делителей. С каждой квадратной матрицей связаны два многочлена: характеристический и минимальный. Эти многочлены играют большую роль в различных вопросах теории матриц. Во II главе рассматриваются свойства характеристического и минимального многочлена. Этому исследованию предпосылаются основные сведения о многочленах с матричными коэффициентами и о действиях над ними. Последний параграф II главы посвящен нормальной форме Жордана. В главе III, посвященной функциям от матрицы, рассмотрены вопросы матричного исчисления, для решения которых используется возможность приведения матриц к нормальной жордановой форме. В §5, этой же главы, рассмотрен вопрос о сходимости степенного ряда от матрицы А. Приведенный материал иллюстрируется в решениях различных примеров. Глава I. Многочленные матрицы. §1. Элементарные преобразования многочленной матрицы. Многочленной матрицей или λ-матрицей называется прямоугольная (в частности, квадратная) матрица А (λ) = ║аі ј (λ) ║, где i=1, 2, ., m; j=1, 2, ., n; элементы которой являются многочленами от одного переменного λ с числовыми коэффициентами из основного поля К. Элементарными преобразованиями λ-матрицы А(λ) называются преобразования следующих типов: I. Перестановка двух строк. II. Умножение строки на число с Є К, с ≠ 0. III. Прибавление к одной строке другой строки, умноженной на любой многочлен f (λ), и аналогичные преобразования столбцов. Элементарные преобразования I, II, III равносильны умножению многочленной матрицы А(λ) слева соответственно на следующие квадратные матрицы порядка m: (i) (j) (i) (j)
1 0 1 0 1 0 S' = с , S" = 1. b(λ) ..(i), S"' = 0 1 , 0 1 0 1 1 0 0 1 т.е. в результате применения преобразований I, II, III матрица А(λ) преобразуется соответственно в матрицы S'· А(λ), S"· А(λ), S"'· А(λ). Поэтому преобразования типа I, II, III называются левыми элементарными операциями. Совершенно аналогично определяются правые элементарные операции над многочленной матрицей (эти операции производятся не над строками, а над столбцами) и соответствующие им матрицы порядка n:
1 0 1 0 1 0 Т' = с (i), Т" = 1.....(i), Т"' = 0.1....(i) b(λ)....(j) 1.0....(j) . 0 1 0 1 0 1 В результате применения правой элементарной операции матрица А(λ) умножается справа на соответствующую матрицу Т. Заметим, что матрица Т' совпадает с матрицей S', а матрицы Т", Т"' совпадают с матрицами S", S"', если в последних поменять местами индексы i и j. Матрицы типа S', S", S"' (или, что то же, типа Т', Т", Т"') называются элементарными. Две λ-матрицы А(λ) и B(λ) одинаковых размеров m x n называются эквивалентными, А(λ) ~ B(λ), если от матрицы А(λ) к B(λ) можно перейти при помощи цепочки из конечного числа элементарных преобразований. Отношение эквивалентности обладает тремя основными свойствами: 1) рефлексивность: каждая матрица эквивалентна сама себе А(λ) ~ B(λ); 2) симметрия: если А(λ) ~ B(λ), то B(λ) ~ А(λ); 3) транзитивность: если А(λ) ~ B(λ), и B(λ) ~ С(λ), то А(λ) ~ С(λ). §2. Канонический вид λ-матрицы Выше было показано, что отношение эквивалентности транзитивно, симметрично и рефлексивно. Отсюда следует, что совокупность всех λ-матриц данных размеров m x n разбивается на непересекающиеся классы эквивалентных матриц, т.е. на такие классы, что любые две матрицы из одного класса эквивалентны, а из разных классов - не эквивалентны между собой. Возникает вопрос о канонической форме λ-матрицы, характеризующей данный класс эквивалентных λ-матриц. Канонической диагональной λ-матрицей размеров m x n называется λ-матрица, у которой на главной диагонали стоят многочлены Е1(λ), Е2(λ), ., Ер(λ), где р - меньшее из чисел m и n, причем не равные нулю среди этих многочленов имеют старшие коэффициенты, равные единице, и каждый следующий многочлен делится на предыдущий, все же элементы вне главной диагонали равны нулю. Т е о р е м а 1. Всякая λ-матрица конечным числом элементарных преобразований может быть приведена к канонической диагональной форме. Доказательство. Пусть А(λ) - прямоугольная многочленная матрица. Применяя к А(λ) как левые, так и правые элементарные операции приведем к канонической диагональной форме. Среди всех не равных нулю элементов аіј(λ) матрицы А(λ) возьмем тот элемент, который имеет наименьшую степень относительно λ, и путем соответствующей перестановки строк и столбцов сделаем его элементом а11(λ). После этого найдем частные и остатки от деления многочленов аі1(λ) и а1ј(λ) на а11(λ): аі1(λ) = а11(λ) qі1(λ) + rі1 (λ), а1ј(λ) = а11(λ) q1ј(λ) + r1ј(λ) (i = 2, 3, ., m; j = 2, 3, ., n). Если хотя бы один из остатков rі1(λ), r1ј(λ) (i = 2, ., m; j = 2, ., n), например r1ј (λ), не равен тождественно нулю, то, вычитая из j-го столбца первый столбец, предварительно помноженный на q1ј(λ), мы заменим элемент а1ј(λ) остатком r1ј(λ), который имеет меньшую степень, нежели а11(λ). Тогда мы имеем возможность снова уменьшить степень элемента, стоящего в левом верхнем углу матрицы, поместив на это место элемент с наименьшей степенью относительно λ. Если же все остатки r21(λ), . rm1(λ); r12(λ), ., r1n(λ) равны тождественно нулю, то, вычитая из i-ой строки первую, помноженную предварительно на qі1(λ) (i = 2, ., m), а из j-го столбца - первый, предварительно помноженный на q1ј(λ) (j = 2, ., n), мы приведем нашу матрицу к виду а11(λ) 0 . 0 0 а22(λ) . а2n(λ) .......... . 0 аm2(λ) . аmn(λ) Если при этом хотя бы один из элементов аіј(λ) (i = 2, ., m; j = 2, ., n) не делится без остатка на а11(λ), то, прибавляя к первому столбцу тот столбец, который содержит этот элемент, мы придем к предыдущему случаю и, следовательно, снова сможем заменить элемент а11(λ) многочленом меньшей степени. Поскольку первоначальный элемент а11(λ) имел определенную степень и процесс уменьшения этой степени не может неограниченно продолжаться, то после конечного числа элементарных операций мы должны получить матрицу вида а1(λ) 0 . 0 (*) 0 b22(λ) . b 2n(λ) .......... , 0 bm2 (λ) .bmn (λ) в которой все элементы bіј(λ) делятся без остатка на а1(λ). Если среди этих элементов bіј(λ) имеются не равные тождественно нулю, то продолжая тот же процесс приведения для строк с номерами 2, ., m и столбцов с номерами 2, ., n, мы матрицу (*) приведем к виду а1 (λ) 0 0 . 0 0 а2(λ) 0 . 0 0 0 с33(λ) . с3n(λ) , .......... 0 0 сm3(λ) . сmn(λ) где а2(λ) делится без остатка на а1(λ), а все многочлены сіј(λ) делятся без остатка на а2(λ). Продолжая этот процесс далее, мы в конце концов придем к матрице вида а1(λ) 0 . 0 0 . 0 0 а2(λ). 0 0 . 0 .............. 0 0 .аs(λ) 0 . 0 , 0 0 . 0 0 . 0 ........... 0 0 . 0 0 . 0 где многочлены а1(λ), а2(λ), ., аs(λ) (s ≤ m, n) не равны тождественно нулю и каждый из них делится на предыдущий без остатка. Помножая первые s строк на соответствующие отличные от нуля числовые множители, мы сможем добиться того, чтобы старшие коэффициенты многочленов а1(λ), а2(λ), ., аs(λ) были равны единице. Таким образом мы доказали, что произвольная прямоугольная многочленная матрица А(λ) эквивалентна некоторой канонической диагональной. Данная теорема утверждает, что каждый класс эквивалентных матриц содержит, по меньшей мере, одну матрицу, имеющую каноническую диагональную форму. Пример Найти каноническую форму и инвариантные множители l-матрицы l l-1 l-1 l-2 l l3+l l-1 l3+l-1 А(l) = l+1 l3+l+2 l+1 l3+2l+3 l-1 l-3 l-3 -6 Третий столбец вычитаем из остальных: 1 0 l-1 -1 1 l3+1 l-1 l3 А(l) ~ 0 l3+1 l+1 l3+l+2 ~ 2 0 l-3 -l-3 (из второй и четвертой строк вычитаем первую, умноженную соответственно на 1 и 2) 1 0 l-1 -1 ~ 0 l3+1 0 l3 +1 ~ 0 l3+1 l+1 l3+l+2 0 0 -l-1 -l-1 (из III и IV столбцов вычитаем I, умноженный соответственно на l-1 и -1) 1 0 0 0 0 l3+1 0 l3+1 ~ 0 l3+1 l+1 l3 +l+2 ~ 0 0 -l-1 -l-1 (II и III столбцы вычитаем из IV) 1 0 0 0 0 l3+1 0 0 ~ 0 l3+1 l+1 0 ~ 0 0 -l-1 0 (из III строки вычитаем II и полученную строку прибавляем к IV) 1 0 0 0 0 l3+1 0 0 ~ 0 0 l+1 0 ~ 0 0 0 0 (переставляем II и III строки, а также второй и третий столбцы) 1 0 0 0 0 l+1 0 0 ~ 0 0 l3+1 0 . 0 0 0 0 Это каноническая матрица; значит матрица А(l) имеет инвариантные множители: Е 1(l)=1, Е2(l)=l+1, Е3(l)=l3+1, Е4 (l)=0. §3. Наибольшие общие делители миноров Пусть F - какая-нибудь λ-матрица порядка n. Составим ее всевозможные миноры порядка к (к = 1, 2, ., n). Эти миноры являются многочленами от λ. Обозначим их наибольший общий делитель через Dк(λ) (наибольшим общим делителем мы условимся называть общий делитель наивысшей степени со старшим коэффициентом 1. Поэтому все не равные нулю многочлены Dк(λ) имеют старший коэффициент 1). Если окажется, что все миноры к-го порядка равны нулю, то по определению будем считать Dк(λ)=0. В частности, D1(λ) есть наибольший общий делитель элементов матрицы F, а Dn(λ) равен определителю F, деленному на свой старший коэффициент. Т е о р е м а 2. Эквивалентные λ-матрицы имеют один и тот же наибольший общий делитель миноров к-го порядка (к = 1, 2, ., n). Доказательство. Пусть F1, F2 - две эквивалентные λ-матрицы. Обозначим наибольшие общие делители их миноров к-го порядка соответственно через Dк1(λ) и Dк2(λ). Требуется показать, что Dк1(λ) = Dк2(λ). Нам известно, что F2 можно получить из F1 цепочкой элементарных преобразований. Предположим сначала, что эта цепочка состоит только из одного элементарного преобразования. Пусть, например, F2 получается из F1 умножением всех элементов i-й строки матрицы F1 на число α≠0. Соответственные миноры F1 и F2 тогда либо совсем не отличаются друг от друга, либо отличаются лишь постоянным множителем α. Однако постоянный множитель не влияет на вычисление наибольшего общего делителя многочленов, поэтому Dк1 = Dк2. То же самое будет и в случае, когда F2 получается из F1 умножением на α элементов какого-либо столбца матрицы F1. Пусть теперь F2 получается из F1 одним из элементарных преобразований; например, пусть F2 возникает в результате прибавления к i-й строке матрицы F1 элементов j-й строки, умноженных на f(λ). Покажем, что Dк2 делится на Dк1. В самом деле, миноры к-го порядка матриц F1, F2 можно разбить на три класса. К первому мы отнесем те из них, которые не содержат элементов i-й строки. Соответственные миноры матриц F1, F2 в этом случае, очевидно, равны друг другу. Ко второму классу отнесем те миноры, которые содержат элементы и i-й и j-й строк. Эти миноры у матриц F1, F2 будут снова равными, так как величина определителя не меняется, если к элементам одной из его строк прибавляются величины пропорциональные элементам какой-либо другой строки. Наконец, к третьему классу мы отнесем миноры, содержащие элементы i- й строки и не содержащие элементов j-й строки. Соответственные миноры этого класса имеют вид: ..... ............. м1 = f iυ1 . f iυk , м2 = f iυ1 + f jυ1 ∙ f(λ) ... f iυk + f jυk ∙ f(λ) , ..... ............. где невыписанные строки у обоих миноров одинаковы. На основании теоремы сложения определителей ..... ...... м2 = f iυ1 . f iυk + f(λ) f jυ1 ... f jυk = м1 ± f(λ)N1 , ..... ....... где N1 - некоторый минор матрицы F1. Все миноры к-го порядка матрицы F1 делятся и на Dк1. Из наших рассуждений видно, что на Dк1 делятся и все миноры к-го порядка матрицы F2, так как они либо совпадают с соответственными минорами матрицы F1, либо выражаются через них линейно. Но в таком случае Dк1 войдет множителем в наибольший общий делитель миноров к- го порядка матрицы F2, т.е. войдет множителем в Dк2. Итак, если F2 получается из F1 одним элементарным преобразованием, то Dк2 делится на Dк1. Совершая над F2 обратное элементарное преобразование, мы получим F1. Поэтому, Dк1 также должен делиться на Dк2. Но если старшие коэффициенты двух многочленов равны и эти многочлены делятся без остатка друг на друга, то они совпадают. Таким образом, Dк1 = Dк2. Пока доказано равенство наибольших общих делителей в предположении, что F2 получается из F1 одним элементарным преобразованием. Однако если Dк(λ) не меняется при каждом отдельном элементарном преобразовании, то, очевидно, Dк(λ) не изменится и при нескольких последовательных преобразованиях. Потому теорема доказана в общем виде. Вычислим многочлены D1(λ), ., Dn(λ) для матрицы, имеющей канонический диагональный вид . d1(λ) d2(λ) . D = . . . dn(λ) Чтобы получить какой-нибудь минор к-го порядка, мы должны из D вычеркнуть n- k строк и n-k столбцов. Если из D вычеркнуть i-ю строку, то в i-м столбце останутся только нули. Поэтому, чтобы получить минор, отличный от нуля, мы должны вычеркнуть все столбцы матрицы D, номера которых равны номерам вычеркнутых строк. Таким образом, отличные от нуля миноры к-й степени должны иметь вид dυ1(λ) dυ2 (λ) . = dυ1 (λ) dυ2 (λ) . dυk(λ). (1) . . dυk(λ) Наибольшим общим делителем этих миноров будет Dк(λ). Из неравенств 1 ≤ υ1 ≤ υ2 ≤ . ≤ υk ≤ n следует, что 1 ≤ υ1, 2 ≤ υ2, . , k ≤ υk. Поэтому dυi(λ) делится на di(λ), а значит, dυ1(λ) . dυk(λ) делится на d1(λ) . dk(λ). Мы видим, следовательно, что все миноры к-го порядка матрицы D делятся на минор d1(λ) . . = d1(λ) . d k(λ) (2) . dk(λ) Если этот минор равен нулю, то и все миноры к-го порядка матрицы D также равны нулю. Согласно определению в этом случае Dк(λ)=0. Если минор (2) отличен от нуля, то многочлены d1(λ), ., d k(λ) отличны от нуля и имеют старший коэффициент 1. Но тогда и минор (2) имеет старший коэффициент 1. Поскольку все миноры (1) делятся на (2), то Dк(λ) совпадает с (2). Следовательно, в обоих случаях имеем Dк(λ) = d1(λ)d2(λ). d k(λ) (к = 1, 2, ., n) (3) Таковы многочлены Dк(λ) канонической диагональной матрицы с диагональными элементами d1(λ), . , d n(λ). Рассмотрим теперь произвольную λ-матрицу F. Обозначим через Dк(λ) наибольший общий делитель миноров степени к этой матрицы. Согласно теореме 1. матрицу F элементарными преобразованиями можно привести к канонической диагональной форме d1(λ) . D = . . . dn(λ) Согласно теореме 1. многочлены Dк(λ), вычисленные для матрицы D, совпадают с соответственными многочленами Dк(λ), вычисленными для F. Таким образом, многочлены Dк(λ) матрицы F и диагональные элементы канонической диагональной матрицы D, к которой можно привести F связаны соотношениями (3). Пусть D1(λ), . , Dr(λ) отличны от нуля, а остальные многочлены Dr+1(λ), . , Dn(λ), если они есть, равны нулю. Тогда из (3) имеем D1(λ) = d1(λ), d1(λ) = D1(λ), D2(λ) = d1(λ)d2(λ), d2(λ) = D2(λ) : D1(λ), ........ ......... Dr(λ) = d1(λ)d2(λ) . dr(λ), dr(λ) = Dr(λ) : Dr-1(λ), Dr+1(λ) = d1(λ)d2(λ) . dr(λ)dr+1(λ), dr+1(λ) = Dr+1(λ) : Dr(λ). Поскольку dr+1(λ)=0, то dr+2(λ), . , dn(λ) также должны быть равны нулю, и мы имеем окончательно d1(λ) = D1(λ), d2(λ) = D2(λ) : D1(λ), . , dr(λ) = Dr(λ) : Dr-1(λ), dr+1(λ) = . = dn(λ) = 0 (4) Тем самым мы получили следующую теорему: Т е о р е м а 3. Если наибольшие общие делители Dк(λ) миноров порядка к λ-матрицы F при к=1, 2, . r отличны от нуля, а Dr+1(λ)=0, то диагональные элементы dк(λ) канонической диагональной матрицы, к которой F приводится элементарными преобразованиями, выражаются через Dк(λ) по формулам (4) и определяются, таким образом, матрицей F однозначно. Многочлены d1(λ), . , dn(λ) называются инвариантными множителями матрицы F. Число r, в равенствах (4) это - ранг матрицы F. Действительно, ранг матрицы F есть порядок наивысшего минора F,отличного от нуля. Если этот порядок равен r, то, следовательно, Dr(λ)≠0, а Dr+1(λ)=0. Обратно, условия Dr(λ)≠0, Dr+1(λ)=0 означают, что некоторый минор порядка r отличен от нуля, а все миноры порядка r+1 равны нулю. Следовательно, ранг F равен r. §4. Условия эквивалентности λ-матриц. Первое условие эквивалентности. Для того чтобы многочленные матрицы порядка n были эквивалентны , необходимо и достаточно, чтобы наибольшие общие делители их миноров к-го порядка совпадали при к=1, 2, . , n. Поскольку совпадение наибольших общих делителей миноров равносильно совпадению соответствующих инвариантных множителей, то первое условие эквивалентности можно сформулировать и в следующем виде: для эквивалентности λ-матриц необходимо и достаточно, чтобы их соответствующие инвариантные множители были равны. Доказательство очевидно. В самом деле, если две λ-матрицы F, G эквивалентны, то их наибольшие общие делители Dк(λ) одинаковы (теорема 2). Обратно, если многочлены Dк(λ) у F и G равны, то F и G элементарными преобразованиями приводятся к одной и той же канонической диагональной матрице (теорема 3). Но две матрицы, эквивалентные третьей эквивалентны между собой; следовательно, F эквивалентна G, что и требовалось. Второе условие эквивалентности. Для того чтобы многочленные матрицы F и G порядка n были эквивалентны, необходимо и достаточно, чтобы они удовлетворяли соотношению G = PFQ, где P, Q - некоторые многочленные матрицы с постоянными отличными от нуля определителями. Прежде чем перейти к доказательству этого утверждения, сделаем несколько замечаний. Пусть 1 . . . 1 А = α i-я строка, 1 . . . 1 где α- некоторое число, отличное от нуля. Умножая произвольную матрицу F на A слева, мы увидим, что все элементы матрицы F останутся без изменения, кроме элементов i-й строки которые умножаются на α. Таким образом, каждое элементарное преобразование типа II, совершаемое над матрицей F, равносильно умножению F на подходящую матрицу А слева. Аналогично, если умножить матрицу F слева на матрицу 1. 0 . 0 . 0 0 .1. f(λ).0 i-я строка В = 0. 0 . 1 . 0 j-я строка, 0 . 0 . 0 . 1 где все диагональные элементы равны единице, элемент, стоящий в i-й строке и j-м столбце, равен f(λ), а остальные элементы равны нулю, то в результате к элементам i-й строки F прибавятся элементы ее j-й строки, умноженные на f(λ). Следовательно, каждое элементарное преобразование типа III равносильно умножению F слева на соответственную матрицу В. Наконец, тем же способом легко убедиться, что элементарные преобразования матрицы F типа II, III столбцов равносильны умножению F на соответственные матрицы А, В справа. Перейдем теперь к доказательству второго условия эквивалентности. Необходимость. Пусть матрица G эквивалентна матрице F. Это означает, что G может быть получена из F цепочкой последовательных элементарных преобразований. Каждое элементарное преобразование мы можем заменить умножением на матрицу вида А или В соответственно слева или справа. В результате получим равенство G = P1P2 . Pp FQ1Q2 . Qq, (5) где каждая из матриц Pi , Qj имеет вид А или В (i, j = 1, 2, .). Положим Р = P1P2 . Pp , Q = Q1Q2 . Qq. Поскольку определители матриц В равны единице, а определители матриц А являются постоянными, отличными от нуля числами то определители матриц P и Q будут также постоянными, отличными от нуля числами. Соотношение (5) дает G = PFQ, и необходимость доказана. Достаточность. Допустим обратно, что G = PFQ, (6) где P и Q - многочленные матрицы с постоянными, отличными от нуля определителями. Наибольший общий делитель Dn(λ) всех миноров порядка n матрицы Р равен определителю Р, деленному на его старший коэффициент. Так как этот определитель есть постоянное число, то Dn(λ)=1. Из соотношения (3) параграфа 3 при к=n получаем Dn(λ) = d1(λ)d2(λ) . dn(λ) = 1, откуда d1(λ) = d2(λ) = . = dn(λ) = 1, где d1(λ), . dn(λ) - инвариантные множители матрицы Р. Но инвариантные множители единичной матрицы E также все равны единице, ибо Е имеет канонический диагональный вид. Согласно первому признаку эквивалентности отсюда следует, что матрица Р эквивалентна Е и, значит, может быть получена из Е цепочкой элементарных преобразований. Каждое элементарное преобразование можно заменить умножением на матрицу типа А или В. В результате Р будет представлена следующим образом: Р = P1 . Pк Е Q1 . Qℓ = P1 . Pк Q1 . Qℓ где Pi, Qj - матрицы типов А, В. Применяя те же самые рассуждения к матрице Q, получим аналогичное разложение Q = M1 . MsN1 .Nt. Подставляя эти разложения в (6) придем к равенству G = P1 . Pк Q1 . QℓFM1 . MsN1 .Nt, (7) Из которого видно, что G получается последовательным умножением матрицы F на матрицы Pi, Qi , Mi, Ni типа А или В. Но каждое такое умножение равносильно некоторому элементарному преобразованию. Следовательно, G эквивалентна F, и доказательство закончено. §5. Элементарные делители многочленной матрицы. Рассмотрим произвольную λ-матрицу F, элементами которой являются многочлены от λ с коэффициентами из основного поля К. Обозначим через d1(λ), d2(λ), . , dn(λ) инвариантные множители матрицы F. Часть этих множителей может равняться нулю, поэтому предположим для определенности, что d1(λ), ., dr(λ) нулю не равны, а dr+1(λ)= .=dn(λ)= 0. Число r, есть ранг матрицы F. Разложим каждый из многочленов d1(λ), ., dr(λ) на множители, неприводимые в поле К. Пусть, например, di(λ) = [ε1(λ)]ⁿ¹ [ε2(λ)]ⁿ² . [εk(λ)]ⁿk , где ε1(λ), ., εk(λ) -различные неприводимые многочлены со старшим коэффициентом 1. Выражения [ε1(λ)]ⁿ¹, ., [εk(λ)]ⁿk ,называются элементарными делителями инвариантного множителя di(λ). Элементарные делители всех непостоянных инвариантных множителей матрицы F называются ее элементарными делителями. Например, если инвариантные множители матрицы F равны соответственно 1, λ, λ²(λ+1), λ²(λ+1)², то элементарные делители будут λ, λ², λ², λ+1, (λ+1)². Элементарный делитель вида [ε(λ)]k, где ε(λ) неприводимый многочлен, называют принадлежащим многочлену ε(λ). В рассмотренном примере элементарные делители λ, λ², λ²принадлежат λ, а λ+1 и (λ+1)² принадлежат λ+1. Т е о р е м а 4. Порядок, ранг и система элементарных делителей λ-матрицы F вполне определяют ее инвариантные множители и, следовательно, определяют F с точностью до эквивалентности. Доказательство легко уясняется из следующего примера. Пусть порядок F равен 6, ранг 4, элементарные делители λ, λ², λ², λ+1, (λ+1)³, λ-1, λ-1. Поскольку порядок F есть 6, то F имеет шесть инвариантных множителей d1(λ), ., d6(λ). Из них d5(λ) = d6(λ) = 0, так как ранг F должен быть 4. Если разложить d1(λ), ., d4(λ) на множители, то должны получиться указанные семь элементарных делителей. Поскольку, однако, d4(λ) делится на d3(λ), d2(λ) и d1(λ), то в d4(λ) входят элементарные делители, принадлежащие ко всем неприводимым многочленам, и притом в высших степенях. Поэтому d4(λ) = λ²(λ+1)³(λ-1). Среди оставшихся элементарных делителей λ, λ², λ+1, λ-1 высшие должны войти в d3(λ); следовательно, d3(λ) = λ²(λ+1)(λ-1). В свою очередь высшие из оставшихся должны составить d2(λ), т.е. d2(λ)=λ. Поскольку все элементарные делители уже распределены, то d1(λ)=1. Ясно, что этот способ применим в любом случае, что и доказывает теорему. Элементарные делители зависят от основного поля К. Например, пусть инвариантные множители некоторой λ-матрицы F равны λ²+1, (λ²+1)². Если основное поле есть поле вещественных чисел, то многочлен λ²+1 неприводим и элементарными делителями матрицы F являются λ²+1 и (λ²+1)². Однако, если основное поле - поле всех комплексных чисел, то λ²+1= (λ - i)(λ + i) и элементарными делителями F будут λ+i, (λ+i)², λ-i, (λ-i)². Чтобы получить элементарные делители λ-матрицы, имеющей каноническую диагональную форму, достаточно взять, согласно определению, все элементарные делители ее диагональных элементов. Это же правило имеет место и для произвольной диагональной λ-матрицы. Лемма. Система элементарных делителей произвольной диагональной λ-матрицы F есть объединение элементарных делителей ее диагональных элементов. Примеры: 1) Найдем элементарные делители l-матрицы 1 0 0 0 0 l2+2 0 0 А(l) = 0 0 (l4-4)2 0 0 0 0 0 над полями рациональных, действительных и комплексных чисел. Матрица А(l) имеет каноническую форму. Поэтому на ее главной диагонали стоят инвариантные множители Е1(l)=1, Е2(l)=l2+2, Е 3(l)=(l4-4)2, Е4(l)=0. Разлагая Е2(l) и Е3(l) на неприводимые множители над каждым из указанных полей, получим системы элементарных делителей: 1. над полем рациональных чисел: е1(l)=l2+2, е2(l)=(l2+2)2, е3(l)=(l2-2)2 ; 2. над полем действительных чисел: е1(l)=l2+2, е2(l)=(l2+2)2, е3(l)=(l+Ö2)2, е4(l)=(l-Ö2)2 ; 3. над полем комплексных чисел е1(l)=l+2i, е2(l)=l-2i, е3(l)=(l+2i)2, е4(l)=(l-2i)2, е5(l)=(l+Ö2)2, е6(l)=(l-Ö2)2 . 2) Найти инвариантные множители l-матрицы А(l) размеров 6х8 и ранга 5, если даны ее элементарные делители l+1, l+1, (l+1)3, (l-1)2, (l-1)2. Так как число всех инвариантных множителей равно меньшему из размеров матрицы, т.е. 6, а число инвариантных множителей, отличных от нуля, равно рангу, т.е. 5, то Е6(l)=0. Далее Е5(l)= (l+1)3(l-1)2 , Е4(l)=(l+1)(l-1)2, Е3 (l)=l+1. Так как мы уже использовали все элементарные делители, то других инвариантных множителей положительной степени не может быть. Поэтому Е2(l)= Е 1(l)=1. 3) Найти каноническую форму диагональной l-матрицы А(l)=íl2, 0, l2+l, l 2-1, l4-l2ý. Здесь порядок матрицы равен 5, а ранг равен 4. Поэтому Е5(l) =0. Элементарные делители диагональных элементов: l2, l, l+1, l+1, l-1, l2, l+1, l-1. Значит, Е4(l)= l2(l+1)( l-1)= l4-l2, Е3(l)= l2(l+1)( l-1)= l4-l2, Е2(l)= l(l+1)= l2+l, Е1(l)=1. Каноническая форма имеет вид í1, l2+l, l4-l2, l4-l2, 0ý. Глава II. Матричные многочлены. §1. Деление матричных многочленов. Обобщенная теорема Безу. Матричным многочленом от переменной l называется выражение вида F(λ) = Ао λm + А1 λm-1 + А2 λm-2 + . + Аm , (1) где Ао, ., Аm - квадратные матрицы одного и того же порядка с элементами из основного поля К. Число m называется степенью многочлена, если Ао≠0. Два многочлена называются равными, если равны матрицы, стоящие в этих многочленах при одинаковых степенях переменной λ. Складываются и перемножаются матричные λ-многочлены по обычным правилам. Ясно, что каждый λ-многочлен можно записать в виде одной матрицы, элементами которой являются обыкновенные многочлены от λ, и обратно. Например, 1 2 + 5 6 λ + 1 0 λ² = λ² +5λ + 1 6l+ 2 0 3 7 -2 0 1 7λ λ²-2λ + 3 . Поэтому матричные λ-многочлены являются лишь особым видом записи λ-матриц . Многочлен F(λ) называется регулярным, если матрица Ао обратима. Сумма (разность) двух матричных многочленов одного и того же порядка может быть представлена в виде многочлена, степень которого не превосходит наибольшей из степеней данных многочленов. Произведение двух матричных многочленов равно многочлену, степень которого меньше или равна сумме степеней сомножителей. Если хотя бы один их двух сомножителей регулярный многочлен, то в этом случае степень произведения всегда равна сумме степеней сомножителей. Пусть даны два матричных многочлена А(λ) и В(λ) одного и того же порядка n, причем В(λ) - регулярный многочлен: А(λ) = Аоλm + А1λm-1 + . + Аm (Ао≠0), В(λ) = Воλр + В1λр-1 + . + Вр (|Во|≠0). Будем говорить, что матричные многочлены Q(λ) и R(λ) являются соответственно правым частным и правым остатком при делении А(λ) на В(λ), если А(λ) = Q(λ)В(λ) + R(λ) (2) и степень R(λ) меньше степени В(λ). Совершенно аналогично будем называть многочлены ^Q(λ) и ^R(λ) соответственно левым частным и левым остатком при делении А(λ) на В(λ), если А(λ) = В(λ) ^Q(λ) + ^R(λ) (3) и степень ^R(λ) меньше степени В(λ). В общем случае многочлены Q(λ) и R(λ) не совпадают с ^Q(λ) и ^R(λ). Покажем, что как правое, так и левое деление матричных многочленов одного и того же порядка всегда выполнимо и однозначно, если делитель - регулярный многочлен. Рассмотрим правое деление А(λ) на В(λ). Если m<p, можно положить Q(λ)=0, R(λ)=А(λ). В случае m≥p для нахождения частного Q(λ) и остатка R(λ) применим обычную схему деления многочлена на многочлен. "Разделим" старший член делимого Аоλm на старший член делителя Воλр . Получим старший член искомого частного АоВо -1λm-p. Умножим этот член справа на делитель В(λ) и полученное произведение вычтем из А(λ). Найдем "первый остаток" А(1)(λ): А(λ)= АоВо -1λm-pВ(λ) + А(1)(λ). (4) Степень m(1) многочлена А(1)(λ) меньше m: А(1)(λ) = Ао(1) λm(1) + . (Ао(1)≠0, m(1)<m). (5) Если m(1)≥p, то повторя этот процесс, получаем: А(1)(λ) = Ао(1)Во -1 λm(1)-р В(λ) + А(2)(λ), (6) А(2)(λ) = А(2)λm(2) + . (m(2)<m(1)) и т.д. Так как степени многочленов А(λ), А(1)(λ), А(2) (λ), . убывают, то на некотором этапе мы придем к остатку R(λ), степень которого меньше р. Тогда из (4), (6) будет следовать: А(λ) = Q(λ) В(λ) + R(λ), где Q(λ) = АоВо-1 λm -р + Ао(1)Во-1 λm (1)-р + . (7) Докажем теперь однозначность правого деления. Пусть одновременно А(λ) = Q(λ) В(λ) + R(λ) (8) и А(λ) = Q*(λ) В(λ) + R*(λ), (9) где степени многочленов R(λ) и R*(λ) меньше степени В(λ), т.е. меньше р. Вычитая почленно (8) из (9) получим: [Q(λ) - Q*(λ)] В(λ) = R*(λ) - R(λ). (10) Если бы Q(λ) - Q*(λ) ≡ 0, то поскольку |Во|≠0, степень левой части равенства (10) равнялось бы сумме степеней В(λ) и Q(λ) - Q*(λ) и потому была бы ≥р. Это невозможно, так как степень многочлена, стоящего в правой части равенства (10), меньше р. Таким образом, Q(λ) - Q*(λ)≡0, а тогда из (10) R*(λ) - R(λ)≡0, т.е. Q(λ) = Q*(λ), R(λ) = R*(λ). Совершенно аналогично устанавливаются существование и единственность левого частного и левого остатка. Т е о р е м а 1. (Обобщенная теорема Безу). При правом (левом) делении матричного многочлена F(λ) на бином λЕ-А остаток от деления равен F(А)(соответственно ^F(A)). Доказательство. Рассмотрим произвольный матричный многочлен n-го порядка F(λ) = Fо λm + F1 λm- 1 + . + Fm (Fо ≠0) (11) Этот многочлен может быть записан и так: F(λ) = λm Fо + λm-1 F1 + . + Fm (12) Обе записи при скалярном λ дают один и тот же результат. Однако если вместо скалярного аргумента λ подставить квадратную матрицу n-го порядка А, то результаты подстановки в (11) и (12) будут различны, так как степени матрицы А могут не быть перестановочными с матричными коэффициентами Fо, F1, ., Fm. Положим F(А) = Fо Аm+ F1 Am-1 + . + Fm (13) и ^F(А) = Аm Fо + Am-1 F1 + . + Fm (14) и будем называть F(А) правым, а ^F(А) - левым значением многочлена F(λ) при подстановке вместо λ матрицы А. Разделим многочлен F(λ) на бином λЕ-А. В данном случае правый остаток R и левый остаток ^R не будут зависеть от λ. Для определения правого остатка рассмотрим обычную схему деления: F(λ) = Fо λm + F1 λm- 1 + . + Fm = Fо λm-1 (λЕ-А) + (Fо А + F1) λm-1 + F2 λm-2 + .= = [Fо λm-1 + (Fо А + F1 ) λm-2] (λЕ-А) + (Fо А2 + F1А1+ F2) λm-2 + F3 λm-3 + . = . . = [Fо λm-1 + (Fо А + F1 ) λm-2 + . + (Fо Аm-1 + F1Аm-2 + . + Fm-1)] (λЕ-А) + + Fо Аm + F1Аm-1 + . + Fm Мы нашли, что R = Fо Аm+ F1Аm-1 + . + Fm = F(А). (15) Совершенно аналогично ^R =^F(А). (16) Из доказанной теоремы следует, что многочлен F(λ) делится без остатка справа (слева) на бином λЕ-А тогда и только тогда, когда F(А)=0 (соответственно ^F(А)=0). Пример Проверить, что А(l)=Q(l)В(l) + R(l).
l3+l 2l3+l2 А(l)= -l3-2l2+1 3l3+l =
= 1 2 l3+ 0 1 l2+ 1 0 l+ 0 0 -1 3 -2 0 0 1 1 0 . 2l2+3 -l2+1 2 -1 l2+ 3 1 В(l)= -l2-1 l2+2 = -1 1 -1 2 , 1 1 3 5 l 2+4 2l2+13 |Bo| = 1, Bo-1 = 1 2 , А0 B0-1 = 2 5 , А0B0 -1В(l) = -l2+1 3l2+12 , l3+l 2l3+l2 l3 +4l 2l3+13l -3l l2 -13l А(1)(l)= -l3-2l2+1 3l3+l - -l3+l 3l3+12l = -2l 2-l+1 -11l ,
0 1 l2+ -3 -13 l+ 0 0 А(1)(l)= -2 0 -1 -11 1 0 ,
Страницы: 1, 2 |
|
© 2010 |
|