0
1
2
3
4
5
6
7
|
–2,00
–2,00
–1,75
–1,75
–1,75
–1,75
–1,75
–1,74
|
–1,00
–1,50
–1,50
–1,63
–1,69
–1,72
–1,73
–1,73
|
–1,50
–1,75
–1,63
–1,69
–1,72
–1,73
–1,74
|
–3,5625
0,3633
–1,8140
–0,7981
–0,2363
–0,0406
0,1592
|
Javob: x1 ≈ –1,73. Ikkinchi ildizni ham xuddi shunday topish mumkin.
misol. Kesmani teng ikkiga bo‘lish usulidan foydalanib, x3+3x2–3=0 tenglamaning [–3;–2] kesmadagi ildizini ε = 0,1 aniqlik bilan hisoblang.
Yechish. Yuqorida keltirilgan algoritga asoslanib, tenglamani yechish jarayonini quyidagi hisob jadvali ko‘rinishida yozamiz:
n
|
an
|
bn
|
f(an)
|
f(bn)
|
xn
|
f(xn)
|
(bn–an)/2
|
0
|
–3
|
–2
|
–3
|
1
|
–2,5
|
0,125
|
0,5
|
1
|
–3
|
–2,5
|
–3
|
0,125
|
–2,75
|
–1,11
|
0,25
|
2
|
–2,75
|
–2,5
|
–1,11
|
0,125
|
–2,625
|
–0,42
|
0,125
|
3
|
–2,625
|
–2,5
|
–0,42
|
0,125
|
–2,5625
|
–0,129
|
0,0625
|
Jadvalga ko‘ra ildiz x = –2,5625 0,0625 yoki natijani yaxlitlasak, u holda x = –2,6 0,1.
Mashqlar
Quyida berilgan f(x)=0 ko‘rinishdagi tenglamalarni kesmani teng
ikkiga bo‘lish usuli bilan yeching (bunda a, b, c, d, parametrlarni o‘zingiz har xil tanlash orqali turli variantlar hosil qilishing‘iz mumkin):
1. a /b x2 cx 0 , a = 1.05; b = 0.1; c = 2.03; = 10-3.
2. a x2 bx2 0, a =1.23; b = –3.14; = 410-5.
3. bx a3 cx d 0 , a = 0.1; b = 2.23; c = 2; d = –1.03; = 210-4.
4. (x a)5 bx 0 , a =0.29; b = 2; = 310-4.
(x a)2 ebx
0 , a = –0.4; b = 0.53; = 10 -4.
a /
cosbx 0 , a =2.07; b =1.19; c =1.13; =210 -5.
Izoh: Dastlab f( x) funksiyaning grafigini matematik paketlardan birida (Maple, Matlab, Mathcad, Mathematica va boshqa) yoki MS Excel dasturida chizing, haqiqiy ildizlar yotgan oraliqlarni aniqlab oling, hisoblashlarni qo‘lda va dastur yordamida bajaring.
Vatarlar usuli (proporsional bo‘laklar usuli) Usulning mazmuni.
Quyidagi shartlarning bajarilishini talab qilamiz:
f(x) funksiya o‘zining f ( x ) va f ( x ) hosilalari bilan [a,b] kesmada uzluksiz;
funksiyaning f(a) va f(b) qiymatlari kesmaning oxirgi nuqtalarida har xil ishorali, ya’ni f(a) · f(b) < 0;
har ikkala f ( x ) va f ( x ) hosilalar [a,b] kesmaning barcha nuqtalari- da o‘z ishorasini saqlab qoladi (berilgan [a,b] kesma f(x) funksiya hosilasining o‘z ishorasini saqlashi bu shu funksiya monotonligining yetarli sharti).
Bulardan foydalanib, 1.14, a-rasmga asosan, dastlab A( a, f( a)) nuqta qo‘zg‘almas, x0=b – nolinchi yaqinlashish, A( a, f( a)) va B( b, f( b)) nuqtalarni tutashtiruvchi AB vatarning absissa o‘qi bilan kesishish nuqtasini x1 – birinchi yaqilashish deb qabul qilamiz.
Keyingi yaqinlashishlarni hisoblash uchun f(x1) qiymatni hisoblaymiz va uni f(a) va f(b) qiymatlar bilan taqqoslaymiz. Hosil bo‘lgan [a, x1] va [x1, b] intervallardan chetlarida f(x) funksiya har xil ishorali bo‘lganini tanlaymiz, chunki aynan ana shu intervalda izlanayotgan x ildiz yotadi. Yuqorida aytilgan uslubni ana shu intervalga qo‘llab, keyingi yaqinlashishni (x2 nuqtani) topamiz. Keyingi yaqinlashishlarda
funksiyaning birinchi va ikkinchi hosilalari intervallarda o‘z ishorasini saqlaydi, deb o‘suvchi ketma-ketlikni tashkil etuvchi va yuqoridan x qiymat bilan chegaralangan barcha x1, x2, ... yaqinlashishlarni topamiz.
Natijada
lim x
n
n
x . Ketma-ket yaqinlashishning formulasini chiqarish
uchun xn dan xn+1 ga o‘tishni qaraylik. Bu holda Bn va B nuqtalardan o‘tuvchi BnB vatar tenglamasini tuzamiz:
y f (xn )
x xn .
f (a) f (xn )
a xn
Agar bu tenglamada y( xn+1) = 0 desak, u holda undan xn+1 had topilad.
Bularga asosan umumlashgan quyidagi to‘rtta holat bo‘ladi:
Agar [a,b] kesmada A(a,f(a)) nuqta qo‘zg‘almas va f(x) funksiya
botiq va kamayuvchi ( f (x) 0, f (x) 0 ) bo‘lsa (1.14, a-rasm), u holda ushbu ketma-ketlik
x0=b ,
n1
n f ( x )
f ( a) n
(n=0,1,2,…) (5.2)
chegaralangan va monoton kamayuvchi: a< x <…< xn+1< xn<… x1< x0= b.
n
Agar [a,b] kesmada B(b,f(b)) nuqta qo‘zg‘almas va f(x) funksiya
botiq va o‘suvchi ( f (x) 0, f (x) 0 ) bo‘lsa (1.14, b-rasm), u holda ushbu ketma-ketlik
x0=a ,
n 1
n f (b)
f (xn )
n (n=0,1,2,…) (5.1)
chegaralangan va monoton o‘suvchi: a= x0< x1<…< xn< xn+1<…< x < b.
Agar [a,b] kesmada A(a,f(a)) nuqta qo‘zg‘almas va f(x) funksiya
qovariq va o‘suvchi ( f (x) 0, f (x) 0 ) bo‘lsa (1.14, a-rasm), u holda ushbu ketma-ketlik
x0=b ,
n1
n f ( x )
f ( a) n
(n=0,1,2,…) (5.2)
chegaralangan va monoton kamayuvchi: a< x <…< xn+1< xn<… x1< x0= b.
n
Agar [a,b] kesmada B(b,f(b)) nuqta qo‘zg‘almas va f(x) funksiya
qovariq va kamayuvchi ( f (x) 0, f (x) 0 ) bo‘lsa (1.14, b-rasm), u holda ushbu ketma-ketlik
x0=a ,
n 1
n f (b)
f (xn )
n (n=0,1,2,…) (5.1)
chegaralangan va monoton o‘suvchi: a= x0< x1<…<xn<xn+1<…< x <b.
1.14-rasm. Proporsional bo‘laklar usuli (vatarlar usuli)ning har xil hollari uchun sxemalar.
Endi bu to‘rtta holatni umumlashtiramiz:
agar kesmaning qaysi bir chetida f(x) funksiya va uning f(x) ikkinchi hosilasi bir xil ishoraga ega bo‘lsa, o‘sha chetki nuqta qo‘zg‘almas deb olinadi;
agar x ildizning qaysi tarafida f(x) funksiya o‘zining f(x) ikkinchi hosilasiga qarama qarshi ishoraga ega bo‘lsa, xn ketma-ket yaqinlashishlar o‘sha tomondan x ildizga yaqinlashadi.
Iteratsion jarayonning tugallanishi ikkita qo‘shni xn va xn-1 iteratsiyalarning hisob hatijalari bo‘yicha hisobni tugallash kriteriyasini beradi, ya’ni bu iteratsion jarayon ushbu
xn+1–xn<
shart bajarilgunga qadar davom ettiriladi va xn+1= x yoki xn = x yechim deb olinadi, bu yerda – berilgan limitik absolyut xato.
Usulning qulayliklari: usulning yaqinlashuvchanligi kafolatlangan; oraliqni teng ikkiga bo‘lish usuliga qaraganda kamida ikki yoki uch marta tezroq yaqinlashishni beradi;
Usulning kamchiliklari: agar a dan b gacha bo‘lgan kesmada umuman ildiz mavjud bo‘lmasa yoki unda bir nechta ildizlar mavjud bo‘lsa, u yechimni izlash vaqti cheksizga yaqinlashishi mumkin; agar f(x) funksiya grafigi [a,b] kesmada yetarlicha yotiq bo‘lsa, u holda f(a) – f(b) farq katta bo‘ladi va hisoblashlarda xatolik ko‘payadi, bunday holda keyingi hisoblashlarda oraliqni teng ikkiga bo‘lish usuliga o‘tgan ma’qul.
Usulning hisob algoritmi:
[a,b] kesmani va aniqlikni berish.
Agar f(a) va f(b) lar bir xil ishorali yoki f '(a) va f '(b) lar har xil ishorali bo‘lsa, ildizni topish mumkin emasligini bildirish.
Boshlang‘ich yaqinlashishni va navbatdagi yaqinlashishning iteratsion hisob formulasini yuqoridagi to‘rtta holatdan biri bo‘yicha tanlash.
Hisoblashlarni tanlangan iteratsion hisob formulasida bajarish.
Aniqlikni baholash:
xn1 xn
.
Agar bu shart bajarilsa, ildiz deb x =xn+1 ni qabul qilish, aks holda 4- qadamga o‘tish.
Usulning blok-sxemasi 1.15-rasmda tasvirlangan.
Dasturda cheksiz takrorlanishlar ku-zatilmasligi uchun funksiya grafigining qovariq yoki botiqligini (1.16-rasm) va iteratsiyalar sonini nazorat qilish maqsadga muvofiq.
Shunday qilib, vatarlar usulidan foydalanishda ushbu qoidaga amal qilish maqsadga muvofiq: kesmaning qaysi chetida funksiyaning ishorasi uning ikkinchi tartibli hosilasi ishorasi bilan bir xil bo‘lsa, o‘sha chet qo‘zg‘almas qilib olinadi.
1-misol. Ushbu qoidani (x 1) ln(x) 1 0 tenglamaning [2;3] kesmadagi izolyatsiyalangan ildizini topishga qo‘llang.
Yechish. Bu yerda
f(x)= (x 1) ln(x) 1;
f (x) ln(x) x 1 ; f (x) 1 1 .
x x x 2
Kesmaning b = 3 chetki nuqtasida funksiyaning qiymati musbat f(3)>0 va kesmada esa ikkinchi tartibli hosila ham musbat f (x) 0 ,ya’ni f (b) f (x) 0 . Shunday qilib, ushbu misolda berilgan tenglamaning izolyatsiyalangan ildizini vatarlar usuli bilan topish uchun ushbu
x x f (xn )(b xn )
n1 n f (b) f (x )
n
formuladan foydalanish tavsiya etiladi.
|
1.15.-rasm. Vatarlar usulining blok-sxemasi.
|
a)
b)
1.16-rasm. Vatarlar usulining geometrik interpretatsiyasi:
a) f (x) f (x) 0 ; b) f (x) f (x) 0 .
2-misol. Ushbu f(x) = x3 – 0.2x2 – 0.2x – 1.2 = 0 tenglamaning [1;1.5] kesmadagi x ildizini 0.002 aniqlik bilan hisoblang.
Yechish. Berilgan tenglama [1;1.5] kesmada yagona x ildizga ega (buni talabaning o‘zi Maple matematik paketida yoki MS Excel dasturida mustaqil tekshirib ko‘rsin). Usulning formulasiga ko‘ra
f '(x) = 3x2 – 0.4x – 0.2 ; f ''(x) = 6x – 0.4 ;
f(1) = –0.6 < 0 ; f(1.5) = 1.425 > 0,
demak x[1;1.5] da 2.4 f '(x) 5.95 ; 5.6 f ''(x) 8.6 .
Bu yerdan ko‘rinadiki, x[1;1.5] da f '(x) f ''(x) > 0. Shuning uchun
x0=a ,
n 1
n f (b)
f (xn )
n (n=0,1,2,…)
formuladan foydalanib (bunda x0 = a = 1, b = 1.5), ketma-ket quyidagilarni topamiz:
x 0.6(1.5 1) 1 1.15;
1 1.425 0.6
f ( x1
) 0.173;
x 1.15 0.173
1.5 1.15
1.19;
f (x
) 0.036;
2 1.425 0.173 2
x3 1.19 0.036
1.5 1.19
1.425 0.036
1.198;
f ( x3
) 0.0172;
x 1.198 0.0172
1.5 1.198
1.199;
f (x
) 0.0061;
4 1.425 0.0172 3
Bu yerdan ko‘rinadiki,
x4 x3
1.199 1.198 0.002
ekanligidan taqribiy
ildiz sifatida
x x4 1.199 ni qabul qilishimiz mumkin (aniq ildiz x =1.2).
Do'stlaringiz bilan baham: |