Әдістемелік жинақ



бет40/40
Дата06.01.2022
өлшемі2,2 Mb.
#16295
1   ...   32   33   34   35   36   37   38   39   40
xi

yi

Xi-x

Li-1,i

Li-2,i-1,i

Li-3,i-2,i-1,i



X0

Y0

X0-x













X1

Y1

X1-x

L01(x)










X2

Y2

X2-x

L12(x)

L012(x)







X3

Y3

X3-x

L23(x)

L123(x)

L0123(x)



X4

Y4

X4-x

L34(x)

L234(x)

L1234(x)



Эйткен схемасын есептеуді көршілес L0123…n(x), L0123…n,n+1(x) мәндері берілген дәлдік маңайында бір бірімен беттессе тоқтатуға болады.

Xi нүктелерінде yi мәндерін қабылдайтын n-ші дәрежелі интерполяциялық көпмүшелік келсі түрде де жазылады:



. (4.7)

1-Мысал: Кестемен берілген функция үшін Лагранж көпмүшелігін құру.

I

0

1

2


(4.8)
3

xi

0

0.1

0.3

0.5

yi

-0.5

0

0.2

1

Шешімі: (4.4)-формула бойынша n=3, i=0,1,2,3 болғандағы өрнекті анықтаймыз:

L13(x) мүшесін есептемейміз, себебі y1=0. Бәрін бір біріне қосамыз да көпмүшеліктің соңғы түрін аламыз:






(4.9)
2-мысал: Кесте мен берілген функцияның x=0.45 нүктесіндегі мәнін анықтау керек.

X

0.05

0.15

0.20

0.25

0.35

0.40

0.50

0.55

y

0.9512

0.8607

0.8187

0.7788

0.7047

0.6703

0.6065

0.5769

Шешімі:

92

Есептеуді жеңілдету үшін x=0.05t деп алайық. X-тердің мәні белгілі болғанда t-лардың мәндерін тауып алуға болады, олар: 1, 3, 4, 5, 7, 8, 10, 11. Және x=0.45 болғандағы t=9 болады. Есептеу қадамдары 2-кесте да келтірілген.



2-кесте. (4.9)-есептің есептелу қадамдары.

i


ti-tj

(i<>j)


Di


yi






0

8

-2

-3

-4

-6

-7

-8

-10

-725 760

0.9512

-0.0131*10-4

1

2

6

-1

-2

-4

-5

-7

-8

26 880

0.8607

0.3202*10-4

2

3

1

5

-1

-3

-4

-6

-7

-7 560

0.8187

-1.0829*10-4

3

4

2

1

4

-2

-3

-5

-6

5 760

0.7788

1.3520*10-4

4

6

4

3

2

2

-1

-3

-4

-3 456

0.7047

-2.0390*10-4

5

7

5

4

3

3

1

-2

-3

2 520

0.6703

2.6599*10-4

6

9

7

6

5

5

2

-1

-1

11 340

0.6065

0.5348*10-4

7

10

8

7

6

6

3

1

-2

-80 640

0.5769

-0.0715*10-4






Сонымен y(0.45)= *1.6604*10-4=0.6376.



Егер есепте керісінше функция мәні белгілі болып сол мәнге сәйкес абсцисса мәнін табу керек болса, ондай есепті кері интерполяциялау деп атайды. Кері интерполяциялау формуласы:


(4.10)


Ньютоннның интерполяциялық формулалары

Егер интерполяциялық түйіндердің бір бірінен ара қашықтығы тұрақты болса, практикада Ньютонның интерполяциялық формулалары қолданылады. Бұл формулалар екіге бөлінеді:



  1. Алдыға қарай интерполяциялау

  2. Кері интерполяциялау

Егер берілген х нүктесінің мәні кестенің бас жағында жатса, 1-формуласы қолданылады:


(4.11)
.

Мұндағы

Егер берілген х нүктесінің мәні кестенің соңғы жағянда жатса, 2-формула қолданылады:


(4.12)

93

Формулалардағы , , т.с. сияқтылар шектік айырымдар деп аталады және 3-кестені толтыру арқылы анықталады. Кестедегі мысал үшін 6 интерполяциялық түйін және шектік айырымдардың 4-ші дәрежесіне дейінгі мәндер қарастырылған. 1-формула үшін кестенің бірінші жолындағы мәндер, 2-формула үшін кестенің соңғы жолындағы мәндер қолданылады.

3-кесте . Шектік айырымдар кестесі.



x

Y









X0

Y0









X1

Y1









X2

Y2










X3

Y3











X4

Y4












X5

Y5













Егер интерполяциялық түйіндер саны 1 немесе 2-ге тең болса сызықты интерполяциялық формуланы қолдануға болады: .

Қателіктерін бағалау:



1-формула үшін мына формула қолданылады:

,

немесе

2-формула үшін мына формула қолданылады:



,

Мысал: y=lg(x) функциясының мәндері кесте да берілген, lg1001 мәнін табу керек.

4-кесте . y=lg(x) функциясының мәндері және шектік айырымдары кестесі



x

Y







1000

3.0000000

43 214

-426

8

1010

3.0043214

42 788

-418

9

1020

3.0086002

42 370

-409

8

1030

3.0128372

41 961

-401




1040

3.0170333

41 560







1050

3.0211893









3-ретті шектік айырымдар тұрақты бола бастағандықтан кесте ны толтыруды тоқтатамыз. Формулада n=3 деп аламыз. Q=0,1. x=1001. Ньютонның бірінші формуласын қолданамыз, себебі х-тің мәні кестенің бас жағында жатыр, сонда lg1001=3.00043417+0.5*10-9 болатынын қалдық мүшенің формуласын қолдану арқылы анықтаймыз.


