Тор құрудың әдістері

Жұмыс түрі: Дипломдық жұмыс
Тегін: Антиплагиат
Көлемі: 48 бет
Таңдаулыға:

Жоспар

Кіріспе3

1 Бір айнымалы функцияларды интерполяциялау5

1. 1 Бір айнымалы функцияны кубтық сплайнның көмегімен интерполяциялау5

1. 2 Тегістеу арқылы құрама-кубтық интерполяциялау8

1. 3 Тегіс толықтыру11

1. 4 Екі немесе бірнеше айнымалылар функциясын интерполяциялау12

1. 5 Көп айнымалылар функциясының r - тегіс жуықтауы15

2 Есептік торлар құрудың әдістері17

2. 1 Алгебралық әдістер. Дифференциалдық теңдеулерді19

шешуге негізделген әдістер19

2. 2 Бір байланысты аймақ26

2. 3 Көп байланысты аймақ30

2. 4 Көпбайланысты аймақ үшін эквиүлестіру әдісі32

2. 5 Дербес туындының теңдеуін шешу негізінде құрылатын торлар35

2. 6 Интегралдау аймағында тор құрудың әдістері37

2. 7 Бейімді торлар. Тор құрудың әдістері39

3 Торларды құру49

3. 1 Кездейсоқ қисық сызықты координаталарда, ортогональ қисық сызықты тор құру51

Қорытынды60

Кіріспе

Физикалық процесстерді математикалық моделдегенде, шекаралары қисық сызықты болатын физикалық аймақ жиі кездеседі. Мұндай аймақтарда бейімді қисық сызықты тор құрған, санау уақытын қысқартуға, бағдарламалауды автоматтандыруға мүмкіндік береді. Соңғы жылдары қисық сызықты тор құру әдістері және ерікті қисық сызықты координаталарда гидродинамика есептерін сандық шешу жылдам дамыды. Бұл мәселе С. К. Годунова, Г. П. Прокопова [1], Н. Н. Яненко, В. Д. Лисейкина, Н. Т. Данаева [2-3], П. Н. Вабищевича [4], J. F. Thompson, Z. U. A. Warsi, C. W. Mastin [5] бұл жұмыстарда негізінен қисық сызықты және динамикалық бейімді торларды, бір байланысты, екі өлшемді және үш өлшемді аймақтарда құру қарастырылған.

Торда берілген функцияны барлық облыстарын үздіксіз функциялармен толықтыру кезінде интерполяция мәселесі туындап отырады. Бұл жайлы Н. С. Бахвалов, А. А. Самарский, Н. Н. Калиткин, Г. И. Марчук, А. Л. Самарскийдің жұмыстарында қарастырылған. Интерполяция есебі жоба-конструкторлық жұмысын автоматтандыру жүйесінде маңызды, негізгі бөлігі болып табылады. Бұл интерполяциялау мәселесі математикалық әдебиеттерінде классикалық әдістерде толығымен баяндалғандықтан жаңадан мәселе болып табылмайды. Соңғы он жылдықта интерполяция теориясының бағыты сплайн интерполяцияны қолдану болып табылады.

Сплайн интерполяция қисық сызықты торлы функцияны тегіс толықтыруды тұрғызуда ең жақсы амал. Сплайнның тиімділігі оның арнайы экстремальді қасиетімен байланысты. Сплайндық жуықтау барлық ғылым мен техника ортасында кеңірек қолданылып келеді.

Жалпы дискретті нүктелер жиынында берілген шаманы интерполяциялау мәселесі варияциялы-айырымдық сұлбаны және айырымдық есептердің шешімін үзіліссіз көрсетілімін құрумен тығыз байланысты. Біз есептерді басқа тұрғыда қарастырамыз. (a, b) кесіндісінің n+1 нүктесінде (интерполяция түйінінде) функцияның -ң мәні берілген.

2. барлық қалған мәндер үшін

қосалқы функцияны (интерполяцияланған функцияны) жинақтау керек. Басты сұрақ мынада - интерполяцияны қалай таңдау керек және кесіндісінде функцияның жуықтау дәлдігін қалай бағалау керек.

