Как найти многочлен если задана матрица

Определение 14. Пусть дан многочлен

и пусть

– квадратная матрица, тогда значением
многочлена

от матрицы

называется матрица

,
где

– единичная матрица,

– матрица, получающаяся при умножении
матрицы

на себя

раз.

827 (П).
Найти значение многочлена

от матрицы


.

Р е ш е н и е.

Найдем

.


;


.

Ответ:

.

    1. Обратная матрица

Определение 15. Матрица

называется обратной к квадратной
матрице

,
если

.

Определение 16. Квадратная матрица

называется невырожденной, если
она имеет единственную обратную матрицу


.
В противном случае

вырожденная матрица.

Утверждение. Квадратная матрица

порядка

является невырож-денной в том и только
том случае, если определитель этой
матрицы отличен от нуля.

Для отыскания обратной матрицы
существуют два способа.

  1. Припишем
    к матрице

    справа единичную матрицу и, применяя
    метод Гаусса (см. §5), преобразуем
    расширенную матрицу так, чтобы слева
    стояла единичная матрица, тогда справа
    будет находиться обратная матрица

    :


.


.

Обоснование этого способа состоит в
следующем.

Пусть нам дана невырожденная квадратная
матрица. Задачу нахождения обратной
матрицы можно рассматривать как задачу
решения матричного уравнения

,
которое эквивалентно системе

уравнений с

неизвестными.

Эта система является объединением

систем уравнений, каждая из которых
содержит

неизвестных. Умножая поочередно строки
матрицы

на 1-й столбец матрицы

и приравнивая к 1-му столбцу матрицы

,
получим систему уравнений, матричная
форма записи которой имеет вид

(1.2.1)

С помощью элементарных операций над
строками матрицы систему уравнений
можно привести к виду

Умножая
поочередно строки матрицы

на второй столбец матрицы

и приравняв ко второму столбцу матрицы

,
получим систему уравнений


.
(1.2.2)

С помощью тех же элементарных операций,
что применялись для решения системы
(1.2.1), мы приведем систему (1.2.2) к виду

и т.д.

Поэтому
для нахождения обратной матрицы и был
предложен описанный выше способ.


  1. ,
    где

    – алгебраические дополнения к элементу

    ,

    – определитель матрицы

    (см. §2).

840 (П).
Найти обратную матрицу для матрицы

.

Р е ш е н и е.

I
способ.


.

Ответ:

.

II
способ.


;


;


;

;


;


;


;


;


;


.

Таким образом,


.

Ответ:

.

861 (П).
Решить матричное уравнение

.

Р е ш е н и е.

1 вариант.

Пусть

,
тогда


.


.

Ответ:

.

2 вариант.

Очевидно, что

.
Найдем матрицу, обратную к матрице

.

I
способ:


.

II
способ:


.

Таким образом,


.

Ответ:

.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]

  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #

Многочлены от матриц

Напомним определение многочлена от матрицы. Пусть заданы многочлен (степени m) переменной lambda:

f(lambda)=a_mcdot lambda^m+a_{m-1}cdot lambda^{m-1}+ldots+a_1cdot lambda+a_0,

(7.40)

где A — квадратная матрица n-го порядка. Выражение вида

f(A)=a_mcdot A^m+a_{m-1}cdot A^{m-1}+ldots+a_1cdot A+a_0cdot E,

(7.41)

называется многочленом от матрицы A.

При больших значениях m и n вычисление выражения (7.41) затруднительно из-за операции возведения матрицы в натуральную степень. Поэтому требуется найти другие, эквивалентные определению (7.41), формы записи и алгоритмы эффективного вычисления многочлена от матрицы. Для упрощения (7.41) имеются две возможности. Во-первых, можно упростить матрицу A так, чтобы многочлен (7.40) от упрощенной матрицы уже вычислялся сравнительно просто. Например, выражение (7.41) легко вычисляется, если матрица A диагональная. Во-вторых, можно понизить степень m многочлена, тогда самая трудоемкая операция — возведение матрицы в степень — упрощается.


Использование жордановой формы для нахождения многочлена от матрицы

Использование жордановой формы матрицы для нахождения многочлена от матрицы основано на трех свойствах.

1. Многочлены от подобных матриц подобны.

Действительно, пусть при помощи преобразования подобия матрица A приведена к жордановой форме J_A=S^{-1}AS. Подставим A=SJ_AS^{-1} в правую часть (7.41):

f(A)=a_mcdot(SJ_AS^{-1})^m+ a_{m-1}cdot(SJ_AS^{-1})^{m-1}+ldots+ a_1cdot(SJ_AS^{-1})+a_0cdot E

Учитывая, что (SJ_AS^{-1})^k= underbrace{SJ_AS^{-1}cdot SJ_AS^{-1} cdotldotscdot SJ_AS^{-1}}_{k}= SJ_A^kS^{-1} для любого натурального k, получаем

