logo

To’lqin tenglamasi uchun qo’yilgan umumiy masalani chekli ayirmali usuli bilan yechish va algoritmini tuzish

Загружено в:

20.11.2024

Скачано:

0

Размер:

524.4365234375 KB
Mavzu: To’lqin tenglamasi uchun qo’yilgan umumiy masalani chekli
ayirmali usuli bilan yechish va  algoritmini  tuzish .
Reja:
1.Kirish
2.Asosiy qismi:
 Tor   tebranish   tenglamasi   uchun   umumiy   boshlang`ich-chegaraviy
masalaning  qo` yilishi.
 Bir parametrli ayirmali approksimatsiya.
O(τ2+h2)  approksimatsiyali masala;
 To’lqin tenglamasini oshkor sxemada aproksematsiya qilish
 Natijalar
3.Xulosa
4.Foydalanilgan adabiyotlar.
    Ilova Kirish
Chekli   ayirmalar   usulining   asosiy   g ‘ oyasiga   ko ‘ ra   biror   muhitda
berilgan   funksiya   to ‘ r   vektor   bilan   ifodalanadi ,   differensial   operatorlar   esa
hech   bo ‘ lmaganda   fazoviy   o ‘ zgaruvchilar   va   vaqt   bo ‘ yicha   to ‘ rlarda
ularning   ayirmali   analoglarini   approksimatsiyalaydi .   Vaqt   bo‘yicha
masalalarda   EHMning   real   vaqtiga   to‘g‘ri   keluvchi   yechimni   olishga
majburmiz va vaqt bo‘yicha hosilalarni hisoblash oson kechmaydi, chunki
yangi vaqt momentiga to‘g‘ri keluvchi yechim noma’lum. Shuning uchun
vaqtdan   bog‘liq   hadlarni   o‘z   ichiga   olgan   masalalarda   vaqt   to‘tlarida
integrallash   operatsiyasini   bajarish   maxsus   ta’riflarni   talab   qiladi.   Biz   bu
yerda   vaqt   bo‘yicha   masalani   qaraganimizda   bitta   nuqtada   berilgan
chegaraviy   shartli   masalalarni   tushunamiz.   Shuning   uchun   real   vaqt
momentida   yechiladigan   masalalarning   o‘ziga   xos   jihatlari   va   muhim
ahamiyati mavjud.
Chekli   ayirmalar   usulining   asosiy   g‘oyasiga   ko‘ra   biror   muhitda
berilgan funksiya to‘r vektor bilan ifodalanadi, differensial operatorlar esa
hech   bo‘lmaganda   fazoviy   o‘zgaruvchilar   va   vaqt   bo‘yicha   to‘rlarda
ularning   ayirmali   analoglarini   approksimatsiyalaydi.   Vaqt   bo‘yicha
masalalarda   EHMning   real   vaqtiga   to‘g‘ri   keluvchi   yechimni   olishga
majburmiz va vaqt bo‘yicha hosilalarni hisoblash oson kechmaydi, chunki
yangi vaqt momentiga to‘g‘ri keluvchi yechim noma’lum. Shuning uchun
vaqtdan   bog‘liq   hadlarni   o‘z   ichiga   olgan   masalalarda   vaqt   to‘tlarida
integrallash   operatsiyasini   bajarish   maxsus   ta’riflarni   talab   qiladi.   Biz   bu
yerda   vaqt   bo‘yicha   masalani   qaraganimizda   bitta   nuqtada   berilgan
chegaraviy   shartli   masalalarni   tushunamiz.   Shuning   uchun   real   vaqt
momentida   yechiladigan   masalalarning   o‘ziga   xos   jihatlari   va   muhim
ahamiyati mavjud.  Tor tebranish tenglamasi uchun umumiy boshlang`ich-
chegaraviy masalaning  qo` yilishi.
Birjinsli tor tebranish tenglamasini qaraymiz∂2u	
∂t1
2=	a2∂2u	
∂	x1
2+	f(x1,t1),   	0<	x1<l,  t1>	0
.	
x=	x1/l,  t=	at	1/l
  o`lchovsiz   kattaliklarni   kiritib   bu   tenglamani
quyidagicha yozamiz	
∂2u	
∂t2=	
∂2u	
∂	x2+	f(x,t),       	0<	x<1	,      0	<t≤	T
. (1)
Boshlang`ich momentda	
u(x,0)=	u0(x),    	
∂u(x,0)	
∂t	
=	¯u0(x)
(2)
shartlar   berilgan,   bu   erda     u
0 ( x )   –   boshlang`ich   chetlashish   va    	
¯u0(x)     -
boshlang`ich tezlik.
Tor oxirlari quyidagi berilgan qonun bo`yicha harakatlansin 	
u(0,t)=	μ1(t),    	u(1,t)=	μ2(t)
. (3)	
¯D	=	(0≤	x≤	1	,    	0≤	t≤	T	)
  sohada   isssiqlik   o`tkazuvchanlik
