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

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

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

геоцентрического сечения хорды, а их максимальный угол раство­ ра у>90°, часть синхронных наблюдений ИСЗ можно и целесооб­ разно выполнить и в центре зоны, не нарушая условия (2.115). Это объясняется тем, что векториальные суммы

^2>1р = Р \ { Ы , ) + Р'2,2(А,+Ѵ) + Р 0>2 ^ і+ 2_j

будут равны нулю, если веса наблюденных точек на краях зоны равны между собой (рі=Р2 =р), а в центре зоны ро= —2cosy/?, так как

cos 2А1+ cos 2 (Лг +

у) — 2 cos у cos (2Ах -f у) •■=0 1

2 }}

sin 2Ах + sin 2 (i4j +

у) — 2 cos у. sin (2Лх + у) = 0 )

 

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

п = Ра

mw, О *

(2.117)

крі

2 cos упі\Ѵ

 

где k = —2cosy.

 

измерений

Третья задача — проектирование дополнительных

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

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

С

2q. (2.118)

1\ѵ 2А ,

Для этого потребуются наблюдения в северной половине зоны (Ир >-0°), а если условия наблюдений данного спутника не позво­ лят этого сделать, то следует обратиться к наблюдениям ИСЗ с другими параметрами орбиты.

Для повышения точности отдельно взятых углов і|з и Л до­ полнительные наблюдения следует проектировать так, чтобы под­ спутниковые точки располагались на изолинии Лр = 90°, для умень­ шения ошибки угла или Ар= 0° — для уменьшения ошибки угла Л.

ПО

Г л а в а 4

ОПРЕДЕЛЕНИЕ ДЛИН ХОРД

§ 12. МЕТОД ПЕРЕСЕЧЕНИЙ

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

ным путем,устанавливая

в створе

 

 

 

измеряемой

линии

промежуточ­

 

 

 

ные пункты. В отдельном слу­

 

 

 

чае

такими

 

пунктами

служат

 

 

 

рефлекторы

 

антенн

передвиж­

 

 

 

ных

станций, устанавливаемых

 

 

 

на некоторой высоте над поверх­

 

 

 

ностью Земли. Новейшие радио­

 

 

 

дальномеры типа РЕМ-2 способ­

 

 

 

ны

измерять,

таким

 

образом,

 

 

 

расстояния до 150 км с ошибкой

 

 

 

меньшей 3-10_G.

 

длин

хорд

 

 

 

Для

определения

 

 

 

от 60 до 400 км в геодезической

 

 

 

практике нашел применение ме­

 

 

 

тод

пересечения

створа.

Этот

 

 

 

метод,

разработанный

вскоре

 

 

 

после войны канадскими учены­

 

 

 

ми [16], предусматривает ис­

 

 

 

пользование

в качестве

проме­

 

 

 

жуточного

пункта

передвижную

 

 

 

приемо-передаточную

 

станцию

Рис.

44. Геометрия

пересечения

(ПС),

устанавливаемую

на

са­

створа

 

 

молете

в

отдельных

случаях

 

створ наземных пунктов

и на ИСЗ), которая должна пересекать

на постоянной высоте Н0 и под некоторым углом а (рис.

44).

Для определения длины хорды Ll2 необходимо по мере при­ ближения передвижной станции к створу равномерно и синхронно (через 1—2 сек) измерять дальности Dи и D2i от наземных стан­

ций до ПС

и составлять суммы измеренных расстояний Si = Du +

+ D2i. Установив по

результатам измерений функциональную за­

висимость

S; = /(/,•)

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

тории /, пройденным ПС от начальной точки отсчета, вычисляют минимальную сумму s0, которая соответствует моменту пересече­ ния ПС створа измеряемой линии.

В процессе полета высоту ПС над поверхностью Земли изме­ ряют периодически с помощью радиовысотомера. В дальнейшем, по известным значениям Dl0, D2о и Н0 вычисляют длину хорды,

111

соединяющей наблюдательные станции или длину соответствую­ щей геодезической линии.

