Участник:SKirill/Метод Ньютона решения систем нелинейных уравнений: различия между версиями
Строка 32: | Строка 32: | ||
Основная вычислительная нагрузка алгоритма заключается в решении СЛАУ: | Основная вычислительная нагрузка алгоритма заключается в решении СЛАУ: | ||
− | <math>\frac{\partial F(x^{(k)})}{\partial x}\Delta x^{(k)} = -F(x^{(k)})</math> | + | <math>\frac{\partial F(x^{(k)})}{\partial x}\Delta x^{(k)} = -F(x^{(k)}) (*)</math> |
Для нахождения значения <math>\Delta x^{(k)}</math>, по которому вычисляется значение вектора <math>\overline{x}</math>, на очередной итерации вычисляется значение: <math>x^{(k+1)} = x^{(k)} + \Delta x^{(k)}</math> | Для нахождения значения <math>\Delta x^{(k)}</math>, по которому вычисляется значение вектора <math>\overline{x}</math>, на очередной итерации вычисляется значение: <math>x^{(k+1)} = x^{(k)} + \Delta x^{(k)}</math> |
Версия 12:24, 14 октября 2016
Метод Ньютона решения систем нелинейных уравнений | |
Последовательный алгоритм | |
Объём выходных данных | n-мерный вектор |
Авторы: Шохин К.О., Лебедев А.А.
Содержание
- 1 Свойства и структура алгоритма
- 1.1 Общее описание алгоритма
- 1.2 Математическое описание алгоритма
- 1.3 Вычислительное ядро алгоритма
- 1.4 Макроструктура алгоритма
- 1.5 Схема реализации последовательного алгоритма
- 1.6 Последовательная сложность алгоритма
- 1.7 Информационный граф
- 1.8 Ресурс параллелизма алгоритма
- 1.9 Входные и выходные данные алгоритма
- 1.10 Свойства алгоритма
- 2 Программная реализация алгоритма
- 3 Литература
1 Свойства и структура алгоритма
1.1 Общее описание алгоритма
Метод Ньютона решения систем нелинейных уравнений является обобщением метода Ньютона решения нелинейных уравнений, который основан на идеи линеаризации. Пусть F(x) : R^1 \to R^1 - дифференцируемая функция и необходимо решить уравнение F(x) = 0.
Взяв некоторое x_0 в качестве начального приближения решения, мы можем построить линейную аппроксимацию
F(x) в окрестности x_0 : F(x_0+h) \approx F(x_0)+F^'(x_0)h и решить получающееся линейное уравнение F(x_0 )+F^' (x_0 )h =0.
Таким образом получаем итеративный метод :
x_{k+1} = x_k - {F^'(x_k)}^{-1}F(x_k) , k = 0,1,\ldots
Данный метод был предложен Ньютоном в 1669 году. Более точно, Ньютон оперировал только с полиномами; в выражении для F(x+h) он отбрасывал члены более высокого порядка по h , чем линейные. Ученик Ньютона Рафсон в 1690 г. предложил общую форму метода (т. е. не предполагалось что F(x) обязательно полином и использовалось понятие производной), поэтому часто говорят о методе Ньютона—Рафсона.
Дальнейшее развитие исследований связано с именами таких известных математиков, как Фурье, Коши и другие. Например, Фурье доказал в 1818 г., что метод сходится квадратично в окрестности корня, а Коши (1829, 1847) предложил многомерное обобщение метода и использовал метод для доказательства существования решения уравнения.
1.2 Математическое описание алгоритма
1.3 Вычислительное ядро алгоритма
Основная вычислительная нагрузка алгоритма заключается в решении СЛАУ:
\frac{\partial F(x^{(k)})}{\partial x}\Delta x^{(k)} = -F(x^{(k)}) (*)
Для нахождения значения \Delta x^{(k)}, по которому вычисляется значение вектора \overline{x}, на очередной итерации вычисляется значение: x^{(k+1)} = x^{(k)} + \Delta x^{(k)}
1.4 Макроструктура алгоритма
1.5 Схема реализации последовательного алгоритма
1.6 Последовательная сложность алгоритма
1.7 Информационный граф
1.8 Ресурс параллелизма алгоритма
1.9 Входные и выходные данные алгоритма
1.10 Свойства алгоритма
2 Программная реализация алгоритма
2.1 Масштабируемость алгоритма и его реализации
2.2 Существующие реализации алгоритма
3 Литература
- Тыртышников Е. Е. Методы численного анализа — М., Академия, 2007. - 320 c.
- Бахвалов Н. С., Жидков Н. П., Кобельков. Г. М. — 6-е изд. — М. : БИНОМ. Лаборатория знаний, 2008. — 636 с.
<references \>