Метод касательных в excel как делать

»
В отличие от метода хорд, в методе
касательных вместо хорды на каждом шаге
проводится касательная к кривой y=F(x)
при x=x_n
и
ищется точка пересечения касательной
с осью абсцисс:

Формула
для (n+1) приближения имеет вид:

Если
F(a)*F»(a)>0,
x₀=a,
в противном случае x₀=b.

Итерационный
процесс продолжается до тех пор, пока
не будет обнаружено, что:

.

Пример:

Пусть
дана задача следующего характера:
Уточнить корни уравнения cos(2x)+x-5=0
методом касательных с точностью до
0,00001.

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

Изначально
необходимо определиться с тем, чему
равно x0: либо a, либо b. Для этого необходимо
выполнить следующие действия:

Найти
производную первого порядка от функции
f(x)=cos(2x)+x-5. Она будет выглядеть следующим
образом: f1(x)=-2sin(2x)+1.

Найти
производную второго порядка от функции
f(x)=cos(2x)+x-5. Она будет выглядеть следующим
образом: f2(x)=-4cos(2x).

Заполнить
ячейки следующим образом (обратить
внимание на названия и номера столбцов
при заполнении — они должны быть такими
же, как на рисунке):

В
итоге получается следующее:

Так
как x0=b, то необходимо выполнить следующие
действия:

Заполнить
ячейки следующим образом (обратить
внимание на названия и номера столбцов
при заполнении — они должны быть такими
же, как на рисунке):

В
ячейку A6 ввести формулу =D5.

Выделить
диапазон ячеек B5:E5 и методом протягивания
заполнить диапазон ячеек B6:E6.

Выделить
диапазон ячеек A6:E5 и методом протягивания
заполнить диапазон нижерасположенных
ячеек до получения в одной из ячеек
столбца E результата (диапазон ячеек
A6:E9).

В
итоге получаем следующее:

Ответ:
Корень уравнения cos(2x)+x-5=0 равен 5,32976.

4. Комбинированный метод хорд и касательных

Для
того чтобы достичь наиболее точной
погрешности, нужно одновременно
использовать методы хорд и касательных.
«По формуле хорд находят x_n+1
, а по формуле касательных —
z_n+1
. Процесс нахождения приближенного
корня прекращается, как только:

В
качестве приближенного корня берут
значение, равное (11):»[2]

Пример:

Пусть
требуется уточнить корни уравнения
cos(2x)+x-5=0 комбинированным методом с
точностью до 0,00001.

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

1. Так
   как в комбинированном методе необходимо
   использовать одну из формул хорд и
   формулу касательных, то для упрощения
   следует ввести следующие обозначения:

• Для
  формул хорд обозначить:

—
xn как mn.

—
Переменная c будет играть роль a или b в
зависимости от ситуации.

—
Остальные обозначения аналогичны
приведенным в формулах хорд, только
учитывая выше введенные переменные.

• Для
  формулы касательных обозначить:

—
xn как nn.

—
Остальные обозначения аналогичны
приведенным в формуле касательных,
только учитывая выше введенные переменные.

1. Найти
   производную первого порядка от функции
   f(x)=cos(2x)+x-5. Она будет выглядеть следующим
   образом: f1(x)=-2sin(2x)+1.

2. Найти
   производную второго порядка от функции
   f(x)=cos(2x)+x-5. Она будет выглядеть следующим
   образом: f2(x)=-4cos(2x).

3. Заполнить
   ячейки следующим образом (обратить
   внимание на названия и номера столбцов
   при заполнении — они должны быть такими
   же, как на рисунке):

1. В
   итоге получается следующее:

1. В
   ячейку G1 ввести e, а в G2 ввести число
   0,00001.

2. В
   ячейку H1 ввести c, а в H2 ввести число 6,
   так как c=b (см. ячейку F2).

3. В
   ячейку I1 ввести f(c), а в I2 ввести формулу
   =COS(2*H2)+H2-5.

