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

0-LALect2

.pdf
Скачиваний:
15
Добавлен:
10.02.2015
Размер:
382.92 Кб
Скачать

что из 6= 0 следует существование обратной матрицы X 1, которая является матрицей перехода от базиса (xi) к базису (ei).

Как изменяется матрица оператора при переходе к другому базису.

Замечание (существенное). Теперь мы можем понять, почему матрицы умножаются по правилу (строка на столбец). Это умножение соответствует умножению линейных операторов. Сначала мы применяем оператор A с матрицей (aij)i n:j n. Применяя A к вектору ei, мы получа-

ем A(ei) = aijej. Аналогично, применяя B к вектору ej, мы получаем

B(ej) =

Pk bjjkek. Что же будет, если мы сначала применим A к ei, а

потом

B к A(ei)?

 

 

 

 

 

 

 

P

 

aijej! = Xj

 

 

 

 

 

B(A(ei)) = B

Xj

aijB(ej) = Xj

aij Xk

bjkek = Xk

Xj

aijbjk

!

ek:

Рассмотрим линейный оператор A, матрицу которого в базисе (ei) мы также обозначим A:

0

 

 

 

 

 

 

1

 

a11

a12

a1n

C

 

a

 

a

 

 

a

 

 

A = B a

21

a

22

 

a

2n

:

B

n1

 

n2

 

 

nn

C

 

@

 

 

 

 

 

A

 

(напомним, что A(ei) = (ai1; ai2; ; ain) для всех i n). Теперь мы хотим вычислить матрицу этого же оператора в базисе (xi) (эту матрицу мы обозначим B), то есть

0

 

 

 

 

 

 

1

 

b11

b12

b1n

C

 

b

 

b

 

 

b

 

 

B = B b

21

b

22

 

b

2n

:

B

n1

 

n2

 

 

nn

C

 

@

 

 

 

 

 

A

 

(где A(xi) = (bi1; bi2; ; bin) для всех i n). Здесь доказательства не будет, только объяснение: сначала мы перейдем от базиса (xi) к базису (ei) с помощью матрицы X 1, потом применим оператор A с помощью матрицы A, после этого обратно вернемся к базису (xi) с помощью матрицы X. В результате мы получаем представление

B = XAX 1:

51

Все таки, для желающих во всем разобраться я написал доказательство (но учить это не обязательно). Введем обозначение для матрицы X 1:

X 1 =

0 y21

y22

 

y2n

1:

 

y11

y12

 

y1n

 

 

B y

y

 

y

C

 

B n1

n2

 

nn

C

 

@

 

 

A

По формуле умножения матриц нам нужно получить следующее представление для произвольного элемента матрицы B:

n

n

 

Xk

X

( )

bij =

xikaklylj:

=1 l=1

 

 

 

n

Сначала мы разложим вектор xi

по базису (ei): xi =

kP

xikek. Теперь

 

 

=1

применим к этой линейной комбинации оператор A. Так как линейный оператор, примененный к линейной комбинации, является линейной комбинацией, мы имеем:

n

n n

X

XX

A(xi) = xikA(ek) =

xikaklel:

k=1

k=1 l=1

n

P

А теперь мы разложим каждый el по базису (xj): el = yljxj: Подста-

j=1

вим это разложение в предыдущую формулу и получим:

n

n

n

XXl

X

A(xi) =

 

xikaklyljxj:

k=1

=1 j=1

Таким образом, коэффициент при элементе xj в разложении A(xi) по

базису (xj) равен bij из равенства (**).

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

Пример. Если вы сидите в среднем ряду, то понятия ’левый ряд’ и ’правый ряд’ у вас и преподавателя противоположны. У вас с преподавателем разные системы координат, но ряды по бокам не остается прежними.

52

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

Собственные векторы и собственные значения

Во многих задачах математического анализа, физики, механики, геометрии (применения к геометрии мы скоро рассмотрим), теории вероятностей, математической статистики и в других областях возникает важная и естественная задача: для каких векторов x и чисел имеет место следующее соотношение:

A(x) = x:

( )

