- •Введение
- •Глава 1. Погрешность результата численного решения задачи
- •1.1. Источники и классификация погрешностей
- •1.2. Абсолютная и относительная погрешности. Формы записи данных.
- •1.3. Вычислительная погрешность
- •Глава 2. Решение нелинейных уравнений
- •2.1. Отделение корней уравнения
- •2.1.1. Аналитический метод отделения корней
- •2.1.2. Графический метод отделения корней
- •2.2. Уточнение приближенных корней
- •2.2.1. Метод половинного деления
- •2.2.2 Метод хорд
- •2.2.3. Метод Ньютона – метод касательных
- •2.2.4. Модифицированный метод Ньютона
- •2.2.5. Метод секущих
- •2.2.6. Метод итераций
- •Глава 3. Решения систем линейных алгебраических уравнений
- •3.1. Метод итераций
- •3.1.1. Оценка погрешности приближений процесса итераций
- •3.1.2. Приведение линейной системы к виду, удобному для итерации:
- •3.2. Метод Зейделя
- •3.3. Метод релаксаций
- •Глава 4. Решение систем нелинейных уравнений
- •4.1. Метод Ньютона для систем нелинейных уравнений
- •4.2. Распространение метода Ньютона на системы из n уравнений с n неизвестными
- •4.3. Метод итераций для систем нелинейных уравнений
- •4.4. Распространение метода итераций на системы из n уравнений с n неизвестными
- •Глава 5. Интерполяция
- •5.1. Постановка задачи интерполирования
- •5.2. Конечные разности
- •5.3. Интерполяционная формула Ньютона №1
- •5.4. Интерполяционная формула Ньютона №2
- •5.5. Интерполяционный многочлен Лагранжа
- •5.5.1. Вычисление лагранжевых коэффициентов
- •5.5.2. Схема Эйткина
- •5.5.3. Остаточный член формулы Лагранжа
- •5.6. Обратное интерполирование
- •5.6.1 Итерационные методы для обратного интерполирования
- •Глава 6. Аппроксимация функций с помощью сплайнов
- •6.1. Кубические сплайны
- •Глава 7. Методы обработки экспериментальных данных
- •7.1 Построение эмпирической формулы.
- •7.2. Метод выбранных точек (метод натянутой нити)
- •7.3 Метод средних
- •7.4. Метод наименьших квадратов
- •7.5. Метод выравнивания
- •7.6. Метод наименьших квадратов для полиномов
- •Глава 8. Численное интегрирование
- •8.1. Квадратурные формулы Ньютона-Котеса
- •8.1. Формула трапеций и ее остаточный член
- •8.2. Общая формула трапеций и ее остаточный член
- •8 .3 Формула Симпсона и ее остаточный член
- •8.4. Общая формула Симпсона и ее остаточный член
- •8.5. Формулы Ньютона-Котеса высших порядков
- •8.6. Квадратурная формула Чебышева
- •8.7. Квадратурная формула Гаусса
- •Глава 9. Приближенное решение обыкновенных дифференциальных уравнений
- •9.1. Аналитические методы
- •9.1.1. Метод последовательного дифференцирования
- •9.1.2. Метод последовательных приближений.
- •9.1.3 Метод неопределенных коэффициентов.
- •9.2. Численные методы
- •9.2.1. Метод Эйлера
- •9.2.2. Модифицированные методы Эйлера Первый улучшенный метод Эйлера
- •Второй улучшенный метод Эйлера
- •Третий улучшенный метод Эйлера
- •9.2.3. Метод Рунге-Кутта для уравнений первого порядка
- •Список литературы
2.1.2. Графический метод отделения корней
Действительные корни уравнения f(x)=0 можно определить как абсциссы точек пересечения графика функции с осью Ох. Если уравнение не имеет близких между собой корней, то этим способом корни легко определяются. На практике часто удобно тождественно преобразовать уравнение к виду , где и - более простые функции, чем функция . Тогда, построив графики и , искомые корни получаются как абсциссы точек пересечения этих графиков.
Пример 2.2
Отделить графически корни уравнения x·ln(x) -1=0. Преобразуем уравнение к виду 1/x=ln(x) и построим графики.
Рис 2.1. Графический метод отделения корней
Из графика видно, что .
2.2. Уточнение приближенных корней
2.2.1. Метод половинного деления
Пусть уравнение (2.1) имеет на отрезке [a,b] один корень, а функция f(x) на данном отрезке непрерывна и f(a)·f(b)<0. Разделим отрезок [a,b] пополам точкой x1=(а+b)/2. Если f(x1)≠0 , то для продолжения вычислений выберем ту часть промежутка, где знаки функции различны. Концы полученного отрезка обозначим [a1,b1] и снова разделим отрезок [a1,b1] пополам точкой x2=( a1+ b1)/2 и т. д. В результате на каком-то этапе получим или точный корень уравнения(2.1) или бесконечную последовательность вложенных друг в друга отрезков [a1,b1], [a2,b2],… [an,bn],…таких, что
f(an)·f(bn)<0, (n=1,2,…) (2.2)
bn - an=2 -n·(b-a). (2.3)
Так как левые концы a1, a2,… ,an образуют монотонную неубывающую ограниченную последовательность, а правые концы b1 b2,…,bn образуют монотонную невозрастающую ограниченную последовательность, а расстояние между ними в силу (2.2) стремится к нулю, то у последовательностей существует общий предел . Число ξ, которое является общим пределом последовательностей {an} и {bn}, точный корень уравнения (2.1). Оценим погрешность решения на n-м шаге. Считаем до тех пор, пока длина промежутка не станет меньше заданной точности ε.
В качестве ответа возьмем середину отрезка [an,bn].
.
Рис.2.2. К объяснению метода половинного деления
Метод половинного деления практически удобно применять для грубого нахождения корня данного уравнения, так как при увеличении точности значительно возрастает объём вычислений.
Пример 2.3. Найти, используя пакет Matchcad, методом половинного деления корень уравнения x4-x3-2x2+3x-3=0 на промежутке [1,2]
Функция koren(a,b,ε) возвращает длину отрезка, который будет меньше заданной точности ε, и значение корня на этом промежутке, если на концах отрезка [a,b] функция имеет противоположные знаки, или сообщение об отсутствии корня, в противном случае.
Метод легко реализуется на компьютере. Далее приводится листинг программы, написанной на языке, встроенном в систему Mathcad.
Рис. 2.3. Листинг программы в Mathcad, реализующей метод половинного деления для примера 2.3
2.2.2 Метод хорд
В этом методе нелинейная функция f(x) на отделенном промежутке
[a,b] заменяется хордой, проходящей через точки (a,f(a))и (b,f(b))
Рис.2.4. Метод хорд. Неподвижен правый конец промежутка b
Уравнение хорды: . Найдем точку пересечения хорды с горизонтальной осью. Полагая и , получим
.
Точку x1 принимаем за новую границу отрезка, где содержится корень. Через эту точку с координатами (x1,f(x1)) и соответствующую границу предыдущего интервала (b,f(b)) опять проводим хорду, находим и т.д., получая последовательность x1,x2,x3,…xn,…, сходящуюся к корню уравнения.
Вторая производная сохраняет постоянный знак на . Следовательно, возможны два случая. Если f(b)·f "(b)>0, то хорда имеет правый фиксированный конец, причем последовательность x0,x1,…xn приближается к корню слева. За начальное приближение x0, естественно, берут a
; ; ;
.
Рис.2.5. Метод хорд. Неподвижен левый конец промежутка a
Если f(a)·f "(a)>0, то хорда имеет левый фиксированный конец, причем последовательность x0,x1,…xn … приближается к корню справа. За начальное приближение x0, берут b
; ; ;
.
Для оценки точности можно воспользоваться формулой
,
где -точный корень, - приближенный корень, , на промежутке [a,b]. Считаем до тех пор пока, не выполнится условие . Если имеет место неравенство , то счет можно прекратить, когда.
Пример 2.4. Найти методом хорд корень уравнения x4-x-1=0
Решение находим, используя пакет Mathcad.
Функция монотонна на промежутках (-∞, 0.63), (0.63, ∞) и меняет на концах промежутков знак. Уравнение имеет два корня. Сузим промежутки отделения корней методом проб, т.е. подстановкой.
Первый корень принадлежит промежутку (-1,-0.5)
Будем находить корень на промежутке (-1,-0.5)
Вторая производная всюду положительна, функция положительна в точке a = -1, значит, этот конец неподвижен.
-максимальное, a -минимальное значение модуля производной
на промежутке
нужно учитывать при оценке точности решения,
Нашли корень исходного уравнения с точностью .
Рис. 2.6. Вычисления в Mathcad, реализующие метод хорд для примера 2.4