Файл: Ермолаева Э.Н. Элементы численного анализа учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 07.07.2024

Просмотров: 139

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

 

г

 

 

4

4

0 =

& ёгч

 

 

 

0 = 4 4 4

4

Контролем правильности вычислений коэффициентов этой таблицы будет служить выполнение следующих равенств:

Ьхь

=

b i 2

-

Ь г ъ

+

Ь

н — \ ;

4 г ,

= — b-,0

b23

- f

Ъ,Л ;

*вз

=

- ь г г

-

& з з

+

Ь

и .

Если Ь 2 2 Ф 0 , то, опять сделав один шаг жордановых исклю­ чений с этим разрешающим элементом, получим табл. 3.15.

Таблица 3.15 Таблица 3.1в

 

~Х<

1

 

 

х-

с ,»

C,s

1

 

 

0 =

Css

4 4

Контрольные формулы правильности вычислений, коэффи­ циентов этой таблицы будут иметь вид:

С15 =

С13 + С Н

1 ">

с2-, = • • С-23 -(- с,4

1 ;

СЗо — .

С 33 "Г С З І •

 

Над табл. 3.15 произведем шаг жордановых исключений с разрешающим элементом с13Ф0 и получим табл. 3.16.


Б этой таблице должно выполняться условие:

^25 = du

1 . d.,:> du

1

,

d3r> =

d3i - -

1 .

Решением

системы

(3.6)

будут

 

числа

X | = £ / M ,

^2 = ^/24,

л:з = ^з4, а контрольные

числа di5

= xi—1, dys — Хг—1. #г5 = *з—1

дадут решение контрольной системы.

 

 

 

 

 

 

 

.

 

 

 

Таблица 3.17

Решение

 

 

 

 

1

о -

 

 

а

У,

0 =

<*„

а

Уг

 

 

 

гз

0=.

а»

 

Уз

 

 

 

- Х Г

 

-*з

1

 

Зи

 

 

 

 

О -

 

 

4 ,

 

 

 

 

 

0 = 4

 

 

 

 

 

 

 

/

 

X

=

 

г

 

1

 

 

 

 

 

 

 

с

 

0

-

 

Г

 

 

 

 

 

 

 

 

і

 

 

 

 

d l 4

 

 

 

^3 =

< 4

 

Контроль

 

/

б; —

агГагГагіг

 

1

1

- c „ +

c H - 1 )

f)

 

1

 

 

 

\

 

 

 

 

 

 

d

l

4

-

1 )

 

d

i

4

-

0

d 3 S

( = < 4 "

1 )

Схема решения системы (3.5) приведена в табл. 3.17 в пред­ положении, что система имеет единственное решение.


Пример 3.5. Решить систему

 

 

 

 

 

 

10 я,, — 5,80.al ~ 21,74.00 л, -

0,33 = 0 ,

 

-

5,80.я0 4- 21,74 • 00 ал

-

42,1 • 66 а,

-

6,5.26 =

0 ,

 

21,74.00 а0 -

42,1.66 а,

4- 121.5.01 а,

-

17,1.55 =

0 .

а,, Ф 0,. Ь.„ Ф 0,

с; ) 3 Ф 0 — разрешающие элементы. Решени­

ем СИСТеМЫ будут ЧИСЛа X\=d\^

Х2 = ^24, Х 3 = іІз4.

 

значащие

Замечания. 1. В условии точкой отделены верные

цифры от остальных.

 

 

 

 

 

 

 

2. Напомним лравило подсчета цифр; при умножении и де­

лении следует сохранять в результате столько значащих

цифр,

сколько их имеет приближенное данное с наименьшим

числом

верных значащих цифр.

 

 

 

 

 

 

 

В промежуточных результатах каждый раз оставляют две

запасные цифры.

 

 

 

 

 

 

 

 

Решение оформлено в виде табл. 3.18. Ответом будут слу­

