Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Разумов, О. С. Пространственная геодезическая векторная сеть

.pdf
Скачиваний:
5
Добавлен:
19.10.2023
Размер:
5.69 Mб
Скачать

то переход от предварительных координат каждого пункта урав­ ниваемого хода (вычисленных по результатам непосредственных измерений от исходного пункта) к их уравновешенным значениям осуществляется по таким формулам

Рис. 51. Схема сети из векторных ходов

(2.250)

Если ход близок к равностороннему, то тогда

1

і ;

Y ,= ■ - ¥ ) - fy J -

 

п

п

 

 

Z/ = Z'j f z ~

(2.251)

 

1

п

 

При этом ошибка

единицы веса, опреде­

ляемая по

материалам

уравнивания,

будет

равна

 

 

 

 

 

'у9 ' г-

(2.252)

 

 

3П п

 

Из

рассмотрения

соотношении

(2.246) —(2.251)

видно, что

при соблюдении принципа равного влияния (2.230)

становится пра­

вомерным

и строгим

раздельное

распределение

невязок

в при­

ращениях

координат

сети. Это

обстоятельство

значительно

облегчает задачу уравнивания и допускает возможность

приме­

нения

способа полигонов В. В. Попова [52].

 

 

 

Так, если придать каждому ходу векторной сети, изображен­

ной на

рис. 51, вес, равный р , = ---------- , то можно

будет,

поль­

з ъ зуясь правилами способа полигонов, написать по чертежу сети

систему нормальных уравнений. Для случая, изображенного иа рис. 51, эта система предстанет в виде

и т. д.

140

Решение этой системы дает соотношение

 

, Кг* \

Q u

Q l 2 Q13 • • • Q 15 \ / f ,A

 

 

I Kir

'

Q l 2

Q

2 2

Q - 2 3

. Q 2 5

\ / „ A

(2.254)

 

к ях

 

Q i 3

Q

2 3

Q 33

- Q

35

/ . I .

 

\ к , : )

Q i B Q 25 Q 3 5 • ■ . Q 5 5 / \ K ,

 

где квадратная матрица Q есть обратная

матрица

системы нор­

мальных уравнений

коррелат.

 

 

 

 

 

и другие невязки Uy и

Аналогичным

образом

распределяются

}іг в полигонах сети с одной

н

той

же

системой

коэффициентов

в нормальных уравнениях коррелат

(2.253).

 

 

 

Поправки и приращениям координат по отдельным ходам сети

определяются на основании

(2.253)

по формулам

 

 

