
- •0. Указатель/вопросы.
- •I. Основные понятия.
- •II. Численные методы линейной алгебры.
- •III. Решение нелинейных уравнений и систем нелинейных уравнений.
- •IV. Интерполирование функций.
- •V. Численное дифференцирование и интегрирование функций.
- •VI. Численные методы решения задачи коши для обыкновенных дифференциальных уравнений.
- •VII. Приближение функций и смежные вопросы.
- •Основные понятия.
- •1. Понятие метрики и метрического пространства, примеры.
- •2. Сходящиеся и фундаментальные последовательности, открытые и замкнутые шары, полные метрические пространства.
- •3. Линейные нормированные пространства. Аксиомы нормы. Примеры определения норм в различных пространствах.
- •4. Элементы общей теории погрешностей. Основные определения, утверждения (абсолютная и относительная погрешности, погрешности основных арифметических операций).
- •5. Особенности представления чисел в эвм и компьютерной арифметики. Машинное эпсилон и алгоритм его вычисления.
- •6. Обусловленность вычислительной задачи и вычислительного алгоритма. Абсолютное и относительное число обусловленности.
- •7. Понятие оператора и неподвижной точки оператора. Принцип сжимающих отображений.
- •II. Численные методы линейной алгебры.
- •8. Нормы векторов и матриц. Относительная и абсолютная погрешность вектора. Связь погрешности решения систем линейных алгебраических уравнений и невязки.
- •9. Характеристика и отличительные особенности прямых и итерационных методов численного решения систем линейных алгебраических уравнений. Примеры.
- •11. Алгоритм lu разложения матрицы.
- •12. Решение систем линейных алгебраических уравнений с помощью lu разложения матрицы системы.
- •13. Алгоритм qr разложения матрицы.
- •15. Использование lu и qr разложения матрицы системы для решения систем с одинаковой матрицей и различными правыми частями (на примере вычисления обратных матриц).
- •16. Вычисление определителя матрицы на основе ее lu или qr разложения.
- •17. Метод прогонки решения систем линейных алгебраических уравнений с матрицей специального вида.
- •18. Итерационный метод Якоби решения слау, оценка погрешности полученного решения. Геометрическая иллюстрация метода для двух уравнений.
- •19. Итерационный метод Зейделя решения слау, оценка погрешности полученного решения. Геометрическая иллюстрация метода для двух уравнений.
- •III. Решение нелинейных уравнений и систем нелинейных уравнений.
- •IV. Интерполирование функций.
- •32. Задача интерполяции. Обобщенный многочлен по системе функций. Глобальная и локальная интерполяции. Преимущества и недостатки. Примеры.
- •33. Интерполяция алгебраическими многочленами. Системы функций Чебышёва. Существование и единственность интерполяционного алгебраического многочлена.
- •34. Интерполяционный многочлен Лагранжа, особенности применения.
- •35. Вычислительная схема Эйткена.
- •36. Оценка погрешности интерполяции многочленом Лагранжа.
- •37. Разделенные разности, основные свойства. Интерполяционный многочлен Ньютона, особенности применения, оценка погрешности.
- •40. Определение сплайна. Гладкая кусочно-полиномиальная интерполяция. Определение интерполяционного сплайна.
- •42. Глобальные и локальные кубические интерполяционные сплайны, Определения и способы построения.
II. Численные методы линейной алгебры.
8. Нормы векторов и матриц. Относительная и абсолютная погрешность вектора. Связь погрешности решения систем линейных алгебраических уравнений и невязки.
Норма — число, отличное от нуля. Предназначена для сравнения векторов, матриц, и т. п. (какой вектор/матрица «больше», какой «меньше»).
Нормы векторов.
Норма вектора — число, характеризующее величину многомерного объекта.
Норма вектора — число, обозначаемое и удовлетворяющее следующим условиям — аксиомам нормы (обобщающим свойства длины 3-мерного в-ра):
1. , причем .
2. , — действительное число.
3. — неравенство треугольника.
Часто норму
обозначают в виде
.
Нормированное пространство — векторное пространство с нормой, аксиомы нормы в таком случае — аксиомы нормированного пространства.
Нормированный
вектор (нормальный) — с единичной
нормой (
).
Любой ненулевой
вектор можно нормировать, т. е.
разделить его на свою норму: вектор
имеет единичную норму. С геометрической
точки зрения это значит, что мы берем
сонаправленный вектор единичной длины.
Векторные нормы в различных пространствах: см. вопрос 3.
Абсолютная
погрешность вектора:
.
Относительная
погрешность вектора:
.
Нормы матриц.
Норма матрицы
A — вещественное число
,
удовлетворяющее условиям:
1.
,
причем
.
2.
,
— действительное число.
3.
— неравенство треугольника.
Виды матричных норм:
Если выполняется 4-ое условие:
, норма называется мультипликативной.
Операторная норма, порожденная векторной нормой — число
.
Матричная норма, составленная как операторная, называется подчинённой по отношению к норме, использованной в пространствах векторов.
m-норма:
.
l-норма:
.
Норма Фробениуса:
Невязка —
вектор
.
— система уравнений.
Связь между
погрешностью решения СЛАУ и величиной
невязки определяется числом
обусловленности матрицы —
характеристикой, не связанной с
каким-либо численным алгоритмом и
отражающей только свойства системы.
Это число обозначают Cond(A) и определяют
по формуле:
.
9. Характеристика и отличительные особенности прямых и итерационных методов численного решения систем линейных алгебраических уравнений. Примеры.
Численное решение уравнений и их систем — приближённое определение корней уравнения или системы уравнений. Применяется, когда точное значение вычислить невозможно или очень трудоёмко.
Методы численного решения СЛАУ делятся на прямые и итерационные.
Прямые методы — позволяющие получить решение системы за конечное число арифметических операций. Метод Крамера, метод Гаусса, LU-метод и т.д.
Итерационные
методы (методы последовательных
приближений) — решение системы находится
как предел последовательных приближений
при номере итерации
.
При использовании методов итерации
обычно задается некоторое малое число
и вычисления проводятся до тех пор,
пока не будет выполнена оценка
.
Метод Зейделя, Якоби, метод верхних
релаксаций и т.д.
Реализация прямых методов на компьютере приводит к решению с погрешностью, т. к. все арифметические операции над переменными с плавающей точкой выполняются с округлением. В зависимости от свойств матрицы исходной системы эти погрешности могут достигать значительных величин.
1. Метод Жордана-Гаусса (метод последовательного исключения неизвестных).
Алгоритм. Система в развернутом виде:
Метод Гаусса
состоит в последовательном исключении
неизвестных из этой системы. Предположим,
что
.
исключается из всех уравнений, кроме
первого, последовательно умножая первое
уравнение на
и складывая с i-м уравнением:
Прямой ход Гаусса. Аналогичным образом исключаются остальные неизвестные и получается система треугольного вида:
Можно получить все значения неизвестных, выполняя последовательные подстановки в последней системе (начиная с последнего уравнения). Обратный ход Гаусса.
Характеристики.
Метод легко реализуем на компьютере.
При выполнении вычислений, как правило,
не интересуют промежуточные значения
матрицы
.
Поэтому численная реализация метода
сводится к преобразованию элементов
массива размерности
,
где
столбец содержит элементы правой части
системы.
Один из основных недостатков: накапливается вычислительная погрешность.
2. LU-метод — см. вопросы 11-12.
3. QR-метод — см. вопросы 13-14.