4. Заполнить
   ячейки последовательно следующим
   образом (обратить внимание на названия
   и номера столбцов при заполнении — они
   должны быть такими же, как на рисунке):

1. В
   ячейку A6 ввести формулу =E5.

2. В
   ячейку F6 ввести формулу =I5.

3. Выделить
   диапазон ячеек B5:E5 и маркером автозаполнения
   заполнить диапазон ячеек B6:E6.

4. Выделить
   диапазон ячеек G5:K5 и маркером автозаполнения
   заполнить диапазон ячеек G6:K6.

5. Выделить
   диапазон ячеек A6:K6 и методом протягивания
   заполнить все нижестоящие ячейки до
   получения ответа в одной из ячеек
   столбца K (диапазон ячеек A6:K9).

В
итоге получаем следующее:

Ответ:
Корень уравнения cos(2x)+x-5=0 равен 5,32976.

Соседние файлы в папке Численные_методы

• #
• #
• #
• #
• #
• #

Решите нелинейное уравнение методом касательных в экселе

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

Для включения режима циклических вычислений в Excel 2003 в меню Сервис/Параметры/вкладка Вычисления следует поставить флажок Итерации и флажок выбора вида ведения вычислений: автоматически. В MS Excel 2010 следует зайти в меню Файл/Параметры/Формулы и поставить флажок в поле «Включить итеративные вычисления» :

Найдем производную функции f(x)=x-x 3 +1

f’(x)=1-3x 2
В ячейку А3 введем значение а =1, ячейку В3 введем формулу расчета текущего значения х: =ЕСЛИ(B3=0;A3;B3-(B3-СТЕПЕНЬ(B3;3)+1)/(1-3*СТЕПЕНЬ(B3;2)))
В ячейку С3 введем формулу для контроля значения f(x): =B3-СТЕПЕНЬ(B3;3)+1.
Получим корень уравнения в ячейке В3 х=1,325.

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

Получим тот же результат, значит корень на данном промежутке один.

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

При прохождении темы численные методы учащиеся уже умеют работать с электронными таблицами и составлять программы на языке паскаль. Работа комбинированного характера.Расчитана на 40 минут. Цель работы повторить и закрепить навыки паботы с программами EXCEL, ABCPascal. Материал содержит 2 файла. Один содержит теоретический материал, так как он и предлагается ученику . Во 2-м файле пример работы ученика Иванова Ивана.

Скачать:

Вложение	Размер
материал для ученика	57.5 КБ
работа ученика	27 КБ

Предварительный просмотр:

Аналитическое решение некоторых уравнений, содержащих, например тригонометрические функции может быть получено лишь для единичных частных случаев. Так, например, нет способа решить аналитически даже такое простое уравнение, как cos x=x

Численные методы позволяют найти приближенное значение корня с любой заданной точностью.

Приближённое нахождение обычно состоит из двух этапов:

1) отделение корней, т.е. установление возможно точных промежутков [a,b], в которых содержится только один корень уравнения;

2) уточнение приближённых корней, т.е. доведение их до заданной степени точности.

Мы будем рассматривать решения уравнений вида f(x)=0. Функция f(x) определена и непрерывна на отрезке [а.Ь]. Значение х 0 называется корнем уравнения если f(х 0 )=0

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

• Если f(a)* f(b] a, b существует, по крайней мере, один корень
• Если функция y = f(x) непрерывна на отрезке [a, b], и f(a)*f(b) и f ‘(x) на интервале (a, b) сохраняет знак, то внутри отрезка [а, b] существует единственный корень уравнения

Приближённое отделение корней можно провести и графически. Для этого уравнение (1) заменяют равносильным ему уравнением р(х) = ф(х), где функции р(х) и ф(х] более простые, чем функция f(x). Тогда, построив графики функций у = р(х) и у = ф(х), искомые корни получим, как абсциссы точек пересечения этих графиков

