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

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

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

После введения поправок результаты измерений оказы­ ваются 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

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