Бақылау сұрақтары:

  1. Лагранждың интерполяциялық көпмүшелігі?

  2. Эйткен схемасы?

  3. Ньютоннның интерполяциялық формулалары?

  4. Шектік айырымдар кестесі қалай есептеледі?

94

Тақырыбы: Гаусс квадратурасын қолданып, Рунге принципінің негізінде интегралды есептеу

Мақсаты: Функцияны интерполяциялау есептерін шығаруды үйрену

Тапсырма:

E=10-5 дєлдікпен тµмендегі интегралдарды сєйкес єдісті тањдау арќылы есептеңіздер.

1. ;

2.

3.

4.

5.

6.

7.

8.

9.

10.


11.

12.

13.

14.

15. , аралыќты 2-ге бµлу арќылы есептеу.
95

16.

17.

18.

19.

20.



Тапсырманы орындауға әдістемелік нұсқаулар:

Канторовичтің интеграл астындағы функцияның ерекшелігін айшыќтау әдісі.

Б±л єдістіњ негізгі идеясы f(x) функциясынан єлдебір g(x) функциясын бµліп алады да f(x)-g(x) айырымынан интеграл алады, м±ндаѓы f(x)-элементарлы интегралданатын болуы керек, ал g(x) – сандыќ єдістердіњ біреуімен интегралданатын болуы керек:



(5.7)

Б±л єдіс интеграл астындаѓы функция келесі т‰рде берілген жаѓдайда ќолданылады:



М±ндаѓы функциясы [a,b] аралыѓында ‰зіліссіз дифференциалданатын функция болсын. Онда б±л функцияны Тейлор ќатарына жіктеуге болады. Тейлор ќатарына жіктелген функция т‰рін деп белгілесек:



болады. Сонда айырманы деп ќарастыруѓа болады.

Мысалы: ;



Интеграл астындаѓы функияны келесі т‰рде жазып алуѓа болады:

. Онда болады. Б±л функцияны Тейлор ќатарына жіктеуге болады:

. Онда: , болады.



Дәрежелік ќатарлар кµмегімен интегралдау єдісі.

Интеграл астындаѓы функция [a,b] аралыѓы жатќан (-R,+R) интервалында жинаќталатын дєрежелік ќатарѓа жіктелетін болсын:

96

(5.8)

Дєрежелік ќатарды м‰шелеп интегралдауѓа болады, сонда:



(5.9)

Егер (5.9)-ќатар жинаќты болса, онда (5.10)

Єдіс ќателігі ќатар ќалдыѓы мен дµњгелектеу ќателігінен т±рады.

Егер ќатардыњ тањбасы ауыспалы жєне абсолютті шамасы бойынша монотонды кемімелі болса, онда ќатар ќалдыѓы тасталатын (отбрасываемый) ќатар м‰шелерініњ ењ біріншісініњ абсолют ќателігінен аспайды.

Басќа жаѓдайларда ќатар ќалдыѓын баѓалау ‰шін ќатарлары жењіл баѓаланатын сандыќ ќатарлармен мажорлайды.

Мысал:

Интеграл астындаѓы функция дєрежелік ќатарѓа жіктеледі жєне кез келген х ‰шін жинаќталады:







Ќатардыњ ‰шінші м‰шесінен бастап ќалѓан м‰шелерін тастап кетуге болады, себебі олар нµлге µте жуыќ.



Гаусстыњ квадратуралыќ формулалары.

Б±л єдісте х айнымалысын µрнекпен алмастырады: . Сонда Гаусс формуласы былай жазылады:



(5.11)

м±ндаѓы , ал ,

Б±л єдістіњ ыњѓайсыздыѓы, абсциссадаѓы мєндер мен коэффициенттер иррационал сандар. Біраќ соѓан ќарамай интегралдау т‰йіндерініњ саны аз болса да, дєлдігі жоѓары.

Абсциссадаѓы мєндер мен коэффиценттер м‰мкін, жоѓары дєрежелі кµпм‰шеліктердіњ барлыѓы ‰шін (5.11)-формула д±рыс болатындай етіп тањдалынып алуы керек. сандары N=2n-1 болѓанда бірмєнді болатыны дєлелденген.

мєндерін аныќтау ‰шін Лежандр кµпм‰шелігі ќолданылады:

97

Ал -лерді табу ‰шін тµмендегі ж‰йе шешіледі:

(5.11.1)

Б±л ж‰йені шешу барысында коэффициенттерді тµмендегі 1-кесте кµмегімен де аныќтауѓа болады.

1-кесте . (5.11.1)-жүйені шешу барысында анықталған коэффициенттер.


n

I

ti

Ai

1

1

0

2

2

1;2



1

3

1;3

2


0






4

1;4

2;3




0.34785484

0.65214516



5

1;5

2;4


3



0.23692688

0.47862868

0.56888889



Мысалы:

N=5 болѓанда интегралын Гаусс формуласымен есептеу керек болсын. Айнымалыны ауыстырамыз: . Сонда интеграл мына т‰рге келеді:



; (5.12)

;

x=0 болѓанда t=-1;

x=1 болѓанда t=1;

Интеграл астындаѓы функцияныњ мєндер кесте сын ќ±рамыз:



98

Алынѓан мєндерді (4)-ќою арќылы интегралды есептейміз:



.

Интегралдау ќадамдарын тањдау.