Одно тривиальное решение мы знаем если x = 0, то равенство (***) справедливо для всех . Но интерес представляют ненулевые решения этого уравнения. Число в этом уравнении называется собственным значением оператора A, вектор x собственным вектором оператора

A. Ввиду линейности оператора A собственный вектор определяется с точностью до множителя если для x равенство (***) справедливо, то оно справедливо и для x, умноженного на любое число, например, на 2, 3, 4, 5 и т. д. Напомним, что единичный линейный оператор E имеет свойство: E(x) = x. Поэтому соотношение (***) можно переписать в следующем виде:

(A E)(x) = 0:

( )

Перепишем это равенство в виде уравнения для координат вектора x, для этого сначала запишем матрицу оператора A E:

A E =

0 a21

a22

 

 

a2n

1

 

0

0

0

1 =

 

 

a11

a12

 

 

a1n

 

 

 

 

0

 

0

 

 

B a

a

a

C

 

B

0

0

 

C

 

B

n1

n2

 

 

nn

C

 

B

 

 

 

 

C

 

@

 

 

 

 

A

 

@

 

1:

 

A

 

= 0 a21

 

a22

a2n

 

 

 

 

 

 

 

a11

a12

 

 

 

a1n

 

 

 

 

 

 

 

B a

 

a

 

 

 

a

 

 

C

 

 

 

 

B

n1

 

 

n2

 

 

nn

 

 

C

 

 

 

 

@

 

 

 

 

 

 

 

A

 

 

Для координат вектора x мы получаем следующую систему уравнений:

8 a21x1 + (a22

 

)x2

+

 

+ a2nx3

= 0

(a11

)x1

+ a12x2

+ + a1nxn = 0

>

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

<

 

 

 

 

 

 

 

 

>

>

: ann1x1 + an2x2 + + (ann )xn = 0

53

Мы уже знаем, что ненулевое решение этой системы существует лишь если определитель системы равен нулю:

 

a21

a22

 

 

a2n

 

 

= 0:

 

a11

a12

 

 

a1n

 

 

 

 

 

 

 

 

 

 

 

 

 

a

a

 

 

a

 

 

 

 

n1

n2

 

 

nn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получающееся уравнение является уравнением n-го порядка и называется характеристическим уравнением для A. Некоторые факты из теории уравнений n-го порядка будут изложены позже.

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

Доказательство несложно. Напомним, что в новом базисе матрица оператора будет иметь вид B = XAX 1, где X матрица перехода к новому базису. Напомним, что XX 1 = E, а так как для всех n-мерных матриц T справедливо T E = ET = E, то E = XEX 1. Мы будем также использовать свойство определителя произведения матриц: если S и Tматрицы n-го порядка, а определитель матрицы T мы обозначаем jT j, то

jT Sj = jT j jSj;

поэтому jXj jX 1j = jEj = 1. Итак,

jB Ej = jXAX 1 XEX 1j = jX(A E)X 1j = jXj jA Ej jX 1j =

jXj jX 1j jA Ej = jA Ej:

Рассмотрим конкретный пример из геометрии: линейный оператор A : R2 ! R2 задается матрицей

A =

5 2

2 1

Тогда мы получаем уравнение:

5 2

= (5 )(1 ) 4 = 0;

2 1

54

p p

2 6 + 1 = 0; 1;2 = 3 9 1 = 3 2 2:

Для каждого из двух собственных значений мы вычисляем собственный

p

 

 

p

 

 

вектор. Пусть 1 = 3 + 2 p2. Тогда из двух уравнений [5 (3 + 2

 

2)]x1 +

2x2 = 0, 2x1 + [1 (3 + 2 2)]x2 = 0 мы можем взять лишь одно.

Таким образом, собственный вектор, соответствующий собственному значению 1, имеет вид:

 

 

 

 

1

 

x1 = C

1

;

p

; где C произвольное число, отличное от 0:

2

 

 

 

 

 

 

2

2 2

Аналогично находится собственный вектор x2.

Мы увидим в дальнейшем, что эти векторы определяют направления главных осей эллипса, задаваемого уравнением 5x2 + 4xy + y2 = 1. Итак, мы уже подошли к квадратичным формам, с помощью которых в геометрии задаются уравнения различные фигуры и тела второго порядка на плокости и в пространстве, частными случаями которых являются окружности и сферы.

