
книги из ГПНТБ / Болдырев, В. С. Методы математической статистики в гидрографии и кораблевождении учеб. пособие
.pdf/w2= х г+ у 2 ;
А?2“ г 2(А2COS2 + ß ZSir?2<р) j
Al2=A2tß2;
A 2+ ß2=-v2(A 2co s2cf ) + ß 2s i n 2(f >') .
Разделим |
на |
>1 |
|
|
|
|
|
|
||
|
/ |
в * |
- |
|
|
+ —у sin |
|
|
||
|
/ f |
A * |
|
|
А |
|
|
|
||
|
|
д2 |
|
|
|
|
І+ К‘ |
-7 -5 — вероят |
||
Пусть |
= К 2 • Тогда |
г г= |
||||||||
7 |
||||||||||
|
А г |
|
|
|
|
cos (ft«sin <А |
||||
ность |
попадания |
в |
круг |
радиуса М будет |
|
|
||||
|
|
|
|
гт |
|
|
і+к2 |
|
|
|
|
Л |
|
/ |
Г |
Г |
|
d |
f |
||
|
p" ~ t ' w |
Y * |
F |
cos |
f i - К s in |
|||||
|
f_ |
|
||||||||
Вероятность |
попадания в |
круг |
радиуса г =А |
ыояет |
||||||
быть |
рассчитана |
так: |
|
|
|
|
|
|||
|
|
|
|
|
, 2 |
|
|
)s . |
|
|
|
A 2= x 2+ y = 't2( A c o $ (P + ß 2s i n 2cf |
|
Отсюда
г1
7 =
|
COS7<P + K 2s i n 2(f |
2S |
|
Р*~ ‘ ~2Ж j |
d ? , |
е* ? 1 г (с М’г , / r W f >. |
167
И, наконец, если
R ^ ] f Ä ß 7
то
ъ*(А2cos2<f + t zs i n f ) ■
Здесь
2 ________К_______
Г = cos2f t K 2s i n 2f
2Я
к
2 [cos2f t К 2sin 2f \
о
Таким образом, вероятность попадания в круг в за висимости от способа вычисления радиуса является функ цией К
|
|
|
|
|
Таблица 3 .I |
|
||
|
^ A zt ß 2 |
|
А |
|
]fi)ß |
|
||
0 |
0 , 6 8 |
|
0 , 6 8 |
|
0 , 0 0 |
|
|
|
I |
0,63 |
|
0,39 |
|
0,39 |
|
|
|
Из таблицы |
вероятностей |
(табл. 3.1) |
в и н е , что |
|
||||
при изменении К от |
0 до I |
вероятность |
нахождения |
точки |
||||
в круге меньше |
всего |
изменяется |
при |
первом способе |
|
|||
(М = ^ А г+ Зг ) |
вычисления радиуса |
Можно показать,что |
||||||
и при промежуточных |
значениях К |
вероятность |
нахожде |
|||||
ния точки в пределах круга |
радиуса А1 |
также |
меняется |
|||||
незначительно и не выходит |
за пределы, |
указанные |
в |
|||||
таблице. |
|
|
|
|
|
|
|
|
168
Вероятность нахождения точки в пределах круга радиу
са 2 Мдостаточно |
велика и составляет |
|
|||
А |
, 1 |
0 , 2 |
0,5 |
0,7 |
1 , 0 |
0 |
|||||
& |
|
|
|
|
|
м |
6 8 |
6 8 |
6 6 |
64 |
63 |
2 М |
98 |
97 |
97 |
96 |
95 |
В официальных руководствах по навигационному обе спечению в качестве предельной ошибки принята величина 2,58М. Вероятность нахождения точки в пределах такого круга равна 99$. Численно эта вероятность совпадает с
вероятностью линейной ошибки, равной 2,58(3" . При из-
д
ѵенении —z- вероятность варьируется в пределах 2 £.
5
При известных элементах эллипса величина Мрассчиты вается преете;
|
|
|
M ^ A zt ß 2' |
■ |
|
|
Па основании теоремы Апполония можно получить такое |
||||||
же простое |
выражение для |
М, используя величммм сопряжен |
||||
ных полудиаметров |
эллипса |
(векториальные |
оии<ши) |
|||
-& ‘ Ѵ г + ѵ ; |
9 |
|||||
|
|
А ‘ |
|
|
|
|
M e у , . |
Ѵ г - |
сопряженные |
полудиаметры; |
|||
А , |
В |
- |
главные |
полудиаметры эллипса. |
||
|
|
|
|
|
Отсюда
М = 1ІАг+Вг •
169
§ 23. Эквивалентные преобразования систем уравнений погрешностей
При наличии избыточных измерений линии положения, как правило, не пересекаются в одной точке и образуют фигуру погрешности. Для приведения результатов наблюде ний в согласие необходимо выполнить уравнивание и найти вероятнейшие значения координат. Геометрически это со ответствует такой точке, сумма квадратов расстояний от которой до каждой линии положения становится минималь ной.
Пусть исходные линии положения представлены в общем виде:
А f t - ё;А иУ + = ca. ; |
|
I * 1 , 2 ,3 , • • •, п ? |
|
п - число измерений. |
|
Каждое уравнение ошибок имеет свой вес |
. Соответ |
ствующая исходная система нормальных уравнений будет иметь вид:
\paa~^A<f>+ [ра # ] д US + [ра<Г]~Оі
[ p a é ] A f+ \p é é ] A us + [ p t t b o .
Заменим исходную систему уравнений погрешностей эквивалентной ей системой, сосгоящей из двух уравнений. Согласно поставленному условию решение по методу наи меньших квадратов исходной и эквивалентной систем должно совпадать.
Пусть уравнения эквивалентных линий положения будут представлены в нормальной форме:
170
c o sT Д<р + sinTA w - i-A N /^V i |
у |
S in T Д ^ - c o sT A u f + A Ng=V2 |
• |
Веса этих линий положения £> и Р2 |
90+ Т - азимут, |
определявдий направление первой линии; Т - азимут вто
рой линии.
Уравнения поправок эквивалентных линий полохения приводят к следу«вей системе нормальных уравнений«
(j>l cos2T+P2s in T ) A P+(Pl ~P2 ) s in TcosTA as ~
-(P'Co sTa N * P 2 sin T A N 2) =0 у
(jfj-P^)Sin Teas TA(f >i'(Pt S in 2T +P2cos2T ) a u S -
-(Pt s in T A N -P 2cosTAN 2 ) = 0 .
Нормальные уравнения исходной и эквивалентной систе мы могут различаться только по форме. Сравнивая коэф
фициенты при неизвестных, получим, |
что |
(P'Cos2T+P2 s i n T ) = [ p а а J |
j |
(P 'SirfT* P2 cos2T )= [p ß É J |
* |
(Pt ~P2 ) s in Tcos T = \p a ë ] .
Слоят и вычтем два первых равенства и перепишем третье:
Р,+Рг ~ \p a a ] + [ p é é J ;
(Pr P2)c o s 2 T = [ p a a ] -[ p é é ]
I7I
|
( e ;- p 2 ) s i « 2 T - 2 [ f * t J • |
|
||||
Сравнение свободен* членов дает |
|
|||||
-(P'Cos |
TAN'* Рг sin Ta N2)~ [pa £ J |
, |
||||
- |
(P' siгг Та Nr Pz cos Та n2) —\ f i £ ] |
- |
||||
Из этих систем уравнений мои !м найти |
величиям |
|||||
|
|
р, , |
Р, |
« |
т |
|
затем |
найдем |
A N t |
и |
A N 2 |
• |
|
Таким образом, |
задача определении кеиэаеетимк и |
|||||
оценки их точности полностью решается: |
|
|||||
|
|
[fcrt]cts Т+ [ f it ] Sin Т |
|
|||
|
a n = - |
|
Pr |
|
|
|
|
|
|
|
|
|
|
|
a N , |
[patjsin T- [pg£J cos T |
|
|||
|
|
. |
|
|
|
|
|
$ |
|
\?а а ] - { р и ] |
|
Это уравнение показывает, что экстремальные овибки расположены по направлениям двух эквивалентных линий полохения, модули экстремальных овибок равны Pt и Р0 •
Как обычно, в качестве весов приняты величины, обрат ные дисперсиям, то есть
Р = ~В2
где А и В - полуоси результирующего эллияса овибок.
Нахождение элементов результирующего эллипса и вероят
1 7 2
нейших значений координат (или переноса A N f , ДЫг )
можно упростить, применив графоаналитический прием ре шения.
Если уравнения исходных линий положения даны в нор мальной Форме, то
а і = cos Ti ; |
= s in |
= - д п ; ■ |
Теперь получим:
р,+р2 = 1 °] ;
(P-Pg )cos2T=[f>cosZZ]
(Pr Pz )sin2T= [P sin2 7] ;
P ,A N t ^ [ p A n c o s ( T - T ) \= P cr )
P2AN = [ fA « s in ( 7 - T ) ] = P v ;
второе и третье уравнения легко представляются геомет рически. Для их решения нужно построить полигон из
в е к т о р о в , направив эти векторы под углом Z Z к из
бранной системе координат. Геометрическая сумма этих векторов по абсолютной величине
Угол между |
замыкающей полигона и осью Х равен 2Т. |
Таким образом, |
подучаем следующее репение: |
ITS
>{[Л -?] :
Для определения |
переносов |
a N, |
и a Nz построим |
||||
полигон переносов, |
составленный из |
векторов |
|
|
• |
||
Дирекционные направления этих |
векторов |
Т2- |
. |
Проектируя |
|||
замыкающую полигона |
переносов |
на оси & |
и |
V |
, |
имеющие |
|
направления Т и Т+900, получим проекции |
Рц |
и р |
: |
Дальне по переносам ANt и a N2 надо построить
эквивалентные линии полохения и в их пересечении по лучить место. Рассмотренный выне способ эквивалентного преобразования исходной системы уравнений погреиностей достаточно удобен для ручных вычислений. Однако схема решения в этом случае все-таки останется громоздкой. В ряде задач оолее целесообразно использовать другие эк вивалентные преобразования. Поэтому рассмотрим общие вопросы эквивалентных преобразований. Согласно опре делению эквивалентными системами называются такие, у которых ревения, полученные по методу найменьиих квад ратов, одинаковы. А для этого достаточно, чтобы исход ные системы имели одинаковые нормальные уравнения.
Известио, что в самом общем виде систему уравнений погрешностей можно записать так»
AX- L=V - |
(ЗЛО) |
174
Запилен эквивалентную ей систему в такой хе форме:
|
a ' X - L ' - V ' 1 |
( з . Н ) |
|
где a '-Q A |
I L‘=QL я |
К = Q V . |
|
Здесь Q - матрица преобразования. Обозначим через
р ' весовую матрицу эквивалентной системы и напиием систему нормальных уравнений:
A rQ Tp ,C ( A X - A Tq TP lQL = 0 .
Tax как системы уравнений погрешностей по условию эквивалентны, то системы нормальных уравнений должны быть одинаковыми. Отсюда получим условие эквивалентнос ти двух систем:
A TQ TP ' Q = A TP . |
(3.12) |
Подбирая матрицу преобразования Q |
так, чтобы ус |
ловие (3.12) выполнялось, найдем способ эквивалентного преобразования исходной системы.
Определим зависимость между корреляционными матрица ми исходной и эквивалентной систем уравнений погрешнос
тей. Для этого обе |
части |
(3.12) |
умножим справа |
на K Q r9 |
|||
где к - корреляционная |
матрица исходной системы. |
Р |
|||||
результате |
получим |
в правой части (3.12) A |
Q T |
, |
а в |
||
левой части |
A TQTP 'Q К 6)Т |
|
|
|
|
||
Для того |
чтобы |
равенство не |
наруиалось, |
необходимо |
K ,= (P ')~ '= Q K Q T.
Отсюда следует, что при эквивалентном преобразовании корреляционная матрица системы преобразуется так же,как
175
преобразуется корреляционная матрица случайного векто ра при его линейном преобразовании.
Большой интерес при использовании вычислительной техники представляет собой эквивалентная замена систе
мы уравнений |
погрешностей системой, состоящей из п |
|||||
уравнений |
с |
п |
неизвестными. Рассмотрим этот случай |
|||
подробнее. |
|
|
|
|
|
|
Полоним |
Q= C~‘A ‘P |
. Здесь |
С |
- матрица коэффициен |
||
тов нормальных уравнений исходной системы. |
||||||
Тогда |
|
|
|
|
|
|
|
А '- Q A -C 'W P A = C~‘C - J |
|||||
|
l!= ql = c 'a tp l = x . |
|
|
|||
X - матрица оценок |
неизвестных, |
найденных по методу |
||||
наименьших |
квадратов. |
|
|
|
||
Таким образом, эквивалентная |
система приобретает в |
|||||
данном случае |
следующий вид: |
|
|
х- х = ѵ ‘ ,
аее корреляционная матрица
к'= Ц Р ~'О іТ= С ~1 .
Пример. Пусть дана система уравнений линий положе ния:
а, |
; |
а2Х |
Уг ’ |
Х" тCr>jt ~^гГ |
’ |
каждое из которых имеет вес |
Р і |
Л
176