- •«Численные методы»
- •Оглавление
- •§1. Теоретические основы численных методов 10
- •§2. Особенности математических вычислений на эвм. Погрешности вычислений 13
- •§3. Численные методы решения алгебраических и трансцендентных уравнений 25
- •§4. Методы решения систем уравнений 38
- •Введение
- •Из истории вычислительной математики
- •§1. Теоретические основы численных методов
- •§2. Особенности математических вычислений на эвм. Погрешности вычислений
- •Задания для самостоятельного решения
- •Практическая работа №1
- •Примеры выполнения заданий работы
- •Контрольные вопросы.
- •§3. Численные методы решения алгебраических и трансцендентных уравнений
- •3.1 Задача решения алгебраических и трансцендентных уравнений
- •3.2 Локализация корней
- •3.3 Метод деления отрезка пополам (метод бисекции, метод дихотомии)
- •3.4 Метод простой итерации
- •Задания для самостоятельного решения
- •Практическая работа №2
- •Примеры выполнения заданий работы
- •Контрольные вопросы.
- •3.5 Методы Ньютона
- •3.6. Решение уравнений с помощью табличного процессора Excel
- •Задания для самостоятельного решения
- •Практическая работа №3
- •Примеры выполнения заданий работы
- •Контрольные вопросы.
- •§4. Методы решения систем уравнений
- •4.1 Система линейных уравнений
- •4.1.1 Прямые методы решения систем линейных уравнений
- •Задания для самостоятельного решения
- •Практическая работа №4
- •Примеры выполнения заданий работы
- •4.1.2 Вычисление определителей и обратной матрицы
- •Задания для самостоятельного решения
- •Контрольные вопросы.
- •4.1.3 Итерационные методы решения систем линейных уравнений. Метод простой итерации
- •4.2. Решение системы уравнений и вычисление определителя с помощью табличного процессора Excel
- •Задания для самостоятельного решения
- •Практическая работа №5
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •§5. Методы приближения и аппроксимации функций
- •5.1 Понятия интерполяции и экстраполяции
- •5.2 Интерполяционный многочлен Лагранжа
- •Задания для самостоятельного решения
- •5.3 Приближение функций с помощью табличного процессора Excel
- •Задания для самостоятельного решения
- •Практическая работа №6
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •5.3 Интерполяционные формулы Ньютона
- •Задания для самостоятельного решения
- •§6. Численное интегрирование
- •6.1 Задача численного интегрирования
- •6.2 Методы прямоугольников и трапеций
- •6.3 Метод Симпсона (метод парабол)
- •Задания для самостоятельного решения
- •Практическая работа №7
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •6.4 Квадратурная формула Гаусса
- •6.5. Вычисление интеграла с использованием табличного процессора Excel.
- •Задания для самостоятельного решения
- •Практическая работа №8
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •§7. Решение обыкновенных дифференциальных уравнений
- •7.1. Задача численного решения обыкновенных дифференциальных уравнений
- •7.2. Методы Эйлера
- •7.3 Метод Рунге – Кутта
- •Задания для самостоятельного решения
- •Практическая работа №9
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •§8. Методы оптимизации
- •8.1 Методы одномерной оптимизации
- •Задания для самостоятельного решения
- •8.2 Методы многомерной оптимизации
- •8.3. Решение задач оптимизации с помощью табличного процессора Excel
- •Задания для самостоятельного решения
- •Практическая работа №10
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •Литература
Примеры выполнения заданий работы
Задание 1. Дана система трех линейных уравнений с тремя неизвестными:
решите систему линейных уравнений методом простой итерации.
Решение:
Исходную систему приводим к системе с преобладающими диагональными коэффициентами и к виду Х = АХ.
Переставим уравнения в системе по преобладанию коэффициентов: в первом уравнении преобладает коэффициент перед х3, во втором – перед х1, сложив третье уравнение с первым, мы получим уравнение с преобладанием второго коэффициента:
Разделим каждое уравнение на диагональный коэффициент и выразим диагональное неизвестное:
(итерационная формула)
Проверяем условия сходимости:
условие выполняется
α = 0,78
Заполняем итерационную таблицу:
№ итерации |
х(к) |
х(к + 1) |
|
0 |
– 0,801 – 5,735 – 1,241 |
2,963 – 5,8158 0,679 |
1 3,35 0,29 13,35 – 6,8 1 |
1 |
2,963 – 5,8158 0,679 |
2,9504 – 4,3956 1,46866 |
0 ,05 5,04 5,04 2,7 998 |
2 |
2,9504 – 4,3956 1,46866 |
2,0018 – 4,6121 0,9506 |
3 ,37 0,77 3,37 1,84 |
На второй итерации с погрешностью 3,37 получили результат: х1= 2,0018, х2 = – 4,6121, х3 = 0,9506. Продолжая итерационный процесс, можно получить результат с меньшей погрешностью.
Задание 2. Дана система трех линейных уравнений с тремя неизвестными:
решите систему линейных уравнений методом Зейделя.
Решение:
Все вычисления до итерационного процесса те же что и в методе простой итерации.(см. задание 1).
Расчеты в нулевой итерации производятся следующим образом:
№ итерации |
х(к) |
х(к + 1) |
|
0 |
– 0,801 – 5,735 – 1,241 |
2,963 – 3,8848 0, 7678 |
13,33
|
1 |
2,963 – 3,8848 0, 7678 |
1,694 – 5,0702 1,198 |
4,5 |
2 |
1,694 – 5,0702 1,198 |
2,449 – 4,797 0,995 |
2,7 |
На второй итерации с погрешностью 2,7 получили результат:
х1= 2,449, х2 = – 4,797, х3 = 0,995.
Контрольные вопросы
1. Как привести исходную систему линейных уравнений к системе с преобладающими диагональными коэффициентами?
2. В чем отличие метода Зейделя от метода простой итерации?
3. В чем заключается решение системы линейных уравнений методом Ньютона?
4. Какие методы положены в основу решения системы линейных уравнений в математических программных средствах (табличный процессор Excel, MATCAD, MATLAB)?
§5. Методы приближения и аппроксимации функций
5.1 Понятия интерполяции и экстраполяции
Замена одной функции на другую более простую и удобную называется аппроксимацией.
Интерполирование – аппроксимация функции с помощью алгебраического многочлена n(x) степени n, значения которого в заданных узлах xi (i=0,1,…,n) равны значениям функции в этих же узлах n(xi)=yi=f(xi), при этом полагается, что среди значений xi нет одинаковых, т.е. xixk при ik. Многочлен n(x) называется интерполяционным. Геометрическая интерпретация интерполирования представлена на рис. 5.1.
Чаще всего задача интерполяции функции одной переменной состоит в замене дискретной зависимости y(xi), т.е. N пар чисел (xi,yi), или, по-другому, узлов, некоторой непрерывной функцией (x). При этом основным условием является то, что функция y(x) должна проходить через точки (xi,yi), т. е. y(xi)=yi ,i=1...N, а также возможность вычислить значение (x) в любой точке, находящейся между узлов.
Когда искомое значение (x) вычисляется в точке x, которая находится между каких-либо из узлов xi, говорят об интерполяции, а когда точка x лежит вне границ интервала, включающего все xi - об экстраполяции функции (x).