Решение полиномиальных матричных уравнений
Приведенная ниже функция вычисляет собственные значения матричного полинома.
- [Х.е] = polyeig(AO,Al,...Ap) — решает задачу собственных значений для матричного полинома степени р вида:
где степень полинома р — целое неотрицательное число, а А 0 , А 1 ,..., А p — входные матрицы порядка п. Выходная матрица X размера nхnр содержит собственные векторы в столбцах. Вектор е размером пр содержит собственные значения.
Пример:
» А=[1:4:5:8:9:12:13:16]
А =
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-0.5426 |
-0.7594
|
|
|
|
|
-0.1314
|
|
|
|
|
-0.1314
|
|
|
|
|
0.3771
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|