1. 1Differensial tenglamalar haqida boshlangich tushunchalar

Ta’rif. F(x,y,y’,....,y⁽ⁿ⁾)=0 ko’rinishdagi tenglamaga n - tartibli differensial tenglama deyiladi.
Ta’rif. n - tartibli differensial tenglamaning umumiy yechimi deb n ta с₁, с₂, .... с_n - ixtiyoriy o’zgarmas miqdorlarga bog’liq bo’lgan
y= (x, с₁, с₂, .... с_n)
funksiyaga aytiladi. Bu funksiya:
с₁,...,с_n larning ixtiyoriy qiymatlarida tenglamani qanoatlantiradi;berilgan y(x₀)=y₀, (x₀)=y₁,..., y^(n-1)(x₀)=y_n-1 boshlang’ich shartda с₁, с₂, .... с_n larni shunday tanlash mumkinki,
y= (x, с₁, с₂, .... с_n) funksiya bu boshlang’ich shartni qanoatlantiradi.
Ta’rif. Umumiy yechimdan с₁, с₂, .... с_n miqdorlarning tayin qiymatlarida hosil bo’ladigan funksiya xususiy yechim deyiladi. Birinchi tartibli differensial tenglamalarni yechishning sonli usullari. Differensial tenglamalarni yechishning sonli usullari. Oddiy differensial tenglamalarni yechish
Raqamli yechim differensial tenglamalar
Fan va texnikaning ko'pgina muammolari oddiy differensial tenglamalarni (ODE) echishga qisqartiriladi. ODElar - kerakli funktsiyaning bir yoki bir nechta hosilalarini o'z ichiga olgan tenglamalar. Umuman olganda, ODE quyidagicha yozilishi mumkin:Bu erda x mustaqil o'zgaruvchi, talab qilinadigan funktsiyaning i-hosilasi. n - tenglamaning tartibi. n-tartibli ODE ning umumiy yechimi n ta ixtiyoriy konstantadan iborat, ya’ni umumiy yechim.Bitta yechim tanlash uchun n ta qo'shimcha shartni ko'rsatish kerak. Qo'shimcha shartlarni belgilash usuliga qarab ikki xil turdagi masalalar mavjud: Koshi muammosi va chegaraviy masala. Agar bir nuqtada qo'shimcha shartlar ko'rsatilgan bo'lsa, unda bunday muammo Koshi muammosi deb ataladi. Koshi masalasidagi qo'shimcha shartlar boshlang'ich shartlar deyiladi. Agar bir nechta nuqtada qo'shimcha shartlar ko'rsatilgan bo'lsa, ya'ni. mustaqil o'zgaruvchining turli qiymatlari uchun bunday muammo chegaraviy masala deb ataladi. Qo'shimcha shartlarning o'zi chegara yoki chegara shartlari deb ataladi.

Ko'rinib turibdiki, n = 1 uchun biz faqat Koshi muammosi haqida gapirishimiz mumkin.

