Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Практические занятия по высшей математике. Численное решение алгебраических и трансцендентных уравнений, матричное исчисление, векторный анализ и интегрирование линейных дифференци

.pdf
Скачиваний:
11
Добавлен:
12.11.2023
Размер:
15.42 Mб
Скачать

а 22 — ^ ^ 1

= ^ (1

а 23^32

а 2 4 ^ « )

— я [ 1 - и ( - т ) - ( - т ) < - 13)

Ё “2*С*1

»- 2

а 21

си

 

^22^21 ’Н" а 23с 31 + а 24С41 =

6*11

 

 

 

-А-»+5-*+(-4)*

5

 

 

 

4

 

 

 

 

й

:

 

 

 

 

 

 

 

 

 

а 12 = —

X

&1*а *2 = —

1 2 * 2 2 + ^13*32 +

Ь ц

Л к ) =

 

= — [ т ( ~ й ) + т ( — ® ) + 2 • ° | = й ;

 

 

4

 

 

 

 

 

 

а 11 =

~ | 1

а 1*С**] =

^ ----а 12^21

а 13^31 -----а 14^411 =

 

- ■ г [ ' - я - ® - ( - 8 ) - » - 1

- * ] —

А -

Итак,

обратная матрица

 

 

 

 

 

 

 

 

I

5

23

1

 

 

 

 

15

54

“ Я

У

 

 

 

А~1=

5

 

67

3

 

 

 

Я

48

45

т

 

 

 

 

 

7^

13

7

 

 

 

 

 

 

24

45

48

4

 

 

 

 

 

О

О

- 2

1

 

 

 

 

 

— 4

10

—46

24"

 

 

 

_1

1 0 — 1

67

-36

 

 

 

 

48

— 1 4 — 13

7

12

 

 

 

 

 

О

0

—96

48

 

 

П р о в е р к а .

Должно быть

А~1 А1=

Е. Дейетвите

 

—4

10

— 46

24

1

Ю

— 1

 

67

—36

48

— 14

— 13

7

12

 

0

9

 

—96

48_

 

 

48

0

0

0 "

 

1

0

48

0

0

 

= 45

0

0

48

0

 

_

0

0

0

48

"2

7

1

 

4"

 

5

2

0

 

— 1

 

3

4

2

 

1

 

о>

00

 

 

со

 

1

 

 

__1

 

1

0

0

0

 

 

0

1

0

0

=

Е.

0

0

1

0

 

 

„0

0

0

1

_

 

Повторяем, что все вычисления были выполнены в простых дробях для более простой проверки произведенных действий. На практике вычисления обыкновенно ведутся в десятичных дробях, но вносимые округления влекут за собой неизбежные погрешности.

Задача 6,13 (для самостоятельного решения). Пользуясь фор­ мулами (6,26) — (6,28) и указанной схемой получения обратной матрицы, обратить следующие матрицы:

 

л*

 

О

1

А

 

 

0,888

1,059

0,667

1) А =

1

0

 

5

—3

 

2) А

 

:

1,059

1,512

1,031 .

4

 

2

1

1

 

 

 

 

0,667

1,031

0,837

 

О

 

Л

О

С

 

 

Вычисления производить в десятичных дробях с помощью ариф­ мометра или клавишной вычислительной машины. Полученные результаты проверить: произведение А~1 • А должно быть равно

единичной матрице

А~1- А - Е.

О т в е т .

А"1=

1) Для

контроля

 

 

 

 

 

2

0

 

0

0

 

С =

1

- 1 . 5

 

0

0

ф

4

—4

— 15

0

9

 

 

 

_3

- 0 ,5

—9

 

2,4667_

 

 

~ 1

1,5

2

-

1

 

В =

0

1

- 2

 

1,3333

 

0

0

1

—0,6889

9

 

_0

0

0

 

1

 

—0,0538

—0,1168

+0.4414

-0,1804

+0,4052

—0.2884

—0,1441

0,0181

—0,2163

+0,5314

—0,2342

0,2793

0,3784

+0,5134

—0,2432

0,4054

2)

8,08977

—7,93495

3,32745'

А- 1 =

-7,93495

11,91472

—8,35304

л

 

3,32745

-8,35304

8,83220

 

