Koshi masalasining ko’p qadamli ayirmali metodlar


Download 0.5 Mb.
bet6/7
Sana25.04.2023
Hajmi0.5 Mb.
#1397539
1   2   3   4   5   6   7
Bog'liq
Ro\'zmetova Maryamjon 151-matematika.pdf.


sm
s εs

i0
i nmis

jm
nmj j j

Qarayotgan D sohaning aniqlanishiga ko’ra
ls L
tengsizlik bajariladi, bundan


m
tashqari, ma’lumki
Г m  1.

Demak,n
oldidagi koeffisentni pastdan quyidagicha baholash mumkin:
1  hb0ln 1  hL b0 .

Shuning uchun ham h yetarlicha kichik


h 1

bo’lganda


 oldidagi





L b0

n
koeffisentni doim musbat qilish mumkin, bu esa (1.63) tenglikni quyidagicha yozishga imkon beradi:

m1
s
minm,ns
m

  1
Гn
s0

  • hls


ims
Гnmis s


(1.64)

n  


1  hb0ln
n1
minm,ns
m
n1
m

hlss
bs Гnmi Гnmj hrj
 

sm i0 j0
Oxirgi tenglik xn nuqtadagi taqribiy yechimning n xatoligi (1.52) hisoblash

formulasini parametrlari
ε0 , ε1,...,εm1
dastlabki xatoliklar, hisoblashning hamma

pog’onalaridagi formulaning xatoligi, yaxlitlash xatoligi hamda xm , xm1..., xn1
nuqtalardagi taqribiy yechimning xatoliklari oqrali ifodalanadi. Bizni n da n
xatolikning raftori qiziqtiradi. (1.64) formulaga ko’ra n xatoning raftori xususiy

hlda
Г i i  0,1,..., m 1
yoki bunga teng kuchli bo’lgan


n

i i
nnj 0,1,...,k
1;i 1,2,...,q
funksiyalarning raftoriga bog’liqdir. Oldingi

bandda (1.52) formulaning foeffisentlari
m

A1  1  aj  0
j1
Shartni qanoatlantirishini ko’rgan edik, bu esa (1.60) tenglama doimo

 1



yechimga ega ekanligini ko’rsatadi. Shuning uchun ham eng qulay holdan n da

n
Г i i 0,1,..., m 1funksiyalar moduli bo’yicha chegaralangan bo’lishiga umid qilish mumkin.
Agar A 0 xarakteristik tenglamaning barcha , ,2 ,...,q ildizlari moduli
bilan birdan oshmasa va moduli birga teng bo’lganlari karrali bo’lmasa, u holda
ildizlar sharti bajarilgan deymiz.
Ko’rinib turibdiki, nn j j 0,1,...,k 1;i 1,2,...qfunksiyalar
i i
chegralangan bo’lishi uchun ildizlar shartining bajarilishi zarur va yetarlidir.

  1. ta’rif. Agar ildizlar sharti bajarilsa, u holda (1.52) hisoblash metodlari turg’un deyiladi. Biz turg’unlikning ikki xil ta’rifini keltirdik, bu ta’riflarning teng kuchliligini qralayotgan masalalarning xususiy holi bo’lgan aniqmas integrallarni taqribiy hisoblash qoidasi uchun ko’rsatgan edik.

Biz bu yerda 2-ta’rifdan 1-ta’rifning kelib chiqishini, ya’ni ildizlar sharti bajarilganda (1.52) formula yordamida topilgan taqribiy yechim n da (1.1) tenglamaning yechimiga tekis yaqinlashishini ko’rsatamiz.
Aytaylik, ildizlar sharti bajarilsin, u holda


Г
i
n
max
0im1
0nN
i

Г
n
Г h Г
(1.65)

bo’lib, shu bilan birga
h  0 da Г noldan farqli chekli limitga intiladi.

Faraz qilaylik, i i  0,1,...,m 1, rj
bo’lsin:
va  j lar uchun quyidagi baholar o’rinli

n   i  0, m 1; rj


r;  j  ; j m, N.

Bu baholarni va (1.65)ni hisobga olib, (4.64)dan quyidagi bahoga ega bo’lamiz:
n1

bunda
n P  hQ s
sm
T n mhr ,
(1.66)

m m



P
1  hL bi
j1 ,
1  hL b0
ГL bi
Q j1 ,
T Г .

Ko’rinib turibdiki,
n  0 da
P,Q va T miqdorlar chekli limitlarga intiladi.

Endi
n m
X x0 ligini hisobga olib va ushbu
h

P  a,
ha b,T hr X

  • x0 c h

belgilashlarni kiritib, (1.66) tengsizlikni quyidagicha yozib olamiz:
n1

n a b s
sm

  • c,

bunda
m n N. Bu formuladan ketma-ket quyidagilarga ega bo’lamiz:

m a c,

m1
a b m
c a c ba c a c1  b,

m2
a c b m
  a c1  b2 ,


m1
m3
a c1  b3

..............................................................................

n
Barcha n larda n
Demak,
a c1 bnm .
ni baholash uchun yuqoridagi bahoni ko’proq qilib olamiz:

