Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Коршунов / М07ОцнкФнкц.doc
Скачиваний:
37
Добавлен:
26.04.2015
Размер:
653.82 Кб
Скачать

Вектор-функция случайного вектора

Пример: при вычислении координат точек в системе снимка мы должны учесть влияние ошибок наблюдений координатных меток, а не только ошибок наблюдений точек. Отсчеты на координатные метки используются при вычислении координат и параллаксов всех точек снимка, значит, все вычисленные координаты х , у снимка (а для стереопары еще р и q) суть функции отсчетов Vx и Vу на точку и отсчетов +1, +2 +3 +4 на координатные метки, поэтому они будут коррелированны:

Столбец слева от знака равенства есть вектор функций случайных величин. Для вектора функций согласно (9.2) мы можем записать вектор A, состоящий из градиентов Fi . Образуется матрица. Это матрица первых производных, которую называют матрицей Якоби ее определитель - якобиан). Зная матрицу Якоби и ковариационную матрицу Kx случайного вектора , размер которой будет 2(n+4) на 2(n+4), мы вычисляем по (9.2) ковариационную матрицу вектор-функции ( 9.4) .Для наглядности раскроем матрицы ( 9.4)

Запись для вектора функций в алгебраическом виде получается крайне сложной:

(f/x1)2 S12 +(f/x2)2S22+. . . +(f/xn)2Sn2 + . . . +2f/x1 f/x2 r12S1S2+. . .+2f/xn-1f/xnrn-1 nSn-1Sn (это частица формулы)

Формулы (9.2) или (9.3) служат для вычисления дисперсии функций случайного вектора, элементы которого коррелированны. Они действительны при рассмотренных выше условиях: (1) значение функции вычисляется не в особой точке функции, и (2) на рассматриваемом интервале функция близка к линейной.

Если корреляция между элементами вектора Х равна нулю, т.е. все , то ковариационная матрица KF вырождается в диагональную, содержащую только дисперсии или же, что то же самое, корреляционная - в единичную. Тогда выражение для вычисления дисперсии функции ( 9.2) упрощается . (9.5) , а ковариационная матрица вектора функций ( 9.4) –станет диагональной (9.6)

Вычисление значения оценки

Пользуясь конкретными численными значениями аргументов и параметров, входящих в формулы, вычисляют численные значения частных производных . Если стандартSi аргумента xi соизмерим с используемым значением аргумента xi, то следует подсчитать значение соответствующих производных при , если данная величина находится в числителе, и при , если в знаменателе, и сравнить с ранее вычисленными. При отклонениях, более 10% от значения производной, в формулу оценки дисперсии функции подставляют большее значение производной.

Числовой пример. (1) Оценка точности определения превышения по разности продольных параллаксов.

Высота фотографирования H= 2100м определена с СКО SH= 200 м, продольный параллакс на снимке р = 70 мм измерен с СКО Sp= 1 мм, а разность продольных параллаксов р= I мм - с СКО Sp= 0.1 мм. Находим превышение h= 2I00 1/70 = 30м. Корреляция отсутствует. Вычисляем =

= 302 (I/I0.52+I/702+I/I02) = 8.I6 +0.I8 +9 = I7.34. Отсюда СКО определения превышения в 30 м будет Sh=√ I7.34= 4.16 м.

Так как относительная ошибкапри этих данных постоянна, то меньшее превышение определится точнее и наоборот.

(2) Так, если при тех же условиях разность продольных параллаксов составит 10мм, что дает превышение в 300 м, а точность ее измерения о Sp остается та же самая, то превышение в 300 м определится с СКО = (8I6+18+9)1/2 = 8431/2=29м. Мы здесь не учли того, что высота фотографирования и продольный параллакс изменились: Н уменьшилась, а р увеличился. Предполагается пересчитать самостоятельно.

Сравнивая составляющие ошибки в примерах (1) и (2), видим, что во втором влияние ошибки определения высоты фотографирования является преобладающим.

Рассмотрим, каково будет ее влияние при H- SH =2100-200 м. (нижний предел, ибо Н в производной стоит в знаменателе). Влияние составит 99% общей ошибки, которая возросла на 10%, т.е. общее СКО - 32 м. Эту оценку принимаем за оценку СКО функции. Отметим, что в примере (1) влияние ошибки измерения продольного параллакса, а в (2) - еще и ошибки разности продольных параллаксов пренебрегаемо малы.

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

Соседние файлы в папке Коршунов