Дифференциалдық теңдеулерді сандық шешу



бет1/3
Дата06.01.2022
өлшемі95,5 Kb.
#15887
  1   2   3

Дифференциалдық теңдеулерді сандық шешу

1. Ақырлы – айырымдық әдіс.

2. Қуалау әдісі.

F(x,y,y,y’’ )=0 (1)

2-ретті дифференциалдық теңдеу берілсін.

Біз осыған дейінгі есептерде Коши есебі деп дифференциалдық теңдеуге қосымша функцияның бір тәуелсіз айнымалыдағы мәні берілген жағдайды айтып жүрдік. Ал функцияның екі тәуелсіз айнымалыдағы мәндері берілсе, есепті ҚДТ үшін шектік есеп деп айтады. Қосымша шарттарды шекаралық шарттар дейді. Коши есебінде тәуелсіз айнымалы ретінде уақытты қарастырғанбыз, сонда есептің физикалық мағынасы серіппеге ілінген дененің еркін қозғалысын зерттеу болған. Ал шектік есепте тәуелсіз айнымалы ретінде ұзындық қарастырылады, яғни есептің физикалық мағынасы қатты серіппенің деформациясын зерттеу немесе, электр желілерін басқару, есептеу, химиялық реакциялардың нәтижелерін есептеу, снаряд қозғалысының теңдеуін құру сияқты практикада қарапайым классикалық әдістермен шешілмейтін «қатты» есептерді зерттеу болып табылады.

(1) теңдеу үшін екінүктелі шектік есеп келесідей қойлады: [a,b] аралығында (1) теңдеуін және осы аралықтың екі шеткі нүктелерінде

1[y(a), y /(a)]=0 (2)

2[y(b),y /(b)]=0

шекаралық шарттарын қанағаттандыратын y=y(x) шешімін табу керек.

(1) теңдеу және (2)- шекаралық шарт сызықты болған жағдайды қарастырайық. Мұндай шектік есептер сызықты шектік есептер деп аталады және төмендегідей түрде жазылады:

y // +p(x) y / + q(x) y = f(x) (3)



(4)

мұндағы p(x), q(x), f(x) – функциялары [a,b] аралығында үзіліссіз функциялар, ал 0,1,0,1, A,B –берілген тұрақтылар және 0+10, 010.

Егер A=B=0 болса, онда (4) – шекаралық шарттар біртекті шекаралық шарттар деп аталады.

Қарапайым дифференциалдық теңдеулер үшін қойылған шектік есепті шешудің жалпы 2 әдісі бар:



  1. Дифференциалдық теңдеудің ақырлы – айырымдық немесе шектік – айырымдық түрін қолдану – сандық әдіс.

  2. Шектік есепті бірнеше Коши есебін шешуге келтіру аналитикалық әдіс.

1. Ақырлы – айырымдық әдіс.

Бұл әдістің негізгі идеясы - шектік есепті шешуді алгебралық теңдеулер жүйесін шешуге келтіру.

Берілген аралықты бірнеше бөлікке бөлу арқылы бірдей қашықтықта орналасқан түйіндер жүйесін құрамыз: x0=a, xn=b, xi=x0+ih (i=1,2,…,n-1), қадамы h=b-a/n және pi=p(xi), qi=q(xi), fi=f(xi) болсын.

xi түйіндерінде ізделінді функцияның жуық мәнін y(x) және y /(x) ,y //(x) туындыларын сәйкесінше yi ,yi/ ,yi // деп белгілейік. Құрылған түйіндер жүйесінің әрбір ішкі түйінінде yi/(xi), yi//(xi) туындыларын сәйкес ақырлы – айырымдық қатынастарымен алмастырамыз:



(5)

Шекаралық шарттар үшін:



(5!)

қатынастарын қолданылады.