Берілген интегралды сандыќ єдістердіњ ќайсысымен болса да шешу уаќытында берілген дєлдікті ќамтамасыз ететін ќадам тањдау керек. Кейде интегралдау аралыѓын бірнеше бµлікке бµлу барысында дµњгелектеу, есептеу ќателіктері µсуі м‰мкін. М±ндай ыњѓайсыздыќќа ±шырамас ‰шін интегралдау ќадамын д±рыс тањдау керек. Практикада интегралдау ќадамын 2 тєсілмен тањдайды:



  1. Ќалдыќ м‰шені баѓалау арќылы

  2. Екілік есептеу арќылы.

1-тєсілде берілген интегралды шешуге тиімді бір сандыќ єдісті тањдап алып, сол єдістіњ ќалдыќ м‰шесініњ формуласын

(5.13)

баѓалау арќылы h –ты аныќтайды.



Мысалы:

Трапеция әдісін таңдайық. Ол әдістің қалдық мүшесінің формуласы:



.

;

Табылѓан мєндерді формуласына ќойып (5.13)-баѓалауды баќылаймыз, сонда екендігі табылады. Енді ќадамды табуѓа болады: . Трапеция формуласына ќойсаќ: .

2-тәсілде берілген интеграл h ќадаммен аралыќты n рет бµледі жєне ќадаммен аралыќты 2n рет бµледі де екі рет есептеледі. алынѓан интегралдарды сєйкесінше жєне деп белгілесек, шарты орындалса, онда деп есептеуге болады, м±ндаѓы I – интегралдыњ дєл мєні. Егер б±л шарт орындалмаса, онда ќадамды таѓы да 2-ге кішірейтеді.

Интегралдау ќадамы берілмеген есептерде, алѓашќы ќадамды санына жуыќ сан ретінде алуѓа болады, м±ндаѓы m=2 трапеция формуласы ‰шін, m=4 – Симпсон формуласы ‰шін. Б±л тєсіл есепті ЭЕМ кµмегімен шешу кезінде ќадамды автоматты т‰рде компьютер тањдайтындай жаѓдай тудырады жєне бірмезгілде есептеу ќадамдары да баќыланып отырады.
Бақылау сұрақтары:


  1. Канторовичтің интеграл астындағы функцияның ерекшелігін айшыќтау әдісі?

  2. Дәрежелік ќатарлар кµмегімен интегралдау єдісі?

  3. Гаусстыњ квадратуралыќ формулалары?

  4. Интегралдау ќадамдарын тањдау қалай жүзеге асады?

99

Тақырыбы: Монте-Карло әдісімен еселі интегралды есептеу

Мақсаты: Интегралдап жуықтап есептеу есептерін шығаруды үйрену

Тапсырма:

E=10-5 дєлдікпен тµмендегі интегралдарды сєйкес єдісті тањдау арќылы есептеңіздер.



1. ;

2.

3.

4.

5.

6.

7.

8.

9.

10.

11.

12.

13.

14.

15. , аралыќты 2-ге бөлу арқылы есептеу.

100


16.

17.

18.

19.

20.


Тапсырманы орындауға әдістемелік нұсқаулар:

Сандық интегралдау есебі

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

Сандық интегралдауды сандық квадратура деп те атайды. Ал қолданылатын формулалар квадратуралық формулалар деп аталады.

Сандық интегралдау да дәл және жуықтау болып екіге бөлінеді.

Егер абсцисса өсі бойынан алынатын нүктелер бірқалыпты орналасатын болса, онда Ньютон – Котестің дәл квадратуралық формулалары қолданылады, басқа жағдайда жуықтау – Гаусс формулалары қолданылады.

Сандық интегралдаудың негізгі идеясы - интеграл астындағы функцияны [a,b] аралығында интерполяциялық полиномға жіктеу және полиномның әр мүшесін интегралдау арқылы есептеу процесін жеңілдету.

Интегралдың қателігін төмендету үшін интеграл астындағы функция анықталған [a,b] аралығы h қадаммен бірнеше аралыққа бөлу керек: xi+1-xi=h, i=1,2,…,n-1. Қадам тұрақты болған жағдайды қарастырайық.

(5.1)

түрдегі интеграл берілсін. Дәл әдістерге Ньютон-Котес квадратуралық формулалары жататыны жоғарыда айтылған.



Трапеция әдісі

Егер n=1 болса квадратуралық формула трапеция әдісі деп аталады. Әдіс бойынша; интегралдық қисық пен ох өсі аралығындағы фигура ауданын табу үшін сол фигураны трапециямен толықтырып, ауданын табуға болады:



(5.2)

Қателікті азайту үшін аралықты бірнеше бөлікке бөліп әр трапецияның ауданын тауып барлығының қосындысы берілген интегралдың мәні деуге болады:



(5.3)

Мұндағы



, (5.4)

.

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



1-мысал:



(5.3)- формулаға қоямыз:



101











Симпсон әдісі.

N=2 болса Ньютон – Котес формуласы Симпсон єдісін аныќтайды. [a,b] аралыѓын екі симметриялы бµлікке бµледі: н‰ктелері болса, аралыќ ж±п болады, есептеу формуласы:



(5.5)

Онда есептеу ќателігі 16-ѓа азаяды. Ал бµлу аралыѓы таќ болса, онда [a,b] аралыќтыњ алѓашќы ‰ш бµлігінен ‰шінші дєрежелі парабола ж‰ргізнміз, б±л жаѓдайда Симпсоннныњ ‰штен сегіздік формуласы ќолданылады:





Мысал:

a=0; h=0.25; b=1;



102




.

Тіктөртбұрыштар әдісі.