f(A)=ScdotBigl(a_mJ_A^m+a_{m-1}J_A^{m-1}+ldots+a_1J_A+a_0EBigr)cdot S^{-1}=Scdot f(J_A)cdot S^{-1}.

Таким образом, многочлены f(A) и f(J_A) подобны (с той же самой преобразующей матрицей S):

A=Scdot J_Acdot S^{-1}qquad Rightarrowqquad f(A)=Scdot f(J_A)cdot S^{-1}.

2. Многочлен от блочно-диагоналъной матрицы является блочно-диагоналъной матрицей.

Пусть A=begin{pmatrix}A_1!!&vline!!&O\hline O!!&vline!!&A_2end{pmatrix}, где A_1 и A_2 — квадратные матрицы, а O — нулевые матрицы соответствующих размеров. Для блочно-диагональных матриц справедливы равенства (они следуют из операций над блочными матрицами):

A^k=begin{pmatrix} A_1!!&vline!!&O\hline O!!&vline!!&A_2end{pmatrix}^k= begin{pmatrix} A_1^k!!&vline!!&O\hline O!!&vline!!&A_2^kend{pmatrix}!,quad a_kcdot A_k=begin{pmatrix} a_kA_1!!&vline!!&O\hline O!!&vline!!&a_kA_2end{pmatrix}!, где kinmathbb{N}

Поэтому f(A)=begin{pmatrix} f(A_1)!!&vline!!&O\hline O!!&vline!!& f(A_2) end{pmatrix}.Для большего числа блоков доказательство V. о 1Аг)) аналогичное.

3. Многочлен (7.41) от жордановой клетки J_r(lambda_0) имеет вид

f(J_r(lambda_0))= begin{pmatrix}f(lambda_0)& dfrac{1}{1!}f'(lambda_0)& cdots& dfrac{1}{(r-1)!}f^{(r-1)}(lambda_0)\[8pt] 0&f(lambda_0)& cdots& dfrac{1}{(r-2)!}f^{(r-2)}(lambda_0)\ vdots&vdots&ddots&vdots\ 0&0&cdots&f(lambda_0) end{pmatrix}!.

(7.42)

Это верхняя треугольная матрица r-го порядка, на главной диагонали которой стоят значения функции f(lambda) в точке lambda_0, над диагональю — значения первой производной в этой же точке и т.д., т.е. коэффициенты ряда Тейлора для функции f(lambda).

Действительно, разложим многочлен (7.40) по формуле Тейлора в окрестности точки lambda=lambda_0:

f(lambda)= f(lambda_0)+ frac{1}{1!},f'(lambda_0)(lambda-lambda_0)+ frac{1}{2!}, f''(lambda_0)(lambda-lambda_0)^2+ldots+ frac{1}{m!},f^{(m)} (lambda_0)(lambda-lambda_0)^m.

Остаточный член в данном случае равен нулю, так как все производные более высокого порядка, чем m, тождественно равны нулю. При вычислении f(J_r(lambda_0)) линейный двучлен (lambda-lambda_0) заменяется матрицей

I=Bigl(J_r(lambda_0)-lambda_0EBigr)= begin{pmatrix}lambda_0&1&cdots&0&0\ 0&lambda_0&ddots&0&0\ vdots&vdots&ddots&ddots&vdots\ 0&0&cdots&lambda_0&1\ 0& 0&cdots&0&lambda_0end{pmatrix}- begin{pmatrix}lambda_0&0&cdots&0&0\ 0&lambda_0 &ddots&0&0\ vdots&vdots&ddots&ddots&vdots\ 0&0&cdots&lambda_0&0\ 0&0&cdots& 0&lambda_0end{pmatrix}= begin{pmatrix}0&1&cdots&0&0\ 0&0&ddots&0&0\ vdots&vdots& ddots&ddots&vdots\ 0&0&cdots&0&1\ 0&0&cdots&0&0 end{pmatrix}!,

у которой элементы над главной диагональю равны единице, а остальные элементы равны нулю, т.е. I=begin{pmatrix}0&e_1&e_2&cdots&e_{r-1}end{pmatrix}, где e_i — i-й столбец единичной матрицы r-го порядка.

Можно показать, что при возведении в степень единичные элементы матрицы I смещаются вверх:

I^2=begin{pmatrix}0&0&e_1&e_2&cdots&e_{r-2}end{pmatrix},quad I^3=begin{pmatrix}0&0&0&e_1&cdots&e_{r-3}end{pmatrix} и т.д.

причем I^k — нулевая матрица при kgeqslant r. Подставляя эти матрицы в формулу Тейлора, получаем

f(A)= f(lambda_0)cdot E+frac{1}{1!},f'(lambda_0)cdot I+ frac{1}{2!}, f''(lambda_0)cdot I^2+ frac{1}{m!},f'^{(m)}(lambda_0)cdot I^m,

