10.07.2015 Views

С.П. Шарый - Институт вычислительных технологий СО РАН

С.П. Шарый - Институт вычислительных технологий СО РАН

С.П. Шарый - Институт вычислительных технологий СО РАН

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

334 3. Численные методы линейной алгебрычастей на некоторую матрицу Λ. <strong>С</strong>ама эта матрица Λ называетсяпредобуславливающей матрицей или, коротко, предобуславливателем.Цель предобуславливания — изменение (вообще говоря, улучшение)свойств матрицы A исходной системы Ax = b, вместо которой мы получаемсистему(ΛA)x = Λb.<strong>П</strong>родуманный выбор предобуславливателя может, к примеру, изменитьвыгодным нам образом расположение спектра матрицы A, так необходимоедля организации сходящихся итерационных процессов.Естественно выполнить предобуславливание до перехода к системе(3.96), т. е. до прибавления вектора неизвестных x к обеим частям исходной<strong>С</strong>ЛАУ. <strong>П</strong>оскольку тогда вместоAx = b будем иметь(ΛA)x = Λb,то далее получаемx = (I −ΛA)x+Λb.Теперь в этом рекуррентном виде с помощью подходящего выбора Λможно добиваться требуемых свойств матрицы (I −ΛA).Каким образом следует выбирать предобуславливатели? <strong>С</strong>овершеннообщего рецепта на этот счёт не существует, и теория разбиваетсяздесь на набор рекомендаций для ряда более или менее конкретныхважных случаев.Например, если в качестве предобуславливающей матрицы взятьΛ = A −1 или хотя бы приближённо равную обратной к A, то вместосистемы Ax = b получим (A −1 A)x = A −1 b, т. е. систему уравненийIx = A −1 bили близкую к ней, матрица которой обладает всеми возможными достоинствами(хорошим диагональным преобладанием, малой обусловленностьюи т. п.). Ясно, что нахождение подобного предобуславливателяне менее трудно, чем решение исходной системы, но сама идеяпримера весьма плодотворна. На практике в качестве предобуславливателейчасто берут несложно вычисляемые обратные матрицы к тойили иной «существенной» части матрицы A. Например, к главной диагоналиматрицы или же к трём диагоналям — главной и двум побочным.Другой способ приведения <strong>С</strong>ЛАУ к рекуррентному виду основан нарасщеплении матрицы системы.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!