
книги из ГПНТБ / Скворцов, М. И. Теория и практика решения задач кораблевождения с учетом влияния систематических ошибок учебник
.pdfПосле введения поправок результаты измерений оказы ваются t отягощенными остаточными систематическими ошибками. Величина r-й остаточной систематической ошиб ки в результате /-го измерения равна
* , , = (С, - С г Э /,(а„ |
& . . . ) = *,/,(«/, |
h ...)• |
(1-3) |
|
Будем |
называть величину гг = Сг — Сг |
амплитудой г-й |
||
остаточной |
систематической |
ошибки. Она |
является |
слу |
чайной величиной уже потому, что оценка £г как функция результатов некоторых предыдущих измерений, неизбеж но отягощенных случайными ошибками, есть случайная ве личина. Естественным является предъявить требование,
чтобы величина £г была несмещенной1 |
оценкой |
амплиту |
ды Сг компенсируемой систематической |
ошибки. В |
обычной |
практике назначения поправок, служащих для компенса ции систематических ошибок, это требование, как правило, удовлетворяется. Для отыскания поправок применяется способ наименьших квадратов, а он, как известно, дает несмещенные оценки искомых величин. Но тогда, как не
посредственно |
следует из выражения (1.3), |
математиче |
||
ские ожидания |
амплитуды |
остаточной |
систематической |
|
ошибки и величины, которую эта ошибка |
примет в любом |
|||
измерении, следует полагать |
равными |
нулю. |
Дисперсия |
же амплитуды остаточной систематической ошибки в об щем случае отлична от нуля, поскольку наблюдения, из
которых |
определяется |
величина |
Сг, неизбежно |
отягощены |
||||
некоторыми |
ошибками. |
|
|
|
|
|
||
Введем |
обозначения: |
|
|
|
|
|
||
|
а |
— число |
измерений |
в |
рассматриваемой |
сово |
||
|
|
купности измерений; |
|
|
|
|||
Zr = \\zir\\nl—вектор |
значений |
г-и остаточной системати |
||||||
|
|
ческой |
ошибки в |
этих |
измерениях; |
|
||
Fr = ЦД. Dm — вектор |
значений |
|
/;г, |
которые |
функция |
|||
|
|
/г (а, {3 ... ) приняла |
в этих измерениях. |
|||||
Тогда, как следует из формул |
(1.3), (3.10) |
и из |
опре |
|||||
деления |
(3.29) корреляционной |
матрицы, вектор ZT |
и его |
1 Несмещенной называется такая оценка искомой величины, мате матическое ожидание которой равно истинному значению искомой величины.
10
корреляционная матрица |
Kz |
могут быть |
записаны в |
виде |
|
|
|
Zr |
= |
zrFry |
(1.4) |
Kz=M(ZrZj) |
|
= oiFrFj, |
(1.5) |
где о2—дисперсия амплитуды zT рассматриваемой оста точной систематической ошибки.
В некоторых случаях для вычисления элементов, ма трицы Кг удобно пользоваться правилом, сформулирован- 'ным В. Ф. Лукьяновым [49]: корреляционный момент двух случайных величин равен сумме дисперсий их общих сла гаемых.
Из выражения (1.5) видно, что коэффициент корреля ции между величинами, которые принимает остаточная систематическая ошибка в двух любых измерениях одной совокупности, равен единице ( + 1 или — 1) . В этом со стоит наиболее существен«ое отлимие систематических оши бок от случайных (коэффициент корреляции между слу чайными ошибками двух любых измерений одной сово купности, как видно из приведенного выше определения случайных ошибок, равен нулю). Расположив элементар ные ошибки измерений в порядке возрастания коэффи циента корреляции между их величинами в двух любых измерениях одной совокупности, мы увидим, что случай ные и систематические ошибки являются предельными чле
нами |
этого ряда, |
соответствующими |
предельным |
значе |
||
ниям |
коэффициента |
корреляции (нулю |
или |
единице).. Те |
||
ошибки измерений, |
которым |
соответствуют |
промежуточ |
|||
ные |
значения коэффициента |
корреляции, отличные |
и от |
нуля, и от единицы, принято называть зависимыми. Чет ких границ, которые отделяли бы одни ошибки от других, не существует: случайные и систематические ошибки яв ляются предельными частными случаями зависимых.
Говорить о том, что та или иная ошибка является слу чайной, зависимой или систематической имеет смысл толь ко тогда, когда одновременно указывается, в какой со вокупности измерений рассматриваются свойства этой ошибки. Одна и та же ошибка в одной совокупности изме
рений может |
проявлять свойства систематической, а в дру |
|
гой — случайной. Это обстоятельство |
отмечалось еще Гаус |
|
сом [15, стр. |
17]. Приведем пример. |
Если корабль выпу |
скает ракеты по цели залпом, то ошибка поправки |
систе |
мы курсоуказания вызывает одинаковые отклонения |
всех |
11
ракет от цели и выступает как систематическая ошибка стрельбы. Если несколько кораблей независимо один от другого ведут стрельбу по одной цели, то ошибки курсоуказания каждого из них ведут к различным, взаимно не зависимым отклонениям точек падения ракет от цели и выступают как случайные ошибки.
Нельзя отождествлять термины «случайная ошибка» и , «ошибка измерений, являющаяся случайной величиной». Случайными называются не всякие ошибки, являющиеся Случайными величинами, а лишь те из них, которые в рас сматриваемой совокупности измерений взаимно независи мы (некоррелированы). Систематические ошибки — слу чайные величины, хотя случайными ошибками не являют
ся. Применение |
к |
их |
изучению математического |
аппарата |
теории вероятностей |
и математической статистики |
не толь |
||
кодопустимо, |
но |
и |
необходимо. |
|
Можно предположить, |Что встречающееся во многих по собиях по теории ошибок противопоставление случайных ошибок, как обладающих свойствами случайных величин, систематическим ошибкам, как якобы этими свойствами не обладающим, явилось следствием несовершенства термино логии, приводившего к невольной подмене понятий (ото ждествлению понятий «случайная ошибка» и «ошибка, яв ляющаяся случайной величиной»). С этой точки зрения термины, которыми пользовался Гаусс, были более удач ными. Он подразделял ошибки измерений на правильные (regulare, regelmaJHge) и неправильные (irregulare, unregelniafige). Однако мы будем и впредь применять ставшие уже привычными современные термины.
Каноническое представление ошибок измерений. Из выражений (1.1) и (1.3) видно, что систематические ошиб ки измерений являются элементарными случайными функ циями (§ 3.1) параметров, характеризующих условия из мерений. Сопоставим это обстоятельство с гипотезой об аддитивной структуре ошибок измерений [27, стр. 271—274], [32, стр. 9], [50, стр. 116]. Согласно этой гипотезе ошибки измерений представляют собой суммы большого числа элементарных ошибок разного происхождения, чем объ ясняется близость распределений реально наблюдаемых ошибок к нормальному распределению. Обратим внима ние на другой аспект проблемы. Обозначим истинную ошибку г'-го измерения
д; = £ / ; - £ / , „ „ , |
' |
(1.6) |
12
где |
U] — отсчет прибора |
(инструмента) |
в/-м |
измерении; |
Ui |
ист — истинное значение измеряемой |
величины ,в мо |
||
|
мент i'-го измерения. |
|
|
|
Будем рассматривать |
величину А'. |
как |
случайную |
функцию параметров, характеризующих условия измере ний. В соответствии с теоремой о каноническом представ лении случайной функции [61, гл. 6], [11, гл. 16] она может быть сколь угодно большим числом способов представлена
в виде суммы эл-ементарных случайных |
функции |
(т. |
е. |
||
систематических ошибок): |
|
|
|
|
|
|
|
s |
|
|
|
А;(в / , Р / . - . . ) = |
2 ^ ( я « . Р / - ) |
+ 8/. |
С1 -7 ) |
||
где 8г — остаточный |
член |
канонического |
разложения. |
|
|
Если дисперсия |
случайной функции |
конечна, |
то, |
вы |
брав достаточно большое число s, можно сколь угодно большим числом способов представить ее в виде (1.7) та ким образом, чтобы дисперсия остаточного члена oi была меньше любого наперед заданного положительного числа. Это приводит к следующей формулировке утверждения об аддитивной структуре ошибок измерений. Если дисперсия истинной ошибки в рассматриваемой совокупности изме рений конечна, то с любой наперед заданной точностью эта ошибка может быть представлена в виде конечного числа систематических ошибок.
Обычно нас удовлетворяют не любые представления ошибок измерений в виде (1.7), а лишь те из них, которые
отвечают некоторым |
дополнительным требованиям. |
Пусть |
|||
Сг—значение |
коэффициента £г в |
разложении |
(1.7) |
истин |
|
ных ошибок |
одной |
совокупности |
измерений, |
Сг —значе |
ние того же коэффициента в разложении ошибок другой совокупности измерений, отделенной от первой некоторым промежутком времени. Если можно приближенно считать, что справедливо равенство
с; = с;, |
(1.8)- |
причем дисперсия ошибки, с какой оно удовлетворяется, доступна оцениванию, то это означает принципиальную возможность, найдя из результатов первых измерений
оценку Сг этого коэффициента, компенсировать r-ю систе-
13
матическую |
ошибку |
в измерениях |
второй совокупности |
введением |
поправок |
вида (1.2). |
|
Отыскание таких |
функций fr(a, |
(3 . . . ) , чтобы дисперсии |
ошибок, с какими удовлетворяются равенства (1.8), были минимальными, составляет одну из важнейших задач тео рии и практики измерений. Практически оправдано введе ние в результаты измерений лишь тех поправок, в отно шении которых с удовлетворяющей нас точностью требо
вание (1.8) |
можно считать выполненным. Но, ограничи |
вая таким |
образом число членов канонического разложе |
ния (1.7), мы уже не можем считать остаточный член Зг пренебрежимо малым.
Учитывая обозначения (1.1) — (1.3) и (1.7), можно представить в виде канонического разложения и истинную
остаточную ошибку |
Ас |
|
|
|
|
|
д, = u t - |
и , |
,(СТ =* (£/; + |
2 |
- и , |
и с т = |
|
|
= |
2 * г / г ( а / . & . . 0 |
+ 8„ |
|
(1.9) |
|
|
|
/=1 |
|
|
|
|
где Ui—результат |
г-го измерения, |
исправленный |
всеми |
|||
учитываемыми |
поправками. |
|
|
|
|
|
Величину о,- остаточного члена |
разложения |
(1.9) |
мож |
но было бы назвать остаточной ошибкой. Но этот термин уже применяется для обозначения другого понятия. По этому будем называть ее пост-остаточной ошибкой /-го
измерения. |
|
|
|
|
|
|
|
Введем обозначения: |
|
|
|
|
|
||
о = || о\ ||я] |
— вектор |
пост-остаточных |
ошибок |
совокуп |
|||
|
|
ности |
измерений; |
|
|
|
|
F—\\fir\\„s— |
|
блочная матрица, |
столбцами которой яв |
||||
|
|
ляются |
векторы |
/•'[, F2, |
Fr, |
Ft\ |
|
Z—\\zr\\s\ |
— вектор |
амплитуд |
zu |
z2, |
zr, |
zs оста |
|
|
|
точных |
систематических |
ошибок. |
|
||
Тогда вектор истинных остаточных ошибок совокупно |
|||||||
сти измерений |
может быть записан в виде |
|
|
||||
|
|
Д = 11-Ми = ^ |
|
+ 8. |
|
(1Л0) |
Будем считать, что при любых г и i случайные вели чины zT и S; взаимно независимы. Тогда в соответствии с
14
выражением (3.32) корреляционная матрица вектора Д равна
|
К, |
= М ( М т ) = FKZFT |
+ Kv |
(1.11) |
|
где /fz = |
Ж ( Z Z T ) — корреляционная |
матрица |
вектора |
Z |
|
|
|
амплитуд остаточных систематических |
|||
|
|
ошибок; |
|
|
|
Къ= |
М (ойт ; — корреляционная |
матрица |
вектора |
8 |
|
|
|
пост-остаточных ошибок совокупности |
|||
|
|
измерении. |
|
|
|
Если корреляция между величинами 8i в двух любых |
|||||
измерениях одной |
совокупности отсутствует (матрица |
К& |
днагональна), то их можно называть случайными ошиб ками измерений, в противном случае — зависимыми ошиб ками. Но в обоих случаях величина В; является случайной функцией параметров, характеризующих условия измере ний, и-.также может быть представлена в виде (Г.7). Та ким образом, называть случайные и зависимые ошибки
элементарными можно только |
условно, подобно тому, как |
в физике объекты, состоящие |
из протонов, нейтронов и |
электронов, условно называются атомами (неделимыми). Столь же условно и выражение «измерения, свободные от систематических ошибок». Его можно понимать только в переносном смысле для обозначения предположения, что компоненты векгора Д в рассматриваемой совокупности из мерении обладают свойствами, мало отличающимися от свойств случайных ошибок.
Сложная структура случайных ошибок любых измере ний проявляется ощутимо каждый раз, как обнаруживает ся новый, дотоле неизвестный источник систематических ошибок. Тогда оказывается, что вновь открытые система тические ошибки являются составными частями тех оши бок, которые ранее считались случайными, т. е. элемен тарными. Этот процесс выделения систематических оши бок из случайных является одной из необходимых предпо сылок повышения точности измерений. Чем больших успе хов мы добиваемся на этом пути, тем большую практи ческую значимость приобретает изучение свойств систе матических ошибок.
Влияние систематических ошибок измерений на точ
ность оценок искомых величин. Пусть вектор X = || jc |j„a
15
оценок искомых величин отыскивается как произведение матриц:
|
|
, |
. |
X |
= |
GL, |
|
|
|
|
(1.12) |
|
где |
G = |
Gmn |
—матрица |
линейного |
преобразования; |
|
||||||
|
^ — |
II h |
llm — случайный |
вектор. |
|
|
|
|
|
|||
|
Предположим, что корреляционная матрица вектора L |
|||||||||||
идентична |
корреляционной |
|
матрице |
ЛГД |
вектора |
А = |
||||||
= |
II d < II ni |
истинных |
остаточных |
ошибок |
совокупности |
из |
||||||
мерений. Тогда, как следует из правила |
(3.32) |
отыскания |
||||||||||
корреляционной матрицы |
линейной |
векторной |
функции |
|||||||||
случайного |
вектора, |
корреляционная |
матрица |
К~ |
|
век- |
||||||
тора /Y равна |
|
|
|
|
|
|
х |
|
|
|||
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
К~ = |
Gl<fi\ |
|
|
|
(1.13) |
|
||
|
|
|
|
х |
|
|
|
|
|
|
|
|
|
Для |
оценивания |
корреляционной |
матрицы |
|
суще |
||||||
ствуют два |
пути. |
|
|
|
|
|
|
|
|
|
||
|
Первый |
путь — статистическое |
исследование |
истинных |
||||||||
ошибок |
измерений, непосредственно |
наблюденных |
в |
спе |
циально организованном эксперименте. Истинная ошибка каждого измерения оценивается как разность его резуль тата и эталонного значения измеряемой величины (ре зультата более точного измерения), которое принимается за истинное. Получив таким образом ряд наблюденных значений истинных ошибок, вычисляют их корреляцион ные моменты. Это осуществляется способами, общеприня тыми в математической статистике и описанными во мно гих пособиях, например [11, § 14.6], [39, § 10]. Считая, что корреляционные матрицы однотипных измерений в сход ных условиях примерно одинаковы, распространяют затем
результаты такого исследования на те |
измерения, |
кото |
||
рые предстоит выполнять в будущем. |
|
|
||
Второй |
путь — косвенное |
оценивание |
корреляционной |
|
матрицы /с~д основанное на |
представлении истинной оста |
|||
точной ошибки каждого измерения в виде (1.9) |
суммы |
|||
остаточных |
систематических |
ошибок и |
пост-остаточной |
ошибки, что приводит к выражению (1.11). Этот путь несравненно более прост и удобен, поскольку не требует проведения специальных сложных экспериментов и экс траполяции их результатов на все измерения, какие пред стоит выполнять в будущем. Чтобы получить возможность
16
идти этим путем, надо уметь оценивать корреляционные матрицы Kz и /(•.. Ниже будет описан способ обработки результатов измерений, который позволяет получать оцен ки не только искомых величин, но и корреляционной ма трицы Кг- Что же касается корреляционной матрицы Кй, то в нередко встречающихся случаях, когда ее можно по лагать диагональной, она может оцениваться по разно стям уравненных и измеренных значений измеряемых ве личин, как это принято при классическом применении спо
соба наименьших |
квадратов.. |
|
|
|
||
Если матрицы |
Kz |
и/7<"5 известны, то, как |
видно |
из |
||
выражений (1.11) |
и |
(1.13), |
|
|
|
|
К~ |
= |
GFKZ |
(GF)T + GKfiT- |
(1.14) |
||
Если рассматривается влияние только одной, г-н систе |
||||||
матической ошибки, то каноническое разложение |
(1.9) |
ис |
||||
тинной остаточной ошибки |
Д; можно |
представить в виде |
||||
|
Л/ = г,Л(«/ . Р/ - ) + |
8|. |
(1Л5) |
т. е. считать, все остальные систематические ошибки со ставными частями пост-остаточной ошибки. Тогда равен ства (1.11) и (1.14) примут вид
|
К~ = °*rGF, |
(GFry |
+ |
GK&\ |
(1.17) |
|
где of — дисперсия |
амплитуды |
r-й |
остаточной |
системати |
||
ческой ошибки. |
|
|
|
|
|
|
Введем |
обозначение |
|
|
|
|
|
|
|
К~ = |
a-rGFT |
(GFr)r. |
|
(1.18) |
|
|
X ( r ) |
|
|
|
|
Матрица |
АГ~ |
есть |
корреляционная матрица вектора |
Х(г)
ошибок в оценках искомых величин, обусловленных влия нием r-й остаточной систематической ошибки измерений. Поскольку ранг матрицы Fr равен единице, ясно, что ранг матрицы К~ также равен единице, т. е. ошибки оценок
искомых величин, происходящие от влияния рассматри ваемой r-й систематической ошибки измерений, линейно зависимы. По аналогии с терминами, применяемыми в от ношении двухмерной случайной величины, можно сказать,
17
что |
если ошибки оценок искомых величин, |
происходящие |
от |
влияния случайных ошибок измерений, |
характеризуют |
ся /71-мерным средним квадратическим эллипсоидом оши бок (или, что то же самое, системой из т взаимно неза висимых векториальных ошибок), то ошибки оценок иско
мых |
величин, |
происходящие от |
влияния |
г-\\ систематиче |
|||
ской |
ошибки |
измерений,. характеризуются |
одной |
векто |
|||
риальной ошибкой системы оценок искомых величин. |
|||||||
Сравнение |
выражений |
(1.18) |
и (1.12) |
приводит к сле |
|||
дующему правилу оценивания |
корреляционной |
матрицы |
|||||
/<"- . |
Пусть |
А - = II Л , II л1 —вектор |
значений |
/,> — |
|||
А'(г) |
|
|
|
|
|
|
|
=/г ( а л Р* • • •)> которые |
приняла функция |
/,-(а, р ... ) в из |
мерениях рассматриваемой совокупности. Подвергнув его тому же преобразованию, каким из вектора L отыскивает
ся вектор |
X оценок искомых |
величин, |
найдем вектор |
|||
|
|
.Cr = |
||c,|U = C/v |
|
(1.19) |
|
Тогда искомая |
матрица |
К~ |
будет равна |
|||
|
|
К~ |
=°'rCrCj. |
|
(1.20) |
|
|
|
Х(г) |
|
|
|
|
Величины |
3rCj, |
агсъ ..., |
orCj, |
. . . , <згст |
представляют собой |
|
1-ю, 2-ю, |
|
/n-ю компоненты искомой |
векториальной |
|||
ошибки. |
Второе решение |
составляет |
система величин |
{—агС}) •
*
§1.2. СПОСОБ НАИМЕНЬШИХ КВАДРАТОВ И АЛГОРИТМ
ПО С Л Е Д О В А Т Е Л Ь Н О Г О УТОЧНЕНИЯ ОЦЕНОК ИСКОМЫХ
ВЕЛИЧИН
Способ наименьших квадратов. При обработке наблю дений часто встречаются случаи, когда измеряются не сами искомые величины непосредственно, а некоторые другие величины, являющиеся функциями искомых. На пример, при определении места корабля искомыми явля ются его географические координаты, измеряются же на вигационные параметры — пеленги, расстояния до ориен тиров и т. д., зависящие от места корабля на земной по верхности и являющиеся функциями географических ко ординат.
Пусть зависимость между истинными значениями £ ь
. . . , 5ji • • •> 5m искомых величин и истинным значением
18
fr'tmvr |
измеряемой |
величины |
аппроксимируется |
(прибли |
|||
женно выражается) |
некоторой функцией ф{(£ь |
£j, |
|||||
• • •> |
$ m) ^ |
|
|
|
|
|
|
|
|
|
iJt |
. . , |
g = t / , „ C I + A; , |
|
(1.21) |
где |
Д(° — истинная |
ошибка |
аппроксимации. |
|
|
||
С другой стороны, |
в |
соответствии с формулой |
|
(1.9) за |
висимость между исправленным всеми учитываемыми по правками результатом 1-го измерения С,- и истинным зна
чением измеряемой |
величины выражается |
уравнением |
|||
|
|
Ut = UlH„+A'r |
' |
(1.22) |
|
где А.' — истинная ошибка |
измерения. |
|
|
||
Таким |
образом, |
|
|
|
|
Как бы |
ни были |
точны |
измерения, они |
всегда |
отяго |
щены неизбежными ошибками. Найти из результатов из мерений истинные значения искомых величин невозможно,
поэтому мы довольствуемся |
отысканием |
приближенных |
||||
значений — оценок |
искомых |
величин, |
стремясь |
к тому, |
||
чтобы |
их отличия |
от истинных значений были минималь |
||||
ными. |
Зависимость |
между оценками |
£l t |
. . . , |
Хт |
искомых величин и результатом /-го измерения £Д- выра
жается |
уравнением |
|
|
|
|
|
|
|||
|
|
|
*?ЛхЛъ |
. . . , i , . . . , f J - ^ - = ^ . |
|
(1.24) |
||||
|
Величина |
<Ь (£}>•••> |
£,•>•••> |
£ т ) |
называется |
у р а в |
||||
н е н н ы м |
(уравновешенным) |
значением |
измеряемой вели |
|||||||
чины, а |
величина |
ог- — о т к л о н е н и е м |
уравновешенного |
|||||||
значения измеряемой величины от результата |
измерения. |
|||||||||
|
Поскольку в общем случае функции |
ф;(...) |
не |
являют |
||||||
ся |
линейными, непосредственно |
пользоваться |
уравнения |
|||||||
ми |
(1.23) и |
(1.24) |
было |
бы |
неудобно. Поэтому |
каждую |
искомую величину Sj выражают как сумму ее произвольного
приближенного значения $jnp и искомой |
поправки ху.' |
^ = ^ п Р + А- ; . |
(1.25) |
19