Әдістемелік жинақ
ҚДТ-ді шешудің көпқадамды сандық әдістері
жүктеу/скачать 2,2 Mb.
|
| ҚДТ-ді шешудің көпқадамды сандық әдістері.
Адамс әдісі Милн әдісі Дәріс тезисі: (6.11) (6.12) (6.11)-теңдеу екі өлшемді қарапайым дифференциалдық теңдеу (ҚДТ) және (6.12)-бастапқы шарт берілсін. [x0, xn] аралығында у-тің мәндерін анықтап, функция графигін сызу керек болсын. Бұл есепті шешудің көпқадамды әдістері: Адамс және Милн әдістері деп аталады. Адамс әдісі Адамс әдісінің идеясы бірқадамды әдіспен табылған мәндер кесте сын толықтыру немесе жалғастыру. Сондықтан есептеің берілгенінде бастапқы шартпен бірге бірнеше нүктедегі функция мәндері табылған болады. Бастапқы шартты пайдаланып, функция өсімшесінің мәндерін анықтаймыз, оларды qi (i=0,1,2,3) деп белгілейік: 32
, Енді функцияның мәндерін есептеу үшін Адамс формуласын қолданамыз. Ол екі түрлі: 1 Экстраполяциялық формула: , , (6.13) Бұл формуламен табылған мәндерін алдын ала анықталған функция мәндері деп атаймыз және , деп белгілейміз. (6.13)-формуламен табылған мәндерді тереңірек анықтау үшін интерполяциялық формуланы қолданамыз. 2 Интерполяциялық формула : , , (6.14) Бұл формуламен табылған мәндерін жөнделген немесе түзетілген функция мәндері деп атаймыз және , деп белгілейміз. Сосын (6.13) және (6.14)-формулалармен алынған мәндерді бір бірімен салыстырамыз. Егер төмендегі шарт орындалса: онда әдіс өзінің жалғыз шешіміне жинақталады, орындалмаса – қадамды кішірейтіп есептеуді қайта жүргізу керек. Практикада есептеуді жеілдету үшін Адамстың басқа формулалары да қолданылады: 1-формуласы: ), (6.15) 2-формуласы: ), (6.16) Адамс әдісі дифференциалдық теңдеулер жүйесіне де қолданылады: жүйесі берілсе, оған қолданылатын Адамс формулалары төмендегідей болады: (6.17) (6.18) Мұндағы: , Милн әдісі Бұл әдіс те Адамс әдісі сияқты мәндер кесте сын жалғастыруға мүмкіндік береді. Теңдеу, бастапқы шарт, және қандай да бір әдіспен табылған функцияның бірнеше мәндері берілсін. Функцияның қалған мәндерін анықтау керек. 33 (6.19) (6.20) , мәндерді анықтау үшін Милн формулаларын қолданамыз: Алдын ала анықтау: , Осы мәндерді қолданып Милннің 2-ші формуласымен алдында табылған мәндерді түзетеміз немесе дәлдейміз: . Табылған мәндердің қателігін бағалаймыз: . Бұл формула есептеудің әр қадамында алынған мәннің дәлдігін тексеріп отырады. Егер дәлдік берілсе және болса, онда деп алып yi+1-лерді есептеуге болады, кері жағдайда қадамды кішірейту керек. Милн әдісін жүйені шешуге де қолдануға болады. Егер жоғарғы ретті теңдеу берілсе оны 1-ші ретті теңдеуге келтіру керек. 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-кесте . теңдеуін шешудің алгоритмі.
формуласымен 3-бағандағы сәйкес мәндерді анықтаймыз. 5-6 – бағандардағы белгілеулері (6.15)- (6.16)- формулалардағы жақша ішіндегі қосындыны білдіреді. K=3 34 болғанда мәнін 5-бағанның сәйкес жолына жазамыз. K=4 болғанда (6.15)-формуламен мәнін тауып, кестеде өз орнына жазамыз. Осы табылған х4, у4 мәндерін қолданып табамыз. мәнін есептеп өз орнына жазамыз. Әрі қарай K=3 болғанда мәнін есептеп өз орнына жазамыз. Енді табылған у4 мәнін (6.16)- формуламен түзетеміз: . Табылған екі мән бір біріне өте жуық болғандықтан кестедегі алдыңғы табылған мәнін мәнімен түзетеміз. Осы әдіспен кестенің келесі жолын толтыруға болады. жүктеу/скачать 2,2 Mb. Достарыңызбен бөлісу:
|