II.BOB. INTERPOLYATSION KUBATUR FORMULALAR.
2.1. Interpolyatsion kvadratur formulalar uchun algoritmva dasturlar.
Buyuk matematik Gauss kvadratura nazariyasiga butunlay yangi va juda muhim g’oya’ni kiritdiki, u amaliy analizning tub sohalari rivojlanishi uchun asos bo’lib qoldi. Faraz qilaylik, ba’zi bir y = f ( x) integrallanuvchi funksiya x o’zgaruvchining uzliksiz oraliqni xar bir nuqtasida emas balkim, shu oraliqda yotuvchi maxsus tanlangan x1 , x 2 , x 3 ,..., xn nuqtalarda berilgan bo’lsin. Biz bu yerda faqat chekli oraliqni qaraymiz. Shuning uchun uni darhol normalab
qo’yamiz.
Oraliqni
1 < x < 1 (2.1)
ga keltiramiz va x1 , x 2 , x 3 ,..., xn nuqtalar ham qaysikim, y = f ( x) funksiya berilgan oraliqdategishli bo’lsin. Umuman olganda n ning kattabo’lishidan qat’iy
nazar,
y 1 = f ( x 1 ) , y 2 = f ( x 2 ) , …, yn = f ( xn ) (2.2)
ordinatalar f ( x) funksiyani aniqlash uchun yetarli emas. Lekin biz f ( x) funksiyani oraliq nuqtalari uchun integrallashga harakat qilamiz. Shu maqsadda x ning darajali funksiyalaridan foydalanamiz. Biz shunday n 1 darajali Pn 1 (x) ko’phad topishimiz mumkinki , u ham xn nuqtalarda yn qiymatga ega bo’ladi. Odatdachekliayirmalarni hisoblashda berilgan x = xn nuqtalartengtaqsimlangan
qilib taqsimlanadi.
Gaussning g’oyasi shundan iboratki nuqtalarning holatini oldindan belgilamasdan o’shanday sondagi ordinatalar bilan yuqori aniqlikka erishish mumkinligi kabi, bu yerda nuqtalar shunday joylashtiriladiki, natijada eng yaxshi natijalar olinadi. Bu yo’lda Gauss kvadratur formulalarning nafaqat eng yuqori aniqlikka erishdi, balkim bu jarayon ko’phadlar bilan teng taqsimli
interpolyatsiyalashda xavfdan ham xolidir. Qaysikim bu xavf u davrda ham ma’lumemasdi. Faraz qilaylik x = xk interpolyatsiyalashnuqtalaritamomanerkin bo’lsin va biz bu nuqtalarda y 1 , y 2 , y 3 ,..., yn qiymatlarni qabul qiladigan
U = Pn- 1 ( x) ko’phadni topamiz. Bu masalani hal qiladigan formula Lagranjning
interpolyatsion formulasi sifatidama’lum. U
Fn ( x) = ( x - x 1 )( x - x 2 )...( x - xn ) .
fundamental ko’phadni qurishga va uni ketma-ket xar bir
bo’lishga asoslangandir.
Shunday qilib biz quyidagixossalarga egabo’lgan
Fn ( x)
Q
(i=1,2,…,n),
'
i ( x) =
Fn ( xi )( x - xi )
x
ko’phadni oldik . Qi ( x)
= x i nuqtadan tashqari barcha
(2.3)
n ta ikki hadliga
(2.4)
x = xk nuqtalarda
nolga teng,
ya’ni
x = xi da esa birga teng. Agar
(1, агар Qi ( x k ) = f ik =〈
l0 , агар
fik - Kroneker simvolini kiritsak,
(
i = k
i 丰 k
,
2.5)
ko’phad qo’yilgan shartni qanoatlantiradi: ya’ni x = xk nuqtalarda y
(k = 1, 2 ,..., n) qiymatlarni qabul qiladi .
Pn- 1 ( x) - ko’phadning yagonaligi shu dalildan kelib chiqadiki ,
y = y k
P
-
n 1 ( x)
ko’phad bilan ikkinchi gipotetik Pn- 1 ( x) ko’phad o’rtasidagi ayirma birga x = x k
nuqtalarda nolga aylanadi . Lekin Pn- 1 ( x) - Pn- 1 ( x) ayirma ham yana n - 1
darajali ko’phad bo’lib, u esa aynannolga aylanmasdan n - 1 tadantub ildizga ega
bo’lmaydi: bu esa Pn- 1 ( x) = Pn- 1 ( x)
ekanligini bildiradi.
Endi agar biz Pn- 1 ( x) ni y = f ( x) funksiyaga yetarlicha yaqinlashgan deb
hisoblasak,
A = ∫ Pn- 1 ( x)dx = Σ y k ∫ Q k ( x)dx , (2.7)
- 1 k = 1 - 1
hisoblasak, amaliyotda noma’lum f ( x) egrilik ostidagi yuzaga ega bo’lamiz. Berilgan ayrim taqsimlangan x = xk nuqtalar uchun Q k (x) ko’phadlar bir
qiymatli aniqlangan va shuninguchun ham
1
∫ Q k ( x)dx = 负k , (2.8)
- 1
aniq integrallar ba’zi bir sonli qiymatlarga ega bo’ladiki, qaysikim ular uchun
jadvallartuzish mumkin.
Bizni qiziqtiruvchi yuza uchun bu qiymatlar
funksiyaningtabiatiga bog’liq emas.
Oldingi x = x i nuqtalarni o’zgartirmasdan yangi
nuqtani qo’shamiz. Qo’shimcha x - xn+1 ikki hadnikiritib,
tamoman y = f ( x)
x = xn+1 qo’shimcha
Q n+1 (x) - qo’shimcha
ko’phadni xosil qilamiz. (2.4) ta’rifdan Qi ( x) uchun kelib chiqadiki, Q n+1 (x)
ko’phad Fn ( x) ko’phadga proporsionaldir, qaysikim ko’paytuvchi qisqarib ketadi. Xuddi shunday yangi
ko’paytiriladigan vaznli 负n+1 vaznli ko’paytuvchi
1
∫ Fn ( x)dx , - 1
aniq integralga proporsionaldir. Shunga o’xshash, agar yangi
x = xn+1 , xn+ 2 ,..., xn+ m
nuqtalarni ularni ordinatalari bilan kiritsak, u holdaularga mos
负n+1 , 负n+ 2 , ... , 负n+ m
1
负
vaznlar
n+ i = ∫ Fn ( x) ξm - 1 ( x) dx
integral bilan aniqlanadi, bu yerda ξm - 1 ( x) ayrim m - 1 darajali
ko’phadlardir. Ixtiyoriy ξm +1 ( x) ko’phad, x 0 , x 1 , x 2 , ..., xm - 1 darajali funksiyalarning chiziqli superpozitsiyasidan iborat ekanligidan, agar Fn ( x) quyidagi integrallarni qanoatlantirsa, bu hamma vaznlar avtomatik ravishda nolga
aylanadi.
1 1
∫ Fn ( x)dx = 0 , ..., ∫ Fn ( x) x m - 1 dx = 0 , (2.12)
- 1 - 1
haqiqatdan ham bizning talablarimiz m = n gacha borib,
1
∫ Fn ( x) xa dx = 0 (a = 0 ,1, 2 ,..., n - 1) , (2.13) - 1
integral shartiningbajarilishidir.
Natijada bizning boshida berilgan n ta nuqtani ixtiyoriy ravishda qo’shsak
ham baribir xech bir yangi ordinata oldinginatijalarni o’zgartirmaydi.
Oldingi natija shundan iboratki, xuddi biz 2 n ta ordinata bilan ish ko’rib, haqiqatdan esa biz n ta ordinatadan foydalanamiz, yangi qurilgan ordinatalar esa
hisoblanayotganyuzaga hech nima tushmaydi.
Bu jarayonda biz
A = y k 负k
yigindiga n ta hadni tejaymiz. Bu fikrlashlar yuqoridagi muloxazalar uchun
yetarlicha emasdir. To’liqroqbo’lishi uchun quyidagi muloxazanitavsiya etamiz.
Haqiqatdan ham yangi x n+1 , xn + 2 ,..., x2 n nuqtalarning berilishi nafaqat
Q n+ m ( x) (m = 1,2 ,..., n) yangi ko’phadlarni qo’shadi, xatto
oldingiQi ( x) (i = 1,2 ,..., n) ko’phadlar ham o’zgaradi: xar bir yangi x n+ m nuqta
x - xn+ m
Dostları ilə paylaş: |