Сандық әдісте ең негізгі есептердің бірі - функцияны интерполяциялау болып табылады. Интерполяция кесіндісінде барлық -тің мәні үшін функциясын қалпына келтіруден тұрады, егер кесіндінің нүктесінде ақырғы санының мәні белгілі болса. Бұл мәндер қандай да бір жасалған тәжірибе бақылауында (өлшемінде) немесе есептеу нәтижесінде табылуы мүмкін.

кесіндісіне тор берілсін делік, оның тоғысында тең болатын функциясының мәні берілген болсын. Бұл жерде жұмыстың өзектілігі - тор түйіндерінде функциясына сәйкес келетін функциясын интерполяциялау, яғни интерполяцияны тұргызуды қажет етеді.

Жоғарыда айтып кеткендей, соңғы жылдары интерполяция жаңа әдістермен байи түсті. Олардың бірі сплайндық интерполяция.

Негізгі мақсаты - қарастырылатын қисық сызықты облыстың шекарасынын интерполяциялауда үзіліссіздігін қамтамасыз ете отырып, интерполяциялаудың бағдарламасын құру.

Жұмыстың міндеттері :

Бір айнымалы функцияларды интерполяциялау.
Есептік тор құру әдістері.
Тор салудың әдістерін қарастыру.
Бір байланысты аймақта тор құру.
Кездейсоқ қисық сызықты координаталарда, ортогональ қисық сызықты тор құру.
Қисық сызықты облыстың аймағын интерполяциялау.

Жұмыстың құрылымында: кіріспе, негізгі бөлім, қорытынды бөлім, әдебиеттер тізімі, формулалар, суреттер қолданылады.

1 Бір айнымалы функцияларды интерполяциялау

1. 1 Бір айнымалы функцияны кубтық сплайнның көмегімен интерполяциялау

нақты осінің кесіндісінде анықталған, түйінінде функциясының мәні берілген тор берілсін делік. Онда құрама-кубтық интерполяция есебі келесі түрдегідей қойылады. [a, b] кесіндісінде келесі талаптарды қанағаттандыратындай функциясын табу керек:

С²класына жатады, яғни екінші ретке дейін өзінің үздіксіз туындыларымен бірге.
Әрбіркесіндідекубтық көпмүшелік түрінде

(1. 1. 1)

болып табылады.

{xn}k=0 тордың түйініде теңдік орындалады.

(1. 1. 2)

g´´(x) шекаралық шартты қанағаттандырады g´´(a) = g´´(b) =0 (1. 1. 3)

Алдағы уақытта анықталған функциясының жалғыз экстремальді қасиетін орнатқан кезде біз таңдап отырған интерполяцияның артықшылығы анықталады.

Интерполяцияланған құрама-кубтық функциясын табуға берілген есептің жалғыз шешімі бар екенін көрсетейік. Ол үшін жоғарыда көрсетілген (1) -(4) талаптарды қолданамыз.

функцияның екінші ретті туындысы әрбір тордың кесіндісінде үздіксіз және сызықты болғандықтан, біз болғанда

(1. 1. 4)

деп жазамыз. Мұнда .

Теңдіктің екі бөлігін де екі рет интегралдасақ

(1. 1. 5)

аламыз, мұнда және интегралдау тұрақтылары. Олар g(xi-1) =fi-1 шартынан шығады. (1. 1. 5) x =xi және x =xi-1 қоя отырып,

аламыз.

Қорытындылай келе

(1. 1. 6)

(1. 1. 7)

(1. 1. 7) формуласынан туындының x1, x2, …, xn-1 нүктелерінде бір жақты шектеулерді табамыз.

g´´(x) және g´(x) функциясы үздіксіз шартына сәйкес нүктесіндегі g´(x) үздіксіз шартынан n-1 теңдігін аламыз.

(1. 1. 8)

(1. 1. 3) шартынан осы теңдікті толықтыра отырып, белгісіз табу үшін сызықты алгебралық жүйені аламыз.

Am=Hf (1. 1. 9)

А квадратты матрица

(1. 1. 10)

түрде болады. m және f векторлары және H тікбұрышты матрицасы былай болады:

(1. 1. 11)

