Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математика / Шилкин Лекции.PDF
Скачиваний:
31
Добавлен:
11.05.2015
Размер:
7.25 Mб
Скачать

ЛЕКЦИЯ 3

Метод Гаусса. Решение системы m линейных уравнений с n неизвестными (метод последовательного исключения неизвестных).

Элементарными преобразованиями системы (9) будем называть:

1.Умножение какой-либо строки системы (9) на число α 0

2.Прибавление к какой-либо строке системы (9) другой строки, умноженной на какое-то число α

3.Замена местами двух строк

Справедлива следующая теорема.

Если мы приходим от данной системы к другой системе с помощью элементарных преобразований, то новая система эквивалентна данной системе (эквивалентна – значит, все решения одной системы являются решениями другой и наоборот).

Перейдем теперь непосредственно к алгоритму.

Пусть дана система (9) и пусть среди коэффициентов a11, a21,K, am1 при неизвестном x1 есть отличные от нуля, например, a11 0 (если a11 = 0 , выбираем коэффициент ai1 0 и ставим i -ю строку первой, перенумеровав коэффициенты, получим a11 0).

Возьмем первое уравнение, до множим его почленно на ai1 и прибавим к

a11

его i -му уравнению для каждого i = 2,K, m . В результате получим систему

a11x1

+ a12 x2

+K+ a1n xn = b1

 

 

~

 

 

 

 

 

~

 

 

 

 

~

 

 

 

 

a22 x2

+K+ a2n xn

= b2

(14)

вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

x

 

 

 

 

~

 

x

 

~

 

 

 

a

m2

2

+K+ a

mn

n

= b

 

 

 

 

 

 

 

 

 

 

 

m

 

~

 

 

aij ai1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где aij = aij

 

 

 

 

 

, i = 2, m,

j = 2, n .

 

 

a11

 

 

 

 

 

 

 

 

 

исключаем x1 из всех уравнений, кроме первого.

Таким

образом,

 

мы

 

Система (14) эквивалентна системе (9), но более проста по виду. Может получиться следующее:

1.Какая-то строка системы целиком состоит из нулей, тогда ее просто вычеркиваем.

2.В какой-то строке слева от равенства стоят нули, а справа – число, отличное от нуля, тогда система (9) несовместна.

~= ~ =

3.Пропали все неизвестные xj , т.е. a2 j 0,K, amj 0 . В этом случае

система не изменяется, ведь xj входит в первое уравнение.

Далее поступаем аналогичным образом с x2 . После конечного числа таких

операций (не более

(m 1)m ) может получиться следующая система:

 

2

a11x1 + a12 x2

+K+ a1k xk +K+ a1n xn = b1

 

 

 

 

 

a22 x2

+K+ a2k xk

+K+ a2n xn

= b2

(15)

LLLLLLLLLLLL

 

 

 

 

 

 

 

′ ′

 

 

 

 

akk xk +K+ akn xn = bk

 

 

где a11 a22′ K akk 0 .

На этом завершен прямой ход метода Гаусса. Отметим так же, что k m, k n . В этом случае система (9) совместна. Она будет определенной

при k = n и неопределенной при

k < n. Действительно, при k = n система

a11x1 + a12 x2

+K+ a1n xn = b1

 

 

 

(15) имеет вид

a22 x2 +K+ a2n xn

= b2

(16)

LLLLLLLL

 

 

 

 

 

 

ann xn = bn

 

 

xn =

bn

,

 

Т.к. a11 a22

K ann 0, то из последнего уравнения находим

 

 

 

ann

 

подставляя xn в предпоследнее, находим xn1 и т.д., пока не дойдем до первого уравнения и не определим x1 . Это называется обратным ходом метода Гаусса.

Если k < n, то для “свободных” переменных xn+1,K, xn выберем

произвольные числовые значения (на практике обычно обозначают их буквами и получают решение в виде зависимости от параметров), после чего, перенося их в правые части уравнений, получим систему вида (16), у которой уже будет k уравнений, и аналогично отыщем xk , xk 1,K, x1 . Т.к. значения для

свободных неизвестных xk +1,K, xk можно выбрать бесконечным числом

различных способов, то система (9) будет совместной, но неопределенной. Итак, метод Гаусса применим к любой системе линейных уравнений. При

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

2x

+ x

 

3x + x

 

 

= 0

 

23

 

2

 

12

 

8

T

 

 

1

 

2

1

4

 

 

 

 

 

 

 

 

x2 2x4 = 2

 

 

 

x =

 

 

,

 

,

 

 

,

 

 

 

 

 

 

 

7

 

7

7

7

Пример 1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

x2 2x3 + x4

 

= −1

x <

(33,14,20,8)T

x

