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