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

книги из ГПНТБ / Дроздов Н.Д. Линейная алгебра в теории уравнивания измерений

.pdf
Скачиваний:
2
Добавлен:
24.10.2023
Размер:
9.33 Mб
Скачать

н . д . ДРОЗДОВ

Л И Н Е Й Н А Я

А Л Г Е Б Р А

ВТ Е О Р И И

У Р А В Н И В А Н И Я

И З М Е Р Е Н И Й

ИЗДАТЕЛЬСТВО «НЕДРА»

Мо с к в а - 1 9 7 3

У Д К 512.8

Д р о з д о в Н. Д . Линейная алгебра в теории уравнивания измерений. М., «Недра», 1973. 216 с.

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

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

Иллюстраций 11. Список литературы — 44 названия.

© и з д а т е л ь с т в о . . Н Е Д Р А " , 1 9 7 3

ПРЕДИСЛОВИЕ

Теория математической обработки результатов измерений, пройдя большой путь развития, продолжает и в настоящее время постоянно обогащаться новыми идеями и взглядами на свои задачи. Но всякие более глубокие идеи с неизбежностью требуют более выразительных и тонких математических средств. В частности, язык и методы линейной алгебры стали неотъемлемым средством во многих работах по математической статистике и теории обра­ ботки измерений (см., например, Н. И. Идельсон [16], А. Н. Кол­ могоров [18], Ю. В. Линник [20], Т. Андерсон [1], С. Р. Рао [28], А. И. Мазмишвили [23], Ю. В. Кемниц [17] и др.).

В книге последовательно рассматриваются некоторые фраг­ менты линейной алгебры и показывается, как они могут быть при­ менены к теории уравнивания измерений. При этом упор делается на то, чтобы показать, что методы линейной алгебры открывают

пути к новым

идеям и возможностям, а не только лишь

ведут

к упрощению

формы изложения общеизвестных фактов.

Так,

в книге приведено несколько современных методов решения систем нормальных уравнений, а также непосредственно условных уравнений и уравнений поправок. Среди них — метод, позволя­ ющий получить любое число произвольных элементов решения, не располагая уравнения и неизвестные специальным образом. Обсуждаются вопросы обусловленности систем уравнений и воз­ можности улучшения обусловленности. Но нужно заметить, что сами новые идет теории уравнивания измерений могут предъяв­ лять свои требования к средствам линейной алгебры. Так, необхо­ димость рассмотрения в книге вопросов алгебры каркасов и ал­ гебры псевдообратных матриц, а также геометрии подпространств в большой степени вызвана рядом вопросов теории уравнивания измерений. С помощью этих средств стало возможным рассматри­ вать уравнивание по методу наименьших квадратов как частный случай решения общей задачи уравнивания, включающей в себя все «приближенные» уравнивания. В связи с этим появляются выразительные средства сравнения разных решений задачи урав­ нивания между собой по степени обусловленности, по вероятност­ ным характеристикам и т. д. Удается так видоизменить обоснова­ ние метода наименьших квадратов, что получается простая

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

Все сказанное определяет монографический характер книги. Ее читателями предполагаются творчески активные люди — науч­ ные работники, аспиранты, инженеры и студенты старших курсов, имеющие интерес к обработке измерений. Впрочем, первые не­ сколько глав написаны весьма подробно и вполне могут быть использованы для начального ознакомления с предметом. В после­ дующих главах трудность чтения заметно возрастает.

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

На разных этапах работы над книгой автор получал поддержку и замечания от проф. Ю. В. Кемница, проф. А. И. Мазмишвили, доц. Ю. М. Неймана, доц. Ю. И. Маркузе, доц. А. П. Тищенко, доц. 3. С. Хаимова и других лиц. Очень большой и кропотливый труд в эту книгу внес ее редактор проф. А. В. Гордеев. Всем им автор выражает свою искреннюю признательность.

Дроздов Н. Д

Г л а в а 1

МАТРИЦЫ И ЭЛЕМЕНТАРНЫЕ ДЕЙСТВИЯ НАД НИМИ

§1. МАТРИЦЫ

Представим себе, что нужно записать систему 400 линейных алгебраических уравнений с 400 неизвестными с тем, чтобы что-то сказать о ней. Ясно, что обычными полиграфическими средствами это сделать невозможно. Поэтому в книгах всегда прибегают к каким-нибудь «хитростям». Одна из них состоит в том, что вместо всей системы пишут символ

 

+ а12Х2 + • • • + аХП

=

 

а21хг

+ а22х2 + . . . +

а2пхп

= b2,

(I)

amixi

чт2х2 -{-...-{-

атпхп

Ьт>

 

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

В самом частном случае система (1) может состоять из одного уравнения с одним неизвестным, т. е. иметь вид

ах= Ъ.

В связи с этим естественна мысль всякую систему (1) изобразить

символически в виде

 

А х = Ь .

(2)

Теперь символам А, а; и & следует придать вполне определен­ ный смысл, причем так, чтобы они в собирательной форме харак­ теризовали соответственно все коэффициенты, неизвестные и сво­ бодные члены системы (1). Символам же Ах, который формально будем называть произведением А на х, и Ах = Ь, который фор­ мально называется равенством А г и Ь, необходимо также придать такой смысл, чтобы равенство (2) действительно выражало си­ стему (1). Реально это можно оформить по-разному (имеется, на­ пример, краковяновая система записи [39, 36, 26]). Одним из самых распространенных оформлений является такое, когда

символ А понимается как

следующая

таблица, составленная

из коэффициентов системы (1),

 

 

 

а11 Я 1 2

• • •

а

А = 1 Ы =

#21

• • •

&2п

символы х и Ь — как таблицы-столбцы, составленные из неизвест­ ных системы (1) и ее свободных членов:

 

 

 

хх

 

 

 

 

X =

Х-2 ; 6 = h

 

 

 

 

Все записанные таблицы называют матрицами

[6, 7,

19, 24,

33,

34].

 

 

 

 

 

 

В общем случае матрицы состоят из тп чисел (і =

1 , 2 , .

. ., т;

/ =

1, 2, . . ., п) — элементов матрицы — и

содержат т строк

и п столбцов. Числа тип

определяют размеры

матрицы. Говорят,

что

матрица А является т X тг-матрицей, понимая при этом, что

она

состоит из т строк

и

п столбцов.

 

 

 

Иногда, когда это необходимо для понимания, размеры мат­ рицы записываются в виде нижнего индекса при символе матрицы, например АтХп. В общем случае размеры матрицы не указы­ ваются. В символе аГп как правило, первый индекс означает номер строки, а второй — номер столбца, на пересечении которых стоит этот элемент.

При т = п матрица называется квадратной матрицей п-то порядка. Элементы квадратной матрицы с одинаковыми номерами

строк и столбцов называют диагональными.

Множество всех

диа­

гональных элементов матрицы х1,

а 2 2 ,

. . ., апп)

образуют ее

диагональ.

 

 

 

 

 

 

 

 

 

 

Когда нужно подчеркнуть, что рассматриваются произвольные

по размерам т X гс-матрицы, применяют

термин

«прямоуголь­

ные

матрицы».

При этом

матрица

называется

горизонтальной,

если

т _п, невертикальной,

если

 

 

п, вертикальной

при

т >

п и негоризонтальной

при т ^

п.

 

 

 

 

 

Матрицы х и Ъ в выражении (2) имеют частный вид: они содер­

жат лишь один

столбец

(п = 1). Такие

матрицы называют

мат­

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

просто

столбцами. Говорят,

что столбец

имеет т-Ш порядок, если число строк в нем равно

т.

 

Матрицы, имеющие только одну строку =

1) и п столбцов,

называются матрицами-строками

или

просто

строками

п-го

порядка.В обозначении элементов строки или столбца один индекс обычно опускают.

В дальнейшем матрицы-строки и матрицы-столбцы обозна­ чаются малыми латинскими буквами, а матрицы других видов —

большими латинскими буквами.

 

§

2. ЭЛЕМЕНТАРНЫЕ ДЕЙСТВИЯ НАД МАТРИЦАМИ

 

После

того,

как символам А, х и Ъ в выражении (2)

придан

смысл матриц,

легко определить смысл произведения

матриц

Ах так, чтобы выражение (2) означало систему (1). Для этого

нужно просто положить,

что

 

 

 

ап

ах2

 

 

 

Ах = ^21 ^22

а2п

 

 

а11х112

х2 +

 

 

^21 Х1

~Ь ^22

~f" • • • "Т &2п

ХП

(3)

 

 

 

 

amlxl

~\~ &т2х2 "f" • • • ~Ь атпхп

 

Иначе, произведение матрицы А на столбец х нужно понимать как столбец, любой г-й элемент которого равен сумме произведе­ ний элементов Ї-Й строки матрицы А на соответствующие по номеру

