Библиотека Eigen - метод наименьших квадратов

Я хочу использовать задачу наименьших квадратов с использованием библиотеки Eigen. Мои варианты 2,

  • sysAAA.jacobiSvd( Eigen::ComputeThinU | Eigen::ComputeThinV ).solve( sysBBB )
  • sysAAA.colPivHouseholderQr().solve( sysBBB );

Сначала я использовал первый, но он оказался очень медленным (1)(2).

Поэтому я пошел ко второму решению (другие методы не подходят для моего случая, потому что они требуют специальных матриц (2))

Дает ли colPivHouseholderQr().solve решение методом наименьших квадратов?

У меня сложилось впечатление, что это не так (3), но я хочу убедиться, прежде чем искать «обходной путь».

  1. http://forum.kde.org/viewtopic.php?f=74&t=102088
  2. http://eigen.tuxfamily.org/dox/TopicLinearAlgebraDecompositions.html
  3. http://eigen.tuxfamily.org/dox/TutorialLinearAlgebra.html#TutorialLinAlgLeastsquares

person dim_tz    schedule 03.11.2012    source источник


Ответы (1)


Да, ColPivHouseholderQr::solve() вычисляет решение методом наименьших квадратов.

person ggael    schedule 15.11.2012
comment
Спасибо большое ggael! Может быть, у вас есть справочная ссылка, в которой это указано, чтобы использовать ее в будущем для подобных вопросов? Я не мог найти что-то, что говорит об этом ясно. - person dim_tz; 22.11.2012
comment
Просто для протокола: GGael предоставил доказательство здесь forum.kde. org/ Спасибо, ggael! - person dim_tz; 15.01.2013