tenglamasini approksimatsiyalashda qo`llangan to`rga o`xshash 	
¯ωhτ  to`g`ri
to`rtburchakli to`r kiritamiz	
y=	y	j,    	y	
^¿=	yj+1,    	y¿=	yj−	1,    	yt=	y^¿−y
τ	,     	y¯t=	y−	y¿	
2τ	,¿
¿	
Λy	=	y¯xx,   	y¯tt=	
yt−	y¯t	
τ	=	^y	
^¿−2y+y¿	
τ2	,    	y
t0=	yt+	y¯t	
2	=	^y
^¿−	y
¿	
2	τ	¿¿
.
(1) da hosilalarni quyidagi formulalarga almashtiramiz	
∂2u	
∂t2	~	u¯tt,   	∂2u	
∂	x2~	u¯xx=	Λu	,   	f	~	ϕ
.
Quyidagi vaznli sxemalar oilasini qaraymiz y¯tt=Λ(σ	^y	
^¿+(1−2σ)y+σy
¿
)+ϕ,    ϕ=f(x,tj),
¿y0=μ1(t),   yI=μ2(t),    	y(x,0)=u0(x),    	yt(x,0)=~u0(x),(4)
bu yerda  	
~u(x)  ni keyinroq aniqlaymiz.
Chegaraviy   shartlar   va   birinchi   u ( x,0 ) =u
0 ( x )   boshlang`ich   shart  	
¯ωhτ
to`rda   aniq   bajariladi.  	
~u(x)   ni   shunday   tanlaymizki.	
~u0(x)−	
∂u(x,0)	
∂t	
=	~u0(x)−	¯u0(x)
  approksimatsiya   xatoligi  	
O	(τ2)   kattalik
bo`lsin.  Quyidagi 	
ut(x	,0)=	˙u	(x	,0	)+	0,5	τ	¨u(x	,0	)+	O	(τ2)=	
=	¯u0(x	)+	0,5	τ(u''(x	,0)+	f(x	,0))+	O	(τ2)=	
=	¯u0(x)+	0,5	τ(u0
''(x	,0)+	f(x	,0))+	O	(τ2)
formuladan   ko`rinadiki  	
~u(x)−	ut(x,0)=	O	(τ2)   belgilashni   qo`ysak
quyidagini hosil qilamiz 	
~u(x)=	¯u0(x)+0,5	τ(¯u0
''(x)+	f(x,0)
. (5)
 Bir parametrli ayirmali approksimatsiya.
Shunday qilib, (4), (5) masala qo`yildi. (4) dan  
y	
^¿=yj+1
¿   ni aniqlash uchun
progonka usuli bilan yechiladigan chegaraviy masalani hosil qilamiz	
σγ	2(yi+1
j+1+	yi−1
j+1)−	(1+2	σγ	2)yi
j+1=	−	F	i,   	0<	i<	I,   	y0=	μ1,  	yI=	μ2
,	
γ=	τ/h	,   	F	i=	(2	yi
j−	yi
j−1)+τ2(1−	2σ	)Λy	j+στ	2Λy	j−1+	τ2ϕ
.
Bunda 	
σ>	0 bo`lganda progonka usuli turg`un bo`ladi.
	
O(τ2+h2)   approksimatsiyali masala	
ϕ=	f(x,tj)
  da (4) sxema approksimatsiya xatoligini hisoblaymiz.   y   – (4),
(5) masalaning, 	
u=	u(x,t)  - esa (1)-(3) masalaning yechimlari bo`lsin. (4)
ga 	
y=	z+	u qo`yib quyidagini hosil qilamiz z¯tt=	Λ	(  	σ	z	
^¿+(1−2σ)z+σz¿)+ψ	
¿,  	
z0=	zI=	0,   	z(x,0)=	0,   	zt(x,0)=	ν(x)
, (6)
bu   erda  	
ψ	=	Λ	(  	σ	^u	
^¿+(1−	2σ)u+σu¿)+ϕ−	u¯tt
¿     -   (4)   sxemaning   u=u ( x , t )
yechimdagi   approksimatsiya   xatoligi ,  	
ν=	~u0(x)−	ut(x,0) -   esa  	yt=	~u0(x)
ikkinchi   boshlang`ich   shart   uchun   approksimatsiya   xatoligi.
Yuqoridagilardan ayonki, 	
ν=	O	(τ2) .	
u	
^¿=	u+τut
¿
, 	u
¿
=	u−	τ	u¯t   lardan foydalanib quyidagiga ega  bo`lamiz  	
σ	u	
^¿+(1+2σ)u+σu¿=	u+στ	2u¯tt
¿
,   (7)
ya`ni har qanday  	
   (	 , 	   va  h  dan bog`liq emas) 	
