Файл: Скворцов, М. И. Теория и практика решения задач кораблевождения с учетом влияния систематических ошибок учебник.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.10.2024

Просмотров: 77

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

Итак, в уравнения поправок в качестве неизвестных помимо поправок к приближенным значениям «основных»

искомых величин включаются оценки —Z\, —z2, ..., —zr,

..., —zs поправок к произвольным приближенным значе­

ниям амплитуд поправок,- предназначенных для

компен­

сации

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

ошибок

измерений.

Если

хи ...,

л-,-,

...., Хт — искомые оценки

поправок

к

приближенным

значениям «основных»

искомых

величин,

то, как

следует

из

выражений

(1.2),

(1.3) и

(1.26), t'-e

уравнение

попра­

вок

принимает

вид

 

 

 

 

 

 

 

 

 

 

Я / i * i

+

••• +

aijxl

+

... + almxm

 

+

 

 

 

+

а 1 (m-f 1)Хт+\^~

 

••• а 1 {т+г)Хт+г

+ ••• +

U i

 

(m+s)Xm+s

 

 

 

 

 

 

 

=

 

 

 

 

(1-95)

 

 

 

 

 

Я»+, =

 

 

 

 

(1-97)

 

Последующее

решение осуществляется

в

соответствии

с обычными предписаниями способа наименьших квадра­

тов. Как и прежде,

введем

обозначения;

 

 

