25. Решение полиномиальных матричных уравнений

 

Решение полиномиальных матричных уравнений

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

  • [Х.е] = polyeig(AO,Al,...Ap) — решает задачу собственных значений для матричного полинома степени р вида:

где степень полинома р — целое неотрицательное число, а А 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



 

11.gif

Изображение: