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

↑

▷ Содержание

⇐ ПредыдущаяСтр 3 из 7Следующая ⇒

▷ Похожие статьи

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

Два этапа решения:

1. На первом этапе локализуют корни (т.е. прикидывают, в какие интервалы облшасти определения попадают корни, чтобы далее сужать поиск, производя уточнение не на всей области определения, а только в пределах нескольких узких отрезков – это ускоряет процесс поиска решения). Области локализации (отрезки) стараются выбрать минимальными, т.е. так, чтобы в в каждый отрезок локализации попал ТОЛЬКО один корень уравнения f (x) = 0 (иначе есть риск найти не все корни уравнения). Локализация корней может производиться (1) графически или (2) аналитически.

(1) В первом случае строят график функции f (x) и выделяют те промежутки 0Х, где график пересекает ось Ox. Если построение графика y = f (x) затруднительно, но f (x) = 0 представляется в виде f1(x) = f2 (x), то строят графики y = f1(x), y = f2 (x) и за начальное Х берут промежутки/точки пересечения этих графиков.

(2) Для аналитического отделения корней находят все критические точки функции f (x), т.е. точки, в которых производные равны нулю или не существуют. Для этого f (x) дифференцируют, приравнивают производную к нулю и решают полученное уравнение относительно x. Кроме того, определяют все точки, где производная может не существовать. В этих (критических) точках или в непосредственной близости от них определяют знак функции f (xi). Затем строят ряд знаков функции в критических точках, включая в рассмотрение и . Анализируя этот ряд, и по числу смен знаков определяют количество корней и интервалы локализации корней: на левой и правой границах такого интервала функция f (x) должна иметь разные знаки. Одна из проблем, в которую упирается решение задачи локализации корней, – это практическая невозможность точного вычисления значения функций.

1. Графический способ: Построить график f(x) и прикинуть отрезки, где график пересекает ось ОХ (см. рисунок справа – 3 пересечения графика f(x) с осью ОХ говорят о наличии 3 действительных корней у уравнения f(x)=0. Правда, мнимые корни таким графиком выявить невозможно).

2. Если график построить затруднительно, но можно представить f(x) в виде суммы 2 более простых функционалов f(x)=g(x)+h(x), то можно построить графики g(x) и h(x) и прикинуть точки (или отрезки) оси х, где графики пересекаются – там есть корни (Рис. Второй справа ниже, 2 линии – синяя и красная – это 2 функции, в сумму которых оказалось с возможным разложить f(x) (если такая возможность вообще есть) – здесь корнями являются 3 точки пересечения этих графиков друг с другом. Их и надо уточнить до заданного уровня погрешности).

3. Аналитический способ: подставить несколько значений Х в f(x) (включая оба края области определения), прикидывая, какой знак в этих точках имеет (т.е. f(x)>0 или f(x)<0). Если на соседних таких точках знаки различны – значит где-то между ними есть корень (или корни) f(x)=0.

4. Найти стационарныеточки (т.е. те х, где производная f ^/ (x)=0) - они тоже могут быть краями интервалов, где есть корни.

2. На втором этапе численного решения нелинейного уравнения итерационно уточняют корни с заданной точностью. Способы уточнения различны:

1) путем перебора всех возможных значений аргумента (с малым шагом) из отрезка локализации с проверкой наличия решения;

2) либо замена нелинейной функции f более простой функцией (линейной, параболической), близкой к исходной, и поиск ее корня итерационными процедурами;

3) либо сведение нелинейного уравнения f (x) = 0 к виду f1(x) = f2 (x) и обеспечение в нем равенства итерационными процедурами.

Условием окончания процесса решения уравнения (т.е. получения корня x* с заданной погрешностью) может быть одно из двух возможных условий:

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

Численные методы уточнения корней (2 этапа):

