Решал задачу извлечения матричного корня 5-ой степени из положительной стохастической матрицы
P 4х4 стандартным методом через жорданову нормальную форму, к-рая оказалось тривиально диагональной с положительными числами. Ошибка полученного решения, т.е. 2-норма разности (
P -
), имеет порядок
. Теоретически тот же корень можно искать как решение задачи минимизации этой нормы с соотв. ограничениями, однако
fmincon НЕ находит нужное решение даже в в том случае, когда в качестве стартовой точки алгоритма выбиралось найденное ранее значение корня (уходит в сторону от истинного решения)!
Кто-нибудь знает разгадку ?