2021 Mundarija
Download 0.78 Mb.
|
327798 (1)
- Bu sahifa navigatsiya:
- Ayirmali tenglama hosil qilish uchun aniqmas koeff isientlar metodi.
- Puasson
1.2. Elliptik tenglamalarni to’r metodi bilan yechish
Quyidagi L(u) = a + c + d + e + gu = f (1.2.1) Elliptik tenglamani (1.1.7) ayirmali tenglama bilan almashtirganda hosil bo’ladigan xatoliknibaholashniko’rib chiqamiz. Bu yerda hisoblashlar soda bo’lishi uchun aralash hosilaningoldidagikoeffisientni b(x1 , x2 ) = 0 deb oldik. (1.2.1) differensialtenglamaning u(x,y) yechimini to’rtinchitartiblixususiy hosilalarga ega deb faraz qilib va Teylor formulasidan foydalanib, (1.1.3)- (1.1.6) taqribiy tengliklar o’rnida quyidagilarni hosil qilamiz: = (x1i ,x2 k ) + (ξ,x2 k ) y 2h2 |( δx2 )|(x1j ,x2 k ) 6 |( δx2 )|(ξ1 ,x2 k ) i ,k+1 - yi ,k-1 ,i u ξ < x1, 1 1 )| ( 2,k-1 2, +1 < x x < η k ) = (x1j ,x2 k ) + (ξ1 ,x2 k ) = (x1j ,x2 k ) + (x1j ,η1 ) (1.2.2) (1.2.3) (1.2.4) (1.2.5) L ( ?2 u ?2 u ?u ?u ) hyik =〈aik 2 + cik 2 + dik l ?x1 ?x2 ?x1 ?x2 J(i ,k) + 12〈 aik (||(?x (ξ1 ,x2 k ) +a2 cik (||(?x (x1i ,η1 ) + 2dik (||(?x (ξ,x2 k ) + 2alik (||(?x (x1i ,η) = = [L(u)](i ,k) + Ri, k, Bunda h = h1 ,a = h1 bo’lib, Rik -qoldiq had. Agarushbu m
M3 = ax (| )| max (| )| G |l J| 4 G |l J| Belgilashlarni kritsak, qoldiq had uchun Rik < 12 {( aik +a2 bik )M4 + 2( dik +al ik )M3 } (1.2.6) Baho o’rinli bo’ladi. Demak , Lhyik - fik = {L(u)- f}(i ,k) + Rik = Rik Bundanko’ramizki, (1.2.1) differensialtenglamani (1.1.7) ayirmali tenglama bilan almashtirganda Rik xatolikhosilbo’lib , uninghqadamganisbatan tartibi h 2 dir. Agar Rik qoldiq hadnitashlasak, to’rustidagi yik funksiyauchun Lhyik = fik (1.2.7) Tenglamalar sistemasiga egabo’lamiz. Xususiy holda ushbu + = f (x1 , x2 ) (1.2.8) Puasson tenglamasiuchun h1 = h2 = h kvadrat to’rniqarasak , u holda (1.2.7) tenglamalar sistemasi yi+1,k + yi-1,k + yi ,k+1 + yi ,k-1 - 4yik = h2 fik Ko’rinishga egabo’lib, (1.2.6) dan qoldiq had uchun Rik < M4 bahoga egabo’lamiz. (1.2.9) (1.2.10) Ayirmali tenglama hosil qilish uchun aniqmas koeffisientlar metodi. Yuqoridagi differensialtenglamani (i, k) nuqtada ayirmali tenglama bilan almashtirganda har bir xususiy hosilali alohida-alohida bo’lingan ayirmalar bilan almashtirgan edik. Differensialtenglamani to’laligicha ayirmali tenglama bilan almashtirish mumkin. Hozirqaraladigan metodda to’r sohato’g’rito’rtburchakdan iborat bo’lishi shart emas, to’ruchburchaklar , parallelogramlardan iborat yoki umuman notekis bo’lishi ham mumkin. Differensialtenglamani (i, k) tugunda ayirmalisxema bilan almashtirish uchun (i, k) tugunatrofinima’lumtartibda joylashgan r tatugunniqaraymiz. Qulay bo’lishi uchun (i, k) tugunni 0 orqali belgilab , qolgantugunlarni 1,2,......., P orqalibelgilaymiz. Endi с i aniqmas koeffisientlarbilan ushbu cj uj (1.2.11) chiziqli kombinatsiyanituzamiz, bunda uj miqdor u ning j tugundagi qiymati. Faraz qilaylik , u funksiya (n +1) tartibli hosilalarga egabo’lsin, u holda uj larni 0 tugun atrofida Teylor qatorigayoyamiz: u k1 +k2 =n 1 j = u(x1j , x2j )= Σ - - 0 + R(j), (1.2.12) j = 1,2,.....,... . Bu ifodalarni (1.2.11) ga qo’yib , u funksoyaning bir xil hosilalari oldidagi koeffisientlarni qo’shib chiqamiz, natijada cj uj = 0< naik 0 + cjR(j) (1.2.13) Bu yerda aik koeffisientlar cj lar orqalichiziqli ravishda ifodalanadi. Qoldiq hadesa θhn+1KMn+1 ko’rinishga egabo’ladi, bunda θ < 1 , K qandaydir son bo’lib, h bog’liq emas; h ningo’zi esa 0 tugun va j(j = 1,2,....., P) tugunlar koordinatalari ayirmalarining moduli bo’yicha eng kichigi hamda M = n-1
Endi G sohada (n +1) tartibli uzluksizhosilaga egabo’lgan har qanday u(x1 , x2 ) funksiyauchun i< nik 0 = [L(u)]0 (1.2.14) tenglikningbajarilishinitalab qilamiz. Buninguchun cj koeffisientlarni shunday tanlashimiz kerakki , 0 < i + k < nshartni qanoatlantiruvchibarcha i va k uchun (1.2.14) tenglikning chap va o’ng tomonlarida 0 oldidagi koeffisientlar ustma-ust tushsin. Bu esa c1 , c2 ,..., cp noma’lum koeffisientlarga nisbatan quyidagi chiziqli algebraik tenglamalar sistemasiga olib keladi: a00 = g0 (i + k = 0), a10 = d0 ,a01 = l0 (i + k = 1), a20 = a0 ,a02 = c0 ,a11 = 0(i + k = 2), a30 = a21 = a12 = a03 = 0(i + k = 3), an0 = an-1,1 = ...a1,n-1 = a0n = 0(i + k = n). Agar bu sistema yechimgaega bo’lib, yechim cj (j=0,1,….,P) bo’lsa, u holda cj uj = [L(u)]0 + θKhn+1Mn+1 (1.2.15) Endiqoldiq hadnitashlabyuborib, uj ning to’rustidagitaqribiy qiymati yj uchun ushbu cj yj = f 0 (1.2.16) Ayirmalitenglamaga egabo’lamiz. Bu tenglama (1.2.1) differensialtenglamani 0 tugunda O(hn+1) aniqlikdaalmashtiradi. Chegaradan uzoqroq ichkitugunlaruchun ayirmali tenglamanituzishda qatnashadigantugunlarining joylanishini (1.2.16) dagideksaqlashmaqsadga muvofiq bo’ladi. Chegaragayaqin tugunlaruchun buholatni saqlash har doim ham mumkin bo’lavermaydi. Ammo qaralayotgametoddatugunlarni biroz boshqacha joylashtirib , differensialtenglamanikerakli aniqlikda ayirmali tenglama bilan almashtirish mumkin. Bu metod chegaraviy shartlarni approksimatsiya qilish uchun yaxshinatujaga olib keladi. Puasson tenglamasi uchun aniqmas koeffisientlar metodi asosida 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