Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Alg i metodi vichisl / Teorija / course_359.ppt
Скачиваний:
75
Добавлен:
23.02.2016
Размер:
695.3 Кб
Скачать

Постановка задачи.

Пусть дана система линейных алгебраических уравнений

а11 x1

+ a12 x2 + … + a1n xn = b1

а21 x1

+ a22 x2 + … + a2n xn = b2

……………………………………

аn1 x1

+ an2 x2 + … + ann xn = bn

Используя понятия матрицы:

a11

A ...

an1

и вектор столбца:

b1 B ...

bn

a12

...

an2

è

...

...

...

X

a1n

...

ann

x1

...

xn

систему уравнений можно переписать в виде:

AX = B.

Правило Крамера состоит в вычислении главного определителя матрицы А:

 

a11

a12

a13

...

a1n

D

...

 

... ... ...

 

an1

an2

an3

...

ann

и побочных определителей:

 

b1

a12

...

a1n

 

a11

b1

a13

...

a1n

 

D1

...

...

...

...

, D2

...

 

...

...

...

,...

 

bn

an2

...

ann

 

an1

bn

an3

...

ann

 

И тогда вычисление корней:

x1 = D1 /D, x2 = D2 /D, … xn = Dn /D.

Число требуемых для вычисления определителей операций:

N n(n! + 1).

2. Метод обратной матрицы

Отталкиваясь от АХ = В, умножим обе части на обратную матрицу A-1 :

A-1 А Х = A-1 В

или

Х = A-1 В.

Таким образом, решение системы сводится к вычислению A-1, но её вычисление также трудоемко.

3. Метод Гаусса

Рассмотрим на примере системы трех уравнений:

а11 x1 + a12 x2 + a13 x3 = b1 а21 x1 + a22 x2 + a23 x3 = b2 а31 x1 + a32 x2 + a33 x3 = b3

Прямой ход метода Гаусса состоит в следующем:

Умножим первое уравнение на – a21 /a11 и сложим со вторым:

a21x1 a21a12 / a11x2 a21a13 / a11x3 a21b1 / a11 a21x1 a22 x2 a23x3 b2

_______________________________________

(a22 a21a12 / a11)x2 (a23 a21a13 / a11)x3 b2 a21b1 / a11

Если обозначить первую скобку за a22 , а вторую за a’23 и полученную правую часть за b2 , то получим:

a11 x1 a12 x2 a13 x3 b1

a x

2

a x

b

22

23

3

2

Аналогично: умножим первое уравнение на –a31/a11 , сложим с третьим уравнением

и обозначим полученное уравнение так:

a‘32 x2 + a‘33 x3 = b‘3 .

Обобщенная формула для «штриховых» коэффициентов.

a'

a

ai1a1 j

;

b' b ai1b1

,

i, j 2,3

 

ij

ij

a11

i i

a11

 

 

 

 

 

 

 

Исключим теперь x2 из 3-го уравнения. Таким образом, система сведется к виду:

a11 x1 a12 x2 a13 x3 b1

 

 

 

 

a22 x2

a23 x3

b2

a x

b

 

33

3

3

 

Для «дважды штриховых» коэффициентов общие формулы имеют вид:

a''

a'

 

a'i2a'2 j

 

b'' b'

 

a' b'

 

 

 

;

 

i2 2

,

i, j 2,3

 

 

ij

ij

 

a'22

 

i

i

 

a'22

 

 

 

 

 

 

 

 

 

Обратный ход метода Гаусса состоит в следующем:

Находим x3 из 3-го уравнения:

x3 = b'‘3 /a'‘33

Потом находим

x2 = (b‘2 – a‘23 x3 )/a‘22 ,

x1 = (b1 – a12 x2 – a13 x3 )/a11 .

Аналогично строится процесс метода Гаусса для системы из N уравнений. Обобщим формулы «штриховых» коэффициентов.

a(k) a(k 1)

 

a(k 1)ika(k 1)kj ;

b(k) b(k 1)

 

a(k 1)ikb(k 1)k ,

ij

ij

 

 

a(k 1)kk

 

i

i

a(k 1)kk

 

 

 

 

 

 

i 2,3,..., n,

j i 1,i,..., n 1, k

i 1,i 2,..., n.

Пример. Методом Гаусса решить систему:

 

10x1 7x2 7

 

 

 

 

 

 

3x1 2x2 6x3 4

 

5x

x

5x

6

 

1

2

3

 

4. Итерационный метод Гаусса—Зейделя

Рассмотрим метод Гаусса—Зейделя (на примере трех уравнений)

a11 x1a21 x1a31x1

a12 x2 a13 x3 b1

a22 x2 a23 x3 b2

a32 x2 a33 x3 b3

Выразим x1, x2 и x3 соответственно из 1- го, 2-го и 3-го уравнений:

x1 1 a11 (b1 a12 x2 a13 x3 )

 

 

1/ a22 (b2 a21x1 a23 x3 )

x2

x 1/ a (b a x a x

2

)

 

3

33

3

31

1

32

 

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