Решение полиномиальных матричных уравнений
Приведенная ниже функция вычисляет собственные значения матричного полинома.
где степень полинома р — целое неотрицательное число, а А 0 , А 1 ,..., А p — входные матрицы порядка п. Выходная матрица X размера nхnр содержит собственные векторы в столбцах. Вектор е размером пр содержит собственные значения.
Пример:
» А=[1:4:5:8:9:12:13:16]
А =
|
|
|
|
|
|
1 2 |
3 4 |
|
|
|
5 6 |
7 8 |
|
|
|
9 10 |
11 12 |
|
|
|
1314 |
15 16 |
|
|
|
» В=[4:7 |
;2:5;10:13;23:26] |
|
|
|
3 - |
|
|
|
|
4 5 |
6 7 |
|
|
|
2 3 |
4 5 |
|
|
|
1011 |
12 13 |
|
|
|
2324 |
25 26 |
|
|
|
» [F.a]= |
polyeig(A.B) |
|
|
|
F = |
|
|
|
|
0.4373 |
0.0689 |
-0.5426 |
-0.7594 |
|
-0.3372 |
-0.4969 |
0.6675 |
-0.1314 |
|
-0.6375 |
0.7870 |
0.2927 |
-0.1314 |
|
0.5374 |
-0.3591 |
- 0.4176 |
0.3771 |
|
a = |
|
|
|
|
4.4048 |
|
|
|
|
0.4425 |
|
|
|
|
-0.3229 |
|
|
|
|
-1.0000 |
|
|
|
|
|
|