Отделение корней В Excel

Лабораторная работа

Отделение корней нелинейного уравнения

Пусть имеется нелинейное уравнение .

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

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

Этап отделения корня уравнения алгоритмизирован только для некоторых классов уравнений (наиболее известным из которых является класс алгебраических уравнений), поэтому отделение корней нелинейных уравнений, обычно, выполняется «вручную» с использованием всей возможной информации о функции . Часто применяется графический метод отделения действительных корней, обладающий большой наглядностью.

Методы отделения корней

Отделение корней во многих случая можно произвести графически. Учитывая, что действительные корни уравнения F(x)=0 – это есть точки пересечения графика функции y=F(x) с осью абсцисс y=0, нужно построить график функции y=F(x) и на оси OX отметить отрезки, содержащие по одному корню. Но часто для упрощения построения графика функции y=F(x) исходное уравнение заменяют равносильным ему уравнением f₁(x)=f₂(x). Далее строятся графики функций y₁=f₁(x) и y₂=f₂(x), а затем по оси OX отмечаются отрезки, локализующие абсциссы точек пересечения двух графиков.

На практике данный способ реализуется следующим образом: например, требуется отделить корни уравнения cos(2x)+x-5=0 графически на отрезке    [–10;10], используя Excel.

1 способ

Построим график функции f(x)=cos(2x)+x-5 в декартовой системе координат. Для этого нужно:

1.                 Ввести в ячейку A1 текст х.

2.                 Ввести в ячейку B1 текст y=cos(2x)+x-5.

3.                 Ввести в ячейку А2 число -10, а в ячейку А3 число -9.

4.                 Выделить ячейки А2 и А3.

5.                 Навести указатель «мыши» на маркер заполнения в правом нижнем углу рамки, охватывающий выделенный диапазон. Нажать левую кнопку «мыши» и перетащить маркер так, чтобы рамка охватила диапазон ячеек А2:А22.

6.                 Ячейки автоматически заполняются цифрами :

7.                 Ввести в ячейку В2 формулу =COS(2*A2)+A2-5.

8.                 Методом протягивания заполнить диапазон ячеек В3:В22.

9.                 Вызвать "Мастер диаграмм" и выбрать диаграмму график (первый вид), нажать «далее».

10.            Указать диапазон данных, для этого щелкнуть кнопку в поле «Диапазон» и выбрать диапазон данных В2:В22.

11.            Выбрать вкладку ряд, указать имя ряда, щелкнув кнопку в поле «ряд» и выбрав В1.

12.            В поле «подписи по оси Х», щелкнуть кнопку и выбрать диапазон А2:А22, нажать «далее».

13.            Подписать названия осей x и y соответственно, нажать «далее».

14.            Вывести диаграмму на том же листе, что и таблица, нажать кнопку «готово».

В итоге получаем следующее (рисунок 1):

Рисунок 1 – Локализация корня

Анализируя полученное изображение графика, можно сказать, что уравнение cos(2x)+x-5=0 имеет один корень – это видно из пересечения графика функции y=cos(2x)+x-5 с осью OX. Можно выбрать отрезок, содержащий данный корень:[5;6] – отрезок локализации.

2 способ

Для подтверждения полученных данных, можно решить эту же задачу вторым способом. Для этого необходимо уравнение cos(2x)+x-5=0 преобразовать к виду: cos(2x)=5-x. Затем следует каждую часть уравнения рассмотреть как отдельную функцию. Т. е. y₁=cos(2x) и y₂=5-x. Для решения этой задачи в  Excel необходимо выполнить следующие действия:

1.                 Вести в ячейки А1:C1 соответственно текст: «x», «y₁=cos(2_x)», «y₂=5-x».

2.                 A2:A22 заполнить так же как при решении задачи первым способом.

3.                 В В2 ввести формулу =COS(2*A2).

4.                 Методом протягивания заполнить диапазон ячеек В3:В22.

5.                 В С2 ввести =5-A2.

6.                 Методом протягивания заполнить диапазон ячеек С3:С22.

7.                 С помощью Мастера диаграмм выбрать график (первый вид).

8.                 В данном случае диапазон данных следует указывать для построения двух графиков. Для этого нужно нажать кнопку в поле «Диапазон» и выделить ячейки В2:В22, затем нажать Ctrl (на клавиатуре) и выделить следующий диапазон C2:C22.

9.                 Перейти на вкладку ряд, где выбрать именем ряда 1 ячейку В1, а именем ряда 2 ячейку С2.

