2021 Mundarija


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

1.2. Elliptik tenglamalarni tor 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)

Endi (1.2.2)- (1.2.5) lardan foydalanib, (1.1.7) dan quyidagiga egabo’lamiz:

L
( ?2 u ?2 u ?u ?u )
hyik =〈aik 2 + cik 2 + dik + lik + gik u +
l ?x1 ?x2 ?x1 ?x2 J(i ,k)
+ 12aik (||(?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



? 4 u

? 4 u

?x1 ,

?x2







?3u

?3u

?x1 ,

?x2




M3 = 〈 〉, M = 〈 〉
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 xatolikhosilbolib , 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 tor sohato’g’ritortburchakdan
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 korinishga egaboladi, 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

max max

i + k = n + 1 G


δ i+ku

δx δx1 2




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 boladi. 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:
1   2   3   4   5   6   7   8   9   10




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