А матрицасы қатаң диагональдық басымдылықпен симметриялы. Гершгоринның өзіндік мәнін таратпаушылық жайлы теоремасы бойынша айқындалған және де ерекше емес. Сонымен, m1, m2, …, mn-1 коэффициенттері (1. 1. 9) жүйесінен бір мәнінде анықталады. Сондықтан, g(x) сплайн функциясы да (1. 1. 1- 1. 1. 6) формуласымен бір мәнде қалпына келтіріледі және g(x) құрама-кубтық функцияны табудағы есеп жалғыз шешімнен тұрады.

Құрама-кубтық функциясы сплайн-интерполяциясының жоғары тиімділігін қамтамасыз ететіндіктен өте маңызды қасиеттерді иеленеді. Атап айтқанда, екінші ретті туындының квадратты қосындылаулары бар функциядан тұратын класын [a, b] кесіндісінде қарастырайық. Есепті класында функционалды минимизациялайтын ізделінді интерполяцияланған функцияға қоямыз:

, u(xk) =fk, k=0, 1, …, n (1. 1. 12)

(1. 1. 13)

Мұндай функционалдың минимумы жәнегі біз салған g(x) құрама-кубтық сплайн-функциясына жететіндігі анықталады. Шынымен де мына шаманы қарастырайық:

(1. 1. 14) g мен функцияларының қасиеттерін қолдана отырып және бөлшек бойынша интегралдай отырып

Бірақ = Ck=const [xk-1, xk] кесіндісінде болғандықтан

Осыдан және (1. 1. 14) формуласынан кез келген функциясы үшін

(1. 1. 12), (1. 1. 13) формулаларына негізделе отырып, құрама-кубтық сплайн-функциясына: тор түйінінен берілген мәнді қабылдайтын, классынан және (1. 1. 13) функционалын минимилизациялайтын функция деп эквиваленттік анықтама беруге болады.

Кубтық сплайн-функциясы жақсы жуықтау қасиетін иеленеді. Егер иетерполяцияланатын функция класына жататын болса, онда функция қателігі үшін

(1. 1. 16)

теңсіздігі болғаны жөн, мұндағы - торына тәуелсіз теріс емес тұрақты.

1. 2 Тегістеу арқылы құрама-кубтық интерполяциялау

Бұл бөлімде біз тағы да торда анықталған тегіс толықтыру жайлы есебін қарастырамыз. Бірақ та функцияның мәні тор түйінінде біршама қателіктерге ие болады. Бұл кезде түйіндіде, дәлдікте берілген мәнмен сәйкес келетін интерполяциянған функцияны салудың қажеті жоқ. Оның үстінде интерполяцияланғаннан гөрі берілген мәнге жақынырақ өтетін функцияны салу керек болады. Ал мұндай функциялар енді интерполяциялық емес, тегістелген деп аталады.

Ізделінді тегістелген функциясы классында

(1. 2. 1)

функционалын минимизациялауын талап етеміз, мұнда - кейбір оң сандар. функционалында берілген қисық мәннің жанынан және «бүгу» функциясының минимал шартынан өтетін интерполяциялық шарттар араласқан. зілдеме коэффициенті көп болған сайын, функционалға интерполяциялық шарттарды қосады, берілген мәнге жақын тегістеу функциясы өтеді.

(1. 2. 1) вариациялық есептердің шешімі кубтық сплайн болып табылатынын қарастырайық, яғни алдыңғы бөлімшедегі талаптарды қанағаттандыратын функция. - есептің шешімі болсын. болатын сплайнын тұрғызамыз. (1. 2. 1) екінші қосылғыш пен функциясы үшін бірдей болғандықтан

(1. 2. 2)

Бірақ, алдыңғы бөлімшеде көрсетілгендей, интерполяциялау кезінде минимум өрнегін беретін жалғыз функция. Сондықтан .

Сонымен функционалының минимумын кубтық сплайн классында табу жеткілікті. түйінінде қабылдайтын кубтық сплайн бір мәнде оның мәнінде жиынмен анықталғандықтан, минимизациясы айнымалылардан функция минимумын табуға апарады.

- құрама-кубтық функция екенін білеміз, сонымен ,

, ,

(1. 2. 3)

Болғандықтан

(1. 2. 4)

(1. 2. 4) формулада интегралды туындай келе

(1. 2. 5)

аламыз.

Мұнда - белгілі матрица. Сонымен

(1. 2. 6)