[a,b] аралығынан х0 бір түйін алатын болсақ, чғни f(x)=const болады, онда қарастырып отырған аралықта деуге болады. Х0 нүктесін аралықтың тура ортаңғы нүктесі деп алсақ формуласы шығады, оны тіктөртбұрыштар формуласы дейді, әдістің қателігін азайту мақсатында аралықты бірнеше бөлікке бөліп, әр аралықты тіктөртбұрышпен толтырып, ауданын тауып, барлық аудандарды бір біріне қосады:



(5.6)





Мысалы:




Табылған мәндерді (5.6)-формулаға қойсақ:

Жоғарыда келтірілген мысалдармен салыстыру арқылы қай әдістің қателігі аз екенін анықтай аласыздар.


Бақылау сұрақтары:

  1. Сандық интегралдау есебіне мысал келтіріңіз?

  2. Трапеция әдісі қалай жүзеге асады?

  3. Симпсон әдісі қалай жүзеге асады?

  4. Тіктөртбұрыштар әдісі қалай жүзеге асады?

103


Тақырыбы: Болжау-түзеу схемасын пайдаланып КДТ сандық шешімін табу.

Мақсаты: Қарапайым дифференциалдық теңдеулерді шешудің бірқадамды сандық әдістерімен есептер шығаруды үйрену

Тапсырма:

Төмендегі қарапайым дифференциалдық теңдеулерді [0.2; 1.2] аралығында 0.1 қадаммен у(0.2)=0.25 бастапқы шартты қанағаттандыратын шешімін Эйлер және Эйлер-Коши әдістерімен тауып, қателіктерін бағалау. Есептеуді үтірден кейін 4 орынмен жүргізу.

№ 1.

№ 2.


№ 3


№ 4


№ 5


№ 6


№ 7


№ 8


№ 9


№ 10


2. Берілген аралықта h=0,2 қадаммен теңдеулерді және теңдеулер жүйелерін Рунге-Кутта әдісімен шешу.

a) y(0)=1.5; a=0; b=1;

b) , y(1)=1; a=1; b=2;

c) , y(0)=1; z(0)=1; a=0; b=1;

d) , y(0)=1, табу керек y(0.5)

e) , y(0)=1, табу керек y(1)

f) , y(0)=1, табу керек y(0.5)

g) y(0)=2, z(0)=-2, x=0.5 болғандағы y,z-тің мәндерін табу керек
i) y(0)=2, z(0)=-1, x=0.5 болғандағы y,z-тің мәндерін табу керек

3. Рунге-Кутта әдісін қолданып берілген аралықта 10-4 дәлдікпен төмендегі ҚДТ –ді шешу.

a) , y(0)=1, a=0, b=1

104


b) , y(0)=1, a=0, b=1

c) , y(1)=0, a=1, b=2

4. Рунге-Кутта әдісін қолданып бастапқы шарты х(0)=0 болғанда төмендегі теңдеулердің шешімдерін анықтау. һ=0,1 болсын.

a) , a=1+0.4n, n=0,1,…, 5, b=1+0.8k, k=0,1,2

b) , a=1+0.4n, n=0,1,…, 5, b=1+0.4k, k=0,1,2,3,4,5

c) , a=1+0.4n, n=0,1,…, 5, b=1+0.4k, k=0,1,2,3,4,5

d) , a=1+0.4n, n=0,1,…, 5

e) , a=1+0.4n, n=0,1,…, 5, b=1+0.8k, k=0,1,2

f) , a=1+0.4n, n=0,1,…, 5

g) , a=1+0.4n, n=0,1,…, 5

5. Рунге-Кутта әдісін қолданып 10-4 дәлдікке дейін төмендегі теңдеулерді шешу.

a) , y(1)=1, a=1, b=2,

b) , y(0)=1, a=0, b=1,

Тапсырманы орындауға әдістемелік нұсқаулар:

Эйлер әдісі

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

(6.1)-(6.2) Коши есебі берілсін.

Геометриялық мағынасы: һ қадам таңдап алып берілген аралықта бірдей қадаммен нүктелер жиынын құраймыз:

xi=x0+ih (i=0,1,2,…). (6.4)

M0(x0,y0) нүктесінен өтетін ізделінді y=y(x) интегралдық қисықты төбелері Mi (xi,yi) (i=0,1,2,…) болатын Эйлер сынықтарымен M0M1M2… алмастырылады:



(6.5)

Әрбір MiMi+1 сынықтары бағыты Mi нүктелерінен өтетін (6.1)-теңдеумен берілген интегралдық қисықтың бағытымен беттеседі.

Сонда есептеу формуласы келесі түрде жазылады:

Yi+1=yi+yi, (6.6)

yi=hf(xi,yi) (i=0,1,2,…)

Әдіс теңдеулер жүйесіне де қолданылады. Ол 2-мысалмен келтірілген.



1-Мысал:

Эйлер әдісін қолданып [0,1] аралығында h=0,2 қадаммен теңдеуін және у(0)=1 бастапқы шартты қанағаттандыратын мәндер кестесін құру керек болсын.

Есептеу нәтижелері 1-кестеде келтірілген. Кестенің толтырылуы: Бірінші жолға i=0 болғанда бастапқы мәндер жазылады: x0=0; y0=1,0000. Үтірден кейін мәнді цифрларды

жоғалтпау үшін 4 орын сақтап отырайық. Осы мәндер және (6)- формула бойынша f(x0,y0)=1, сосын , у1=1+0,2=1,2 мәндері есептеледі.



105

Екінші жолға i=1 болғанда x1=0.2, y1=1.2000 мәндері жазылады. Осы мәндерді қолданып f(x1,y1)=0.8667 мәні, мәні есептеледі. Сосын мәнін анықтауға болады. Дәл осылайша i=2,3,4,5 болғандағы есептеулерді анықтауға болады.