ψ	=	Λu	+στ	2Λu	¯tt+ϕ−	u¯tt	=	Lu	+	στ	2L	¨u+	f−	¨u+	O	(τ2+	h2)
,	
ψ	=	O	(τ2+h2)
.
3-tur chegaraviy shartlar 	
∂u(0,t)	
∂x	
=	β1u(0,t)−	μ1(t),     	−	
∂u(1,t)	
∂x	
=	β2u(1,t)−	μ2(t)
quyidagi yoyilmalarni qo`llab approksimatsiya qilinadi	
∂u
∂	x
=	ux−	h
2	
u''−	h2
6	
u'''−	h3	
24	
u(IV	)+O	(h4)
va 	
∂u
∂	x
=	u	x+	h
2	
u''−	h2
6	
u'''+	h3	
24	
u(IV	)+O	(h4)
.	
O	(τ2+	h2)
 approksimatsiya tartibini hosil qilish uchun
i =0 da   	
ux−	h
2	
u''=	β1u−	μ1(t)+O	(h2)  
yoki	
ux−	β1u	
0,5	h	
+	
μ1(t)	
0,5	h	
=	u''+O	(h2)
ifodalarni qo`llaymiz. 
Tenglamaning o`zi dan  quyidagiga ega bo`lamiz	
uxx
''=	utt
''−	f
. y   ayirmali funktsiya uchun quyidagiga egamizu¯tt=	Λ−(σ	y+(1−	2σ	)y+	σ	y
¿
)y+	ϕ−	,    	i=	0
,
bu erda	
Λ−	y=	
yx−	β1y	
0,5	h	
,    	ϕ−=	ϕ+	
μ1(t)	
0,5	h	
,   	ϕ=	f.
Shunga o`xshash	
u	¯t	t	=	Λ	
+	
¿	¿
,
bu yerda 	
Λ+=	−	
y¯x−	β2y	
0,5	h	
,    	ϕ+=	ϕ+	
μ2(t)	
0,5	h	
,   	ϕ=	f.
Bundan   tashqari  	
∂u
∂x   yoyilmada   yanada   yuqori   tartibli   hosilalarni
qo`llab 	
O	(τ2+	h4)  aniqlik bilan sxemalar hosil qilish mumkin. 
 To’lqin   tenglamasini   oshkor   sxemada   approksimatsiya
