
книги из ГПНТБ / Скворцов, М. И. Теория и практика решения задач кораблевождения с учетом влияния систематических ошибок учебник
.pdfУчитывая влияние как случайных, так и систематиче ских ошибок измерении, среднюю квадратическую ошибку обсервованного места следует оценивать формулой
/s
где s — число систематических ошибок.
Во многих пособиях по навигации рекомендуется оце нивать точность обсервованного места, исходя из предпо ложения, что влияние остаточных систематических ошибок измерений существенно не отличается от влияния случай
ных ошибок, т. е. пользуясь |
формулой (2.125) и |
подстав |
|
ляя в нее вместо величин |
<з\, |
02 средние квадратическне |
|
величины сумм случайных |
и |
систематических |
ошибок. |
Представление о погрешностях, которые возникают при таком оценивании, дает следующий пример. Как известно,
при определении |
места корабля по визуальным |
пеленгам |
|
1 ориентиров преобладают постоянные |
систематические |
||
ошибки (ошибка |
определения поправки |
компаса, |
измене |
ние среднего значения ошибки курсоуказания за время, прошедшее после ее определения, случайное отклонение мгновенного значения ошибки курсоуказания от средне го). Примем среднюю квадратическую величину их суммы
равной |
<Jr = 0,6°, среднюю квадратическую |
величину |
суммы |
случайных ошибок измерения и прокладки |
пеленгов |
о с л = |
|
= 0,3° |
(преобладание по величине систематических ошибок |
измерений над случайными нередко бывает еще более зна
чительным). |
Оценим |
среднюю квадратическую |
|
ошибку |
||||||
определения места корабля при расстояниях до |
|
ориенти |
||||||||
ров Д 1 = Д 2 = Ю0 |
каб |
сначала |
|
формулой |
(2.125), |
полагая |
||||
а\ — <z1 = |
Yа1 |
+ о 2 |
л = 0,67°, |
а |
затем |
формулами |
(2.125), |
|||
(2.127) и |
(2.129), |
полагая а г |
= 0,6°; |
о С л = 0,3°. |
Результат |
|||||
приведен |
в табл. |
2.3. |
|
|
|
|
|
|
|
|
Приступим |
к |
рассмотрению |
общего |
случая |
определе |
ния места корабля, когда число линий положения более двух. Применяемые при этом способы отыскания обсер вованного места по линиям положения, образующим фи гуру погрешностей, будем обозначать так же, как рассмо тренные в § 1.3 аналитические способы обработки резуль татов измерений, которым они эквивалентны.
Способ А. Обсервованное место отыскивается в пред положении, что причиной образования фигуры погрешно-
100
Т а б л и ц а 2.3
Оценки средней квадратической ошибки обсервации по двум
пеленгам |
в кабельтовых при |
= Д.г — 100 каб; аг — 0,6°; |
|||||||
|
|
|
я с л |
- 0,3° |
|
|
|
|
|
|
|
|
Угол (] между направлениями на ориентиры, |
град |
|||||
Формула |
|
|
|
|
|
|
|
|
|
|
|
|
30 |
60 |
|
90 |
120 |
|
150 |
(2.125) |
|
|
3,3 |
1.9 |
|
1,7 • |
1,9 |
|
3,3 |
(2.125), |
(2.127), |
1.8 |
1.5 |
|
1.7 |
2.3 |
|
4,3 |
|
(2.129) |
|
|
|
|
|
|
|
|
|
стей могут быть только случайные |
ошибки |
(систематиче |
|||||||
ские ошибки |
пренебрежимо |
малы). |
Графоаналитические |
||||||
приемы |
отыскания |
обсервованного |
места, |
соответствую |
|||||
щие этому предположению |
(центриграфнческий, |
построе |
|||||||
нием полигонов весов и переносов, |
лроведением |
противо- |
|||||||
медиан треугольника погрешностей), и оценивания |
влия |
||||||||
ния случайных ошибок измерений |
на его точность |
описаны |
в любом современном учебнике навигации и мореходной
астрономии, |
поэтому |
останавливаться |
на этих |
приемах |
нет необходимости. |
Аналитическое оценивание координат |
|||
обсервованного места и корреляционной матрицы |
вектора |
|||
их ошибок, |
обусловленных влиянием |
случайных |
ошибок |
измерений, осуществляется в соответствии с общими пред писаниями способа наименьших квадратов (§ 3: 4). Влияние систематических ошибок измерений на точность обсерво
ванного - места |
может быть |
оценено |
формулами (1.19), |
||
(1.20), |
(1.94) |
(пример — в § |
3.6). |
Но |
значительноt более |
удобны |
таблицы готовых ответов |
(§ 3.11). |
Иногда может оказаться полезной следующая интер претация формулы (1.20). Пусть каким-либо графическим или графоаналитическим способом по нескольким линиям положения найдена' обсервованная точка А, причем по правки, введенные в результаты. измерений для компенса
ции r-й систематической |
ошибки, отыскивались |
в |
виде |
AUir = ZrCjr (а , Р ( .... ); |
амплитуда ZrC известна |
со |
сред |
ней квадратической ошибкой аг с . Изменим каждую из этих
поправок на величину о,с ] г (а., |
£)..'. .), проложим линии |
положения, • соответствующие |
полученным значениям на- |
101
вигационных параметров, и применим к ним тот же алго ритм, каким из первоначальных линий положения отыски
валась точка А. Отстояние найденной |
таким образом |
точ |
|
ки от точки |
А будет равно (в избранном для прокладки |
||
масштабе) |
векториальной ошибке тг |
обсервованного |
ме |
ста, выражающей влияние рассматриваемой r-п система
тической |
ошибки |
измерений. |
|
|
В некоторых частных случаях этот прием приводит к |
||||
простым |
и удобным |
выражениям. Так, если место |
кораб |
|
ля определяется |
по |
радиомаякам, расположенным |
на од |
ной прямой линии, линия курса корабля ей параллельна (близкие ситуации часты в практике кораблевождения) и рассматривается влияние ошибки коэффициента D радио
девиации, величина тг в милях |
|
|
|||
|
|
|
|
|
( 2 - 1 3 0 ) |
где d—расстояние |
от |
корабля до прямой, на которой |
|||
расположены |
радиомаяки, |
мили; |
|
||
aD—средняя |
квадратическая |
величина |
ошибки в |
||
коэффициенте |
D |
радиодевиации, град. |
|||
|
|
<--> |
|
|
|
Векториальная |
ошибка |
тг |
перпендикулярна |
линии кур |
са корабля. Все обсервованные точки располагаются мо ристее фактической линии пути корабля при положитель ной ошибке в коэффициенте D и ближе к берегу при от рицательной ошибке.
При рассмотрении способов отыскания обсервованного места корабля нас будет интересовать точность не только
самой |
обсервации |
непосредственно, |
но и доставляемых |
этим |
способом |
оценок амплитуд |
поправок, которыми |
впредь должны исправляться результаты измерений для компенсации систематических ошибок. Когда применяется способ А, эти амплитуды в число искомых величин не включаются. Следовательно, и после его применения мы будем исправлять результаты измерений поправками, имеющими прежние счислимые амплитуды, средняя ква дратическая величина амплитуды остаточной систематиче
ской |
ошибки |
не изменится: |
|
|
||
|
zr(A) = |
0> |
Z r W — ZrZ |
+ z r ( A ) = |
2ГгС; |
о г ( А ) = з,с . (2.131) |
Способ |
В — способ исключения |
систематических оши |
||||
бок. |
Обсервованное место |
отыскивается |
в предположении, |
102
что поправки, вводимые в результаты измерении для ком пенсации некоторых систематических ошибок, известны с ничтожно малой точностью, и что той информацией, ко торая в них содержится, можно пренебречь.
В практике кораблевождения "нередко бывает, что одна или несколько систематических ошибок заметно преобла дают по величине над случайными. Если к тому же каж
дый |
навигационный |
параметр измеряется |
по |
нескольку |
|
раз |
и к прокладке |
принимается |
среднее |
арифметическое |
|
из результатов измерений, то это преобладание |
становится |
||||
еще |
более заметным. Признаком, |
его подтверждающим, |
является сохранение конфигурации и размеров фигуры погрешностей в нескольких выполненных одна за другой обсервациях при полной уверенности в правильности опо
знания ориентиров. В таких ситуациях |
применение спосо |
ба В может оказаться оправданным. |
Кроме того, он |
имеет и вспомогательное значение как промежуточная сту пень к отысканию оценок искомых величин способом D.
Если исключению подвергается только одна система тическая ошибка, то практическое применение способа В сводится к очень простым геометрическим построениям, вследствие чего он часто применяется на практике. Та
ковы, |
например, прием отыскания |
обсервованного места |
||
в точке пересечения |
биссектрис треугольника |
• погрешно |
||
стей, |
образованного |
равновесными |
линиями |
положения |
(при определении места по высотам светил исключается постоянная систематическая ошибка измерения высот), и аналогичный прием отыскания обсервованного места по визуальным пеленгам трех ориентиров (исключается по стоянная ошибка поправки компаса). Ряд других частных случаев рассмотрен Н. Н. Матусевичем [55] и М. М. Леско вым [47].
В любом частном случае может быть применен общий прием исключения систематических ошибок, являющийся
геометрической |
интерпретацией |
решения |
уравнений |
|||
(1.95) — (1.104) |
. Пусть надо исключить |
некоторую систе |
||||
матическую |
ошибку |
Cr /r (а, р . . . ) . |
Для |
этого следует сна |
||
чала исправить результаты измерений поправкой |
&Ujr = |
|||||
— Zrfr(a., |
р / |
. . . ) , |
задавшись |
произвольным |
(например, |
счислимым) значением амплитуды Zr исключаемой систе матической ошибки, и проложить линии положения /, //,
/// и т. д., соответствующие полученным значениям нави гационных параметров (рис. 2.4), затем исправить резуль-
103
таты измерении поправкой A'Ulr = (Zr+Ar) /, (а,, Р( ....). где А,-— произвольное число, и проложить смещенные линии
положения /', |
//', |
III' и т. д., соответствующие |
изменив |
|
шимся значениям |
навигационных параметров. |
|
|
|
Прямые, соединяющие одноименные точки |
(13 |
и 13', 23 |
||
и 23', 24 и 24') |
пересечения первоначальных и |
смещенных |
линий положения, соответствуют навигационным параме трам вида U.i,=frl,Ui—jnUi,. Поэтому будем назы-
|
|
Рнс. 2.4. Графическим прием исключения система |
|
|
тической ошибки из фигуры погрешностей |
вать |
их |
изоразностными линиями положения. Очевидно, |
что |
они |
свободны от влияния исключаемой систематиче |
ской ошибки (лишь исключаемой, но могут быть отягоще ны другими систематическими ошибками, иным образом зависящими от параметров, характеризующих условия из мерений).
При отыскании обсервованной точки В внутри фигуры погрешностей, образованной изоразностными линиями по ложения, следует считать, что они отягощены только слу чайными ошибками. Амплитуда поправки, соответствую щая условию полной компенсации исключаемой система
тической ошибки, может быть |
найдена по |
формулам |
|
^ S - 2 |
" ^ - " |
( 2 Л 3 2 ) |
|
2 л |
/ = i |
|
|
i=I |
|
|
|
^ < * , = |
^ |
( 2 . |
1 3 3 ) . |
104
где
Pi =
UI (B) = ^i\xB> Ув) |
— значение |
1-го навигационного пара |
||||||
|
метра, вычисленное исходя из ко |
|||||||
|
ординат |
точки |
5 ; |
|
|
|
||
п.—число |
линий |
|
положения; |
|
|
|||
U, — измеренное-значение |
/-го |
навига |
||||||
|
ционного |
параметра, |
исправленное |
|||||
|
поправками |
со |
счислимыми |
значе |
||||
а2(1) |
ниями |
амплитуд; |
|
|
|
|||
— дисперсия |
ошибки |
измерения, |
вес |
|||||
|
которого |
принят равным единице; - |
||||||
g[ |
— градиент |
i-.ro |
навигационного пара |
|||||
|
метра; |
|
|
|
|
|
|
|
о] — дисперсия |
случайной |
ошибки |
1-го |
измерения.
Способ В далеко не всегда приводит к меньшим дис персиям оценок искомых величин, нежели способ А. В слу чаях когда приходится делать выбор между этими спосо бами, следует основываться на сравнении точностей до ставляемых ими результатов. При этом можно руковод ствоваться общими выражениями (1.19), (1.20), (1.14). Если пользоваться только логарифмической и арифмети ческой линейками, то вычисления окажутся трудоемкими. Па помощь могут прийти таблицы готовых ответов. Такая таблица, составленная для частых в практике корабле вождения обсерваций по трем равновесным линиям поло жения (по высотам трех светил, по визуальным пеленгам трех равноудаленных от корабля ориентиров и т.. д . ), и пример пользования ею приведены в § 3.11.
Косвенными признаками, позволяющими судить о до пустимости .применения способа В, являются углы пере сечения изоразностных линий положения и величина об разованной ими фигуры погрешностей. Если углы пере сечения изоразностных линий положения малы, то от при менения способа В лучше воздержаться. Большая величи на образованной изоразностными линиями положения фи гуры погрешностей свидетельствует о том, что предполо жение о пренебрежимой малости случайных ошибок изме рений и всех других систематических ошибок,' кроме ис-
105
ключаемой, лежащее в основе способа, не подтвердилось, следовательно, применение способа В нельзя считать оправданным.
Заметим, что если п — число исходных линий положе ния, то число взаимно независимых изоразностных линий положения будет равно п—1. В частности, если число ис ходных линий положения равно трем, то изоразностные линии положения всегда пересекутся в одной точке; воз можность контроля правильности измерений и геометри
ческих построений |
по величине фигуры погрешностей бу |
дет отсутствовать. |
Поэтому результаты, полученные из |
трех исходных линий положения способом В, следует под вергать особо тщательной проверке (лучше всего — по вторными обсервациями по другим комбинациям ориен тиров, применяя для вычисления поправок только что найденное значение амплитуды Zr В ) ) .
Как отмечалось в § 1.3, включение амплитуды Zr в число искомых величин может приводить к неопределен
ности, к тому, что система нормальных |
уравнений не |
бу |
|
дет иметь единственного конечного |
решения (типичный |
||
пример — случай неопределенности |
при |
определении |
ме |
ста по пеленгам нескольких ориентиров, когда все ориен
тиры и корабль находятся на одной окружности |
и делает |
||
ся |
попытка исключить |
постоянную систематическую ошиб |
|
ку |
пеленгования). Если |
точка В отыскивается |
описанным |
выше графическим приемом (рис. 2.4), признаком, указы
вающим |
на ситуацию |
неопределенности, |
является |
взаим |
||
ная параллельность или приблизительная |
параллельность |
|||||
изоразностных |
линий |
положения. |
|
|
||
В § 1.3 было показано также, что если в ситуации не |
||||||
определенности |
амплитуда |
Z r це будет включена |
в число |
|||
искомых |
величин, то |
эта' |
систематическая |
ошибка |
не мо |
жет приводить ко взаимному несогласию результатов из мерений (в графической интерпретации — к образованию фигуры погрешностей исходными линиями положения). В случаях когда измерения отягощены несколькими си стематическими ошибками, зависимости которых от пара метров, характеризующих условия измерений, различны, этим обстоятельством можно пользоваться для последо вательного уточнения оценок амплитуд поправок, служа щих для компенсации этих систематических ошибок. На пример, из систематических ошибок измерений пеленгов корабельным радиопеленгатором часто преобладают по ве личине ошибки коэффициентов А и D радиодевиации [22], 106
[72], [73]. Когда ошибка в коэффициенте А не может вызы вать несогласия измерений, можно способом В уточнить
оценку коэффициента D радиодевиации, а когда |
возник |
нет ситуация неопределенности для исключения |
коэффи |
циента D, уточнить оценку коэффициента А. |
|
Способ D — способ последовательного уточнения ам плитуд поправок, предназначенных для компенсации систе матических ошибок. Если измерения отягощены и случай
ными, и систематическими ошибками, |
способ D всегда |
при |
|
водит к оценкам искомых величин, |
обладающим меньши |
||
ми дисперсиями, чем оценки, доставляемые |
способом А |
||
или В. Случая неопределённости при его применении |
воз |
||
никнуть не может. Следует помнить, |
что, когда применяет |
||
ся способ D, результаты измерений |
должны |
исправляться |
всеми учитываемыми поправками с наибольшей тщатель ностью. Пользоваться произвольными приближенными зна чениями амплитуд, как в способе В, недопустимо.
В часто встречающихся случаях, когда одна из систе матических ошибок заметно преобладает по величине над остальными, а также в ситуациях, когда другие система тические ошибки не могут приводить к внутреннему несо гласию результатов измерений (образованию фигуры по грешностей), можно пользоваться простым приемом оты скания обсервованного места, идея которого была пред ложена профессором В. В. Каврайским, а обоснование выражено формулой (1.130). Практическое применение этого приема сводится к следующим действиям: найти об-
сервованные точки способами А и В, |
измерить |
расстоя |
||
ние d между ними, от точки А отложить по |
прямой АВ |
|||
отрезок |
Ы. |
|
|
|
Как |
видно из выражений (1.130), |
(2.131) |
и |
(2.133), |
амплитуда поправки, предназначенной |
для |
компенсации |
г-й систематической ошибки, которая в способе В исклю чалась бы, может быть найдена по формуле
|
Я( 0 ) |
= ^ с + ^ ( Л ( |
В ) - ^ с ) , |
(2.134) |
|
где |
^-г{В) — оценка |
r-й амплитуды, |
найденная |
||
|
ZrC |
способом |
В |
по формуле |
(2.133); |
|
— счислимое |
значение амплитуды. |
|||
Введем |
обозначения: |
|
|
|
107
A^gj |
='Cr -\- Zr^ |
•—истинная |
ошибка |
оценки амплиту |
|||||||||||
|
|
|
|
|
ды г-н поправки, доставляемой спо |
||||||||||
ДгС = l r + |
|
|
собом |
В; |
|
|
|
|
|
|
|
||||
Z*rlг С |
|
истинная |
ошибка |
счислимого зна |
|||||||||||
|
|
|
|
|
чения амплитуды |
этой поправки. |
|
||||||||
Рассмотрим |
корреляционную |
матрицу |
|
вектора |
|
||||||||||
|
|
|
|
|
|
д |
= |
r(D) |
|
|
|
|
(2.135) |
||
|
|
|
|
|
|
Д.п |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Обратим |
внимание |
|
на то, |
что |
формула |
|
(2.134) может |
||||||||
быть представлена |
в |
|
виде |
|
|
|
|
|
|
|
|
||||
|
|
|
Z r(D) |
|
(l-X)ZrC |
|
|
+ |
XZr(B). |
|
(2.136) |
||||
Следовательно, |
если |
учесть, |
что |
ошибки |
амплитуд ZrC |
и |
|||||||||
Z r l B ) |
являются |
взаимно |
независимыми |
|
случайными |
ве |
|||||||||
личинами, то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A* [A rU>Aej = |
( l - * R . |
|
|
С 2 ' 1 3 7 ) |
||||||||
Кроме |
того, |
можно |
|
показать, |
что |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.138) |
|
Таким |
образом, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
КЛ = М(ДДТ ) |
= |
(!-*)<£ |
|
|
(1-Х) |
|
of |
|
||||||
|
|
|
|
|
|
|
|
(2.139) |
|||||||
При |
неоднократном |
применении способа |
D следует учи |
||||||||||||
тывать, |
что |
амплитуда |
практически |
любой |
систематиче |
ской ошибки подвержена случайным изменениям и являет
ся |
случайной функцией времени. Введем |
обозначения: |
||
|
_ С,-(О — мгновенное |
значение амплитуды |
||
|
Сг в момент |
времени |
|
|
|
С,0 — математическое |
ожидание |
сред |
|
|
него значения |
случайной |
функ |
|
|
ции С,- ( 0 ; |
|
|
|
о2 |
= М { [Сг ( 0 — С,0 ]2 j —дисперсия |
случайной функции |
||
|
М О ; |
|
|
|
108
Zru = — C^o— амплитуда поправки, которой мы компенсируем r-ю систематиче скую ошибку, если не распола гаем оценкой ее мгновенного зна чения (взятая с обратным знаком оценка величины do, полученная из результатов предыдущих на блюдений, например среднее зна чение поправки гирокомпаса, оп ределенное из длительных наблю дений, табличные значения по правок высоты светила за накло нение горизонта и астрономиче скую рефракцию и т. д . ) ;
a2r0 = М [(с^о -\- ZrQJ ] —дисперсия оценки Zr n- Поставим себе такую задачу. Пусть нам известны пе
речисленные |
характеристики стационарной |
случайной |
|
функции £г(/) |
и |
ее структурная функция S.(-c); |
кроме |
того, в некоторый момент времени t способом |
D опреде |
||
лена амплитуда |
Zr(D) = — ? г ( 0 — в з я т а я с |
обратным |
знаком оценка мгновенного значения амплитуды /--Доста точной систематической ошибки. Какова оценка Zr^ ам плитуды поправки, которой должна компенсироваться эта систематическая ошибка в момент времени t + x, чтобы ди сперсии остаточных систематических ошибок были мини мальны? Если изменение по времени амплитуды t.r(t) мо жет быть описано как стационарный случайный процесс, то решение окажется подобным изложенному в § 2.1 для поправки компаса. Поэтому рассмотрим случай, когда
—нестационарная случайная функция.
Для простоты будем считать, что при первом приме нении способа D за счислимое значение амплитуды по правки, включенной в число искомых величин, принята
оценка ZTQ\
Zrz = Zr0. |
.. (2.140) |
109 .