Koshi muammosini o'rnatishga misollar:
Chegaraviy masalalarga misollar:
Bunday masalalarni faqat ayrim maxsus turdagi tenglamalar uchun analitik tarzda yechish mumkin.Birinchi tartibli ODElar uchun Koshi masalasini echishning raqamli usullari Muammoni shakllantirish... Birinchi tartibli ODE yechimini toping
Taqdim etilgan segmentda Taxminiy yechim topilganda, hisob-kitoblar hisoblangan qadam bilan amalga oshirilgan deb hisoblaymiz, hisoblangan tugunlar oraliq nuqtalari [ x 0 , x n ].Shu nuqtada haqiqatdan foydalanib x 0 ma'lum yechim y(x 0)= y 0 va uning hosilasi qiymati, siz nuqtada kerakli funktsiya grafigiga teginish tenglamasini yozishingiz mumkin:. Etarlicha kichik qadam bilan h qiymatning o'ng tomoniga almashtirish orqali olingan bu tangensning ordinatasi ordinatadan ozgina farq qilishi kerak y(x 1) yechimlar y(x) Koshi muammosi. Shuning uchun tangensning to'g'ri chiziq bilan kesishish nuqtasi x = x 1 taxminan yangi boshlanish nuqtasi sifatida qabul qilinishi mumkin. Ushbu nuqta orqali yana to'g'ri chiziq torting, bu nuqtaga tegishning harakatini taxminan aks ettiradi. Bu yerni almashtirish (ya'ni, chiziq bilan kesishish x = x 2), biz taxminiy qiymatni olamiz y(x) nuqtada x 2: va boshqalar. Natijada, uchun i- Shu nuqtada Eyler formulasini Bu usul deyiladi yashirin Eyler usuli, chunki ma'lum qiymatdan noma'lum qiymatni hisoblash uchun umumiy holatda chiziqli bo'lmagan tenglamani echish kerak.
Yashirin Eyler usuli aniqlik yoki yaqinlashishning birinchi darajasiga ega.
Ushbu usulda hisoblash ikki bosqichdan iborat:
Bu sxema bashoratchi-tuzatuvchi usul (bashoratchi-tuzatuvchi) deb ham ataladi. Birinchi bosqichda taxminiy qiymat past aniqlik (h) bilan bashorat qilinadi, ikkinchi bosqichda esa bu bashorat tuzatiladi, natijada olingan qiymat ikkinchi aniqlik tartibiga ega bo'ladi.Runge - Kutta usullari: aniq Runge-Kutta usullarini yaratish g'oyasi p- Bu buyurtma qiymatlarga yaqinliklarni olishdir y(x i+1) shakl formulasi bo'yicha Bu yerda a n , b nj , p n, - ba'zi sobit raqamlar (parametrlar).Runge – Kutta usullarini qurishda funksiya parametrlari ( a n , b nj , p n) kerakli yaqinlashish tartibini oladigan tarzda tanlanadi.Runge - to'rtinchi darajadagi aniqlikning Kutta sxemasi:Misol... Koshi muammosini hal qiling:Uchta usulni ko'rib chiqing: aniq Eyler usuli, modifikatsiyalangan Eyler usuli va Runge - Kutta usuli.Aniq yechim:Ushbu misol uchun aniq Eyler usulidan foydalangan holda hisoblash formulalari:O'zgartirilgan Eyler usulining hisoblash formulalari:Runge - Kutta usulini hisoblash formulalari:

y1 - Eyler usuli, y2 - modifikatsiyalangan Eyler usuli, y3 - Runge Kutta usuli.

Ko'rinib turibdiki, eng aniq Runge - Kutta usuli hisoblanadi.

Birinchi tartibli ODE tizimlarini echishning raqamli usullari

Ko'rib chiqilgan usullardan birinchi tartibli differensial tenglamalar tizimini yechishda ham foydalanish mumkin.

Keling, buni ikkita birinchi tartibli tenglamalar tizimi uchun ko'rsatamiz:

Aniq Eyler usuli:

O'zgartirilgan Eyler usuli:

Runge - to'rtinchi darajadagi aniqlikdagi Kutta sxemasi:

Yuqori tartibli tenglamalar uchun Koshi masalalari ODE tenglamalar tizimini yechish uchun ham qisqartiriladi. Masalan, ko'rib chiqing ikkinchi tartibli tenglama uchun Koshi muammosi

Ikkinchi noma'lum funktsiyani kiritamiz. Keyin Cauchy muammosi quyidagilar bilan almashtiriladi:

Bular. oldingi vazifa bo'yicha:.

Misol. Koshi muammosiga yechim toping:

Segmentda.

Aniq yechim:

Haqiqatan ham:

Eyler va Runge-Kutta usulida h = 0,2 qadam bilan o'zgartirilgan aniq Eyler usuli yordamida masalani hal qilaylik.

Funktsiyani tanishtiramiz.

Keyin ikkita birinchi tartibli ODE tizimi uchun quyidagi Koshi muammosini olamiz:

Aniq Eyler usuli:

O'zgartirilgan Eyler usuli:

Runge-Kutta usuli:

Eyler sxemasi:

O'zgartirilgan Eyler usuli:

Runge - Kutta sxemasi:

Maks (y-y nazariyasi) = 4 * 10 -5

ODE uchun chegaraviy masalalarni echishning chekli farq usuli

Muammoni shakllantirish: chiziqli differensial tenglamaning yechimini toping

chegara shartlarini qondirish: (2)

Teorema. Mayli. Keyin muammoning o'ziga xos echimi bor.

Bu muammo, masalan, uchlarida mentli bo'lgan nurning egilishlarini aniqlash muammosi kamayadi.

Cheklangan farq usulining asosiy bosqichlari:

1) argumentning uzluksiz o'zgaruvchanligi mintaqasi () tugunlar deb ataladigan diskret nuqtalar to'plami bilan almashtiriladi:.

