Система уравнений не нарушится, если левую и правую часть каждого из уравнений системы (59) или (62) умножить на одно и то же число.
В матричном виде эта операция соответствует умножению системы слева на матрицу вида:
, (68)
элементами которой являются величины, обратные СКП измерений (1/m). После выполнения такого умножения система (62) будет иметь вид:
(69)
Последовательность вывода ничем не отличается от вывода изложенного в предыдущем параграфе.
Выражение для Q, (64) ,перепишется так:
(70)
Обозначим , (71)
- весовая матрица, обратная ковариационной матрице погрешностей
(72)
Система нормальных уравнений запишется так:
(73)
а решение будет иметь вид:
(74)
Если элементы матрицы (69) равны, это означает, что измерения равноточные и (74) превращается в (67). Формула (73) носит более общий характер, чем (66).
Вероятнейшая точка, полученная по формуле (74) - является средневзвешенным значением измеренных параметров.
После того как система нормальных уравнений будет решена, определяют обсервованные координаты. При необходимости применяется итерационная процедура Ньютона. Итерационная процедура Ньютона применяется для уменьшения погрешностей линеаризации изолиний. Обычно, определив место, необходимо провести анализ точности измерений.
- измерить навигационные параметры Uo1, и Uo2„
- найти счислимые значения Uc1 и Uc2 этих же навигационных параметров по счислимым координатам (jс и lс) на момент измерений;
- вычислить приращения навигационных параметров DU1 и DU2, рассчитать модули градиентов g1, g2 ,
- вычислить коэффициенты системы уравнений;
- сформировать матрицу частных производных;
- сформировать ковариационную матрицу погрешностей измерений;
- найти матрицу коэффициентов нормальных уравнений;
- найти ковариационную матрицу погрешностей координат;
- найти псевдообратную матрицу
- найти вектор приращения координат;
- вычислить географические координаты обсервованного места:
jo = jc +Dj
lj = lc + Dw secjт
Если приращения географических координат отвечают требованиям точности, то решение задачи заканчивается, если нет, то полученные координаты записывают как счислимые, и цикл счета повторяется. Процедуру выполняют до тех пор, пока требования не будут удовлетворены. Обычно во всех алгоритмах определения места судна предусмотрена данная циклическая схема расчета. Как правило, итерационный процесс включает не более трех циклов.
Изучив основные технические приемы использования МНК в обработке навигационных измерений, полезно сделать анализ основных свойств метода.
Во-первых, необходимо различать статистические и нестатистические свойства. Ясно, что в какой-то эпизодической обсервации по избыточному числу линий положения МНК выступает, как формальный и наиболее простой математический аппарат, используя который можно однозначно фиксировать точку относительно фигуры погрешностей при условии, что систематических погрешностей в измерениях нет или они включены в определяемые переменные, т. е. в векторDХ. В этом случае фигура погрешностей с точностью до случайных погрешностей измеряемых параметров находится в окрестности истинной точки. Поэтому разовая обсервация, привязанная к фигуре погрешностей, с точностью до случайных погрешностей измеряемых навигационных параметров будет находиться в окрестности истинной точки (рис. 19).
Рис. 19. Возможное расположение треугольника, погрешностей измерений относительно истинной точки Z: D 1, 2, 3 - истинные случайные погрешности измерений линий положения |
В пределах доверительного интервала, являющегося суммой доверительных интервалов вершин многоугольника погрешностей, размеры фигуры погрешностей определяются случайными погрешностями измерений, а поэтому ее размеры не могут служить характеристикой точности измерений. В этой области малая фигура не обязательно должна быть ближе к истинной точке. Систематические же погрешности измерений приводят к четкой аналитической зависимости размеров фигуры и расстояния до истинной точки.
Теперь есть смысл говорить о вероятностно-статистических характеристиках МНК, если результаты измерений удовлетворяют следующим условиям:
значения коэффициентов линий положения известны точно;
результаты измерений содержат лишь случайные составляющие погрешностей;
погрешности измерений имеют нормальный закон распределения.
Свойства МНК таковы.
1. Несмещенность оценки, т. е. математическое ожидание оценки, равна истинному значению. Это свойство имеет явно асимптотический характер и может наблюдаться в большой совокупности измерений.
2. Состоятельность оценки, т. е. оценка при увеличении п,, стремится к истинному значению параметра по вероятности:
Р{|Х-Х|>x}® 0, е® 0.
3. Эффективность оценки (оценка называется эффективной, если она в данном классе оценок имеет минимальную дисперсию).
4. Линейная процедура при нормальном законе распределения погрешностей измерений определяет и нормальность оценки. 180
5. При нормальном законе распределения МНК является частным случаем метода максимального правдоподобия.
6. Линейность обработки делает МНК очень чувствительным к наличию промахов, систематических погрешностей измерений и к неравноточности навигационных параметров.
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.