жить

числа: а0 = 1,0.215, а\ = 2,0.07, а 2 = 1,0.204.

 

 

 

§ 3.5.

М Е Т О Д ГАУССА

 

 

 

 

 

 

 

Часто описываемый в литературе метод Гаусса

отличается

от предыдущего

лишь тем, что после каждого

шага

жордано-

ьых исключений вычеркивают из таблицы разрешающую стро­ ку, выписывая ее отдельно.

На последнем шаге начодят значение одного из неизвест­ ных, а значения остальных вычисляются последовательной подстановкой уже полученных значений неизвестных в выпи­

санные выражения для

х,.

 

 

 

 

Пример

3.6. Решить

систему

примера

3.2 методом

Гаусса,

осуществляя контроль.

 

 

 

 

 

Р е ш е н и е .

Составим табл. 3.19.

 

 

 

 

 

 

 

'

• •"' ' Т а б л и

ц а 3.19

 

 

 

 

'•2 9

і

Нон тр олЬ

 

0=

г

г

-/

га

ч

 

 

 

0 =

Ч

3

- /

г

в

 

 

 

0=

8

5

ч

-2 =

/2~8-5-(-3)-Ч

0=

3

3

г

6

о=б-з-з-ы)-г

 

)


 

 

і'о=

юав

-

5,80, а, + 21,74.00аг-

 

о_зз.

 

 

 

ТйБПиЦб

3J8

 

 

 

 

 

 

 

 

 

 

1 0=

-5jg0,a0

+21J9.Q0CL, -

42,1.66

а г

-

6,5.26

 

 

 

 

 

 

 

 

о=21,74.зоас-

42,і.бба(ч-121,501

 

аг

-

17,1.55

 

 

 

 

 

 

 

 

 

-а.

 

 

 

і

 

~ а <

 

 

 

 

С 8. чл.

 

 

\і

0=1

 

Е Ш 2

 

 

 

!

 

5,80.

 

-21ЦОС

-о,зз.

25, 6/. 00

 

|

 

5, SO

 

 

[

- гі,

г</. оо

 

4£1.66

-6,5.26

-32, 75.2

!

о-- І

-гі,7Ч.оо

 

|

 

42, 1. 66

 

-12

5,01

-17,1.55

85,3.20

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Св. чл.

 

 

 

 

<*о =

 

 

5,80. (-Ю)

 

=-0,58.0

 

\ (-21,74.00):(-ю)=2,174.00

(-0,35): но) ~ о,озз.

25,6/,00:(-/о)= -2, 56/. оо

о

 

 

3J£4G

 

 

 

\

 

 

 

 

-гге.,09

 

-is. іч _

 

I4B..S4

 

 

-21,79.00-15,80. • 5,80)^-/0)=-13,37.601

42,166-[(-21,74.00).5Mf-/0)^2S,5-.6,SІ'•.26-0,33) • 5,80]:(-/0)= -6,717

 

 

 

 

 

 

 

^1)36,40

 

I

 

 

 

 

f2,6.0S

 

0,19.14

 

-32/5,2-(2S,t>/.00-5JO):(-10}r-rZ!.n

 

 

 

-/26.,09

 

 

'

 

 

 

 

чгг^е

 

Г,1.74

 

 

 

 

 

 

 

 

 

-121,5.оі-[(-гі7Щ(-гід^о)$гЩ -м**$і5з)НШ«№Ф-/б. чм

 

 

 

 

 

42, Ш-[5,80 {-21,74.00)];Н0)=2Э,5,57

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12,6.09

 

 

 

 

 

 

-4~Z26\22

Ґ

-0,71.74

 

&3?20-[25.Ш0-(-2/, 74.00)):(-Ю)=282.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-а,

 

 

 

 

 

 

 

 

 

С

чл.

 

 

zz

 

 

 

 

 

J 4,2.5?

 

 

 

 

 

 

 