элементы столбца х.

Если

воспользоваться

символом сложения

 

 

 

 

т

Ч

ацхі +

ai2x2

+ ...

+ аіпхп

= 2

a-ip],

то (3) запишется так:

 

 

 

J - 1

 

 

 

 

 

 

а11 а 1 2 • • • %п * 1

2

 

&210,22

• •

а2п

х% =

2

(3е )

ат\атг

• "ял

 

2

 

Нужно заметить, что число строк второго сомножителя в про­ изведении Ах должно быть равно числу столбцов первого сомно­ жителя. Число строк первого сомножителя при этом может быть произвольным. В частности, можно говорить о произведении одной строки ?г-го порядка на столбец п-го порядка

II «її яг, • • • а„ х2 = а1х1 + а2х2 + ... + апхп =

= 2 яр}.

п

Так что символ 2 aixj можно считать произведением строки /-і

на столбец.

С этой точки зрения произведение матрицы А на стол­

бец х есть

столбец,

і-ж элемент

которого равен произведению

і-й строки матрицы А

на столбец

х.

Из соображения отождествления выражений (1) и (2) опреде­

ляется также смысл равенства

двух столбцов

2

<hjxj

ъ,

2

0,2 jX j

=

2

amjxj

ьт

Для этого нужно лишь потребовать, чтобы соответствующие

по номеру

элементы

равных

столбцов были равны между собой,

 

п

 

 

т. е. чтобы

было 2

аахі

°і П Р И любом і. Важно отметить, что

 

/ - 1

 

 

порядки столбцов должны быть при этом одинаковыми.

Понятие произведения матрицы на столбец естественно обоб­ щается в понятие произведения двух матриц. Пусть число столбцов матрицы А равно числу строк матрицы В. Тогда произведением АВ матрицы А на матрицу В называется матрица С, элемент с,у

которой

получается

как

произведение

г-й

строки

матрицы

А на

;-й

столбец

матрицы

В, т. е.

 

 

 

 

 

 

 

 

 

а11 Я 1 2 • • а1п її || ь п ъ п

• • • V

 

С 11 С 12 • • •

с

 

 

 

а 2 1 а 2 2

• «2/1

1

^21&2 2

"

 

^21 ^22

• •

с

 

 

 

ffimlffim2

• О-тп II

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

71

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сц

= аЪ1}

+

ai2b2j

+ . . + ainbni

= 2

ЩФкі

 

 

 

 

 

 

 

 

 

 

 

 

 

fe-i

 

 

 

при

і =

1,

2,

. . ., т; /

=

1, 2, . . ., г.

следующим

образом

обоб­

Понятие

равенства

двух

столбцов

щается на случай двух матриц. Пусть матрицы А =

 

|| а,;-|| и В =

=|| fyyll имеют одинаковые размеры. Тогда говорят, что они

равны между собой и пишут А = В или || at!\\ = || если все элементы одной из них равны соответствующим по индексам

элементам

другой, т. е. atj =

bt-r

Поэтому можно, например,

 

 

п

 

записать,

что \\alk\\ \\ bkS\\ =

||2

aikbki\\-

Помимо действия умножения над матрицами, вводится еще несколько элементарных операций. Как известно из элементар­ ной алгебры, если с = ах и d = Ьх, то

с + d = ах + Ъх = (а + Ь) х.

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

A z + B a ; = ( A + B ) x .

 

(4)

Это требует установления понятия суммы матриц.

Покажем,

что если под суммой двух матриц А =

|| а1}\\ и

В =

|| ЪГ1\\ оди­

наковых размеров понимать матрицу

А + В тех

же

размеров,

любой элемент которой

получается

как сумма

соответствующих

по индексам элементов матриц А и В, т. е.

 

А + В = || аи

М

= К/+&//II,

 

то будет выполняться равенство (4). Действительно,

II л

и

||

п

||

|| п

п

I h-1

II

II

ft=l

• II

II h=l

h=l

= 2

 

+

 

= ( A + B ) x .

 

Ilft-i

 

 

 

 

 

 

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

S = A + B + - . . + K

5)

образуются по правилу

*І/ = <**/-|-ЬІН

1-А//.

В частном случае, когда в выражении (5) стоит к равных слагаемых, пишут

S = A + A - )

[-А-кА = Ак.

Иначе

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

Соседние файлы в папке книги из ГПНТБ