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

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

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

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

Добавлен: 23.10.2024

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

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

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

Тогда система из п уравнении поправок (1.26) может быть записана в виде одного матричного уравнения

AX-L

=

V,

(1.34)

а выражениям (1.25), (1.29) и

(1.31)

будут соответство­

вать матричные равенства:

 

 

 

Г = 5 п р

+

й

(1.35)

М®

= ^р

+ Х;

(1.36)

L =

AX+\.

(1.37)

Рассмотрим выражение (1.29). В левой

его

части

стоят

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

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

ожидание

любой

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

является случайной

вели­

чиной, не является таковой и величина ^ П р, которую при обработке наблюдений выбирают произвольно, исходя из

удобства вычислений). Следовательно, не

случайна

и ве­

личина

Xj, не

является

 

случайным

вектор ,Y=||.^||.

Учтем, кроме того, что не

является

случайной

и

матри­

ца А. Значит, из равенства

(1.37)

следует

важный

вывод

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

матрице

KL

вектора L свободных

чле­

нов

 

уравнений

поправок:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

KL

= K„

 

-

 

(1.38)

где

КА—корреляционная

 

матрица

вектора А.

 

 

Дальнейшие выводы основываются' на предположении,

что

выполнены

следующие

условия:

 

 

 

 

 

 

— математическое ожидание любой остаточной

ошиб­

ки

равно

нулю, т. е.

 

 

 

 

 

 

 

 

 

 

 

 

 

М ( Д ) =

[|М(Л; )|и =

О п 1 Г

 

 

(1.39)

матрица

КА

— неособенная

и

известна

хотя

бы с

точностью до

постоянного

множителя;

 

 

 

 

— ранг матрицы КА

не меньше числа

т искомых ве­

личин (для этого необходимо, чтобы

число п

измерений

или,

что

равнозначно, число

п уравнений

поправок,

со­

ставленных по их результатам, было не. меньше числа m искомых величин);

— ранг матрицы А коэффициентов при неизвестных в уравнениях поправок равен числу m искомых величин.

22


Если второе условие выполнено, то может быть най­

дена матрица

 

 

 

 

Р

= \\Ри.\\„п = ^

- К , у 1 ,

(1.40)

где

—дисперсия ошибки

измерения,

вес которого

принят

равным

единице.

 

 

К вектору $ оценок искомых величин

предъявляются

следующие требования:

 

 

— компоненты вектора должны быть несмещенными оценками искомых величин, т. е. должно выполняться ра­ венство

 

 

М$) = М{Ь)\

 

(1.41)

дисперсии оценок искомых

величин

(т. е. диагональ­

ные

элементы

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

матрицы

К~

вектора £)

должны быть минимальны (такие оценки

принято назы­

вать

эффективными).

 

 

 

До сих пор мы не ограничивали себя какими-либо пред­

положениями о

распределениях

вероятностей

рассматри­

ваемых нами ошибок. Все выводы этого и предыдущего параграфов справедливы при любых распределениях этих ошибок, если только их дисперсии конечны. Теперь же возникла необходимость в таком предположении, посколь­ ку способ удовлетворения второго из перечисленных тре­ бованийсущественно зависит от того, каковы распреде­ ления вероятностей ошибок hi. В теории способа наимень­ ших квадратов доказывается, что в частном (но имеющем наибольшее практическое значение) случае, когда эти рас­ пределения нормальны, второе требование удовлетворяет­

ся, если вектор X отыскивается из матричного

уравнения

ArPAX = ArPL*

(1.42)

Нетрудно показать, что при этом удовлетворяется и первое требование. Действительно,, из выражения (1.42) следует, что

X={ArPA)-lATPL,

(1.43)

23


откуда, учитывая равенства (1.37), (1.39), (1.35) и (1.36), получим

М (X) =

( Л Т А 4 ) - 1 Л Т Я Ж

(L) = (Ат РА)'1

Я

PAX

= X;

(1.44)

 

^ ( Г ) = 6 „ р +

ЛГ(^) = 6п р +

^

=

^ ( 6 ) .

' (1.45)

Из

выражений (3.32),

(1.43),

(1.38)

и

(1.40),

учиты­

вая также, что матрицы

Р, Р~\

А^ РА,

тРА)'1

и КА

симметрические, а вектор £пр не случаен, следует, что кор­ реляционная матрица К~ вектора % оценок искомых ве­ личин равна

К~ = К- =

ТРА)~1 АТРКА

[ ( / Г Л 4 ) - 1 АТР\Г'=

 

= [А'РАУ'А'РА

{£РА)-1РКЬ

=

Т А4)~"1 РР-1

=

 

= а ? 1 ) ( Д 7 М ) - ' .

 