Первые работы [16] по измерению расстояний методом пере­ сечений были невысокой точности. Но последние опытные работы, выполненные с системами Аэродист, Шорам и самолетным радио­ дальномером ЦНИНГАиК [14] показали, что линии длиной 100 и более километров могут быть измерены в благоприятных условиях с точностью порядка 2-н5-10~6 L. В последние годы проводились также опытные работы по определению длин хорд до 4000 км методом пересечений с помощью системы SECOR, с ПС, устанав­ ливаемой на ИСЗ [43]. Результаты этих работ характеризуются мощностью Зч-8-10-1і.

Геометрическую сущность метода пересечений целесообразно рассмотреть на примере кругового маршрута ПС над шарообраз­ ной Землей1 (см. рис. 44). Предположим, что траектория ПС пе­ ресекает геоцентрический створ наблюдательных станций 1 и 2 под некоторым углом а, при постоянной скорости, высоте и на­ правлении полета. Обозначим: С£— начальную точку отсчета рас­ стояний /, пройденных ПС, Іо— расстояние, соответствующее мо­ менту пересечения створа; s, v, е — центральные углы, соответст­ вующие отрезкам Д-, D0 и (/—/0). Из треугольника ОС{1 находим

D\j = (R + Ях)2 + (R + Я0)2 - 2 (R + Нх) (R + Я0) cos Sl. (2.119)

Но так как

cos = cos vx cos e -f- sin vx sin e cos a,

(2.120)

и в силу малости углов е

8 ЯЗ sin 8 =

и - ІО .

COS 8

R + Rо

 

 

(h - /о)2

( 2. 121)

2 (# + Я 0)*

то подставляя последние выражения в (2.119), получим

Du — D\o — 2 (/£— /0) (R + Ях) sin vj cos а + (It — l0f

cos v2,

 

(R + no)

 

(2.122)

где Dl0— наклонная дальность, соответствующая

моменту пере­

сечения створа.

 

Аналогично

 

D\i = Dl о+ 2 — /0) (R 4- Я2) sin ѵ2 cos а -f (/£— /0)2

+ н Cos ѵ2.

 

(R + Ң0)

 

(2.123)

1 Такое представление является приближенным для действительного поле­ та ПС над геоидом или эллипсоидом. Но при современной точности непосредст­ венных измерений решений для сферы пока удовлетворяет практические запросы.

112

После извлечения корня из уравнений (2.122 и 2.123), найдем

Dv = Dio j 1 — - (li~— [R + Hi) sin V! cos X +

D 10

i k - l o Yü -

D\o

D„(- — Do0 f1 -1- 2

‘ 1 '

(R-l-Wi)

V i *

v

cos ѵЛ

(W+ H„)

 

. (2.124)

~ l°^ (R + Я,) sin v0 cos a Dfo

, (/f - /о)2

(fi +

cos vx

Ö-20

(Я + Wo)

 

Если использовать биномиальные ряды для разложения пра­

вых частей уравнений (2.124), то

в общем виде можно

записать

D j ^ D j b + SEjbiklbY;

( / =1, 2) ,

[k = 1 , 2 , 3 , . .

., n),

где Е — коэффициенты разложения.

 

 

 

 

Суммы измеренных дальностей будут равны

 

 

si =

Dio + Do0 +

2 (Elk +

E2k) (II /„)*.

 

(2.125)

Последнее уравнение решает задачу по определению вида

функции s= /(/);

сумма измеренных

дальностей

здесь

аппрок­

симируется степенным полиномом

k-то порядка.

Последующая

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

Заметим, что наибольшее расстояние D, которое можно опре­

делить методом пересечения створа, примерно равно

 

•^тах(км) = 4, \ ( У Н г - Н

+ Ѵ Н 2- Н + 2 ] / Я 0

— Я ), (2.126)

где Н — среднее значение

высоты земной поверхности вдоль из­

меряемой линии в м.

 

 

Накопленный практический опыт показал, что при существую­ щей точности измерений степенной ряд (2.125) достаточно огра­

ничить квадратным трехчленом вида

 

Si = P + q ( U - Іо) + г [lt - l0f ,

(2.127)

который представляет собой уравнение параболы

второго порядка.

Для равноотстоящих значений Іі+І—Ц = Е1

можно записать

/f = ( t - l) A /,

(2.128)

и уравнение параболы привести к виду

 

Pi + qii + r-jY = sh

(2.129)

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

/о, А/, a, DJ0, Dio-

113

Составив по результатам измерений уравнения поправок

Р о + Ч і 4~ r i + Asi — ѵ і

 

P o "г 2£7і + 4/у As2= t>2

(2.130)

 

Po + «Д + «2/'i + As« = °n

где As = s0si :p0 = pl—s0; s0—-меньшая из сумм измеренных рас­ стояний, решают эти уравнения под условием [u2]= min и опреде­ ляют искомые коэффициенты уравнения (2.129).

Отыскание наименьшего значения функции S;=/(/) выполняется

стандартным способом. Так, при — =0

аі

q1+ 2i\i — 0,

откуда

і =

2гI .

 

Тогда

2

 

2

 

 

 

 

smin = D10 + Dw =- Pl -

= s0+

p0 - -p-

(2.131)

 

4rl

 

4rx

 

Для оценки точности измерений вычисляется ошибка единицы веса

(2.132)

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

пи

=

Р2 Qp,4

QqI +

 

 

 

 

4rj

16

Яі

о

I

Qp,r,

(2.133)

----Ѵг,?,

Г\

 

 

 

 

где Q — элементы обратной матрицы коэффициентов нормальных уравнений. Для приближенного определения этой ошибки можно использовать формулу

< “ CJ / - £ ■

(2-134)

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

114

влияния навигационных ошибок на отдельные и минимальные суммы измеренных расстояний.

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

валов времени измерений т определяется формулой

 

 

2

0,26-ІО“ 4 ѵ2 км/час

ТСск2

sina X

 

 

 

 

2L км

 

 

 

 

. .

Г .

/п 2-\-2п +

45\

т±

 

(2.135)

X у (П —

( -

 

2 I

)

 

 

где

(2.136)

Средние квадратические ошибки минимальной суммы, завися­ щие от влияния изменения высоты полета Я0 и отклонений от

заданного направления

(рысканий р),

соответственно

равны

mSH= ( - ^

+ 1,2-10-j L km) — ,

(2.137)

\

Т

 

J уГп

 

mSp = 2,2- IO-4

• sin 2сб

n ^ 1 mp.

(2.138)

Из всех навигационных ошибок наибольшее влияние на значе­ ние минимальной суммы оказывает изменение высоты полета.

Так, например, при

п=21; 1 = 200 км,

Я0 = 3 км, т н =15 м,

mSH = 0,3 м. Если

же

предположить, что передвижная станция

будет установлена

на

ИСЗ, высота полета

которого сопоставима

с длиной порды L, то изменение величины Я за время измерений будет сказываться еще сильнее на искомой минимальной сумме расстояний.

Поэтому в методе пересечений возможно использовать только круговые орбиты ИСЗ или орбиты с весьма малым эксцентриси­ тетом. В последнем случае, при известной траектории спутника, в измеренные суммы дальностей нужно вводить поправки за от­ клонение высоты полета ИСЗ от средней высоты полета Я0 в мо­ мент пересечения.

Для этого можно воспользоваться формулами Б. Б. Серапи-

наса [69].

 

ZS3,

 

 

(2.139)

 

As, = 4 ------ -----Ah,

 

sj 4x°-L°-

 

 

где

 

 

 

 

X « (/ — /0) cos a -{-

2

D-°;

z = Hn

L2

ң

__

ff14-

(R + Hcp) ’

8

1cp

 

2

115

Для уменьшения влияния навигационных ошибок па искомое расстояние измерение линии методом пересечения створа делают многократно, причем чтобы ослабить воздействие возможных си­ стематических ошибок измерений — инструментальных пли обу­ словленных влиянием внешней среды, — маршруты пересечения створа обычно видоизменяют по специально разработанной про­ грамме [16], [26].

§ 13. ВЫЧИСЛЕНИЕ ДЛИНЫ ХОРДЫ

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

мую длину хорды L,

если все

необходимые

данные

для

расчета

1, # 2, Но, Dю, D20)

известны.

Из решения

треугольника

1, 2 С0

(рис. 45) получим

 

 

 

 

 

Lr — D]о + Dio — 2Dl0D.,0 cos (cpj

ф2).

 

(2.140)

Так как

 

 

 

 

 

hi = H0 — H1 = D10cos Фі — 2 (R -f

Я,) sin2

,

(2.141)

TO

 

 

 

 

 

 

ih +

2(R + Hl) sin"- Y

 

 

 

COS ф х --=

(2.142)

Dio

В силу того, что

sin

4’i

 

D

 

 

(2.143)

= ------------- ,

 

 

2

2

(R +

Hi)

 

 

 

 

 

,2

-4- h\ -T

О. hi

D1 о ~ D1 -}- h\ -j- 2D| hi sin —„— = Di

{R +

(2.144)

 

 

2

 

 

Hi)

имеем

(2.145)

Возвращаясь к уравнению (2.143), получим

,2 Ч’і _

öfo-Af

(2.146)

sin*

i ( R + Hi){R + На)

2

 

Подставляя найденное соотношение в (2.142), будем иметьг

COS ф !

=п

hi

I

Dio

Af

(2.147)

Du

 

2 (R + H0)

2D10 (R -j- H0)

 

 

 

 

Аналогично найдем

 

 

 

 

 

COS ф ,

/

I

Djo

h\

(2.148)

= (

-------------------D20

---------------------------

2 (R +- -------------Ho)

2D20 (R + Ho)

 

'

 

 

116

Рассмотрим вопрос о точности определения длины хорды под влиянием ошибок измеренных величин D и Я.

Дифференцируя уравнение (2.140) с учетом (2.141) и глав­ ного члена уравнении (2.147—148), получим после некоторых пре­ образований

dL = Т ' DuAo

2 (R + До)

\D.о — D,0cos (фі -Ь Ф2) 4“

о

 

sin ( < p i + ф о )

 

 

j dD10 + -j-

°20 — D10 COS(фі + Фе) +

sin cp!

Sin ( ф і + ф о )

^10^20

sin (фі

— 1dDiü +

sin Cp2

2 (R + До)

sin cp.,

Sin ( ф і . + ф о )

dHl + Blü.. sin (tpt +

Фа)

sin Cp!

L

sin ф о

 

 

 

 

 

1

( D j o D o p S i n

( ф і

+

 

Dio

I

Dso \

 

 

 

 

 

L {

2 (« +

//„)*

 

Sin Cp!

 

sinCp2 J

 

 

 

 

 

- sin (Фі +

cp,) I

Dso

I

D10

 

dH0.

(2.149)

 

 

 

 

 

 

 

 

 

sin Cp!

sin cp.,

J \

 

 

На основании последней формулы молено получить выралсение

для ■

средней

квадратической

 

 

 

 

 

ошибки искомой. длины хорды.

 

 

 

 

 

Этой же формулой (2.149) мож­

 

 

 

 

 

но

воспользоваться

для

 

приве­

 

 

 

 

 

дения длины найденной хорды к

 

 

 

 

 

фиксированным

точкам

 

земной

 

 

 

 

 

поверхности, учитывая высоты Я

 

 

 

 

 

приемных

 

антенн

наземных

 

 

 

 

 

станций.

 

 

 

 

 

 

что

 

 

 

 

 

 

Принимая во внимание,

 

 

 

 

 

условия измерений,

близкие

к

 

 

 

 

 

действительным,

предполагают

 

 

 

 

 

использование симметричной

или

 

 

 

 

 

примерно

 

симметричной

 

схемы

 

 

 

 

 

измерений

относительно

 

норма­

 

 

 

 

 

ли точки С (см. рис. 45), поло­

 

 

 

 

 

жим для

 

простоты,

что

фі = ср,=

 

 

 

 

 

= Ф

и

2Д э~ Д

Тогда уравнение

Рис.

45.

 

Схема определения длины

(2.149)

примет

вид

 

 

 

 

 

хорды

 

 

 

 

dL =

I sin2 ф -j- cos ф

Г Л,ср

 

 

Дс р

 

 

j (dDl0 -j- dO20) -f-

 

 

Dср

 

2 (R +

H0)

 

 

 

+

cos ф dH2-L cos <pdH2

2

 

D.с р

соэф^Яц.

(2.150)

 

 

 

 

 

 

 

 

 

 

 

(R + До)2

 

 

Заменив в этом уравнении sin ф и соэф их функциями

(2.147),

опустим

малые

члены,

существенно

не

влияющие на

точность

117

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

2 (R + я о)

Тогда

 

 

 

 

 

dL = dD10+ dD20 +

(^ D^ o)- (dH, + dH2) -

-j— щ - dH0

и

 

 

 

 

 

 

D2

 

 

4°cp

 

mi = m(DrhD2)

cp

К / , +

m%t) +

пн.

(* + н 0у

(R + w„)*

 

 

 

 

(2.151)

(2.152)

Высоты Hi и #2 земной поверхности почти всегда могут быть найдены с высокой точностью, и их ошибки, как видно из (2.152), не окажут большого влияния на точность определения длины хор­ ды. Что же касается высоты На задающей станции, то ошибка в ее определении оказывает вдвое большее влияние на точность

окончательного результата

и потому

она

должна определяться,

по возможности, точнее и

особенно

при

измерении длин хорд

большой протяженности; особо это следует учитывать при приме­ нении ИСЗ.

На

основании

(2.149) и (2.152) можно сделать заключение и

о том,

что при

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

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

С этой целью воспользуемся формулой, предложенной А. А. Изотовым [24], для определения величины радиуса-вектора пространственной точки. Тогда

 

H\ = N( 1— е2 sin2 В)

N2 е2sin2 В cos2 В

Hi R,

(2.153)

 

 

2(W+ Hi)

 

 

где

Hi — высота точки над

поверхностью референц-эллипсоида.

Применив биноминальные ряды для преобразования выраже­

ния

(2.153), получим приближенную формулу

для вычисления

высот Н'і.

 

 

 

 

Н'і = Hi 4- (ае - R ) + -~ аее2sin2 В

 

 

-----—а2е4 sin2 В (sin2 В — 4 cos2 В),

(2.154)

 

8

 

 

 

118

обеспечивающую точность вычислений до 10 см (опущены члены порядка е6). Геодезическую широту точки здесь также достаточно

знать приближенно

(с точностью до

 

 

5-10").

 

 

 

 

 

 

 

Вместе с тем нужно заметить, что

 

 

при переходе от нормального сечения

 

 

к центральному, и в том случае, ког­

 

 

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

 

 

меридиана,

при определении высоты

 

 

Но задающей станции может возник­

 

 

нуть погрешность из-за расхождения

 

 

плоскостей

центрального

и нормаль­

 

 

ного

сечения. Эту

погрешность

(на

 

 

рис. 46 дано положение хорды L2 в

 

 

плоскости

первого

вертикала)

при­

 

 

ближенно

можно оценить

формулой

 

 

АЯп

 

(Во-Ф0) 1

sin А

Рис. 46. Расхождение

между

2 1sin2

 

 

 

нормальным

п центральным

 

 

 

 

 

 

сечениями

 

 

 

 

 

 

(2.155)

 

 

где

^ _ D io D-io Sin (ф! + фо)

высота

треугольника

1С02,

А —

 

 

L

 

 

 

 

 

 

азимут хорды.

Как показывает расчет, эта погрешность мала по сравнению с ошибками, которые допускаются сейчас при определении высо­ ты полета (1—2 м), и потому она заметно не скажется на точ­ ности определения длины хорды. Например, при D|o= D20 = 400 км,

Но=12,5 км, sin(срі + фг) =0,06, /4=90° и (ВФ)0:2 = 35О", АН0 —

=0,1 м.

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

определяется по формуле общей арифметической средины.

Г л а в а 5

СОВМЕСТНОЕ ОПРЕДЕЛЕНИЕ ДЛИНЫ

ИНАПРАВЛЕНИЯ ХОРДЫ

§14. КОМБИНИРОВАННЫЕ НАБЛЮДЕНИЯ НА СТАНЦИЯХ

Для совместного определения длины и направления вектора, соединяющего наблюдательные станции, необходимо и достаточно в некоторый момент Ѳ, звездного гринвичского времени одновре­ менно измерить на обеих станциях топоцентрические экватори­

119

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