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

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

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

ошибок. Если применяется способ С, то для этого требу­

ются

дополнительные вычисления.

 

 

 

 

 

 

Во-вторых, для применения способа С,

как

мы уже

говорили, должныбыть известны

(хотя бы

с

точностью

до

постоянного множителя)

матрица

ЛГЛ

или

матри­

цы Д"3

и Кг- Однако рекомендаций,

как их оценить,

способ

не дает. Применяя способ D, мы, оценивая

корреляцион­

ную матрицу вектора оценок искомых

величин,

как один

из

ее блоков оцениваем и

корреляционную

матрицу К%

вектора оценок амплитуд поправок, включенных в число искомых величин. Что же касается корреляционной ма­

трицы Кд

вектора 8

пост-остаточных ошибок, то в тех

случаях,

когда можно

считать ее диагональной, она мо­

жет оцениваться

по

внутреннему согласию результатов

измерений. Таким

образом, мы получаем все сведения, не­

обходимые для того, чтобы, выполнив новые измерения, опять применить к их обработке способ D.

В-третьих, применяя способ

С, приходится

обращать

матрицу

(\'• а~{1)) К±,

порядок

которой

равен числу уравне­

ний поправок. Применение способа D требует обращения

матрицы

( l : o f i ) ) - ^ z ,

порядок

которой

равен

числу ам­

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

дит способ D, не более сложны,

чем в способе В.

3. Как видно из выражений

(1.34) — (1.48), (1.95) —

(1.105) и (1.120) — (1.124), способы А и В являются пре­ дельными частными случаями способа D. Первый соот­

ветствует предположению, что все

ог = 0 и все

х т + г = 0,

т. е. условию, когда, все поправки,

вводимые в

результаты

измерений, считаются известными со столь высокой точно­ стью, что любое последующее уточнение их оценок апри­ ори полагается лишенным смысла. Способ В, наоборот, со­ ответствует условию, что все о г = + °°; Р\ = 0, т. е. случаю, когда поправки,' вводимые в результаты измерений для компенсации систематических ошибок, - считаются извест­ ными со столь малой точностью, что всеми измерениями,

50

выполненными ранее для их определения, можно пре­ небречь.

4. В случае когда рассматривается влияние только слу­ чайных и одной систематической ошибок, вектор оценок «основных» искомых величин, доставляемых способом D, является линейной комбинацией векторов оценок искомых величин, доставляемых способами А и В:

x D = x A + \ { х в - х Х

(ЫЗО)

причем 0< ! I < . + 1 , т. е. справедливо предположение, вы­ сказанное профессором В. В. Кавранским. Формула (1.130) является следствием предыдущего утверждения и равенств (1.124) — (1.128), если учесть также, что в этом случае матрицы FTP2F и Рх превращаются в положи­ тельные числа.

5. В отличие от способа В, при пользовании которым включение амплитуд поправок в число искомых величин может повести к ухудшению точности, с какой отыскива­ ются оценки остальных искомых величин, или даже к не­ определенности, применение способа D, при котором вклю­ чение новых неизвестных в |Число искомых величин сопро­ вождается присоединением к системе уравнений поправок новых уравнений вида (1.82), не может вызвать ухудше­ ния точности оценок искомых величин (дисперсии их оши­ бок уменьшаются или, по крайней мере, не увеличива­ ются).

Г л а в а

2

С Ч И С Л Е Н ИЕ ПУТИ И О П Р Е Д Е Л Е Н И Е

МЕСТА КОРАБЛЯ

§ 2.1. УЧЕТ ПОПРАВКИ ГИРОКОМПАСА

Ошибки курсоуказания будут

нами рассматриваться

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

ИК—истинный

курс

корабля (угол

между плоско­

стью истинного

меридиана и

диаметральной

плоскостью корабля);

 

КК— компасный курс (отсчет курса по репитеру ги­ рокомпаса);

ИП—истинный

 

пеленг (азимут) на ориентир (угол

между плоскостью истинного меридиана и пло­

скостью

вертикала

ориентира);

КП—компасный

 

пеленг

(результат измерения пе­

ленга);

 

 

 

М<—учитываемая

поправка компаса;

i — номер

измерения.

 

Истинная ошибка

i-ro измерения курса

 

 

 

(2.1)

 

 

 

(2.2)

 

 

 

(2.3;

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

52

равновесия ошибки измерения курса и пеленга испыты­ вают изменяющиеся по времени отклонения от некоторых

средних

значений, являются

случайными

функциями

времени:

 

 

 

 

 

 

 

 

 

 

ск =

«ко +

* i

(0;

сп =

£по + * i (0;

 

( 2 4 )

 

г к —

ZK0 ~\~ z\ (0>

Za

~ zn0 "Г" %\ (01

 

где Ск 0 (С„о)среднее значение истинной ошибки изме­

 

 

рения

курса

(пеленга);

 

 

г ко

(~ло) среднее

значение

остаточной

оши'бкн

из­

 

 

мерения

курса

(пеленга).

 

 

Сечение случайной

функции

Z\(t),

соответствующее

мо­

менту времени t, будем впредь называть отклонением мгновенного значения ошибки курсоуказания от среднего. Как показало исследование [69], случайная функция Z\ (t) является практически стационарной. Ее корреляционная

функция хорошо

аппроксимируется

-выражением

 

 

/ ? 2 W =

^ H ) « - e N c o s p , ,

(2.5)

где т — разность

моментов

времени,

в которые

рассматри­

ваются мгновенные значения ошибки курсоуказания. У не­ апериодических гирокомпасов типа «Курс» при плавании в средних широтах без больших изменений курса и скоро­ сти при умеренной качке значения параметров этой корре­

ляционной функции равны в среднем:

в к ( м ) =

0,5°; а =

==0,75 ч - 1 ; Р = 5,3 ч - 1 .

 

 

В свою очередь, среднее значение ошибки курсоуказа;

ния постепенно изменяется с течением

времени

вследствие

влияния таких факторов, ка-к неравномерность износа под­ шипников роторов, изменение положения центра тяжести чувствительного элемента относительно центра величины и т. д. У исправного гирокомпаса такие изменения стано­ вятся ощутимыми только в течение длительного срока (ме­ сяцы, годы). Нет. оснований ожидать, что этот случайный процесс окажется стационарным. Поэтому удобнее харак­

теризовать его

не корреляционной, а структурной функ­

цией (§ 3.1).

Видимо,

с удовлетворительной

тччностью

можно аппроксимировать

ее выражением

 

 

• W ) = T l ' l .

( 2 - 6 )

где t — разность моментов времени, в которые рассматри­ ваются средние значения истинной ошибки измерения кур-

53

са корабля. Параметр у может оцениваться как -среднее значение коэффициента пропорциональности между про­ межутком времени от одного до другого определения сред­ него значения поправки компаса и квадратом разности оценок поправок компаса, к которым эти наблюдения при­ вели. Например, если определения поправки гирокомпаса примерно равноточны, то можно принять

2^(ДА'( --ДА',._02

 

 

7 = -

(-2-7)

/=1

Если ось 0—180° азимутального кольца репитера ги­ рокомпаса, установленного на пелорусе, с которого изме­ ряются пеленги при определении поправки компаса, со­ ставляет с диаметральной плоскостью корабля угол s, то исправление компасного курса поправкой компаса, опре­ деленной по пеленгам, ведет к постоянной остаточной ошибке измерения курса корабля, равной Е. Е С Л И не па­ раллельны одна другой осп 0—180° двух репитеров, уста­ новленных на разных пелорусах, причем пеленги, изме­ ренные по одному из них, исправляются поправкой компа­

са, определенной по другому, то возникает

систематиче­

ская ошибка измерения пеленга, равная углу

между ося­

ми 0—180° этих репитеров. Эксцентриситет оси вращения пеленгатора относительно центра картушки грубого от­ счета ведет к дополнительной систематической ошибке из­ мерения пеленгов, имеющей полукрутовой характер.

Описания поверок и регулировок, которые должны про­ водиться для уменьшения этих погрешностей, в инструк­ циях по эксплуатации гирокомпасов отсутствуют, а опу­ бликованные в периодической печати [71], [78] неполны. Поэтому такое описание приведено в § 3.7. Если поверки выполнены тщательно, средние квадратические величины амплитуд остаточных погрешностей составляют в сред­ нем 0,1—0,2°.

Отметим некоторые особенности применяемой ныне терминологии.

• Термином «истинный курс» в современном кораблево­ ждении обозначается несколько разных величин, в общем случае не. совпадающих: угол между плоскостью истин­ ного меридиана и диаметральной плоскостью корабля,

54

угол между плоскостью истинного меридиана и направле­

нием

перемещения корабля относительно

водной

среды

при

штиле, сумма

компасного курса и учитываемой по­

правки компаса.

Впредь будем применять

этот

термин

только

в первом из перечисленных значений, термин «ис­

тинный

пеленг» — только в значении,, указанном в начале

параграфа. Величины КК+АК

и КП + АП удобнее на­

зывать

оценками истинных курса

и пеленга.

Термин «поправка компаса» в современном кораблево­ ждении применяется для обозначения как учитываемой по­ правки компаса, так и. величины, противоположной истин­ ной ошибке измерения курса пли пеленга (например, в выражениях «штурман должен следить за постоянством поправки компаса», «обнаружив изменение поправки ком­ паса, надо ее определить» и т. д . ) . Если условиться для обозначения первого понятия применять только термин «учитываемая поправка компаса», то применение слово­ сочетания «поправка компаса» во втором его значении можно считать допустимым.

Поскольку истинная ошибка курсоуказания испыты­ вает непрерывные случайные отклонения относительно среднего значения, рекомендуется исправлять показания гирокомпаса средним значением его поправки. Оно. опре­ деляется, как правило, в базе. В те-ченне 1,5—3 часов через равные промежутки времени (обычно через 10 или 15 мин) определяются мгновенные значения поправки компаса.

Каждое

из них вычисляется как среднее арифметическое

из серии

(обычно пяти) разностей

эталонного

(принимае­

мого за

истинный) н измеренного

пеленгов

отдаленного

предмета или небесного светила. Среднее значение по­

правки компаса

вычисляется как

среднее

арифметиче­

ское из оцененных таким образом

мгновенных

ее

зна­

чений.

 

 

 

 

 

Обозначим символами: т—число

измерений

в серии;

п — число серий

измерений; т — промежуток

времени

ме­

жду средними моментами двух следующих одна за другой серий измерении.

Если курс корабля и пеленг на ориентир за время на­ блюдений существенно не изменяются, то ошибка за экс­ центриситет пеленгатора во всех измерениях будет прак­ тически одинаковой. Постоянными останутся также ошиб­ ка эталонного пеленга и ошибка, возникающая вследст­ вие непараллельности оси 0—180° репитера и диаметраль­ ной плоскости корабля.

55

Обозначим С И М В О Л О М

О 2 дисперсию суммы

этих по­

стоянных систематических

ошибок, символом о2

диспер­

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

но

независимы

(к ним относятся ошибки наведения пе­

ленгатора, отсчета по шкале репитера

и т. д . ), символом

о'

дисперсию

ошибкиопределения

среднего значения

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

чения

поправки гирокомпаса

равна

 

 

о2

=

а2 +

— а 2

+

а2 .

(2.8)

 

О.П

 

П

1

/ц/1

с

1

м*

\"'^J

Оценивая величину

о2 ,

 

будем

 

считать, что

истинная

ошибка

курсоуказания

в

течение

 

промежутка

времени,

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

ошибки курсоуказания от среднего, можно

считать

рав­

ной

 

 

п

 

 

соответственно

 

 

+ 2 / ? , [ ( я - 1 ) т ] } ,

(2.10)

или

 

 

где ki(n, т)коэффициент, зависящий от

числа п

серий

измерений, от промежутка времени х между средними мо­

ментами соседних серии и от параметров

корреляционной

функции (2.5). Величины

k\{n, т) и

У kx

(п., t)

при т =

= 10 мин, а = 0,75 ч - 1 ,

Р =

5 1 3 ч _ )

приведены в табл. 2.1.

Заметим, что основными

источниками

ошибок

в опре­

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

56

Т а б л и ц а 2.1

Значения коэффициентов

А, (/г,

т) и

j

/ ^ (,г,

т )

п р и ,

- ю мин;

 

 

 

в =

0.75

ч - 1 ;

р =

5,3

ч - 1

 

 

 

 

п

1

2

3

.4

5

6

 

8

10

13

16

19

 

1,00

0,78

0,55

0.35

0,20

0,12

0,07

0,06

0,04

0,03

0,03

 

1,00

0,88

0,7-1

0,59

0,45

0,34

0,27

0,25

0,20

0,17

0,16

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

величине

над этими ошибками и приняты должные меры

к предупреждению промахов в наблюдениях,

то увеличе­

ние числа

т измерений в каждой серии сверх

одного-двух

не имеет смысла: точности определения среднего значения поправки компаса оно практически не повысит.

Слагаемыми истинной остаточной ошибки результата однократного измерения курса (пеленга), исправленного средним значением поправки компаса, являются ошибка

определения

этой поправки, происшедшее

затем

измене­

ние среднего значения

истинной

ошибки

курсоуказания,

отклонение

мгновенного

значения

этой ошибки

от сред­

него и (при измерении пеленга) ошибки эксцентриситета, наведения пеленгатора, отсчета по шкале. Их можно счи­ тать взаимно независимыми случайными величинами. Если

после

определения

поправки компаса

прошло

более 1 —

1,5 ч,

то

взаимно

независимы и

отклонения

мгновенных

значений

ошибки

курсоуказания

от

среднего

во время

определения поправки компаса и в момент измерения, ис­ правляемого этой поправкой. При этом условии диспер- "сия истинной остаточной ошибки результата однократного измерения курса (пеленга), исправленного средней по­

правкой компаса,

равна

 

 

^ =

+ S « о ( 9 + < ^ > + ( < « + <#.

С 2 - 1 2 )

где

t — промежуток времени, прошедший после

опре­

 

деления

поправки компаса:

 

а2 п дисперсия ошибки за эксцентриситет пеленга­ тора.

57

При оценивании ошибок измерения курса два

послед­

них члена (в скобках) формулы (2.12) учитываться

не дол­

жны.

 

Исправляя результаты измерений средним значением поправки компаса, надо следить за ее постоянством, для чего сличать с результатами контрольных определений. Пусть Д/0—-учитываемое среднее значение поправки ком­

паса,

— результат ее

контрольного

определения, вы­

полненного

через промежуток

времени

t

после первого

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

Обозначим

ДЛ*2— А /С i = /.

 

Перед

нами возникают два

вопроса:

 

1. Не

свидетельствует

ли величина

•/ о

существенном

изменении среднего значения истинной ошибки курсоука­ зания? Если да, то определенной ранее средней поправке компаса доверять нельзя. При первой же возможности сле­ дует определить ее заново. В простейшем случае, когда ошибки величин АК\ и Д/Сг можно считать практически независимыми, а их распределения нормальными, для от­ вета на этот вопрос следует, пользуясь формулой (2.8), ап­

риорно оценить дисперсии

о| п j и п 2

ошибок опреде­

ления поправок АК\, ДЯа

и дисперсию

их разности:

а также сравнить величину / с величиной 2ad или 3od (в зависимости от того, какая доверительная вероятность по­ ложена в основу критерия аномально большого отклоне­ ния).

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

+

( 2 Л 4 )

 

где

 

ffo.nl

+ ^ко (О

c

iо. п2

 

Дисперсия

остаточной

ошибки

этой

оценки'

 

 

0 2

=

С ( » •

 

 

(2.16).

 

к

- с

Р\ + Рг

 

 

 

 

 

 

2. Какова

эффективная

оценка

мгновенного

значения

поправки компаса в момент времени t + x? Будем

считать,

58

что в момент времени t выполнены единичное определение поправки компаса или серия измерений, общая длитель­

ность

которой

невелика, так что мгновенное значение ис­

тинной

ошибки

курсоуказания за время измерений замет­

но измениться

не

могло.

 

 

 

 

Воспользуемся выражением линейной регрессии двух

случайных величин [И, § 14.8]. Пусть Мх, Му

математи­

ческие

ожидания

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

X и У; а 2 Д ) а\ — и х

дисперсии;

RXY

— корреляционный момент;

ГХУ~^ХУ'-

(ал-ау)коэффициент корреляции.

Если известно зна­

чение

у,

которое

в некотором

опыте

приняла

случайная

величина

Y,

то

эффективная априорная

оценка

соответст­

вующего

значения

случайной

величины

А' равна

 

х = Мх

+

 

(у - Му) =МХ

+ ^-

rXY

(у - My), (2.17)

 

 

 

G~y

 

°у

 

 

 

еедисперсия

 

a i =

М [ [ х -

xf\

=

(1 -

r\Y) 4 .

(2.1-8)

Примем

 

 

 

 

 

 

 

 

 

 

x =

AK(t

+ 'J;

у =

АК2;

МХ

=

МУ=ЬК»

(2.19)

^

=

«y = W

^

y =

^ i W ;

 

г* =

' (2.20)

Это приведет

к

оценке

 

 

 

 

 

 

 

\К (t + х) = (1 -

,;) Д/<с +

гАКъ

(2.21)

дисперсию ошибки которой можно приближенно считать

равной

-

^ = ° к . с + ( 1 - 0 ^ < » г - -

; (2-22)

Подобным образом могут решаться и вопросы опреде­ ления, учета и-корректировки поправки гироазимута. Од­

нако следует

иметь в виду, /что изменение его поправки

по времени

не является стационарной случайной функ­

цией. Но отклонение скорости ухода гироазимута от сред­ ней можно описывать как стационарный случайный про­ цесс.

59

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