Файл: Скворцов, М. И. Теория и практика решения задач кораблевождения с учетом влияния систематических ошибок учебник.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