С о д е р ж а н и е . Матричная запись системы линейных алгебраических урав­ нений. Численное решение линейных алгебраических уравнений способом исключения.

СВЕДЕНИЯ ИЗ ТЕОРИИ

Система п линейных алгебраических

уравнений

с

п неизвест­

ными имеет такой

вид:

 

 

 

 

 

 

0 Ц * 1

+

^ 12*2 +

0 1 3 *8 +

’ •1 +

а 1пХ п

^ 1,

 

 

021*1 + 022*2 + 023*3 +

' '*+ а2п^п = <2>

 

031*1 +

032*2 +

а з3*3 +

• ‘‘ +

03л*л =

^3»

(7,1)

0 Я1*1 +

0/12*2 +

0/13*3 +

’ • ' + 0 /1 л * л =

« V

 

Мы будем

рассматривать только такие системы линейных алге­

браических уравнений,

у

которых

число

уравнений

равно числу

неизвестных.

В

системе

уравнений (7,1)

неизвестными являются

*1, хг, х3, ...

х„,

а

коэффициенты

а{/ (*,

^ =

\, 2,

3 , . . . п) при

неизвестных

и свободные

 

члены

й{ (1 =

1 ,

2,

3 , . . . п) — действи­

тельные числа.

 

 

А из коэффициентов о/;- при неизвестных

Составим матрицу

 

 

 

 

012

013

 

 

 

 

 

 

 

 

^21

^22

^23

 

 

 

 

 

 

х

 

- ап\

ап2

апЗ • • •

а1

 

 

матрицу-столбец

из

неизвестных

 

 

 

 

 

 

 

 

 

 

X =

 

хв

 

 

 

 

и матрицу-столбец й из свободных членов

С^2

й =

Учитывая правило

умножения

матриц

и

условие равенства

двух матриц, систему (7,1) можно записать

в виде

0ц

а 12

а13. . . аХп

XI

 

~ёГ

021

022

^23 • • • а2п

хг

ё2

031

032

0зз • • • 0зп

х?

й?

—^т

&п2

0/»з • • • ^/1я—

 

 

- Л _

или, учитывая введенные обозначения, всю систему уравнений (7,2) можно записать компактно в виде одного матричного уравнения

Ах = ё.

(7,3)

Такая запись большого числа линейных уравнений в виде одного уравнения является одним из достоинств матричных обо­ значений.

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

раических уравнений.

 

 

 

 

 

 

 

Если матрица

А — неособенная, то

она

имеет

обратную мат­

рицу А~У Умножим слева обе части

уравнения

(7,3) на

А~1 и,

заметив, что

А~* ■А = Е— единичная матрица,

получим

столбец

неизвестных

х из равенства

 

 

 

 

 

 

 

 

 

 

 

х = А~1й.

 

 

 

 

 

(7,4)

В этой формуле х может быть

не только

матрицей-столбцом, но

и матрицей

размера

п х т. В этом случае

и

матрица ё состав­

ленная

из свободных

членов, должна

тоже иметь размер п х т.

Такой

случай имеет место, например,

в задаче (7,3).

 

 

Таким образом,

нам стоит

только определить

элементы а,-/

(|, / — 1, 2

, ,

п) обратной матрицы

А~1,

и задачу

определения

неизвестных систем (7,1) можно считать решенной. Несколько таких упражнений мы выполним в начале этого практического занятия.

Но здесь же следует заметить, что в случае большого числа л неизвестных вычисление элементов обратной матрицы становится громоздким и затруднительным. Поэтому формула (7,4) имеет больше теоретическое значение, чем практическое, так как по сравнению с формулами Крамера для решения системы линейных алгебраических уравнений она никаких преимуществ в вычисле­ ниях не дает.

Эта формула оказывается безусловно полезной тогда, когда рассматриваются такие системы уравнений, у которых матрица коэффициентов при неизвестных одна и та же (таксе системы уравнений встречаются, например, в строительной механике). Вы-

числение обратной матрицы в таком случае приносит большую экономию в вычислительной работе (см., например, задачу (7 ,3).

Это практическое занятие проведем в таком порядке: сначала будем решать системы уравнений с небольшим числом неизвестных (два, трй) по формуле (7,4) (задачи 7,1—7,7), а потом укажем удобную компактную схему решения системы (7,1) методом исклю­ чения (алгоритм Гаусса), которому дадим матричную трактовку. Теория этого метода, формулы, сюда относящиеся и схема его применения указаны после задачи 7,7.

Задача 7,1. Записать систему уравнений

2^ + 3*2 = 9) Зд.| 4*2 = 5

ввиде одного матричного уравнения и решить ее по формуле (7,4).

Ре ш е н и е . Систему представим в виде

Здесь матрица

Найдем обратную ей матрицу по формуле (5,2), которая в развернутом виде выглядит так:

(7,5)

где | — определитель матрицы Л, а элементы А{1— алгебраиче­ ские дополнения ее элементов а„, причем матрица в правой части

формулы (7,5) есть союзная матрица А для А. Но здесь еще раз подчеркнем, что когда п > 3, то формулой (7,5) для обращения матрицы обыкновенно не пользуются, а применяют методы, рас­ смотренные на пятом и шестом практических занятиях.

Воспользуемся формулой (5,2) для определения обратной матрицы

\А\=

I

 

 

= - 1 7 ;

 

 

Аи = (— 1)1+1 • ( - 4 ) =

- 4

;

Аи =

( - 1)1+*

3 =

—3;

Ап = (-1)*+* 3 =

-

3;

 

Ап =

(-!)»<■*

(2) =

2.

Составим

матрицу,

союзную матрице А

 

 

 

 

 

 

 

11

^21

 

Г -4

- 3

 

 

 

 

 

 

 

^12

^22

I - 3

 

 

 

 

 

 

 

 

 

 

 

 

 

‘ 4

3 ‘

 

 

 

 

=

1

 

- 4

 

- 3 1

_

17

17

 

 

 

 

 

—17 —3

 

2 I—

3

2

 

 

Теперь по формуле

(7,4)

 

 

 

 

 

-Т7

Т7.

 

 

х= А ~ Ч .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'4

 

3 '

 

’9‘

 

 

 

 

 

 

М

-

17

Т7

 

 

 

 

 

 

 

3

 

—2

 

5

 

 

 

 

 

 

N

 

-17

17.

 

 

 

 

 

 

 

 

 

 

 

51

 

 

 

 

 

 

 

 

 

 

н

 

17

 

 

 

 

 

 

 

 

 

 

 

- -17Г7.

 

 

 

 

 

 

 

Х\ _51 __ о.

^2 17 .1-

 

 

 

Задача 7,2 (для

самостоятельного

решения).

Записать систему

уравнений

 

 

 

 

5*1 +

 

хг =

7;

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*1

 

3*2 =

8

|

 

 

 

в виде одного матричного уравнения и решить

ее по формуле (7 ,4)

У к а з а н и е .

Обратная

матрица по формуле

(5,2)

 

 

 

 

А~* =

ТВ

4

 

 

 

Л

*, =

29

 

 

 

 

33

 

 

 

О т в е т .

]$;

*2 = — -ре.

 

 

 

 

 

Задача 7,3.

Решить системы уравнений:

 

 

 

 

1) 3*, — 2л:2 = 2;

 

| 2)