+ 2x

2

+ x + 3x

4

=1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 5x

8x + x

=

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

3

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3x1 + x2

3x3 5x4 =1

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2:

x

 

7x + 2x

 

= −5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11x

+ 20x 9x

= 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2x1 +

3x2 x3 + x4 =1

 

 

 

1

 

c1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 3 x =

 

 

c1

 

c2 +

 

 

 

 

 

 

Пример 3:

 

 

 

 

 

 

 

 

 

 

10

2

5

c ,c

 

R

8x1 +12x2 9x3 +8x4

 

 

 

 

 

 

 

2

 

 

 

 

 

4x1 +6x2 +3x3 2x4 = 3

 

 

 

 

 

1 4

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2x

 

+3x +9x

7x

 

= 3

 

 

 

 

 

5

+

5 c1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

3

 

 

4

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисление обратной матрицы методом Гаусса

 

 

 

Пусть

 

матрица

A

 

невырожденная.

 

 

Составим

 

 

матрицу

 

 

 

a

 

a

 

L a

 

1

0

 

L 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

12

 

 

1n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[A

 

E]= a21

 

a22 L a2n

 

0

1

 

L 0

 

,

 

которая

получается,

если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

L L L

L L L L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an2 L ann

 

0 0 L 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

справа к

матрице

 

A

приписать

единичную

матрицу

E .

 

Применяя

метод

Гаусса,

матрицу

[A

 

E]

 

 

 

можно

 

 

 

преобразовать

к

 

 

виду

 

 

 

 

 

 

 

 

 

1 0 L 0

 

b11

b12

L b1n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

L

0

 

 

b

b

L b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

22

 

 

2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L L L L

 

L L L L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

bn1

bn2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0 L

 

L bnn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(элементарные преобразования делаем только над строками матрицы). Отсюда

b

b

11

12

B = b21

b22

L L

 

bn2

bn1

Пример 4:

L b

 

 

 

 

 

 

 

 

 

 

1n

 

 

 

 

 

 

 

 

 

 

L b2n

= A1 , т.к. A1[A

 

E]= [A1 A

 

A1 E]= [E

 

A1 ]

 

 

 

 

L L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L bnn

1

1

0

 

 

 

 

 

 

 

 

Найти A1 для A = 2

1

 

2 .

 

 

 

 

 

 

 

 

 

 

3

0

 

 

 

0

 

 

 

1

1

0

 

1

0

0

 

 

 

 

 

[A

 

E]= 2

1

2

 

0

1

0

. После элементарных преобразований

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

3

0

 

0

0

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

1 0 0

 

1 0

3

 

 

1 0

3

 

 

 

0

0

1

 

1

 

0

1

 

получаем 0 1 0

 

3

. Отсюда A

 

= 0

3

.

 

 

 

 

1

 

 

 

 

1

 

 

 

1

1

 

 

 

1

 

0 0 1

 

2

2

 

 

1

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Векторы и линейные операции над ними

В математике различают скалярные и векторные величины. Скалярные величины полностью определяются своим значением, например, температура тела, масса, объем и т.д. Но для силы, скорости, ускорения, кроме величины надо указать еще и направление в пространстве. Такие величины называют векторными.

Вектором называется направленный отрезок AB с начальной точкой A и конечной точкой B , обозначается AB (или a,b ,K, x ). Будем считать, что два вектора AB и CD равны, если они одинаково направлены и длины отрезков AB и CD равны AB = CD .

Мы будем рассматривать свободные векторы, т.е такие векторы, которые можно перемещать параллельно самим себе. Далее рассматриваем только свободные векторы.

Нулевым называется вектор, начало, и конец которого совпадают (обозначается o ).

Длина нулевого вектора равна нулю, и ему можно принимать любое направление.

Длиной или нормой вектора AB называется длина отрезка AB и

обозначается AB или AB .

Вектор a , длина которого равна 1, называется единичным или ортом. Векторы, параллельные одной прямой, называются коллинеарными.

Векторы, параллельные одной плоскости, называются компланарными. Суммой двух векторов a и b называется вектор a +b , являющийся

диагональю параллелограмма, построенного на векторах a и b как на сторонах:

 

 

a +

 

b

b

a

Суммой конечного числа n векторов a1, a2 ,K, an , называется замыкающий

вектор OAn :

a2

 

 

a4

a1

a3

O

 

 

An

Вектор

(

 

) называется противоположным вектору

 

, если (

 

)

 

 

,

b

b

b

 

b

b = b и его направление противоположно вектору b .

Разностью двух векторов a и b называется вектор a b , являющийся суммой векторов a и (b )

a

 

a +

 

b

b

b

a

Произведением вектора a на действительное число α называется вектор

 

 

