Файл: Скворцов, М. И. Теория и практика решения задач кораблевождения с учетом влияния систематических ошибок учебник.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 23.10.2024
Просмотров: 76
Скачиваний: 0
выдающиеся советские ученые — моряки и геодезисты про фессора В. В. Кавранский, А. С. Чеботарев, Н. Н. Мату-
севич, доцент |
В. Н. Зимовнов [26] |
и другие. |
|
A. С. Чеботарев относил эту проблему к |
числу наибо |
||
лее насущных |
в области обработки |
результатов наблюде |
ний [79, стр. 31]. В течение последних лет проблеме систе матических ошибок был посвящен ряд серьезных исследо
ваний, но считать ее окончательно решенной еще |
нельзя |
[33]. |
|
B. В. Каврайский в одной из своих последних работ, к |
|
сожалению оставшейся неопубликованной, высказал |
пред |
положение, что в частном случае, когда число искомых
величин |
равно |
двум (координаты |
обсервованного места ко |
||
рабля), |
число |
измерений — трем |
(три |
астрономические ли |
|
нии. положения), причем измерения |
отягощены |
постоян |
|||
ной (повторяющейся) систематической ошибкой, |
вектор |
оценок искомых величин, обладающих наименьшими, ди сперсиями, должен представлять собой линейную комби нацию векторов оценок искомых величин, доставляемых способом наименьших квадратов и способом исключения систематической ошибки. Впоследствии эта идея, допол ненная соображением М. М. Лескова о том, что эффектив ная оценка исключаемой -систематической ошибки должна отыскиваться как среднее весовое из двух ее оценок — найденной из предыдущих наблюдений и доставляемой спо
собом исключения [47, стр. 44—45], легла в |
основу способа |
|||
отыскания |
обсервованного |
места, |
предложенного |
|
В. |
В. Вейхманом [10]. Решение проблемы, |
основанное на |
||
ее |
сведении |
к задаче уравнивания зависимых результатов |
измерений, было предложено В. Т. Кондрашихиным [37], [38], В. А. Коугия [41], В. М. Полянским [58] и другими авторами.
Способы обработки результатов измерений, отягощен ных систематическими ошибками, можно разделить на че тыре основные группы.
Способ А — способ наименьших квадратов в его клас сическом виде (§ 1.2). При его применении мы ограничи
ваемся минимальным числом |
«основных» искомых вели |
чин, ради которых, собственно, |
и выполняются измерения, |
и полагаем, что измерения отягощены только случайными ошибками, т. е. матрица КА = М (ДДт ) вектора истинных остаточных ошибок диагональна.
Нельзя -не видеть произвольности и условности этих
36
предположений. Так, при определении места корабля по результатам практически одновременных измерений двух навигационных параметров мы считаем «основными» иско мыми величинами географические координаты, при опре делении места по ряду разновременных измерений навига ционных параметров — также . и проекции вектора скоро сти течения, а иногда и их производные по времени, при определении маневренных элементов корабля на мерной линии — еще и поправку лага и аванс. Однако будем счи тать, что в каждом конкретном случае молено конкретно указать смысл, который придается понятию «основные ис комые величины».
Не меньше произвола и в предпололсении, что корреля ционная матрица АГД является диагональной. Как видно
из выражения (1.32), каледая из компонентов At вектора Д слагается из собственно ошибки измерения Д^ ошибки аппроксимации Дг* и ошибки Д™( возникающей вслед ствие случайных отклонений искомых величин от их ма тематических ожиданий. Стремясь уменьшить влияние слу чайных ошибок, мы нередко выполняем каждое измерение по нескольку раз и подставляем в выражение (1.28) сво бодного члена уравнения поправок среднее арифметиче ское результатов измерений. При этом дисперсии случай ных ошибок свободных членов уравнений поправок умень шаются пропорционально числу осредияемых результатов измерений, систематические л<е ошибки остаются неизмен ными и начинают преобладать по величине над случайны ми даже в тех случаях, когда такое преобладание вначале отсутствовало. Что же касается ошибок Д/ и Д " , то основании для предположения о диагональное™ их кор реляционных, матриц имеется еще меньше.
Обратимся |
теперь |
к |
каноническому |
представлению |
||
(1.10) |
вектора |
Д: |
|
|
|
|
|
|
|
|
+ |
|
(1.88) |
Поскольку |
при Z Ф 0 |
корреляционная |
матрица |
векто |
||
ра FZ |
недиагональна |
(§ |
1.1), ясно, что лелощее в |
основе |
способа наименьших квадратов предположение о диаго-
нальности корреляционной |
матрицы Кь |
может выпол |
няться лишь при условии, |
что |
|
2 |
= 0, |
(1.89) |
37
следствием которого являются равенства |
|
|
КА = Кй; |
Р=(-т-кУ., |
(1.90) |
где Ка—корреляционная |
матрица вектора |
8. |
Рассмотрим, к каким искажениям вектора \ оценок ис комых величин будут вести эти предположения в случаях, когда в действительности условие (1.89) не выполняется. Если поправки, вводимые в результаты измерений, явля ются несмещенными оценками компенсируемых ими си стематических ошибок, то в соответствии со сказанным в § 1.1 о свойствах остаточных систематических ошибок
М ( Д ) = 0 , |
(1.91) |
откуда, повторив вывод (1.43) — (1.45), |
получим |
М(Г) = ЛГ(«). |
(1.92) |
Этот результат имеетпринципиальное значение, по скольку показывает, что способ наименьших квадратов, по праву считающийся одним из основных способов обработ ки результатов измерений и познания количественных за кономерностей окружающего нас реального мира, способен доставлять несмещенные, не искаженные систематическими погрешностями сведения об этих закономерностях не только в тех случаях, когда выполнены условия Гаусса — Колмо горова, но и в тех, когда считать остаточные систематиче ские ошибки пренебрежимо малыми нельзя.
Чтобы судить о том, как влияет пренебрежение систе матическими ошибками на оценку корреляционной матри цы вектора \ оценок искомых величин, обратимся к выра жению (1.14). Учитывая соотношения (1.35), (1.37) и то •обстоятельство, что матрица £п р приближенных значений искомых величин — не случайная матрица, а также то, что, как видно из выражений (1.12) и (1.43), при пользовании способом наименьших квадратов
G = |
(АТРА)~1АТР, |
(1.93) |
получим |
|
|
К~ = GFKZ |
(GF)T+ 4 ) (АГРАГ\ |
(1.94) |
£ |
|
|
38
Здесь матрица GFK7,(GFy отражает то влияние оста точных систематических ошибок, которым, производя урав нивание результатов наблюдений способом наименьших квадратов, мы пренебрегаем. Если надо оценить влияние, которое оказывает на точность оценок искомых величин некоторая конкретная систематическая ошибка, удобно пользоваться выражениями (1.19) и (1.20). Пример при
веден в § |
3,6. |
|
|
|
|
|
Способ |
В — способ |
исключения систематических |
оши |
|||
б о к — упоминается уже |
Ф. Р. Гельмертом |
[17, стр. 19]. Его |
||||
разновидностями |
являются |
предложенные |
Гауссом |
спосо |
||
бы организации |
измерений |
и обработки |
их результатов, |
иногда именуемые способами замещения, компенсации по грешностей по знаку и противопоставления [51, стр. 253— 257], [52, стр. 167—168]. Пользоваться способом В для об работки наблюдений, отягощенных систематическими ошиб ками, рекомендовали А. П. Ющенко [82, стр. 108], В. Ф. Дья конов [25] и другие авторы. Единственное его отличие от способа А заключается в том, что амплитуды поправок, предназначенных для компенсации систематических оши бок, включаются в число искомых величин, причем воз можные значения этих амплитуд какими-либо предположе
ниями не |
ограничиваются. |
|
|
|
|
Вспомним, что оценки искомых величин, доставляемые- |
|||||
способом |
наименьших |
квадратов, |
с точностью до вели |
||
чин второго порядка малости не |
зависят |
от тех |
произ |
||
вольных |
приближенных |
значений |
Е 1 п р , |
^ п Р , |
••• . Ет п Р |
искомых величин, которые при вычислении свободных чле нов уравнений поправок подставляются в формулу (1.28). Это полностью справедливо и по отношению к амплиту дам поправок, предназначенных для компенсации систе матических ошибок, включаемым в число искомых вели чин. Какие бы значения этим амплитудам ни были при даны при вычислении поправок, которыми исправляются результаты измерений, вводимые в выражение -(1.28), оценки искомых величин, доставляемые способом В, изме няться не будут. Следовательно, поправки, которыми ис правляются результаты измерений, могут вычисляться ис ходя из любых произвольных приближенных значений этих амплитуд. Например, если место корабля определяется по
высотам трех светил и исключается |
постоянная (повто |
|||||
ряющаяся) |
ошибка |
измерения |
высот |
(обсервованное |
ме |
|
сто отыскивается |
в точке |
пересечения биссектрис |
тре |
|||
угольника |
погрешностей), |
то |
отсчеты |
секстана могут |
ис- |
39