З* 3 — 2* 4 =

3; 1

 

 

4*1 — * 2 = 1 ;

]

4*3 — *4 =

4;

[

 

 

 

 

3)

 

3*з — 2*е ~

5;

 

 

 

 

 

 

 

 

 

4*» — *в = 7.

 

 

 

Р е ш е н и е .

Матрица

 

А коэффициентов

при

неизвестных во

всех трех системах одна

и та

же:

 

 

 

 

 

Ч < 4

Это тот случай, когда применение формулы (7,4) является выгод­ ным, т. е. выгодно решить все три системы путем отыскания обрат­ ной матрицы А~ 1 для матрицы А. Запишем все три системы в виде одного матричного уравнения

3

- 2 1

ж,

х3

х8]

Г2

3

5]

4

— 1|

Х-2

х8] — [1

4

7]

(отметим, что матрица, составленная из неизвестных, должна иметь тот же размер, что и матрица, составленная из свободных членов). Легко видеть, что для матрицы А обратная матрица

 

 

 

 

I

2 '

 

 

 

 

 

~ Т

Т

 

 

 

 

 

4

3

а потому по формуле

(7,4)

 

т Т-

 

1

2 "1

 

 

г*1

ха

3

5]

2

[х2

Х 4

хв

I ■[1

4

7\-

 

 

 

. т

т.1

 

 

Выполняя умножения в правой части равенства, получим

*1

хя

х*

 

х2

Х«

хв

 

откуда

X! =

0;

х2 = — 1 .

 

 

х3 =

1;

х4 = 0.

 

 

9

1

Решение этих трех систем уравнений по формуле (7,4) путем отыскания обратной матрицы потребовало значительно меньше труда, чем решение нх по известным читателю формулам Крамера.

Задача 7,4 (для самостоятельного решения). Решить системы уравнений:

1)

2Х] "4* Зх2 =