Кестенің ең соңғы бағанында салыстыру үшін теңдеудің дәл шешімінің мәндері келтірілген. Кестеден абсолютті қателіктің мәні е=0,0917, салыстырмалы қателігінің 5% екендігі көрінеді.



1-кесте . теңдеуін есептеу алгоритмі.

I

xi

yi

-ді есептеу
























0

1

2



3

4

5



0

0,2


0,4

0,6


0,8

1,0


1,0000

1,2000


1,3733

1,5294


1,6786

1,8237


0

0,3333


0,5928

0,7846


0,9532

1,0000

0,8667


0,7805

0,7458


0,7254

0,2000

0,1733


0,1561

0,1492


0,1451


1,0000

1,1832


1,3416

1,4832


1,6124

1,7320



2-мысал:

Бастапқы шарттары: .

Аралық [1, 1.5], қадам h=0.1 болсын. Мұндай жағдайда алмастыру қолдану арқылы 1-ретті теңдеулер жүйесін құрып алуға болады:

және бастапқы шарттары: .

Мұндай жүйелер үшін Эйлер формуласы келесі түрде жазылады:

Yi+1=yi+hf1(xi,yi,zi)

zi+1=zi+hf2(xi,yi,zi) , i=0,1,2,… (6.7)

(6.7)-ні қолдансақ және f1(x,y,z)=z; f2(x,y,z)= болатынын ескерсек:

i=0; x0=1.00; y0=0.77; z0=-0.44; f1(x0,y0,z0)=-0.44; f2(x0,y0,z0)= =-0.33;

y1=y0+=0.726; z1=z0+=-0.473 екендігі шығады. Дәл осылайша келесі мәндерді табу оңайға түседі.

Эйлер – Коши әдісі

Бастапқы нүктедегі ақиқат қисыққа жанама көлбеуі бұрышының тангенсі белгілі және -ға тең болса да, ол тәуелсіз айнымалының өзгеруіне байланысты өзгеріп отырады. Сондықтан x0+h нүктесінде жанама көлбеуі x0 нүктесіндегі жанама көлбеуіндей болмайды. Осыдан, h интервалында бастапқы жанама көлбеуін сақтай отырып есептеу барысында қателік пайда болады. Эйлер әдісінің дәлдігін арттыру үшін туынды аппроксимациясын жақсарту керек, яғни интервалдың бастапқы және соңғы нүктелерінде

туындының орта мәнін алуға болады. Бұл әдісті Эйлер – Коши әдісі дейді. Бұл әдісте алдымен Эйлер формуласы қолданылады:

106


, i=0,1,2,… Сосын осы мәнді интервал соңындағы туынды мәнін жуықтап есептеуге қолданады. Табылған екі мәннің ортасын анықтап, дәл мәнге өте жуық мән аламыз, бұл формула Эйлер – Коши формуласы деп аталады:

(6.8)

Рунге-Кутта әдісі

Бұл әдіс те бірқадамды әдіске жатады.



(6.9) (6.10)

(6.9)-теңдеу екі өлшемді қарапайым дифференциалдық теңдеу (ҚДТ) және (6.10)-бастапқы шарт берілсін. [x0, xn] аралығында у-тің мәндерін анықтап, функция графигін сызу керек болсын.



нүктелерінде мәндерін Рунге-Кутта формуласымен табамыз:

Мұндағы К аралық сандары төмендегідей табылады:









Практикада есептеу қадамдарын 2-кестеге толтырған тиімді.



2-кесте. Рунге-Кутта әдісінің схемасы.

i

X

Y

K=hf(x,y)



0


























1




























2





Т.с.с.




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

Бұл әдістің дәлдігі 4-ші ретті. Қателігін бағалау өте күрделі, сондықтан екілік есептеу арқылы алынған мәндер бір бірімен салыстырылып, айырмасы дәлдіктен асып кетпеуі

тексеріледі, егер айырма мәндері көп ауытқыса, есептеудің келесі қадамындағы х мәнін екі еселенген қадаммен алады, әйтпесе қадамның жартысын алады.

107


1-мысал: Рунге-Кутта әдісін қолданып теңдеуін шешу. Бастапқы шарты y(0)=-1, h=0.1, аралық [0,0.5] болсын.

Шешімі:


Есептеу қадамдарын 3-кестеге толтырып отырған дұрыс. Кестені толтыру ережесі:

  1. 1-ші жолға бастапқы шарт мәндерін жазамыз. Х0=0, у0=-1.

  2. f(x0, y0)=0.25, K10=0.1*0.25=0.025 мәндерін тауып аламыз.

  3. 2-ші жолға мәндерін жазамыз.

  4. , мәндерін есептейміз.

  5. 3-ші жолға мәндерін жазамыз.

  6. , мәндерін есептейміз.

  7. 4-ші жолға , мәндерін жазамыз.

  8. , мәндерін есептейміз.

  9. 5-ші жолға мәнін жазамыз.

  10. 6-шы жолға

  11. Осы алгоритмді қайталау арқылы кестенің қалған жолдары толтырылады.

Рунге-Кутта әдісін ҚДТ жүйесіне де қолдануға болады. Егер жүйе жоғарғы ретті болса, теңдеуді 1-ретті түрге келтіріп алу керек.

3-кесте . теңдеуін шешудің кестелік алгоритмі.



i

X

Y







0

0

0.05


0.05

0.1


-1

-0.98750


-0.98769

-0.97536













0.024

1

0.1

0.15


0.15

0.2


-0.97528

-0.96289


-0.96257

-0.94987












0.025

2

0.2

0.25


