Приближённые методы решения алгебраического уравнения

Реферат по курсу численных методов выполнил студент группы РЭ–01-1

Днепропетровский Национальный Университет

Радиофизический факультет

Кафедра физики СВЧ

Днепропетровск 2002

1. Численное решение уравнений с одним неизвестным

В данной работе рассматриваются метода приближённого вычисления действительных корней алгебраического или трансцендентного уравнения

f(x)=0 (1.1)

на заданном отрезке [a, b].

Уравнение называется алгебраическим, если заданная функция есть полином n-ой степени:

f(x) = P(x) = a₀ xn + a₁ xn- 1 + … + a_{n -1} x + a_n = 0, a₀ ¹ 0

Требование a₀ ¹ 0 обязательно, так как при невыполнении этого условия данное уравнение будет на порядок ниже.

Всякое уравнение (1.1) называется трансцендентным, если в нём невозможно явным образом найти неизвестное, а можно лишь приближённо.

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

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

Корнем уравнения (1.1) называется такое число x, где f(x)=0.

При определении приближённых корней уравнения (1.1) необходимо решить две задачи:

отделение корней, т. е. определение достаточно малых промежутков, в каждом из которых заключён один и только один корень уравнения (простой и кратный);

уточнение корней с заданной точностью (верным числом знаков до или после запятой);

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

Вторая задача решается непосредственно в методах рассмотренных ниже.

При графическом отделении корней уравнения (1.1) нужно последнее преобразовать к виду:

j₁ (x)=j₂ (x) (2.1)

и построить графики функций y₁ =j₁ (x), y₂ =j₂ (x).

Действительно, корнями уравнения (1.1)

f(x) = j₁ (x) - j₂ (x) = 0

являются абсциссы точек пересечения этих графиков (и только они).

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

Условия, наложенные на функцию f(x) на отрезке [a, b].

Будем предполагать, что функция f(x) непрерывна на отрезке [a, b] (для метода хорд можно потребовать на интервале) и имеет на этом интервале первую и вторую производные, причём обе они знакопостоянны (в частности отличны от нуля). Будем также предполагать, что функция f(x) принимает на концах отрезка значения разного знака. В силу знакопостоянства первой производной функция f(x) строго монотонна, поэтому при сделанных предположениях уравнение (1.1) имеет в точности один корень на интервале (a, b).

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

Этот метод ещё называется методом вилки.

Нам необходимо найти корень уравнения (1.1) на отрезке [a, b]. Рассмотрим отрезок [x₀ , x₁ ]: [x₀ , x₁ ]Ì[a, b]. Пусть мы нашли такие точки х₀ , х₁ , что f (х₀ ) f(х₁ ) £ 0, т. е. на отрезке [х₀ , х₁ ] лежит не менее одного корня уравнения. Найдём середину отрезка х₂ =(х₀ +х₁ )/2 и вычислим f(х₂ ). Из двух половин отрезка выберем ту, для которой выполняется условие

f (х₂ ) f(х_гран .) £ 0, так как один из корней лежит на этой половине. Затем новый отрезок делим пополам и выберем ту половину, на концах которой функция имеет разные знаки, и т. д. (рис 1.2).

рис. 1.2

Если требуется найти корень с точностью Е, то продолжаем деление пополам до тех пор, пока длина отрезка не станет меньше 2Е. Тогда середина последнего отрезка даст значение корня с требуемой точностью. Дихотомия проста и очень надёжна. К простому корню она сходится для любых непрерывных функций в том числе и не дифференцируемых; при этом она устойчива к ошибкам округления. Скорость сходимости не велика; за одну итерацию точность увеличивается примерно вдвое, т. е. уточнение трёх цифр требует 10 итераций. Зато точность ответа гарантируется.

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

Предположим для определённости, что функция f(x) принимает на левом конце отрезка [a, b] отрицательное значение, а на правом – положительное:

f(a) < 0, f(b) > 0.

Возьмём среднюю точку отрезка [a, b], h=(a+b)/2 и вычислим значение в ней функции f(

Внимание, отключите Adblock

Вы посетили наш сайт со включенным блокировщиком рекламы!
Ссылка для скачивания станет доступной сразу после отключения Adblock!

Скачать