Складывая матрицы в правой части, получаем квадратную матрицу r-го порядка, у которой элементы главной диагонали равны f(lambda_0), элементы над главной диагональю равны — frac{1}{1!},f'(lambda_0) и т.д., т.е. матрицу вида (7.42).


Пример 7.16. Найти многочлен f(lambda)= lambda^2+ lambda+1 от матриц

A=begin{pmatrix}2&1&0\0&2&1\0&0&2end{pmatrix}!,qquad B=begin{pmatrix} 2&1&0\ 0&2&0\ 0&0&2 end{pmatrix}!,qquad C=begin{pmatrix}2&1&0\ 0&2&0\ 0&0&3 end{pmatrix}!.

Решение. а) Матрица A– это жорданова клетка 3-го порядка, соответствующая собственному значению 2: A=J_3(2). Находим значения функции и ее производных в точке lambda=2: f(2)=7, f'(2)=5, f''(2)=2. Составляем матрицу вида (7.42), учитывая, что r=3:

f(A)=f(J_3(2))= begin{pmatrix}f(2)&dfrac{1}{1!},f'(2)&dfrac{1}{2!},f''(2)\[8pt] 0&f(2)&dfrac{1}{1!},f'(2)\[4pt] 0&0&f(2) end{pmatrix}= begin{pmatrix}7&5&1\ 0&7&5\ 0&0&7 end{pmatrix}!.

Матрица B имеет жорданову форму B=operatorname{diag} (J_2(2),,J_1(2)), т.е. является блочно-диагональной. По свойству 2 многочлен от матрицы B является блочно-диагональной матрицей. Записываем многочлен f(lambda) от каждой жордановой клетки по формуле (7.42):

f(J_2(2))=begin{pmatrix}7&5\0&7end{pmatrix}!,qquad f(J_1(2))=f(2)=7.

Здесь число 7 рассматривается как квадратная матрица 1-го порядка. Составляем из этих квадратных матриц искомую блочно-диагональную матрицу

f(B)= operatorname{diag}!left(begin{pmatrix}7&5\0&7end{pmatrix}!,,(7)right)= begin{pmatrix} 7&5!!&vline!!&0\ 0&7!!&vline!!&0\hline 0&0!!&vline!!&7 end{pmatrix}!.

Матрица C имеет жорданову форму C=operatorname{diag} (J_2(2),J_1(3)), т.е. является блочно-диагональной. Записываем многочлен f(lambda) от каждой жордановой клетки по формуле (7.42):

f(J_2(2))=begin{pmatrix}7&5\0&7end{pmatrix}!,qquad f(J_1(3))=f(3)=13.

Составляем из этих квадратных матриц искомую блочно-диагональную матрицу

f(C)= operatorname{diag}!left(begin{pmatrix}7&5\0&7end{pmatrix}!,,(13)right)= begin{pmatrix} 7&5!!&vline!!&0\ 0&7!!&vline!!&0\hline 0&0!!&vline!!&13 end{pmatrix}!.


Первый способ нахождения многочлена от матрицы

1. Привести матрицу A к жордановой форме J_A=S^{-1}AS, т.е. определить жорданову форму J_A и преобразующую матрицу S.

2. Составить блочно-диагональную матрицу f(J_A), размещая на ее диагонали многочлены от жордановых клеток (7.42).

3. Найти многочлен от матрицы А по формуле f(A)=Scdot f(J_A)cdot S^{-1}.

Пример 7.17. Найти многочлен f(lambda)=lambda^m (при m>1) от матриц:

A=begin{pmatrix} 4&4\-1&0 end{pmatrix}!;quad B=begin{pmatrix} 1&0&1\ 0&1&-1\ -1&-1&1 end{pmatrix}!;quad C=begin{pmatrix} 1&1&1\ 1&1&1\ 1&1&1 end{pmatrix}!;quad D=begin{pmatrix} -2&5&-3\ -2&5&-3\ 1&-1&0end{pmatrix}!.

Решение. Матрица A. 1. Жорданова форма J_A и преобразующая матрица S были найдены в примере 7.15:

J_A=begin{pmatrix}2&1\0&2end{pmatrix}!,qquad S=begin{pmatrix}2&1\-1&0 end{pmatrix}!,qquad S^{-1}=begin{pmatrix}0&-1\1&2 end{pmatrix}!.

2. Жорданова форма J_A состоит из одной жордановой клетки J_A=J_2(2) 2-го порядка, соответствующей собственному значению lambda=2. Найдем значения функции f(lambda)=lambda^m и ее производной f'(lambda)=m lambda^{m-1} при lambda=2: f(2)=2^m, f'(2)=m2^{m-1}. Запишем многочлен от жордановой формы (блочно-диагональную матрицу с одним блоком): f(J_A)= begin{pmatrix} 2^m&m2^{m-1}\ 0&2^mend{pmatrix}.

