xn(n= 1,2,3...) sonlarini (2) formula yordamida ketma-ket hisoblash jarayoni iteratsiya metodi deyiladi. 9-rasmda iteratsiya metodining geometrik interpretatsiyasi keltirilgan.
Agar [a, p] oraliqda |
shart bajarilsa, u holda jarayon yaqinlashuvchi bo‘ladi.
9-rasm. Iteratsiya metodining geometrik та 'nosi.
Iteratsiya jarayonini ikki qo'shni bo'lgan ketma-ket yaqin-
I
lashishlar: xn_t va xn uchun quyidagi tengsizlik
Я
bajarilmaguncha davom ettiramiz. Bu yerda £ - berilgan x’ ildizning absolyut xatoligi.
Agar q<0,5 bo‘lsa, u holda |л:я — jc„_1 | <£ oddiy ifodadan foyda- lansa bo‘ladi. Agar f'{x)< 0, u holda /(x) = 0 tenglama o‘rniga -fix) = 0 tenglama ko‘riladi.
(1) ni unga ekvivalent bo‘lgan tenglama x = x- A f{x) Я >- 0 bilan almashtiriladi. я parametri shunday tanlab olinadiki, bunda 0 < cp'{x) = 1 - A/'(*) < q < 1 shart bajarilishi kerak.
Misol:
[—0,5;0] oraliqda joylashgan (3) tenglamaning ildizini £ = 10"
absolyut xatolik bilan iteratsiya metodi orqali toping. Bundan tashqa- ri, ildizni hisoblashga kerak bo‘ladigan qadamlar sonini ham toping.
Yechish:
(3) tenglamani x = cp(x) ko‘rinishga keltiramiz:
arcsin(2x +1) = x2
Sin(arcsm(2x +1)) = Sinx2
2x +1 = Sinx2
x = 0.5{Sinx2 -1)
Bu yerdan q>(x) = 0,5{Sinx2 -1) (4)
q>'(x) = xcosx2 topamiz. Ko'rinib turibdiki, hamma 0,5 0 lar uchun. Shuning uchun q=0,5 va iteratsiya jarayoni yaqinlashuvchi boiadi.
Algoritmi:
Boshi
X=-0.4
e =0,0001
n=i
|
|
|
Л
|
f
|
;
|
Y=0t5(Sinx2-l)
S=\y-x\
|
I
I
|
X
|
= V
n+1
|
Boshlang‘ich yaqinlashuv sifatida [—0,5;0] oraliqning xohlagan nuqtasini tanlab olish mumkin, masalan: x0=-0,4. (3) tenglamani iteratsiya metodi orqali yechishning blok-sxemasi 8-rasmda ko'rsatilgan.
x= -0,4, s = 1*1(И,и = 0
XYo'q
Bosmaga
x,n
Keyingi yaqinlashish xn+| ni x„+1 = 0,5(sinx2-1) formula orqali hi- soblanadi.
5 = xn+x-xn farqni hisoblab, n bittaga oshiriladi.
|5| > £ shart tekshiriladi. Agar u
Tamom
10-rasm. Iteratsiya metodi blok-sxemasi.
bajarilsa, keyingi yaqinlashish hisobla- nadi. Agar u bajarilmasa, u holda xn+| ni ildiz deb va hisoblash tugatiladi, bunda n bajarilgan iteratsiyalar soniga teng bo'ladi.
Dasturi:
10 X=-0,4 N=0 E=lE-4 20 Y=0,5*(SIN(XA2)-1)
30 D=ABS(Y-X) X=Y N=N+1
40 IF D>E GOTO 20
50 PRINT «Х=»; X, «N=»; N
60 STOP
RUN
Javobi:
X=-0.414539 N=7
Nyuton metodi
[a,p] oraliqda f(x) = 0 tenglama bitta ildizga ega, hamda f'(x)
va f"{x) ushbu oraliqda aniqlangan va uzluksiz bo‘lsin.
Nyuton metodining geometrik ma’nosi 11-rasmda keltirilgan.
[a,P] oraliqdan ixtiyoriy x0 nuqta olinadi va y=f(x) egri chizig‘ini
P0(x0,f(x0)) nuqtasidan OX o‘qi bilan kesishguncha urinma o‘tkazila- di. Abssissa bilan kesishish nuqtasi x, ni ildizga yaqinlashishni birinchi nuqtasi deb olish mumkin. Yangi Р^Дх,)) nuqtadan urinma o‘tka- zib, abssissa bilan kesishish nuqtasini - X2 ildizga yaqinlashishning ikkinchi nuqtasi sifatida olinadi. Keyingi yaqinlashishlarni shu tariqa topish mumkin.
Ildizga ketma-ket yaqinlashish formulasini keltirib chiqaramiz. PO nuqtadan o‘tuvchi urinma tenglamasi quyidagi ko‘rinishda bo‘ladi:
y=f(x0)+f(x0)(x-x0) (5)
y=0 deb, x, abssissa nuqtasidan o‘tkazilgan urinmani OX o‘qi bilan kesishish nuqtasini topamiz:
fixo)
'■'■/hi (6) Keyingi yaqinlashishlar mos ravishda quyidagi formulalardan topiladi:
fix,)
X, = X, --
fix,)
(7)
Yaqinlashish jarayonini hisoblashni quyidagi shart bajarilsa to‘xtatamiz:
2m. £
-ПГ (8)
Bu yerda m, - |/'0)| ning [a,/3] oraliqdagi eng kichik qiymati,
M2 - |/"(x)| ning [a,p] oraliqdagi eng ka tta qiymati.
Bu yerda \x*-xn\<£ tengsizlik o‘rinli bo‘ladi, bu yerda £ - x* ildizning berilgan chegaraviy absolyut xatoligi. Agar M2 ning qiymati
2w, 7 I—
0 bo‘lsa, u holda (8) tengsizlikni |*я -x„_,| < 10 ’V£ ko‘rinish- da yozish mumkin
Tamom
Bosmaga x,n
T
12-rasm. Nyuton metodining blok- sxemasi.
Misol uchun £=10-6 bo‘lsa, u holda quyidagi oddiy shartdan foydalanish
munkin: |jc„ — | с 10-4
xQ boshlang‘ich yaqinlashishni shun- day tanlab olish maqsadga muvofiqki, bunda quyidagi (9) shart bajarilishi kerak:
/(*o)/"W>0 (9)
Aks holda metodning yaqinlashishi kafolatlanmaydi.
Ko'pincha (9) shart bajarilishiga qarab x0 =a yoki x0 = p deb tanlab olinadi. Nyuton metodi f(x)=0 tengla- malarni yechishda f(x) funksiyasi grafigi berilgan ildizning atrofida katta egrilikka ega bo‘lgan hoi uchun yaxshi samara beradi.
Misol:
f(x)=sinx-x+0.15=0 tenglamaning oraliqdagi ildizini absolyut xatolik bilan Nyuton metodi yordamida toping. Bundan tashqari, ildizni hisoblashga kerak bo'ladigan qadamlar sonini ham toping.
Yechish:
f'(x) = cosx-1 topamiz.
Nyuton metodining hisoblash formulasi quyidagi ko'rinishga ega bo‘ladi:
sinx.
' xn-\ +0,15
cos X , -1
, bu yerda n= 1,2,3...
Boshlang‘ich yaqinlashish nuqtasi qilib X0=l tanlab olamiz. Bu nuqtada yaqinlashish sharti f(x0)f"(x0) > 0 bajarilyapti. 12-rasmda ushbu misolni yechishning blok-sxemasi keltirilgan.
Dostları ilə paylaş: |