Министерство образования и науки РФ

Государственное образовательное учреждение

высшего профессионального образования

«Омский государственный технический университет»

Численные методы

Учебное пособие

для студентов второго и третьего курсов всех форм обучения

Омск - 2010

Составитель: Котюргина Александра Станиславовна, доцент

Печатается по решению редакционно-издательского отдела

Омского государственного технического университета.

Решение нелинейных уравнений

1. Постановка задачи

Пусть дана некоторая функция и требуется найти все или некоторые значения, для которых.

Значение , при котором, называетсякорнем (или решением) уравнения. Относительно функции часто предполагается, чтодважды непрерывно дифференцируема в окрестности корня.

Корень уравнения называетсяпростым, если первая производная функции в точкене равна нулю, т. е.. Если же, то кореньназываетсякратным корнем.

Геометрически корень уравнения есть точка пересечения графика функции с осью абсцисс. На рис. 1 изображен график функции, имеющей четыре корня: два простыхи два кратных.

Рис. 1.

Большинство методов решения уравнения ориентировано на отыскание простых корней.

2. Основные этапы отыскания решения

В процессе приближенного отыскания корней уравнения обычно выделяют два этапа: локализация (или отделение) корня и уточнение корня.

Локализация корня заключается в определении отрезка , содержащего один и только один корень. Не существует универсального алгоритма локализации корня. Иногда удобно бывает локализовать корень с помощью построения графика или таблицы значений функции. На наличие корня на отрезкеуказывает различие знаков функции на концах отрезка. Основанием для этого служит следующая теорема.

Теорема. Если функция непрерывна на отрезкеи принимает на

его концах значения разных знаков, так, что , то отрезоксодержит, по крайней мере, один корень уравнения.

Однако корень четной кратности таким образом локализовать нельзя, так как в окрестности такого корня функция имеет постоянный знак. На этапе уточнения корня вычисляют приближенное значение корня с заданной точностью. Приближенное значение корня уточняют с помощью различных итерационных методов. Суть этих методов состоит в последовательном вычислении значений, которые являются приближениями к корню.

3. Метод половинного деления

Метод половинного является самым простым и надежным способом решения нелинейного уравнения. Пусть из предварительного анализа известно, что корень уравнения находится на отрезке , т. е., так, что. Пусть функциянепрерывна на отрезкеи принимает на концах отрезка значения разных знаков, т.е..

Разделим отрезок пополам. Получим точку. Вычислим значение функции в этой точке:. Если, то- искомый корень, и задача решена. Если, то- число определённого знака:, либо. Тогда либо на концах отрезка, либо на концах отрезказначения функцииимеют разные знаки. Обозначим такой отрезок. Очевидно, что, и длина отрезкав два раза меньше, чем длина отрезка. Поступим аналогично с отрезком. В результате получим либо корень, либо новый отрезоки т.д. (рис. 2).

Рис. 2.

Середина -го отрезка. Очевидно, что длина отрезкабудет равна, а так как, то

. (1)

Критерий окончания. Из соотношения (1) следует, что при заданной точности приближения вычисления заканчиваются, когда будет выполнено неравенство или неравенство . Таким образом, количество итераций можно определить заранее. За приближенное значение корня берется величина.

Пример. Найдем приближенно с точностью. Эта задача эквивалентна решению уравнения, или нахождению нуля функции. В качестве начального отрезкавозьмем отрезок. На концах этого отрезка функция принимает значения с разными знаками:. Найдем числоделений отрезка, необходимых для достижения требуемой точности. Имеем:

.

Следовательно, не позднее 6-го деления найдем с требуемой точностью,. Результаты вычислений представлены в таблице 1.

Таблица 1

	0	1	2	3	4	5	6
	1,0000	1,0000	1,0000	1,1250	1,1250	1,1406	1,1406
	2,0000	1,5000	1,2500	1,2500	1,1875	1,1875	1,1562
	1,5000	1,2500	1,1250	1,1875	1,1406	1,1562	1,1484
Зн	-	-	-	-	-	-	-
Зн	+	+	+	+	+	+	+
	5,5938	0,7585	-0,2959	0,1812	-0,0691	0,0532	-0,0078
–	1,0000	0,5000	0,2500	0,1250	0,0625	0,0312	0,0156