15-ma’ruza chiziqli algebraik tenglamalar sistemasini taqribiy yechish usullari: Chiziqli algebraik tenglamalar sistemasini taqribiy yechishning oddiy iteratsiya va Zeydel usullari. Iteratsion jarayon yaqinlashishining zaruriy shartlari



Download 0,59 Mb.
bet2/4
Sana30.05.2023
Hajmi0,59 Mb.
#945754
1   2   3   4
Bog'liq
15- мавзу

ODDIY ITERATSION USUL

Faraz kilaylik, tizim biror usul bilan
Ax = b (2.1)


x + Cx + f (2.2)

ko`rinishga keltirilgan bo`lsin, bu erda S — qandaydir matritsa, f - vektor ustun. Dastlabki yaqinlashish vektori x(0) biror usul bilan (masalan, x(0) = 0) topilgan bo`lsin. Agar keyingi yaqinlashishlar
x(k+1) = Cx(k) + f (k=0,1,2, …)
rekkurent formula yordamida topilsa, bunday usul oddiy iteratsiya usuli deyiladi.
Agarda S matritsa elementlari
n

Cij i1
va
n
Cij i1
a  1
   1
i  1,2,..., n
j  1,2,..., n
(2.3)

(2.4)


shartlardan birortasini kanoatlantirsa, u xolda iteratsion jarayon berilgan tenglamaning x echimiga ixtiyoriy boshlangich x(0) vektorda yaqinlashishi isbotlangan, ya`ni

x  lim
k 
xk

Shunday kilib, tizimning aniq echimi cheksiz qadamlar natijasida -hosil qilinadi va hosil kilingan ketma-ketlikning ixtiyoriy vektori taqribiy echimni beradi. Bu taqribiy echimning xatoligini quyidagi formulalardan biri orqali ifodalash mumkin:

x x k | 
max
x k x k 1
(2.5)

i i

1
j1,2,... n j j

agarda (2.3) shart bajarilsa, yoki

x x k
n x k x k 1
(2.6)

i i 1 
j j j1

agarda (2.4) shart bajarilsa. Bu baxolarni moc ravishda quyidagicha kuchaytirish mumkin:

mxi x k  max x k x k 1



eki
i 1 i i


n x x k
n x k x k 1

i i j1
1 
j j j1

Iteratsion jarayonlarni yuqoridagi baxolar oldindan berilgan aniqlikni kanoatlantirganda tugallaydilar.


Boshlangich x(0) vektor, umuman olganda, ixtiyoriy tanlanishi mumkin. Ba`zan x(0) = f deb olishadi. Ammo x(0) vektorning komponentlari sifatida noma`lumlarning ko`pol taxminlarda aniqla-ngan qiymatlari olinadi.
(2.1) tizimni (2.2) ko`rinishga keltirishni bir necha xil usullarda amalga oshirish mumkin. Faqat (2.3) yoki (2.4) shartlardan birortasining bajarilishi lozim. Shunday usullardan ikkitasiga tuxtalamiz.
"Birinchi usul. Agarda A matritsaning diagonal elementlari noldan farqli bo`lsa, ya`ni

u xolda berilgan tizimni


aii 0 (I=1,2,…, n)



x1
1
1
a11
b1 a12 x2  ...  a1n xn

x2


a22
b2 a22 x1 a23 x3  ...  a2n xn
(2.7)





xn

... .................................
1 bn an1 x1  ...  an,n1 xn1
ann

ko`rinishda yozish mumkin. Bu xolda S matritsa elementlari quyida-gicha aniqlanadi:

Cij
aij
aii
i
j,
Cii  0

hamda (2.3)va (2.4)shartlar mos ravishda quyidagi ko`rinishni qabul kiladi:



n

j1
j i
a  1
i  1,2,..., n

(2.8)


n


i1
   1
j  1,2,..., n

(2.9)




(2.8)va (2.9)tengsizliklar A matritsaning diagonal elementlari

aii

  • aij j i

i  1,2,..., n
(2.10)