2) Uzluksiz argument x ning talab qilinadigan funktsiyasi taxminan berilgan to'rdagi diskret argumentning funktsiyasi bilan almashtiriladi, ya'ni. ... Funktsiya grid deb ataladi.

3) Asl differensial tenglama panjara funksiyasiga nisbatan ayirma tenglama bilan almashtiriladi. Bu o'zgarish farqni yaqinlashish deyiladi.

Shunday qilib, differensial tenglamaning yechimi algebraik tenglamalar yechimidan topilgan to'r tugunlarida to'r funksiyasining qiymatlarini topishga qisqartiriladi.

Hosilalarni yaqinlashtirish.

Birinchi hosilani taxmin qilish (almashtirish) uchun siz quyidagi formulalardan foydalanishingiz mumkin:

- to'g'ri farq hosilasi,

- chap farq hosilasi,

Markaziy farq hosilasi.

ya'ni hosilani taxminan aniqlashning ko'plab usullari mavjud.

Bu ta'riflarning barchasi lotinning chegara sifatidagi tushunchasidan kelib chiqadi: .

Birinchi hosilaning ayirma yaqinlashuviga asoslanib, ikkinchi hosilaning ayirma yaqinlashuvini qurish mumkin:

Xuddi shunday, yuqori tartibli hosilalar uchun taxminiy ma'lumotlarni olish mumkin.

Ta'rif. Farqi n-chi hosilaning yaqinlashish xatosi deyiladi:.

Teylor kengayishi yaqinlashish tartibini aniqlash uchun ishlatiladi.

Birinchi hosilaning o'ng tomonidagi farqning yaqinlashuvini ko'rib chiqing:

Bular. to'g'ri farq hosilasi bor birinchi bo'lib h yaqinlashish tartibi.

Xuddi shunday chap farq hosilasi uchun.

Markaziy farq hosilasi bor ikkinchi tartibli yaqinlashish.

(3) formula bo'yicha ikkinchi hosilaning yaqinlashuvi ham ikkinchi darajali yaqinlashish tartibiga ega.

Differensial tenglamani yaqinlashtirish uchun barcha hosilalarni ularning yaqinliklari bilan almashtirish kerak. (1), (2) muammoni ko'rib chiqing va (1) dagi hosilalarni almashtiring:

Natijada biz quyidagilarni olamiz:

(4)

Dastlabki masalani yaqinlashtirish tartibi 2 ga teng, chunki ikkinchi va birinchi hosilalar 2-tartib bilan almashtiriladi, qolganlari esa aynan.

Shunday qilib, (1), (2) differensial tenglamalar o'rniga biz tizimni oldik chiziqli tenglamalar panjara nuqtalarida aniqlash uchun.

Sxema quyidagicha ifodalanishi mumkin:

ya'ni matritsali chiziqli tenglamalar tizimini oldik:

Ushbu matritsa tridiagonal, ya'ni. asosiy diagonalda va ikkita qo'shni diagonalda joylashgan bo'lmagan barcha elementlar nolga teng.

Hosil bo'lgan tenglamalar tizimini yechish orqali biz dastlabki masala yechimini olamiz.

Differensial tenglamalar noma’lum funksiya hosila belgisi ostida kiradigan tenglamalardir. Differensial tenglamalar nazariyasining asosiy vazifasi shunday tenglamalarning yechimi bo'lgan funksiyalarni o'rganishdir.

Differensial tenglamalarni noma’lum funksiyalari bitta o‘zgaruvchining funksiyasi bo‘lgan oddiy differensial tenglamalarga va noma’lum funksiyalari ikki va bir o‘zgaruvchining funksiyasi bo‘lgan qisman differentsial tenglamalarga bo‘lish mumkin. Ko'proq o'zgaruvchilar.

Qisman differensial tenglamalar nazariyasi murakkabroq boʻlib, u toʻliqroq yoki maxsus matematika kurslarida koʻrib chiqiladi.

Differensial tenglamalarni o'rganishni eng oddiy tenglama - birinchi tartibli tenglama bilan boshlaylik.

Shakl tenglamasi

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

bu erda x - mustaqil o'zgaruvchi; y - talab qilinadigan funktsiya; y "- uning hosilasi, birinchi tartibli differensial tenglama deyiladi.

