Навигация
Усовершенствованный метод хорд
31486
знаков
0
таблиц
15
изображений
9. Усовершенствованный метод хорд
Если итерационная последовательность, полученная методом хорд, сходится, то скорость сходимости будет такой же, как и у метода итераций, - погрешность значения корня убывает, как геометрическая прогрессия. Существует усовершенствование способа хорд, дающее гораздо более быструю сходимость. В обычном методе хорд мы на каждом шагу используем один из концов отрезка [a, b] последнее получившееся приближение. Вместо этого можно использовать два последних приближения – ведь они ближе к искомому корню, чем концы отрезка [a, b].
рис.1.9 а) б)
Формула, при которой мы используем два последних приближения, имеет вид:
(1.9)
При этом а1 вычисляется по формуле:
а а2 в зависимости от знаков f(a), f(b), f(a1), если f(a)<0, f(b)>0,
, f(a1)<0
, f(a1)>0
Если случайно окажется, что точка а3, вычисленная по формуле (1.9), лежит за пределами отрезка [a, b], то на следующем шаге надо вместо этой точки взять ближайший к ней конец этого отрезка (рис. 1.9, б). Оказывается, что сходимость усовершенствованного метода хорд гораздо быстрее, чем у обычного. Именно, если x - корень уравнения f(x)=0, то:
|an+ 1|<C×|an-x| S, где 
10. Комбинированный метод решения уравнений
При решении уравнений часто комбинируют методы хорд и Ньютона. Если график функции y=f(x) обращён вогнутостью вверх, то находят точки а1 и х1 по формулам:
(1.10)
(2.10)
Если же график функции y=f(x) обращён вогнутостью вниз, то точку а1 находят по формуле (1.10), а точку х1 – по формуле:
(3.10)
Как видно из рис.1.10 а) и б), корень x уравнения f(x)=0 лежит обычно между полученными точками а1 и х1. Применяя снова к этим точкам формулы метода хорд и метода Ньютона, получают новую пару точек а2 и х2 и т. д.
Таким путём получают две последовательности точек а1, а2, а3, …, an, … и x1, x2, x3, … , xn, …, приближаются с разных сторон к искомому корню x. Преимущество описанного метода состоит в том, что при нём получаются приближённые значения как с избытком так и с достатком.
рис.1.10
а) б)
11. Заключительные замечанияСитуация, когда одну и ту же задачу можно решить многими способами, является довольно типичной. В таких случаях естественно возникает необходимость сравнения их между собой.
При оценке эффективности численных методов существенное значение имеют различные свойства:
универсальность;
простота организации вычислительного процесса и контроля над точностью;
скорость сходимости.
Наиболее универсальным является метод деления пополам (дихотомии): он только требует непрерывности функции. Остальные методы накладывают более сильные ограничения. Во многих случаях это преимущество метода вилки может оказаться существенным.
С точки зрения организации вычислительного процесса все виды численного нахождения корней уравнения очень просты. Однако и здесь метод деления пополам обладает некоторым преимуществом. Вычисления можно начинать с любого отрезка [a, b], на концах которого непрерывная функция f(x) принимает значения разных знаков. Процесс будет сходится к корню уравнения f(x)=0, причём на каждом шаге он даёт для корня двустороннюю оценку, по которой легко определить достигнутую точность. Сходимость же метода итераций или касательных зависит от того, насколько удачно выбрано нулевое приближение.
Наибольшей скоростью сходимости обладает метод касательных. В случае, когда подсчёт значений функции f(x) сложен и требует больших затрат машинного времени, это преимущество становится определяющим. На вопрос о том, какой метод – метод итераций или дихотомия даёт большую скорость сходимости, однозначно ответить нельзя. При методе дихотомии знаменатель геометрической прогрессии убывания погрешности равен q=0.5, а при методе хорд он может принимать значения 0<q<1.
Из вышесказанного следует, что ответ на вопрос о наилучшем численном методе решения уравнения не однозначен. Он существенно зависит от того, какую дополнительную информацию о данной функции мы имеем, в соответствие с этим, каким свойствам метода придаём большее значение.
При обосновании метода итераций и метода Ньютона на функции j(x) и f(x), а также на выбор начального приближения х0 накладывались определённые ограничения. Однако при решении конкретных задач проверить их выполнение часто бывает трудно и даже практически не возможно. Функция может не задаваться в виде простой формулы, а находится в результате численного решения некоторой математической задачи, получаться из измерений и проверять «экспериментально»: начинают расчёт и следят за поведением первых членов последовательности {xn}. Если по ним видно, что процесс сходится, то расчёт продолжают, пока не достигнут нужной точности. В противном случае вычисления прекращают и анализируют полученные данные, пытаясь установить причину рассходимости и, в соответствии с ней выбрать другой метод решения задач.
Список литературыА. Н. Тихонов, Д. П. Костомаров «Вводные лекции по прикладной математике»
М. «Наука» 1984
Л. Д. Кудрявцев «Математический анализ т. 2» М. 1984 «Наука»
П. Ф. Фильчаков «Справочник по высшей математике» К. 1973 «Наукова Думка»
Н. Н. Калиткин «Численные методы» М. «Наука» 1978
Н. Я. Виленкин «Итерационные методы» М. «Наука» 1984
Похожие работы
... Мы видим, что одна и та же математическая модель может содержать в себе различные аспекты приложений, что вполне характерно для духа прикладной математики. 3. Методы решения алгебраических уравнений Большинство задач физики, экономики, социологии, биологии и других областей знания приводят к решению алгебраических уравнений или систем уравнений. Несмотря на наличие множества приближённых ...
... решения от численных методов расчёта. Для определения корней уравнения не требуется знания теорий групп Абеля, Галуа, Ли и пр. и применения специальной математической терминологии: колец, полей, идеалов, изоморфизмов и т.д. Для решения алгебраического уравнения n - ой степени нужно только умение решать квадратные уравнения и извлекать корни из комплексного числа. Корни могут быть определены с ...
... «проявляется» лишь в процессе преобразований. Очевидность и «завуалированность» новой переменной мы рассмотрим на конкретных примерах во второй главе данной работы. 2. Возможности применения метода замены неизвестного при решении алгебраических уравнений В этой главе выявим возможности применения метода замены неизвестного при решении алгебраических уравнений в стандартных и нестандартных ...
... D2 ∙ c = p. Отсюда получаем, что: c = D2 ∙ ( p - D1∙ u ) Таким образом, искомые константы найдены. Далее показано как применять этот метод для решения «жестких» краевых задач. Запишем V∙ K(1←0) ∙ ∙ = p. совместно с K(1←0) = K(1←x2) ∙ K(x2←x1) ∙ K(x1←0) и получим: V∙ K(1←x2) ∙ K(x2&# ...
0 комментариев