10.            Подписать ось x , выбрав диапазон А2:А22.

11.            Подписать соответственно оси x и y.

12.            Поместить диаграмму на имеющемся листе.

Результат представлен на рисунке 2: Анализируя полученный результат, можно сказать, что точка пересечения двух графиков попадает на тот же самый отрезок локализации [5;6], что и при решении задачи первым способом.

Рисунок 2 – Локализация корня

Аналитический способ отделения корней

Аналитический способ отделения корней основан на следующей теореме, известной из курса математического анализа.

ТЕОРЕМА: Если непрерывная на  функция , определяющая уравнение , на концах отрезка  принимает значения разных знаков, т.е. , то на этом отрезке содержится, по крайней мере, один корень уравнения. Если же функция  непрерывна и дифференцируема и ее производная сохраняет знак внутри отрезка , то на этом отрезке находится только один корень уравнения.

В случае, когда на концах интервала функция имеет одинаковые знаки, на этом интервале корни либо отсутствуют, либо их четное число.

Для отделения корней аналитическим способом выбирается отрезок , на котором находятся все интересующие вычислителя корни уравнения. Причем на отрезке  функция F(x) определена, непрерывна и  F(a)*F(b)<0. Требуется указать все частичные отрезки , содержащие по одному корню.

Будем вычислять значение функции F(x), начиная с точки x=a, двигаясь вправо с некоторым шагом h. Если F(x)*F(x+h)<0, то на отрезке [x;x+h] существует корень (рисунок 3).

Рисунок 3 – Аналитический способ локализации корней

Если F(x_k)=0, x_k-точный корень. 

Доказательство существования и единственности корня на отрезке.

В качестве примера рассмотрим функцию f(x)=cos(2x)+x-5.

1.      Ввести в ячейки А1, В1 и С1 соответственно «x», «y=cos(2x)+x-5» и «ответ».

2.      В А2 и А3 ввести граничные значения отрезка изоляции.

3.      В В2 ввести формулу =COS(2*A2)+A2-5 и методом протягивания заполнить В3.

4.      В С2 ввести формулу =ЕСЛИ(B2*B3<0;"корень существует";"корень не существует").

Таким образом, на отрезке изоляции корень существует:

Рисунок 4 – Проверка существования корня на отрезке

Для доказательства единственности корня на отрезке изоляции необходимо выполнить следующие действия:

1.                 Продолжить работу в том же документе MS Excel.

2.                 Заполнить D1 и E1 соответственно: «y'=-sin(2x)*2+1» и «ответ» (причем выражение y'=-sin(2x)*2+1 – это производная первого порядка от функции y=cos(2x)+x-5).

3.                 Ввести в D2 формулу =-SIN(2*A2)*2+1 и методом протягивания заполнить D3.

4.                 Ввести в E2 =ЕСЛИ(D2*D3>0;"корень на данном отрезке единственный";"Корень не единственный").

В результате получаем (рисунок 5):

Рисунок 5 – Доказательство единственности корня на отрезке

Таким образом доказано существование и единственность корня на отрезке изоляции.

Рассмотрим решение задачи отделения корней уравнения 
cos(2x)+x-5=0 аналитическим способом с шагом 1 на отрезке [-10;10].

Чтобы отделить корни уравнения аналитическим способом с помощью Excel, необходимо выполнить следующее:

1.                 Заполнить ячейки A1:D1 соответственно: «x», «y=cos(2x)+x-5», «h», «ответ».

2.                 В С2 ввести значение 1.

3.                 Ввести в А2 значение -10.

4.                 Ввести в А3 =A2+$C$2 и методом протягивания заполнить ячейки А4:А22.

5.                 В В2 ввести =COS(2*A2)+A2-5 и методом протягивания заполнить диапазон В3:В22.

6.                

В С3 ввести формулу =ЕСЛИ(B2*B3<0;"Корень на отрезке существует";ЕСЛИ(B3=0;"точный корень";"-")) и методом протягивания заполнить диапазон ячеек С4:С22.

В результате получаем следующее (рисунок 6):

Рисунок 6 – Отделение корня

Следующий пример (рисунок 7) демонстрирует отделение нескольких корней. Пусть исследуется функция cos(x)=0,1x на интервале [–10;10] с шагом 1.

Табулирование функции и построение графика осуществляется как в предыдущих примерах. Видно, что на заданном отрезке имеем 7 корней, находящихся внутри отрезков: [-10;-9]; [-9;-8]; [-5;-4]; [-2;-1]; [1;2]; [5;6]; [7;8].