Для уточнения корня разделим отрезок [а, b] пополам и вычислим значение функции f(х) в точке x sr =(a+b)/2. Выбираем ту из половин [a, x sr ] или [x sr ,b], на концах которых функция f(x) имеет противоположные знаки.. Продолжаем процесс деления отрезка пополам и проводим то же рассмотрение до тех пор, пока. длина [a,b] станет меньше заданной точности . В последнем случае за приближённое значение корня можно принять любую точку отрезка [a,b] (как правило, берут его середину). Алгоритм высокоэффективен, так как на каждом витке (итерации) интервал поиска сокращается вдвое; следовательно, 10 итераций сократят его в тысячу раз. Сложности могут возникнуть с отделением корня у сложных функций.

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

ПРИМЕР : Определим графически корень уравнения . Пусть f1(х) = х , a и построим графики этих функций. (График). Корень находится на интервале от 1 до 2. Здесь же уточним значение корня с точностью 0,001(на доске шапка таблицы)

Алгоритм для программной реализации

1. а:=левая граница b:= правая граница
2. m:= (a+b)/2 середина
3. определяем f(a) и f(m)
4. если f(a)*f(m)
5. если (a-b)/2>e повторяем , начиная с пункта2

Точки графика функции на концах интервала соединяются хордой. Точка пересечения хорды и оси Ох (х*) и используется в качестве пробной. Далее рассуждаем так же, как и в предыдущем методе: если f(x a ) и f(х*) одного знака на интервале , нижняя граница переносится в точку х*; в противном случае – переносим верхнюю границу. Далее проводим новую хорду и т.д.

Осталось только уточнить, как найти х*. По сути, задача сводится к следующей: через 2 точки с неизвестными координатами (х 1 , у 1 ) и (х 2 , у 2 ) проведена прямая; найти точку пересечения этой прямой и оси Ох.

Запишем уравнение прямой по двум точках:

В точке пересечения этой прямой и оси Ох у=0, а х=х*, то есть

, откуда

процесс вычисления приближённых значений продолжается до тех пор, пока для двух последовательных приближений корня х„ и х п _1 не будет выполняться условие abs(xn-x n-1 ) е — заданная точность

Сходимость метода гораздо выше предыдущего

Алгоритм различается только в пункте вычисления серединной точки- пересечения хорды с осью абсцисс и условия останова (разность между двумя соседними точками пересечения)

Уравнения для самостоятельного решения: (отрезок в excel ищем самостоятельно)

Метод Ньютона в Excel

Как видно, процесс нахождения корней нелинейного уравнения методом Ньютона состоит из следующих этапов:

1. Получения шаблона.
2. Уточнение интервалов в ячейках B2 , B3 .
3. Замена в формуле ЕСЛИ запятую ( , ) на точку с запятой ( ; ).
4. Копирование строки итераций до требуемой точности (столбец E ).

Примечание: столбец A — номер итерации, столбец B — корень уравнения X , столбец C — значение функции F(X) , столбец D — значение первой производной dF(X) , столбец E — точность eps .

Решение уравнений

Аналитическое решение некоторых уравнений, содержащих, например тригонометрические функции может быть получено лишь для единичных частных случаев. Так, например, нет способа решить аналитически даже такое простое уравнение, как cos x=x

Численные методы позволяют найти приближенное значение корня с любой заданной точностью.

Приближённое нахождение обычно состоит из двух этапов:

1) отделение корней, т.е. установление возможно точных промежутков [a,b], в которых содержится только один корень уравнения;

 2) уточнение приближённых корней, т.е. доведение их до заданной степени точности.

Мы будем рассматривать решения уравнений вида f(x)=0. Функция f(x)определена и непрерывна на отрезке [а.Ь]. Значение х0  называется корнем уравнения если  f(х0)=0

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

• Если  f(a)* f(b] < 0 , то внутри отрезка a, b существует, по крайней мере, один корень
• Если функция y = f(x) непрерывна на отрезке [a, b], и  f(a)*f(b)<0 и f‘(x) на интервале (a, b) сохраняет знак, то внутри отрезка [а, b] существует единственный корень уравнения