7;

1

2)

2х8 -)• Зх4 =

 

1:

1

 

Зхх — 2х2 =

1;

|

 

Зх3 — 2х« =

0; }

3)

2х* +

Зх, =

9; 4

4) 2х, + Зх8 =

1;

)

Зх* — 2х* =

3;

|

 

Зх7 2х8 =

3.

\

У к а з а н и е .

В матричной записи все четыре системы запишут­

ся так (проверьте, используя правило

умножения матриц и условие

их равенства):

 

Гх1

 

 

 

 

 

 

 

2

3]

х3

х$ х?

(7 1 9 1]

3

— 2

1*2

х4

х8

х8

[ 1 0

3

з ]‘

Д л я матрицы

 

 

 

 

а = \ : ;

 

2^

 

 

 

 

 

 

 

1

 

 

обратная

м атрица

 

 

 

 

 

 

"2

3 ‘

 

 

 

1 Г -2 —з]

 

А~* =*

13

13

 

[ - 3

 

2Г

3

2

 

 

 

~ 1 3

 

.13

13

 

 

 

 

 

17

 

 

 

О т в е т .

 

 

 

 

 

 

*2

19

 

 

 

2 .

 

 

Тз;

 

 

 

 

 

 

 

 

3 #

 

 

 

* 3 = = Т5;

 

 

* 4 — 13 ;

 

 

 

 

 

27

 

 

 

д:в ~

21

 

 

 

 

=

1 з ;

 

 

тз;

 

 

 

к

 

11

 

 

 

д8 =*

3

 

 

 

и

’|Ь

 

 

 

 

 

 

 

 

 

 

13

 

 

 

 

 

 

 

 

 

 

Задача

7,5.

Решить систему уравнений

 

 

 

 

2х + У

— 2г = 5; I

 

 

 

 

Ъх— 2у +

3г =

4;

 

 

 

 

— Зу + 5 г =

1,1

 

пользуясь

формулой

(7,4).

 

 

 

 

 

 

Р е ш е н и е .

Запишем

систему

в виде одного матричного урав­

нения

 

" 2

1

2 '

 

 

н

 

 

 

 

 

 

3

—2

 

 

 

 

 

 

2

—3

 

 

 

 

Здесь

матрица коэффициентов

 

 

 

 

А =

 

2

 

1

— 2

 

 

 

 

 

3

— 2

 

3

 

 

 

 

 

 

2

—3

 

5

 

а неизвестны е

х. у и

г найдутся

из

ф орм улы

 

 

 

X

 

= А- 1

'

5

(а)

 

 

 

У

 

 

4

 

 

 

 

г

 

 

 

 

1

 

 

 

 

 

для

чего применим

 

 

 

 

А =

'

2

1

1— 2

'

 

 

 

 

 

3

— 2 !

з

 

 

 

 

 

 

 

2

- 3

 

5 .

 

17В

Здесь

«11

Г

2

1 1

 

 

а12

- 2

1

 

«21 = |

 

2 ;

3 |

«22 *

 

 

1

3

- 2

}

 

 

 

3

] ’

 

 

 

 

 

а^ 1

находим

непосредственно

 

'2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

—1

 

1

- 2

-

1

1

 

У

 

 

 

 

 

 

 

 

 

 

а11

=

 

7

- 3

 

2

 

 

Л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

2)

р = а22 — «21ац «12 = 5 — [2;

—з] •

 

 

 

 

[1 ]-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 5 —

 

7