(5) – қатынастарды (1) – ге қойсақ және шекаралық шарттарды ескерсек, сызықты алгебралық теңдеулер жүйесін немесе сызықты емес теңдеулер жүйесін аламыз. Бұл жерде (1)-теңдеудің сипатына көңіл аудару арқылы оның қандай жүйеге келетінін алдын ала анықтауға болады. Егер (1)-теңдеу сызықты болса, онда есеп сызықты алгебралық теңдеулер жүйесін шешуге келеді, ал сызықты емес болса, онда сызықты емес теңдеулер жуйесін шешуге кейде трансцендентті теңдеуді шешуге келеді.

Егер құрылған жүйені шешу қиынға соқса, онда шектік есепті екі Коши есебін шешуге келтіруге болады.



(3), (4)- шектік есепті қарастырайық. (4)-шекаралық шартта деп есептесек,

(6)

екенін көреміз. Мұндағы , - төмендегі түрге келтірілген Коши есептерінің шешімдері болады:

(7)
(8)
(7)-есептің шешімі , (8)-есептің шешімі , ал (6)-теңдеудегі с – тұрақты, белгісіз. Оны анықтау үшін (4)-шекаралық шартты ескереміз:

(9)

(7), (8) - шектік есептердегі туындыларды сәйкес айырымдық қатынастармен алмастырамыз:



бұдан келесі есепті аламыз:



(10)

және
(11)

(10)-(11)-жүйелерді шешіп болған соң (9)-формуламен с тұрақтысы табылса теңдеуін шешу арқылы ізделінді мәндерді табуға болады.

Егер шекаралық шарттар түрінде берілсе, болады да есептеу формулалары жеңілдейді:





болады.

2. Қуалау әдісі.

(3)-(4)-шектік есептерді ақырлы-айырымдық қатынастармен алмастырған соң алынған есепті келесі түрде қарастырайық:



Мұндағы


mi=-2+hpi , ki=1-hpi+h2qi (i=0,1,2,…,n-2) (16)
Сонда жүйе мынадай түрге келеді:

Yi+1=ci (di-yi+2) (i=0,1,2,…,n-2) (17)



сi,di сандары i=0 болғанда

(18)

i=1,2,…,n-2 болғанда біртіндеп

ci=1/(mi-kici-1) (19)

di=fih2-kici-1di-1

формулалармен есептеледі.

Есептеу этаптары:



Тура жол. (16)-формуламен mi,ki мәндерін табамыз. Сосын (18)-формуламен c0, d0 –дарды және (19)-формуламен i=1,2,…,n-2 үшін біртіндеп ci, di мәндерін табамыз.

Кері жол. (17)-формуламен i=n-2 болғанда және (3)-(4)-теңдеулер жүйесінің соңғы теңдеуінен

Yn-1=cn-2(dn-2-yn)



екенін анықтаймыз. Бұл жүйені yn бойынша шешеміз:

(20)

Алдын ала табылған cn-2, dn-2 мәндерін қолданып yn мәнін табамыз. Сосын біртіндеп (i=n-1,…,1) үшін (17)- формуланы қолданып yi мәндерін есептейміз:



(21)

Сосын y0 мәнін (3)-(4)-теңдеулердің ең соңғысының алдындағы теңдеуінен табамыз:



(22)

Сонымен, барлық есептеулер екі рет қуаланып шығады. Тура жолда i индексінің өсу ретімен көмекші ci, di сандары алынады. Бұл арада c0, d0 сандарын есептеу үшін интегралдау аралығының сол жақ шетіндегі шекаралық шарт қолданылады. Сосын кері жолдың алғашқы қадамында есептелген cn-2, dn-2 сандарының интегралдау аралығының оң жақ шетіндегі шекаралық шартпен үйлестірілуін қадағалайды. Осыдан барып I индексінің кему ретінмен біртіндеп yi ізделінді мәндер есептеледі.





Достарыңызбен бөлісу:
  1   2   3




©emirsaba.org 2024
әкімшілігінің қараңыз

    Басты бет