shartlartlarni kanoatlantirganda urinli bo`ladi.
Ikkinchi usul. Bu usulni quyidagi misol orqali namoyish qilamiz.
Umuman olganda, har qanday keltirilmagan matritsali tizim uchun yaqinlashuvchi iteratsion usullar mavjud, ammo ularning barchasi kisoblash uchun qulay emas.
Agarda iteratsiya usuli yaqinlashuvchi bo`lsa, u xolda bu usul yuko-rida kurilgan usullardan quyidagi afzalliklarga ega bo`ladi:

  1. Iteratsion jarayon tezrok yaqinlashsa, ya`ni tizimning echimini aniqlash uchun p dan kamrok iteratsiya talab kilinsa, u xolda vaktdan yutiladi, chunki arifmetik emallar soni p2 ga mutanosib (proportsional) (Gauss usuli uchun esa bu son p3 ga mutanosib).

  2. Yaxlitlash xatoliklari iteratsiya usulida natijaga kamrok ta`-sir etadi. Bundan tashqari iteratsiya usuli o`z xatoligini to`g’rilab boruvchi usuldir.

  3. Iteratsiya usuli tizimning muayyan koeffitsientlari nolga teng bo`lgan kolda juda ham qulaylashadi. Bunday tizimlar xususiy hosilali differentsial tenglamalarni echganda ko`prok uchraydi.

  4. Iteratsiya jarayonida bir xil turdagi amallar bajariladi, bu esa eX.M uchun programmalashtirishni osonlashtiradi.

1- misol. Quyidagi tizim oddiy iteratsiya usuli bilan yechilsin:

10x1 x2 3x3 2x4 x5 6




x1  25x2 x3  5x4  2x5  11
2x1 x2  20x3  2x4  3x5  19


x2 x3 10x4  5x5  10
x1  2x2 x3  2x4  20x5  32

E c h i s h . Birinchi usulda aytilganidek, bu tizimning tenglamalarini mos ravishda 10, 25, - 20, 10, 20 larga bo`lib, quyidagi ko`rinishda yozib olamiz:



x1  0,6
 0,1x2
 0,3x3  0,2x4
 0,1x5

x2
 0,44  0,04x1  0,04x3  0,2x4  0,08x5


x3
 0,95  0,1x1  0,05x2  0,1x4  0,15x5


4


x  1 0,1x2  0,1x3  0,5x5

x5  1,6  0,05x1  0,1x2  0,05x3
 0,1x4

bu erda (2.8)shart bajariladi. Xaqiqatdan ham,





5
C1 j j1
5
C3 j j1
5
C5 j j1
 0,3  1;


 0,41  1;


 0,3  1;
5
C2 j j1
5
C4 j j1
 0,28  1;


 0,5  1;

Dastlabki yaqinlashish x(0) sifatida ozod xadlar ustuni (0,6; 0,44; 0,95; 1; 1,6) ni olib keyingi yaqinlashishlarni topamiz:
x1  0,6  0,1x0  0,3x0  0,2x0  0,1x0 =
1 2 3 4 5
0,6 – 0,1  0,44 + 0,3  0,95 + 0,2  1 – 0,1  1,6 = 0.881

2
х 1= 0,44 + 0.04  0,6 – 0,04  0,95 + 0,2  1 + 0,08  1,6 = 0,754

Shunga o`xshash
х 1= 0,892;
х 1= 1,851;
х 1= 1,72. Hisoblashlarning davomini


3

4

5
1- jadvalda keltiramiz:
1-jadval


k

х k
1

х k
2

х k
3

х k
4

х k
5

0

0,6

0,44

0,95

1

1,6

1

0.881

0,754

0.892

1,851

1,72

2

0.9884

0.9482

1,0029

1,9147

1,9859

3

0,9904

0,9814

0,9908

1,9939

1,9854

4

0,99944

0.99753

0,99789

1,99364

1.99897

5

0,99839

0,99865

0,99929

1,99954

1,99970

6

0.99986

0,99989

0,99977

1,99976

1.99960

7

0,999934

0,999920

1,000018

1,999788

1,999947

8

0.999974

0,999951

0,999976

2,000042

1,999978

Yuqoridagi 3.4- jadvaldan ko`ramizki, 8-iteratsiya x1= 0,999974; x2= 0,99951; x3= 0,99998; x4 = 2,00004; x5= 1,99998 echimdan iborat. Bu topilgan taqribiy echim aniq echim
x1* = x2* = x3* = 1; x4* = x5* = 2 dan beshinchi xonaning birliklari buyichagina farqlanadi.

Download 0,59 Mb.

Do'stlaringiz bilan baham:
1   2   3   4




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

kiriting | ro'yxatdan o'tish
    Bosh sahifa
юртда тантана
Боғда битган
Бугун юртда
Эшитганлар жилманглар
Эшитмадим деманглар
битган бодомлар
Yangiariq tumani
qitish marakazi
Raqamli texnologiyalar
ilishida muhokamadan
tasdiqqa tavsiya
tavsiya etilgan
iqtisodiyot kafedrasi
steiermarkischen landesregierung
asarlaringizni yuboring
o'zingizning asarlaringizni
Iltimos faqat
faqat o'zingizning
steierm rkischen
landesregierung fachabteilung
rkischen landesregierung
hamshira loyihasi
loyihasi mavsum
faolyatining oqibatlari
asosiy adabiyotlar
fakulteti ahborot
ahborot havfsizligi
havfsizligi kafedrasi
fanidan bo’yicha
fakulteti iqtisodiyot
boshqaruv fakulteti
chiqarishda boshqaruv
ishlab chiqarishda
iqtisodiyot fakultet
multiservis tarmoqlari
fanidan asosiy
Uzbek fanidan
mavzulari potok
asosidagi multiservis
'aliyyil a'ziym
billahil 'aliyyil
illaa billahil
quvvata illaa
falah' deganida
Kompyuter savodxonligi
bo’yicha mustaqil
'alal falah'
Hayya 'alal
'alas soloh
Hayya 'alas
mavsum boyicha


yuklab olish