^ [ Д а ] , — [ М х і I ^ 1 л - >

 

И [ Д у ] , —

[ P z l 1

 

 

У [ Д г ] ,

 

 

|

 

Ѵ[&Х]. ~

[Pzla

 

^2л-)>

У[Ду]„ =

[Рд]г

 

^2J')> I

 

 

 

 

 

- * *

. )

 

 

 

[ '

<2 -255>

~

X lJ’

 

и[й!/1,

 

lW .ll ^ 23"

У[й2],

=

и ч із

^32 j

 

а поправки на отдельную сторону хода 1

составят

 

 

 

^Д.с(. =

f j |

К і х ,

ѴА х . =

fl[

Кіу,

ѴА г . =

_Кхг.

 

Переход от поправок приращений координат к поправкам непо­ средственно измеренных величин можно осуществить по формулам, аналогичным (2.248), а именно

vL. =

р|

(cos і|); cos Aj I\u. -f cos ip; sin Ai Kly + sin іід/фи)

 

 

Оф. =

P"

РІ (— sin ljj(- COS Л; Klx sin l|); sin A;ATly -j- COS

K\z)

 

I

 

Li

 

 

 

 

 

 

n,.

 

 

 

%

=

P" Lt сбзЧТ (~

Sin Л/ Kl r + C0S Лі

(2.256)

 

 

 

 

 

Ошибка единицы веса, определяемая по результатам

уравнива-

ния,

будет равна

 

 

 

 

 

 

2

\Pl ѵі \ + [р'і.% I + [ р > л ]

(2.257)

 

 

 

і-1 =

з7

 

 

 

 

 

где г — число полигонов в сети.

141

В соответствии с (2.254) можно будет написать и общую для всей сети матрицу Q весовых коэффициентов итоговой (утроен­ ной) системы нормальных уравнений коррелат. Элементы этой матрицы расположатся симметрично друг относительно друга и будут равняться соответствующим элементам матрицы (2.254), причем

<2.ѴЛ = Q yt --= 0 .гі2 і = Qii

Qxixj = QuiUj = Qzi2j = Qj

(2.258)

 

Q x y = Qx-Z = Q yz = 0

 

Имея полную матрицу Q, вес любой функции уравновешенных

элементов системы можно вычислить по формуле

 

 

X

 

(2.259)

где F — частные производные искомой функции по результатам

измерений.

 

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

В этом способе, после вычисления предварительных координат пунктов сети, уравнения поправок к измеренным величинам име­ ют тот же вид, что и приведенные для уравнивания космической триангуляции (1.79—1.81).

Чтобы упростить вид уравнений поправок и соответствующих им нормальных уравнений, системе уравнений поправок измерен­ ных величин противопоставим эквивалентную систему

cos

cos A°k(ѵхк vxi) -f cos ф!/;sin A°ik\vyk — vyi) ■

 

 

+ sin

(vzk -

vzi) + h ik = eL[k,

PL - ~

 

 

— sin f ik cos A°k(vxk — vxi) — sin

sin A°k (oyb vyi)

 

cos v ik{vzk— vzi) +

 

C

= 8

 

cp-

 

L{ 4 ik

 

(2.260)

 

 

 

ill

 

'Ьь’

p*'~

Ц ,» ’

 

 

 

 

 

 

 

 

4)

 

• sin A ° k

(vxk vxi) -f cos A °u (V

k U

i ) +

L°ik cos \pik ik

 

 

 

 

 

 

 

yU

 

 

 

 

=

 

P*

 

cp1

 

 

 

 

 

^ifcCOS21|Hk ml

 

 

 

 

 

 

 

 

142

Рис. 52. Геодези­ ческий векторный четырехугольник

решаемую под условием

 

 

 

 

L- cos'-1|

+

А

min.

(2.261)

 

 

.

 

J

 

В такой записи коэффициенты при неизвестных в уравнениях

поправок представлены направляющими

косинусами градиентов

измеренных величин, свободные

члены — уклонениями

прибли­

женных позиций точек от соответствующих поверхностен положе­ ния, а веса — весами поверхностей положения.

Обозначив в уравнениях поправок коэффициенты при неизвест­

ных буквами а,

b и с,

систему

(2.260)

запишем так:

a U ІРхк Ѵх і ) +

Ь ц {Ѵук Vyi) +

с ц (vzk v zi) +

I ц

== 8и

a M (Vx/i

— 1

 

) +

b Aj

(

Vyk

Vyi)

+

l.\f

■= 8д/

(2.262)

Ux i

 

 

 

 

 

 

Ö1|V (Vx i

Ux i ) byV (Vyk

Vyi)

+

Cy,j (Vzk

Vzj) -\- ly,j

где ] — индекс стороны векторной сети.

 

 

 

В том

случае,

когда

точность

результатов

 

непосредственных

измерений

 

согласована

по

 

принципу равенства

 

влияний

(2.230),

и вследст­

 

вие этого Р2 = р'^=р'х = р,

матрица

нормальных

 

уравнений приобретает квазидцагональный вид. В этой матрице квадратичные коэффициенты при неизвестных будут равны сумме весов [р] сторон векторной сети, сходящихся в данном пункте, коэффициенты типа — [а,-аА,, — [bjbj

окажутся равными весу р, стороны, взятому со знаком минус, а прочие неквадратичиые коэффи­ циенты [аД;], [öjcj, [bjCj] будут равны нулю.

Это обстоятельство позволяет написать матри­ цу коэффициентов нормальных уравнений прямо по чертежу сети (рис. 52).

В силу равенства и полной симметрии коэффициентов при неизвестных ., Ѵу., ѵг. , а также из-за того, что корреляционные

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

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

/(Рі + P i ~+ Р з )

— P -г

—1

— Р з

Q' =

— Ра

(Ра + Р з + Р6)

— Рз

V

~ Р з

Р з '

(Р з + Р і + Рй ) і

ИЗ

и искомые неизвестные определить из соотношении

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

[Pl 41 + К 8ф~] +ІРлел]

(2.263)

п k

 

Найденные весовые коэффициенты

как составляющие

элементы обратной матрицы нормальных уравнений, в итоговой сводной таблице расположатся симметрично, в силу равенств

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

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

При использовании радиодалыюмерных систем итоговую си­ стематическую ошибку в длинах сторон можно представить как сумму двух слагаемых

sL — oL-f- ЯВ,

(2.265)

из которых первое не зависит от величины измеряемого расстоя­ ния, а второе связано с коэффициентом систематического влия­ ния Я; учет таких ошибок выполняется исходя из характера их воздействия на результаты измерений (в зависимости от типа

144

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

В общем случае учет систематических ошибок sL может быть осуществлен путем введения дополнительных неизвестных в ус­ ловные уравнения поправок или уравнения ошибок измерений.

Так, если у нас L j — измеренная длина хорды, sL— ее система­

тическая ошибка, 6j — случайная ошибка наблюдения, L\

— вычис­

ленная по предварительным координатам пунктов

длина

хорды,

dLj — изменение длины линии в зависимости от изменения

коор­

динат ее концов,

то для окончательного значения

длины хорды Т'.

можно написать уравнение

 

 

 

 

Lj = Lj Si -}- еj — Lj + dLj,

 

 

(2.266)

откуда

 

 

 

 

 

ej — Sl T dLj + {Lj Lj).

 

 

(2.267)

Соотношение

(2.267) представляет собой уравнение

поправок

сучетом систематической ошибки sL.

Взависимости от характера влияния этой ошибки на измеряе­ мое расстояние в уравнения поправок может быть включено то

или иное слагаемое из (2.265), или, в общем случае, оба.

 

С учетом

равенства (2.267)

первое уравнение поправок си­

стемы (2.260)

приобретает вид

 

 

cos ф°*cos A°ik(ухк — ѵхі) +

cos фі* sin A ik(vyk — vyi) +

 

+ sin ф°* (vzk vzi) + oL — ILik + h ik = e,Lik.

(2.268)

Если при измерении сторон сети использовалась разная

аппа­

ратура, со свойственными только ей систематическими ошибками уравнения поправок (2.268) будут содержать в каждом отдельном

случае свои значения sL.

измерений

по способу условий

При уравнивании результатов

в условные уравнения поправок

вводятся те

же дополнительные

неизвестные а и А. с коэффициентами, отражающими влияние этих неизвестных на невязки условных уравнений.

’ Так, если исправленное значение длины хорды найти по фор­

муле

 

Lj ==Lj CTL LLj,

 

то приращения координат между пунктами сети

определяются

из уравнений

 

Ax'. = AXj aLcos фу cos Ay АДдгу

 

Ay'. = At/у oLcos фуsin Ay KAt/j ..

(2.269)

Аг. = Azj cLsin фу AAzy

 

145

Поэтому в условные уравнения поправок (2.238) или (2.241) нужно ввести дополнительные неизвестные aL и А с теми же ко­ эффициентами (или суммой коэффициентов), что и в соотношениях

(2.269).

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

Что касается учета ошибок положения исходных пунктов при уравнивании векторных сетей, то он может быть выполнен по схемам (1.167) или (1.168).

§ 18. ПРОЕКТЫ ВЕКТОРНЫХ СЕТЕЙ

Рассматривая представленные и осуществляемые на практике проекты построения пространственных опорных геодезических се­ тей, можно сделать несколько замечаний.

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

Так как точность определения направлений хорд земной поверх­ ности при надлежащем выборе орбит наблюдаемых спутников (или

ракет) практически не зависит

от

длины этих хорд,

а короткие

расстояния между пунктами земной поверхности пока

определя­

ются точнее, чем длинные, то

при

многоступенчатой

системе по­

строения опорной сети возникает опасность искажения сетей сгу­ щения за счет ошибок сети первого порядка.

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

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

личие

одного пункта

на участке

поверхности Земли

размером

5X5°,

а где позволяют условия, то и 1X Г '.

переходу

2.

Наблюдаемая

в настоящее

время тенденция к

только на радиотехнические средства наблюдения ИСЗ, обеспечи­ вающие проведение работ* независимо от погодных условий, ка­ жется нежелательной, так как при отсутствии фотографических наблюдений спутников мы одновременно отказываемся и от опре­ деления направлений сторон в пространстве, т. е. от усиления ори­ ентировки сети.

1 Один из первых проектов космической триангуляции США предусматри­ вал создание такой густой сети спутниковых пунктов на территории США. Реа­ лизация этого проекта в настоящее время задержана.

146

Лучшим, наиболее надежным и гибким методом построения глобальных геодезических сетей является векторный метод, кото­ рый, по мере возможности, и следует применять.

В 1969 г. И. Д. Жоиголович представил первый проект геоде­ зического векторного хода, связывающего Арктику и Антарктику

118J. Этот ход, соединяющий пункты Баренбцург

(о. Шпицбер­

ген)— Звенигород — Каир — Могадишо

— Реюньон

— Мирный,

имеет общую протяженность 17 000 км со средними

длинами хорд

2800 км.

ф, к), соединяющих на­

В этом ходе элементы векторов L(L,

блюдательные станции, предполагается получить по результатам синхронных наблюдений ИСЗ с точностью 2 -10-6. По предвари­ тельным расчетам И. Д. Жонголовича для определения с задан­ ной точностью длины и направления каждого вектора L потре­ буется около 20—30 пар синхронных наблюдений ИСЗ, причем точность определения топоцентрических координат должна быть не ниже 1" по направлению и 1—2 м по дальности. При успешном выполнении этих работ длина стягивающей хорды L=12 400 км будет получена с точностью 10~6. В этой связи автор проекта ука­ зывает, что никакие триангуляционные ряды не смогут обеспечить указанную точность для столь большой длины [18]. Результаты таких измерений могли бы послужить базой и для уточнения мас­ штаба солнечной системы (АЕ). В проекте [18] отмечается полез­ ность проложения такого же хода в западном полушарии для по­ лучения замкнутого мериднального кольца и целесообразность соединения векторного хода с экваториальной сетью станций на­ блюдения SAO и другими геодезическими сетями.

Работы по проложению этого хода, начатые в рамках между­ народной программы 1SAGEX, в настоящее время продолжаются.

При рассмотрении возможности сгущения сети спутниковых пунктов на территории СССР возникает мысль о целесообразно­ сти создания такого геодезического построения, которое могло бы одновременно с астрономо-геодезической сетью служить главной геодезической основой страны.

Такой основой могла бы стать континентальная пространст­ венная векторная сеть, построенная в форме геодезических четы­

рехугольников

с вершинами в

узловых пунктах

пересечения

рядов триангуляции

I класса.

 

 

Если в этой

сети

длины хорд

измерить, с ошибками порядка

3 -10~бч-2 ■10_б,

а направления

хорд определить

с точностью

0,5—0,7", то мы получим высокоточную геодезическую основу для уравнивания сетей сгущения и для решения всех других задач геодезии. Главным достоинством такой сети будет ее однородная точность и абсолютная ориентировка в пространстве.

Построение опорной сети по предлагаемой схеме позволяет в кратчайшие сроки решить задачу по распространению единой системы координат на территории СССР, а в случае необходимо­ сти,— и по оперативному обновлению этой сети. Состав наблюде-

147

ыип для определения направления н длины отдельно взятой хорды

в этой сети таков, что

при наличии соответствующей аппаратуры

и погодных условий, он

может быть выполнен в течение нескольких

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

Практическое определение длин сторон в этой сети можно будет осуществить методом пересечений створа. Что касается оп­ ределения направлений сторон, то для этой цели придается выпол­ нять синхронные наблюдения вспышек ламп, установленных на ракетах или шарах-пилотах, которые можно будет запускать в нужные моменты н в заданном направлении. Использование ИСЗ для этих целей будет затруднительным, так как наблюдение высо­ ких ИСЗ не обеспечит необходимую точность из-за малого парал­ лактического угла, под которым будет виден спутник, а низкие спутники трудно наблюдать, так как они обладают слишком боль­

шой угловой скоростью движения.

Опытные работы по определению направлений коротких хорд, выполненные в Финляндии и Чехословакии, дали обнадеживаю­ щие результаты. Работы Вяйселя по определению направлений

хорды Турку — Хельсинки

по 26 парам синхронных

наблюдений

ПВЦ, поднимавшихся на

высоты 17—22

км, дали

результат

с точностью ±1,2" [871. Новые опыты [88] подтвердили

возмож­

ность определения

топоцентрического

направления

 

на

вспышку

лампы, установленную на шаре-зонде,

с точностью

1".

 

рас­

чет

В предварительном порядке нами выполнен сравнительный

точности пространственной векторной

сети,

состоящей

из

16

геодезических

четырехугольников

со сторонами

в

200 км, и

аналогичной сети, состоящей из рядов государственной триангу­

ляции I

класса.

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

Модель опорной сети,

I класса,

исследовалась

в ЦНИИГАиКе А. Ш. Татевяном [72].

По материалам этих исследований найдены ожидаемые погрешно­

сти положения узловых пунктов и уравновешенных значении

длин

звеньев сети (цг = 0,8 м). Аналогичный расчет

выполнен

и для

пространственной векторной сети ((.іг. = 0,7 м).

Результаты

этих

вычислений представлены на рис. 53, из которого видны заметные преимущества новой системы.

Сделаем для примера предварительный расчет состава наблю­ дений на станциях для построения векторной сети, изображенной на рис. 53. При выполнении этого расчета будем считать, что точность определения топоцентрического направления в зените со­ ставляет цо=1>5", а необходимая точность направления хорды

то = 0,7".

Предположим вначале, что направления хорд в сети (L = 200—

240 км)

определяются автономно. Тогда, для достижения равно­

точного

результата т^ = т к cos г|/= т 0 потребуется применить та­

кую схему наблюдений, в которой углы у между плоскостями син­ хронизации были бы не менее 90°.

148

Поставив это предварительное

условие, выполним расчет по

формулам (2.103—2.108)

и найдем, что для симметричной

схемы

наблюдений

со станций,

когда е= 45° (у= 2е = 90°)

и 2тах=70°,

ракеты или

шары-пилоты

следует

запускать на

высоту

Н =

= 45—50 км, справа и слева от хорды на удалении 60—70 км от ее середины.

Рис. 53. Сравнительная точность семейства полигонов астрономо-геоде­ зическом сети СССР и векторной сети (в кругах показаны ошибки положения пунктов, по линиям сети— ошибки их взаимного положения)

 

При поставленных условиях параллактический угол при ПВЦ

ßc— 120°,

а ошибка m~w

углового положения

плоскости

синхро­

низации

будет

определяться

формулой,

согласованной

с

(2.39)

и

(2.40)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 cos z sin2 —

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Подставляя

численные

значения, находим,

что,

здесь

т\ѵ=,

= 4".

 

 

 

 

точности направления

хорды

іщ —

=

Для достижения заданной

0,7", суммарный вес всех

наблюдений

должен

определяться

соотношением

(2.115).

 

 

 

 

 

находим,

что [р] =

 

Полагая в нашем случае p = cjm\v = 16/4= 4,

32:0,5 = 64. Иначе говоря,

для достижения

заданной

точности /?г0

=

0,7" потребуется выполнить 16 синхронных

наблюдений

ПВЦ

(по 8 наблюдений с одной и другой стороны хорды).

 

хорды на

 

При раздельном определении направлений каждой

участке,

изображенном

на

рис. 53, для

отыскания

направлений

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

149

Соседние файлы в папке книги из ГПНТБ