0.25

0.3


-0.94978

-0.93650


-0.93569

-0.92164


0.026553

0.028176


0.028138

0.030236


0.02824

0.023

3

0.3

0.35


0.35

0.4


-0.92154

-0.90642


-0.90514

-0.88881


0.030231

0.032790


0.032732

0.035743


0.03284

0.023

4

0.4

0.45


0.45

0.5


-0.88870

-0.87083


-0.86910

-0.84957


0.035745

0.039209


0.039134

0.04307


0.03925

0.022

5

0.5

-0.84945










108

Тақырыбы: Адамс схемасын екі теңдеу жүйесінің сандық шешімін табу Мақсаты: Қарапайым дифференцилдық теңдеулерді шешудің көпқадамды сандық әдістерімен есептер шығаруды үйрену

Тапсырма:

  1. Адамс әдісін қолданып 10-2 дәрежесіне дейінгі дәлдікпен төмендегі қарапайым дифференциалдық теңдеулер мен теңдеулер жүйелерін шешу. Бастапқы мәндерді Рунге-Кутта әдісімен есептеу.

a) , y(0)=1, табу керек y(0.5)

b) , y(0)=1, табу керек y(1)

c) , y(0)=1, табу керек y(0.5)

d) y(0)=2, z(0)=-2, x=0.5 болғандағы y,z-тің мәндерін табу керек

e) y(0)=2, z(0)=-1, x=0.5 болғандағы y,z-тің мәндерін табу керек

2. Адамс әдісін қолданып берілген аралықта 10-4 дәлдікпен төмендегі ҚДТ –ді шешу. Бастапқы мәндерді Рунге-Кутта әдісімен анықтау.

a) , y(0)=1, a=0, b=1

b) , y(0)=1, a=0, b=1

c) , y(1)=0, a=1, b=2

3. Адамс әдісін қолданып бастапқы шарты х(0)=0 болғанда төмендегі теңдеулердің шешімдердің мәндер кесте сын екі қадамға жалғастыру. һ=0,1 болсын. Бастапқы мәндерді Рунге-Кутта әдісімен анықтау.

a) , a=1+0.4n, n=0,1,…, 5, b=1+0.8k, k=0,1,2

b) , a=1+0.4n, n=0,1,…, 5, b=1+0.4k, k=0,1,2,3,4,5

c) , a=1+0.4n, n=0,1,…, 5, b=1+0.4k, k=0,1,2,3,4,5

d) , a=1+0.4n, n=0,1,…, 5

e) , a=1+0.4n, n=0,1,…, 5, b=1+0.8k, k=0,1,2

f) , a=1+0.4n, n=0,1,…, 5

g) , a=1+0.4n, n=0,1,…, 5

4. Милн әдісін қолданып 10-4 дәлдікке дейін төмендегі теңдеулерді шешу. Бастапқы мәндерді бірқадамды әдістердің біреуімен анықтау.

a) , y(1)=1, a=1, b=2,

b) , y(0)=1, a=0, b=1,



5. Адамс әдісін қолданып дифференциалдық теңдеудің бастапқы шартты қанағаттандыратын шешімдерін аралығында қадаммен анықтау. Есептеуді төрт ондық таңбамен жүргізу. Бастапқы аралықтағы мәндерді Рунге-Кутта әдісімен анықтау.

109























Тапсырманы орындауға әдістемелік нұсқаулар:

(6.11)

(6.12)

(6.11)-теңдеу екі өлшемді қарапайым дифференциалдық теңдеу (ҚДТ) және (6.12)-бастапқы шарт берілсін. [x0, xn] аралығында у-тің мәндерін анықтап, функция графигін сызу керек болсын. Бұл есепті шешудің көпқадамды әдістері: Адамс және Милн әдістері деп аталады.



Адамс әдісі

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



Бастапқы шартты пайдаланып, функция өсімшесінің мәндерін анықтаймыз, оларды qi (i=0,1,2,3) деп белгілейік:







Енді осы мәндердің шектік айырымдарын табамыз:







,

Енді функцияның мәндерін есептеу үшін Адамс формуласын қолданамыз. Ол екі түрлі:



  1. Экстраполяциялық формула:

,

, (6.13)

110


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

  1. Интерполяциялық формула :

,

, (6.14)

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



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

Практикада есептеуді жеілдету үшін Адамстың басқа формулалары да қолданылады:

1-формуласы: ), (6.15)

2-формуласы: ), (6.16)



Адамс әдісі дифференциалдық теңдеулер жүйесіне де қолданылады:



жүйесі берілсе, оған қолданылатын Адамс формулалары төмендегідей болады:

(6.17)

(6.18)

Мұндағы: ,





Милн әдісі.

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

(6.19)

(6.20)



,

мәндерді анықтау үшін Милн формулаларын қолданамыз:

  1. Алдын ала анықтау: ,

  2. Осы мәндерді қолданып




  1. Милннің 2-ші формуласымен алдында табылған мәндерді түзетеміз немесе дәлдейміз: .

111

  1. Табылған мәндердің қателігін бағалаймыз: . Бұл формула есептеудің әр қадамында алынған мәннің дәлдігін тексеріп отырады. Егер дәлдік берілсе және болса, онда деп алып 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

Есептеу қадамдарын 1-кестеге жазуға болады.

Кестені толтыру ережесіне тоқталайық:

1,2- бағандарға белгілі мәндерді толтырамыз. 3-бағанда уk-дің (k=0,1,2,3)

белгілі мәндерін толтырамыз. Осы мәндерді қолданып,

1-кесте . теңдеуін шешудің алгоритмі.



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)- формуламен түзетеміз:

