2021 Mundarija


Download 0.78 Mb.
bet5/10
Sana05.01.2023
Hajmi0.78 Mb.
#1079323
1   2   3   4   5   6   7   8   9   10
Bog'liq
327798 (1)

ayirmali sxema qurish. Faraz qilaylik , G sohakvadrat bo’lib , shu sohada ushbu
Δu = + = f (x1 , x2 ) (1.2.17)
Puasson tenglamasiuchun ayirmali sxema qurishtalab qilinsin.



3















2



4 1
1.2-chizma. Tugunnuqtalar.
Bunday sxemani ikki xil to’r ustida bajaramiz. Avvalo, qadami h ga teng bo’lgan kvadrat to’rni qaraymiz, 1.2-chizmada ko’rsatilganidek, 0 tugun atrofida 1,2,3,4 bilan belgilangan tugunlarni olamiz. Bu yerda x1 va x2 teng huquqli bo’lganligi hamda tugunlar simmetrik ravishda joylashganligi sababli
ayirmali approksimatsiyani quyidagiko’rinishdaizlash mumkin:
Lhu0 = c0 u0 + c1 (u1 + u2 + u3 + u4 ). (1.2.18)


Lhu0 = c0 u0 + 4c1u0 + c1 h - u + - 2 u + - 3 u +
Qaralayotgan sohada (1.2.17) tenglamaning yechimi to’rtinchi tartibli uzluksiz
hosilalarga ega deb faraz qilib (1.2.18) ifoda uchun quyidagiga egabo’lamiz:
h + u + + 2 u - + 3 u 0 +
+ + 4 u (ξ1 ,η1 ) +
| 0 | + h 4 δ - δ )||4 u + 4! |L(δx1 δx2 ) 」|(ξ2 ,η2 ) + h(||(- + u + - + 2 u + - + 3 u 0
Bunda
x10 - h < ξj < x10 + h, x20 - h j < x20 + h, j = 1,2,3,4.
Bu ifodani soddalashtirib quyidagini olamiz:
Lhu0 = (c0 + 4c1 )u0 + 4c1 + 0 + R(h), (1.2.19)
Bunda R(h)- qoldiq had. (1.2.19) ifoda (1.2.17) tenglamani approksimatsiya qilishi
uchun
с0 + 4с1 = 0,2c1h2 = 1
shartlar bajarilishikerak. Bulardan esa
c0 = - , c1 =
h 2h
kelib chiqadi. Shunday qilib, natijada quyidagiga egabo’ldik:
Lhu0 = (u1 + u2 + u3 + u4 - 4u0 ) = (Δu)0 + R(h) (1.2.20)
Agar M4 orqali tortinchi hosilalarning G dagi maksimumi modulini
belgilasak , u holda R(h) qoldiq had uchun ushbu bahoga egabo’lamiz:
R(h) < 4 c 1 24 M4 = M4 (1.2.21)
Yuqoridagi (1.2.20) ifodada (Δu)0 ni f 0 orqali almashtirib , R(h) qoldiq
hadnitashlabyuborsak , natijada uj ning to’rdagitaqribiy qiymati yj uchun ushbu
y1 + y2 + y3 + y4 - 4y0 = 2h2 f0
ayirmali tenglamaga ega bo’lamiz. Bundan approksimatsiyaning xatoligi h 2 M4
dan oshmaydi.
Izoh. (1.2.10) va (1.2.21) baholarni solishtirish shuni ko’rsatadiki , (1.2.10)
baho (1.2.21) ganisbatan 8 marta kichik.
Ayirmali sxemaning turg’unligi. Biz yozuvni qisqaroq qilish maqsadida ayirmali sxemaning turg’unligini Puasson tenglamasi uchun Dirixle masalasini
yechishmisolida ko’rib chiqamiz.
Faraz qilaylik, G soha to’g’ri burchakli to’rtburchak G {0 < x1 < a,0 < x2 < b} bo’lib , Γ uning chegarasi bo’lsin. Shunday u(x1 , x2 )
funksiyanitopishkerakki , u G da
+ = f (x1 , x2 ) (1.2.22)
Tenglamani qanoatlantirib, Γ chegarada Dirixle shartini qanoatlantirsin:
u Γ = Q(x1 , x2 ) (1.2.23)
Bunda Q(x1 , x2 ) ma’lum funksiya. Faraz qilaylik, (1.2.22)- (1.2.23) chegaraviy

masala G = GUГ sohadayagona yechimga egava buyechim G da 4 va 4
δ 4 u δ 4 u
δx1 δx2
uzluksiz hosilalarga egabo’lsin.
Biz quyidagi to’g’ri burchakli tortburchaklardan iborat bo’lgan to’rni
qaraymiz:
x
.x2k = kh2 , k = 0,1,...., N, Nh2 = b.J
1i = ih1 , i = 0,1,....., M , Mh1 = a; ) (1.2.24)


Endi G da yotuvchibarcha tugunlarni Gh deb olib chegaraviy nuqtalar Γh
sifatida Γ dayotuvchitugunlarni olamiz. Keyin
δ 2 u δ 2 u
Δu = 2 + 2
δx1 δx2
Laplas operatorini Gh ga tegishli nuqtalarda 3-chizmadagi besh nuqtali andaza



yordamida
Δ hyik = + = fik i = 1,2,...., M - 1, k = 1,2,......, N - 1
Ayirmali sxema bilan approksimatsiya qilamiz. Agar
bo’lsa,u holda (1.2.26) approksimatsiyaningxatoligi
Rik (h) < 12 (1 + c2 )M4
daniborat. (1.2.23) shartni quyidagilarga almashtiramiz:
yik Γh = Q(ih1 , kh2 ), (ih1 , kh2 )e Γh


(1.2.25)


h1 = h, h2 = ch, (c= const)


(1.2.26)

Qaralayotgan soha to’g’ri burchakli to’rtburchak bo’lganligi tufayli (1.2.26)approksimatsiyaning xatoligi nolga teng. Chegaradagi (1.2.26) qiymatlar ma’lum bo’lganligi uchunularni (1.2.25) tenglamaga qo’yib, keyin ma’lum hadlarni o’ng tomonga o’tqazib, quyidagi chiziqli algebraik tenglamalar
sistemasiga egabo’lamiz:
Lhyik = fik , i = 1,2,...., M - 1; k = 1,2,....., N - 1. (1.2.27)
Ravshanki , (1.2.27) tenglamalar faqat chegara yaqinidagi tugunlarda (1.2.25) tenglamalardan farq qiladi. Masalan , (i,1) ko’rinishdagi tugunlarda
(1.2.27) tenglama quyidagiko’rinishga egabo’ladi:
+ = fi1 - =Ψi1
(1.2.27) sistemada tenglamalarning soni noma’lumlarning soniga teng. Shuning uchun ham (1.2.25) sistemaning matritsasini Gh to’r ustidagi funksiyani o’ziga
akslantiradigan chiziqli operatordek qarash mumkin.

Endi (1.2.25), (1.2.26) tenglamalar
mavjudliginiko’rsatamiz.
1-lemma. Faraz qilaylik , 9(h) = {9ik }
aniqlangan qandaydir funksiya bo’lsin. Agar


sistemasining yagona yechimi

miqdorlar Gh = GhUΓh to’r ustida
Gh sohaning tugunlarida Δ h9(h) 0

shart bajarilsa , u holda 9(h) o’zining eng katta qiymatini Gh ning chegarasida ,
yani Γh daqabul qiladi.
Isboti. Teskarisini faraz qilamiz. Aytaylik , 9(h) o’zining eng katta
qiymatini ichki nuqtada qabul qilsin. Umuman aytganda , shunday nuqtalar ko’p
bolishi mumkin. Ular orasida shunday (i, k)eGh tugunni tanlaymizki,
9i+1,k ,9i-1,k ,9i,k+1 ,9i,k-1 qiymatlarning birortasida 9ik dan qat’iyan kichik , masalan,
9i+1,k < 9ik bo’lsin. U holda (i, k) tugunda quyidagiga egabo’lamiz:
Δ h9+ik = (9i-1,k - 29ik +9i+1,k )+ (9i,k+1 - 29ik +9i,k-1 )=
= [(9i-1,k - 9ik )+ (9i+1,k9ik )]+ (1.2.28)
+ [(9i,k+1 - 9ik )+ (9i,k-1 - 9ik )]< 0,
chunki 9i+1,k -9ik < 0 bo’lib, qolgan kichik qavslar ichidagi ifoda musbat emas, (1.2.28) tengsizlik esa lemma shartiga ziddir. Demak , bizning farazimiz noto’g’ri
ekan. Shu bilan lemma isbotlandi.

Download 0.78 Mb.

Do'stlaringiz bilan baham:
1   2   3   4   5   6   7   8   9   10




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