Приближённое отделение корней можно провести и графически. Для этого уравнение (1) заменяют равносильным ему уравнением  р(х) = ф(х), где функции р(х) и ф(х] более простые, чем функция f(x). Тогда, построив графики функций у = р(х) и у = ф(х), искомые корни получим, как абсциссы точек пересечения этих графиков

Метод дихотомии

Для уточнения корня разделим отрезок [а, b] пополам и вычислим значение функции f(х) в точке xsr=(a+b)/2. Выбираем ту из половин [a, xsr ] или [xsr ,b], на концах которых функция f(x) имеет противоположные знаки.. Продолжаем процесс деления отрезка пополам и проводим то же рассмотрение до тех пор, пока. длина [a,b] станет меньше заданной точности. В последнем случае за приближённое значение корня можно принять любую точку отрезка [a,b] (как правило, берут его середину). Алгоритм высокоэффективен, так как на каждом витке (итерации) интервал поиска сокращается вдвое; следовательно, 10 итераций сократят его в тысячу раз. Сложности могут возникнуть с отделением корня у сложных функций.

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

ПРИМЕР: Определим графически корень уравнения    . Пусть f1(х) = х,  a и построим графики этих функций. (График). Корень находится на интервале от 1 до 2. Здесь же уточним значение корня с точностью 0,001(на доске шапка таблицы)

 Решение в Excel,

Алгоритм для программной реализации

1. а:=левая граница b:= правая граница  
2.  m:= (a+b)/2    середина
3. определяем f(a) и f(m)
4. если f(a)*f(m)<0  то b:=m иначе a:=m
5. если (a-b)/2>e повторяем , начиная с пункта2
   

Метод хорд.

Точки графика функции на концах интервала соединяются хордой. Точка пересечения хорды и оси Ох (х*) и используется в качестве пробной. Далее рассуждаем так же, как и в предыдущем методе: если f(xa) и f(х*) одного знака на интервале , нижняя граница переносится в точку х*; в противном случае – переносим верхнюю границу. Далее проводим новую хорду и т.д.

Осталось только уточнить, как найти х*. По сути, задача сводится к следующей: через 2 точки с неизвестными координатами (х1, у1) и (х2, у2) проведена прямая; найти точку пересечения этой прямой и оси Ох.

Запишем уравнение прямой по двум точках:

        

В точке пересечения этой прямой и оси Ох у=0, а х=х*, то есть

, откуда

процесс вычисления приближённых значений продолжается до тех пор, пока для двух последовательных приближений корня х„ и хп_1 не будет выполняться условие  abs(xn-xn-1)е — заданная точность

Сходимость метода гораздо выше предыдущего

Алгоритм различается только в пункте вычисления серединной точки- пересечения хорды с осью абсцисс и условия останова (разность между двумя соседними точками пересечения)

Решение в Excel

Уравнения для самостоятельного решения: (отрезок в excel ищем самостоятельно)

1.  (х=1,261)
2.  (х=?)

1. sin(x/2)+1=x^2   (х=1,26)

1. y=sin3x*cos5x  (х=?)
2.     (х=0,756)

1. x-cosx=0  (х=0,739)

1. x^2+4sinx=0  (х=-1,933)

1. x=(x+1)3  (х=-2,325)

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

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

Чаще всего речь идет о решении нелинейных уравнений различного типа. Сделать это максимально быстро, особенно с использованием ЭВМ, позволяют численные методы. Они хорошо изучены и давно доказали свою эффективность. К их числу относится и метод касательных Ньютона, которым посвящена эта статья.

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

В данном случае имеется функция g, которая задана на отрезке (a, b) и принимает на нем определенные значения, т. е. каждому x, принадлежащему (a, b) возможно сопоставить конкретное число g(x).

Требуется установить все корни уравнения из промежутка между точками a и b (включая концы), для которых функция обнуляется. Очевидно, что это будут точки пересечения y = g(x) с ОХ.

