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.
...
Невтона ему похотелось, блин-н, с эпсилонами...
Второе в первое подставить, х+2(0,725-0,725)^ =1,05 получить.
К чему ε=0.01, если сразу х=1,05 вытекает.
Во второе ур-е подставить, "у" найти:
1,05-у=0,725, откуда и козе понятно: x,y>0.
...
Невтона ему похотелось, блин-н, с эпсилонами...
Метод линеаризации Ньютона:
Есть система уравнений.
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.
А насчет реализации именно на маткаде... кроме простых алгебраических операций и циклов вам тут ничего и не понадобится. Как писать циклы, складывать и умножать в маткаде, легко найдете. Удачи =)
Есть система уравнений.
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). К третьей итерации уже достигается нужная точность. И кроме последнего начального приближения ответ полностью точный. В последнем до точного сходится еще несколько итераций.
Похожие вопросы
- Найти решения системы уравнений матричным способом.
- Помогите с решением дифференциальных уравнений методом Бернулли
- решение системы уравнения. x^2-2x+y^2=0 y-lnx=0 Выразим y: y=корень квадратный из (x^2-2x) а как выразить x?
- Доброго времени суток срочно нужна ваша помощь! В университете дали вопросы по ИРК ( история русской культуры) ...
- Доброго времени суток! Дорогие,уважаемые,помогите решить проблему. Не поленитесь заглянуть!Заранее спасибо!
- Доброго времени суток! Очень срочно нужна любая сказка (более-менее известная) , переделанная на бизнес лад.
- Доброго времени суток, категория! Влияет ли портфолио на поступление в ВУЗ? И как его правильно оформить?
- доброго времени суток) уважаемые знатоки административного права,прошу помощи)
- всем доброго времени суток! Нишите в сокращенном варианте про "Принципы трудового права", за раннее спасибо)
- Всем доброго времени суток, у кого хорошо с английским помогите перевести текст, заранее большое спасибо тем, кто поможет!
Размышляю примерно так же...
то, что при этих раскладах все так, как Вы, Николай, написали, мне тоже понятно.