2w, 7 I— 63
1лА=т1 70
+t(e -*A, 98
* = 1
ega bo‘lamiz.
Det [*,*"'J*0 (Vandermond aniqlovchisi), demak, bundan (25)
tizim yagona yechimga ega.
Quyidagi ko‘rinishda yozilgan interpolyatsiya ko‘phadi:
gn « = к (x) = £ fix, )Пт7^ (26)
/=1 j*i X, Xj
Lagranj interpolyatsiya polinomi deyiladi.
n
a>„(x) = ^(x-xJ) belgilashlami kiritamiz. U holda:
>i
м=X П (* - xj y&n (xi)=П (x, - xi)
£=1 j*k j*i
Demak, (26) ni quyidagi ko‘rinishda yozish mumkin:
Lagranj interpolyatsiya polinomi qoldiq hadini baholash
x) = /(*) -g„(x)~K(0n (x)
deb olamiz.
К ni
shartdan tanlaymiz, bu yerda x* - xatolikni baholovchi nuqta. Unda:
f(x*)-g„(x*)
K = -
«„(**)
К ni bunday tanlashda
xl,x2,...,xn,x* nolga aylanadi. Roll teoremasiga ko‘ra uning hosilasi q> (x) bo‘lmaganda bitta nuqtada nolga aylanadi: \yx,y2\\yx=min{xx,x2,...,xn,x*y, уг=т&х(хх,х2,...,хп,х*). (P(n)(x) = f{n)(x) — nK, bo‘lgani uchun
M{£,) = 0, bundan K= f(n)(%)/n\
«.[*•*]- (28)
n\
(26) Lagranj interpolyatsiya formulasining qoldiq hadi.
Eytken hisoblash sxemasi. Eytken sxemasi (24) polinom a:, / = 1, n,
koeffitsientlarini (26) hisobga olib, hisoblashni ikkinchi darajali funksional aniqlovchilarga keltirishga yordam beradi.
Belgilashlar kiritamiz:
,i = \,n-\\ (28)
F
va h.k.
(30)
U holda Px,2,...,n(x), funksional aniqlovchini hisoblab, (26) ni (24) ko‘rinishida hosil qilamiz, ya’ni а,,а2,...,<з„ koeffltsientlaming qiymatini hisoblaymiz.
Ayniqsa, fiksirlangan £ e[x,.,x(+l] nuqtalarda Lagranj polinomini
hisoblashda bu sxema katta iqtisod beradi.
(28) aniqlovchilarni hisoblash uchun quyidagiga ega bo‘lamiz:
pu+1 W = au+ a2.ix’ i = 1’" -1 , (31)
Г, = 1 /(*i+, - x,); au = r (xMyt - xtyM); a2i = r, (yM - y.t),
bkJ,k = 1,3 koeffitsientlarni quyidagi formulalar orqali hisoblaymiz:
r, =1 - x,~\); К = n (auxi+i - ei.i+ixi);
^2,i = ri(a2,i-lXi+l ~a2,iXi-I + a2,i+l — °l,i+l)>
К,=гЛа2- «2./-1); (32)
^-U.U W = + *2./^ + йз,,* ; ' = 2, П -1.
akj = bk,i> к = 1,3 qayta belgilab, (32) ko'rinishdagi formuladan umumiy ko‘rinishdagi (30) polinom koeffitsientlarini hisoblash mumkin:
(33)
r, = 1 /(x,.+t+1 - xf); = r, {auxi+M - a, ,.+Ix,.);
bj+H = ri(ai+uxn
x-2
|
4
|
1 1
|
x-5
|
3,5
|
|
3,5
|
= /,45(JC)=7 0 1
|
x - 6
|
x-5
|
6
|
P„W = t
Ых) = -
x +1 4x +1 x — 2 l,5x + l
= ^(6 + 19x-5x2);
Aj+\.iAi+k*\ aj*\.i+\Xi aj,i + aj+\,i+\ )>
к = 1, n - 2; i = 1 ,k +1; j = \,n-k-\.
Misol:
XI
|
-1
|
0
|
2
|
5
|
6
|
У1
|
-3
|
1
|
4
|
3,5
|
6
|
1
|
x+1 -3
|
1
|
x 1
|
|
№) = y
|
x 1
|
= l + 4x; P2}(x) = -
|
x-2 4
|
= 1 + l,5x;
|
-9 + 2,5x;
= i(6 + 13x-2x2);
= — (66 - 29x + 4x2); 6
x 1,5л: +1 x-5 1/6 (26-x)
x-2 1/6(26 — jc)
X —6 2,5x-9
^.2.3.4 М — '
х + 1 1/6(6 + \9х-5х2) х —5 1/6(6 +13х-2х2)
= — (12 + 36х-11х 2+х 3); 12
х 1/6(6 + 13х-2х2) х — 6 1/6(66-29х + 4х2)
= — (6 + 23х-9х2 + х3);
6
х + 1 1/12 (12 + 36х — 11х 2 +х 3) х-6 1/6(6+ 23х-9х 2 +х 3) ' Shunday qilib:
4 (*) = ^1.2.3,4.5 W = 1 + 3,119х - 0,881х2 + 0,0119х3 (1 + х). Hisoblashlarni tekshirish uchun topamiz:
L«(0) = 1; L 4(2) = 4; L 4(3) = 3,5.
Dostları ilə paylaş: |