В некоторых случаях удобнее заменить g(x)=0 на аналогичное, вида g₁(x) = g₂(x). В таком случае в качестве корней выступают абсциссы (значение x) точек пересечения графиков g₁(x) и g₂(x).

Решение нелинейного уравнения важно и для задач оптимизации, для которых условие локального экстремума — обращение в 0 производной функции. Иными словами, такая задача может свестись к поиску корней уравнения p(x) = 0, где p(x) тождественна g'(x).

Методы решения

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

Способы извлечения корней нелинейных уравнений принято делить на аналитические (прямые) и итерационные. В первом случае искомое решение имеет вид формулы, используя которую за некоторое число арифметических операций можно найти значение искомых корней. Подобные методы разработаны для показательных, тригонометрических, логарифмических и простейших алгебраических уравнений. Для остальных же приходится использовать специальные численные методы. Их легко реализовать с помощью ЭВМ, которые позволяют найти корни с требуемой точностью.

К их числу относится и так называемый численный метод касательных. Последний был предложен великим ученым Исааком Ньютоном в конце XVII века. В последующие столетия метод неоднократно совершенствовался.

Локализация

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

Рассмотрим отрезок [a,b]. Если g(x) на нем не имеет разрывов и принимает в концевых точках значения разных знаков, то между a и b или в них самих расположен по крайней мере 1 корень уравнения g(x) = 0. Чтобы он был единственным, требуется, чтобы g(x) на [a,b] была монотонной. Как известно, таким свойством она будет обладать при условии знакопостоянства g’(x).

Говоря иначе, если на [a,b] g(x) не имеет разрывов и монотонно растет или убывает, а ее значения в концевых точках имеют не одинаковые знаки, то на [a, b] существует 1 и только 1 корень g(x).

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

Решение уравнения делением пополам

Прежде чем рассматривать более сложные численные методы (метод касательных и его разновидности) стоит познакомиться с наиболее простым способом выявления корней. Он называется дихотомией и относится к интуитивным методам. Алгоритм нахождения корней основан на теореме о том, что если для g(x), непрерывной на [x₀, x₁] выполняется условие разнознаковости, то на рассматриваемом отрезке есть хотя бы 1 корень g(x) = 0.

Для его обнаружения нужно поделить отрезок [x₀, x₁] пополам и обозначить среднюю точку как x₂. Тогда возможны два варианта: g(x₀) * g(x₂) либо g(x₂) * g(x₁) равны или меньше 0. Выбираем тот, для которого верно одно из этих неравенств. Повторяем процедуру, описанную выше, пока длина [x₀, x₁] не станет меньше некой, заранее выбранной величины, определяющей точность определения корня уравнения на [x₀, x₁].

К достоинствам метода относится его надежность и простота, а недостаток — необходимость изначально выявить точки, в которых g(x) принимает разные знаки, поэтому его нельзя применять для корней, обладающих четной кратностью. Кроме того, он не обобщается на случай системы уравнений или если речь идет о комплексных корнях.

Пример 1

Пусть мы хотим решить уравнение g(x) = 2x⁵ + x — 1 = 0. Чтобы долго не искать подходящий отрезок, строим график, используя, например, известную программу «Эксель». Мы видим, что в качестве отрезка для локализации корня лучше брать значения из промежутка [0,1]. Мы можем быть уверены, что хотя бы один корень искомого уравнения на нем есть.

g'(x) = 10x⁴ + 1, т. е. это монотонно возрастающая функция, поэтому на выбранном отрезке есть только 1 корень.

Подставляем концевые точки в уравнение. Имеем 0 и 1 соответственно. На первом шаге за решение берем точку 0,5. Тогда g(0,5) = -0,4375. Значит ,следующий отрезок для деления пополам будет [0,5, 1]. Его серединная точка — 0,75. В ней значение функции равно 0,226. Берем для рассмотрения отрезок [0,5, 0,75] и его середину, которая находится в точке 0,625. Вычисляем значение g(x) в 0,625. Оно равно -0,11, т. е. отрицательное. Опираясь на этот результат, выбираем отрезок [0,625, 0,75]. Получаем x = 0,6875. Тогда g(x) = -0,00532. Если точность решения 0,01, то можем считать, что искомый результат равен 0,6875.