1. Метод половинного деления – с простым способом выбора пробной точки путем деления промежутка существования корня пополам. Алгоритм такой: На каждом шаге итерации делят отрезок локализации пополам, за оценку корня берут середину отрезка: . Если значение f(xn) в точности =0, то корень считается найденным, итерации прерывают. Если же значение f(xn) не 0, то смотрят знак. (т.к. f(xn)<0 или f(xn)>0), и берут за новые, более узкие границы интервала локализации те, где f(x) принимает разные знаки. Затем весь блок повторяется. Критерий окончания итераций: , где - заданная погрешность.

2. Метод касательных (Ньютона). . (Начальное приближение решения x₀ берут таким, чтобы , иначе метод может не сойтись к решению (сходимость зависит от близости начального приближения x₀ к корню).) Критерий окончания итерационного процесса уточнения корня: . (Точность задается изначально).

3. Метод простых итераций. f(x) разлагают в форму: x=g(x). - В таком случан решение ищут итерационно: x_n+1=g(x_n). За начальное х₀ берут какое-нибудь значение из интервала локализации корня. Известно, что метод сходится со скоростью геометрической последовательности и справедлива оценка погрешности: . Если можно подобрать такую константу q, чтобы |g^/(x)|<=q<1, то метод гарантированно сходится (иногда может сойтись, даже если условие |g^/(x)|<1 не выполняется).

Можно взять за критерий сходимости итерационного процесса к решению : на отрезке локализации. Критерием окончания итераций считают выполнение условия: где - заданная погрешность.

4. Метод хорд. расчетная формула для уточняемого корня на очередной итерации выглядит так: Метод хорд сходится медленнее, может вовсе не сойтись, во избежание чего начальное значение x₀ (и x₁) берут таким, чтобы . Сходимость зависит от близости начального приближения x0 к корню. Критерий окончания итераций тот же .

5. Метод ложного положения. . Критерий окончания

Метод дихотомии (бисекций).

x1

f(x)

x

f1

x0

f0

x2

f2

Иначе называется методом половинного деления. Пусть задан начальный интервал [ x ₀, x ₁], на котором f (x ₀) f (x ₁) ≤ 0 (т.е. внутри имеется не менее чем один корень). Найдем x ₂ = ½ (x ₀ + x ₁) и вычислим f (x ₂). Если f (x ₀) f (x ₂) ≤ 0, используем для дальнейшего деления отрезок [ x ₀, x ₂], если > 0 – используем для дальнейшего деления отрезок [ x ₁, x ₂], и продолжаем деление пополам.

Итерации продолжаются, пока длина отрезка не станет меньше 2ξ ­– заданной точности. Тогда середина последнего отрезка даст значение корня с требуемой точностью. В качестве иного критерия можно взять | f (x)| ≤ ξ_y.

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

Если на заданном интервале предполагается несколько корней, то существует возможность последовательно исключать найденные корни из рассмотрения. Для этого воспользуемся вспомогательной функцией , где – только что найденный корень. Для функций f (x) и g (x) совпадают все корни, за исключением (в этой точке полюс функции g (x)). Для достижения требуемой точности рекомендуется грубо приблизиться к корню по функции g (x), а затем уточнить его, используя f (x).

Метод хорд.

Идея метода проиллюстрирована рисунком. Задается интервал [ x ₀, x ₁], на котором f (x ₀) f (x ₁) ≤ 0, между точками x ₀ и x ₁ строится хорда, стягивающая f (x). Очередное приближение берется в точке x ₂, где хорда пересекает ось абсцисс. В качестве нового интервала для продолжения итерационного процесса выбирается тот, на концах которого функция имеет разные знаки. Условия выхода из итерационного цикла: или

x1

f(x)

x

f1

x0

f0

x2

f2

| f (x)| ≤ ξ _y.

Для вывода итерационной формулы процесса найдем точку пересечения хорды (описываемой уравнением прямой) с осью абсцисс: ax ₂ + b = 0, где ; b = f (x ₀) ­­– ax ₀.

Отсюда легко выразить .

Метод хорд в большинстве случаев работает быстрее, чем метод дихотомии. Недостатки метода те же, что и в предыдущем случае.

Познавательные статьи:



Техника прыжка в длину с разбега

Организация работы процедурного кабинета

Области применения синхронных машин

Оптимизация по Винеру и Калману