Файл: Разумов, О. С. Пространственная геодезическая векторная сеть.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.10.2024
Просмотров: 57
Скачиваний: 0
Следуя правилу (1.85),. получим выражение для среднего квадратического тензора ошибок положения конечной точки век тора г относительно его начала
M l = CMl,ÖJCT |
(2.159) |
где Л4агбд — средний квадратический тензор |
ошибок непосредст- |
венных измерений. |
|
Если результаты измерений независимы друг от друга, то со
ставляющие элементы тензора (2.159) равны |
|
|
||
т2к — cos2 бcosЧтг2+ г2 sin2 бcos2 tm2 + |
г2 cos2 бsin2 tm2 |
|
||
т2у — cos2 бsin2 tin2-j- г2 sin2 бsin2 tin2 + |
г2 cos2 6cos2 tin2 |
|
||
in2 |
— sin2 6in2-f- r2 cos2 in2 |
|
|
|
1112 |
— cos26sin^cos//?r2-j- r2 sin2 6 sin ^ COS Й722 |
-f- |
. . (2.160) |
|
|
+ г2 cos2 бsin t cos tm2 |
|
|
|
in2 |
— cos 6sin 6 cos tin2-— r2 sin 6 cos 6 cos tin2 |
|
|
|
xz |
г |
о |
|
|
tn2z = sin 6cos sin tin2— r2 sin 6 cos 6 sin im2
Геодезический вектор L определяется как сумма (разность) астрономических векторов г, и при строгой синхронности наблю дений на станциях средний квадратический тензор ошибок поло жения конечной точки вектора L относительно начала будет равен сумме тензоров (2.159)
(2.161)
Если синхронность работы двух станций нарушается некоторой ошибкой А0, то чтобы учесть влияние этой ошибки на общую погрешность вектора L, необходимо знать возможный сдвиг спут ника по орбите за время ДѲ. В соответствии с (1.149) имеем
(2.162)
а совокупное влияние ошибки несинхронности наблюдений найдем как
М \%= Se(ЛѲ)2 Si « 2SoinlSl, |
(2.163) |
где тв — средняя квадратическая ошибка фиксации |
моментов |
наблюдений на станциях.
В итоге, для одной серии наблюдений будем ішеть
(2.164)
121
При этом средняя квадратическая ошибка конечной точки вектора окажется равной
т~г = т;-х + т~ии+ т;: \т~гХ\ + [/'2 /гг5І Т+ ]r2cos26m?]r +
-|- 2сгп2т%, |
(2.165) |
а при условии, что т&=m tcos 6= р.
nij- — \піг\ \ -f- 2 [rp2]f + 2ö2/i2mo. |
(2.166) |
В формулах (2.165 и 2.166) сдвиг спутника определялся по сред ней скорости его движения п; а — большая полуось орбиты ІІСЗ.
Имея тензор ошибок взаимного положения наблюдательных станций, можно получить, используя (1.100 и 1.103), средние квад ратические ошибки составляющих элементов вектора L. Заметим, что задача оценки точности в значительной степени упростится, если ошибка несинхрониостн пренебрегаемо мала и результаты непосредственных измерений отвечают условиям
тг _ |
' "'б _ |
т "іcos 6 |
|
|
и |
|
|
|
(2.167) |
|
ГПа |
< 0,3j.ir |
|
|
|
ап —у |
|
||
Тогда |
|
|
|
|
Л'/г = (м-і И + |
М-2r\)Ew |
(2.168) |
||
(£33 — единичная матрица) |
|
|
|
|
и |
|
|
|
|
ml = L2— |
= U cos2 1 — |
= (р? г2\ р? г2). |
(2.169) |
|
Р2 |
|
Р2 |
|
Соотношения (2.168) и (2.169) показывают, что при прочих равных условиях итоговая ошибка положения конечной точки векторов L относительно его начала возрастает с увеличением топоцентрических дальностей г. И если не принимать во внимание искажающих воздействии рефракции, то оказывается, что для уменьшения ошибки т - спутник выгоднее всего наблюдать на
малых высотах в моменты пересечения его траекторией геоцентри ческого сечения хорды (е= 0°).
При равноточных измерениях минимальная ошибка конечной точки вектора L составит
т = m2tx + т1у nVzz = 6рѴ2, если rx — r2 — r, (2.170)
122
а принимая во внимание зависимость точности непосредственных измерений от зенитного расстояния (2.40), можно записать
О |
6Ро |
(2.171) |
COS 2 |
|
|
|
|
|
где ß — угол между хордой |
и направлением на НСЗ. |
|
Минимум этой функции, как показал Н. Л. Макаренко в [38], |
||
наступит при условии, когда |
|
|
t g ? ^ t g f + |
^ j / 9 t g s f + 8 |
(2.172) |
(например, для хорды длиной L = 2800 км, -^-=12°,6, минималь
ная средняя квадратическая ошибка будет при ß= 41°,8 и 2 = 60°,8). Опираясь на уравнения (2.172), (2.103) и (2.108), можно рас считать оптимальную высоту ИСЗ над поверхностью Земли, при использовании его в качестве визирной цели в рассматриваемой программе наблюдений, и при заданной точности окончательного
результата подсчитать необходнмое_число измерений.
Для оценки точности вектора L по результатам нескольких серий синхронных измерений необходимо и достаточно подсчитать весовые тензоры в каждой из серий
Р і = ( М г ) ~ \ |
(2.173) |
п найдя их сумму сделать обратный переход
Р ц = [Р,-]; М ц = р ц . |
(2.174) |
Заканчивая разбор программы 1, следует заметить, что нару шение синхронности наблюдений может привести и к системати ческим ошибкам в элементах вектора L. Тогда, чтобы уменьшить значения этих ошибок, спутники целесообразно наблюдать при различных прохождениях через хорду. Из геометрической сущ ности метода следует, что если траектория ИСЗ будет пересекать хорду примерно посредине и под прямым углом к ней, то при нарушении синхронности наблюдений более надежно будет опре деляться длина хорды, а при пересечении под острым углом — ее направление.
Кроме рассмотренной программы наблюдений 1, для опреде ления длины и направления хорды, соединяющей наблюдательные станции, могут применяться программы 2 п 3; в этих программах распределение ошибок подчинено более сложным^ закономерно стям, так как положение конечной точки вектора L определяется линейной или угловой засечкой относительно позиции ИСЗ, за данных топоцентрическими векторами (см. рис. 27 и 28). Надеж ность определения элементов вектора L такими способами зависит
123
от геометрической формы засечек и точности непосредственных измерений. Для расчета погрешностей положения конечной точки здесь могут применяться формулы оценки точности засечек с уче том ошибок положения промежуточностью позиций ИСЗ.
Вместе с тем некоторые исследователи [18], [91] справедливо полагая, что проведение фотографических наблюдений ИСЗ проще организации лазерных измерений, предлагают определять длину хорды уже после нахождения ее направления в пространстве. Тогда искомую длину отыскивают из решения треугольника (см. рис. 13), в котором измерена одна или две стороны, а приближен ную оценку точности этого результата осуществляют по формулам: для программы 1
т1 = |
m-r cos2 ßx |
-г tri1 cos2 ß2 -1- рj r\ sin2 ßx -)- p,2 r1sin2 ß2; |
||
для программы 2 |
|
|
||
m-L |
f |
pJL2 ctg2 (ß, + ß2) + p |L 2 |
sin2 ßi |
|
sin2 ß, sin2 (ßL-|- ß2) ’ |
||||
|
|
|
||
для программы 3 |
|
|
||
tri1. = |
m; sec2 ßx -f- tri1 cos2(ßx -f ß2)sec2 ßx + p2 L2tg2 ßb (2.175) |
|||
|
|
(p = ni8 — mt cos 6). |
|
Обозначая в этих уравнениях символами А и В члены, зави сящие от формы наблюденных треугольников, для общего случая
равнбточных измерений получим |
|
|
|
|
|
|
|
|
т1 = Ат1+ В |
|
р"2 . |
|
|
(2.176) |
|||
L |
г |
4 -Ют ^ |
|
|
|
ѵ |
’ |
|
Изменение коэффициентов А |
н В |
в |
зависимости |
от |
формы |
|||
треугольников показано на рис. 47 и |
48, |
для |
фигур |
с |
углами |
|||
ß>30° (лазерные наблюдения ведутся |
на |
зенитных расстояниях |
||||||
до 60°). Но так как при |
современной |
точности |
непосредственных |
измерений (пгг—1—2 м; р = ± 1 —2") влияние второго члена на итоговую погрешность длины хорды значительно больше влияния первого члена, то для более точного определения длины L выгод нее наблюдать низкие спутники в плоскости геоцентрического се чения хорды.
Формулы (2.175) |
не |
учитывают |
падение точности непосред |
ственных измерений |
с |
увеличением |
зенитного расстояния. Если |
же принять во внимание зависимость |
(2.40), то наилучшие условия |
для наблюдений по программе 1 будут определяться соотноше ниями (2.172) и (2.103), при Гі = г2 и е= 0°.
Выгоднейшие условия для определения длины хорды по про грамме 3 создаются при прохождении ИСЗ через зенит второго пункта (на котором определяется направление на ИСЗ а2) на
124
предельно допустимом зенитном расстоянии первого пункта (т. е.
при zmax= 60° и е= 0°).
В этом случае
ßi = 90° — ^znl!X---- ^ — 30° + |
и ßi — 90° + ■ |
Наилучшпе условия для выполнения программы 2 создаются при прохождении 11С3 через середину хорды, при е= 0° и -ßi= ßs= = 45°.
Рис. 47. Номограмма для определения ко эффициентов А и В для программы 1
Рис. 48. Номограмма для определения ко эффициентов А и В для программы 2
В одном цикле наблюдений получить длину хорды с высокой точностью нельзя. Расчет по формулам (2.175) показывает, что для хорды длиной 3000 км при ßi = ß2= 45°, т,=2 м, |.і= 2", ошиб
ка в |
определяемой длине хорды |
по |
первой программе составит |
21 м, |
а по второй и третьей ~30 |
м. |
Поэтому для достижения ре |
125
зультата высокой точности требуются многократные наблюдения. В этой связи может исчезнуть целесообразность раздельного оп ределения длины хорды и ее направления и встанет вопрос о не обходимости совместного уравнивания угловых и линейных изме рений, выполненных для определения геодезического вектора L.
§ 15. УРАВНИВАНИЕ НАБЛЮДЕНИЙ
Рассмотрим общий случай математической обработки резуль татов комбинированных наблюдений по методу условий с дополни тельными неизвестными, когда совместно определяется длина и направление хорды. В качестве дополнительных' неизвестных бу дем отыскивать поправки к приближенным значениям элементов вектора L, соединяющего наблюдательные станции
cos a|j cosA\ |
|
||
(COSO], sin A . |
(2.177) |
||
|
sm ф |
J o |
|
Условные уравнения, отвечающие программам 1—3, |
имеют вид |
||
F i = /yij 4- г2о2— Ь 0Е 0 — 0, |
(2.178) |
||
Foa = (Ql Оа Ea) = |
0, |
(2.179) |
|
Fob — ri sin (ßi + ßa) — Ц sin ß., = 0, |
(2.180) |
||
F3= r\ — r\ — Lo + |
2гоL0cos ß2 = 0, |
(2.181) |
|
где |
|
|
|
sin (ßi f ß 2) = I |
ауХОгІ; |
|
sin ßj = I E0 X ax I:
COS ßo — E0Go.
Уравнения составляются на моменты 0j синхронных наблю дений.
Условные уравнения поправок, получаемые путем линеариза ции уравнений (2.178—2.181), в общем случае имеют вид
а |
и |
-fa . |
щ -1------ '-г-ѵ' |
+ |
a, |
ѵ |
-f a. |
v. -f |
г 1 |
Г 1 |
1 б , |
б, 1 COS б і О |
1 |
r 2 |
r 2 |
1 b , |
6 . |
|
|
v t„ ' b |
b LE L 0 - f |
|
b |
AA0 -f U? = 0, (2.182) |
||
+ cos |
g., |
|
cos |
1|)0 |
||||
|
|
|
|
|
|
|
где а и b — значения частных производных функций F по изме ренным величинам и искомым параметрам, at// = t>(cos б и АА'о=
= A A qCOS фо-
126