- SS. 7.S<

 

 

€1,4.02

 

 

 

 

-0,5,80 -(29,557- 2,17400)

.(-74,2.38)=0,285.55

0 033. -[(-16,4.38)•<,

74,00):(-74,2.38)=- 0,448.Зг -2,56/. 00-(28,2, 44-2,174.00): (-74^.38;= -£733.90

 

 

 

 

 

-0J6S.S5

 

 

 

 

 

-485.,86

 

 

 

-0,127/0

 

 

 

 

 

6ГЛ.,62

 

 

 

 

 

 

 

 

 

 

S34..8I

 

 

0 = -18,37.60-(29. 557-29, 557): (- 7fJ. 38) =(-

6,Ш

 

-<$, 7l7-[(-/6,4.38}29 iS7]:(~

74,2.38)=-15,2.62

-17,898-(28,2 44-25,5,5?)/(-74,2.38)=

-6,6.55

 

 

 

 

 

-117.68

 

 

 

 

 

 

 

 

 

 

-112.45

 

 

°*'

29, S. 57; (-74, 2.38)=

- 0,398.14

 

 

 

 

(-16, Ч. 38): (- 74,2 58) = 0,221. *2

28,2.44:(- 74,2.38) = -0,380.45

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

70

 

С 8.

чл.

 

 

 

 

-11.16

 

 

 

 

 

 

 

-S,

?в,

 

 

 

 

 

 

 

 

 

 

 

Ц А

 

-0,448.37-[(-13,2.62)•0,285.55):(-6,6.08) =

-1,02/5

-/,753.90-1(-6,6.53У 0,285.55J/(- 6,6.08) =-2,00.78

 

 

I

 

 

 

 

 

 

o7sr.s/~

 

 

 

 

 

 

 

0,27.39

 

 

41

 

 

(-15,262): (-6,6.08) =

2,о.ог

 

 

(- 6, 6.53) : (- 6, 6.08)-

1,0.07

 

 

 

 

 

 

SJH.Q1

 

 

 

 

 

 

 

 

 

2,

6.49

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Q221.42-lc/J,2.62X-QJS8,r4)J^6,6.08)= і 02.04

(- 0, J80.45)-[(-6,6.5JJ-(- ft398.ту. (-ЄЛ08)=0,02.04

 

 

 

 

 

 

-0,?9.90

 

 

 

 

 

 

 

 

-0,40.09

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 


Сделав шаг жордановых исключений с разрешающим эле­ ментом ЙІ4=1 и вычеркнув разрешающую строку, получим табл. 3.20.

 

 

 

 

 

 

Таблица 3.20

 

0

-/

 

Кон

т

роль

 

 

-г=о-(-о-1+ы)

0 = 0

- J

 

-V - 2 =

 

0-(-Ъ)-1+(-Ч)

0 -

-/

-/ О

1

С S• o =

 

-(-t)-(-l)-0+(-2)

 

 

 

 

 

 

Записав

отдельно

разрешающую

строку 0 = 2х\—2хч +

+ х3—Х4

+ 4, из нее получаем

 

 

 

 

 

 

2 х. 2 х2

+ х3-\- 4

Следующий шаг произведем с разрешающими первой стро­ кой и третьим столбцом. После вычеркивания разрешающей строки 0 = х2x-i—2 получим табл. 3.21.

 

 

 

 

Таблица 3.21

-X, ~У-г

J

 

Контроль

і

 

о - 2

0-

 

о - ( - 2 ) + ( - 2 )

0 = -/E7J - 2

0

=

-(-/)-(-t)+(-2)

Из разрешающей строки имеем

х>) - 2 •

Кпоследней таблице применим шаг жордановых исключе­ ний с разрешающими второй строкой и вторым столбцом. Окончательно получим табл. 3.22.

Таблица 3.22

 

-Xj

1

Контроль

о

2

I 2

ІО--2+2