Коррелатный способ уравнивания по МНК.



Мы поможем в написании ваших работ!


Мы поможем в написании ваших работ!



Мы поможем в написании ваших работ!


ЗНАЕТЕ ЛИ ВЫ?

Коррелатный способ уравнивания по МНК.



Пусть имеем n – измерений, из которых t – необходимые:

y1, y2, y3…yn – количество выполненных измерений

Y1, Y2, Y3…Yn – истинные значения найденных величин результатов измерений.

Связь между ними может быть выражена:

Yi + Δi = yi (1)

Устанавливаем систему весов измеренных величин:

P1, P2, P3…Pn

Вычисляем обратный вес:

П1, П2, П3…Пn Пi = 1/Pi (2)

Подсчитываем число избыточных измерений:

R = n-t (3)

Каждое избыточное измерение приводят к математическому соотношению между истинными значениями измеряемых величин, т.е. в геодезической сети возникает r – условий.

Φj (y1, y2…yn) = 0 (4) – исходная система условных уравнений связи.

J = 1, 2, 3…r.

Вследствие неизбежных ошибок в измерениях в те же функции, но от измеренных величин примут вид:

Φj (y11, y22, …yn + Δn) = Φj (y1, y2…yn) = Wj (5)

Это система условных уравнений связи от измеренных значений

Wj – невязки, [Δi] = Wj (6)

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

[Vi]j = -Wj (7)

Vi – поправки в результаты измерений, I = 1, 2…n.

J = 1, 2…r. – по числу условий.

Yiур = yi + Vi (8)

Тогда система φj:

Φj(y1ур, y2ур…ynур) = Φj (y1+V1, y2+V2…yn+Vn) = 0 (9)

Это система условных уравнений связи от уравненных значений.

В правой части получились 0.

Систему (9) приведем к линейному виду, раскладывая в ряд Тейлора, пренебригая при этом малыми нелинейными членами разложения.

Φj (y1+V1, y2+V2…yn+Vn) = Φj (y1, y2…yn) + (dφj/dy1)V1 + (dφj/dy2)V2 +…+ (dφj/dyn)Vn = 0 (10)

На основании (5) получим:

(dφj/dy1)V1 + (dφj/dy2)V2 +…+ (dφj/dyn)Vn + Wj = 0 (11)

Обозначим частные производные от первой функции по каждому результату измерений через ai.

1/dyi = ai, dφ2/dyi = bi, dφ3/dyi = ci

1/dy1 = a1 1/dy2 = a2 1/dyn = an
2/dy1 = b1 2/dy2 = b2 2/dyn = bn
r/dy1 = r1 r/dy2 = r2 r/dyn = rn

Система (11) с учетом (12) примет вид:

A1V1 + a2V2 +…+ anVn + W1 = 0

b1V1 + b2V2 +…+ bnVn + W2 = 0 (13)

r1V1 + r2V2 +…+ rnVn + Wr = 0

(13) – это система условных уравнений поправок.

Wj – невязки, Ai, bi, ci – коэф – ты при поправках,Vi – неизвестные поправки, которые надо найти, решив (13).

Но в (13) число поправок равно n, а число уравнений r, т.е. число уравнений r меньше числа неизвестных поправок (r<n). Такая система имеет множество решений, т.е. не решается однозначно. Чтобы выбрать из этого множества решений один наилучший вариант, который дает вероятнейшее решение нужно поставить дополнительные условия, которые являются принципом МНК (14): [pV2] = min (14)

Тогда решая систему (13) под условием (14) переходят к системам нормальных уравнений коррелат.

[Паа]К1 + [Пав]К2 + [Пас]К3 +…+ [Паr]Кr + W1 = 0

[Пав]К1 + [Пвв]К2 + [Пвс]К3 +…+ [Пвr]Кr + W2 = 0 (15)

[Паr]К1 + [Пвr]К2 + [Пcr]К3 +…+ [Пrr]Кr + Wr = 0

В этой системе К1, К2…Кr – неизвестные коррелаты или дополнительные множители Лагранжа, которые вводятся искусственно в процессе перехода от системы (13) к (15). Их число r как число уравнений такая система решается однозначно. А способы решения могут быть различны:

1. По схеме Гаусса.

2. По спец. Разработанным программам.

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

3. ачно. жа, которые вводятся искусственно в процессе перехода от системы (13) к (15).

Из решения (15) находятся неизвестные коррелаты, а по ним поправки Vi.

Vi = Пi (aiK1 + biK2 +…+ riKr) (16) – коррелаты уравнений поправок.

Контроль:

[Vi]j = -Wj (17)

Затем yiур = yi + Vi (18)

Контроль:

Φj (y1ур, y2ур…ynур) = 0 (19)

Оценка точности по результатам уравнений выполняется по поправкам по формуле:

μ = √[pVV]/r (20)

Чтобы оценить конкретный элемент сети, например, отметку, коор-ту, угол необходимо этот элемент выразить математически и тогда ошибка функции будет:

mF = μ√1/pF (21)

Уклонение отвесных линий.

УОЛ – угол, образованный нормалью к поверхности эллипсоида вращения и отвесной линией.

УО бывают:

Абсолютные – угол, образованный отвесной линией и нормалью к общеземному эллипсоиду. Зависят от неравномерного распределения масс в теле земли.

Относительные – угол, образованный отвесной линией и нормалью к референц-эллипсоиду. Зависят от размеров и ориентировки поверхности относимости.

1 для изучения фигуры земли;

2 для редуцирования измеренных величин с физической поверхности земли на поверхность эллипсоида вращения;

3 чтобы иметь возможность связывать астрономические и геодезические координаты: j , l « B, L;

4 чтобы связать астрономические и геодезические азимуты: a* « А.

Астрономическая широта – острый угол, образованный отвесной линией и плоскостью экватора.

Астрономическая долгота – двугранный угол, образованный плоскостью Гринвичского меридиана и плоскостью астрономического меридиана данной точки.

Плоскость астрономического меридиана – плоскость проходящая через отвесную линию данной точки ^ плоскости земного экватора.

Астрономический азимут – двугранный угол, образованный плоскостью астрономического меридиана точки и вертикальной плоскостью, содержащей отвесную линию и визирную цель.

U – Уклонение отвеса

ZГ – геодезический зенит

ZА – астрономический зенит

Р – полюс мира

1) U, q 2) x, h

È ZГ G – проекция полного уклонения на геодезический меридиан или составляющая в плоскости геодезического меридиана x

È ZА G – проекция полного уклонения отвеса на 1-й вертикал или составляющая в плоскости 1-го вертикала h

x = U cos q tg q = x / h h = U sin q - прямая связь U = Ö x2 + h2 – обратная связь

Способы определения составляющих уконений отвесных линий.

По результатам наземных измерений можно определить составляющие отвесных линий: 1 Астрономо-геодезический

xАГ = j - В hАГ = (l - L) cos j

Достоинство: - позволяет определить абсолютные и относительные уклонения отвеса

Недостаток: - очень дорогой способ.

2 Гравиметрический g, g

D g = g - g - аномалия силы тяжести

xГР = -1/2П ò0ò0 D g Q cos A a A d j;

hГР = -1/2П ò02Пò2П0 D g Q cos A a A d j

Достоинства:- очень дешевый способ;

- можно получить абсолютное уклонение отвеса

Недостаток:- нельзя на прямую получить относительное уклонение отвеса.

3 Астрономо-гравиметрический

Этим способом можно получить астрономо-геодезические (относительные) уклонения отвесной линии. Чтобы их вычислить нужно, чтобы была создана геодезическая сеть редких астропунктов и гравиметрическая съёмка территорий. А, В, С – пункты геодезической сети, на которой выполняются астрономические определения, расстояние между ними » 200 км

Д – рядовой пункт геодезической сети

Для точки А, В, С известны B, L, j , l

xАГ = j - В; hАГ = (l - L) cos j

Вокруг точки Д выполняем гравиметрическую съёмку

Rs » 300 km Dx = xАГ - xГР Dh = hАГ - hГР

Внутри s выделим ещё 1 область меньшего радиуса s0

Rs » 1/2 Rs

s0 – область в которой Dx и Dh изменяются по линейному закону при переходе от одной точке к другой

Dx = a1x + a2y + a0 Dh = b1x + b2y + b0 a1, a2, a0 b1, b2, b0

Когда коэффициенты найдены можно найти уклонение отвеса для любого пункта.

Для Д известно: DxД, DhД, xАГ, xГР, hГР

xДАГ = xДГР + xДГР hДАГ = hДГР + hДГР

Нормальная и геодезическая высоты. Области их применения.

Нормальная высота – высота, которая вычисляется по формуле

HKg = W0 – WK / gKm (1) HKg = ò0K gdh / gKm (2)

gK0 = gэкв0(1+ b sin2 Bk) gKm =gK0 - g / R*Hk dg / dH = R

В геометрическом смысле нормальная высота это отрезок нормали к эллипсоиду, заключенный между физической поверхностью земли и поверхностью квазигеоида. HKg = HKg + zк

Основные условия выбора системы счета высот:

1 высота не зависит от выбора трассы нивелирного хода;

2 высота вычисляется только по результатам измерений на физической поверхности земли без привлечения гипотез;

3 поправки в измеренные превышения за переход к выбранной системе счета высот малы;

4 высоты точек, лежащих на одной уровенной поверхности не одинаковы;

5 принятой системе счета высот соответствует достаточно простой способ вычисления геодезической высоты.

Геодезической высотой точки называется расстояние НГ = mM между поверхностью эллипсоида и данной точкой, отсчитываемое по нормали к нему в этой точке. Точки, находящиеся выше поверхности эллипсоида, имеют положительные высоты, а ниже – отрицательные. Точки, лежащие на поверхности эллипсоида, имеют геодезические высоты равные нулю (НГ = 0). Геодезические координаты непосредственно измерены быть не могут, их получают путем вычисления и суммирования приращений координат относительно исходного пункта государственной геодезической сети, имеющего координаты B0, L0, H0.

Билет 16

Параметрический способ уравнивания. Это способ абсолютного экстремума, при котором все измеренные величины представляют в виде функций некоторых независимых параметров.



Последнее изменение этой страницы: 2016-04-07; Нарушение авторского права страницы; Мы поможем в написании вашей работы!

infopedia.su Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Обратная связь - 3.236.231.61 (0.016 с.)