- •«Численные методы»
- •Оглавление
- •§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. Дана система трех линейных уравнений с тремя неизвестными:
Решите систему методом Гаусса, используя схему единственного деления; расчеты выполняйте с тремя знаками после запятой; подставьте найденные решения в исходную систему, вычислите невязки и сделайте вывод о полученном решении.
Решение:
Расчеты будем производить с использованием расчетной таблицы:
Таблица 4.1
раздел |
х1 |
х2 |
х3 |
свободные члены |
S1 |
S2 |
А |
2,34 8,04 3,92 |
– 4,21 5,22 – 7,49 |
–11,61 0,27 8,37 |
14,41 6,44 55,56 |
0,93 7,09 59,86 |
|
|
1 |
– 1,799 |
– 4,962 |
6,158 |
0,397 |
0,397 |
А1 |
|
19,685 – 0,938 |
40,161 27,819 |
– 55,951 31,420 |
3,895 58,302 |
3,894 58,302 |
|
|
1 |
2,040 |
– 2,842 |
0,198 |
0,198 |
А2 |
|
|
29,732 |
28,756 |
58,488 |
58,487 |
|
|
|
1 |
0,967 |
1,967 |
|
В |
1 |
1 |
1 |
0,967 – 4,816 2,293 |
1,967 – 3,816 3,293 |
|
Раздел А – из заданной системы коэффициенты и свободные члены.
Для исключения случайных ошибок предусматривается контроль правильности вычислений, с этой целью в таблицу включаются столбец контрольных сумм S1 и столбец строчных сумм S2. Числа в этих столбцах должны практически совпадать.
Значения в разделе А в столбце S1 получаются как сумма чисел в строке.
Строка под разделом А получается путем деления 1-й строчки раздела А на соответствующий коэффициент х1, т.е. на 2,34.
Значение в столбике S1 получается путем тех же преобразований, что и в строчке (аналогично весь столбик S1). Значения в столбике S2 получаются как сумма чисел в строке левой части таблицы.
Раздел А1. Значения первой строчки получаются путем преобразования второй строки раздела А: 19,685 = 5,22– 8,04∙ (– 1,799); 40,161= 0,27– 8,04∙ (– 4,962); – 55,951= 6,44– 8,04∙ 6,158. Значения второй строчки получаются путем аналогичных преобразований третьей строки раздела А
Раздел А2 – аналогично.
Раздел В – обратный ход. Последняя строка прямого хода переносится, т.е.
х3 = 0,967. Вторая строка раздела В получается из строчки под разделом А1:
строчка под разделом А1 записывается в виде уравнения х2 + 2,040х3 – 2,842 = 0, в него подставляется значение х3 и выражается х2. Аналогично находится х1. Значения в столбике S1 должны отличаться от значений неизвестных на 1
х1 = 2,293; х2 = – 4,816; х3 = 0,967.
Если подставить найденные значения в систему:
тогда разности между правой и левой частями, называемые невязками, будут следующими:
ε1 = 14,41– 14,411= – 0,001;
ε2 = – 6,44 – (– 6,441) = 0,001;
ε3 = 55,56 – 55,561 = 0,001.