Студопедия

КАТЕГОРИИ:


Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748)

Априорная оценка точности рассчитанных координат





Доверь свою работу кандидату наук!
1500+ квалифицированных специалистов готовы вам помочь

Априорная оценка точности измеряемых навигационных параметров основана на многочисленных статистических исследованиях, которые происходили до конкретной обсервации. В качестве основного параметра для априорной оценки точности измеряемых величин применяется средняя квадратическая погрешность измеряемого навигационного параметра т.

Особенностью определения координат является тот Аакт, что измерения -косвенные, т.е. измеряются навигационные параметры и допущенные погрешности затем переносятся в погрешности координат. Рассмотрим процедуру переноса погрешностей измерений в погрешности координат на примере ОМС по двум измерениям.

В этом случае линеаризованная система принимает вид: ;

Так как измерения имеют погрешности, то перепишем систему в виде

Тогда

откуда

Формирование ковариационной матрицы погрешности измерений выполняется по формуле

где D - обозначение ковариационной.

Для двумерного случая выражение (2.24) выглядит так:

а операция математического ожидания, выполненная с матрицей 8U6U , превраща­ет ее в ковариационную матрицу D,


На главной диагонали D находятся дисперсии .измеряемых навигационных параметров, а вне диагонали - ковариационные моменты, которые характеризуют статистическую связь между измерениями.

Аналогично определим ковариационную матрицу погрешностей искомых параметров, используя правила матричного исчисления (ABYl=B~\A~\ и (Д-1)-^^)-1.

Получим

N = D(A^) = D(SXSX7) = D^A-'SU^A^SUY] = (A'D^A)^. (2.25)

В дальнейшем при написании ковариационных матриц, где это не вносит двузначности, будем опускать аргумент при D.

Для двумерного случая ковариационная матрица N имеет вид:

 

где п 11 - дисперсия погрешностей широты, n-г.г - дисперсия погрешностей отшест-вия, ni2 = пг\ - ковариационные моменты.

Вся информация о погрешностях содержится в матрице N. В судовождении часто используется ее геометрическая интерпретация в виде эллипса погрешностей. Установим- связь между элементами матрицы N и параметрами эллипса:



полуосями и углом ориентации.

В общем случае такая задача рассматривалась Г. Хоттелингом в 1933 г. Ученым было доказано, что для ковариационной матрицы существуют . векторы, направлениям которых соответствуют максимальные и минимальные значения рассеивания (погрешностей). Эти значения соответствуют собственным числам матрицы. Направления собственных векторов, указывающие на направление максимального и минимального рассеивания (дисперсии), соответствуют направлениям полуосей эллипса. Собственные числа - это экстремальные значения дисперсий. Для перехода к линейным величинам - полуосям эллипса (гиперэллипса для /7-мерного пространства), необходимо извлечь квадратный корень.

Рассмотрим эту задачу для двумерного случая, т.е. для плоскости. Физика и геометрия собственных чисел и векторов заключается в том, что результатом

"<

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

Поставим численный эксперимент, который прояснит эту запись. Выполним умножение Nz, где в качестве z будем выбирать единичный вектор с направлением У от 0 до 360°. Формирование компонент единичного вектора



выполним по формуле:


В качестве примера возьмем матрицу



Рис. 2.4. Результат умножения N -z


В результате умножения матрицы N на z конец вектора р опишет эллипс (рис. 2.4). Процедуру умножения можно рас­сматривать как оператор, преобразую­щий единичный вектор z. После пере­множения вектор изменит направление и длину. Результаты такого перемно­жения с дискретностью в один градус

 


- компоненты вектора р (значения Х и У) приведены в табл. 2.1.

Таблица 2.1

У ^1 Y X R
8,13 21,000 3,000 21,213
8,48 21,049 3,139 21,282
8,83 21,092 3,277 21,345
9,18 21,128 3,415 21,402
9,53 21,158 3,551 21,454
9,87 21,182 3,686 21,500
10,22 21,199 3,820 21,540
10,55 21,209 3,953 21,574
10,90 21,213 4,084 21,603
11,24 21,211 4,215 21,625
11,58 21,202 4,344 21,642
11,92 21,187 4,471 21,653
12,26 21,165 4,598 21,658
12,60 21,137 4,723 21,658
12,93 21,102 4,846 21,651
13,27 21,061 4,968 21,639
...     ... ... ...
95,38 -0,692 7,358 7,390
98,30 -1,062 7,281 7,358
101,24 -1,432 7,201 7,342
104,19 -1,801 7,120 7,344
107,13 -2,169 7,037 7,363
110,05 -2,537 6,951 7,400
112,94 -2,905 6,863 7,453
115,72 -3,271 6,113 7,522

 


В графе «У» указано направление единичного вектора z, в графе «Ti» -направление уже преобразованного вектора р. В графе «R» приведены значения длины вектора р. Из табличных данных видно, что расхождение в направлении вектора z и вектора р - величина переменная, но в районе 12° и 102° эти направления совпадают. Кроме того, им соответствуют максимальное и минимальное значение длины R. Таким образом, направления собственных векторов 12° и 102° - ортогональны. Собственные значения равны приблизительно 21,658 и 7,342 соответственно.

Для двумерного случая можно получить простые формулы. Согласно выражению (2.26), запишем:

или

а также представим в матричном виде:

(N-AE)z=0


где


- нулевая матрица.


Формально получаем



Следовательно,


Так как Z - произвольный вектор и, в общем случае ненулевой, то

Запишем для двумерного случая



- квадратное уравнение. Решая его относительно X и,




принимая во внимание, что п-ц = пц (т.к. матрица N симметрическая), получим

Подставив значения из матрицы (2.27), получим а] == 21,659; Л^ = 7,341. Эти значения практически совпали с максимальным и минимальным значениями, приведенными в табл. 2.1.

Определим ориентацию собственных векторов, соответствующих найденным собственным значениям. Считая ^ известным, подставим это значение в уравнения системы (2.28) и решим ее относительно zi и z:, учитывая, что z/ == cosfF), z-г = = sinCF).

Первое уравнение системы (2.28) будет выглядеть так:


Разделив в первом уравнении левую и правую часть на cos(T), получим:


откуда

Подставив числовые значения, получим ^F = 12,388.

Таким образом, фактически получено направление большой полуоси эллипса У относительно норда. Если в уравнение (2.31) подставить другое значение Л, то получим направление малой полуоси, но так как они ортогональны, то практически это не требуется.

Для отыскания полуосей необходимо извлечь квадратные корни из собственных чисел

Когда говорят об оценке точности, то обычно добавляют слова априорная или апостериорная. Априорная - это оценка точности, выполненная по информации о погрешностях измерений, полученной ранее. Как правило, такая информация о точности измеряемых навигационных параметров основывается на многочисленных статистических исследованиях, которые происходили до конкретной обсервации в каких-либо осредненных условиях. Именно такая информация обычно содержится в ковариационной матрице погрешностей измерений, используемой при расчете координат. В формуле (2.18) она обозначена как D. Если погрешности измерений статистически независимы, то внедиагональные элементы равны нулю и матрица имеет вид:

Эти погрешности, в соответствии с правилом их переноса, формируют априорную ковариационную матрицу определяемых параметров.

Процедура построения эллипса погрешностей по ковариационной матрице сводится к следую­щим операциям:

- расчету собственных значений ^ по формуле (2.29);

- определению угла ориентации У по формуле (2.31);

- расчету полуоси по формулам (2.32).

Рис. 2.5. Связь матрицы N с эллипсом погрешностей

На рис. 2.5 показана связь между элементами ковариационной матрицы и эллипсом. Отрезок, заключенный между касательной к эллипсу параллельной оси Y и самой осью, соответствует

СКП по широте:

Отрезок на оси Y,


отсекаемый вертикальной касательной, соответствует СКП по отшествию:

На рис. 2.5 также показана СКП обсервации М, которая рассчитывается как корень квадратный из следа ковариационной матрицы либо с помощью полуосей эллипса:

Поможем в написании учебной работы
Поможем с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой




Дата добавления: 2014-11-29; Просмотров: 1762; Нарушение авторских прав?;


Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет



ПОИСК ПО САЙТУ:


Читайте также:
studopedia.su - Студопедия (2013 - 2022) год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! Последнее добавление




Генерация страницы за: 0.02 сек.