= αa , длина которого

 

 

 

=

 

αa

 

=

 

α

 

 

 

a

 

,

а направление совпадает с a , если

 

b

 

b

 

 

 

 

 

 

 

λ > 0

и противоположно

 

a , если α < 0.

Если a 0 , то единичный вектор

ao =

 

1

 

ar, называется ортом вектора a .

 

 

a

 

 

 

 

 

 

 

 

 

 

 

Введенные операции сложения векторов и умножения на число называются линейными и обладают свойствами коммутативности, ассоциативности и дистрибутивности:

1.

a +

 

 

=

 

 

+ a a,

 

 

 

 

 

 

b

b

b

2.

a +(

 

 

+c) = (a +

 

 

) +c a,

 

,c

b

b

b

3.

a + o = a a

4.

(a +(a)) = o a

5.

α(βa) = (αβ)a α, β R, a

6.

1 a = a a

7.

(α + β)a = αa + βa α, β R, a

8.

α(a +

 

) = αa +α

 

α R, a,

 

 

b

b

b

Множество векторов пространства, удовлетворяющих свойствам 1-8, называется линейным или векторным пространством (ЛВП).

Заметим, что в определение ЛВП не входит никакое умножение векторов. Отметим также, что с математической точки зрения вектором можно назвать любой объект, а не только направленный отрезок. Лишь бы были определены

операции их сложения и умножения на число. Примером ЛВП может быть:

1.множество обычных векторов в R2 или R3

2.множество матриц размерами m×n с введенными ранее операциями сложения и умножения на число

3.множество функций, непрерывных на отрезке [a,b] и т.д.

Линейной комбинацией n векторов e1, e2 ,K, en , называется вектор

x = α1e1 +α2e2 +K+αnen , где α1,α2 ,K,αn R называются коэффициентами комбинации.

Если вектор x представлен в виде линейной комбинации, то говорят, что он разложен по векторам e1, e2 ,K, en .

Векторы e1, e2 ,K, en называются линейно-зависимыми, если существуют такие постоянные числа λ1, λ2 ,K, λn , не все одновременно равные нулю, что

выполняется равенство λ1e1 + λ2e2 +K+ λnen = 0

(1)

Если же равенство (1) выполняется только при всех λ1 = λ2 =K= λn = 0 , то векторы e1, e2 ,K, en называются линейно-независимыми.

Примером линейно-зависимых векторов могут быть любые два вектора a и

b , т.к. в этом случае имеет место соотношение b =αa , где α R .

Из элементарной математики известно, что если e1 и e2 два неколлинеарных

вектора, то всякий компланарный им вектор a однозначным образом линейно выражается через векторы e1 и e2 : a = α1e1 +α2e2 .

Отсюда следует, что любые три компланарных вектора всегда линейнозависимы. Аналогично, если e1, e2 , e3 - три некомпланарных вектора, то любой

вектор a пространства единственным образом разлагается в их линейную комбинацию: a = α1e1 +α2e2 +α3e3 , т.е. любые четыре вектора в пространстве

линейно-зависимы.

Базисом пространства называется совокупность линейно-независимых векторов, по которым можно разложить любой вектор этого пространства. Векторы, составляющие базис называются базисными.

Базис на прямой образует любой ненулевой вектор, на плоскости – любая пара неколлинеарных векторов плоскости, в пространстве – любая тройка некомпланарных векторов трехмерного пространства.

Координатами вектора a

в базисе e1, e2 , e3 называются коэффициенты

линейной комбинации базисных векторов, представляющих вектор a

a = α1e1 +α2e2 +α3e3

(2)

Основное значение базиса состоит в том, что линейные операции над векторами в заданном базисе сводятся к обычным линейным операциям над числами-координатами векторов. Скажем, при сложении двух векторов их

соответствующие координаты суммируются, при умножении вектора на число все его координаты умножаются на это число.

Среди множества систем координат в пространстве R3 особенно удобна декартовая совокупность трех взятых перпендикулярно осей x, y, z

называемых осями координат и точки O - начала координат. Единичные векторы i , j, k , направленные вдоль соответствующих осей x, y, z образуют

ортонормированный прямоугольный базис. z M (x, y, z)

 

 

 

 

 

 

 

 

 

k

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

j

x

Известно, что каждый вектор r можно однозначным образом разложить по

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

, k

 

 

 

 

 

 

 

 

 

 

 

векторам i

j

 

 

 

 

 

 

 

 

 

 

 

rr = xi + yj + zk

(3)

Числа x, y, z называются координатами вектора

 

 

 

 

 

 

 

 

 

 

 

 

r

в базисе i

,

j

, k

, и

записывается, что

r

= (x, y, z) , что равносильно (3).