Agar (1) tenglamani y "ga nisbatan yechish mumkin bo'lsa, u holda u shaklni oladi

va hosilaga nisbatan yechilgan birinchi tartibli tenglama deyiladi.

Ayrim hollarda (2) tenglamani f (x, y) dx - dy = 0 ko'rinishida yozish qulayroqdir, bu umumiyroq tenglamaning alohida holatidir.

P (x, y) dx + Q (x, y) dy = O, (3)

Bu erda P (x, y) va Q (x, y) ma'lum funktsiyalardir. Simmetrik shakldagi tenglama (3) qulay, chunki undagi x va y o'zgaruvchilar tengdir, ya'ni ularning har biri ikkinchisining funksiyasi sifatida qaralishi mumkin.

Keling, tenglamaning umumiy va xususiy yechimining ikkita asosiy ta'rifini beraylik.

(2) tenglamaning umumiy yechimi Oksi tekislikning ba'zi G hududida x va ixtiyoriy doimiy C ga bog'liq bo'lgan y = q (x, C) funktsiya, agar u (2) tenglamaning yechimi bo'lsa. C doimiysining istalgan qiymati va agar (x 0; y 0) = G bo'lgan har qanday boshlang'ich shartlar uchun yx = x0 = y 0 bo'lsa, C = C 0 doimiysining yagona qiymati mavjud bo'lib, funktsiya y = q bo'ladi. (x, C 0) berilgan dastlabki shartlarni y = q (x 0, C) qanoatlantiradi.

(2) tenglamaning G sohasidagi qisman yechimi y = q (x, C 0) funksiya bo lib, u S = C konstantasining ma lum qiymatida y = q (x, C) umumiy yechimdan olinadi. 0.

Geometrik jihatdan umumiy yechim y = q (x, C) Oksi tekisligidagi bir ixtiyoriy doimiy C ga bog‘liq bo‘lgan integral egri chiziqlar turkumidir va xususiy yechim y = q (x, C 0) buning bir integral egri chizig‘idir. oila o'tadi belgilash nuqtasi(x 0; y 0).

Birinchi tartibli differensial tenglamalarni Eyler usulida taqribiy yechish. Bu usulning mohiyati shundan iboratki, ma'lum bir yechimning grafigi bo'lgan talab qilinadigan integral egri chiziq taxminan siniq chiziq bilan almashtiriladi. Differensial tenglama berilsin

va dastlabki shartlar y |x = x0 = y 0.

Berilgan dastlabki shartlarni qanoatlantirgan holda [x 0, b] oraliqda tenglamaning taqribiy yechimi topilsin.

[ x₀ b] segmentini x 0 nuqtalarga ajratamiz<х 1 ,<х 2 <...<х n =b нa n paвных чacтeй. Пycть х 1 --х 0 =х 2 -- x 1 = ... =x n -- x n-1 = ?x. Oбoзнaчиm чepeз y i пpиближeнныe знaчeния иckomoгo peшeния в тoчkaх х i (i=1, 2, ..., n). Пpoвeдem чepeз тoчkи paзбиeния х i - пpяmыe, пapaллeльныe ocи Oy, и пocлeдoвaтeльнo пpoдeлaem cлeдyющиe oднoтипныe oпepaции.

y "= f (x, y)) tenglamaning o'ng tomoniga x 0 va y 0 qiymatlarini qo'ying va integral egri chiziqqa teginishning y" = f (x 0, y 0) qiyaligini hisoblang. nuqta (x 0; y 0). Kerakli yechimning taxminiy qiymati y 1 ni topish uchun [x 0, x 1,] segmentidagi integral egri chiziqni (x 0; y 0) nuqtadagi tangensining segmentiga almashtiring. Bunday holda biz olamiz

y 1 - y 0 = f (x 0; y 0) (x 1 - x 0),

shuning uchun x 0, x 1, y 0 ma'lum bo'lgani uchun topamiz

y1 = y0 + f (x0; y0) (x1 - x0).

y "= f (x, y) tenglamasining o'ng tomoniga x 1 va y 1 qiymatlarini qo'yib, biz integral egri chiziqqa teginishning y" = f (x 1, y 1) qiyaligini hisoblaymiz. nuqta (x 1; y 1). Keyinchalik, segmentdagi integral egri chiziqni tangens segment bilan almashtirib, x 2 nuqtasida y 2 yechimning taxminiy qiymatini topamiz:

y 2 = y 1 + f (x 1; y 1) (x 2 - x 1)

Bu tenglikda x 1, y 1, x 2 ma'lum va y 2 ular orqali ifodalanadi.

Xuddi shunday, biz topamiz

y 3 = y 2 + f (x 2; y 2)? x,…, y n = y n-1 + f (x n-1; y n-1)? x

Shunday qilib, kerakli integral egri chiziq taxminan siniq chiziq shaklida tuziladi va x i nuqtalarda kerakli eritmaning y i ning taxminiy qiymatlari olinadi. Bunday holda, y i ning qiymatlari formula bo'yicha hisoblanadi

y i = y i-1 + f (x i-1; y i-1)?x (i = 1,2, ..., n).

Formula Eyler usulining asosiy hisoblash formulasi hisoblanadi. Uning aniqligi qanchalik baland bo'lsa, farq shunchalik kichikmi? X.

Eyler usuli kerakli y (x) funktsiyasining taxminiy qiymatlari jadvali ko'rinishida yechim beradigan raqamli usullarni anglatadi. Bu nisbatan qo'pol va birinchi navbatda qo'pol hisoblar uchun ishlatiladi. Biroq, Eyler usuli asosidagi g'oyalar bir qator boshqa usullar uchun boshlang'ich nuqtadir.

Umuman olganda, Eyler usulining aniqlik darajasi past. Differensial tenglamalarni taxminiy yechish uchun ancha aniqroq usullar mavjud.

ODDIY DIFFERENTSIAL TENGLAMALARNI YECHISH

Muammoni shakllantirish

Ilmiy va muhandislik masalalarini hal qilishda ko'pincha dinamik tizimni matematik tarzda tasvirlash kerak bo'ladi. Bu eng yaxshi differentsial tenglamalar shaklida amalga oshiriladi ( DU) yoki differentsial tenglamalar sistemasi. Ko'pincha bunday muammo kimyoviy reaktsiyalarning kinetikasini va turli xil uzatish hodisalarini (issiqlik, massa, impuls) modellashtirish bilan bog'liq muammolarni hal qilishda paydo bo'ladi - issiqlik uzatish, aralashtirish, quritish, adsorbsiya, makro va mikrozarrachalarning harakatini tavsiflashda.

Oddiy differentsial tenglama n tartibli (ODE) quyidagi tenglama bo'lib, u istalgan y (x) funksiyaning bir yoki bir nechta hosilalarini o'z ichiga oladi:

Bu yerda y (n) baʼzi y (x) funksiyaning n tartibli hosilasini bildiradi, x mustaqil oʻzgaruvchidir.

Ayrim hollarda differensial tenglamani eng yuqori hosila aniq shaklda ifodalangan shaklga aylantirish mumkin. Belgilanishning bu shakli tenglama deb ataladi, eng yuqori hosilaga nisbatan ruxsat etiladi(bu holda, eng yuqori hosila tenglamaning o'ng tomonida yo'q):

Aynan shu yozuv shakli sifatida qabul qilinadi standart ODElarni echishning raqamli usullarini ko'rib chiqishda.

Chiziqli differentsial tenglama y (x) funksiya va uning barcha hosilalariga nisbatan chiziqli tenglamadir.

Masalan, quyida birinchi va ikkinchi darajali chiziqli ODElar keltirilgan

Oddiy differensial tenglamani yechish orqali har qanday x uchun ma'lum chekli yoki cheksiz oraliqda bu tenglamani qanoatlantiradigan y (x) funksiya. Differensial tenglamani yechish jarayoni deyiladi differensial tenglamani integrallash orqali.

Umumiy ODE yechimi n-tartibda n ta ixtiyoriy C 1, C 2, ..., C n konstantalar mavjud.

Bu shunisi aniqki, noaniq integral integrandning anti hosilasi va integrasiya doimiysiga teng.

n-tartibli DE ni yechish uchun n ta integrallashni amalga oshirish zarurligi sababli umumiy yechimda n ta integrallash konstantasi paydo bo ladi.

Shaxsiy yechim Agar biz integratsiya konstantalariga ba'zi qo'shimcha shartlarni belgilash orqali ba'zi qiymatlarni belgilasak, ularning soni barcha aniqlanmagan integratsiya konstantalarini hisoblash imkonini beradigan umumiy qiymatdan ODE olinadi.