1. 1. 9 нәтижесінде арқылы сызықтық көрсетіледі, сондықтан - да - ден анықталатын оң форм болуы мүжетті шарты болып табылатын оның экстремумы болып тек қана минимум болуы мүмкін. Бірақ матрицасы - ге тәуелді емес.

Сондықтан (1. 1. 9) байланысты

(1. 2. 7)

Мұнда (1. 1. 11) формуласымен анықталады. Осыдан минимум шарты векторлық формада мынандай түрде (1. 2. 7) шығады, мұнда , ал - диагональді матрица:

(1. 2. 8)

1. 2. 7-нің сол жағын көбейте отырып, , немесе

(1. 1. 9) ескере отырып, түпкілікті

(1. 2. 9)

аламыз.

(1. 2. 9) жүйенің матрицасы бес диагональды, симметриялы және дұрыс анықталған. (1. 2. 9) жүйесін мысалы, Гаусс әдісімен шығаруға болады. Осыдан кейін m векторы анықталғандықтан, 1. 2. 7 формуласынан жалғасатын формуласымен тегістеу сплайнның торлы мәнінің векторын іздеу қажет. Одан 1. 1. 6 формуласы бойынша сплайны қалыптасады.

Кубтық сплайнмен тегістеу - тегістеудің жалғыз ғана әдісі емес.

1. 3 Тегіс толықтыру

Енді торлы функцияның тегіс толықтырудың басқа әдісін айтамыз, оның сплайн-функция әдістер теориясынан бірнеше айырмашылығы бар, сонымен алгоритмдік көзқараста өте тиімді. тұрақты тегістеу классының интерполяциялық функциясын құру процессін өте қысқаша сипаттайық. Түйінде функцияларының мәні белгілі - кейбір бекітілген тор болсын. жетерліктей үлкен деп болжайық. бүтін санын белгілейік. Алдымен, тор кесінді де тегістеу классының интерполяциялық функциясын құрайық. Тор түйінінде берілген мәнімен сәйкес дәрежесінен жоғары емес келетін Лагранж көпмүшелігін арқылы белгілейміз, ал арқылы түйінде берілген мәндерді қабылдайтын дәрежесінен жоғары емес Лагранж көпмүшелігін белгілейміз. Ары қарай мына шартты қанағаттандыратын

(1. 3. 1)

Қанағаттандыратын, 2p+1 дәрежесінен жоғары емес көпмүшелігін тұрғызамыз. Мұндай көпмүшелігі бар және шартта бірмәнде анықталғаны айқын.

кесіндіде интерполяциялық функция тең болсын. Сурет ретінде p=1 жағдайын қарастырумен шектелеміз. Онда кесіндіде 3-ші дәрежелі көпмүшелігін:

түрде жазылады. Олардың коэффициенттері

(1. 3. 2)

формуласы бойынша есептеледі.

Егер бастапқы үшін санауда нүктесін қабылдасақ, онда кесіндіде сәйкес интерполяциялық функциясы салынады. Мұндайд салу процессін интервалына дейін жалғастыруға болатыны айқын. Салынған интерполяция тегістеу классының құрама-полиномиалдық функциясында бар; тор интервалында полиномның дәрежесі 2p+1-ге тең.

Егер p=1 деп қойсақ (кубтық сплайн тегістеуіне сәйкес келеді), онда бес дәрежелі көпмүшелігімен жұмыс жасаймыз. Көпмүшелік дәрежесі 2-ге жоғары, бірақ тұрғызу үшін сызықтық алгебралық жүйені шешудің қажеті жоқ.

Тегіс толықтыру жақсы жуықтама қасиеттерін иеленеді. Әсіресе, классының толықтырумен көп айнымалылар функциясын интерполяциялау кезінде, болғанда, функцияның өзі үшін және оның туындылары үшін

(1. 3. 3)

теңсіздік болады, мұнда мультииндекс

-дан тәуелді және h тор қадамы мен f функциясынан тәуелді емес.

Тегіс толықтыру тағы да бір қасиетті иеленеді. Егер бір түйінді бірге және қалғандарында нольге тең функциясын, яғни базистік функциясын тұрғызсақ, онда кез-келген p үшін финитті болады, сонымен тасушының «радиусы» p-ны асырмайды. Егер де базистік құрама-кубтық сплайн-функциясын салатын болсақ, онда ол мүлдем финитті емес. Бұл жағдай ақырғы элемент әдісімен шешілетін вариациялық есептерде тегіс толықтыруды өте ыңғайлы етеді.

