на каждой сетке отличается от мажорантной. Убедиться, что фактическая погрешность не имеет вида Rah.

5. Найти погрешность нелинейной квадратурной формулы (33) на равномерной сетке.

6. Нцйти погрешность квадратурной формулы Филона (38), аналогичной формуле трапеций.

7. Для слабо меняющихся функций формулы средних и трапеций близки по точности. Почему их аналоги для быстро осциллирующих функций (35) и (37) имеют существенно разную точность?

8. Вывести формулу Филона, соответствующую квадратичной аппроксимации амплитуды по трем соседним узлам. Сравнить ее с формулой Симпсона.

9. Найти асимптотическое выражение погрешности квадратурной формулы (41).

10. Найти формулу для определения числа л способом Бюффона (§ 4, п. 4) и дисперсию этого способа; для этого удобно свести бросания к вычислению интеграла статистическими методами. Сделать то же для иголок, скрепленных крестом и снежинкой.

ГЛАВА V

СИСТЕМЫ УРАВНЕНИЙ

В главе V рассмотрены методы решения систем алгебраических уравнений. В § 1 изложено решение линейных систем методом исключения Гаусса, а также вычисление определителя и обращение матрицы; дан обзор других методов решения этих задач. В § 2 приведены различные методы нахождения корня одного трансцендентного уравнения. В §3 некоторые из этих методов обобщены на системы нелинейных уравнений.

§ 1. Линейные системы

1. Задачи линейной алгебры. Выделяют четыре основные задачи линейной алгебры: решение системы линейных уравнений Ах = Ь, где А--квадратная матрица и х, b - векторы; вычисление определителя; нахождение обратной матрицы; определение собственных значений и собственных векторов матрицы. В этом параграфе мы подробно рассмотрим первую задачу и попутно решим вторую и третью. Четвертая задача суще-.0 ственно сложней, и ей посвящена следующая глава.

Известно, что если det Л = О, то система линейных уравнений или не имеет решения, или имеет бесчисленное множество решений. Если же det Л О, то система имеет решение, притом единственное. Дальше мы будем рассматривать только последний случай..

Все эти случаи хорошо иллюстрируются геометрически на системе двух уравнений (рис. 24). Каждому уравнению соответствует прямая в плоскости х, у, а точка пересечения этих прямых есть решение системы (для п уравнений решение есть точка пересечения всех п гиперплоскостей в п-мерном пространстве). Если det Л = О, то наклоны прямых равны, и они либо параллельны, либо совпадают. В противном случае прямые имеют единственную точку пересечения.

Рис. 24.

На практике кроме существования и единственности решения важна еще устойчивость относительно погрешностей правой части и элементов матрицы. Формально перепишем линейную систему в виде х = А~Ь. Варьируя это равенство и определяя вариацию обратной матрицы из соотношения 6Е = 6{АА")=АЬА~ -\- = 0, получим

бл; = Л-1 (бй-бЛ-х).

Формально устойчивость есть, ибо при det Л =7 О обратная матрица существует. Но если матрица А- имеет большие элементы, то можно указать такой вид погрешности исходных данных, который сильно изменит решение. В этом случае систему называют плохо обусловленной (по-видимому, плохая обусловленность была известна еще Гауссу). Очевидно, у плохо обусловленных систем det Л я« fO; однако заметим, что,этот признак плохой обусловленности является необходимым, но недостаточным.

Плохо обусловленная система геометрически соответствует почти параллельным прямым. При этом небольшое изменение наклона или сдвиг одной прямой сильно меняют положение точки пересечения (рис. 24, пунктир). В многомерном случае геометрическая картина может быть более сложной. Так, для трех переменных возможен случай плохой обусловленности, когда соответствующие трем уравнениям плоскости пересекаются под большими углами (т. е. далеки от параллельности), но линии их попарного пересечения почти параллельны.

В теоретических исследованиях обусловленность часто характеризуют числом х= /4 • Это число зависит от того, какая норма-матриц

выбрана, но при любой норме xSl. Чем больше это число,,тем хуже обусловленность системы; обычно х~103-10 уже означает плохую обусловленность.

В практических расчетах этим определением плохой обусловленности пользуются редко, ибо для его проверки надо находить обратную матрицу, что при плохо обусловленной матрице А нелегко сделать. Чаще ограничиваются проверкой условия det=aO, хотя оно является необходимым, но недостаточным, что видно из простого примера. Положим А - еЕ, где «Б - единичная матрица; тогда det Л = 8", и даже при не очень малых е детерминант высокого порядка п очень мал. Но система с диагональной матрицей хорошо обуслов,гена, и для нее критерий х= Л • = 1 наиболее благоприятен.

Методы решения линейных систем делятся на прямые и итерационные. Прямые методы дают решение за конечное число действий, просты и наиболее универсальны; они рассматриваются в этом параграфе. Для систем небольшого порядка п<1200 применяются практически только прямые методы. Итерационные методы приведены в § 3; они выгодны для систем специального вида, со слабо заполненной матрицей очень большого порядка пя«10-т-10. Сравнительно недавно для решения плохо обусловленных систем стали применять методы регуляризации.