- •«Численные методы»
- •Оглавление
- •§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. Локализуйте корни заданного уравнения, на заданном отрезке. Укажите количество корней и отрезки, в каждом из которых заключен корень уравнения.
Для уравнения sin 2x – ln x = 0 имеется единственный отрезок, содержащий корень данного уравнения – [1,3; 1,4].
Задание 2. Найдите один из кореней заданного уравнения, с погрешностью ε = 0,001 методом половинного деления;
sin 2x – ln x = 0, [1,3; 1,4]
Решение:
F(x) = sin 2x – ln x
а = 1,3; b = 1,4
х0 = (1,3 + 1,4)/2 = 1,35; ∆x0 = |1,3 – 1,4|/2 = 0,05
получаем два отрезка: [1,3; 1,35] , [1,35; 1,4]
F(1,3) = sin (2∙1,3) – ln 1,3 = 0,25; F(1,35) = 0,13; F(1,4) = - 0,001
на отрезке [1,35; 1,4] функция меняет знак
х1= (1,35 + 1,4)/2 = 1,375; ∆ = |1,35 – 1,4|/2 = 0,025
[1,35; 1,375] , [1,375; 1,4]
F(1,35) = 0,13; F(1,375) = 0,06; F(1,4) = - 0,001
[1,375; 1,4]
х2=(1,375 + 1,4)/2 = 1,3875; ∆ = |1,375 – 1,4|/2 = 0,0125
[1,375; 1,3875] , [1,3875; 1,4]
F(1,375) = 0,06; F(1,3875) = 0,03; F(1,4) = - 0,001
[1,3875; 1,4]
х3= (1,3875 + 1,4)/2 = 1,39375; ∆ = |1,3875 – 1,4|/2 = 0,00625 0,01
Ответ: х = 1,39 ± 0,01
Задание 3. Найдите один из кореней заданного уравнения, с погрешностью ε = 0,001 методом простой итерации;
sin 2x – ln x = 0, [1,3; 1,4]
Решение:
F(x) = sin 2x – ln x
1. Приводим уравнение к виду x = f(x)
x = x – μ(sin 2x – ln x)
2. Подберем константу μ, μ = –1/׀М׀
F'(x) = 2cos 2x – 1/ x
F'(1,3) = – 2,48; F'(1,4) = – 2,6
׀М׀ = 2,6, μ = – 1/2,6 0,38
3. Получаем итерационную формулу:
xn+1= xn + 0,38(sin 2xn – ln xn)
4. Найдем q
5. Выберем х0, например х0 = 1,35 (середина отрезка)
Итерационный процесс удобно оформлять в виде таблицы:
n |
xn |
xn+1= xn + 0,38(sin 2xn – ln xn) |
|
0 |
1,35 |
1,3983646 |
0,00059 |
Как видим, уже значение х1 = 1,3983646 соответствует заданной точности 0,001, округляя результат до верных цифр, получаем
Ответ: х = 1,398 ± 0,001
Контрольные вопросы.
1. Из каких этапов состоит решение нелинейного уравнения?
2. В чем заключается задача отделения корней?
3. В чем состоит основная идея метода половинного деления?
4. Каков алгоритм решения уравнений методом простой итерации?
5. В чем недостатки и достоинства методов половинного деления и простой итерации?
3.5 Методы Ньютона
Метод касательных
Пусть корень локализован на отрезке [a,b].
В
ыбирается
точка x0
на этом отрезке и в этой точке строится
касательная к графику функции y=f(x)
(рис.3.3). За
новое приближение x1
принимается точка, в которой касательная
пересекает ось 0x.
Далее процесс повторяется. Таким образом,
на каждом этапе решение исходного
нелинейного уравнения заменяется на
решение линейного уравнения (уравнения
касательной), т.е. проводится линеаризация
исходного уравнения.
Формула метода
касательных
имеет вид:
(3.3)
Оценка
погрешности производится по формуле
(3.4),
где
,
ε – заданная погрешность.
Метод хорд (метод секущих)
Пусть задан отрезок [a,b], на концах которого функция f (x) имеет значения разных знаков – f (a) f (b) < 0 .
Алгоритм метода хорд:
1. с – конец отрезка [a; b] для которого выполняется F(x)·F''(x) >0, противоположный конец принимают за х0.
2.
Итерационная формула:
(3.5),
3.
Оценка погрешности производится по
формуле:
(3.6),
где
,
ε – заданная погрешность.
Графическая интерпритация метода хорд представлена на рисунке 3.4.
Рисунок 3.5. Приближения к корню уравнения по методу хорд