. Табылған екі мән бір біріне өте жуық болғандықтан кестедегі алдыңғы табылған мәнін мәнімен түзетеміз. Осы әдіспен кестенің келесі жолын толтыруға болады.

Бақылау сұрақтары:

  1. Қарапайым дифференциалдық теңдеулерді шешудің көпқадамды сандық әдістерін атаңыз?

  2. Адамс әдісі?

  3. Милн әдісі?

  4. Қарапайым екіөлшемді дифференциалдық теңдеу қалай жазылады?

112

МАТЕРИАЛДЫ МЕҢГЕРГЕНДІГІН БАҒАЛАУҒА АРНАЛҒАН СҰРАҚТАР


  1. Абсолютті және салыстырмалы қателіктер

  2. Қателіктерді бағалау әдістері

  3. Қателіктер теориясының тура есебі

  4. Қателіктер теориясының кері есебі

  5. Қателіктерге қолданылатын арифметикалық амалдар

  6. Функциялардың қателіктері

  7. Сызықты емес теңдеулердің түрлері

  8. Сызықты және сызықты емес теңдеулер туралы ұғым.

  9. Аралықты қақ бөлу әдісі

  10. Хорда әдісі

  11. Біріктірілген хорда және жанама әдісі

  12. Сызықты және сызықты емес теңдеулер жүйесі

  13. Теңдеулердің итерациялық түрі

  14. Теңдеулер жүйесінің итерациялық түрі

  15. Теңдеулерді және теңдеулер жүйесін итерациялық түрге келтіру әдістері

  16. Алгебралық және трансцендентті теңдеулер туралы ұғым.

  17. Сызықты алгебралық теңдеулер жүйесі

  18. Матрица және вектор ұғымы.

  19. Меншікті мәндер және меншікті векторлар

  20. Норма ұғымы. Метрикалық кеңістіктер.

  21. Сандық әдістердің нормада, метрикалық кеңістікте жинақтылығы

  22. Итерациялық әдістердің жинақтылығы

  23. Итерациялық процестерді құру

  24. Сандық интегралдау есебі.

  25. Интегралдау қадамын таңдау.

  26. Интегралдау қателіктерін бағалау

  27. Канторович әдісі және қателігі

  28. Гаусс әдісі және қателігі

  29. Қарапайым дифференциалдық теңдеулерді шешудің бірқадамды әдісі

  30. Қарапайым дифференциалдық теңдеулерді шешудің көпқадамды әдісі

  31. ҚДТ-ны шешу әдістерінің жинақтылығы

  32. Бастапқы және шекаралық шарттардың берілуі және айырмашылықтары

  33. Липшиц теоремасы және оның дифференциалдық теңдеулер шешіміне қатысы

  34. Шектік есептер ұғымы.

  35. ҚДТ үшін шектік есептердің қойылуы

  36. ҚДТ үшін шектік есептердің шешімінің жинақтылық шарттары

  37. ҚДТ үшін шектік есептерді шешудің итерациялық әдістері

  38. ҚДТ үшін шектік есептерді шешудің вариациялық әдістері

  39. Дербес туындылы дифференциалдық теңдеулердің классификациясы

  40. Гиперболалық типті дербес туындылы дифференциалдық теңдеулер және оларды шешу әдістері

  41. Параболалық типті дербес туындылы дифференциалдық теңдеулер және оларды шешу әдістері

  42. Эллипстік типті дербес туындылы дифференциалдық теңдеулер және оларды шешу әдістері

  43. Дербес туындылы дифференциалдық теңдеулерді шешу әдістерінің жинақтылығы мен орнықтылығы

113


Қорытынды

Сандық әдістер тарихы. Математикалық моделдеу мең есептеу(Қолданбалы есептерді шығару кезеңдері. Есептеу эксперимеңті сүлбесі. Есептеу алгоритмі. Сандық әдістер мең математикалық моделдеуге қойылатын негізгі талаптар.) Есептеу математика пәні. Есептеу информатикасы жайлы. Есеп шешімі қателерін жіктеу. Үнемді, орнықты алгоритмдер. Орындылық Аппроксимация (дискреттеу). Есептеу алгоритмдерін жүзеге асырудың инструмеңтальдық кұралдары. Абсолют және салыстырмалы қателіктер. Берілгеңдерді жазу түрлері. Функция кателігі. Түбірлерді бөлу есебі. Дихотомия ( қақ болу) әдісі. Жай итерация әдісі. Ньютон ( жанамалар әдісі) әдісі. Кесу сызыктар әдісі. Үздіксіз Ньютон әдісі. Итерациялык әдістердің геометриялық түсініктемесі..



Біртіндеп Гаусс әдісі. Әдістің көбейткіштерге жіктеумең байланысы. Бас элемеңтті таңдау арқылы орындалатын Гаусс әдісі. Бірнеше теңдеулер жүйесін шешу. Матрицаның анықтауышын есептеу және кері матрицаны есептеу. Басқа дәл әдістер жайында. Үш диогональды матрицамең берілгең сызықтық жүйеңі шешуге арналған қуу әдісі. Квадрат түбір әдісі. Сызықтық алгебрадан косымша мәлімет. Матрица нормасы. Матрицалық геометриялық прогрессияның жинақтылығы (қажетті және жеткілікті шарттар). Қарапайым итерация әдісі. Қарапайым итерация әдісінің жинақтылығының жеткілікті шартының теоремаеы. Зейдель әдісі. Зейдель әдісінің жинақтылығының қажетті және жеткілікті шарттары.Зейдель әдісінің жинактылығының жеткілікті шартының теоремасы. Нашар шарттасқан жүйелер туралы. Матрицаның шарттасқан өлшемі. Нашар шарттасқан жүйелер мысалы. Есептің қойылуы. Дәл әдістер туралы. Өзіндік мәндердің дербес мәселелері. Өзіндік мәндердің абсолют шамасы бойынша ең үлкеңін және оған сәйкес келетін өзіндік вектор табу. Симметриялык матрицадағы жағдай. Кейбір күрделі жағдайлар. Еселі басып шығатын өзіндік мәннің жағдайы. Екі ең үлкең таңбалары бөлек өзіндік мәндердің жағдайы. Басқа әдістер туралы. Симметриялык матрицалардын айналдыру әдісі үшін өзіндік мәндердің толык мәселелерін шешу.