1. 4 Екі немесе бірнеше айнымалылар функциясын интерполяциялау

Құрама - бикубтық функция көмегімен екі өлшемді интерполяциялау проблемасы көптеген автордық зерттеу пәніне айналды. Біз келесі модельді есебін қысқаша қарастырамыз.

- жазықтығының кейбір тікбұрышы. ға тор құрайық.

нүктесінде берілген, осындай құрама - бикубтық функциясының интерполяциялау есебінің болжамы келесі шарттарды қанағаттандыратын функциясын тұрғызумен аяқталады:

1) (1. 4. 1)

2) торының әрбір ұяшығында бикубтық көпмүшелік түрі болып

(1. 4. 2)

3) торында берілген мәнді қабылдайды.

(1. 4. 3)

4) функциясы (1. 4. 4) шартын қанағаттандырады (мұнда v-D облыстың Г шекараға ішкі нормаль) .

Мұндай функцияны қалай құруға болады? Былайда, құру бір өлшемдіден еш айырмашылығы жоқ. 1. 1. 6 қарапайым формула бойынша кез келген нүктеде бірөлшемді сплайн-функцияны есептеу үшін функцияның мәнін және оның түйінді нүктесінде екінші ретті туындысын білу керек екендігін еске түсірейік.

Екінші ретті туындыны табу үшін бір рет үшөлшемді матрицамен сызықтық алгебралық жүйені шешу керек.

Нақты формула бойынша екіөлшемді жағдайда кез келген нүктеде функцияны есептеуге болатындай қандай алғашқы есептеулер жасауға болады?

Алдымен тор сызығындағы кубтық сплайн интерполяциямен бірөлшемді есептерді қарастырайық. Ол үшін 1. 1. 9 типінің m+1 сызықтық алгебралық жүйесін шешеміз. Нәтижесінде тор түйіндегі функция мәнін табамыз. Одан кейін сәйкес сызығындағы бірөлшемді сплайн-интерполяция есебін шешеміз және торында функция мәнін табамыз. Енді - ң кейбір нүктелерде мәнін есептеу керек деп болжайық. , болсын.

(1. 1. 6) сәйкес , нүктелерінде мәнін мына формула бойынша:

, (1. 4. 5 )

(1. 4. 6)

табуға болады. Мұнда ары қарай , , ,

Егер және мәндері белгілі болса, онда 1. 1. 6 типінің формуласы бойынша мәнін табуға болады. бойынша функциясы құрама-кубтық болып табылатынын белгілейміз. Тор мәндері белгілі, функциясы үшін сызықта m+1 бірөлшемді есебін шешеміз. Нәтижесінде торында функциясын табамыз.

белгілесек, онда

(1. 4. 7)

(1. 4. 8)

Тағы да 1. 6 формуласын пайдаланып мәндерін табамыз:

(1. 4. 9)

Жоғарыда айтылғандардың бәрін қорытындылай келе, бикубтық сплайнды есептеу үшін қажетті есептеулердің көлемін бағалауға болады. керекті нүктеде есептемес бұрын бір рет (1. 1. 9) типінің сызықтық алгебралық жүйесін (n+1) +(m+1) +(m+1) =2m+n+3 шешу керек. Ары қарай облысты бір нүктесінде функциясын есептеуі үшін (1. 1. 6) типінің, әсіресе: (1. 4. 5), (1. 4. 6), (1. 4. 7), (1. 4. 8), (1. 4. 9) формулалары бойынша 5 рет есептеулерді орындау керек. Осы формулалардың тізбегін тордың әрбір ұяшығында полиномиальді айқын өрнегін алу үшін қолдануға болады. Бірақ та ЭЕМ оперативті жадта көпмүшеліктің массивтердің коэффициенттерін сақтау үшін 16пт жуық ұяшықты талап етеді. Бұл алгоритмде айнымалыларды орындарымен ауыстыруға болады және сонда 2m+n+3 бірөлшемді есептерді шешу керек болады, бірақ соңғы нәтиже өзгермейді. Сипатталған алгоритм параллелипед типінің көпөлшемді облыстарында оңай жалпыланған болуы да мүмкін.