3. Найдем многочлен от матрицы A:

begin{aligned}f(A)&= begin{pmatrix}2&1\-1&0end{pmatrix}!cdot! begin{pmatrix} 2^m&m2^{m-1}\ 0&2^m end{pmatrix}!cdot! begin{pmatrix}0&-1\1&2end{pmatrix}= begin{pmatrix} 2&1\-1&0 end{pmatrix}!cdot! begin{pmatrix}m2^{m-1}&(m-1)2^{m}\ 2^m&2^{m+1} end{pmatrix}=\[2pt] &= begin{pmatrix}(m+1)2^m&m2^{m+1}\ -m2^{m-1}&-(m-1)2^{m} end{pmatrix}!.end{aligned}

Матрица В. 1. Жорданова форма J_B и преобразующая матрица S были найдены в примере 7.15:

J_B=begin{pmatrix}1&1&0\ 0&1&1\ 0&0&1end{pmatrix}!,qquad S=begin{pmatrix}-1&0&1\ 1&0&0\ 0&-1&0end{pmatrix}!,qquad S^{-1}= begin{pmatrix}0&1&1\ 0&0&-1\ 1&1&0 end{pmatrix}!.

2. Жорданова форма J_B состоит из одной жордановой клетки J_B=J_3(1) 3-го порядка, соответствующей собственному значению lambda=1. Найдем значения функции f(lambda) и ее производных f'(lambda),,f''(lambda) при lambda=1: f(1)=1, f'(1)=m, f''(1)=m(m-1). Запишем многочлен от жордановой формы (r=3)colon

f(J_B)= begin{pmatrix} 1&m&frac{m}{2}(m-1)\ 0&1&m\ 0&0&1 end{pmatrix}!.

3. Найдем многочлен от матрицы B:

f(b)= begin{pmatrix}-1&0&1\ 1&0&0\ 0&-1&0end{pmatrix}!cdot! begin{pmatrix} 1&m&frac{m}{2}(m-1)\ 0&1&m\ 0&0&1end{pmatrix}!cdot! begin{pmatrix} 0&1&0\ 0&0&-1\ 1&1&0 end{pmatrix}= begin{pmatrix} 1-dfrac{m}{2}(m-1)&-dfrac{m}{2}(m-1)&m\[9pt] dfrac{m}{2}(m-1)& 1+dfrac{m}{2}(m-1)&-m\[4pt] -m&-m&1 end{pmatrix}!.

Матрица C. 1. Жорданова форма J_C и преобразующая матрица S были найдены в примере 7.15:

J_C=begin{pmatrix}0&0&0\0&0&0\0&0&3end{pmatrix}!,qquad S=begin{pmatrix} 1&0&-1\ 1&1&-1\ -2&-1&-1 end{pmatrix}!,qquad S^{-1}=begin{pmatrix} 2/3&-1/3&-1/3\ -1&1&0\ -1/3&-1/3&-1/3end{pmatrix}!.

2. Жорданова форма J_C состоит из трех жордановых клеток 1-го порядка J_C= operatorname{diag}(J_1(0),,J_1(0),,J_1(3)), соответствующих собственным значениям lambda=0 и lambda=3. Найдем значения функции f(lambda)=lambda^m при lambda=0 и lambda=3: f(0)=0, f(3)=3^m. Запишем многочлен от жордановой формы:

f(J_C)= operatorname{diag}Bigl(f(J_1(0)),,f(J_1(0)),,f(J_1(3))Bigr)= begin{pmatrix}0&0&0\ 0&0&0\ 0&0&3^m end{pmatrix}!.

3. Найдем многочлен от матрицы C:

begin{aligned}f(C)&= begin{pmatrix} 1&0&-1\ 1&1&-1\ -2&-1&-1end{pmatrix}!cdot! begin{pmatrix} 0&0&0\ 0&0&0\ 0&0&3^m end{pmatrix}!cdot! begin{pmatrix} 2/3&-1/3&-1/3\ -1&1&0\ -1/3&-1/3&-1/3end{pmatrix}=\[2pt] &=begin{pmatrix} 3^{m-1}&3^{m-1}&3^{m-1}\3^{m-1}&3^{m-1}&3^{m-1}\3^{m-1}&3^{m-1}&3^{m-1}end{pmatrix}= 3^{m-1}cdot! begin{pmatrix}1&1&1\ 1&1&1\ 1&1&1 end{pmatrix}!.end{aligned}

Результат совпадает с найденным в примерах 7.10, 7.12.

Матрица D. 1. Жорданова форма J_D и преобразующая матрица S были найдены в примере 7.15:

J_D=begin{pmatrix}0&1&0\0&0&0\0&0&3end{pmatrix}!;qquad S= frac{1}{9}!begin{pmatrix} 3&5&9\3&2&9\3&-1&0end{pmatrix}!;qquad S^{-1}=begin{pmatrix} 1&-1&3\3&-3&0\ -1&2&-1end{pmatrix}!.

2. Жорданова форма J_D состоит из двух жордановых клеток 2-го и 1-го порядков J_D=operatorname{diag}(J_2(0),,J_1(3)). соответствующих собственным значениям lambda=0 и lambda=3. Найдем значения функции f(0)=0, f(3)=3^m и производной f'(0)=0 (так как m>1). Запишем многочлен от жордановой формы

f(J_D)= operatorname{diag}Bigl(f(J_2(0)),,f(J_1(3))Bigr)= begin{pmatrix}0&0&0\ 0&0&0\ 0&0&3^m end{pmatrix}!.

3. Найдем многочлен от матрицы D:

f(D)= frac{1}{9}!begin{pmatrix} 3&5&9\3&2&9\ 3&-1&0 end{pmatrix}!cdot! begin{pmatrix} 0&0&0\ 0&0&0\ 0&0&3^m end{pmatrix}! cdot! begin{pmatrix} 1&-1&3\ 3&-3&0\ -1&2&-1 end{pmatrix}= 3^m! begin{pmatrix}-1&2&-1\ -1&2&-1\ 0&0&0 end{pmatrix}!.


Использование аннулирующих многочленов

Для понижения степени многочлена (7.41) можно использовать аннулирующие многочлены матрицы A, например, ее характеристический Delta_A(lambda)=det(A-lambda E) или минимальный mu_A(lambda) многочлены.

Обозначим через nu степень минимального многочлена

mu_A(lambda)= (lambda-lambda_1)^{m_1}cdot (lambda-lambda_2)^{m_2}cdot ldotscdot (lambda-lambda_k)^{m_k}.

Заметим, что nu не превосходит порядка n матрицы A (или, что то же самое, степени n характеристического многочлена Delta_A(lambda)), т.е. nu=m_1+ldots+m_kleqslant n. Разделим заданный многочлен (7.40) на минимальный:

f(lambda)= q(lambda)cdot mu_A(lambda)+r(lambda).

(7.43)

Здесь q(lambda) — частное, а r(lambda) — остаток, степень которого меньше nu:

r(lambda)= r_{nu-1}cdot lambda^{nu-1}+ r_{nu-2}cdot lambda^{nu-2}+ldots+ r_1cdot lambda+r_0.

(7.44)

Подставив в (7.43) вместо переменной lambda матрицу A, получим:

f(A)=r(A).

(7.45)

поскольку минимальный многочлен является аннулирующим (mu_A(lambda)=O).

Таким образом, вместо вычисления многочлена (7.41) степени m можно вычислить многочлен (7.44), степень которого меньше nu. Коэффициенты r_0,r_1,ldots,r_{nu-1} многочлена (7.44) находятся следующим образом.

Если все корни минимального многочлена простые, то, подставляя корень lambda_j в (7.43), получаем f(lambda_j)=r(lambda_j), так как mu_A(lambda_j)=0, т.е.

f(lambda_j)= r_{nu-1}lambda_j^{nu-1}+ r_{nu-2} lambda_j^{nu-2}+ldots+ r_1 lambda_j+r_0,quad j=1,2,ldots,nu.

Если lambda_j — корень минимального многочлена кратности m_j, учитывая, что

mu_A(lambda_j)=0,quad mu'_A(lambda_j)=0,quad ldots,quad mu_{A}^{(m_j-1)}(lambda_j)=0

из (7.43), последовательно дифференцируя, получаем

f(lambda_j)=r(lambda_j),quad left.{left(dfrac{d^i}{d lambda^i} f(lambda)right)}!right|_{lambda=lambda_j}= left.{left(dfrac{d^i}{d lambda^i} r(lambda)right)}!right|_{lambda=lambda_j},quad i=1,2,ldots,m_j-1.

(7.46)

Записывая равенства (7.46) для каждого корня минимального многочлена, получим совместную систему nu линейных уравнений с nu неизвестными r_0,r_1,ldots,r_{nu-1}.


Второй способ нахождения многочлена от матрицы

1. Найти минимальный многочлен mu_A(lambda) матрицы A одним из способов, рассмотренных в разд.7.2.4. Определить его степень mu и записать многочлен (7.44) с неопределенными коэффициентами r_0,r_1,ldots,r_{nu-1}:

r(lambda)= r_{nu-1}cdot lambda^{nu-1}+ r_{nu-2}cdot lambda^{nu-2}+ldots+ r_1cdot lambda+r_0.

2. Для каждого корня lambda_j (кратности m_j) минимального многочлена по формулам (7.46) составить m_j уравнений. Все уравнения объединить в одну систему.

