Квазиминимизация невязки — функция qmr
Метод решения СЛУ с квазиминимизацией невязки реализует функция qmr:
qmr (А, В) — возвращает решение X СЛУ А*Х=Ь. Матрица А должна быть квадратной. Функция qmr начинает итерации от начальной оценки, представляющей собой вектор длиной
п,
состоящий из нулей. Итерации производятся либо до сходимости метода, либо до появления ошибки, либо до достижения максимального числа итераций. Максимальное число итераций — минимум из п и 20. Функция qmr(...) имеет и ряд других форм записи, аналогичных описанным для функции bi eg (...). Пример:
» qmr(A.B)
QMR converged at iteration 4 to a solution
with relative residual 1.1e-014
ans =
1.0000
2.0000
3 .0000
4.0000