Есептің қойылуы. Функцияны интерполяциялау есебі. Лагранж интерполяциялык көпмүшелігі. Лагранж интерполяциялык формуласының қалдық мүшесі және оны бағалау. Чебышыев көпмүшелігі. Эйткең есептеу схемасы. Бөлінетін айырымдар және оның қасиеттері. Ньютон интерполяциялык көпмүшеліктері. Басқа жуықтаулар туралы. Сплайндармең интерполяциялау. Ең кіші квадраттар әдісі. Бірқалыпты жуықтаулар Сандық дифференциалдау есебінің шартты қисынлылыгы. Сандық интегралдау. Ньютон-Котес квадратуралық формулалары. Тік төртбұрыштар, трапециялар, Симнсон формулаларының қалдық мүшелерін бағалау. Гаусстың квадратуралық

формулалары-ең жоғары алгебралық дәлдікті интерполяциялық квадратуралар. Сандык интегралдаудьщ формулаларының бағалау қателіктерінің практикалык қолданылуы. Рунге ережесі. Қарапайым дифференциалдық теңдеулер үшін Коши есебін сандық шешудің Эйлер әдістері. Эйлер әдісінің модификациялау. Рунге-Кутта әдісі. Екінші ретті сызықтық шектік есептерді шешудің қуу әдісі. Айырымдық аппроксимациясы. Жылу өткізгіштік және ішек тербелісінің бір өлшемді тендеулер үшін үнемді айырымдық сұлбаларға мысалдар. Жылу өткізгіштің сызықты емес теңдеулері және олар үшін айрымдық сұлбалар. Айрымдық сұлбапар үшін максимум принципі. Пуассон теңдеуі үшін Дирихле айырымдық есебінің орнықтылығы және жинақтылықгы зерттегендегі айнымалыларды бөлу әдісі. Торлы шектік есептерді шешудің итерациялық әдістері қарастырылған.


114

Пайдаланылған әдебиет тізімі



Негізгі адебиеттер:

  1. Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков. Числеңные методы, М.: Лаборатория базовых знаний, 2002.-632 с.

  2. Н.С. Бахвалов, А.В. Лапин, Е.В.Чижонков. Числеңные методы в задачах и упражнеңиях. М.: Высшая школа, 2000

  3. В.М. Вержбицкий Числеңные методы. Линейная алгебра и нелинейные уравнеңия, М.: Высшая школа, 2000.-266с.

  4. В.М. Вержбицкий Числеңные методы. Математический анализ и обыкновеңные диффереңциальные уравнеңия. М.: Высшая школа, 2001.-382с.

  5. Джон Г.Мэтьюз, Куртис Д. Финк Числеңные методы. Использование MatLab, 3-издание.: Пер. С англ.-М.:Издательский дом «Вильяме», 2001.-720с.

  6. В.И.Киреев, А.В. Пантелеев. Числеңные методы в примерах и задачах, М.: «Высшая школа», 2004.-480с.

  7. М.П.Лапчик, М.И. Рагулина, Е.Л..Хеңнер. Числеңные методы, М.: Издательский цеңтр «Академия», 2004.-384с.

  8. СВ. Поршнев Вычислительная математика. Курс лекций.-СПб.:БХВ- Петербург, 2004.-320с.

  9. С.В Поршнев, ИВ. Белеңкова Числеңные методы на базе MathCAD.- СПб.:БХВ-Петербург, 2005.-464с.

  10. Д.К. Фаддеев, В.Н. Фаддеева Вычислительные методы линейной алгебры. СПб.: Издательство «Лань», 2002-736с.

  11. А.И.Плисс, Н.А.Сливина MathCAD 2000. Математический практикум. М.: Финансы и статистика, 2000.-655с.

  12. Н.С. Бахвалов, А.В. Лапин, Е.В.Чижонков. Числеңные методы в задачах и упражнеңиях, М.: Высшая школа,2000.

Қосымша адебиеттер:

  1. Бахвалов Н.С., Жидков Н.П. Числеңные методы - М., 1987

  2. Бахвалов Н.С. Числеңные методы - М., 1975

  3. Самарский А.А. Введеңие в числеңные методы,- М., 1987

  4. Самарский А.А., Гулин А.В. Числеңные методы - М., 1989

  5. Солтангазин О.М., Атанбаев С.А. Есептеу әдістерінің қысқаша теориясы 1-і кітап.-Алматы,2004.-384

  6. Волков Е.А. Числеңные методы.-М.:Наука,1987.-248 с.

  7. Марчук Г:И: Методы вычислительной математики.-М.:Наука,, 1980.-536 с.

  8. М., 1989.-608 с. М„ 1977.-456 с.

  9. Білімді бағалау туралы анықтама


115

Достарыңызбен бөлісу:
1   ...   32   33   34   35   36   37   38   39   40




©emirsaba.org 2024
әкімшілігінің қараңыз

    Басты бет