qilish
Endi   masalani   umumiyroq   qo`yilishini   qaraymiz.	
¯D	T=	[0≤	x≤	1	]×	[0≤	t≤	T	]
  to`g`ri   to`rtburchakda   giperbolik   tipli
tenglama uchun 1-chegaraviy masalani qaraymiz	
∂2u	
∂t2=	Lu	+	f(x,t),   	Lu	=	∂
∂	x(k(x,t)∂u
∂	x),	
(x,t)∈	DT,
  (8)	
u(x,0)=	u0(x),    	
∂u(x,	0)	
∂t	
=	¯u0(x),
(9)	
u(0,t)=	u1(t),    	u(1,t)=	u2(t),
(10)	
0<¯¯c1≤	k(x,t)≤	¯¯c2,
bu erda 	
D	T=	(0≤	x≤	1)×	(0<	t≤	T	].
Faraz qilaylik masala 	
¯DT  da  yagona  uzluksiz va  etarlicha  hosilalar ga
ega   bo`lgan   echimga   ega   bo`lsin.   k ( x , t )   (va   f ( x , t )   o`ng   taraf     ot     o`qqa parallel   chekli   sondagi   to`g`ri   chiziqlarda   1-   tur   uzilishga   ega   bo`lishi
mumkin   (qo`zg`almaydigan   uzilishlar).   Har   bir   chiziqdax=	ξs,   	s=	1,2	,...,s0
 uzilishlarda qo`shmalik sharti bajariladi:	
[u]=	u(ξs+0,t)−	u(ξs−	0,t)=	0¿}¿¿¿
. (11)
Endi   (8)-(10)   masala   uchun   bir   jinsli   ayirmali   sxemalarni   tuzishga
kelamiz.  	
¯^ω	h=	{xi,i=	0,1	,...,I,   	x0=	0,    	xI=	1}   to`r    	0≤	x≤	1   da   ixtiyoriy
notekis to`r,  	
¯^ω	τ=	{tj=	j⋅τ,  	j=	0,1,2	,...,J}  to`r  	0≤	t≤	T  da ixtiyoriy tekis
to`r,  	
¯ω	hτ=	¯^ω	h×	¯ω	τ   -   to`r   esa  	DT   to`g`ri   to`rtburchakda   berilgan   to`r
bo`lsin.   Avvalo   fiksirlangan  	
t∈ω	τ   da     Lu+f   operatorni
approksimatsiyalaymiz   va  	
Λu	+	ϕ=	(a(x	,t)u¯x)^x+	ϕ   -   ayirmali   operatorga
keltiramiz.
Bularni quyidagicha almashtiramiz 	
∂2u	
∂t2|t=tj~	u¯tt,    	Lu	+	F	=	Λ(tj)u
(σ1,σ2)
+ϕ
,
Bu erda  	
u(σ1,σ2)=	σ1	u	
^¿+(1−	σ1−	σ2)u+	σ2u
¿
,
¿	
Λ	(tj)u=	(a(x,tj)u¯x)
^x^¿,    	u=uj,   	u
¿
=	u	j−1,   	^u
^¿
=	uj+1¿
¿
.
Quyidagi belgilashlarni eslaymiz	
v¯x,i=	
vi−	vi−1	
hi	
,
   	vx,i=	
vi+1−	vi	
hi+1	
=	v¯x,i+1,   	v^x,i=	
vi+1−	vi	
ℏi	
,
bu erda  	
ℏi=	0,5	(hi+hi+1),
 	Lv	=	d2v	
dx	2	,  	
(Lhv)i=	1
ℏi[
vi+1−	vi	
hi+1	
−	
vi−	vi−1	
hi	]−
=	v¯x^x,i=	v¯x^x
, 
va bulardan quyidagi vaznli bir jinsli uch qatlamli sxemani hosil qilamiz	
y¯tt=	Λ	(tj)y(σ1,σ2)+	ϕ
. (12) t=t
j   o`rta qatlamda  a   koeffitsientni olamiz 
Quyidagi  ^y=	y+τy	t+0,5	τ2y¯tt ,  	
˘y=	y−	ε	y
t
0+0,5	τ2y¯tt ,   bu   erda	
yt
∘=	(˘y+	y)/(2τ)
,  	y¯tt=	(y−	2	y+	˘y)/τ2   larni   qo`llab  	
y(σ1,σ2)=	
=	y+(σ1−	σ2)εy	
t
0+	0,5	(σ1+σ2)τ2y¯tt
 ni hosil qilamiz, bulardan keyin  (12)
