10.07.2015 Views

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

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

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

3.8. Метод прогонки 321<strong>П</strong>редложение 3.8.1 Если в системе линейных алгебраических уравненийс трёхдиагональной матрицей (3.84) имеет место диагональноепреобладание, т.е.|b i | > |a i |+|c i |, i = 1,2,...,n,то метод прогонки c выбором начальных значений согласно (3.89) являетсяреализуемым.<strong>П</strong>о поводу формулировки <strong>П</strong>редложения 3.8.1 можно заметить, чтоусловие диагонального преобладания в матрице влечёт её строгую регулярность,как мы видели в §3.6д. <strong>П</strong>оэтому в силу Теоремы 3.6.2 существуетLU-разложение такой матрицы, и оно может быть получено спомощью прямого хода метода Гаусса без перестановки строк и столбцов.Но это и означает реализуемость метода прогонки. Ниже, тем неменее, даётся другое доказательство этого факта, которое позволяетпомимо установления реализуемости дать ещё числовые оценки «запасаустойчивости» прогонки, т. е. того, насколько сильно знаменателивыражений (3.87)–(3.88) для прогоночных коэффициентов отличны отнуля в зависимости от элементов матрицы <strong>С</strong>ЛАУ.Доказательство. <strong>П</strong>окажем по индукции, что в рассматриваемой реализациипрогонки для всех индексовiсправедливо неравенство|ξ i | < 1.<strong>П</strong>режде всего, ξ 1 = 0 и потому база индукции выполнена: |ξ 1 | < 1.Далее, предположим, что для некоторого индекса i уже установленаоценка |ξ i | < 1. Если соответствующее c i = 0, то из (3.87) следуетξ i+1 = 0, и индукционный переход доказан. <strong>П</strong>оэтому пусть c i ≠ 0.Тогда справедлива следующая цепочка соотношений|ξ i+1 | =∣ − c ∣i ∣∣∣=a i ξ i +b i|c i ||a i ξ i +b i |≤|c i |∣ |bi |−|a i |·|ξ i | ∣ из оценки снизу для модуля суммы

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

Saved successfully!

Ooh no, something went wrong!