Snp = ( 1 п р ~

вектор

произвольных

приближенных

 

значений

«основных»

искомых

вели­

 

чин;

 

 

 

 

^ = II-X/ILi вектор

оценок поправок к произволь­

 

ным приближенным

значениям

«ос­

 

новных»

искомых величин;

 

L — || lL ||л1

— вектор

свободных

членов уравнений

 

поправок;

 

 

 

40


А — || a,j \\пт — матрица

коэффициентов

при

«основ­

 

 

 

ных» неизвестных в уравнениях по­

 

 

 

правок;

 

 

 

 

 

 

 

 

Fr—\\fir\\n\ — вектор

значений

/,>, которые функция

 

 

 

fr(a,

 

р ... )

приняла в

1-м, 2-м, ... , i-м,

 

 

 

... , /i-м измерениях.

 

 

 

 

Кроме

того,

введем

следующие

обозначения:-

 

 

 

I hp II ji — вектор

п р о и з в о л ь н ы х п р и б л и ж е н н ы х

 

 

 

з н а ч е н и й

а м п л и т у д

поправок,

вклю ­

 

 

 

ченных в число и с к о м ы х

величин;

 

Х7

х m-j-r

 

вектор

искомых

оценок

поправок

к

 

 

 

 

 

произвольным

приближенным

значе­

 

 

 

ниям этих амплитуд;

 

 

 

 

К6 = М(ЪЪТ)— корреляционная

матрица

вектора

8 =

 

 

 

— II ^

II «1

 

остаточных

ошибок,

не'

 

 

 

включенных

в

число

искомых'

вели­

 

 

- j

чин;

 

 

 

 

 

 

 

 

 

Р 2 =

 

матрица

весов

уравнении

поправок

 

 

Обозначим

блочные

матрицы:

 

 

 

 

 

 

 

 

^ =

I I / „ I L

= II

 

Ft...Fr ЛИ;

 

(1.98)

 

 

 

^2 пр

=

'•пр

 

 

 

 

(1.99)

 

 

 

 

пр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

(1.100)

 

 

 

 

 

 

 

х 7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A.2

=

\\A\-F\\.

 

 

 

(1.101)

Тогда

решение

(1.42) — (1.48)

примет

вид

 

 

 

 

 

 

А\Р2А2Х2

=

А\Р2Ц

 

 

(1.102)

 

 

 

 

 

 

 

 

 

 

Х2 =

 

{А\Р2А2ул

 

А\Р21\

 

 

(1.103)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^2 пр +

Х2\

 

 

 

(1.104)

 

 

 

KiB

=

 

o2{l)(AlP2A2r\

 

 

(1.105)

Отметим наиболее существенные особенности этого ре­ шения.

41


1. При вычислении свободных членов уравнении попра­ вок мы задаемся произвольными приближенными значе­ ниями всех искомых величин, в том числе и амплитуд по­ правок, предназначенных для компенсации систематиче­ ских ошибок, включенных в число искомых величин. С точ­ ностью до величин второго порядка малости вектор 5в оценок искомых величин не зависит от тех произвольных приближенных значений, которые при вычислениях при­ даны компонентам вектора ij2 np. Информация о значениях амплитуд поправок, предназначенных для компенсации си­ стематических ошибок, которые были найдены при преды­ дущих наблюдениях, нами не используется. Она нам не нужна. Таким образом, применение способа В равносиль­ но обычно не оговариваемому предположению, что все на­ блюдения, выполненные ранее для определения этих ам­ плитуд, считаются не заслуживающими доверия.

 

2. Любая из остаточных систематических ошибок, ам­

плитуда которой включена в

число искомых величин,

как

бы

велика она ни была, не оказывает

никакого влияния

на

точность оценок искомых

величин:

все элементы

кор­

реляционной матрицы К~ , которой согласно формулам

(1.18) и (1.20) характеризуются ошибки вектора оценок искомых величин, происходящие от влияния этой система­ тической ошибки, оказываются равными нулю. Действи­ тельно, если применить общее правило (1.20) к оценива­

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

матрицы

К~ ,

получим

следующие

уравнения поправок

и нормальные

уравнения:

 

 

r

— F-, =

V*;

 

(1.106)

 

AlP2A2C/~A\P,Fr.

 

 

(1.107)

Приняв во внимание

обозначение

(1.101) и правило (3.16),

получим

 

 

 

 

 

 

 

С г

- 0

м

у

 

(1.108)

Следовательно,

 

 

 

 

 

 

^

"

^

 

^

w

( 1 Л 0 Э )

3. Включение амплитуды

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

ошибки в

число искомых величин может привести к тому, что систе­ ма нормальных уравнений (1.102) окажется вовсе це

42


имеющей единственного конечного решения. Это происхо­ дит, например, при попытке исключить постоянную систе­ матическую ошибку пеленгования при определении места корабля по пеленгам нескольких (более двух) ориентиров в случае неопределенности, когда все ориентиры и ко­ рабль находятся на одной окружности (§ 3.6). Нетрудно убедиться, что подобные ситуации неопределенности возни­ кают в тех, и только тех случаях, когда рассматриваемая систематическая ошибка до ее включения в число искомых величин не могла вызывать несогласия результатов изме­ рений (т. е. когда любое изменение амплитуды этой систе­ матической ошибки не ведет к изменению отклонений Vi результатов измерений от их уравновешенных значений).

Чтобы доказать правильность этого утверждения, об­

ратимся

к

следующей

теореме

высшей

алгебры

[46,

стр 153].

Пусть

X — || х- || , л 1 — в е к т о р неизвестных

вели­

чин;

А =

|| я , - 1 | п т — м а т р и ц а коэффициентов при

неизве­

стных.

 

 

 

 

 

 

 

 

 

Система

однородных

линейных

уравнений

 

 

имеет

систему

решений,

отличных

от нуля,

тогда,

и толь­

ко тогда, когда ранг матрицы А меньше числа т неизве­ стных величин.

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

мальных уравнений

 

Л Т Р Л ^ = ЛТ Я/.

(1.110)

 

V

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

величин), причем систематическая ошибка zir =

zrfr [ah (3,....)

не могла вызывать

несогласия

результатов

измерений,

т. е. удовлетворялось

матричное

равенство

 

 

АХ + zrFr

=

0.

(1-111)

Согласно приведенной теореме это возможно только тогда, когда ранг блочной матрицы 'Л 2 = |1/4|/г;.|| мень­ ше числа т-\-\, что, в свою очередь, влечет за собой ра­ венство нулю главного определителя системы нормальных уравнений (1.102) и отсутствие конечного ее решения.

43


Справедливо и обратное утверждение:

если . матри­

ца Р2 — неособенная' и ее ранг не

меньше числа искомых

величин, а система' нормальных

уравнений

(1.102) не

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

Действительно,

если п^.т-\-1

и ранг

матрицы Р 2

боль­

ше

числа т + 1, то система уравнений

(1.102)

не

может

иметь конечного

решения

только

тогда,

когда

ранг

блоч­

ной

матрицы /1о меньше

ш + 1, что, в свою

очередь,

озна­

чает

возможность

удовлетворения

условия

(1.111).

 

Поскольку влияние случайных ошибок на точность оце­ нок искомых величин может быть уменьшено увеличением числа измерений с последующим осреднением их резуль­ татов, а также принимая во внимание, что способ В часто сводится к простым и удобным графическим построениям, ему нередко отдают предпочтение. Количественную харак­ теристику относительных достоинств способов А и В дает сравнение корреляционных матриц векторов ошибок до­ ставляемых ими оценок искомых величин. Эти корреля­ ционные матрицы можно оценить, пользуясь общими вы­ ражениями (1.94) и (1.105).

Следует иметь в виду, что в общем случае возможно найти такие оценки искомых величин, которые будут обла­ дать дисперсиями меньшими, нежели оценки, доставляе­ мые любым из рассмотренных способов А и В. Действи­ тельно, пусть $д — вектор оценок «основных» искомых ве­ личин, доставляемых способом A; i j B — вектор оценок тех же искомых величин, доставляемых способом В, если в число искомых величин дополнительно включена ампли-

• туда г-й остаточной систематической ошибки. При 5=5А минимальны дисперсии ошибок в оценках искомых вели­ чин, обусловленных влиянием случайных ошибок измере­ ний; при $ = ?в минимальны (равны нулю) дисперсии оши­ бок в оценках искомых величин, обусловленных влиянием рассматриваемой систематической ошибки. Поскольку эти дисперсии являются непрерывными однозначными функ­ циями оценок искомых величин, должны существовать оценки искомых величин, лежащие между оценками, до­ ставляемыми способами А и В, дисперсии ошибок кото-

44