книги из ГПНТБ / Дроздов Н.Д. Линейная алгебра в теории уравнивания измерений
.pdfн . д . ДРОЗДОВ
Л И Н Е Й Н А Я |
А Л Г Е Б Р А |
ВТ Е О Р И И
У Р А В Н И В А Н И Я
И З М Е Р Е Н И Й
ИЗДАТЕЛЬСТВО «НЕДРА»
Мо с к в а - 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 + • • • + а1ПХП |
= |
|
|
а21хг |
+ а22х2 + . . . + |
а2пхп |
= b2, |
(I) |
amixi |
~Ь чт2х2 -{-...-{- |
атпхп |
— Ьт> |
|
в котором все, что не умещается на странице, заменяется точками. Но почему тогда не пойти на максимально возможные графиче ские сокращения?
В самом частном случае система (1) может состоять из одного уравнения с одним неизвестным, т. е. иметь вид
ах= Ъ.
В связи с этим естественна мысль всякую систему (1) изобразить
символически в виде |
|
А х = Ь . |
(2) |
Теперь символам А, а; и & следует придать вполне определен ный смысл, причем так, чтобы они в собирательной форме харак теризовали соответственно все коэффициенты, неизвестные и сво бодные члены системы (1). Символам же Ах, который формально будем называть произведением А на х, и Ах = Ь, который фор мально называется равенством А г и Ь, необходимо также придать такой смысл, чтобы равенство (2) действительно выражало си стему (1). Реально это можно оформить по-разному (имеется, на пример, краковяновая система записи [39, 36, 26]). Одним из самых распространенных оформлений является такое, когда
символ А понимается как |
следующая |
таблица, составленная |
|
из коэффициентов системы (1), |
|
|
|
|
а11 Я 1 2 |
• • • |
а1П |
А = 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х1+а12 |
х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 • • • |
с1г |
|
|||||||||
|
|
а 2 1 а 2 2 |
• |
• «2/1 |
1 |
^21&2 2 |
" |
|
^21 ^22 |
• • • |
с2г |
|
||||
|
|
ffimlffim2 |
• |
• О-тп II |
|
|
|
|
|
|
|
|
||||
где |
|
|
|
|
|
|
|
|
|
|
|
71 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
сц |
= а1хЪ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 - ) |
[-А-кА = Ак. |
Иначе
Это правило естественно распространить на случай, когда число к не только натуральное, но и любое действительное (или даже комплексное). Так определяется произведение матрицы на число. В дальнейшем, если не будет оговорено противное, рас сматриваются только действительные множители к.