- •Численные методы
- •Содержание
- •Введение
- •Тема №1. Решение уравнений
- •1. Основные сведения
- •1.1 Отделение корней уравнения
- •Решение
- •1.2 Нахождение корней уравнения
- •Решение
- •Метод касательных (метод Ньютона)
- •Решение
- •2. Лабораторная работа №1 «Решение уравнений».
- •3. Контрольные вопросы
- •Тема №2. Решение систем линейных уравнений
- •1. Основные сведения
- •1.1 Прямые методы
- •Решение
- •Решение
- •1.2 Итерационные методы Метод простой итерации
- •Решение
- •Метод Зейделя
- •2. Лабораторная работа №2 «Решение систем линейных уравнений»
- •3. Контрольные вопросы
- •Тема №3. Интерполяция и аппроксимация функций
- •Основные сведения
- •Интерполяция
- •1.1.1. Глобальная полиномиальная интерполяция
- •Интерполяционный полином Лагранжа
- •Решение
- •Интерполяционный многочлен Ньютона
- •Решение
- •Решение
- •Локальная (кусочно-полиномиальная) интерполяция
- •Решение
- •Решение
- •Интерполяция сплайнами
- •1.2. Аппроксимация. Метод наименьших квадратов
- •Решение
- •2. Лабораторная работа № 3 «Интерполяция и аппроксимация функций»
- •3. Контрольные вопросы
- •3.1. Интерполяция
- •Аппроксимация
- •Тема №4. Дифференцирование функций
- •Основные сведения
- •Вычисление первой производной
- •Вычисление второй производной
- •Дифференцирование с использованием интерполяционных формул Дифференцирование на основе интерполяционной формулы Ньютона
- •Решение
- •Дифференцирование на основе многочлена Лагранжа
- •Лабораторная работа №4 «Дифференцирование функций».
- •Контрольные вопросы.
- •Тема №5 Интегрирование функций
- •Основные сведения
- •1.1. Общие принципы численного интегрирования.
- •Простейшие квадратурные формулы прямоугольников, трапеций и Симпсона.
- •Формулы прямоугольников
- •Формула трапеций
- •Формула Симпсона (формула парабол)
- •Решение
- •1.3. Квадратурные формулы Ньютона – Котеса.
- •1.4. Квадратурные формулы Гаусса.
- •Решение
- •2. Лабораторная работа №5 «Интегрирование функций».
- •3. Контрольные вопросы.
- •Тема №6 Решение обыкновенных дифференциальных уравнений
- •1. Основные сведения.
- •1.1. Понятие о численном решении задачи Коши для дифференциального уравнения первого порядка.
- •1.2. Метод Эйлера. Модифицированный метод Эйлера.
- •Решение
- •1.3. Методы Рунге - Кутта.
- •Решение
- •Лабораторная работа №6 «Решение обыкновенных дифференциальных уравнений».
- •3. Контрольные вопросы. Метод Эйлера
- •Решение
- •Решение
2. Лабораторная работа №1 «Решение уравнений».
Задания
Найти все корни уравнения с точностью до 0,0001.
Примечание. Для расчетов можно использовать программы №№1 - 2 (на языке QBasic) из Приложения 3.
1. |
2.
|
3.
|
|
4.
|
5.
|
6.
|
|
7.
|
8.
|
9.
|
|
10.
|
11.
|
12.
|
|
13.
|
14.
|
15.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
3. Контрольные вопросы
Отделение корней уравнения
Что дает отделение корней?
Можно ли произвольно задавать значения по оси
для отделения корней?Какие основные проблемы могут встретиться при аналитическом отделении корней?
Нахождение корней уравнения
Метод бисекции (половинного деления)
В чем заключается геометрический смысл метода бисекции (половинного деления)?
Всегда ли позволяет метод половинного деления вычислить отделенный корень уравнения с заданной погрешностью?
Как выбираются концы отрезка следующего интервала в методе половинного деления?
Какими свойствами должна обладать функция
,
чтобы методом половинного деления
можно было гарантированно решить
уравнение
?Можно ли найти корень методом половинного деления, если он находится на границе интервала?
Метод простой итерации
Какой функцией заменяется левая часть уравнения
в
методе итераций?Что называется сходимостью метода итераций?
С какой стороны может осуществляться приближение к корню в методе итераций — слева или справа?
Если на заданном отрезке имеется два корня, то, что можно сказать о сходимости метода итераций на этом отрезке?
Что означает несходимость процесса итераций?
Метод касательных (метод Ньютона)
В чем заключается геометрическая интерпретация метода Ньютона?
Из каких соображений в методе Ньютона выбирается первое приближение
?Что необходимо для того, чтобы уравнение
решалось
методом Ньютона?В каких случаях применение метода Ньютона не рекомендуется?
Тема №2. Решение систем линейных уравнений
1. Основные сведения
Рассмотрим решение системы n уравнений с n неизвестными:
(2.1)
Система такого вида называется линейной.
Методы решения системы линейных уравнений разделяются на прямые и итерационные.
1.1 Прямые методы
Одним из наиболее распространенных прямых методов является метод Гаусса и его модификации: схема частичного выбора, схема полного выбора.
Простейший вариант метода Гаусса называется схемой единственного деления, и включает два этапа: прямой ход и обратный ход.
Прямой
ход
метода состоит из
шагов исключения, в результате чего
система приводится к виду:
(2.2)
Это преобразование системы осуществляется следующим образом:
Первый шаг
Предположим, что элемент
.
Будем называть этот элемент главным
или ведущим
элементом первого шага.Найдем величины
.Вычтем последовательно из второго, третьего,…, n-ого уравнения системы (2.1) первое уравнение, умноженное соответственно на
.
В результате коэффициенты при x1 во всех уравнениях, кроме первого, обратятся в нуль.
Второй шаг
На втором шаге при помощи a22 (главного или ведущего элемента второго шага) также обратим в нуль коэффициенты при x2 во всех уравнениях, начиная с третьего.
И так далее.
После
таких действий получим систему (2.2),
матрица которой является верхней
треугольной. На этом заканчивается
прямой ход метода.
Обратный
ход
метода начинается с нахождения
из последнего уравнения системы (2.2).
Подставляя найденное
в предпоследнее уравнение, находим
.
И таким образом, продвигаясь от последних
уравнений системы к первому, находим
все значения неизвестных.
Пример №1. Решить систему методом Гаусса
