Жорданова нормальная форма матриц
Эквивалентность матриц, понятие унимодулярных матриц. Связь подобия числовых матриц с эквивалентность их характеристических матриц. Приведение матрицы к жордановой нормальной форме и особенности минимального многочлена. Решение типовых матричных задач.
Рубрика | Математика |
Вид | дипломная работа |
Язык | русский |
Дата добавления | 20.03.2016 |
Размер файла | 685,0 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Министерство образования и науки Российской Федерации
Федеральное Государственное бюджетное образовательное учреждение
высшего профессионального образования
«Оренбургский государственный педагогический университет»
Физико-математический факультет
Кафедра алгебры, геометрии и истории математики
ЖОРДАНОВА НОРМАЛЬНАЯ ФОРМА МАТРИЦ
(выпускная квалификационная работа)
Специальность - 050201.65 Математика
Квалификация - учитель математики
Специализация - преподавание в классах с углубленным изучением математики
Форма обучения - очная
Допущена к защите: зав. кафедрой алгебры, геометрии и истории математики к.т.н., доцент А.Н. Колобов
Научный руководитель: к.п.н., доцент кафедры алгебры, геометрии и истории математики М.И. Черемисина
Студент Л.В. Муслимова
Оренбург 2015
Содержание
Введение
Глава 1. Жорданова нормальная форма
1.1 Эквивалентность л-матриц
1.2 Унимодулярные матрицы
1.3 Связь подобия числовых матриц с эквивалентность их характеристических матриц
1.4 Жорданова нормальная форма
1.5 Приведение матрицы к жордановой нормальной форме
1.6 Минимальный многочлен
Глава 2. Решение задач
Заключение
Список использованной литературы
Введение
Теория матриц давно вошла в число основных инструментов, используемых другими математическими дисциплинами. Она разрабатывалась из-за необходимости изучения таких систем линейных уравнений, у которых число уравнений не равно числу неизвестных. Матрица определяется как система чисел, расположенных в квадратных или прямоугольных таблицах, состоящих из нескольких строк и столбцов. В середине 19 века возникает матричное исчисление, впервые упомянутое в работах Артура Кэли (1821 -1895) и Джеймса Сильвестра (1814 -- 1897) и занявшее в дальнейшем развитии теории матриц одно из главных мест. К концу 19 века оказались созданными важнейшие разделы матричного исчисления: нормальная форма матрицы, элементарные делители, пары квадратичных форм и т.д.
Жорданова нормальная форма матрицы - одно из фундаментальных понятий линейной алгебры, имеющие большое число приложений в различных разделах математике и физики. Термин «Нормальная (жорданова) форма матрицы» связан с именем К. Жордана. Мари Энмон Камиль (Камилл) Жордан (1838 - 1922) - французский математик, известный благодаря своим фундаментальным работам в теории групп и «Курсу анализа». Он родился в Лионе и учился в Политехнической школе. По образованию Жордан был инженером; позже он преподавал в Политехнической школе и Коллеж де Франс. Основными результатами Жордана были жорданова нормальная форма в линейной алгебре; теорема Жордана о кривой. В математическом анализе мера Жордана используется для построения интеграла Римана (1826 -- 1866).
Теория матриц также широко применяется в различных областях математики, механики, теоретической физики, теоретической электротехнике и т.д.
Аттестационная квалификационная работа состоит из двух глав, введения, заключения. В первой главе рассмотрены многочленные матрицы, или так называемые л-матрицы, элементами которых являются многочлены от буквы л с коэффициентами из основного поля, жорданова нормальная форма. Во второй главе работы представлены решения задач и рассмотрена реальная задача «Население страны», где показано применение жордановой формы матрицы для представления изменения численности населения.
Объектом исследования является теория матриц; предмет исследования - нормальная форма матриц.
Цель: изучить понятие жордановой формы матрицы и ее свойства. Для достижения цели в работе были поставлены задачи:
1) проанализировать научную, методическую, учебную литературу по теме «Жорданова нормальная форма матриц»;
2) изучить л-матрицы и их свойства; нормальную форму матрицы;
3) рассмотреть методы нахождения для заданной матрицы ее каноническую и нормальную жорданову формы;
4) подобрать и решить систему задач по теме «Жорданова нормальная форма матриц».
Глава 1. Жорданова нормальная форма
1.1 Эквивалентность л-матриц
Математическим объектом, записываемым в виде прямоугольной таблице элементов кольца или поля, которая представляет собой совокупность строк и столбцов, на пересечении которых находятся её элементы, называется матрицей. Количество строк и столбцов матрицы задают размер матрицы. Хотя исторически рассматривались, например, треугольные матрицы, в настоящее время говорят исключительно о матрицах прямоугольной формы, так как они являются наиболее удобными.
Определение: Матрица называется квадратной n-го порядка, если число её строк равно числу столбцов и равно одному и тому же числу n
[1, с. 6].
Определение: Квадратной матрице A порядка n можно сопоставить число det A (или |А|, или ?), называемое ее определителем (детерминантом), следующим образом:
1. n = 1.
A=; det A = .
2. n = 2.
A=; det A = = .
3. n = 3.
A= ; det A = = [2, с. 14].
С понятием определителя матрицы связаны понятия вырожденной и невырожденной матрицами.
Определение: Квадратная матрица A называется вырожденной (или особенной), если ее определитель равен нулю, |А| = 0, и невырожденной (или неособенной), если |А| ? 0[3, с. 71].
Определение: Пусть дана квадратная матрица A = (aij)nn. Будем называть минором элемента aij матрицы А определитель (n-1)-го порядка, соответствующий матрице, которая получается из матрицы А вычеркиванием i-ой строки и j-го столбца. Минор элемента aij будем обозначать символом Mij [4, с. 15].
С минором связано еще одно понятие.
Определение: Алгебраическим дополнением элемента aij определителя ?, называется Mij этого элемента, взятого со знаком (-1)i+j [5].
Определение: Матрица А-1 называется обратной по отношению к квадратной матрице A, если при умножении этой матрицы на данную матрицу, как справа, так и слева, получается единичная матрица, то есть
AA-1 = A-1A = E [6, с. 33].
Элементами квадратных матриц порядка n служат многочлены произвольных степеней от одного неизвестного л с коэффициентами из поля P. Такие матрицы называются многочленными, полиноминальными матрицами или л-матрицами. Примером л-матрицы служит характеристическая матрица A - лE произвольной квадратной матрицы A с элементами из поля P ; на главной диагонали этой матрицы стоят многочлены первой степени, вне главной диагонали - многочлены нулевой степени или нули. Всякая матрица с элементами из поля P - такие матрицы для кратности будем называть числовыми - также будет частным случаем л-матрицы: её элементы являются многочленами нулевой степени или нулями [7].
Любую л-матрицу n-го порядка можно представить в виде многочлена с матричными коэффициентами:
A(л) = Am лm + Am-1 лm-1 + … + A1 л + A0,
где Am , Am-1, …, A1, A0 - числовые квадратные матрицы n-го порядка, матрица Am ? 0 - старший коэффициент, матрица A0 - свободный член, неотрицательное целое число m - степень многочлена [8].
Пусть дана матрица
A(л) =
Элементарными преобразованиями этой матрицы называются преобразования следующих четырех типов:
1) Умножение любой строки матрицы A(л) на любое число б из поля Р, отличное от нуля;
2) Умножение любого столбца матрицы A(л) на любое число б из поля Р, отличное от нуля;
3) Прибавление к любой i-й строке матрицы A(л) любой ее j-й строки, i?j, притом умноженной на любой многочлен ц(л) из кольца Р[л];
4) Прибавление к любому i-му столбцу матрицы A(л) любого ее j-го столбца, i ? j, притом умноженного на любой многочлен ц(л) из кольца Р[л] [9, с. 365].
Существует обратное преобразование для каждого из элементарных преобразований л-матрицы, также являющееся элементарным. Так, обратным для преобразования 1) будет элементарное преобразование, состоящее в умножении той же строки на число б-1, существующее ввиду условия б?0; обратным для преобразования 3) будет преобразование, состоящее в прибавлении к i-й строке j-й строки, умноженной на - ц(л).
Пусть, например, нужно переставить i-ю и j-ю строки матрицы A(л).Это можно сделать при помощи четырех элементарных преобразований, как показывает следующая схема:
>>>>
Здесь последовательно выполнялись такие преобразования: а) к i-й строке прибавлялась j-я; б) из j-й строки вычиталась новая i-я; в) к новой i-й строке прибавлялась новая j-я; г) новая j-я строка умножалась на -1.
Определение: Две л-матрицы называются эквивалентными, если одна может быть получена из другой некоторой последовательностью элементарных преобразований [10, с. 241].
Говорят, что л-матрицы А(л) и В(л) эквивалентны и записывать это символом А(л) ~ В(л).
Из определения эквивалентности л-матриц непосредственно вытекают следующие свойства:
1. Отношение эквивалентности транзитивно: если A(л) ~ B(л), а B(л) ~ G(л), то A(л) ~ G(л).
В самом деле, по условию B(л) можно получить из G(л), а A(л) из B(л) цепочкой элементарных преобразований; следовательно, A(л) можно получить цепочкой элементарных преобразований из G(л).
2. Отношение эквивалентности симметрично: если A(л) ~ B(л), то B(л) ~ A(л). Другими словами A(л) можно получить из B(л) цепочкой элементарных преобразований, то и B(л) можно получить из A(л) цепочкой элементарных преобразований.
Докажем сначала, что если A(л) можно получить из B(л) каким-нибудь одним элементарным преобразованием, то B(л) также можно получить из A(л) одним элементарным преобразованием. Для этого мы пересмотрим все четыре типа элементарных преобразований. Пусть A(л) получается из B(л) преобразований типа 1, т. е. умножением какой-либо i-й строки B(л) на число б?0. Тогда, умножая i-ю строку A(л) на б-1, получим B(л). Пусть теперь A(л) получится из B(л) преобразованием типа 2, например, путем прибавления к i-й строки матрицы B(л) ее j-й строки, умноженной на f(л). В таком случае, прибавляя к i-й строке матрицы A(л) ее j-ю строку, умноженную на - f(л), получим обратно B(л). То же самое можно сказать и относительно преобразований типа 3, 4. Для каждого элементарного преобразования существует элементарное преобразование, ему обратное, отменяющее результат первого. Поэтому, если матрица A(л) получается из B(л) цепочкой элементарных преобразований, то, совершая обратные преобразования в обратном порядке, из матрицы A(л) получим матрицу B(л), что и требовалось.
3. Отношение эквивалентности рефлексивно: каждая матрица эквивалентна самой себе [11, с. 176].
Например, совершая над A(л) два взаимно обратных преобразования, получим снова A(л).
Все квадратные л-матрицы порядка n над полем Р распадаются на непересекающиеся классы эквивалентных матриц.
Всякий класс л-матриц будет содержать единственную каноническую л -матрицу: каждая л-матрица эквивалентна единственной канонической матрице. Такая матрица называется нормальной формой или канонической формой данной матрицы.
Многочлены, которые стоят на главной диагонали канонической формы данной л-матрицы, называются инвариантными множителями данной матрицы. Одним из методов при вычислении инвариантных множителей является приведение данной л-матрицы к канонической форме.
Определение: Канонической л-матрицей называется л-матрица, обладающая следующими тремя свойствами:
а) Эта матрица диагональная, т.е. имеет вид
; (1)
б) Всякий многочлен ei(л), i=2, 3, …, n, нацело делится на многочлен ei-1(л), т.е. всякий следующий многочлен делится нацело на предыдущий;
в) Старший коэффициент каждого многочлена ei(л), i=1, 2, …, n, равен единице, если этот многочлен отличен от нуля [9, с. 365].
Если среди многочленов ei(л), которые стоят на главной диагонали канонической л-матрицы (1), встречаются равные нулю, то непременно занимают на главной диагонали последние места, ввиду свойства б). С другой стороны, если среди многочленов ei(л) встречаются многочлены нулевой степени, то, по свойству в), все они равны 1 и, по свойству б), занимают на главной диагонали матрицы (1) первые места.
К числу канонических л-матриц принадлежат некоторые числовые матрицы, в том числе матрицы единичная и нулевая.
Определение: Квадратная матрица n-го порядка, у которой все элементы главной диагонали равны единице, а остальные равны нулю, называется единичной матрицей n-го порядка, и она обозначается буквой Е [1, с. 6]. Единичная матрица имеет вид:
E = [1, с. 8].
Определение: Матрица любого размера называется нулевой, или нуль-матрицей, если все её элементы равны нулю [1, с. 8].
Теорема: Всякая л-матрица эквивалентна некоторой канонической л-матрице, т. е., иными словами, она приводится элементарными преобразованиями к каноническому виду [9, с. 366].
Доказательство: (метод математической индукции по порядку n рассматриваемых л-матриц).
При n=1 будет
A(л)=(a(л)).
Если a(л)=0, то матрица уже каноническая.
Если же a(л)? 0, то достаточно разделить многочлен a(л) на его старший коэффициент - это будет элементарное преобразование матрицы - и получим каноническую матрицу.
Пусть теорема уже доказана для л-матриц порядка n - 1. Рассмотрим произвольную л-матриц A(л) порядка n. Если она нулевая, то уже является канонической и доказывать нечего. Будем считать поэтому, что среди элементов матрицы A(л) имеются ненулевые.
Переставляя, если понадобится, строки и столбцы матрицы A(л), можно перевести один из ненулевых элементов в левый верхний угол. Т.о., среди л-матриц, эквивалентных матрице A(л), имеются такие, в верхнем углу которых стоит ненулевой многочлен. Рассмотрим все такие матрицы. Многочлены, стоящие в левом верхнем углу этих матриц, могут иметь разные степени. Степень многочлена является натуральным числом, а во всяком непустом множестве натуральных чисел существует наименьшее число. Можно найти, следовательно, среди всех л-матриц, эквивалентных матрице A(л) и имеющих ненулевой элемент в левом верхнем углу, одну из таких, что многочлен, стоящий в её левом верхнем углу, имеет наименьшую возможную степень. Деля первую строку этой матрицы на старший коэффициент указанного многочлена, получим такую л-матрицу, эквивалентную матрице A(л),
A(л) ~ ,
что e1(л)?0, старший коэффициент этого многочлена равен 1 и никакой комбинацией элементарных преобразований нельзя перейти от полученной матрицы к такой матрице, в верхнем левом углу которой стоял бы ненулевой многочлен меньшей степени.
Докажем, что все элементы первой строки и первого столбца полученной матрицы нацело делятся на e1(л). Пусть, например, для
2 ? j ? n
b1j(л) = e1(л) q(л) + r(л),
где степень r(л) меньше степени e1(л), если r(л) ? 0. Тогда, вычитая из j-го столбца матрицы её первый столбец, умноженный на q(л), а затем переставляя первый и j-й столбцы, приходим к такой матрице, эквивалентной матрице A(л), в левом верхнем углу которой стоит многочлен r(л), т.е. многочлен меньшей степени, чем e1(л), что противоречит выбору этого многочлена. Отсюда следует r(л) = 0, что и требовалось доказать.
Вычитая теперь из j-го столбца нашей матрице её первый столбец, умноженный на q(л), заменим элемент b1j(л) нулем. Делая такие преобразования для j = 2, 3, …, n, заменим нулями все элементы b1j(л). Аналогичным путем заменяются нулями и все элементы bi1(л), i = 2, 3, ..., n. Придем к такой матрице, эквивалентной матрице A(л), в левом верхнем углу которой стоит многочлен e1(л), а все остальные элементы первой строки и первого столбца равны нулю,
A(л) ~ . (2)
По индуктивному предположению, матрица (n - 1)-го порядка, стоящая в правом нижнем углу полученной нами матрицы (2), элементарными преобразованиями приводится к каноническому виду:
~ .
Совершив эти же преобразования над соответствующими строками и столбцами матрицы (2) - при этом первая строка и первый столбец этой матрице останутся, очевидно, без изменения, - получим, что
A(л)~. (3)
Для доказательства того, что матрица (3) является канонической, остается показать, что e2(л) нацело делится на e1(л). Пусть
e2(л) = e1(л) q(л) + r(л),
где r(л) ? 0 и степень r(л) меньше степени e1(л). Прибавляя, однако, ко второму столбцу матрицы (3) ее первый столбец, умноженный на q(л), а затем вычитая из второй строки первую строку, заменим элемент e2(л) элементом r(л). Переставляя, далее, первые две строки и первые два столбца, переместим многочлен r(л) в левый верхний угол матрицы, что противоречит выбору многочлена e1(л).
Теорема о приведении л-матрицы к каноническому виду доказана. Эта теорема должна быть дополнена следующей теоремой единственности:
Всякая л-матрица эквивалентна лишь одной канонической матрице [9, с. 368].
Доказательство: В самом деле, пусть дана произвольная л-матрица A(л) порядка n. Фиксируем некоторое натуральное число k, 1 ? k ?n, и рассмотрим все миноры k-го порядка матрицы A(л). Вычисляя эти миноры, получим конечную систему многочленов от л; наибольший общий делитель этой системы многочленов, взятый со старшим коэффициентом 1, обозначим через dk(л).
Имеем многочлены
d1(л), d2(л),…, dn(л), (4)
однозначно определяемые самой матрицей A(л). При этом d1(л) есть наибольший общий делитель всех элементов матрицы A(л), взятый с коэффициентом 1, а dn(л) равен определителю матрицы A(л), деленному на его старший коэффициент. Заметим также, что если матрица A(л) имеет ранг r, то
dr+1(л) =…= dn(л) =0,
В то время как все остальные многочлены системы (4) отличны от нуля.
Теорема: Наибольший общий делитель dk(л) всех миноров k-го порядка л-матрицы A(л), k = 1, 2, …, n, не меняется при выполнении в матрице A(л) элементарных преобразований [9, с. 369].
Доказательство: Если в матрице выполняются элементарные преобразования типа 1) и 2). Например, если i-я строка матрицы умножается на число б из поля P, то миноры k-ого порядка, через которые i-я строка проходит, будут умножаться на б, все же остальные миноры k-ого порядка останутся без изменения. Но при разыскании наибольшего общего делителя нескольких многочленов любые из этих многочленов можно беспрепятственно умножать на отличные от нуля числа из поля P. Рассмотрим теперь элементарные преобразования типа 3) или 4). Пусть, например, к i-й строке матрицы A(л) прибавляется ее j-я строка, i ? j, умноженная на многочлен ц(л); получающуюся после этого преобразования матрицу обозначим за ?(л), а наибольший общий делитель всех ее миноров k-го порядка, взятый со старшим коэффициентом 1, - через ?k(л).
При указанном преобразовании не будут меняться те миноры, через которые i-я строка не проходит. Не меняются и те миноры, через которые проходят как i-я, так и j-й строки, т. к. определитель не меняется от прибавления к одной его строке кратного другой его строки. Возьмем любой из тех миноров k-го порядка, через которые проходит i-я строка, но не проходит j-я; обозначим его М. Соответствующий минор матрицы ?(л) можно представить как сумму минора М и умноженного на ц(л) минора Мм матрицы A(л), получающегося из минора заменой элементов i-й строки матрицы A(л) соответствующими элементами ее j-й строки. Т. к. и М, и Мм делятся на dk(л), то и М + ц(л) Мм будет делиться на dk(л).
Отсюда следует, что все миноры k-го порядка матрицы ?(л) нацело делятся на dk(л), а поэтому ?k(л) делится на dk(л). Т. к. для рассматриваемого элементарного преобразования того же типа, что и dk(л) делится на ?k(л). Если же учесть, что старшие коэффициенты обоих этих многочленов равны 1, то ?k(л) = dk(л), что и требовалось доказать.
Всем л-матрицам, эквивалентным матрице A(), соответствует один и тот же набор многочленов (4). Это относится к любой канонической матрице, эквивалентной A(л). Пусть (3) будет одна из таких матриц.
Пользуясь матрицей (3) вычислим многочлен dk(л), k = 1, 2, …, n. Ясно, что минор k-го порядка, стоящий в левом верхнем углу этой матрицей равен произведению
e1(л), e2(л), …, ek(л). (5)
Если, далее, берем в матрице (3) минор k-го порядка, стоящий в строках с номерами i1, i2, …, ik, где i1 < i2 < …, < ik, и в столбцах с теми же самыми номерами, то этот минор равен произведению e1(л) ei2(л) … eik(л), которое делится на (5). Действительно, 1 i1 и поэтому ei1(л) делится на e1(л), 2 i2, и поэтому ei2(л) делится на e2(л) и т. д. Наконец, если в матрице (3) взят минор k-го порядка, через который хотя бы для одного i проходит i-я строка этой матрицы, но не проходит её i-й столбец, то этот минор содержит нулевую строку и поэтому равен нулю.
Произведение (5) и будет наибольшим общим делителем всех миноров k-го порядка матрицы (3), и поэтому и исходной матрицы A(л),
dk(л) = e1(л) e2(л) … ek(л), k = 1, 2, …, n. (6)
Многочлены ek(л), k = 1, 2, …, n, определятся однозначным образом самой матрицей A(л). Пусть ранг этой матрице равен r. Тогда, как мы знаем,
dr(л) ?0,
но
dr+1(л) = 0,
а поэтому, ввиду (6),
er+1(л) = 0.
Отсюда, ввиду свойств канонической матрицы, вообще следует, что если ранг r матрицы A(л) меньше n, то
er+1(л) = er+2(л) = … = en(л) = 0. (7)
с другой стороны, для k r из (6) следует, ввиду
dk-1(л) ? 0,
что
. (8)
Этим заканчивается доказательство единственности канонического вида л-матрицы. Одновременно получили способ разыскания многочленов ek(л), называемых инвариантными множителями A(л).
1.2 Унимодулярные матрицы
Критерию эквивалентности л-матрицы можно придать формулировку:
Две л-матрицы тогда и только тогда эквивалентны,
1) если они приводятся к одному и тому же каноническому виду;
2) если они обладают одинаковыми инвариантными множителями.
Выделим еще один критерий, имеющий уже иной характер.
К числу канонических л-матриц принадлежит единичная матрица E. Назовем л-матрицу U(л) унимодулярной, если она имеет матрицу E своим каноническим видом, т.е. если все её инвариантные множители равны единице.
Теорема: л-матрица U(л) тогда и только тогда унимодулярна, если ее определитель отличен от нуля, но не зависти от л , т. е. является отличным от нуля числом из основного поля P [9, с. 371].
Доказательство: Действительно, если U(л) ~ E, то этим двум матрицам соответствует один и тот же многочлен dn(л). Но для единичной матрицы dn(л) = 1. Отсюда следует, что определитель матрицы U(л), отличающийся от d лишь отличным от нуля числовым множителем, будет отличным от нуля числом из поля Р. Обратно, если определитель матрицы U(л) отличен от нуля и не зависит от л, то для этой матрицы многочлен dn(л) будет равен 1, а поэтому по формуле
dk(л) = e1(л) e2(л) … e1k(л),
все инвариантные множители ei(л) матрицы U(л), i= 1, 2, …, n, равны единицы.
А это означает, что всякая невырожденная числовая матрица является унимодулярной л-матрицей. Унимодулярная л-матрица может иметь очень сложный вид. Так, л-матрица
унимодулярна, т. к. ее определитель равен 20, т. е. отличен от нуля и от л не зависит.
Произведение унимодулярных л-матриц само унимодулярно - достаточно вспомнить, что при умножении матриц их определители перемножаются.
Теорема: л-матрица U(л) тогда и только тогда унимодулярна, если для нее существует обратная матрица, также являющаяся л-матрицей [9, с. 372].
Доказательство: Пусть дана невырожденная л-матрица, то для разыскания обратной матрицы обычным способом, разделим алгебраическое дополнение к элементам данной матрицы на определитель этой матрицы, т.е. на некоторый многочлен от л. Поэтому в общем случае элементы обратной матрицы будут рациональными дробями от л, а не многочленами от л, т.е. эта матрица не будут л-матрицей. Если же дана унимодулярная матрица, то делить алгебраическое дополнение придется лишь отличное от нуля число из поля P, т. е. элементы обратной матрицы будут многочленами от л и поэтому обратная матрица сама будет л-матрицей.
Если л-матрица U(л) обладает обратной л-матрицей U-1(л), то определители этих обеих матриц являются многочленами от л, их произведение равно 1, а поэтому оба определителя должны быть многочленами нулевой степени.
Из последнего замечания вытекает добавление к доказанной выше теореме:
л-матрица, обратная к унимодулярной л-матрице, сама унимодулярна.
Рассмотрим новый критерий эквивалентности л-матриц:
Теорема: Две л-матрицы A(л) и B(л) порядка n тогда и только тогда эквивалентны, если существует такие унимодулярные л-матрицы U(л) и V(л) того же порядка n, что
B(л)=U(л) A(л) V(л). (1)
Введем понятие, используемое при доказательстве этого критерия.
Элементарная матрица - это числовая (и, следовательно, л-) матрица, отличающиеся от единичной матрицы лишь тем, что на некотором i-м месте главной диагонали, 1? i ?n, стоит произвольное число б из поля P, отличное от нуля.
(2)
С другой стороны, элементарной матрицей называется л-матрица, которая отличается от единичной матрицы лишь тем, что на пересечении i-й строки и j-го столбца, 1 ? i ? n, 1 ? j ? n, причем i ? j, стоит произвольный многочлен ц(л) из кольца P[л].
(3)
Теорема: Всякая элементарная матрица унимодулярна.
Доказательство: В самом деле, определитель матрицы (2) равен б, но, по условию б ? 0; определитель матрицы (3) равен 1.
Выполнение в л-матрице A(л) любого элементарного преобразования равносильно умножению этой матрицы слева и справа на некоторую элементарную матрицу.
Проверим справедливость следующих четырех утверждений:
1. умножение матрицы A(л) слева на матрицу (2) равносильно умножению i-й строки матрицы A(л) на число б;
2. умножение матрицы A(л) справа на матрицу (2) равносильно умножению i-го столбца матрицы A(л) на число б;
3. умножение матрицы A(л) слева на матрицу (3) равносильно прибавлению к i-й строке матрицы A(л) ее j-й строки, умноженной на ц(л);
4. умножение матрицы A(л) справа на матрицу (3) равносильно прибавлению к j-му столбцу матрицы A(л) ее i-го столбца, умноженного на ц(л);
Доказательство критерия эквивалентности л-матриц.
Если A(л) ~ B(л), то от A(л)можно перейти к B(л) при помощи конечного числа элементарных преобразований. Заменяя каждое из этих преобразований умножением слева или справа на элементарную матрицу, приходим к равенству
B(л) = U1(л) … Uk(л)A(л)V1(л) … Vl(л), (4)
где все матрицы U1(л), … Uk(л), V1(л) … Vl(л) элементарны и, следовательно, унимодулярны. Унимодулярными будут, поэтому и матрицы
U(л) = U1(л) … Uk(л), V(л) = V1(л) … Vl(л), (5)
являющееся произведениями унимодулярных матриц, а равенство (4) перепишется в виде (1). Замечаем, что если k = 0, т. е. элементарные преобразования совершались лишь над столбцами, то полагаем просто U(л) = E.
Из выше приведенных доказательств следует, что л-матрица тогда и только тогда унимодулярна, если она представима в виде произведения элементарных матриц.
Т.о. любая произвольная унимодулярная матрица W(л) будет эквивалентна единичной матрице E. применяя проведенное выше доказательство вместо матриц A(л) и B(л) к матрицам E и W(л), из (4) получим равенство
W(л) = U1(л) … Uk(л)V1(л) … Vl(л),
т.е. матрица W(л) оказалась представленной в виде произведения элементарных матриц.
Доказательство обратного утверждения критерия.
Пусть для матриц A(л) и B(л) существуют такие унимодулярные матрицы U(л) и V(л), что имеет место равенство (1). По доказанному, матрицы U(л) и V(л) можно представить в виде произведений элементарных матриц; пусть это будет представление (5). Равенство (1) перепишется теперь в виде (4) и, заменяя каждое умножение элементарной матрицы соответствующим элементарным преобразованием, получим, что A(л) ~ B(л).
Матричным л-многочленом порядка n над полем P многочлен от л, коэффициенты которого служат квадратные матрицы одного и того же порядка n с элементами из поля P, его общим видом будет
. (6)
Всякий матричный л-многочлен порядка n можно записать в виде л-матрицы порядка n. Так,
.
И обратно, всякая л-матрица порядка n может быть записана в виде матричного л-многочлена порядка n. Так,
.
Соответствие между л-матрицами и матричными л-многочленами является взаимно однозначным. В самом деле, равенство матричного л-многочленов вида (6) как матриц равносильно равенству матричных коэффициентов при одинаковых степенях л, а умножение матрицы на л равносильно умножению ее на скалярную матрицу с л на главной диагонали.
Пусть дана л-матрица A(л), причем
A(л) =,
Где матрица A0 не является нулевой. Число k называется степенью - матрицы A(л).
Теорема: Пусть над полем P даны л-матрицы порядка n
A(л) =,
B(л) =,
причем предположим, что матрица B0 невырожденная, т. е. существует матрица B0-1. Тогда над полем P можно найти такие л-матрицы Q1(л) и R1(л) порядка n, что
A(л) =
причем степень R1(л) меньше степени B(л) или же R1(л) = 0. С другой стороны, над полем P можно найти такие л-матрицы Q2(л) и R2(л) порядка n, что
A(л) =.
Причем степень R2(л) меньше степени B(л) или же R2(л) = 0. Матрицы Q1(л) и R1(л), а также Q2(л) и R2(л), удовлетворяющие этим условиям, определяются однозначно.
1.3 Связь подобия числовых матриц с эквивалентность их характеристических матриц
Начнем с рассмотрения вопроса о подобии числовых матриц A и B, т. е. матрицы с элементами из основного поля P. С другой стороны, их характеристические матрицы A - лE и B - лE являются л-матрицами и вопрос об эквивалентности этих матриц решается вполне эффективно.
Определение: Матрица A называется подобной матрице B, если существует такая невырожденная матрица C, что
A = C-1BC, (1)
где C-1 является обратной или обращением матрицы C, что выполняется следующие
CC-1 = C-1C = E.
В этом случае говорят, что A получается преобразованием матрицы B при помощи C. Умножая равенство (1) слева на C и справа на C-1, получим
B = CAC-1 = (C-1)-1AC-1.
Докажем следующую теорему.
Теорема: Матрицы A и B с элементами из поля P тогда и только тогда подобны, если их характеристические матрицы A - лE и B - лE эквивалентны.
Доказательство: Допустим, что матрицы A и B подобны, т. е. над полем P существует такая невырожденная матрица С, что
B = C-1AC.
Тогда
C-1(A - лE) C = C-1AC - л(C-1EC) = B - лE.
Невырожденные числовые матрицы C-1 и C являются унимодулярными л-матрицами. Заметим, что матрица B - лE получена умножением матрицы A - лE слева и справа на унимодулярные матрицы, т. е. A - лE ~ B - лE.
Доказательство обратного утверждения является наиболее сложным.
Пусть
A - лE ~ B - лE.
Тогда существуют такие унимодулярные матрицы U(л) и V(л), что
U(л)( A - лE) V(л) = B - лE. (2)
Учитывая, что для унимодулярных матриц обратные матрицы существуют и являются л-матрицами, выведем следующие равенства, используемые ниже:
(3)
Т. к. л-матрица B - лE имеет по л степень 1, причем старший коэффициентом соответствующего матричного многочлена служит невырожденная матрица - E, то к матрицам U(л) и B - лE можно применить алгоритм деления с остатком: существуют такие матрицы Q1(л) и R1 - последняя, если она отлична от нуля, должна иметь по л степень 0, т. е. от л не зависит, что
U(л) = (B - лE) Q1(л) + R1. (4)
Аналогично
V(л) = Q2(л) (B - лE) + R2 (5)
Используя (4) и (5), из (3) получаем:
R1 (A - лE) R2 = (B - лE) - U(л) (A - лE) Q2(л) (B - лE) - (B - лE) Q1(л) (A - лE) V(л) + (B - лE) Q1(л) (A - лE) Q2(л) (B - лE)
Или, ввиду (3),
R1 (A - лE) R2 = (B - лE) - (B - лE) V-1(л) Q2(л) (B - лE) - (B - лE) Q1(л) U-1(л) (B - лE) + (B - лE) Q1(л) (A - лE) Q2(л) (B - лE) = (B - лE) {E - [V-1(л) Q2(л) + Q1(л) U-1(л) - Q1(л) (A - лE)Q2(л)](B - лE)}.
Квадратная скобка, стоящая справа, равна в действительности нулю: в противном случае она, являясь л-матрицей, т. к. и V-1(л) и U-1(л) суть л-матрицы, имела бы по меньшей мере степень 0, а тогда степень фигурной скобки была бы не меньше 1 и, следовательно, степень всей правой части была бы не меньше 2. Это, однако, невозможно, т. к. слева стоит - матрица степени 1.
Т.о.,
R1 (A - лE) R2 = B - лE,
откуда, приравнивая матричные коэффициенты при одинаковых степенях л, получаем
R1 A R2 = B, (6)
R1 R2 = E. (7)
Равенство (7) показывает, что числовая матрица R2 не только отлична от нуля, но даже является невырожденной, причем
R2-1 = R1,
а тогда равенство (6) принимает вид
R2-1 A R2 = B,
что доказывает подобие матриц A и B.
Одновременно научились находить ту невырожденную матрицу R2, которая трансформирует матрицу A в матрицу B. Если матрицы A - лE и B - лE эквивалентны, то первая конечным числом элементарных преобразований переводится во вторую. Берем те из этих преобразований, которые относятся к столбцам, и произведение соответствующих элементарных матриц, взятых в том же порядке, обозначаем через V(л). Делим затем V(л) на B - лE, причем так, чтобы частное стояло слева от делителя. Остаток от этого деления и будет матрицей R2.
Указанное деление может на самом деле не выполниться.
Лемма: Пусть
V(л) = V0 лs + V1 лs-1 + … Vs-1 л + Vs, V0 ? 0. (8)
Если
V(л) = (лE - B) Q1(л) + R1, (9)
V(л) = Q2(л) (лE - B) + R2,
то
R1 = BsV0 + Bs-1V1 + … BVs-1 +Vs, (10)
R2 = V0Bs + V1Bs-1 + … Vs-1B + Vs.
Доказательство: Достаточно доказать хотя бы первое из двух утверждений леммы - второе доказывается вполне аналогично. Доказательство состоит в непосредственной проверки справедливости равенства (9), если многочлен V(л) будет заменен его записью (8), вместо R1 будет подставлено (10), а в качестве Q1(л) будет взят многочлен
Q1(л) = V0лs-1 + (BV0 + V1) лs-2 + (B2V0 + BV1 + V2) лs-3 + …+ (Bs-1V0 + Bs-2V1 + … Vs-1).
1.4 Жорданова нормальная форма
В этом параграфе будем рассматривать квадратные матрицы порядка n с элементами из поля P. Выделим один специальный тип таких матриц, так называемые жордановы матрицы, и будет показано, что эти классы служат нормальной формой для весьма широкого класса матриц.
Определение: Характеристической матрицей для матрицы A называется матрица A - лE. Определитель характеристической матрицы представляет собой скалярный многочлен относительно л и называется характеристическим многочленом матрицы A [12, с. 92].
Корни характеристического многочлена матрицы A называются ее характеристическими числами или собственными значениями.
Определение: Матрицы, все характеристические корни которых лежат в основном поле P (и только такие матрицы), подобны некоторым жордановым матрицам, т. е., как говорят, они приводятся к жордановой нормальной форме.
Отсюда будет следовать, если в качестве поля P взято поле комплексных чисел, что всякая матрица с комплексными элементами приводится в поле комплексных чисел к жордановой нормальной форме.
Определение: Жордановой клеткой порядка n, соответствующей собственному значению л0, называется nn-матрица
Jn(л). (1)
называется жордановой клеткой порядка k, относящейся к числу л0.
Иными словами, на ее главной диагонали стоит одно и то же число л0 из поля P; ближайшая к главной диагонали сверху, сплошь занята числом 1; все остальные элементы матрицы равны нулю. Так,
будут соответственно жордановыми клетками первого, второго и третьего порядков.
Определение: Жордановой матрицей порядка n называется матрица порядка n, которая имеет вид
J = (2)
Здесь на главной диагонали будут находиться жордановы клетки J1, J2, …, Js некоторых, не обязательно различных, порядков, относящиеся к некоторым числам, не обязательно различных, из поля Р, все места вне этих клеток заняты нулями.
Можно описать строение жордановой матрицы, не используя понятие жордановой клетки. Очевидно, что матрица J будет жордановой матрицей тогда и только тогда, если она имеет вид
,
где лi, i = 1, 2, …, n, - произвольные числа из поля P, а каждое еj, j = 1, 2, …, n-1, равно единице или нулю, причем, если еj = 1, то
лj = лj+1.
Частным случаем жордановых матриц являются диагональные матрицы, т. е. это будут в точности те жордановы матрицы, у которых жордановы клетки имеет порядок 1.
Найдем канонический вид для характеристической матрицы J - лE произвольной жордановой матрицы J порядка n.
Для начала найдем канонический вид для характеристической матрицы
. (3)
Одной жордановой клетки (1) порядка k. Вычислим определитель этой матрицы и напомним, что старший коэффициента многочлена dk(л) должен равняться 1, отсюда следует, что
dk (л) = (л - л0)k.
С другой стороны, среди миноров (k - 1)-го порядка матрицы (3) имеется минор, равный единице, а именно тот, который получается вычеркиванием первого столбца и последней строки этой матрицы. Поэтому
dk-1 (л) = 1.
А это означает, что каноническим видом для матрицы (3) служит следующая л-матрица порядка k:
. (4)
Лемма: Если многочлены ц1(л), ц2(л), …, цt(л) из кольца P[л] попарно взаимно просты, то имеет место следующая эквивалентность:
.
Доказательство: Достаточно, рассмотреть случай, когда t = 2. Т .к. многочлены ц1 (л) и ц2 (л) взаимно просты, то в кольце P[л] существуют такие многочлены u1 (л) и u2 (л), что
ц1 (л) u1 (л) + ц2 (л) u2 (л) = 1.
Поэтому
~ ~ = ~ ~ ~ ~ ,
что и требовалось доказать.
Прейдем к рассмотрению характеристической матрицы
J - лE = . (5)
Для жордановой матрицы J вида (2); здесь Ei, i = 1, 2, …, s, есть единичная матрица того же порядка, что и клетка Ji. Пусть к жордановым клеткам матрицы J будут относиться следующие различные числа: л1, л2, …, лt, где t ? s. Пусть, далее, к числу лi, i = 1, 2, …, t, относятся qi жордановых клеток, qi ? 1, и пусть порядки этих клеток, расположенные в невозрастающем порядке, будут
ki1? ki2 ? … ? ki qi, (6)
применим элементарные преобразования к тем строкам и столбцам матрицы (5), которые проходят через клетку Ji - лEi этой матрицы, при этом очевидно, не будем затрагивать других диагональных клеток. Отсюда следует, что в матрице (5) можно при помощи элементарных преобразований заменить каждую клетку Ji - лEi, i = 1, 2, …, s, соответствующей клеткой вида (4). Таким образом, матрица J - лE эквивалентна диагональной матрице, на диагонали которой стоят, помимо некоторого числа единиц, также следующие многочлены, соответствующие всем жордановым клеткам матрицы J
(7)
При этом не указываем те места на диагонали, на которых многочлены (7), т. к. в любой диагональной л-матрице диагональные элементы можно произвольно переставлять при помощи перестановок строк и одноименных столбцов.
Пусть q - наибольшее среди чисел qi, i = 1, 2, …, t. Обозначим через en-i+1 (л) произведение многочленов, стоящих в j-м столбце таблицы (7), j = 1, 2, …, q, т.е.
en-i+1 (л) = ; (8)
если при этом в j-м столбце имеются пустые места - для некоторых i может оказаться, что qi < j, - то соответствующие множители в (8) будем считать равными единице. Т.к. числа л1, л2, …, лt по условию различные, то степени линейных двучленов, которые стоят в j-м столбце таблицы (7), попарно взаимно просты. Поэтому, на основании доказанной выше леммы, они могут быть заменены при помощи элементарных преобразований в рассматриваемой матрице их произведением en-i+1 (л) и некоторым числом единиц.
Проделав это для j = 1, 2, …, q, получим, что
J - лE ~ . (9)
Это и есть искомый канонический вид матрицы J - лE. В самом деле, старшие коэффициенты, стоящих в (9) на главной диагонали, всех многочленов, равны единице и каждый из этих многочленов будет нацело делится на предыдущий ввиду условия (6).
Т.о. сразу можно писать канонический вид ее характеристической матрицы по виду данной жордановой матрицы J.
Докажем следующую теорему:
Две жордановы матрицы тогда и только тогда подобны, если они состоят из одних и тех же жордановых клеток, т. е. отличаются, может быть, лишь расположением этих клеток вдоль главной диагонали [9, с. 384].
Доказательство: Таблица многочленов (7) полностью определяется набором жордановых клеток жордановой матрицы J и в ней никак не отражалось расположение жордановых клеток вдоль главной диагонали этой матрицы. Отсюда следует, что если жордановы матрицы J1 и J2 обладают одним и тем же набором жордановых клеток, то им соответствует одна и та же таблица многочленов (7), а поэтому одни и те же многочлены (8). Т.о., характеристические матрицы J1 - лE и J2 - лE обладают одинаковыми инвариантными множителями, т.е. эквивалентны, а поэтому сами матрицы J1 и J2 подобны.
Обратно, если жордановы матрицы J1 и J2 подобны, то их характеристические матрицы будут обладать одинаковыми инвариантными множителями. Пусть многочлены (8) для j = 1, 2, …, q будут те из этих инвариантных множителей, которые отличны от единицы. Но по многочленам (8) восстанавливается таблица многочленов (7). Многочлены (8) разлагаются в произведение степеней линейных множителей, т.к. этим свойством обладают, как уже доказано, инвариантные множители характеристической матрицы для любой жордановой матрицы. Таблица (7) как раз и состоит из всех тех максимальных степеней линейных множителей, на которые разлагаются многочлены (8). Наконец, по таблице (7) восстанавливаются жордановы клетки исходных жордановых матриц: каждому многочлену (л - л)kij из таблицы (7) соответствует жорданова клетка порядка kij, относящаяся к числу лi. Таким образом, доказано, что матрица J1 и J2 состоят из одних и тех же жордановых клеток и отличаются, быть может, лишь их расположением.
Жорданова матрица сама будет диагональной, если она подобна диагональной матрице. Отсюда следует, что две диагональные матрицы тогда и только тогда подобны, если получаются друг от друга перестановкой чисел, стоящих на главной диагонали.
1.5 Приведение матрицы к жордановой нормальной форме
Если матрицу A с элементами из поля P можно привести к жордановой нормальной форме, т. е. она будет подобна жордановой матрице. То жорданова нормальная форма определяется для матрицы A однозначно с точностью до расположения жордановых клеток на главной диагонали, как следует из доказанной выше теоремы. Условие, допускающие такое приведение матрицы A, указывается в следующей теореме, чье доказательство дает практический способ для нахождения жордановой матрицы, подобной матрицы A, если такая существует. При этом заметим, что приводимость в поле P означает, что все элементы трансформирующей матрицы содержатся в поле P.
Теорема: Матрица A с элементами из поля P тогда и только тогда приводится в поле P к жордановой нормальной форме, если все характеристические корни матрицы A лежат в самом основном поле P [9, с. 385].
Доказательство: Если матрица A подобна жордановой матрице J, то у этих матриц будут одними и теми же характеристическими корнями. Характеристические корни матрицы J находятся без всяких затруднения, т.к. определитель матрицы J - лE равен произведению ее элементов, стоящих на главной диагонали, то многочлен | J - лE | разлагается над полем P на линейные множители и его корнями служат числа, стоящие на главной диагонали матрицы J, и только они.
Пусть все характеристические корни матрицы A лежат в самом поле P. Если отличные от 1 инвариантные множители матрицы J - лE будут
en-q+1(л), …, en-1(л), en(л), (1)
то
| J - лE | = (-1)n en-q+1(л) … en-1(л) en(л).
Действительно, определители матрицы J - лE и её канонической матрицы могут отличаться друг от друга лишь постоянным множителем, который на самом деле равен (-1)n, т. к. именно таков старший коэффициент характеристического многочлена | J - лE |. Таким образом, среди многочленов (1) нет равных нулю, сумма степеней этих многочленов равна n и все они разлагаются над полем P на линейные множители, т.к. по условию, многочлен | J - лE | обладает таким разложением.
Пусть
en-j+1 (л) =; (2)
будут разложения многочленов (1) в произведения степеней линейных множителей. Назовем элементарными делителями многочлена en-j+1, j = 1, 2,…, q, отличные от единицы степени различных линейных двучленов, входящих в его разложение (2), т.е.
.
Элементарные делители всех многочленов (1) называют элементарными делителями матрицы A. Выпишем их:
(3)
Возьмем теперь жорданову матрицу J порядка n, составленную из жордановых клеток, определяемых следующим образом: каждому элементарному делителю матрицы A ставим в соответствие жорданову клетку порядка kij, относящуюся к числу лi. Очевидно, что отличными от 1 инвариантными множителями матрицы J - лE будут многочлены (1) и только они. Поэтому матрицы А - лE и J - лE эквивалентны и, следовательно, матрица А подобна жордановой матрице J [14, с. 315].
Рассмотрим необходимое и достаточное условие приводимости матрицы к диагональному виду.
Теорему: Матрица A порядка n с элементами из поля P тогда и только тогда приводится к диагональному виду, если все корни последнего инвариантного множителя en (л) ее характеристической матрицы лежат в поле P, причем среди этих корней нет кратных.
Подобные документы
Основные операции над матрицами и их свойства. Произведение матриц или перемножение матриц. Блочные матрицы. Понятие определителя. Панель инструментов Матрицы. Транспонирование. Умножение. Определитель квадратной матрицы. Модуль вектора.
реферат [109,2 K], добавлен 06.04.2003Применение матриц и их виды (равные, квадратные, диагональные, единичные, нулевые, вектор-строка, вектор-столбец). Примеры действий над матрицами (умножение на число, сложение, вычитание, умножение и транспонирование матриц) и свойства полученных матриц.
презентация [74,7 K], добавлен 21.09.2013Форма записи и методы решения системы алгебраических уравнений с n неизвестными. Умножение и нормы векторов и матриц. Свойства определителей матрицы. Собственные значения и собственные векторы. Примеры использования числовых характеристик матриц.
реферат [203,0 K], добавлен 12.08.2009Понятие, типы и алгебра матриц. Определители квадратной матрицы и их свойства, теоремы Лапласа и аннулирования. Понятие обратной матрицы и ее единственность, алгоритм построения и свойства. Определение единичной матрицы только для квадратных матриц.
реферат [296,6 K], добавлен 12.06.2010Интерпретация ортогональной и унитарной матрицы. Основные детерминанты матриц. Определение комплексных квадратных невырожденных и вырожденных матриц. Методы нахождения определителя. Метод конденсации Доджсона. Кососимметричная полилинейная функция строк.
курсовая работа [620,9 K], добавлен 04.06.2015Расчет денежных расходов предприятия на выпуск изделий, при выражении их стоимости при помощи матриц. Проверка совместимости системы уравнений и их решение по формулам Крамера и с помощью обратной матрицы. Решение алгебраических уравнений методом Гаусса.
контрольная работа [576,6 K], добавлен 28.09.2014Примеры алгебраических групп матриц, классические матричные группы: общая, специальная, симплектическая и ортогональная. Компоненты алгебраической группы. Ранг матрицы, возвращение к уравнениям, совместимость. Линейные отображения, действия с матрицами.
курсовая работа [303,7 K], добавлен 22.09.2009Обратимые матрицы над полем Zp. Формула для подсчета обратимых матриц порядка 2. Формула для подсчета обратимых матриц порядка 3. Общая формула подсчета обратимых матриц над полем Zp. Обратимые матрицы над Zn.
дипломная работа [156,7 K], добавлен 08.08.2007Методика расчета скалярного произведения заданных векторов. Расчет определителей и рангов матриц, нахождение обратных матриц. Разрешение уравнений по методу Крамера, обратной матрицы, а также встроенной функции lsolve. Анализ полученных результатов.
лабораторная работа [86,8 K], добавлен 13.10.2014Базовые действия над матрицами. Решение матричных уравнений с помощью обратной матрицы и с помощью элементарных преобразований. Понятия обратной и транспонированной матриц. Решение матричных уравнений различных видов: АХ=В, ХА=В, АХВ=С, АХ+ХВ=С, АХ=ХА.
курсовая работа [172,0 K], добавлен 09.09.2013