смущает то что по книжке получается что искомые элементы должны быть полиномами от элементов данной диогональной матрицы. Тоесть первый элемент равен полиному от первого элемента данной матрицы, второй - полиному от второго, ... ?
Я книжку не читал (хотя и заглянул), но это какая-то путаница.
Если все диагональные элементы различны, то полином степени
от
даёт общий вид перестановочной матрицы.
Следует это хотя бы из того, что подбором коэффициентов многочлена можно будет получить какие угодно наборы собственных чисел матрицы
(взяв в качестве полинома соответствующий интерполяционный, построенный по собственным числам исходной матрицы
).
Кстати это единственное чем Жарданова форма лучше диогональной?
Она не лучше, а хуже (сложнее) -- просто она более общая.
Добавлено спустя 3 минуты 50 секунд:Крестиком я обозначил произведение Кронекера, "vec" это векторизация матриц. Посмотрите в инете определение этих операций и всё сразу станет понятно.
Да, собственно, и без Кронекера понятно, что система уравнений -- линейная. И этот способ действительно хорош тем, что логически проще и, кроме того, для
конкретной матрицы решение получается за конечное количество шагов (в отличие от диагонализации).
Только он совершенно непрозрачен.