МИНИСТЕРСТВО ОБЩЕГО И ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ НИЖЕГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ им. Н. И. ЛОБАЧЕВСКОГО Факультет вычислительной математики и кибернетики Кафедра математической логики и высшей алгебры ЖОРДАНОВА ФОРМА ЛИНЕЙНОГО ПРЕОБРАЗОВАНИЯ Решение задач Нижний Новгород, 2002 http://vmk.ucoz.net/
24
Embed
New jordan · 2014. 3. 3. · УДК: 512.1 Жорданова форма линейного преобразования. Решение задач: Методи-ческая разработка
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
МИНИСТЕРСТВО ОБЩЕГО И ПРОФЕССИОНАЛЬНОГООБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ
НИЖЕГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТим. Н. И. ЛОБАЧЕВСКОГО
Факультет вычислительной математики и кибернетики
Кафедра математической логики и высшей алгебры
ЖОРДАНОВА ФОРМАЛИНЕЙНОГО ПРЕОБРАЗОВАНИЯ
Решение задач
Нижний Новгород, 2002
http://vmk.ucoz.net/
УДК: 512.1
Жорданова форма линейного преобразования. Решение задач: Методи-ческая разработка по курсу “Геометрия и алгебра” для студентов специ-альностей “Прикладная математика и информатика”, “Прикладная ин-форматика” / Составители: С. И. Веселов, Н. Ю. Золотых, Т. Г. Смирнова,А. Ю. Чирков. – Н. Новгород: Нижегородский государственный универ-ситет, 2002. – 24 с.
Методическая разработка предназначена для студентов факультета ВМКспециальностей “Прикладная математика и информатика”, “Приклад-ная информатика” и содержит теоретическое введение, примеры реше-ния задач и контрольные задания по теме “Жорданова форма линейногопреобразования”.
Составители:
С. И. Веселов, к.ф.-м.н., доц. каф. МЛиВА,Н. Ю. Золотых, к.ф.-м.н., ст. преп. каф. МЛиВА,Т. Г. Смирнова, к.ф.-м.н., ст. преп. каф. МЛиВА,А. Ю. Чирков, к.ф.-м.н., доц. каф. МЛиВА.
Рецензент:
А. В. Баркалов, к.ф.-м.н., доц. каф. МО ЭВМ факультета ВМК.
Нижегородский государственный университет им. Н. И. Лобачевского
Итак, для нахождения жорданова базиса в данном примере следуетнайти H1; H2; H3, выбрать вектор g3 из H3 nH2, затем вычислить g2 и g1по формулам (1).
Пример 2. Известно, что матрица A подобна матрицеG = 0� � 0 00 � 10 0 � 1A :Найти жорданов базис.
Так же как в предыдущем примере рассмотрим подпространстваHk = �x : (A � �E)k = 0 (k = 1; 2; 3):Имеют место следующие соотношения:H1 � H2 = H3; dimH1 = 2; dimH2 = 3; g1; g2 2 H1; g3 2 H2:Так как g1; g2; g3 — базис, то g3 =2 H1.
Итак, для нахождения жорданова базиса в данном примере следуетнайти H1; H2, выбрать вектор g3 из H2 nH1 затем вычислить g2 и в ка-честве g1 выбрать любой вектор из H1 не пропорциональный g2.
Излагаемый далее общий метод, является естественным обобщениемприведенных выше рассуждений.
4
Алгоритм построения жорданова базиса
1) Определить собственные числа �1; : : : ; �s.
2) Выбирая в качестве � поочередно каждое из собственных чисел,выполнить следующие действия.
а) Найти Hi = ker(A� �E)i; (i = 1; : : : ; r), гдеr = min fk j Hk = Hk+1g :б) Составить таблицу по схеме, приведенной на рисунке (1), в
которой первая сверху строка содержит линейно независимыевекторы, для которых выполняется равенствоHr = Hr�1 � L(a1; : : : ; as1); (2)вторая строка содержит линейно независимые векторы, длякоторых Hr�1 = Hr�2 � L(b1; : : : ; bs2) (3)и т. д. Последняя строка содержит базис H1, так, что sr =dimH1. Списки bs1+1; : : : ; bs2 , . . . , sr�1+1; : : : ; sr могут бытьпустыми.
Покажем, что построение таблицы, главной особенностью которойявляются свойства (2, 3), возможно.
1) Так как Hr�1 � Hr, то в верхней cтроке можно выписать любойнабор векторов, дополняющих базис подпространстваHr�1 до ба-зиса подпространства Hr.
2) Векторы b1; : : : ; bs1 определяются ранее выбранными векторами,поэтому формирование второй строки возможно тогда и только то-гда, когда L(b1; : : : ; bs1)\Hr�2 = f0g. Это равенство справедливотак как �1b1 + : : :+ �s1bs1 2 Hr�2 ! �1a1 + : : :+ �s1as1 2 Hr�1! �1 = : : : = �s1 = 0. Таким образом, можно построить наборb1; : : : ; bs1 , для которого выполняется условие (3).
Структура последующих строк таблицы такая же как и у второй, по-этому приведенное в пункте 2 рассуждение применимо и далее (можнодоказать по индукции).
Для того, чтобы получить именно жорданов базис, можно располо-жить их в следующем порядке: 1; : : : ; b1; a1; 2; : : : ; b2; a2; : : : ; sr ;то есть, обходя столбцы таблицы слева направо, а каждый столбец —
снизу вверх. Порядок столбцов не имеет принципиального значения, нов каждом из них векторы следует нумеровать строго снизу вверх.
Жорданова форма содержит s1 клеток размера r, s2� s1 клеток раз-мера r � 1 и т. д.
Замечание 1. Пусть Si обозначает множество свободных неизвестныхсистемы уравнений для нахождения Hi. Так как H1 � H2; : : :, то можновыбирать свободные неизвестные так, чтобы выполнялись соотношенияS1 � S2; : : :. Использование этого замечания может уменьшить объемвычислений при формировании жорданова базиса.
Пример 3. Найти жорданов базис и жорданову форму линейного преоб-разования, заданного в естественном базисе fe1; e2; e3; e4; e5; e6g матри-цей A = 0BBBBBB� 1 2 0 0 0 21 �1 0 1 �1 �20 0 1 0 0 00 0 0 1 0 00 �2 0 �1 1 �2�1 2 0 �1 1 3
1CCCCCCA :Единственным собственным числом матрицы является � = 1. Решая
систему уравнений (A�E)x = 0, находим подпространствоH1 = L8>>>>>><>>>>>>:0BBBBBB� 01000�1
и решаем систему (A�E)2x = 0: (4)Имеем H2 = H1 + L8>>>>>><>>>>>>:
0BBBBBB� 000�4011CCCCCCA ;0BBBBBB� 000010
1CCCCCCA9>>>>>>=>>>>>>; :Так как (A � E)3 = 0, то H3 = R6 и процесс построения корневых под-пространств завершен. Вектор u = (0; 0; 0; 0; 0; 1)> не принадлежит H2,так как не удовлетворяет системе (4). Поскольку размерности H3 и H2различаются на единицу, то вектор u дополняет базис H2 до базиса H3.В соответствии с общей схемой вычисляем произведения (A � E)u =(2;�2; 0; 0;�2; 2)>, (A�E)2u = (0; 4; 0; 0; 0;�4)>, для удобства распо-лагая векторы друг под другом:(0; 0; 0; 0; 0; 1)>;(2;�2; 0; 0;�2; 2)>;(0; 4; 0; 0; 0;�4)>:Во второй строке можем добавить вектор v = (0; 0; 0;�4; 0; 1)>, по-скольку векторы (A � E)u и v дополняют базис H1 до базиса H2. Вы-числяем (A�E)v = (2;�6; 0; 0; 2; 6)> и добавляем его в третью строку(0; 0; 0; 0; 0; 1)>;(2;�2; 0; 0;�2; 2)>; (0; 0; 0;�4; 0; 1)>;(0; 4; 0; 0; 0;�4)>; (2;�6; 0; 0; 2; 6)>:
обратную к ней, вычисляем матрицу A2 = P�1A1P с элементамиa2ij и так далее.
2) Если на некотором шаге оказалось, что akn�k n�k�1 = 0, но 9 j <n � k � 1 такой,что akn�k j 6= 0, то надо переставить местами jи n � k � 1 строки, затем j и n � k � 1 столбцы (эта операцияравносильна перестановке базисных векторов). Получится матри-ца подобная исходной, переходим к 1.
3) Если на некотором шаге akn�k j = 0 8j < n � k, то матрица имеетвид � B C0 D � ;
10
где B и D квадратные матрицы, причем D — фробениусова мат-рица, поэтому задача сводится к вычислению характеристическогомногочлена для B.
Пример 4. Методом Данилевского привести матрицуA = 0BB� 4 1 1 1�1 2 �1 �16 1 �1 1�6 �1 4 2 1CCAк фробениусовой форме и найти характеристическое уравнение этой мат-рицы.
Для того, чтобы обойтись при вычислениях без дробей, переставимместами 3 и 4 столбцы, 3 и 4 строки:0BB� 4 1 1 1�1 2 �1 �16 1 1 �1�6 �1 2 4 1CCA! 0BB� 4 1 1 1�1 2 �1 �1�6 �1 2 46 1 1 �1 1CCA :
Вычитая третий столбец с подходящими множителями из всех осталь-ных столбцов, и затем, выполняя соответствующие преобразования надстроками, получим матрицу, в которой элемент a43 = 1, а остальные эле-менты четвертой строки равны нулю.0BB� �2 0 1 25 3 �1 �2�18 �3 2 60 0 1 0 1CCA! 0BB� �2 0 1 25 3 �1 �2�25 0 6 160 0 1 0 1CCA :
Так как в полученной матрице элемент a32 = 0 , а a31 6= 0, следуетпереставить местами столбцы 1 и 2, и соответственно строки 1 и 2:0BB� 0 �2 1 23 5 �1 �20 �25 6 160 0 1 0 1CCA! 0BB� 3 5 �1 �20 �2 1 20 �25 6 160 0 1 0 1CCA :
Матрица полураспалась на два блока, поэтому характеристическиймногочлен есть (�� 3)(�3 � 4�2 � 3�+ 18).
3. Значение функции от матрицы
Определение 4. Говорят, что функция f(�) определена на спектре
матрицы A, если для каждого собственного числа �0 определена какf(�), так и все ее производные до (k � 1) -ой включительно, где k равнократности корня �0 для минимального многочлена матрицы A.
Теорема 2. Если f(�) определена на спектре матрицы A, то суще-
ствует единственный многочлен If (�) наименьшей степени, зна-
чения которого на спектреA совпадают со значениями f(�). Сте-пень этого многочлена строго меньше степени минимального мно-
гочлена матрицы A.
Определение 5. Если f(�) определена на спектре матрицы A, то мат-рицу If (A) называют значением функции f(�) от матрицы A. Инымисловами, f(A) � If (A).Пример 5. Вычислить значение функции f(�) = e� от матрицыA = 0� 4 2 �56 4 �95 3 �7 1A :
Минимальный многочлен для этой матрицы совпадает с характери-стическим и равен �2(� � 1). Следовательно, степень многочлена If (�)
12
не превышает 2. Представим его в виде If (�) = a0 + a1� + a2�2 и най-дем коэффициенты, приравнивая значения функций e� и If (�) на спектрематрицы A: 8<: If (0) = e0;I 0f (0) = e0;If (1) = e1; 8<: a0 = 1;a1 = 1;a0 + a1 + a2 = e:
Таким образом,eA = E +A+ (e� 2)A2 = 0� 3e� 1 e �3e+ 13e e+ 3 �3e� 33e� 1 e+ 1 �3e 1A :4. Контрольные задания
В примерах 1.1–1.30 построить жорданов базис и найти жордановуформу матрицы A.
В примерах 1.31–1.45 построить жорданов базис и жорданову формуматрицы линейного преобразования ' пространства многочленов с ком-плексными коэффициентами степени не превосходящих n. Рассмотретьслучаи, когда n = 3; 4; 5; 6 .
0BB� �5 1 2 37 0 �3 �31 2 �1 1�11 1 5 6 1CCA :В примерах 3.1–3.30 вычислить значения функций от матриц.
20
3.1. lnA, где A = 0� 3 �1 �4�4 3 82 �1 �3 1A.
3.2. sin ��2A�, где A = 0� �1 3 �3�1 3 �12 �2 4 1A.
3.3. eA, где A = 0� �1 �2 53 4 �91 1 �2 1A.
3.4. os ��2A�, где A = 0� 4 1 29 4 6�9 �3 �5 1A.
3.5. ar sinA, где A = 0� �1 �2 74 5 �161 1 �3 1A.
3.6. ar osA, где A = 0� 2 �1 �4�4 2 82 �1 �4 1A.
3.7. lnA, где A = 0� �2 3 �3�1 2 �12 �2 3 1A.
3.8. sin ��2A�, где A = 0� �1 �2 �22 3 21 1 2 1A.
3.9. eA, где A = 0� 1 0 �13 4 �151 1 �4 1A.
3.10. os ��2A�, где A = 0� 1 1 1�8 �5 �44 2 1 1A.
3.11. ar sinA, где A = 0� 2 �8 �5�1 5 33 �13 �8 1A.
21
3.12. sin ��4A�, где A = 0� 6 �1 �8�4 2 84 �1 �6 1A.
3.13. ar osA, где A = 0� 1 2 5�3 �6 �151 2 5 1A.
3.14. os (�A), где A = 0� 8 �1 �10�4 2 85 �1 �7 1A.
3.15. lnA, где A = 0� 2 �2 �1�1 3 13 �6 �2 1A.
3.16. sin ��2A�, где A = 0� 3 �2 �1�1 4 13 �6 �1 1A.
3.17. eA, где A = 0� 3 �2 51 0 1�1 1 �2 1A.
3.18. os ��2A�, где A = 0� 6 �5 52 �1 2�3 3 �2 1A.
3.19. ar sinA, где A = 0� 5 �8 �14�1 2 33 �4 �8 1A.
3.20. sin ��4A�, где A = 0� 1 3 3�1 1 �1�2 �2 �4 1A.
3.21. ar osA, где A = 0� 5 10 15�4 �8 �121 2 3 1A.
3.22. lnA, где A = 0� 4 �2 13 �1 1�3 2 0 1A.
22
3.23. sin ��2A�, где A = 0� �4 4 6�3 4 3�3 2 5 1A.
3.24. eA, где A = 0� 2 �2 �3�2 2 91 �1 �3 1A.
3.25. os ��2A�, где A = 0� �2 1 14 �5 �4�5 5 4 1A.
3.26. ar sinA, где A = 0� �7 2 �51 0 111 �3 8 1A.
3.27. os�A, где A = 0� �4 �1 2�4 2 8�1 �1 �1 1A.
3.28. ar osA, где A = 0� �4 1 �28 �2 412 �3 6 1A.
3.29. sin ��4A�, где A = 0� �7 3 �9�1 1 �16 �2 8 1A.
3.30. os � �A�, где A = 0� 3 3 6�1 1 �1�4 �2 �7 1A.
Литература
1. Воеводин В. В. Линейная алгебра. — М.: Наука, 1974.2. Беклемишев Д. В. Курс аналитической геометрии и линейной
алгебры. — М.: Наука, 1983.3. Гантмахер Ф. Р. Теория матриц. — М.: Наука, 1988.
23
Содержание
1. Жорданова форма линейного преобразования 3
2. Метод Данилевского 9
3. Значение функции от матрицы 12
4. Контрольные задания 13
Литература 23
Жорданова форма линейного преобразованияРешение задач
(Методическая разработка)
Составители:
С. И. Веселов, к.ф.-м.н., доц. каф. МЛиВА,Н. Ю. Золотых, к.ф.-м.н., ст. преп. каф. МЛиВА,Т. Г. Смирнова, к.ф.-м.н., ст. преп. каф. МЛиВА,А. Ю. Чирков, к.ф.-м.н., доц. каф. МЛиВА.
Подписано в печать Формат 60� 84 1=16.Бумага газетная. Печать офсетная. Усл. печ. л. 1,5. Тираж 250 экз.
Заказ
Нижегородский государственный университет им. Н. И. Лобачевского,603950, ННГУ, Н. Новгород, пр. Гагарина, 23
Типография ННГУ, Н. Новгород, ул. Б. Покровская, 37