1. 5 Көп айнымалылар функциясының r - тегіс жуықтауы

Негізгі бөлімде ретсіз торда қатемен берілген көп айнымалылар функциясын толықтырудың есебіне мүмкін жолдардың ішінен біреуі қарастырылады. Ұсынылып отырған әдіс негізінде талдауда белгілі бірлікті бөліктеу түсінігі жатыр. классының ( - ашық көпмүшелік) функциясының тізбегі көпмүшелігінде бірлік бөліктеуді жасайды деп айтады, егер және барлық үшін.

Әдетте берілген ашық көпмүшеден тыс функциясы нольге бағытталуын талап етеді. ( ) жиынтығы көпмүшелігін бүркейтін болар. ( ) бірлік бөліктеуді осы бүркеуге бағынышты деп атайды. Қолданбалы мақсаттар үшін стандартты көпмүшелікпен: m-өлшемді интервалмен немесе шармен төңіректе ақырғы бүркеу сай. Осы көпмүшеліктің ортасы мен диаметрі функциясының жуықтау қасиеттерімен сәйкес таңдалыну керек. Алдағы уақытта m-өлшемді кубымен бүркеуді қарастырамыз, ал оларға тәуелді бірлік бөліктеуді келесі түрде құрамыз.

кубын ортасы арқылы белгілейміз, ал арқылы - оның қабырғасының ұзындығын белгілейміз. классының егер , егер болатын стандартты функциямен берілеміз. Мысалы ( болғанда)

немесе

қоюға болады. Енді әрбір i индекс мәні үшін

қоямыз. функциясы класының функциялары болып табылады және бүркеуге тәуелді -да бірлікті бөліктеуді жасайды. Енді әрбір функциясынан -ке - да «жақын» класының функциясының тізбегі функциясы үшін белгілі деп болжаймыз. Әрбір үшін бірқалыпты нормада болсын. Онда теңдікпен

(1. 5. 1)

анықталатын функция - да «жақын» болады. Шынымен тепе-теңдікке байланысты әрбір үшін

аламыз, яғни .

Енді функциясын толықтыру үшін бірлікті бөліктеуді қолданамыз; егер кейбір тор торабында берілген болса, онда үшін кубының, оның төңіректік жуықтауы ортасына жақын болатын кейбір жиынтық бойынша торабын құру керек, содан кейін осы жуықтауды (1. 5. 1) формуласы бойынша «жалғау». Алдағы уақытта біз кубының орташасымен түйіні сәйкес келеді деп есептейміз, ал төңіректік жуықтауы ге тең тұрақты функция болып табылады.

Жалғыз жеке жағдайда жуықтау қателігіне баға береміз: тор- нүктесіндегі түйінімен бірқалыпты және шексіз; интервалдардың орташалары тор түйінімен сәйкес келеді, ал , мұнда - кейбір натурал сан. бірлікті бөліктеу класына тиесілі болсын, функциясы жұп. Мұндай жағдайда қалдық өзін тақ, h периодымен периодты функция түрінде береді. Сондықтан да

2) жағдайын қарастыру жеткілікті.

1-ші жағдайды қарастырайық. және функциясы жұп болғандықтан

Осыдан Тейлор формуласы бойынша

аламыз. Осы теңдіктің оң жақ бөліктегі екінші қосынды нольге тең екені айқын, сондықтан

(1. 5. 2)

Ары қарай қосындыға

Абельдің түрлендіруін қолданып

аламыз.

Сондықтан

Мұнда

(1. 5. 3)

онда (1. 5. 4)

(1. 5. 2) (1. 5. 3) формуласынан бағалауын

(1. 5. 5)

аламыз. Мұнда С - 3 ретті функциясының модулінің максималды мәні. Тұрақты C N-нан тәуелді және N бойымен кемиді.

2 жағдайда да осындай нәтижеге келеміз. Сәйкес бағалауларды сызықтық функция үшін және бірнеше айнымалылар жағдайында алуға болады [10] .

2 Есептік торлар құрудың әдістері