:

? ] * [

3

] “

5

( +

74) -

7 *

 

 

 

 

 

 

 

 

 

 

 

Итак,

р =

у1

;

Р- 1

=

7.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3)

7 =

Р_ 1 «21 = 712;

- 3

=

(14;

-21].

 

 

 

 

 

л

 

 

 

 

 

 

 

 

 

 

 

 

32 ]• |14; — 21]

I Т

 

 

4)

Яи1*1Д*и =

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.!

 

 

 

"

1’

 

 

 

 

‘ 2

 

1 "

_ [

- 2

 

3|

Гт

Г

 

 

 

 

7

(14;

- 2 1 ] .

Т

 

7

 

7

 

 

 

 

12

3

 

2

 

[—24

 

36]

3.

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ Т

 

 

 

 

 

 

У

~~Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

8 ‘

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

У

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60

 

96

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

"2

 

1"

 

Г 5

 

8

 

 

 

 

 

 

5)

-1

 

 

 

-I

7

У

+

 

Т “ У

 

= [ 1

- 1) .

 

 

«I, 1 +

<*п а12Тап1

=

3

 

2

 

 

60

 

96

 

[9

- 1 4 ] ’

 

 

 

 

 

 

 

 

 

7

7 .

 

 

т

у .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

"2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

■21)

 

 

 

7

=

- ( - 5 ;

8) = [5;

-

8);

6 )

Та и ’ =

— И 4 :

 

 

 

 

2

 

 

 

 

 

 

" 2

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

— 2

 

 

 

 

 

 

 

 

7) — «11«12р~1= —

 

У

 

У

 

 

у

 

 

 

 

У

и

 

-

 

3

 

2

 

1

з

 

 

 

 

 

 

12

з

 

 

 

 

 

 

 

У “

У

 

 

 

 

 

 

 

~ У

Теперь все элементы обратной матрицы известны

и по фор*

муле

(5,7)

 

 

 

 

 

"

1

-

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

А '1 -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

— 14

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

-

8

 

 

7

 

 

 

 

 

 

 

 

Заметим, что вычисление по формуле (5 ,2) оказалось бы менее громоздким. Следовало бы вычислить 9 определителей второго,

один

определитель

третьего

порядка

и

образовать союзную ма­

трицу

А.

 

А~ 1

 

 

 

 

 

Подставляя матрицу

в формулу

(а),

получим

 

X

'

1

— 1

1

*

'

5 '

 

У

=

9

— 14

12

V

4

 

2

 

5

—8

7

 

 

1

откуда

 

X

'

 

 

 

 

 

 

2

 

 

У = 1

г0

т. е. х = 2; у = 1 ; 2 = 0.

Задача 7,6. Решить системы уравнений

 

 

 

1) 4*, +

 

5х3=

7;

)

2)

4* 4 + 5*в =

1 ;

 

х2 — 6*, = 1 1 ; |

 

*» — 6*. = 2 ;

 

3*1

 

3 =

— 2; ]

 

3*4 -+• 4*4 = 11;

 

 

 

 

 

3)

4* 7 +

5*в =

0; |

 

 

 

 

 

 

 

 

* 86*4 = 5;/

 

 

 

 

 

 

 

 

3*, +

4*9 =

1 .)

 

 

 

Р е ш е н и е .

Запишем все три системы

в виде

одного матрич­

ного уравнения

 

 

 

 

 

 

 

 

 

4

0

 

5 '

X,

*4 * ; |

Г 7

1

0

0 >1 —6 *2

* 8 I = I 11

2

5

со

О

 

 

 

_х3 х, * , |

| — 2 11

1 .

Из этого следует, что матрица неизвестных

 

 

 

*1

*4

*7

 

7

1

0

 

 

 

* 2

 

 

Хй

 

11

2

5

 

(а )

 

, * 3

 

Х в

* 0 _

— 2 11

1

 

 

где А 1 — обратная матрица для матрицы коэффициентов

4 0 5

0 1 — 6

3 0 4

Найдем для этой матрицы обратную. В данном случае, учи­ тывая характер этой матрицы, удобнее воспользоваться общей формулой (5,2) для определения обратной матрицы. Определитель матрицы А

|Л|=1.

Соседние файлы в папке книги