(1.46)

Если в выражение (1.46) подставляется величина (дисперсия ошибки измерения, вес которого принят рав­ ным единице), оцененная при обработке некоторых преды­ дущих измерений, то формула (1.46) выражает действие, которое называется априорным (от латинского а priori — изначально, до опыта) оцениванием корреляционной ма­ трицы К~. Но может быть выполнено и апостериорное (от

латинского а posteriori — из последующего, после опыта)

оценивание этой матрицы. Для этого вычисляется

вели­

чина о 2 ] } апостериорная

оценка

дисперсии измерения,

вес которого принят равным единице:

 

а-

VTPV

 

(1.47)

 

 

Апостериорная оценка

матрицы

К~ считается

рав­

ной

 

 

 

^ г = а ? „ . ( Л т Р / 1 ) - 1 .

(1.48)

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

24


зультатах измерений. Если они невелики (величина мала), говорят о хорошем внутреннем согласии результа­ тов измерений, в противном случае — о плохом согласии.

Поэтому апостериорное оценивание величины о( 1 ) часто называют также оцениванием точности измерений по их внутреннему согласию. Его частным случаем является оценивание точности непосредственных измерений иско­ мой величины по отклонениям их результатов от среднего арифметического значения.

Выше мы говорили, что изучение функций /г (а, |3 . . . ) , т. е. -характера функциональных зависимостей системати­ ческих ошибок от параметров, характеризующих условия измерений, является необходимым условием повышения

точности измерений. Заметим, что единственную

возмож­

ность

экспериментального исследования этих

зависимо­

стей

представляет рассмотрение отклонений о,-. Далее в тех

случаях, когда функции /г (а, (3 ... ) изучаются по сличению

результатов- исследуемых измерений с результатами бо­ лее.точных (эталонных) измерений, разности этих резуль­ татов, условно именуемые наблюденными значениями ис­

тинных ошибок, по существу тоже

являются отклонения­

ми уравновешенных значений от

результатовизмерений

(у эталонных измерений эти отклонения условно можно считать равными нулю).

Формулы (1.40), (1.43), (1.46) — (1.48) представляют собой матричную запись алгоритма так называемого обоб­ щенного способа наименьших квадратов, при нормальном распределении ошибок измерений дающего эффективные оценки искомых величин как в тех случаях, когда эти ошибки являются взаимно независимыми случайными ве­ личинами, так и в тех, когда они взаимно зависимы. При его обосновании используется только часть условий Гаус­ саКолмогорова (условия нормальности распределений, равенства нулю математических ожиданий всех остаточ­ ных ошибок и конечности их дисперсий). Если, кроме того, полагается удовлетворенным и условие взаимной некор­ релированности двух любых остаточных ошибок (предпо­ лагается, что систематические ошибки отсутствуют, т. е. остаточные ошибки случайны и, следовательно, корреля­ ционная матрица АГД днагональна), то этот алгоритм принимает .вид своего частного случая — классического

25


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

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

Коль скоро матрица /Сд диагональна, ее обращение

существенно упрощается. Как следует из выражений (1.40) и (3.23), в этом случае матрица Р также диагональна (она называется матрицей весов уравнений поправок). Ее г'-й

диагональный

(отличный

от нуля)

элемент, находящийся

в i-й строке

и в г'-м столбце (он

называется в е с о м

f-ro уравнения

поправок),

равен

 

 

 

 

(1.49)

Значительно облегчается и составление системы урав­ нений ("1.42). Взамен ранее применявшихся введем новые обозначения коэффициентов при искомых величинах в уравнениях поправок:

 

an

= al;

ai2 = bl;

alm = fii;

 

(1.50)

обозначим также i-ю

контрольную сумму

 

 

 

 

 

a, +

b,

+

... + h l - l l = s l .

 

(1,51)

Воспользуемся обозначением сумм по Гауссу:

 

 

[раа] = р1а\ +

р2а\

+

... + рр]

+ ...

+

 

 

 

[рад] = p1a1bl

+ р2а2Ь2

+

... +

piaibi

+

 

...+рпапЬп;

 

 

. . . .

.

 

 

 

 

 

 

 

 

 

[pah] = pialhi

+ p2a2h2

+

... +

рр^ц

+

... +

pnanhn\

 

 

[pal]

=

+ p2a2l2

+

... +

platlt

+

... +

pnajn;

[pas] = А « 1 5 1

+ / W

2

+

-

 

+

... +

pnans,;,

K '

'

[pbb] = pxb\ +

p2b\

+

... H-Plb\

+ ... +

pjfy

 

 

 

[pbh]

= plbihx

+ p2b2h2

+

... +

plbihl

+

... -b pabahn\

 

 

26