Рисунок 7 – Отделение корней

Обратим внимание на то, что надежность рассмотренного алгоритма отделения корней уравнения зависит как от характера функции F(x), так и от выбранной величины шага h. Для повышения надежности следует выбирать при отделении корней достаточно малые значения h.

Задание

1. Выполнить отделение корней следующих функций:

№ п/п	Уравнение	A	B
1	tg(x) = 1/x	0	n/2
2	e -x = x	0	1
3	ln(x) = 1/x	1	2
4	2 +ln(x) = 1/x	0	1
5	x - x3 + 1 =0	1	2

6

2x + x5 -1 =0

0

1

2. Выполнить индивидуальные задания

Вариант 1

1	ctg(x) = -x2	1,6	4,5
2	2ln(x)+sin(x) =e x	0	2
3	lg(x) = 2 x-x3	0	10
4	cos(x) = 1/x	0	p/2
5	cos(x) = ln(1+x)	0	p/2

Вариант 2

1	tg(x) = 1/x-x2	1,6	4,5
2	2ln(x) =e x	0	2
3	lg(x) = sin(x)	0	10
4	cos(x)+x2 = 1/x	0	p/2
5	cos2(x) = ln(1+x2)	0	p/2

Вариант 3

1	cos2(x) = x	0	p/2
2	1 - 3 x + x3=0	0	1
3	1 - 3x + x4=0	0	1
4	1 - 3 x + x5=0	0	1
5	tg(x) = 1/x2	0	p/2

Вариант 4

1	ln(x) = sin(x)	1	3
2	e - x = sin(x)	0	p/2
3	e x = 1 /sin(x)	0	p/2
4	e - x = x2	0	1
5	2 + ln(x) = 1/x2	0	1

Вариант 5

1	ln(x) = Sin (x)	0	p/2
2	x - x3 + 2 =0	1	2
3	x +5 = x3	1	2
4	x 2-  0,5 x-2=0	0	0,5
5	ln(x -1)+  x-2=0	1	3

Вариант 6

1	ln(x+3) = Sin (x)	0	p/2
2	2x - x3 + 3 =0	1	2
3	x +8 = x3+x2	1	2
4	x -  0,5 x2+4=0	0	0,5
5	(x -1)2-  x=0	1	3

Вариант 7

1	x +0,5 = e -x	0	1
2	2 - x + x3=0	-2	0
3	sin(x) = 1/x	0	p/2
4	sin(x) = x/2	п/2	п
5	ln(x) = e -x	0	2

Вариант 8

1	lg(x) = e - x	0	1
2	cos(x) = x	0	p/2
3	cos(x) = ln(x)	0	p/2
4	cos(x) = tg(x)	0	p/2
5	cos(x) = x3	0	p/2

Вариант 9

1	1 - 5 x + x3=0	0	1
2	1 - 5 x + x4=0	0	1
3	1 - 3 x + x5=0	1	2
4	4cos(x) = x+x2	0	p/2
5	cos(x) = ln(x)-x2	0	p/2

Вариант 10

1	ln(x) = cos (x)	0	p/2
2	x2 - x3 + 2 =0	1	2
3	cosx +5 = x3	1	2
4	x  0,5-   x=0	0	0,5
5	(x -1)2-   x=5	1	3

Вариант 11

1	ln(x+3) =sin (x)-x	0	p/2
2	2x - x3 + 3 =ln(x)	1	2
3	cos(x) +2 = x3+x2	1	2
4	3x -  0,5 x2+cosx=0	0	0,5
5	(x -1)2-  x+tg(x-1)=0	1	3

Вариант 12

1	x2 +0,5 = e -x	0	1
2	2 – sin(x) + x3=0	-2	0
3	sin(x) = 1/x-x2	0	p/2
4	sin(x) = x/2+cos(x)	п/2	п
5	ln(x)-x = e -x	0	2

Вариант 13

1	lg(x)+2sin2(x) = e - x	0	1
2	cos(x+p/2) = x+x2	0	p/2
3	cos(x) = ln(x)-x2	0	p/2
4	cos(x) = tg(x)+sin(2x-p/2)	0	p/2
5	cos(x) = x3-log(x2)	0	p/2

Вариант 14

1	1 - 5 cos(x) + x3=0	0	1
2	1 - 5 tg(x) + x4=0	0	1
3	1 – 3( x-2)2 + x5=0	1	2
4	4cos(x) = cos(x)+x2	0	p/2
5	cos(x) = ln(x)-x2+x	0	p/2