max


mnN
n
a ceQxx0 .

Shunday qilib, taqribiy yechimning xatoligi bo’lamiz:
n uchun quyidagi tekis bahoga ega

P r T X x



eQX x0 .
(1.67)

max n
h  0

mnN
 

Yuqorida ko’rdikki, ildizlar sharti bajarilsa,
h  0 da
P,Q va T chekli



limitga intildai. Shuning uchun ham (1.67)dan ko’ramizki, taqribiy yechim aniq yechimga tekis intilishi uchun

0,
h0
r0,
 0
0
h  0

shartlar bajarilishi kerak.
Yuqoridagi (1.67) baho qo’pol, amaliyotda tatbiqi ham, chunki uning
tarkibidagi P,Q,T miqdorlarni effektiv ravishda baholash qiyin. Ammo bu
bahoning yaxshilik tomoni shundan uboratki, uning yordamida hisoblash jarayonining yaqinlashish shartlarini aniqlash va yaqinlashish tezligini sifat jihatidan baholash mumkin. Xususiy holda bu baho shuni ko’rsatadiki, yechim izlanayotgan

oraliqning uzunligi
X x0
ortishi bilan taqribiy yechimning xatosi tez o’sadi.

Bundan tashqari, (1.67) baho shuni ko’rsatadiki, (1.52) formulaning xatoligi uch qismdan iborat, Xatolikning birinchi qismi dastlabki ma’lumotlarning xatoligi bilan
bog’liq bo’lib, h  0 da  ning raftoriga bog’liq. Odatda, jadvalning bosh qismini
ko’rishda qo’llaniladigan algoritm shunday tanlanadiki, h kichiklashganda 
kichiklashadi. Xatoning ikkinchi qismi dastlabki differensial tenglamani ayirmali
tenglama bilan almashtirishga bog’liq, u h  0 da kamayadi, chunki p -tartibli

approksimasiyaga ega bo’lgan ayirmali sxemalar uchun
r  0hp .

Nihoyat, xatolikning uchinchi qismi (1.52) formula yordamida hisoblash

xatoligiga bog’liq bo’lib,
h  0 da
ning raftoriga bog’liq. Agar hisoblash
h

formulasi tanlangan bo’la, rn
va r larning h ga bog’liqligi qonuni aniq bo’ladi.

Bizning ixtiyorimizda
h, larni tanlash qoladi. Ularni optimal ravishda tanlash

uchun ayrim mulohazalarni aytish mumkin.
Faraz qilaylik, hozircha  belgilangan bo’lsin. Biz h ni kichraytirib,
xatolikning birinchi va ikkinchi qismini kamaytirish hisobidan umumiy xatoni kamaytiramiz. Lekin bu uzoqqa bormaydi, ma’lum paytdan boshlab xato yana o’sib

boradi, chunki h ni kichraytirgan sari
ning miqdori oshib boradi. Shuning uchun
h

ham h ning shunday optimal qiymatini topish mumkinki, bu qiymatda (1.67)ning o’ng tomoni eng kichik bo’ladi. Qo’pol qilib aytganda, bu optimal qiymat shunga olib keladiki, xatolikning uchala qismi bir-biriga teng bo’lishi kerak, Agar biz
(1.67)ning o’ng tomonini yana ham kamaytirmoqchi bo’lsak, u holda hisoblash
aniqligi  ni oshirishimiz kerak   0hm va   0hm bo’lsa, u holda (1.67)

bahoga ko’ra hisoblash jarayoni yaqinlashadi.
hm tartibdagi tezlikda aniq yechimga tekis

Shuni yana bir bor ta’kidlash kerakki, (1.67) baho yaqinlashishini ta’minlash uchun (1.52) ayirmali metod uchun turg’ulik shartlari bajarilishi kerak. Agar bu
shartlar buzilsa, u holda h  0 da (1.67) tengsizlikning o’ng tomoni darajali yoki
ko’rsatkich funksiyadek cheksiz o’sib boradi.
Ko’pincha turg’un hisoblash metodlari orasida qat’iy turg’un metodlari

ajratiladi. Bunday metodlar uchun yana bir qo’shimcha talab qo’yiladi:
  1

aylanada faqat bitta
 1 ildiz yotishi kerak. Tadqiqotlar shuni ko’rsatadiki, qat’iy

turg’un jarayonlarning yaqinlashishi raftori ancha yaxshi bo’ladi.

Yuqorida ko’rilgan Adamsning barcha metodlari qat’iy turg’undir, chunki
A  m1  m  0 xarakteristik tenglama biriga teng bo’lgan bitta tub ildizga va
nolga teng bo’lgan m -karrali ildizga ega.
Ayirmali metodlarni hosil qilish uchun boshqacha yondashish ham mumkin. Misol uchun ushbu

u xn u xn2
xn

Download 0.5 Mb.

Do'stlaringiz bilan baham:
1   2   3   4   5   6   7




Ma'lumotlar bazasi mualliflik huquqi bilan himoyalangan ©fayllar.org 2024
ma'muriyatiga murojaat qiling