Участник:Александр Двойнев/песочница

Материал из MachineLearning.

< Участник:Александр Двойнев(Различия между версиями)
Перейти к: навигация, поиск
(Полностью удалено содержимое страницы)
 
(1 промежуточная версия не показана)
Строка 1: Строка 1:
-
== Введение ==
 
-
Пусть задана функция <tex>f(x)</tex> действительного переменного. Требуется найти корни уравнения
 
-
{{eqno|1}}
 
-
:: <tex>f(x)=0.</tex>
 
-
 
-
Задача нахождения корней уравнения {{eqref|1}} обычно решается в 2 этапа. На первом этапе проводится отделение корней, то есть выделение отрезков, содержащих только один корень. На втором этапе, используя начальное приближение, строится итерационный процесс, позволяющий уточнить значение отыскиваемого корня.
 
-
 
-
== Изложение метода ==
 
-
 
-
Метод секущих получается из [[Метод касательных (Ньютона-Рафсона)|метода касательных]] заменой <tex>f'(x^k)</tex> разностным приближением:
 
-
:: <tex>f'(x^k) \approx \frac{f(x^k)-f(x^{k-1})}{x^k-x^{k-1}}.</tex>
 
-
В результате получим формулу итерационного процесса:
 
-
{{eqno|1}}
 
-
:: <tex>x^{k+1}=x^k-\frac{x^k-x^{k-1}}{f(x^k)-f(x^{k-1})}f(x^k), \;k=1,2,\ldots</tex>
 
-
 
-
Метод секущих является ''двухшаговым'', то есть новое приближение <tex>x^{k+1}</tex> определяется двумя предыдущими итерациями <tex>x^k</tex> и <tex>x^{k-1}.</tex> В методе {{eqref|1}} необходимо задавать ''два'' начальных приближения <tex>x^0</tex> и <tex>x^1.</tex>
 
-
 
-
Скорость сходимости метода будет ''линейной'': <tex>|x^{k+1}-x*|=O(k^k-x*).</tex>
 
-
 
-
=== Геометрическая интерпретация ===
 
-
[[Изображение:sekush.JPG|thumb|200px|Метод секущих]]
 
-
Заметим, что уравнение для секущей, проходящей через точки <tex>M'(x^{k-1},f(x^{k-1}))</tex> и <tex>M''(x^k, f(x^k))</tex>, будет выглядть так:
 
-
:: <tex>\frac{y-f(x^k)}{x-x^k}=\frac{f(x^k)-f(x^{k-1})}{x^k-x^{k-1}}.</tex>
 
-
Положив <tex>y=0</tex> и <tex>x=x^{k+1},</tex> можно получить формулу {{eqref|1}}. Это означает, что <tex>x^{k+1}</tex> — это абсцисса точки пересечения нашей секущей с осью ОХ. Иначе говоря, на отрезке <tex>[x^{k-1},x^k]</tex> функция <tex>f(x)</tex> интерполируется многочленом первой степени и за очередное приближение <tex>x^{k+1}</tex> принимается корень этого многочлена.
 
-
 
-
== Числовой пример ==
 
-
 
-
[[Изображение:kasatex.jpg|thumb|200px|График функции]]
 
-
Рассмотрим функцию <tex>f(x)= cos(x)-x+1.</tex> С помощью метода {{eqref|1}} найдем корень уравнения <tex>f(x)=0.</tex> Исходный код программы, ищущей корень уравнения методом секущих, выложен в разделе «'''Файлы'''».
 
-
 
-
Возьмём в качестве начальных приближений <tex>x^0=4, x^1=3</tex> и точность <tex>\eps=10^{-6}.</tex> В итоге за 8 итераций получим корень <tex>x* \approx 1.283429.</tex>
 
-
 
-
== Рекомендации программисту ==
 
-
 
-
=== Критерий останова ===
 
-
 
-
Как правило, берут один из следующих критериев останова:
 
-
# <tex>f(x^k)< \eps</tex> — значение функции на данной итерации стало меньше заданого ε.
 
-
# <tex>\left|x^k-x^{k-1}\right| < \eps</tex> — изменение х<sup>k</sup> в результате итерации стало меньше заданого ε.
 
-
 
-
=== Ошибки округления ===
 
-
 
-
В методе секущих, как и в других итерационных методах решения уравнений, ошибка округления ''не накапливается''. Общая ошибка округления равна ошибке, возникшей в последней итерации, и не зависит от арифметических операций, выполнявшихся в предыдущих итерациях.
 
-
 
-
== Заключение ==
 
-
 
-
== Файлы ==
 
-
 
-
[[Медиа:Sekush.zip | sekush.zip]]
 
-
 
-
== Ссылки ==
 
-
* [[Практикум ММП ВМК, 4й курс, осень 2008|Практикум ММП ВМК, 4й курс, осень 2008]]
 
-
* [[Метод касательных (Ньютона-Рафсона)]]
 
-
 
-
== Список литературы ==
 
-
* [http://mmphome.1gb.ru/index.php?pid=show&id=79 Н. В. Соснин. Численные методы. Конспект лекций (сост. Д. В. Ховратович, Е. А. Попов)]
 
-
* Самаский А. А., Гулин А. В. Численные Методы. Учеб. пособие для вузов. — М.:Наука, 1989.
 
-
 
-
<!-- {{Stub|}} -->
 

Текущая версия

Личные инструменты