- •«Численные методы»
- •Оглавление
- •§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. Вычислить интеграл от заданной функции f(x) на отрезке [a; b] при делении отрезка на 10 равных частей по формуле прямоугольников.
Решение:
i |
xi |
yi = x2sinx (i =0, 1,…,9, 10) |
0 1 2 3 4 5 6 7 8 9 10 |
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 |
- 0,0009983 0,0079467 0,0265968 0,0623068 0,1198562 0,2032711 0,3156668 0,4591078 0,6344948 0,841471 |
|
|
2,671716 |
= 0,1·2,671716 = 0,2671716
Оценим погрешность: F(x) = x2sinx
F'(x) = 2x sin x + x2cosx
,
= 0,3± 0,1.
Задание 2. Вычислить интеграл от заданной функции f(x) на отрезке [a; b] при делении отрезка на 10 равных частей по формуле трапеций.
Решение:
i |
xi |
yi = x2sinx (i =1,…,9) |
yi / 2 (i = 0;10) |
0 1 2 3 4 5 6 7 8 9 10 |
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 |
0,0009983 0,0079467 0,0265968 0,0623068 0,1198562 0,2032711 0,3156668 0,4591078 0,6344948 |
0
0,4207355 |
|
|
1,8302453 |
0,4207355 |
= 0,1(1,8302453+0,4207355) = 0,225098;
Оценим погрешность: F(x) = x2sinx
F'(x) = 2x sin x + x2cosx, F''(x) = (2– x2) sinx + 4x cosx
,
= 0,23± 0,01.
Задание 3. Вычислить интеграл от заданной функции f(x) на отрезке [a; b] при делении отрезка на 10 равных частей по формуле парабол.
Решение:
i |
xi |
yi / 2 (i = 0;10) |
yi = x2sinx (i =2; 4; 6; 8) |
2yi = 2x2sinx (i =1; 3; 5; 7; 9) |
0 |
0 |
0 |
|
|
1 |
0,1 |
|
|
0,0019966 |
2 |
0,2 |
|
0,0079467 |
|
3 |
0,3 |
|
|
0,0531936 |
4 |
0,4 |
|
0,0623068 |
|
5 |
0,5 |
|
|
0,2397124 |
6 |
0,6 |
|
0,2032711 |
|
7 |
0,7 |
|
|
0,6313333 |
8 |
0,8 |
|
0,4591078 |
|
9 |
0,9 |
|
|
1,2689896 |
10 |
1 |
0,4207355 |
|
|
|
|
0,4207355 |
0,7326324 |
2,1952255 |
= 2∙0,1/3(0,4207355+0,7326324+2,1952255) = 0,2232395
Оценим погрешность: F(4)(x) = (x2– 12) sinx – 8x cosx
= 14;
= 0,22324± 0,00001.