
- •Основы вычислительной математики Учебно-методическое пособие
- •Оглавление
- •Предисловие
- •I. Введение
- •Контрольные вопросы
- •II. Системы линейных алгебраических уравнений
- •II.1. Основные понятия и определения
- •II.2. Прямые методы решения систем линейных уравнений
- •II.3. Метод Гаусса решения слау
- •II.4. Решение слау с помощью -разложения
- •II.5. Метод итераций решения слау
- •II.6. Метод Зейделя решения слау
- •II.7. Сходимость итерационных процессов для систем линейных уравнений
- •III. Методы решения нелинейных скалярных уравнений
- •III.1. Локализация корней
- •III.2. Метод простых итераций
- •III.3. Приведение нелинейного уравнения к виду , допускающему сходящиеся итерации
- •III.4. Метод Ньютона (Метод касательных)
- •III.5. Контрольные вопросы
- •III.6. Задачи для самостоятельного решения
- •IV. Приближенное решение систем нелинейных уравнений
- •IV.1. Метод Ньютона
- •IV.2. Метод итераций
- •IV.3. Контрольные вопросы
- •IV.4. Задачи для самостоятельного решения
- •V. Интерполяция и аппроксимация
- •V.1. Интерполяция
- •V.2. Интерполяционный многочлен в форме Лагранжа
- •V.3. Интерполяционный многочлен в форме Ньютона
- •V.4. Приближение методом наименьших квадратов
- •V.5. Контрольные вопросы
- •V.6. Задания для самостоятельного решения
- •VI. Равномерное приближение функций
- •V.1. Полиномы Чебышева
- •VI.2. Понятие о равномерном приближении функций
- •VI.3. Контрольные вопросы
- •VII. Приближенное дифференцирование
- •VII.1. Постановка вопроса
- •VII.2. Формулы приближенного дифференцирования, основанные на первой интерполяционной формуле Ньютона
- •VII.3. Формулы приближенного дифференцирования, основные на формуле Стирлинга
- •VII.4. Контрольные вопросы
- •VII.5. Задачи для самостоятельного решения
- •VIII. Численное интегрирование функций
- •VIII.1. Общие замечания
- •VIII.2.Квадратурные формулы Ньютона-Котеса
- •VIII.3. Формула трапеций и ее остаточный член
- •Общая формула трапеций (правило трапеций).
- •VIII.4. Формула Симпсона и ее остаточный член
- •Общая формула Симпсона (параболическая формула).
- •VIII.5. Формулы Ньютона-Котеса высших порядков
- •VIII.6. Контрольные вопросы
- •VIII.7. Задачи для самостоятельного решения
- •IX. Приближенное решение обыкновенных дифференциальных уравнений
- •IX.1. Интегрирование дифференциальных уравнений с помощью степенных рядов
- •IX.2. Метод последовательных приближений
- •IX.3. Метод Эйлера
- •IX.4. Метод Эйлера с итерациями
- •IX.5. Метод Рунге-Кутта
- •IX.6. Контрольные вопросы
- •IX.7. Задачи для самостоятельного решения
- •X. Преобразование фурье
- •X.1. Разложение периодических функций в ряд Фурье
- •X.2. Эффект Гиббса
- •X.3. Контрольные вопросы
- •IX.4. Задачи для самостоятельного решения
- •Список литературы
Контрольные вопросы
Виды погрешностей.
Принцип А.Н. Крылова.
Правила действий над приближенными числами.
Примеры неустойчивых задач.
Пример неустойчивости метода.
II. Системы линейных алгебраических уравнений
II.1. Основные понятия и определения
Приведем некоторые сведения из линейной алгебры, которые потребуются в дальнейшем. Рассмотрим прямоугольную матрицу
.
Две
матрицы
и
размерности
равны друг другу, если
для всех
и
.
Сумма
двух матриц
и
размерности
есть матрица размерности
:
.
Произведение
матрицы
на скаляр
есть матрица размерности
:
.
Введем
обозначение
.
Произведение матрицы
размерности
на матрицу
размерности
есть матрица
размерности
:
.
Таким
образом, элемент
матрицы
есть сумма произведений
-й
строки матрицы
на соответствующие элементы
-го
столбца матрицы
,
при этом число столбцов матрицы
должно равняться числу строк матрицы
.
Из существования произведения
вовсе не следует существование
произведения
.
Для квадратных матриц (
)
одного порядка существуют матрицы
и
,
но, вообще говоря,
.
Определитель
(детерминант) квадратной матрицы будем
обозначать
(или
):
.
Определитель
квадратной матрицы с
элементами (действительными или
комплексными числами
)
есть сумма
членов вида
:
,
где
– алгебраическое дополнение элемента
.
Важным частным случаем квадратной матрицы является диагональная матрица
.
При
матрица называется единичной
и обозначается через
.
Матрица называется нижней треугольной, если все ее элементы, расположенные выше главной диагонали, равны нулю:
.
Аналогично
определяется верхняя
треугольная
матрица. Квадратная матрица называется
симметричной,
если ее элементы удовлетворяют соотношению
(
,
).
Если в матрице поменять строки со
столбцами, то получим транспонированную
матрицу
Квадратная
матрица
симметрична, если
.
Матрица
называется обратной
матрице
,
если
.
Необходимым и достаточным условием
существования обратной матрицы
является условие
.
Говорят, что в этом случае матрица
несобственная,
или невырожденная.
Минором
порядка
матрицы
называется определитель
-го
порядка, составленный из элементов,
которые находятся на пересечении
строк и
столбцов матрицы
.
Рангом
матрицы
называется число
такое, что все миноры порядка
и выше равны нулю.
Действительная
матрица называется ортогональной,
если ее транспонированная матрица
совпадает с обратной, т. е.
или
.
Все собственные значения ортогональной
матрицы по модулю равны единице. Строки
(столбцы) ортогональной матрицы попарно
ортогональны, суммы квадратов элементов
каждой строки (столбца) ортогональной
матрицы равны единице, определитель
ортогональной матрицы равен
;
если матрица
ортогональна, то и матрица
тоже ортогональна.
Всякая симметричная действительная матрица может быть приведена к диагональному виду подобным преобразованием
,
где
– ортогональная матрица,
– диагональная матрица. Из свойств
ортогональной матрицы следует, что
.