Адамстың экстраполяция әдісі.
(6.11)
(6.12)
(6.11)-теңдеу екі өлшемді қарапайым дифференциалдық теңдеу (ҚДТ) және (6.12)-бастапқы шарт берілсін. [x0, xn] аралығында у-тің мәндерін анықтап, функция графигін сызу керек болсын. Бұл есепті шешудің көпқадамды әдістері: Адамс және Милн әдістері деп аталады.
Адамс әдісі
Адамс әдісінің идеясы бірқадамды әдіспен табылған мәндер кесте сын толықтыру немесе жалғастыру. Сондықтан есептеің берілгенінде бастапқы шартпен бірге бірнеше нүктедегі функция мәндері табылған болады.
Бастапқы шартты пайдаланып, функция өсімшесінің мәндерін анықтаймыз, оларды qi (i=0,1,2,3) деп белгілейік:
Енді осы мәндердің шектік айырымдарын табамыз:
,
Енді функцияның мәндерін есептеу үшін Адамс формуласын қолданамыз. Ол екі түрлі:
1 Экстраполяциялық формула:
,
, (6.13)
Бұл формуламен табылған мәндерін алдын ала анықталған функция мәндері деп атаймыз және , деп белгілейміз. (6.13)-формуламен табылған мәндерді тереңірек анықтау үшін интерполяциялық формуланы қолданамыз.
2 Интерполяциялық формула :
,
, (6.14)
Бұл формуламен табылған мәндерін жөнделген немесе түзетілген функция мәндері деп атаймыз және , деп белгілейміз. Сосын (6.13) және (6.14)-формулалармен алынған мәндерді бір бірімен салыстырамыз. Егер төмендегі шарт орындалса:
онда әдіс өзінің жалғыз шешіміне жинақталады, орындалмаса – қадамды кішірейтіп есептеуді қайта жүргізу керек.
Практикада есептеуді жеілдету үшін Адамстың басқа формулалары да қолданылады:
1-формуласы: ), (6.15)
2-формуласы: ), (6.16)
Адамс әдісі дифференциалдық теңдеулер жүйесіне де қолданылады:
жүйесі берілсе, оған қолданылатын Адамс формулалары төмендегідей болады:
(6.17)
(6.18)
Мұндағы: ,
1-мысал:
Адамс әдісін қолданып теңдеуін шешу. Бастапқы шарты y(0)=-1
Шешімі:
(6.15)-(6.16)-формулаларды қолданып есептейік. Рунге-Кутта әдісімен алдын ала бірнеше мәндер табылған болсын.
X1=0.1 y1=-0.97528
X2=0.2 y2=-0.94978
X3=0.3 y3=-0.92154
Есептеу қадамдарын 18-кестеге жазуға болады.
Кестені толтыру ережесіне тоқталайық:
1,2- бағандарға белгілі мәндерді толтырамыз. 3-бағанда уk-дің (k=0,1,2,3)
белгілі мәндерін толтырамыз. Осы мәндерді қолданып,
18-кесте . теңдеуін шешудің алгоритмі.
K
|
xk
|
yk
|
Y`k
|
|
|
|
|
0
|
0.0
|
-0
|
0.25
|
|
|
|
|
1
|
0.1
|
-0.97528
|
0.24779
|
|
|
|
|
2
|
0.2
|
-0.94978
|
0.26552
|
|
|
|
|
3
|
0.3
|
-0.92154
|
0.30232
|
0.32834
|
0.32840
|
0.03283
|
0.03284
|
4
|
0.4
|
-0.88871
|
0.35745
|
0.39237
|
0.39246
|
0.03924
|
0.03925
|
|
|
-0.88870
|
|
|
|
|
|
5
|
0.5
|
-0.84946
|
0.43040
|
|
|
|
|
|
|
-0.84946
|
|
|
|
|
|
формуласымен 3-бағандағы сәйкес мәндерді анықтаймыз. 5-6 – бағандардағы белгілеулері (6.15)- (6.16)- формулалардағы жақша ішіндегі қосындыны білдіреді. K=3 болғанда мәнін 5-бағанның сәйкес жолына жазамыз. K=4 болғанда (6.15)-формуламен мәнін тауып, кестеде өз орнына жазамыз. Осы табылған х4, у4 мәндерін қолданып табамыз. мәнін есептеп өз орнына жазамыз. Әрі қарай K=3 болғанда мәнін есептеп өз орнына жазамыз. Енді табылған у4 мәнін (6.16)- формуламен түзетеміз: . Табылған екі мән бір біріне өте жуық болғандықтан кестедегі алдыңғы табылған мәнін мәнімен түзетеміз. Осы әдіспен кестенің келесі жолын толтыруға болады.
Достарыңызбен бөлісу: |