Теоретическая база

Этот способ нахождения корней методом касательных Ньютона пользуется популярностью из-за его очень быстрой сходимости.

Он основан на том доказанном факте, что если x_n — приближение к корню f(x)=0, таком, что f’ C¹, то следующая апроксимация будет в точке, где обнуляется уравнение касательной к f(x), т. е.

Подставляем x = x_n+1 и обнуляем y.

Тогда алгоритм метода касательных выглядит так:

Пример 2

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

Пусть требуется выявить корни для x³ + 4x — 3 = 0 с некоторой точностью, например 0,001. Как известно, график любой функции в виде многочлена нечетной степени должен хотя бы раз пересекать ось ОХ, т. е. сомневаться в существовании корней не приходится.

Прежде чем решить наш пример методом касательных, строим график f(x) = x³ + 4x — 3 поточечно. Это очень легко сделать, например, используя табличный процессор «Эксель». Из полученного графика будет видно, что на [0,1] происходит его пересечение с осью ОХ и функция y = x³ + 4x — 3 монотонно возрастает. Мы можем быть уверены, что на [0,1] уравнения x³ + 4x — 3 = 0 имеет решение и оно единственное.

Алгоритм

Любое решение уравнений методом касательных начинается с вычисления f ‘(x). Имеем:

Тогда вторая производная будет иметь вид x * 6.

Используя эти выражения, можем записать формулу для выявления корней уравнения по методу касательных в виде:

Далее требуется выбрать начальное приближение, т. е. заняться определением, какую точку считать стартовой (об. x₀) для итерационного процесса_. Рассматриваем концы отрезка [0,1]. Нам подойдет тот, для которого верно условие разнознаковости функции и ее 2-ой производной в x₀. Как видим, при подстановке x₀ = 0 оно нарушено, а вот x₀ = 1 вполне подходит.

Так, как

то если нас интересует решение методом касательных с точностью e, то значение x_n можно считать удовлетворяющим требованиям задачи, при условии выполнения неравенства|f(x_n) / f’(x_n)|< e.

На первом шаге решения задачи методом касательных имеем:

• x₁ = x₀ — (x₀³ + 4x₀ — 3) / (3x₀² + 4) = 1- 0,2857 = 0,71429;
• так как условие не выполняется, идем далее;
• получаем новое значение для x₂, которое равно 0,674;
• замечаем, что отношение значения функции к ее производной в x₂ меньше 0,0063, прекращаем процесс.

Метод касательных в Excel

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

Для этого в «Эксель» нужно создать новую страницу и заполнить ее ячейки следующими формулами:

• в C7 записываем «= СТЕПЕНЬ (B7;3) + 4 * B7 — 3»;
• в D7 вписываем «= 4 + 3 * СТЕПЕНЬ (B7;2)»;
• в E7 записываем «= (СТЕПЕНЬ (B7;3)- 3 + 4 * B7) / (3* СТЕПЕНЬ (B7;2) + 4)»;
• в D7 вписываем выражение «=В7 – Е7»;
• в B8 вписываем формулу-условие «= ЕСЛИ(Е7 < 0,001;»Завершение итераций»; D7)».

Далее требуется «растянуть» формулы в столбцах C, D и E сначала на две строки, а после того как в них появились значения, поступить также со столбцом В.

В конкретной задаче уже в ячейке B10 появится надпись «Завершение итераций», и за решение задачи нужно будет взять число, записанное в ячейке, расположенной на одну строку выше. Для него можно выделить и отдельный «растягиваемый» столбец, введя там формулу-условие, согласно которой там будет записан результат, если содержимое в той или иной ячейке столбца B примет вид «Завершение итераций».

Реализация в Pascal

Попробуем получить решение нелинейного уравнения y = х⁴ – 4 – 2 * х методом касательных в Паскале.