3. Решить составленную систему, т.е. найти коэффициенты r_0,r_1,ldots,r_{nu-1} многочлена r(lambda).

4. По формуле (7.45) найти многочлен от матрицы:

f(A)=r(A)= r_{nu-1}cdot A^{nu-1}+ r_{nu-2}cdot A^{nu-2}+ldots+r_1cdot A+r_0cdot E.


Замечания 7.8.

1. Вместо минимального многочлена можно использовать характеристический многочлен матрицы, который также является аннулирующим (см. теорему Гамильтона-Кэли). При этом в пунктах 1,2 алгоритма минимальный многочлен заменяется характеристическим, степень которого равна n.

2. В первом способе нахождения многочлена от матрицы используются все инвариантные множители, так как нужно получить жорданову форму. Во втором способе требуется только один последний инвариантный множитель, который совпадает с минимальным многочленом. Можно сказать, что жорданова форма матрицы излишне информативна для решения поставленной задачи.


Пример 7.18. Найти (вторым способом) многочлен f(lambda)= lambda^m (при m>1) от матриц:

A=begin{pmatrix}4&4\-1&0end{pmatrix}!;quad B=begin{pmatrix}1&0&1\ 0&1&-1\ -1&-1&1 end{pmatrix}!;quad C=begin{pmatrix}1&1&1\ 1&1&1\ 1&1&1 end{pmatrix}!;quad D=begin{pmatrix}-2&5&-3\ -2&5&-3\ 1&-1&0 end{pmatrix}!.

Решение. Матрица A. 1. Для матрицы A в примере 7.15 были найдены инвариантные множители. Минимальный многочлен совпадает с последним инвариантным множителем. Поэтому mu_A(lambda)=e_2(lambda)=(lambda-2)^2. Степень nu минимального многочлена равна двум. Значит, многочлен (7.44) линейный: r(lambda)=r_1 lambda+r_0.

2. Для двойного корня lambda=lambda_1=2~(m_1=2) составляем уравнения (7.46):

begin{cases}f(2)=r(2),\ left.{dfrac{df(lambda )}{d lambda}}right|_{lambda=2}= left.{dfrac{dr(lambda )}{d lambda}}right|_{lambda=2},end{cases}Rightarrowquad begin{cases}2^m=r_1cdot2+r_0,\mcdot2^{m-1}=r_1.end{cases}

3. Решая систему, получаем r_1=m2^{m-1},~r_0=(1-m)2^m и r(lambda)=m2^{m-1}cdot lambda+(1-m)2^m.

4. Находим многочлен от матрицы A:

f(A)=r(A)= m2^{m-1}cdot! begin{pmatrix}4&4\-1&0end{pmatrix}+()2^mcdot! begin{pmatrix}1&0\0&1end{pmatrix}= begin{pmatrix}(m+1)2^m&m2^{m+1}\ -m2^{m-1}& -(m-1)2^mend{pmatrix}!.

Матрица B. 1. Инвариантные множители характеристической матрицы (B-lambda E) найдены в примере 7.15. Минимальный многочлен равен последнему инвариантному множителю: mu_B(lambda)=e_3(lambda)=(lambda-1)^3. Степень nu минимального многочлена равна 3. Значит, многочлен (7.44) — это квадратный трехчлен: r(lambda)=r_2 lambda^2+r_1 lambda+r_0.

2. Для тройного корня lambda=lambda_1=1~(m_1=3) составляем уравнения (7.46):

begin{cases} f(1)=r(1),\[4pt] left.{dfrac{df(lambda )}{d lambda}}right|_{lambda=1}= left.{dfrac{dr(lambda )}{d lambda}}right|_{lambda=1},\[8pt] left.{dfrac{d^2f(lambda )}{d lambda^2}}right|_{lambda=1}= left.{dfrac{d^2r(lambda )}{d lambda^2}}right|_{lambda=1}, end{cases}Rightarrowquad begin{cases}1^m=r_2cdot1^2+r_1cdot1+r_0,\ mcdot1^{m-1}=2r_2cdot1+r_1,\ m(m-1)cdot1^{m-2}=2r_2. end{cases}

3. Решая систему, получаем r_2=frac{m(m-1)}{2},~~r_1=m(2-m),~~r_0=frac{(m-1)(m-2)}{2} и

r(lambda)=frac{m(m-1)}{2}cdot lambda^2+ m(m-1)cdot lambda+frac{(m-1)(m-2)}{2}.

4. Вычисляя B^2=begin{pmatrix}0&-1&2\ 1&2&-2\ -2&-2&1 end{pmatrix}, записываем искомый многочлен:

begin{aligned} f(B)=r(B)&= frac{m(m-1)}{2}cdot! begin{pmatrix}0&-1&2\ 1&2&-2\ -2&-2&1end{pmatrix}+ m(2-m)cdot! begin{pmatrix}1&0&1\ 0&1&-1\ -1&-1&1 end{pmatrix}+ frac{(m-1)(m-2)}{2}cdot! begin{pmatrix}1&0&0\ 0&1&0\ 0&0&1end{pmatrix}=\[2pt] &= begin{pmatrix}1- dfrac{m(m-1)}{2}&-dfrac{m(m-1)}{2}&m\[9pt] dfrac{m(m-1)}{2}& 1+dfrac{m(m-1)}{2}&-m\[4pt] -m&-m&1 end{pmatrix}!.end{aligned}

Матрица C. 1. Минимальный многочлен найден в примере7.15: mu_C(lambda)=lambda(lambda-3). Степень nu многочлена равна 2. Следовательно, многочлен (7.44) имеет первую степень: r(lambda)=r_1 lambda+r_0.

2. Для каждого простого корня lambda=lambda_0=0 и lambda= lambda_2=3 записываем пер вое равенство из (7.46):

begin{cases}f(0)=r(0),\ f(3)=r(3),end{cases},Rightarrowquad begin{cases} 0^m=r_1cdot0+r_0,\ 3^m=r_1cdot3+r_0.end{cases}

3. Решая систему, получаем r_1=3^{m-1},~r_0=0 и r(lambda)=3^{m-1}lambda.

4. Находим многочлен от матрицы C:

f(C)=r(C)=3^{m-1}C= 3^{m-1}cdot! begin{pmatrix}1&1&1\ 1&1&1\ 1&1&1 end{pmatrix}!.

Найдем f(C), используя характеристический многочлен вместо минимального. Согласно пункту 1 замечаний 7.8, выполняем все действия второго способа, заменяя минимальный многочлен характеристическим.

1. Найдем характеристический многочлен матрицы C (см. при мер 7.11): Delta_C(lambda)=lambda^2(3-lambda). Это многочлен 3-й степени. Поэтому многочлен (7.44) будет 2-ой степени: r(lambda)=r_2 lambda^2+r_1 lambda+r_0.

2. Для двойного корня lambda=0 записываем два уравнения из (7.46), а для простого корня lambda=3 одно:

begin{cases} f(0)=r(0),\ left.{dfrac{df(lambda)}{d lambda}}right|_{lambda=0}= left.{dfrac{dr(lambda)}{d lambda}} right|_{lambda=0},\  f(3)=r(3), end{cases},Rightarrowquad begin{cases} 0^m=r_2cdot 0^2+r_1cdot0+r_0,\ mcdot0^{m-1}=2r_2cdot0+r_1,\ 3^m=r_2cdot 3^2+r_1cdot3+r_0.end{cases}

3. Решая систему, получаем r_2=3^{m-2},~ r_1=r_0=0 и r(lambda)=3^{m-2}lambda^2.

4. Вычисляя C^2=begin{pmatrix} 3&3&3\ 3&3&3\ 3&3&3 end{pmatrix}, записываем искомый многочлен:

f(C)=r(C)=3^{m-2}cdot! begin{pmatrix}3&3&3\ 3&3&3\ 3&3&3end{pmatrix}= 3^{m-1}cdot! begin{pmatrix}1&1&1\ 1&1&1\ 1&1&1 end{pmatrix}!.

Поскольку степень характеристического многочлена больше степени минимального многочлена (3>2), его применение менее эффективно.

Матрица D. 1. Инвариантные множители характеристической матрицы (D-lambda E) найдены в примере 7.15. Минимальный многочлен равен последнему инвариантному множителю: mu_D(lambda)=e_3(lambda)=lambda^2(lambda-3). Степень nu минимального многочлена равна 3. Значит, многочлен (7.44) — это квадратный трехчлен: r(lambda)=r_2 lambda^2+r_1 lambda+r_0.

2. Для двойного корня lambda=lambda_1=0~(m_1=2) записываем первые два равенства (7.46), а для простого корня lambda=lambda_2=3 (m_2=1) — первое равенство из (7.46). Получаем систему трех уравнений относительно коэффициентов квадратного трехчлена r(lambda):

begin{cases} f(0)=r(0),\ left.{dfrac{df(lambda)}{d lambda}}right|_{lambda=0}= left.{dfrac{dr(lambda)}{d lambda}}right|_{lambda=0},\ f(3)=r(3), end{cases}Rightarrowquad begin{cases}0^m=r_2cdot0^2+r_1cdot0+r_0,\  mcdot0^{m-1}=2r_2cdot0+r_1,\ 3^m=r_2cdot3^2+r_1cdot3+r_0. end{cases}

3. Решая систему, получаем r_2=3^{m-2},~r_1=r_0=0 и r(lambda)=3^{m-2}lambda^2.

4. Находим многочлен от матрицы D:

f(D)=r(D)=3^{m-2}D^2= 3^{m-2}cdot! begin{pmatrix}-2&5&-3\ -2&5&-3\ 1&-1&0 end{pmatrix}= 3^mcdot! begin{pmatrix}-1&2&-1\ -1&2&-1\ 0&0&0end{pmatrix}!.

Эта формула справедлива при m>1, так как при m=1 или m=0 в системе для нахождения коэффициентов многочлена r(lambda) появляются неопределенные выражения (0^0). Впрочем, для этих показателей степени (m=1,~m=0) многочлен f(D) легко находится по определению D^1=D,~D^0=E.

Все результаты совпадают с полученными в примере 7.17.

Математический форум (помощь с решением задач, обсуждение вопросов по математике).

Кнопка "Поделиться"

Если заметили ошибку, опечатку или есть предложения, напишите в комментариях.


01:26

Многочлен от матрицы

Тема: Линейная алгебра.Калькулятор для вычисления многочлена от матрицы.

Пример. Найти многочлен f(x)= x^2-5x+3 от матрицы

Пусть дан многочлен  и квадратная матрица .

Определение. Многочленом от матрицы  назовем выражение , где  — единичная матрица того же порядка, что и матрица .

Пример 12.

а)  Найти значение многочлена  от матрицы .

Решение.

Многочлен  имеет вид:

.

б)  Найти значение многочлена от матрицы .

Тема: Линейная алгебра.Калькулятор для вычисления многочлена от матрицы. $CUT$ Пример. Найти многочлен f(x)= x^2-5x+3 от матрицы Пусть дан многочлен и квадратная матрица . Определение. Многочленом от матрицы назовем выражение , где — единичная матрица того же порядка, что и матрица . Пример 12. а) Найти значение многочлена от матрицы . Решение. Многочлен имеет вид: . б) Найти значение многочлена от матрицы .

  • 1
  • 2
  • 3
  • 4
  • 5

Категория: Линейная алгебра | Просмотров: 58621 | | Теги: Многочлен от матрицы, Возведение матрицы в степень | Рейтинг: 3.3/14

Всего комментариев: 1

Порядок вывода комментариев:

Как найти многочлен матрицы

На чтение 2 мин Опубликовано 18.02.2022

Многочлен и матрица как аргумент

Расссмотрим в данном материале одну из трудоёмких задач в высшей математике, которая звучит так: Найти чему задан многочлен

если аргумент есть квадратная матрица, то есть

И если сам принцип вычисления понятен, особенно если вы в совершенстве поняли как умножать матрицы, то непосредственное вычисление, для меня лично считается рутиной, которую по возможности нужно избежать.

Сразу хотелось бы сказать, где этот калькулятор пригодится. Для учителей, преподавателей, для создателей учебников, для тех, кому необходимо создавать оригинальные задачи по данной теме.

Также пригодится для студентов или аспирантов которые пишут рефераты, курсовые, дипломы.

Для всех остальных, это легкий способ проверить ошибку в заданном примере, решить, без долгих промежуточных вычислений, поставленную задачу.

Когда калькулятор был написан, оказалось что сайты, которые были посвещенны этой теме, содержали ошибки в промежуточных вычислениях и как как результат были неверные.

Данный калькулятор, я надеюсь избавлен от ошибок и Вы сможете безопасно решать любые примеры.

Как и подавляющее большинство калькуляторов на этом сайте, значениями как коэффициентов полинома, так и элементов матрицы, могут быть комплексные значения.

Такого на конец 2017 года, больше нигде не найдете, не считая конечно специальных созданных математических программ.

Приступим к примерам?

Найти значение полинома от матрицы

Исходный полином f(x) (его коэффициенты)
Аргумент является квадратной матрицей с элементами
Многочлен
Переменная x=
Результат вычислений

Чему равен полином если

Многочлен
Переменная x=
Результат вычислений

Найти значение многочлена от комплексной матрицы

Источник

Характеристический полином
матрицы
A, вычисляется следующим образом:

| Aλ E |

где
E
– единичная матрица, размеры которой совпадают с размерами исходной матрицы
A.

Разберем подробнее приведенную выше формулу. Если матрица
A
задана в виде:

матрица

тогда выражение
Aλ E
имеет вид:

A-λE

Наконец, нам нужно найти определитель:

det(A-λE)

Раскрыв этот определитель, мы получим полином
n-ой степени
(n
– порядок исходной матрицы), зависящий от
λ:

P( λ )
=
cn λn
+
cn1 λn1
+ +
ci λi
+ +
c1 λ
+
c0

Поскольку для вычисления характеристического полинома, требуется нахождение определителя матрицы, то характеристический полином может быть найден только для квадратной матрицы.

Наш онлайн калькулятор находит
характеристический полином матрицы, причем в качестве элементов матрицы, можно вводить не только числа и дроби, но и параметры.

Добавить комментарий