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