2021 Mundarija
Download 0.78 Mb.
|
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 to’rtinchi 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 to’rtburchaklardan 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 bo’lishi 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: |
Ma'lumotlar bazasi mualliflik huquqi bilan himoyalangan ©fayllar.org 2024
ma'muriyatiga murojaat qiling
ma'muriyatiga murojaat qiling