
- •Глава 3 основы численных методов
- •§ 3.1. Основные понятия линейной алгебры
- •3.1.1. Линейное пространство. Евклидово пространство.
- •3.1.2. Арифметическое пространство. Векторы и операции над ними. Линейная зависимость векторов. Базис.
- •3.1.3. Матрицы и операции над ними. Некоторые виды матриц.
- •3.1.4. Скалярное произведение векторов.
- •3.1.5. Норма вектора. Некоторые свойства векторных норм. Абсолютная и относительная погрешности. Сходимость.
- •3.1.6. Норма матрицы. Некоторые свойства матричных норм.
- •3.1.7. Число обусловленности матрицы.
- •§ 3.2. Прямые методы решения систем линейных алгебраических уравнений
- •3.2.1. Общий вид системы линейных алгебраических уравнений (слау). Понятие о прямых методах решения слау.
- •3.2.2. Метод Гаусса.
- •3.2.3. Вычисление определителя матрицы коэффициентов системы. Некоторые свойства определителей.
- •3.2.4. Выбор главного элемента.
- •3.2.6. Итерационное уточнение.
- •3.2.7. Метод Гаусса для ленточной слау.
- •3.2.8. Слау с несколькими правыми частями.
- •3.2.9. Вычисление обратной матрицы .
- •§ 3.3. Итерационные методы решения систем линейных алгебраических уравнений
- •3.3.1. Понятие об итерационных методах решения слау.
- •3.3.2. Метод Зейделя и метод простой итерации.
- •3.3.3. Общий вид итерационного процесса.
- •§ 3.4. Методы вычисления собственных значений и собственных векторов матриц
- •3.4.1. Основные определения.
- •3.4.2. Характеристическое уравнение.
- •3.4.3. Понятие о методах решения проблемы собственных значений.
- •3.4.4. Определение максимального по модулю собственного числа степенным методом.
- •3.4.5. Собственные значения обратной матрицы.
- •3.4.6. Обобщенная проблема собственных значений.
- •3.4.7. Вычисление минимального собственного числа степенным методом.
- •§ 3.5. Методы численного интегрирования
- •3.5.1. Понятие о формулах численного интегрирования.
- •3.5.2. Формула прямоугольников (формула средних).
- •3.5.3. Формула трапеций.
- •3.5.4. Формула Симпсона.
- •3.5.5. Точность квадратурных формул.
- •3.5.6. Вычисление интеграла с заданной точностью.
- •3.5.7. О современных методах численного интегрирования.
- •§ 3.6. Методы решения нелинейных уравнений
- •3.6.1. Уравнения с одним неизвестным (скалярные).
- •3.6.2. Решение систем нелинейных уравнений.
- •§ 3.7. Метод наименьших квадратов
- •3.7.1. Понятие о методе наименьших квадратов (мнк).
- •3.7.2. Применение мнк для решения переопределенных систем линейных алгебраических уравнений.
3.1.4. Скалярное произведение векторов.
Несложно показать, что арифметическое пространство является евклидовым пространством.
Скалярное произведение векторов
и
обозначается как
,
(3.1.30)
т.е. двум векторам и ставится в соответствие число s (скобки – это стандартное обозначение скалярного произведения).
В соответствии с данным выше определением евклидового пространства скалярное произведение векторов должно обладать свойствами:
;
;
;
,
если
и
,
если
.
(3.1.31)
Чаще всего скалярное произведение векторов определяется формулой
.
(3.1.32)
Можно показать, что оно удовлетворяет всем указанным выше свойствам.
В частности, евклидовым пространством является пространство векторов, рассматриваемых в курсе аналитической геометрии.
Положительно определенной матрицей называется вещественная сим-
метричная матрица, если для любого вектора пространства справедливы соотношения:
,
если
и
,
если
.
(3.1.33)
Заметим, что определенная выше операция скалярного произведения векторов является наиболее употребляемой. Вместе с тем, существует более общая формула определения скалярного произведения векторов
,
(3.1.34)
где
– положительно определенная симметричная
матрица. В выше приведенном случае
(формула (3.1.32))
.
Можно показать, что для любой матрицы A справедливо равенство
.
(3.1.35)
3.1.5. Норма вектора. Некоторые свойства векторных норм. Абсолютная и относительная погрешности. Сходимость.
Некоторым обобщением понятия длины вектора (оценкой величины вектора) является величина, называемая норма вектора. Функция
является нормой вектора, если она обладает следующими свойствами:
1)
0
(равенство только при
);
2)
;
3)
– неравенство треугольника.
Чтобы отличить одну норму от другой используются индексы при двойной черте.
Полезный класс векторных норм – это p-нормы, определяемые как
.
(3.1.36)
Наиболее важными из p-норм
являются 1, 2 и
нормы:
(3.1.37)
– норма
по модулю (частный случай (3.1.36) при
);
(3.1.38)
–
«евклидова» норма (частный случай
(3.1.36) при
);
(3.1.39)
– максимум модуля
для элементов вектора (частный случай
(3.1.36) при
).
Очевидно, что
евклидова норма
.
Это соответствует естественному понятию
длины вектора в двумерном и трехмерном
пространстве.
Единичным вектором по отношению к
норме
называется вектор
,
удовлетворяющий равенству
.
Классический результат о p-нормах – неравенство Гельдера
.
(3.1.40)
Очень важным частным случаем этого неравенства является неравенство Коши-Шварца:
.
(3.1.41)
Все нормы в
эквивалентны, т.е. для двух норм
и
в
существуют положительные константы
и
,
такие, что
(3.1.42)
для всех из пространства . Например, для вектора из имеем:
;
;
.
(3.1.43)
В вычислительной математике большую роль играет понятие «близости» двух векторов. Это соответствует понятию «расстояния», определяемому формулой
,
(3.1.44)
т.е. расстояние между векторами равно норме от их разности:
,
где
.
(3.1.45)
Оценка близости двух векторов используется для определения погрешности решения в итерационных методах.
Пусть вектор
пространства
есть аппроксимация к вектору
пространства
.
Для заданной векторной нормы
будем говорить, что
(3.1.46)
есть абсолютная погрешность
,
а при
формула
(3.1.47)
задает относительную погрешность .
Сходимость. Будем говорить,
что последовательность векторов
,
,
…,
,
… сходится к вектору
,
если
.
(3.1.48)
Отметим, что вследствие приведенных выше свойств векторных норм сходимость в -норме влечет сходимость в -норме и наоборот.