Aniq (analitik) yechim (umumiy yoki xususiy) differensial tenglama elementar funksiyalar ifodasi shaklida kerakli yechimni (y (x) funksiya) olishni nazarda tutadi. Bu har doim ham, hatto birinchi tartibli tenglamalar uchun ham mumkin emas.

Raqamli yechim DE (bo'lim) y (x) funksiyani va uning ma'lum bir segmentda yotgan ba'zi berilgan nuqtalarda hosilalarini hisoblashdan iborat. Ya'ni, aslida, shaklning n-tartibining yechimi quyidagi raqamlar jadvali ko'rinishida olinadi (eng yuqori hosilaning qiymatlari ustuni qiymatlarni tenglamaga almashtirish orqali hisoblanadi). ):

Misol uchun, birinchi tartibli differensial tenglama uchun yechim jadvali ikkita ustunga ega bo'ladi - x va y.

Funktsiyaning qiymati aniqlanadigan abscissa qiymatlari to'plami deyiladi to'r, bunda y (x) funksiya aniqlangan. Koordinatalarning o'zi deyiladi mesh tugunlari... Ko'pincha, qulaylik uchun ishlatiladi yagona panjaralar, unda qo'shni tugunlar orasidagi farq doimiy va deyiladi panjara qadami yoki integratsiya bosqichi differensial tenglama

Yoki, i= 1, ..., N

Aniqlash uchun shaxsiy yechim integratsiya konstantalarini hisoblash imkonini beradigan qo'shimcha shartlarni belgilash kerak. Bundan tashqari, bunday shartlar aniq bo'lishi kerak. Birinchi tartibli tenglamalar uchun - bitta, ikkinchisi uchun - 2 va boshqalar. Differensial tenglamalarni yechishda qanday qo‘yilishiga ko‘ra, uch xil masalalar mavjud:

· Koshi muammosi (dastlabki muammo): Bundaylarni topish kerak shaxsiy yechim aniqni qanoatlantiradigan differensial tenglama bir nuqtada berilgan dastlabki shartlar:

ya'ni mustaqil o'zgaruvchining o'ziga xos qiymati (x 0) berilgan va bu nuqtada (n-1) tartibgacha bo'lgan funktsiya va uning barcha hosilalari qiymati. Bu nuqta (x 0) deyiladi boshlang'ich... Masalan, agar 1-tartibdagi DE yechilsa, u holda boshlang'ich shartlar juft sonlar (x 0, y 0) sifatida ifodalanadi.

Bunday muammoni hal qilishda duch keladi ODE Bu, masalan, kimyoviy reaksiyalarning kinetikasini tavsiflaydi. Bunday holda, vaqtning dastlabki momentidagi moddalarning kontsentratsiyasi ma'lum ( t = 0), va ma'lum vaqtdan keyin moddalarning konsentratsiyasini topish kerak ( t). Misol tariqasida issiqlik uzatish yoki massa uzatish (diffuziya), kuchlar ta'sirida moddiy nuqtaning harakat tenglamasi va boshqalarni ham keltirish mumkin.

· Chegara muammosi ... Bunday holda, funktsiya va (yoki) uning hosilalari qiymatlari bir nechta nuqtada, masalan, vaqtning boshlang'ich va oxirgi momentida ma'lum bo'ladi va differentsial tenglamaning ma'lum bir yechimini topish kerak. bu nuqtalar orasida. Bu holatda qo'shimcha shartlarning o'zi deyiladi mintaqaviy (chegara chizig'i) shartlar. Tabiiyki, chegaraviy masala kamida ikkinchi tartibli ODE uchun echilishi mumkin. Quyida chegara shartlariga ega ikkinchi tartibli ODE ga misol keltirilgan (funktsiyaning qiymatlari ikki xil nuqtada berilgan):

· Shturm-Liouvil muammosi (o'ziga xos qiymat muammosi). Ushbu turdagi muammolar chegaraviy masalaga o'xshaydi. Ularni hal qilishda har qanday parametrning qaysi qiymatlarida yechim topish kerak DU parametrning har bir qiymati uchun differensial tenglamaning yechimi bo'lgan chegaraviy shartlarni (o'ziga xos qiymatlar) va funktsiyalarni qanoatlantiradi (o'ziga xos funktsiyalar). Masalan, kvant mexanikasining ko'pgina muammolari xususiy qiymat muammolari.

Birinchi tartibli ODElar uchun Koshi masalasini echishning raqamli usullari

Yechishning bir necha raqamli usullarini ko'rib chiqing Cauchy muammolari Birinchi tartibli oddiy differensial tenglamalar (dastlabki masala). Keling, bu tenglamani hosilaga nisbatan yechilgan umumiy shaklda yozamiz (tenglamaning o'ng tomoni birinchi hosilaga bog'liq emas):

(6.2)

Agar boshlang'ich qiymatlari ma'lum bo'lsa, y (x) funktsiyasining qiymati x 0 boshlang'ich nuqtasida mavjud bo'lsa, to'rning berilgan nuqtalarida y funktsiyasining qiymatlarini topish kerak.

Tenglamani d x ga ko'paytirish orqali o'zgartiring

Va biz chap va o'ng tomonlarni to'rning i va i + 1-tugunlari o'rtasida birlashtiramiz.

(6.3)

Biz to'rning i-tugunidagi x va y qiymatlari bo'yicha i + 1 integratsiya tugunida yechim qurish ifodasini oldik. Biroq, qiyinchilik shundan iboratki, o'ng tomondagi integral aniq berilgan funktsiyaning integrali bo'lib, uni umumiy holatda analitik tarzda topib bo'lmaydi. ODE ni turli yo'llar bilan echishning raqamli usullari ODE ning raqamli integratsiyasi uchun formulalarni qurish uchun ushbu integralning qiymatini taxminiy (taxminan) qiladi.

Birinchi tartibli ODElarni echish uchun ishlab chiqilgan ko'plab usullardan usullarni ko'rib chiqing va. Ular juda oddiy va raqamli yechim doirasida ushbu muammoni hal qilishning yondashuvlari haqida dastlabki fikrni beradi.

Eyler usuli

Tarixiy jihatdan birinchi tartibli ODElar uchun Koshi masalasini raqamli yechishning birinchi va eng oddiy usuli Eyler usuli hisoblanadi. Bu bog'liqning ( y) va mustaqil ( x) yagona tarmoq tugunlari orasidagi o'zgaruvchilar:

bu yerda y i + 1 funksiyaning x i + 1 nuqtadagi kerakli qiymati.

Agar biz ushbu tenglamani o'zgartirsak va integratsiya tarmog'ining bir xilligini hisobga olsak, biz iterativ formulaga ega bo'lamiz, bu orqali biz hisoblashimiz mumkin. y i + 1 agar y i x i nuqtada ma'lum bo'lsa:

Eyler formulasini avval olingan umumiy ifoda bilan solishtirsak, integralni Eyler usulida taqribiy hisoblash uchun eng oddiy integrallash formulasi - segmentning chap qirrasi bo ylab to rtburchaklar formulasi qo llanilishini ko rish mumkin.

Eyler usulining grafik talqini ham oddiy (quyidagi rasmga qarang). Darhaqiqat, echilayotgan tenglamaning () shakliga asoslanib, qiymat y (x) funksiyaning x = xi - nuqtasidagi hosilasining qiymati va shuning uchun tangensga teng ekanligi kelib chiqadi. y (x) funksiya grafigiga x = xi nuqtada chizilgan tangensning qiyaligi.

Rasmdagi to'g'ri burchakli uchburchakdan topishingiz mumkin

bu erdan Eyler formulasi olinadi. Demak, Eyler usulining mohiyati integrallash intervalidagi y (x) funksiyani x = x i nuqtadagi grafaga tangens to g ri chiziq bilan almashtirishdan iborat. Agar kerakli funktsiya chiziqli funktsiyadan integratsiya oralig'ida katta farq qilsa, hisoblash xatosi sezilarli bo'ladi. Eyler usulining xatosi integratsiya bosqichiga to'g'ridan-to'g'ri proportsionaldir:

Xato~ h

Hisoblash jarayoni quyidagicha tuzilgan. Dastlabki shartlarni hisobga olgan holda x 0 va y 0 hisoblash mumkin

Shunday qilib, y (x) funktsiyasi qiymatlari jadvali ma'lum bir qadam bilan tuziladi ( h) yoqilgan x segmentida. Qiymatni aniqlashda xato y (x i) bu holda, u kamroq bo'ladi, qadam uz — Источник: https://muegn.ru/uz/okazanie-pomoschi/chislennye-metody-resheniya-differencialnyh-uravnenii-pervogo-poryadka-chislennye.html © muegn.ru