(7.47)
(7.48)
есебі берілсін. Бұл есептің шешімін келесі түрде іздейміз:
(7.49)
Мұндағы z(x) функциясы - шекаралық шарттары біртекті болатын біртекті емес теңдеудің шешімі:
(7.50)
Ал z1(x) функциясы – шекаралық шарттары біртекті емес болатын біртекті теңдеудің шешімі:
(7.51)
Ал z2(x) функциясы – шекаралық шарттары біртекті емес болатын біртекті теңдеудің шешімі болады:
(7.52)
(7.50)-(7.52)- есептерді шешіп, (7.49)-формула мен (7.48)-шекаралық шарттарды қолданып с1 және с2 коэффициенттерін анықтайтын жүйені құруға болады.
(7.53)
(7.53)-жүйені шығару арқылы у(х) белгісізін тауып, оны (7.47)-(7.48)-есепке қою арқылы шектік есепті шешеміз.
Мысалы:
(7.54)
Шешуі:
1. Біртекті шекаралық шартты біртексіз теңдеу құрамыз:
(7.55)
2.Біртексіз шекаралық шартты біртекті теңдеу құрамыз:
(7.56)
42
3.Біртексіз шекаралық шартты біртекті теңдеу құрамыз:
(7.57)
4. (7.55)-(7.57)-есептерді 1-ретті дифференциалдық теңдеулер жүйесі үшін Коши есебіне келтіреміз:
T(x)=z(x), белгілеулерін енгізейік. Сонда
(7.58)
T1(x)=z1(x), белгілеулерін енгізейік. Сонда
(7.59)
T2(x)=z2(x), белгілеулерін енгізейік. Сонда
(7.60)
есептері құрылады. Бұл есептердегі Т(х), T1(x), T2(x) функциялары (7.55-(7.57)-есептердің шешіимдері болады, ал функциялары осы шешімдердің туындылары болады.
5. Т(х), T1(x), T2(x) функциялары табылған соң (7.53)-формула бойынша с коэффициенттерін табуға арналған жүйе құрамыз:
екенін ескерсек, болады.
екенін ескеріп, (7.53)-ке қойсақ: сонда болады. Келесі теңдеуіне қойсақ:
6. с2-ні табу үшін алдымен (7.58)-(7.60)-есептерді Эйлер әдісімен шешіп Т(хk), T1(xk ), T2(xk ), мәндер кесте сын құру керек. Сосын с2-ні тауып формуласына қою арқылы у(хk) шешімдерін тауып алуға болады.
43
Галеркин әдісі.
Вариациялық, ақырлы-айырымдық әдістерді қолдану барысында шектік есептің шешімі кесте лық мәндер түрінде алынады. Ал Галеркин әдісі аналитикалық әдіске жатады.
y // +p(x) y / + q(x) y = f(x) (7.61)
(7.62)
мұндағы p(x), q(x), f(x) – функциялары [a,b] аралығында үзіліссіз функциялар, ал 0,1,0,1, A,B –берілген тұрақтылар және 0+10, 010 шектік есебі берілсін.
Белгілеу енгізейік:
L[y]=y // +p(x) y / + q(x) y (7.63)
(7.64)
[a,b] аралығында базистік функциялар жүйесі берілсін
(7.65)
және олар келесі шарттарды қанағаттандырсын:
(7.65)-жүйе ортогональды болсын, яғни
(7.66)
(7.65)-жүйе толық жүйе болсын, яғни барлық функцияларына ортогональ нөлден өзгеше басқа бір де бір функция жоқ болсын.
ақырлы базистік функциялар жүйесі төмендегі біртексіз (7.67)-шарттарды қанағаттандыратындай болып, ал функциялары төмендегі біртекті (7.68)-шарттарды қанағаттандыратындай болып таңдалынып алынуы керек:
(7.67)
(7.68)
(7.61)-(7.62)-шектік есептің шешімін
(7.69)
түрінде іздейміз. (7.67), (7.68)-шарттардан бұл функцияның (7.62)-шартты қанағаттандыратыны көрінеді.
Байныспаушылық (невязка) өрнегін қарастырамыз:
(7.70)
сi коэффициенттерін байналыспаушылықтың квадратының интегралының мәні
(7.71)
өте аз шама болатындай етіп таңдап алу керек. Бұл тек қана байланыспаушылық базистік функциялардың барлығына ортогональды болса ғана орындалатыны дәлелденген ([2], [13], [18] қараңыз).
44
Ортогональдылық шарты:
немесе
(7.72)
(7.72)-ден сызықты алгебралық теңдеулер жүйесін аламыз, оны сi-ға қатысты шешу керек болады. Мұндай жүйені шешудің сандық әдістерін жоғарыда келтіргенбіз.
Мысалы:
Галеркин әдісін қолданып
(7.73)
у(0)=у(1)=0 (7.74)
шектік есептің шешімін табу.
Ui базистік функциялар жүйесін таңдап аламыз:
Бұл функциялар сызықты тәуелсіз және нөлдік шекаралық шарттарды қанағаттандырады. Жуық шешімін мына түрде іздейміз:
(7.75)
Бұл теңдеуді (7.73)-теңдеудің сол жағына қою арқылы байланыспаушылықты аламыз:
. (7.76)
R функциясының u1(x) және u2(x) функцияларына ортогональдылығын ескерсек:
Бұл жүйеге (7.77) – шешімді қойсақ және интегралды есептесек, келесі сызықты алгебралық теңдеулер жүйесін аламыз:
Бұл жүйені шешсек . Бұл мәндерді (7.76)-ға қойсақ:
.
Бұл есептің дәл шешімі . Берілген аралықты 0,25 қадаммен бірнеше бөлікке бөліп әр нүктедегі функция мәндерін алып салыстырсақ, төмендегі 22-кестеден қателігінің қаншалықты екенін анықтауға болады.
22-кесте. (7.73-7.74)-есептің мәндер кестесi.
xi
|
0.25
|
0.50
|
0.75
|
yi
|
0.044
|
0.069
|
0.060
|
Y
|
0.044
|
0.070
|
0.050
|
45
Егер (7.61)-теңдеудің коэффициенттерін есептеу қиынға соғатын болса бұл әдістен гөрі ақырлы-айырымдық немесе қуалау әдістерінің бірін қолдануға болады.
Достарыңызбен бөлісу: |