- •4.8. Методы обращения матриц
- •4.10. Методы решения слау с трехдиагональной матрицей Рассмотрим систему трехточечных уравнений
- •5. Полная и частичная проблемы собственных значений.
- •5.1. Постановка задачи
- •5.2. Метод а. Н. Крылова
- •Решение систем нелинейных уравнений
- •6.1. Метод простой итерации
- •6.2. Метод Ньютона
- •Основная литература
6.2. Метод Ньютона
Метод Ньютона применяется к решению систем уравнений вида
fi(х1..., хр) = 0, I = l , ..., р, (6.10)
или, в векторной форме,
F(x) = 0. (6.11)
Введем матрицу Якоби J(х) для функций fi(х), I = l, ..., р, которые будем предполагать непрерывно дифференцируемыми:
…
J(x) = … . (6.12)
……………………………
…
Пусть задано начальное приближение х(0). Вместо нелинейного уравнения (6.11) решаем линейное уравнение
F(x(0)) – J(х(0)) (х-х(0)) = 0. (6.13)
Если det J (х(0)) ≠ 0, то уравнение (6.13) имеет единственное решение, которое обозначим х(1). Здесь удобно решать уравнение (6.13) относительно х(0) = х - х(0), а затем вычислять x(l) = x(0) + x(0). Если найдено x(т) то х(m + 1) вычисляем по формуле x(m+1) = х(m) + х(m), а поправку х(m) = ( х1(m), ..., хp (т)) находим из системы
F(x(m)) – J(х(m)) (х-х(m)) = 0.
которая в координатной форме имеет вид
f1(х(m)) + х1(m) + … + хp(m) = 0,
(6.14)
fp(х(m)) + х1(m) + … + хp(m) = 0.
Для системы второго порядка
f (x, y) = 0,
g (x, y) = 0
последовательные приближения по методу Ньютона вычисляются по формулам
xn+1 = xn - , yn+1 = yn - ,
г де
f {хп, уп) fy '(хп, уп)
Аn =
g(xn, yл) gy'(xn, yn) ,
fx '(хп, уп) f (хп, уп)
Bn =
gx'(xn, yn) g(xn, yл) ,
fx '(хп, уп) fy '(хп, уп)
Jn = ≠ 0.
gx'(xn, yn) gy'(xn, yn)
Метод Ньютона сходится, если начальное приближение выбрано удачно и матрица J(х*) невырожденна. На практике итерации обычно оканчивают, если ||x(n+1)- х(n)|| ≤ ε. Для выбора начального приближения применяют графический метод, метод проб, метод табулирования и т. д.
Рис. 6.2
Пример. Найти решение системы
f (x, y) = х 3 - у2 – 1 = 0,
g (x, y) = ху3 - у - 4 = 0
методом Ньютона с точностью ε =10-3.
Р ешение. Графически находим начальное приближение хо=1,5, y0 = 1,5 (рис. 6.2).
Матрица Якоби имеет вид
3 х2 - 2у
J(x,y) = .
у3 3x у2 - 1
Дальнейшие вычисления отражены в табл. 6.2.
Таблица 6.2
п |
хп yn |
f (хп, уп) g(xn, yл) |
Jп |
Ап |
Вn |
0 |
1,5 1,5 |
0,12500 -0,43750 |
71,71875 |
-0,171875 |
-3,3750 |
1 |
1,502397 1,547059 |
-0,002170 + 0,015844 |
77,73277 |
0,0277998 |
0,1153255 |
2 |
1,5020396 1,545570 |
0,0000017 0,000019 |
|
|
|