ВУЗы и колледжи

Доброго времени суток! Нужно в Маткаде найти решение системы уравнений методом Ньютона:

x+2(x-y-0.752)^2=1.05
x-y=0.752
ε=0.01
x,y>0
А чё Маткад-то?
Второе в первое подставить, х+2(0,725-0,725)^ =1,05 получить.
К чему ε=0.01, если сразу х=1,05 вытекает.
Во второе ур-е подставить, "у" найти:
1,05-у=0,725, откуда и козе понятно: x,y>0.
...
Невтона ему похотелось, блин-н, с эпсилонами...
Анна Грошева
Анна Грошева
35 229
Лучший ответ
Артём Войс Это господа, претензии к тому человеку, который задал эту задачу решить, предъявляйте.
Размышляю примерно так же...
то, что при этих раскладах все так, как Вы, Николай, написали, мне тоже понятно.
Метод линеаризации Ньютона:
Есть система уравнений.
F(x,y) = 0
G(x,y) = 0
Подставим вместо x и y некоторые приближенные значения a и b:
F(a,b) и G(a,b) уже не равны строго нулю. Предположим, что a и b мало отличаются от точного решения x и y, тогда:
F(x,y) = F(a+[x-a], b+[y-b]) = F(a,b) + Fx(a,b) (x-a) + Fy(a,b) (y-b) = 0
Аналогично: G(x,y) = G(a,b) + Gx(a,b) (x-a) + Fy(a,b) (y-b) = 0
обозначим: x-a = da, y-b = db, получим систему:
Fx(a,b) da + Fy(a,b) db = - F(a,b)
Gx(a,b) da + Gy(a,b) db = - G(a,b)
система линейных уравнений для da, db. решаем, находим da, db. Тогда x = a + da, y = b + db
На основании таких рассуждений строится итерационная процедура:
1) Выбираются x, y
2) Составляется и решается слау относительно dx, dy:
Fx(x,y) dx + Fy(x,y) dy = - F(x,y)
Gx(x,y) dx + Gy(x,y) dy = - G(x,y)
3) Уточняется решение:
x = x + dx, y = y + dy
4) Проверка на достижение требуемой точности:
Если max( |dx|, |dy|) < E, то считается, что требуемая точность E достигнута, завершаем процедуру.
Если max( |dx|, |dy|) > E, то считается, что требуемая точность E не достигнута, возвращаемся к пункту 2.
А насчет реализации именно на маткаде... кроме простых алгебраических операций и циклов вам тут ничего и не понадобится. Как писать циклы, складывать и умножать в маткаде, легко найдете. Удачи =)
Алёна Василенко Попробовал осуществить данный алгоритм. Брал начальные приближения (0,0), (1,1), (1,-1), (500, 100). К третьей итерации уже достигается нужная точность. И кроме последнего начального приближения ответ полностью точный. В последнем до точного сходится еще несколько итераций.

Похожие вопросы