6 метод секущих уточнения корней алгебраических и неалгебраических уравнений

Численные методы решения нелинейных уравнений. Метод хорд.

Численные методы решения нелинейных уравнений. Метод хорд.

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

В отличие от метода половинного деления, метод хорд предлагает, что деление рассматриваемого интервала будет выполняться не в его середине, а в точке пересечения хорды с осью абсцисс (ось — Х). Следует отметить, что под хордой понимается отрезок, который проведен через точки рассматриваемой функции по концам рассматриваемого интервала. Рассматриваемый метод обеспечивает более быстрое нахождение корня, чем метод половинного деления, при условии задания одинакового рассматриваемого интервала.

Геометрически метод хорд эквивалентен замене кривой хордой, проходящей через точки и (см. рис.1.).

Рис.1. Построение отрезка (хорды) к функции .

Уравнение прямой (хорды), которая проходит через точки А и В имеет следующий вид:

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

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

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

или .

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

.

Рис.2. Пояснение к определению погрешности расчета.

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

Алгоритм нахождения корня нелинейного уравнения по методу хорд

1. Найти начальный интервал неопределенности одним из методов отделения корней. З адать погрешность расчета (малое положительное число ) и начальный шаг итерации ( ) .

2. Найти точку пересечения хорды с осью абсцисс:

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

— если выполняется условие , то искомый корень находится внутри левого отрезка положить , ;

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

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

4. Проверяем приближенное значение корня уравнения на предмет заданной точности, в случае:

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

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

Пример решения уравнений методом хорд

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

Вариант решения нелинейного уравнения в программном комплексе MathCAD .

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

Рис.1. Результаты расчета по методу хорд

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

Примечание:

Модификацией данного метода является метод ложного положения ( False Position Method ), который отличается от метода секущих только тем, что всякий раз берутся не последние 2 точки, а те точки, которые находятся вокруг корня.

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

Случай №1: 0,

f»(a)>0″ width=»158″ height=»20″ border=»0″ />

Из первого условия получается, что неподвижной стороной отрезка является – сторона a .

Случай №2: 0″ width=»158″ height=»20″ border=»0″ />

Из второго условия получается, что неподвижной стороной отрезка является – сторона b .

В общем виде, для выявления неподвижного конца можно записать следующее условие: 0″ width=»122″ height=»20″ border=»0″ /> , где или .

Рис. 3. Примеры убывающей или возрастающей функции

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

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

, где k =0,1,2,…

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

, где k =0,1,2,…

Случай сводится к рассматриваемому , если уравнение записать в форме: .

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

Метод секущих

Метод секущих — итерационный численный метод приближённого нахождения корня уравнения.

Немного теории о методе секущих под калькулятором.

Метод секущих

Метод секущих

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

Положим, что у нас есть две точки, x0 и x1, в которых значения функции равны соответственно f(x0) и f(x1). Тогда уравнение прямой, проходящей через эти точки, будет

Для точки пересечения с осью абсцисс (у=0) получим уравнение

Это и есть наша итерационная формула. Графическое отображение метода — на рисунке ниже.

Метод работает и в случае, если начальные точки выбраны по одну и ту же сторону от корня (то есть, корня нет на отрезке между начальными приближениями), но при этом возможны случаи, когда метод не сходится.

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

В качестве критерия останова берут один из следующих:

— значение функции на данной итерации стало меньше заданого ε.

— изменение хk в результате итерации стало меньше заданого ε.

Метод секущих

При нахождении нулей функции f, для которой вычисление f'(x) затруднено, часто лучшим выбором, чем метод Ньютона, является метод секущих. В этом алгоритме начинают с двумя исходными числами x1 и х2. На каждом шаге xk+1 получают из xk и xk-1 как единственный нуль линейной функции, принимающей значения f(xk) в xk и f(xk-1) в xk-1. Эта линейная функция представляет секущую к кривой у = f(х), проходящую через ее точки с абсциссами xk и xk-1 — отсюда название метод секущих.

Пусть — абсциссы концов хорды, — уравнение секущей, содержащей хорду. Найдем коэффициенты и из системы уравнений:

Вычтем из первого уравнения второе:

, затем найдем коэффициенты и :

Уравнение принимает вид:

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

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

Повторять операцию следует до тех пор, пока не станет меньше или равно заданному значению погрешности.

Пример решения задачи методом секущих

В Delphi напишем программу для расчета корней уравнений методом секущих:

procedure TForm1.Button1Click(Sender: TObject);


источники:

http://planetcalc.ru/3707/

http://vuzlit.com/980963/metod_sekuschih