
- •«Численные методы»
- •Оглавление
- •§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
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •Литература
6.2 Методы прямоугольников и трапеций
Если функция f (x) на отрезке [xi–1; xi] заменяется полиномом нулевой степени, то приближенное значение интеграла на частичном отрезке будет иметь вид
, где
hi=xi-xi+1
, а fi=f(i)
- значение функции в точке i=(xi+xi-1)/2.
Данная формула означает, что определенный интеграл от функции f(x) на элементарном отрезке [xi-1,xi] приближенно равен площади прямоугольника, у которого основание равно hi, а высота - значению функции, вычисленное в средней точке отрезка.
Рисунок 6.2. Геометрическая интерпритация метода прямоугольников
В зависимости от выбора точки ξ i получаются различные формулы прямоугольников, то есть ξ i = х i–1 (рис.6.2 а) или ξ i = х i (рис. 7.2 б):
(6.1)
или
(6.2)
Формулы (6.1) и (6.2) называются формулами прямоугольников.
Формула
для оценки погрешности:
(6.3),
где
Заменяя в частичном
интеграле
функцию
f (x)
линейным полиномом
получаем
формулу трапеций на частичном отрезке
(рис. 7.3)
Рисунок 6.3 Метод трапеций
Общая формула трапеций получается суммированием частичных интегралов и имеет вид:
(6.4)
формула
для оценки погрешности:
(6.5),
где
6.3 Метод Симпсона (метод парабол)
Аппроксимируя в частичном интеграле функцию f (x) квадратичным полиномом Лагранжа получаем так называемую формулу Симпсона (формулу парабол) для частичного интервала.
Рисунок 6.6. Метод парабол
Для всего отрезка [a, b] соответственно получаем:
(6.6)
Формула
для оценки погрешности:
(6.7),
где
Задания для самостоятельного решения
Приближенно вычислить интеграл с использованием метода трапеций, при этом на всем отрезке интегрирования использовать пять узлов и в результате оставить только верные знаки.
.
Приближенно вычислить интеграл с использованием метода Симпсона, при этом на всем отрезке интегрирования использовать пять узлов и в результате оставить только верные знаки.
.
Приближенно вычислить интеграл с использованием метода прямоугольников, при этом на всем отрезке интегрирования использовать пять узлов и в результате оставить только верные знаки. Оценить погрешность результата
.
Приближенно вычислить интеграл с использованием метода трапеций, при этом на всем отрезке интегрирования использовать пять узлов и в результате оставить только верные знаки. Оценить погрешность результата
.
Приближенно вычислить интеграл с использованием метода Симпсона, при этом на всем отрезке интегрирования использовать пять узлов и в результате оставить только верные знаки. Оценить погрешность результата
.
Практическая работа №7
Тема: «Вычисление интегралов по формулам Ньютона-Котеса»
Цели: освоение вычисления интегралов приближенными методами прямоугольников, трапеций и с помощью формулы Симпсона - методом парабол; сравнение методов.
Задание 1. Вычислить интеграл от заданной функции f(x) на отрезке [a; b] при делении отрезка на 10 равных частей по формуле прямоугольников.
Задание 2. Вычислить интеграл от заданной функции f(x) на отрезке [a; b] при делении отрезка на 10 равных частей по формуле трапеций.
Задание 3. Вычислить интеграл от заданной функции f(x) на отрезке [a; b] при делении отрезка на 10 равных частей по формуле парабол.
Задание 4. Сравнить полученные результаты.
Исходные данные:
Вариант 1.
;
Вариант 2.
Вариант
3.
;
Вариант 4.
Вариант
5.
;
Вариант 6.
Вариант 7.
;
Вариант 8.
Вариант 9.
;
Вариант 10.
Вариант
11.
; Вариант 12.
Вариант 13.
;
Вариант 14.
Вариант 15.
;
Вариант 16.
Вариант 17.
;
Вариант 18.
Вариант 19.
;
Вариант 20.