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

книги из ГПНТБ / Скворцов, М. И. Теория и практика решения задач кораблевождения с учетом влияния систематических ошибок учебник

.pdf
Скачиваний:
15
Добавлен:
19.10.2023
Размер:
6.54 Mб
Скачать

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

/s

где s — число систематических ошибок.

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

ных ошибок, т. е. пользуясь

формулой (2.125) и

подстав­

ляя в нее вместо величин

<з\,

02 средние квадратическне

величины сумм случайных

и

систематических

ошибок.

Представление о погрешностях, которые возникают при таком оценивании, дает следующий пример. Как известно,

при определении

места корабля по визуальным

пеленгам

1 ориентиров преобладают постоянные

систематические

ошибки (ошибка

определения поправки

компаса,

измене­

ние среднего значения ошибки курсоуказания за время, прошедшее после ее определения, случайное отклонение мгновенного значения ошибки курсоуказания от средне­ го). Примем среднюю квадратическую величину их суммы

равной

<Jr = 0,6°, среднюю квадратическую

величину

суммы

случайных ошибок измерения и прокладки

пеленгов

о с л =

= 0,3°

(преобладание по величине систематических ошибок

измерений над случайными нередко бывает еще более зна­

чительным).

Оценим

среднюю квадратическую

 

ошибку

определения места корабля при расстояниях до

 

ориенти­

ров Д 1 = Д 2 = Ю0

каб

сначала

 

формулой

(2.125),

полагая

а\ — <z1 =

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.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 .

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