sxemani quyidagicha yozamiz	
(E	−	0,5	(σ	1+	σ	2)τ2Λ	)y¯tt−	(σ1−	σ	2)τΛy	
t
0=	Λy	+	ϕ
, (13)
bu   erda   E   –   birlik   operator.  	
σ1=	σ2=	σ   da   simmetrik     sxemani   hosil
qilamiz	
(E	−	στ	2Λ	)y¯tt=	Λy	+	ϕ(x,t),   	0<	t=	τj
  (14)
va uni o`rganish bilan cheklanamiz.
(10)   chegaraviy   shartlar   va   (9)     birinchi   boshlang`ich   shart   aniq
qanoatlantir il adi 	
y(0,t)=	u1(t),   	y(1,t)=	u2(t),   	y(x,0)=	u0(x)
. (15)	
∂u/∂t|t=0=	¯u0(x)
  ikkinchi   boshlang`ich   shartni   ikki   usul   bilan
approksimatsiyalash mumkin. Bit t a usul yuqorida ko`rsatildi	
yz(x,0)=	~u0(x),    {	~u0(x)=	¯u0(x)+0,5	τ(Lu	0+	f)t=0¿
. (16)
U 	
   bo`yicha ikkinchi tartibli approksimatsiyaga ega.
Ikkinchi   usul   shundan   iboratki,   y (	
 )   ni   aniqlash   uchun   quyidagi
ayirmali tenglama yoziladi	
(E	−	στ	2Λ(0))yt(x,0)=	¯u0(x)+0,5	τ(Λu	0+	f(x,0)).
(17)
Natijada (14)-(16) (yoki (14), (15), (17)) ayirmali masalani hosil qilamiz.
Bu   sxema   uch   qatlamli   deyiladi.   YAngi   qatlamdagi  	
y	
^¿=yj+1
¿   ni
hisoblash   uchun   avvalgi   ikkita   qatlamdagi  	
yj   va  	yj−1   qiymatlarni   bilish
kerak.   Har   bir  	
t=	tj+1   yangi   qatlamda   chegaraviy   masala  	y	
^¿=yj+1
¿     ga
nisbatan echiladi (progonka usuli bilan): 	
(E	−	στ	2Λ	)	y
^¿=F,
¿
 
 	
0<x=	ih	<1,     {	^y0=	^u1,   {	^y¿I=	^u2¿ , F	(t)=	2	y−	y
¿
−	τ
2
Λ	((2σ−	1)y−	σ	y
¿
)+	τ
2
ϕ	,    	t≥	τ,	
F	(0)=	u0+τ2(0,5	−	σ	)Λ	(0)u0+	τ¯u0(x)+0,5	τ2f(x,0). Natijalar
Testlash jarayoni:
Koordinatalar oralig'ini kiriting: 3 5 
Vaqt oralig'ini kiriting: 0 10
Nuqtalar oralig'i masofasini kiriting: 0.1 
Vaqt oralig'i masofasini kiriting: 0.2 
O'tkazuvchanlikni kiriting: 4 
Koordinatalar oralig'ini kiriting: 3 5
Vaqt oralig'ini kiriting: 5 6
Nuqtalar oralig'i masofasini kiriting: 0.2
Vaqt oralig'i 
masofasini kiriting: 
0.2 
O'tkazuvchanlikni 
kiriting: 4 Koordinatalar oralig'ini kiriting: 1 2 
Vaqt oralig’ini kiriting: 2 4 
Nuqtalar oralig'i masofasini kiriting:0.1 
Vaqt oralig'i masofasini kiriting: 0.02 
O'tkazuvchanlikni kiriting: 2 Xulosa
          Xulosa   qilib   shuni   ta’kidlash   mumkinki,   biz   yashab   kelayotgan
jarayonda,   turli   tuman   muhandislik   va   amaliy   masalalar   yuzaga   keladi,
ularni   hal   etishning   yagona   chorasi   raqamlarga   yuzlanish   bo’ladi.   Biz
qo’llagan   usullar   bilan   siz   hal   etmoqchi   bo’lgan   turli   muhandislik
masalalaringizni osongina hal etishingiz mumkin.
Yuqoridag   python   dasturlash   tilidan   foydalangan   holda   tuzulgan   dastur
natijasini ko’rib buni yaqqol aytishimiz mumkin.
Adabiyotlar
1. Самарский   А.А.   Теория   разно стных   схем.   М:   Наука,   1989.
(407-410  betlar )
2. Х ў жаёров   Б.Х.   Қ урилиш   масалаларини   сонли   ечиш   усуллари.
Тошкент,  “Ў збекистон ” , 1995.  ( 180-189 betlar) 3. Турчак Л.И. Основы численных методов. М.: Наука, 1987.  (272-277
betlar)
4. Калиткин   Н.Н.   Численные   методы   М:   Наука,   1978.   (424-439
betlar )
5. Самарский   А.А.,   Гулин   А.В.   Численные   методы.   М:   Наука,
1989.  (283-286 betlar)
Ilova:
import  matplotlib.pyplot  as  plt
import  numpy  as  np
a,b =  map ( float , input ( "Koordinatalar oralig'ini kiriting: "  ).split())
c,d =  map ( float , input ( "Vaqt oralig'ini kiriting: " ).split())
h =  float ( input ( "Nuqtalar oralig'i masofasini kiriting: " ))
tau =  float ( input ( "Vaqt oralig'i masofasini kiriting: " ))
m =  float ( input ( "O'tkazuvchanlikni kiriting: " ))
# n nuqtalar soni
n =  int ((b-a)/h)
# k vaqtlar soni
k =  int ((d-c)/tau)
l = np.empty((k+ 1 ,n+ 1 ))
# i ustun
i =  0
# j satr
j =  0
for  i  in   range (n+ 1 ):
l[ 0 ][i] =  0 for  j  in   range (k+ 1 ):
l[j][ 0 ] =  20
l[j][n] =  10
for  j  in   range ( 1 ,k):
for  i  in   range ( 1 ,n):
l[j+ 1 ][i] = (tau** 2 *m** 2 *(l[j][i+ 1 ] -2 *l[j][i]+l[j][i -1 ]))/(h*h)+ 2 *l[j][i]+l[j -1 ][i]
for  j  in   range (k+ 1 ):
for  i  in   range (n):
print (l[j][i],end= ' ' )
print ()
for  j  in   range (k+ 1 ):
if  j %  10  ==  0 :
plt.plot(l[j][ 0 :n+ 1 ])

Mavzu: To’lqin tenglamasi uchun qo’yilgan umumiy masalani chekli ayirmali usuli bilan yechish va algoritmini tuzish . Reja: 1.Kirish 2.Asosiy qismi:  Tor tebranish tenglamasi uchun umumiy boshlang`ich-chegaraviy masalaning qo` yilishi.  Bir parametrli ayirmali approksimatsiya. O(τ2+h2) approksimatsiyali masala;  To’lqin tenglamasini oshkor sxemada aproksematsiya qilish  Natijalar 3.Xulosa 4.Foydalanilgan adabiyotlar. Ilova

Kirish Chekli ayirmalar usulining asosiy g ‘ oyasiga ko ‘ ra biror muhitda berilgan funksiya to ‘ r vektor bilan ifodalanadi , differensial operatorlar esa hech bo ‘ lmaganda fazoviy o ‘ zgaruvchilar va vaqt bo ‘ yicha to ‘ rlarda ularning ayirmali analoglarini approksimatsiyalaydi . Vaqt bo‘yicha masalalarda EHMning real vaqtiga to‘g‘ri keluvchi yechimni olishga majburmiz va vaqt bo‘yicha hosilalarni hisoblash oson kechmaydi, chunki yangi vaqt momentiga to‘g‘ri keluvchi yechim noma’lum. Shuning uchun vaqtdan bog‘liq hadlarni o‘z ichiga olgan masalalarda vaqt to‘tlarida integrallash operatsiyasini bajarish maxsus ta’riflarni talab qiladi. Biz bu yerda vaqt bo‘yicha masalani qaraganimizda bitta nuqtada berilgan chegaraviy shartli masalalarni tushunamiz. Shuning uchun real vaqt momentida yechiladigan masalalarning o‘ziga xos jihatlari va muhim ahamiyati mavjud. Chekli ayirmalar usulining asosiy g‘oyasiga ko‘ra biror muhitda berilgan funksiya to‘r vektor bilan ifodalanadi, differensial operatorlar esa hech bo‘lmaganda fazoviy o‘zgaruvchilar va vaqt bo‘yicha to‘rlarda ularning ayirmali analoglarini approksimatsiyalaydi. Vaqt bo‘yicha masalalarda EHMning real vaqtiga to‘g‘ri keluvchi yechimni olishga majburmiz va vaqt bo‘yicha hosilalarni hisoblash oson kechmaydi, chunki yangi vaqt momentiga to‘g‘ri keluvchi yechim noma’lum. Shuning uchun vaqtdan bog‘liq hadlarni o‘z ichiga olgan masalalarda vaqt to‘tlarida integrallash operatsiyasini bajarish maxsus ta’riflarni talab qiladi. Biz bu yerda vaqt bo‘yicha masalani qaraganimizda bitta nuqtada berilgan chegaraviy shartli masalalarni tushunamiz. Shuning uchun real vaqt momentida yechiladigan masalalarning o‘ziga xos jihatlari va muhim ahamiyati mavjud.

 Tor tebranish tenglamasi uchun umumiy boshlang`ich- chegaraviy masalaning qo` yilishi. Birjinsli tor tebranish tenglamasini qaraymiz∂2u ∂t1 2= a2∂2u ∂ x1 2+ f(x1,t1), 0< x1<l, t1> 0 . x= x1/l, t= at 1/l o`lchovsiz kattaliklarni kiritib bu tenglamani quyidagicha yozamiz ∂2u ∂t2= ∂2u ∂ x2+ f(x,t), 0< x<1 , 0 <t≤ T . (1) Boshlang`ich momentda u(x,0)= u0(x), ∂u(x,0) ∂t = ¯u0(x) (2) shartlar berilgan, bu erda u 0 ( x ) – boshlang`ich chetlashish va ¯u0(x) - boshlang`ich tezlik. Tor oxirlari quyidagi berilgan qonun bo`yicha harakatlansin u(0,t)= μ1(t), u(1,t)= μ2(t) . (3) ¯D = (0≤ x≤ 1 , 0≤ t≤ T ) sohada isssiqlik o`tkazuvchanlik tenglamasini approksimatsiyalashda qo`llangan to`rga o`xshash ¯ωhτ to`g`ri to`rtburchakli to`r kiritamiz y= y j, y ^¿= yj+1, y¿= yj− 1, yt= y^¿−y τ , y¯t= y− y¿ 2τ ,¿ ¿ Λy = y¯xx, y¯tt= yt− y¯t τ = ^y ^¿−2y+y¿ τ2 , y t0= yt+ y¯t 2 = ^y ^¿− y ¿ 2 τ ¿¿ . (1) da hosilalarni quyidagi formulalarga almashtiramiz ∂2u ∂t2 ~ u¯tt, ∂2u ∂ x2~ u¯xx= Λu , f ~ ϕ . Quyidagi vaznli sxemalar oilasini qaraymiz

y¯tt=Λ(σ ^y ^¿+(1−2σ)y+σy ¿ )+ϕ, ϕ=f(x,tj), ¿y0=μ1(t), yI=μ2(t), y(x,0)=u0(x), yt(x,0)=~u0(x),(4) bu yerda ~u(x) ni keyinroq aniqlaymiz. Chegaraviy shartlar va birinchi u ( x,0 ) =u 0 ( x ) boshlang`ich shart ¯ωhτ to`rda aniq bajariladi. ~u(x) ni shunday tanlaymizki. ~u0(x)− ∂u(x,0) ∂t = ~u0(x)− ¯u0(x) approksimatsiya xatoligi O (τ2) kattalik bo`lsin. Quyidagi ut(x ,0)= ˙u (x ,0 )+ 0,5 τ ¨u(x ,0 )+ O (τ2)= = ¯u0(x )+ 0,5 τ(u''(x ,0)+ f(x ,0))+ O (τ2)= = ¯u0(x)+ 0,5 τ(u0 ''(x ,0)+ f(x ,0))+ O (τ2) formuladan ko`rinadiki ~u(x)− ut(x,0)= O (τ2) belgilashni qo`ysak quyidagini hosil qilamiz ~u(x)= ¯u0(x)+0,5 τ(¯u0 ''(x)+ f(x,0) . (5)  Bir parametrli ayirmali approksimatsiya. Shunday qilib, (4), (5) masala qo`yildi. (4) dan y ^¿=yj+1 ¿ ni aniqlash uchun progonka usuli bilan yechiladigan chegaraviy masalani hosil qilamiz σγ 2(yi+1 j+1+ yi−1 j+1)− (1+2 σγ 2)yi j+1= − F i, 0< i< I, y0= μ1, yI= μ2 , γ= τ/h , F i= (2 yi j− yi j−1)+τ2(1− 2σ )Λy j+στ 2Λy j−1+ τ2ϕ . Bunda σ> 0 bo`lganda progonka usuli turg`un bo`ladi.  O(τ2+h2) approksimatsiyali masala ϕ= f(x,tj) da (4) sxema approksimatsiya xatoligini hisoblaymiz. y – (4), (5) masalaning, u= u(x,t) - esa (1)-(3) masalaning yechimlari bo`lsin. (4) ga y= z+ u qo`yib quyidagini hosil qilamiz

z¯tt= Λ ( σ z ^¿+(1−2σ)z+σz¿)+ψ ¿, z0= zI= 0, z(x,0)= 0, zt(x,0)= ν(x) , (6) bu erda ψ = Λ ( σ ^u ^¿+(1− 2σ)u+σu¿)+ϕ− u¯tt ¿ - (4) sxemaning u=u ( x , t ) yechimdagi approksimatsiya xatoligi , ν= ~u0(x)− ut(x,0) - esa yt= ~u0(x) ikkinchi boshlang`ich shart uchun approksimatsiya xatoligi. Yuqoridagilardan ayonki, ν= O (τ2) . u ^¿= u+τut ¿ , u ¿ = u− τ u¯t lardan foydalanib quyidagiga ega bo`lamiz σ u ^¿+(1+2σ)u+σu¿= u+στ 2u¯tt ¿ , (7) ya`ni har qanday  (  ,  va h dan bog`liq emas) ψ = Λu +στ 2Λu ¯tt+ϕ− u¯tt = Lu + στ 2L ¨u+ f− ¨u+ O (τ2+ h2) , ψ = O (τ2+h2) . 3-tur chegaraviy shartlar ∂u(0,t) ∂x = β1u(0,t)− μ1(t), − ∂u(1,t) ∂x = β2u(1,t)− μ2(t) quyidagi yoyilmalarni qo`llab approksimatsiya qilinadi ∂u ∂ x = ux− h 2 u''− h2 6 u'''− h3 24 u(IV )+O (h4) va ∂u ∂ x = u x+ h 2 u''− h2 6 u'''+ h3 24 u(IV )+O (h4) . O (τ2+ h2) approksimatsiya tartibini hosil qilish uchun i =0 da ux− h 2 u''= β1u− μ1(t)+O (h2) yoki ux− β1u 0,5 h + μ1(t) 0,5 h = u''+O (h2) ifodalarni qo`llaymiz. Tenglamaning o`zi dan quyidagiga ega bo`lamiz uxx ''= utt ''− f .