Иллюстрированный самоучитель по Matlab

         

Квадратичный метод сопряженных градиентов


Квадратичный метод сопряженных градиентов реализуется в системе MATLAB с помощью функции cgs:

cgs(A.B) — возвращает решение X СЛУ А*Х=В. А — квадратная матрица. Функция cgs начинает итерации от начальной оценки, по умолчанию представляющей собой вектор размера

п,

состоящий из нулей. Итерации производятся либо до сходимости метода, либо до появления ошибки, либо до достижения максимального числа итераций. Сходимость метода достигается, когда относительный остаток norm(B-A*X)/norm(B) меньше или равен погрешности метода (по умолчанию le-6). Функция cgs(...) имеет и ряд других форм записи, аналогичных описанным для функции bieg(...). Пример:

» cgs(A.B)

CGS converged at iteration 4 to a solution

with relative residual 4e-014 

ans =

1.0000

2.0000

3.0000

4.0000



Содержание раздела