Квадратичные формы.

Определение. Квадратичной формой от n переменных называется

функция

n n

XX

A(x1; :::; xn) = aijxixj:

i=1 j=1

В этой сумме при i 6= j коэффициент при xixj имеет вид aij + aji, мы можем менять эти два числа, не меняя их сумму, при этом сама квадратичная форма не изменится. Поэтому мы будем всегд считать, что aij = aji. Например, (x + y + z)2 мы записываем как x2 + xy + xz + yx + y2 +yz +zx+zy +z2. Квадратичная форма определяется своей матрицей

0 a21

a22

 

a2n

1

a11

a12

 

a1n

 

B a

a

 

a

C

B n1

n2

 

nn

C

@

 

 

A

(отметим еще раз, что в этой матрице всегда aij = aji).

Квадратичную форму мы можем считать и будем считать определенной на пространстве Rn. Каждому вектору (x1; :::; xn) мы сопоставляем число A(x1; :::; xn). Если в пространстве Rn введен новый базис, координаты в новом базисе мы обозначаем u1; :::; un, то в новом базисе матрица

55

квадратичной формы меняется и она записывается уже в виде

n n

XX

A(x1; :::; xn) = bijuiuj;

i=1 j=1

где B = (bij) матрица квадратичной формы в новом базисе. Переход к новому базису это замена переменных. Процедуру такой замены мы рассмотрим на примерах.

Заметим также (это важно для математической статистики), что квадратичную форму можно записать как произведение матриц:

 

 

x2 + x2 + ::: + x2

+ x2

= (x1

; x2; :::; xn)

 

0 x2

1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

C

 

 

 

 

 

 

 

1 2

n 1

 

n

 

 

 

 

 

 

 

 

B x

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

 

A

 

 

 

 

 

n

n

aijxixj = (x1; x2; :::; xn)

 

0 a21

a22

 

 

 

a2n

1 0 x2

1

:

Xi

X

 

 

 

 

 

a11

a12

 

 

 

a1n

 

 

 

 

 

x1

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B a

 

a

 

 

 

 

a

 

C B x

 

 

 

 

 

 

 

B

 

n1

 

n2

 

 

 

 

nn

C

 

B

 

n

C

 

=1 j=1

 

 

 

 

@

 

 

 

 

 

 

 

 

 

A @

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Важнейшим примером квадратичной формы является

 

 

 

 

 

 

 

 

 

A(x1; :::; xn) = 1x12 + 2x22 + ::: + n 1xn2

1 + nxn2

 

 

 

 

 

( )

Матрица такой квадратичной формы является диагональной

 

 

 

 

 

 

0 01

2

 

 

0

 

1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

0

 

 

 

 

0

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

n

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

Определение. Если квадратичная форма задана в виде (*****), мы будем говорить, что она задана в каноническом виде. Число ненулевыхi в (*****) называется рангом квадратичной формы. Разность числа положительных и отрицательных i в (*****) называется сигнатурой квадратичной формы.

Центральной задачей теории квадратичных форм является приведение квадратичных форм к каноническому виду с помощью замены переменных xi на новые переменные xi. Кстати, эта техника полезна для

56

вычисления кратных интегралов. Задача состоит в том, чтобы найти новые переменные ui = vi1x1 + ::: + vinxn такие, что имеет место равенство

nn

XX

aijxixj = 1u21 + ::: + nu2n;

i=1 j=1

которое превращается в тождество после подстановки вместо ui их представлеий ui = vi1x1 + ::: + vinxn.

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

Пример. A(x; y; z) = x2 + y2 + z2 + xy + xz 4yz. Приводить к каноническому виду мы начнем с обработки всех слагаемых, содержащих множитель x, при этом все остальные слагаемые также изменятся. Итак,

x2 + xy + xz =

x + 2y +

2z

2

4y2

 

4z2

 

2yz:

 

 

1

 

1

 

1

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вводя новую переменную u = x + 12 y + 12 z, мы получаем

A = u2 + 34y2 + 34z2 92yz:

Теперь мы займемся y:

34y2 92yz = 34(y2 6yz) = 34[(y 3z)2 9z2]:

Итак, введя новую переменную v = y 3z, мы получаем

A = u2 + 34v2 6z2:

Мы привели форму к каноническому виду, причем с двумя положительными коэффициентами и одним отрицательным.

Теперь попробуем сделать то же самое, но начиная со слагаемых, содержащих множитель z:

z2 + xz 4yz =

z + 2x 2y

2

4x2

4y2 + 2xy:

 

 

1

 

1

 

 

 

 

57

 

 

 

 

 

Введем новую переменную u = z + 12 x 2y. Получаем:

A = u2 + 34x2 3y2 + 3xy:

Далее освобождаемся от x:

34x2 + 3xy = 34(x2 + 4xy) = 34[(x + 2y)2 4y2]:

Обозначая v = x + 2y, получаем

A = u2 + 34v2 6y2:

Опять два положительных знака, один отрицательный.

Как приводить к каноническому виду форму без квадратов, например,

A(x; y; z) = xy + xz + yz:

Очень просто, так как коэффициент при xy не равен нулю, мы переходим к новым переменным x = u + v, y = u v. Наша форма принимает вид

A = u2 v2 + 2uz:

Введя новую переменную w = u + z, мы приводим форму к виду

w2 v2 z2:

Ранг равен 3, сигнатура равна 1.

Справедлива следующая теорема, которую доказывать мы не будем:

Теорема. Любая квадратичная форма приводится заменой переменных к каноническому виду. Ранг и сигнатура квадратичной формы в каноническом виде не зависят от способа приведения.

Определение. Рангом и сигнатурой квадратичной формы называются ее ранг и сигнатура после приведения к каноническому виду.

Важное значение имеют определенные квадратичные формы. Определение. Квадратичная форма называется положительно

определенной, если ее значения больше нуля для любого ненулевого вектора (x1; x2; :::; xn). Квадратичная форма называется отрицательно определенной, если ее значения меньше нуля для любого ненулевого вектора (x1; x2; :::; xn). После приведения квадратичной формы к каноническому виду становится ясно, что сигнатура положительно определенной

58

квадратичной формы равна (+n), а отрицательно определенной квадратичной формы ( n).

Теорема. (Критерий Сильвестра положительности квадратичной формы, без доказательства). Квадратичная форма положительно определена тогда и только тогда, когда

a > 0;

 

a11

a12

 

> 0;

a21

a22

a23

 

> 0; :::;

a21

a22

 

a2n

> 0:

11

 

a

 

a

 

 

 

a11

a12

a13

 

 

a11

a12

 

a1n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

22

a

 

a

 

a

 

a

a

 

a

 

 

 

 

 

 

 

 

 

31

 

32

 

33

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

n2

 

nn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что первое неравенство относится к минору первого порядка. Итак, n миноров больше нуля, где i-й минор расположен на пересечении первых i строк и первых i столбцов. Эти миноры называются главными.

Отрицательно определенная квадратичная форма получается из положительно определенной изменением знака, т.е. каждый коэффициент aij заменяется на aij. i-й главный минор для отрицательно определенной квадратичной формы получается из главного минора для положительно квадратичной формы изменением знака всех i строк (или, что то же) изменением знака всех i столбцов. Таким образом, знак определителя тоже меняется i раз и для четного i остается прежним (больше нуля), а для нечетного i меняется (меньше нуля). Итак, для отрицательной определенности квадратичной формы необходимо и достаточно, чтобы знаки главных миноров чередовались отрицательный, положительный, отрицательный, положительный и т.д.

Пример. A(x1; x2x3) = x21 + 12 x1x2 + 12 x1x3 + 12 x2x1 +x22 + 12 x2x3 + 12 x3x1 +

12 x3x2 + x23. Главные миноры имеют вид:

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

1

 

 

 

1

2

2

 

 

 

2

1

 

1

1

1

 

1;

 

11

2

 

;

 

2

2

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проверьте, что все они больше нуля.

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

59

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

Элементы аналитической геометрии на прямой, плоскости и в трехмерном пространстве

Евклидово пространство и квадратичные формы Комплексные числа и многочлены Линейное программирование

60

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]