Дифференциалды теңдеу жүйелерін дербес туындымен шешуде жақсы құрылған есептік торларды қолданып біршама жеңілдетуге болады. Бірақ та қанағаттандыратындай нәтиже бермеуі мүмкін. Кейбір қосымшаларда есептіқ торлардың түйіндерінің орналасуының бірдей емес таңдауы тұрақсыздыққа немесе жинақтылыққа әкелуі мүмкін. Теңдеулерді дербес туындымен сандық шешу кезінде есептік торларды тұрғызу мәселесі туындайды.

Қарапайым теңдеуді шешіп көрейік,

кейбір аймақтарында бастапқы және ақырғы шарттарымен сәйкес. Есептеу әдетте есептеу кеңістікте жүргізіледі. Бастапқы дифференциалдық теңдеуінің дербес туындысы x, y физикалық координатасынан есептік кеңістік координатасына өтуі күрделі дифференциалдық функциялардың ережелерінің көмегімен түрленеді:

(2. 1)

Онда ол мынандай түрде болады:

(2. 2)
Осы теңдеу есептеу кеңістігінде бірқалыпты торда шешіледі. Міндетті түрде физикалық кеңістік пен есептеу кеңістіктегі координаталарының арасына байланыс орнату керек. Бұл байланысты түрленудің метрикалық коэффициенттері береді ( мүшелері қарастырылған теңдеуде дербес туындысымен) .

Сурет 1. Физикалық кеңістіктің есептеуіш кеңістігіндегі кескіні: а) физикалық кеңістік ; b) есептеуіш кеңістік

Есептік торларды құру есебі тор түйіндісін физикалық кеңістіктен (D) есептеуіш кеңістікке (CD) ауыстыратын кескінді табудан тұрады. Бұл кескін кейбір талаптарды қанағаттандыру керек. Солардың ішінен кейбіреуін атап кетуге болады:

1. Кескін бірмәнді болу керек.

2. Тор сызығы тегіс болу керек, ол туындынының үзіліссіздігін қамтамасыз етеді.

3. Тор (D) аймағының үлкен сандық қателіктер туындайтын бөлігіне қалың болу керек.

Бірөлшемді есептеуіш тор құру оңайлау. Бұл торларды тұрғызуда көптеген функциялар немесе әдістер бар. Сондықтан көбінесе екі өлшемді тор тұрғызуда сұрақтар туындайды. Ал үшөлшемді тор тұрғызу өте күрделі есеп болып табылады және қанағаттандыратындай нәтиже беретін әдістер сондай көп емес.

Есептеуіш тор құру әдістерін дөрекі түрде 3 классқа бөлуге болады:

Комплекстік айнымалы функциялар теориясының әдістері.
Алгебралық әдістер.
Дифференциалдық теңдеулерді шешуге негізделген әдістер.

Біз алгебралық әдістер мен дифференциалдық теңдеулерді шешуге негізделген әдістерін қарастырамыз. Жалпы алгебралық әдістер мен дифференциалдық теңдеулерді шешуге негізделген әдістерді күрделі үшөлшемді есептерге қолданылады. Бұл екі әдіс есептеуіш торларды құру әдістерінің ішіндегі тиімді және кең таралған болып табылады.

2. 1 Алгебралық әдістер. Дифференциалдық теңдеулерді

шешуге негізделген әдістер

Сурет 2.

Тор тұрғызу кезінде есептеуіш кеңістік тік бұрышты болу міндетті емес. 2-суретте көрсетілген кескінде тор тұрғызу керек болсын. Оның қабырғасы функциямен беріледі:

(2. 1. 1)

Бұл мысалда x координатасы бойынша тұрақты қадамды таңдай отырып және соплоның осьі мен қабырғасы арасындағы әрбір кесіндіні бірдей бөлікке бөле отырып тор тұрғызу оңай. Бұл процедура келесі тәуелдімен сипатталады:

(2. 1. 2)

мұнда - сопло қабырғасының теңдеуі. Осыдан берілген мен мәндері бойынша x пен y мәндерін табу оңай. Физикалық кеңістікте тұрғызылған тор 2. 3- ші суретте кескінделген.

Түрлендірудің метрикалық коэффициенттерін есептеуде мұқият болу керек. (2. 1. 2) теңдеуі бойынша есептелетін, дербес, пен туындылар келесі түрде болады: