4.1. CHiziqsiz tenglamalar sistemasini yechishda Nyuton usuli
1.CHiziqsiz ikki noma’lumli tenglamalardan tuzilgan
(4.1)
sistema berilgan bo‘lsin.
Bu sistemaning yechimlari yotgan oraliqlarni aniqlashda grafik usuldan foydalanamiz.
F(x,y)=0 va G(x.y)=0 funktsiyalar grafiklari kesishgan nuqtani o‘z ichiga oluvchi kesmani taqriban aniqlaymiz:
D={a x b, c x d}
Bu kesmada yechimga mos keluvchi nuqtaga iloji boricha yaqin bo‘lgan (x0, y0) nuqtani tanlaymiz. Bu x=x0, y=y0 qiymatlardan foydalanib =0.001 aniqlikda hisoblash algoritmini tuzamiz.
n=1,2,3,... lar uchun berilgan sistemadagi funktsiya va ularning xususiy xosilalarini hisoblab sistema yechimini topamiz:
1) F=F(xn-1, yn-1), Fx =Fx(xn-1, yn-1), Fy =Fy(xn-1, yn-1)
G =G(xn-1.yn-1), Gx =Gx(xn-1.yn-1), Gy =Gy(xn-1.yn-1)
2) J= Fx Gy – Gx Fy
1= F Gy – G Fy
2= Fx G – Gx F
3) xn = xn-1 + 1/J , yn = yn-1 + 2/J ,
4) xn – xn=1< , yn – yn-1<,
bo‘lsa, taqribiy yechimni: xn , yn deb olamiz.
4.1-masala. Ushbu
chiziqsiz tenglamalar sistemasining yechimini Nyuton usuli bilan 0.1 aniqlikda toping.
Yechish. Yechim yotgan Kesmani
D={0.4< x <0.5, -0.76< y < -0.73}
deb olsa bo‘ladi (bunga ishonch hosil qilishni o‘quvchining o‘ziga havola qilamiz). U holda, boshlang‘ich yaqinlashishni: x0=0.4, y0=-0.75 deb olsak bo‘ladi.
Fx=2Cos(2x-y)-1.2 , Gx=1.6x
Fy= -Cos(2x-y) , Gy= 3y
F=F(0.4,-0.75)=0.1198 ,
Fx =Fx(0.4, -0.75)=-1.1584,
Fy =Fy(0.4, -0.75) = -0.0208
G=G(0.4,-0.75)=-0.0282,
Gx =Gx(0.4,-0.75) =0.64 ,
Gy =Gy(0.4, -0.75)=-2.25
J==2.6197
1=0.2701 , 2=0.044
x1=x0+1/ =0.5 , y1=y0+2/=-0.733
x1-x0=0.1=0.1, y1-y0=0.02<0.1
Birinchi yaqinlashish: x1=0.5, y1=-0.733 bo‘lganda
F= -0.0131 , Fx=0.8 , Fy=-1.4502
G=0.059, Gx= -2.191, Gy=0.1251
=3.2199
1=-0.0293 , 2=0.0749
x2=x1+1/ =0.491, y2=y1+2/=-0.710
x2-x1=0.009<0.1, y2-y1=0.023<0.1
bo‘lganidan x 0.5, y -0.71
Endi Nyuton usulini n ta noma’lumli n ta chiziqsiz tenglamalar sistemasini yechish uchun qo‘llaymiz.
Buning uchun quyidagi chiziqsiz tenglamalar sistemasini olamiz.
(4.2)
Bu sistemasini yechimini topish uchun ketma-ket yaqinlashish (iteratsiya) usulidan foydalanamiz. Bu ketma-ketlikni yechimga r- yaqinlashishini quyidagicha yozamiz:
x(p+1)= x(p)- W-1(x(p))f(x(p)) (4.3)
bu formulada:
- x(p)=( x1(p), x2(p),…, xn(p))-boshlang‘ich yoki r- yaqinlashishini bildiradi.
-W-1(x(p)) (4.2) sistemaning chap tamonidagi funktsiyalarning har bir argumenti bo‘yicha olingan 1-tartibli xususiy hosilalarning x(p) r-yaqinlashish qiymati bo‘yicha topilgan sonlardan tuzilgan quyidagi Yakobi W
(4.4)
matritsaga teskari matritsa
- f(x(p)) (3.2) sistemaning chap tamonidagi funktsiyalarning x(p) dagi qiymatilaridan tuzilgan matritsa
(4.3) ketma-ketlikni yechimga yaqinlashishining asosiy sharti:
Do'stlaringiz bilan baham: |