Қазақстан республикасының білім және ғылым министрлігі
М.Әуезов атындағы Оңтүстік Қазақстан университеті
Зертханалық жұмыс
Тақырыбы: Қарапайым дифференциальдық теңдеулерге қойылатын Коши есебін шешудің сандық әдістері. Эйлер әдісі.
Орындаған:Турабаев О
Тобы:ИП18-6К2
Қабылдаған: Досанова Г
Лабораториялық жұмыс №6
Тақырыбы: Қарапайым дифференциальдық теңдеулерге қойылатын Коши есебін шешудің сандық әдістері. Эйлер әдісі.
Теориялық материал
Бірінші ретті қарапайым дифференциялдық теңдеулер келесі түрде беріледі:
y′=f(x,y) (1)
Осы теңдікпен байланысқан негізгі тапсырма бізге Коши есебі деп танымал: Бастапқы
y(x0)=y0 (2)
шартын қанағаттандыратын y=y(x) функциясы түрінде (1) теңдеуінің шешімін табу. (1) теңдігі орындалған кезде берілген M(x0, y0) нүктесі арқылы өтетін интегралдық қисық y(x)-ті табу керек. (1) теңдігінің шешімінің табылуы және жалғыз екендігі келесі теоремамен қаматамасыз етіледі.
Пикар теоремасы: кезкелген G-облысында A -функциясы анықталған және үздіксіз болса, теңсіздіктермен анықталады.
[x-x0]≤a, [y-y0]≤b (3)
Және осы облыста шартты қанағаттандыратын Липщица y-бойынша
[f(x1y1)-f(x1 y2)] ≤M[y1-y2]
бұл арақашықтықта [x-x0]≤h, мұндағы h -оң сан. (7.1)- теңдіктің y=y(x) -тің шешімі тек жалғыз y0=y(x0) бастапқы шартын қанағаттандыратын мұнда M - тұрақты, егер f(x,y) G-облыстарындағы туындының шегі A'y(x,y) болады, онда (x,y)€G жатады.
M=max [f'y(x,y)] (4)
Классикалық анализдерде көптеген дифференциялдық теңдеулерді жазудың әдістері бар. (қарапайым немесе күрделі).
Бұған қарамастан осы есепті шығарар кезінде осы әдістер керексіз болып қалады, немесе оларды шешімі жай уақыт кетірумен өткізеді.
Осы себепке байланысты есепті шығару үшін дифференциялдық теңдеудің жуықтап шешу әдісі қарастырылған. Шешімді өрнектеу түріне байланысты бұл әдістер 3 топқа бөлінеді:
Аналитикалық әдіс. Бұл дифференциялдық теңдеудің шешімі аналитикалық түрде беріледі.
Графикалық әдіс. Жуықтап шешудің графикалық түрде берілуі.
Сандық әдіс. Алынған функция кесте түрінде беріледі. Сандық есептің қателігін дифференциялдық теңдеулерде бірінші ретті, яғни (1) формуламен және n-ші ретті дифференциялдық теңдеулер келесідегідей.
Y(n) =f (x,y,y'…y(n-1)),
Коши есебі y=y(x) шешімін табумен тұрады.
Y(x0)=y(x0)=y0, y'(x0)=y'0.., y(n-1)(x0)=y0(n-1)
Мұндағы y0,y0',…y0(n-1) есептелетін яғни белгісіз сан. Көп жағдайларда дифференциялдық теңдеулердің системасы бірінші ретті теңдеумен шешіледі. y''=(x,y,y') Екінші ретті теңдеулер бірінші ретті теңдеулер системасында да жатады.
y'=z y'=f(x,y,z).
Біртіндеп жуықтау әдісі.
Коши есебін қарастырайық. Біртіндеп жуықтау әдісі бойынша шешімі функциясынын тізбектерінін шегі ретінде қарастырылады. Жоғарыда айтылған шарттар қанағаттандырылсын деп ұйғарсақ, келесі рекуренттік формула бойынша табылады
Эйлер әдісі. Коши есебін карастырайык. Эйлер әдісінің негізгі формуласы
мұндағы h-қадам,
.
Эйлер әдісінің бағалау қателігі
мұндағы –теңдеудін шешімінің дәл мәні (болғанда),
-ші қадамдағы алынған жуық мән.
Рунге-кутт әдісі. Коши есебін қарастырайық. арқылы ізделінді шешімнін нүктесіндегі жуык мәнін белгілейік. Рунге-Кутт әдісі бойынша ізделінді функцяның мәндер тізбегін есептеу формуласы
мұндағы
Бағалау қателігі
мұндағы – нүктесіндегі берілген теңдеудін шешімінін дәл мәні,
- мәндері h/2 және h қадаммен алынған жуық мәндер.
№ 9. , h=0.1
1-қадам
2-қадам
3-қадам
4-қадам
5-қадам
6-қадам
7-қадам
8-қадам
9-қадам
10-қадам
Яғни у(2,7)=5.5463 –қа тең болады.
Достарыңызбен бөлісу: |