xi
(0) _
= 0,5; x (0) = -1
x
l
x
x
f2(x1, x2) = x2 +~^ -1 =0;
x 2 +— = 1;
. 2 2
Yechish:
5 f df2
2
a f1( xl(0), x20))
2
= — • x-i
= 1,
0.67,
= xi
dx1 dx2
a f2i x;0), x»0)) dx1
dx-
dx1
dx
'2 3
= 0.5, - f1(x1(0), x (0)) = -0.67, - f2(x|0), x(0)) = 1.87.
2
.(0) V(0K_
ya’ni
J
{
1 • (x1 + 2.36) — 1.53 • (x2 — 2.3) — 3.6;x1 — 1.53x2 — —2.28;
{— 2.36 • (x1 + 2.36) +1 • (x2 — 2.3) = —4.2; — 2.36 x1 + x2 = 3.7;
x((2) =-1.26, x(2) = 0.67; en=2 = max{(-1.26 + 2.36)2,(0.67-2.3)2} = 2.66.
Ushbu misolda berilgan tenglamalar sistemasi ikkita haqiqiy yechimga ega ekanligini quyidagi Maple dasturidan ko‘rish mumkin:
>folve({x-yA2/3 = -0.5, xA2/2+y = 1}, {x,y}); allvalues(%);
{x = 2.895363758, y = -3.191565646} {x = -0.1770315380, y = 0.9843299173}
Takomillashtirilgan Nyuton usuli
Nyuton hisob jarayoni (3.3) ni qurishda har bir qadamda teskari matritsa
W 1 (x(k)) ni hisoblash zarurati noqulaylik tug‘diradi.
Agar W-1(x) matritsa izlanayotgan x* yechimning atrofida uzluksiz va
boshlang‘ich yaqinlashsh x0 izlanayotgan x* yechimga yetarlicha yaqin bo‘lsa, u holda taqriban ushbu
tenglikni o‘rinli deb qabul qilish yoki bu teskari matritsani bir qancha qadamlardan keyin qayta hisoblash mumkin. Bu esa iteratsion jarayonlardagi hisoblashlarni ka- maytirib, quyidagi takomillashtirilgan Nyuton usuli formulasini vujudga keltiradi:
Shuni ta’kidlaymizki, (3.13) va (3.14) jarayonlar uchun dastlabki yaqinlashishlar x(1) va ^(1) o‘zaro mos keladi, ya’ni x(1) = ^(1).
Takomillashtirilgan Nyuton usulining algoritmi (blok-sxemasi 3.4-rasmda tas- virlangan):
x(0) - boshlang‘ich yaqinlashish aniqlanadi.
W-1 (x(0)) matritsani hisoblaymiz.
(3.14) formula yordamida ildizni aniqlashtiramiz.
Agar (3.13) shart bajarilsa, u holda masala yechilgan bo‘ladi va x(k+1) (3.1 ) vektor tenglamaning ildizi deb qabul qilinadi, aks holda 3-qadamga o‘tiladi.
Agar Yakob matritsasidagi hosilalarni hisoblash murakkab yoki uni analitik yo‘l bilan hisoblash mumkin bo‘lmasa, u holda Nyuton usulini qo‘llash murakkablashadi.
117
Dostları ilə paylaş: |