Используем вспомогательную функцию, которая поможет осуществить приближенное вычисление f'(x) = (f(x + delta) — f(x)) / delta. В качестве условия для завершения итерационного процесса выберем выполнение неравенства|x₀-x₁|< некого малого числа. В Паскале его запишем, как abs(x0 — x1)<= epsilon.

Программа примечательна тем, что не требует ручного вычисления производной.

Метод хорд

Рассмотрим еще один способ выявления корней нелинейных уравнений. Процесс итераций заключается в том, что в качестве последовательных приближений к искомому корню для f(x)=0 принимают значения точек пересечения хорды с абсциссами концевых точек a и b с ОХ, обозначаемые, как х₁, …, х_n . Имеем:

Для точки, где хорда пересекается с осью ОХ выражение запишется, как:

Пусть вторая производная положительная при х £ [a,b] (противоположный случай сведется к рассматриваемому, если записать– f(x) = 0). В таком случае график у = f(x) — кривая, выпуклая внизу и расположенная ниже хорды AB. Могут иметь место 2 случая: когда функция имеет положительное значение в точке a или она отрицательное в точке b.

В первом случае в качестве неподвижного выбираем конец a, а за x₀ берем точку b. Тогда последовательные приближения по формуле, представленной выше, образуют последовательность, которая монотонно убывает.

Во втором случае неподвижным является конец b при x₀ = a. Значения х, полученные на каждом шаге итерации, образуют последовательность, которая монотонно возрастает.

Таким образом, можем констатировать, что:

• неподвижным в методе хорд является тот конец отрезка, где не совпадают знаки функции и ее второй производной;
• приближения для корня x — x_m — лежат от него в той стороне, где у f(х) знак, не совпадающий со знаком f» (х).

Итерации можно продолжать, пока не выполнится условия близости корней на этом и предыдущем итерационном шаге по модулю abs(x_m — x_m — ₁)< e.

Модифицированный способ

Комбинированный метод хорд и касательных позволяет устанавливать корни уравнения, приближаясь к ним с разных сторон. Такое значение, при котором график f(x) пересекает OX, позволяет уточнить решение гораздо быстрее, чем по каждому из методов по отдельности.

Предположим, нужно отыскать корни f(x)=0, если они есть на [a, b]. Можно применить любой из описанных выше способов. Однако лучше попробовать их комбинацию, благодаря чему значительно повысится точность корня.

Рассматриваем случай с начальным приближением, соответствующим условию разнознаковости первой и второй производной в конкретной точке х.

В таких условиях решение нелинейных уравнений методом касательных позволяет найти корень с избытком, если x₀=b, а способ с использованием хорд при неподвижном конце b приводит к нахождению приближенного корня с недостатком.

Используются формулы:

Теперь искомый корень х нужно искать в интервале[a₁, b₁]. На следующем шаге нужно применить комбинированный метод уже к этому отрезку. Действуя так далее, получим формулы вида:

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

В качестве условия используется оценочное неравенство| b_n+1 — a_n+1|< e. Иными словами, на практике приходится находить решение при помощи двух методов, но на каждом шаге требуется выяснять, насколько полученные результаты близки друг другу.

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

Комбинированный метод легко реализуется в среде TURBO PASCAL. При большом желании можно попробовать осуществить все вычисления табличным методом в программе «Эксель».

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

При этом каждая строка используется для записи вычислений на конкретном итерационном шаге по двум методам. Затем, в левой части от области решения, на активной рабочей странице выделяется столбец, в котором вписывается результат вычислений модуля разности значений очередного итерационного шага по каждому из методов. Еще один можно использовать для внесения результатов вычислений по формуле расчета логической конструкции «ЕСЛИ», используемой для выяснения, выполняется ли условие или нет.

Теперь вы знаете, как решать сложные уравнения. Метод касательных, как вы уже видели, реализуется достаточно просто, как в Паскале, так и в «Экселе». Поэтому вы всегда сможете установить корни уравнения, которое сложно или невозможно решить посредством формул.