5-lemma. Eyler toʻgʻrilangan usulining lokal xatoligi baholashni qanoatlantiradi:
)
1i
quyidagi
(2)
i1
M h3 , i=0,1,…,N-1, (58)
bu yerda M – chekli oʻzgarmas son boʻlib barcha i lar uchun bir xil va un- ing qiymati M1 – M6 larning qiymatlari orqali topiladi.
Isbot. (58) baholashni chiqarishning gʻoyasi Eyler oshkor usulining
(40) baholashinikiga oʻxshash (bu yerda M oʻzgarmasning qiymati har ikkala usul uchun bir biridan farq qiladi). Qaralayotgan usulning lokal xatoligi uchun ushbu
(2)
i1
y(i) (x
) yi1
(59)
i1
ifodasida yi+1 toʻr yechimni unga teng boʻlgan va hisob ifodasining h ning darajalari boʻyicha Teylor qatoriga yoyilmasi hisob formulasidan topilgan miqdori bilan almashtiramiz, bunda ishonch hosil qilamizki, u h boʻyicha talab qilingan kichiklik darajasiga koʻra cheksiz kichik.
Yuqoridagi (59) formulada
1i
lidagi kabi, ushbu
y (i) ( x )
miqdor, Eylerning oshkor usu-
y (i) (x
i 1 )
f (x, y (i) (x)), (60)
differensial tenglama yordamchi yechimining – Koshi masalasi yechimin- ing xi+1 nuqtadagi qiymati.
Bu qiymatni Teylor formulasi boʻyicha ifodalaymiz, bunda yoyilman- ing uchunchi tartibgacha kichiklikdagi hadlarini saqlab qolamiz:
y( i ) (x
) y(i) (x ) y(i) '(x )h 1 y(i) ''(x )h2 1 y(i) '''(x
h) h3 . (62)
i 1 i
i 2 i 6 i i
(60) va (61) tengliklarga koʻra
y(i) ( x ) y , y (i) '( x ) f ( x , y(i) ( x )) f ( x , y ) . (63)
i 1 i i1 i i i i
Bundan tashqari avval (60) tenglikni differensiallab, keyin esa uni qayta qoʻllab, quyidagi munosabatga kelamiz:
x
y
y(i) ''( x) f '( x, y(i) ( x)) f '( x, y(i) ( x)) y(i) '( x)
x
f '(x, y(i) (x))
f '(x, y(i) (x)) f (x, y(i) (x)) , (64)
y
bu yerda x = x i almashtirish olib, (61) ga koʻra quyidagi tenglikka ega boʻlamiz:
y(i) ''(x) f '(x , y ) f '(x , y ) f (x , y ) . (65)
x i i y i i i i
(63) va (65) lardan foydalanib, (62) yoyilmani quyidagicha yozish mum- kin:
y( i ) (x
) y
f (x , y )h 1 f '(x , y ) f
'(x , y ) f (x , y )h2
i 1 i
i i 2 x i i
y i i i i
i
1 y (i) '''( x
6 i
h)h3
(66)
Lokal xatolik uchun (59) ifodadagi y i+1 miqdor (54) formula orqali beriladi. Bu ifodaning oʻng tarafini ikki oʻzgaruvchili funksiya uchun (xi,yi) nuqtada Teylor formulasi boʻyicha yoysak, quyidagini hosil qilamiz:
y y h f ( x , y )
f ( x
h, y hf (x , y )) y
h f (x , y ) f (x , y )
i 1 i 2 i i
i i i i
i 2 i i i i
f '(x , y ) h f
'(x , y ) h f (x , y ) 1 f
''(~x , ~y ) h2
x i i
y i i
i i 2 xx i i
''( ~x , ~y ) h h f ( x , y ) f
''( ~x , ~y ) h2 f 2 ( x , y ) y
f (x , y ) h
xy i i
i i yy i i
i i i i i
1 f
'(x , y ) f
'(x , y ) f (x , y ) h2 1 f
''(~x , ~y ) 2 f
''(~x , ~y ) f (x , y )
2 x i i
y i i i i
2 xx i i
xy i i i i
y f x y
'h'(~ , ~ ) 2 ( ,
)) 3 , (67)
yy i i i i
bu yerda
~x , ~y
- ikkita (xi,yi) va x, y
(x , y
hf (x , y ))
nuqtalarni
i i i
i 1
i 1 i i i
tutashtiruvchi kesma boʻylab yotuvchi kenglik nuqtalarining koordinata- lari.
(66) va (67) tengliklarni ixchamroq holda yozib, ikki oʻzgaruvchili funksiyalar argumentlarini tashlab yuborib, agar ular x i, y i larda teng boʻlsa, quyidagilarga kelamiz:
y( i ) (x
) y
f h 1 f ' f
' f h2
1 y(i) '''(x h) h3 , (68)
i 1 i
2 x y
6 i i
yi 1
yi
f h 1 f
2 x
' f y
' f h2
h3 , (69)
i
bu yerda quyidagi belgilash qabul qilingan:
1 f ''(~x , ~y ) 2 f
''(~x , ~y ) f f
''(~x , ~y ) f 2 . (70)
i 2 xx i i
xy i i
yy i i
i
(68) va (69) ifodalarni oʻzaro taqqoslab, ularda h boʻyicha nolinchi, birinchi va ikkinchi tartibgacha kichiklikdagi cheksiz kichik miqdorlar bir xil va shuning uchun ularni (59) ifodaga qoʻyganimizda ular oʻzaro qis- qarib ketadi. Demak, Eyler toʻgʻrilangan usulining lokal xatoligi quyidagi- ga teng:
(2)
i 1
1 y(i) '''(x
6 i
i
h) h3
h3 , (71)
shunga koʻra bu h ga nisbatan uchinchi tartibgacha kichiklikdagi cheksiz kichik miqdor. Bunda (71) munosabatning oʻng tarafidagi ikkinchi had (-ih3) ning moduli yuqoridan M 7ˑh3 miqdor bilan baholanadi, bu yerda
(70) ga koʻra M7 ning qiymati quyidagicha:
M7 = (M4+2 M5 M1+ M6(M1)2)/2
Birinchi handing moduli esa xuddi hu tartibli M6h3 cheksiz kichik mi- qdor bilan yuqoridan baholash imkonini beradi, ammo bunda oʻzgarmas MS. Bu oʻzgarmasni topish uchun (64) tenglikning oʻng tarafini differensi- allash lozim va (60) munosabatdan foydalanib, y (i) yechimning uchinchi hosilasini hamda uning ikkinchi tartibli xususiy hosilalarini differensial
tenglamaning oʻng tarafi orqali quyidagicha ifodalash zarur:
yi '''(x)
f xx
''(x, y(i) (x)) f
''(x, y(i) (x)) f (x, y(i) (x))
xy
xy
yy
f ''(x, y(i) (x)) f
''(x, y(i) (x)) f (x, y(i) (x)) f (x, y(i) (x))
x
y
f '(x, y(i) (x))
f '(x, y(i) (x)) f (x, y(i) (x)) f '(x, y(i) (x)), (72)
y
bu yerdan MS oʻzgarmasning qolgan M 1, M2, M3, M4, M5, M6 oʻzgarmaslar orqali ifodasi kelib chiqadi.
Shunday qilib, M = M8 + M7 konstantali (58) baholash oʻrnatildi.
Do'stlaringiz bilan baham: |