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

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

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

хорды. Продифференцировав равенство (2.81) по переменным ф0, Ло II w и перейдя к конечным приращениям, найдем

— sin Wj Awj = (— cos a.\v sin ф>0 cos Л0 — cos ßu--\ sin фп sin Л0 +

+ cos

cos ф0) Дф0 — (cos aw. sin A0 — cos

 

cos A0) cos ф0 AA0.

(2.82)

Члены

уравнения

(2.82),

заключенные

в скобки, равны соот­

ветственно

синусу II

косинусу

угла

Ар

отстояния

плоскости

син­

хронизации Рj

от часового круга

(рис. 38).

А так как с точностью

 

 

 

 

 

 

 

до

малых

величин первого

 

 

 

 

 

 

 

порядка

можно

полагать,

 

 

 

 

 

 

 

ЧТО Sin 01=1, то

 

 

 

 

 

 

 

 

 

 

— AWj = sin А,,. Аф0 —

 

 

 

 

 

 

 

 

— cos ЛрААрСоэфр.

(2.83)

 

 

 

 

 

 

 

 

Перейдя

от равенства

 

 

 

 

 

 

 

(2.81)

к уравнениям попра­

 

 

 

 

 

 

 

вок,

с учетом

равенства

 

 

 

 

 

 

 

(2.83), получим

 

 

 

 

 

 

 

 

 

sin Ар.ѵу — cos Ар~0\' +

IWj =

 

 

 

 

 

 

 

 

= бу с весом Pj

С

 

 

 

 

 

 

 

 

Шш

Р и с 38.

Схема

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

 

 

 

 

 

 

 

 

*7

 

хорды по Вянсялю

 

 

 

 

 

 

 

 

(2.84)

 

 

 

1\ѵ, =

W j — 90°,

 

ѵ л , = ѵ А cos ф0,

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

еі

=

Г д\Ѵ \

, / д\Ѵ \

 

,

/ д\Ѵ \

,

/ д\Ѵ \

 

 

U

r b

Н т г . А '

+

( ч : ) . ”*' +

( ^ : 1 ”-

 

3. Решение уравнений поправок под условием [ре2]= тіп при­ водит к системе нормальных уравнений вида

[psm^Aplv^ — [psinApcosA/;]oA. + [ps\nAplw] = 0

1

Ip sin Apcos A p\ Ü—f- [p cosM/;] uK, — [pcosAp/lv,] =

J, (Z.OÖ)

0 j

из которых находят поправки к приближенным значениям сфери­ ческих координат хорды.

4. По материалам уравнивания определяется ошибка единицы веса

пг = +

/

[ре2]

(2.86)

Уп — 2

100

и тензор ошибок искомых величин

к .

а- = ,щ:

/ЛфЛ' 1 _

 

 

кЩ\' тІ'

 

т-

___ f[pcos2A p] [р sin Ар cos А р]

(2.87)

[P/Pfcsin- (Лр.

у sin Ар cos А р] sin2 Ар]

 

где

 

 

 

 

тХ' — /ил cos2 л\і.

 

Если бы имелся только минимум наблюдений (две плоскости синхронизации), то средние квадратические ошибки сферических

координат хорды на основании равенства

(2.87)

определялись бы

по формулам (2.31) пли (2.32).

 

 

 

 

 

 

Характеризуя погрешность направления хорды в проекции на

сферу эллипсом ошибок, получим

 

 

 

 

 

 

tg2<|’0 =

2 inЧ>Л'

[р sin 2Ар\

 

 

( 2. 88)

 

 

"Л'

1Р cos 2Ap]

 

 

 

 

 

 

 

 

 

 

где сро— угол

между

большой

полуосью

а эллипса ошибок

и ча­

совым кругом,

а

 

 

 

 

 

 

 

 

 

а2

у

[т\ +

т \’ +

Y

т -м)2 +

 

4 ("4л')2) .

(2.89)

б2"

 

 

 

 

 

 

 

 

 

 

 

Согласно исследованию Н. Г. Келл я [28]

 

 

 

т 2(

+

д )

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.90)

 

^

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2q

 

т 2 (

■------LтІ )\ = [Рcos 2 (Ар — фо)] =

 

 

 

\

>-

а2

J

 

 

 

 

 

 

откуда

^Ь'2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2т2

- cos 2

 

 

 

 

 

 

= ГР] —

 

 

 

 

 

 

 

 

 

Поэтому

 

 

= 2 [рsin2 (АР~ср0)].

 

 

(2.91)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[р sin- (ѵ4р

Фо)1

 

 

 

(2.92)

и аналогично

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b2 =

 

 

 

 

 

(2.93)

 

 

 

 

[р cos2 (Ар — сро)]

 

 

 

 

Эллипс ошибок

обратится

в круг

радиуса

/пЕ 2

 

 

— ■ три условии

 

 

 

 

cos 2.Ар]

V

[р]

 

 

 

 

 

= 0.

 

 

(2.94)

 

 

ГрЧ

= (\ [р sin 2А р]

 

 

 

 

 

 

 

 

101

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

тер влияния этих ошибок на направление вектора я) показан на рис. 39.

Если за малый отрезок времени ДѲ спутник С] переместится в точку С'і, то изменение направления вектора я) относительно за-

Рпс. 39. Влияние сдвига спутника на элементы его топоцентрнческого век­ тора

фиксированной синхронной плоскости Pj (т. е. сдвиг линии поло­ жения La, вектора ИД на рис. 35) окажется равным

Ааг = A61sin/4Ci — Aa1cos61cos/4C|,

(2.95)

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

cos ACl =■■cos a\V sin tx— cos

cos tx.

(2.96)

Для вычисления предполагаемого сдвига

спутника за

время

ДѲ необходимо будет знать значения производных топоцентрических координат ГІСЗ по времени. Эти величины могут быть най­ дены или эмпирически — путем соотнесения наблюденных измене­ ний координат к соответствующим отрезкам времени, или анали­ тически, если будут известны на моменты наблюдений элементы оскулирующей орбиты ИСЗ.

Тогда в окончательном виде

 

Ааі= {(Л -) sinАс‘ ~ (lO t)cosöicosЛс,1 АѲъ

^2-97^

а влияние

ошибок регистрации моментов наблюдений на невязки

уравнений

(2.84) получим, подставив (2.97) в (2.38).

 

sin ß2

А \Ѵ в.

' sin (ßi + ß2) sin ßi

sin (ßi + ßo)

д&г

sin Ас,

( lit) C0S 6l C0S Ac■}ЛѲі—

дѲ,

 

dè2

sin Ar

dQj J

cosS2cos Лс,]ДѲ2- (2.98)

00/

 

 

I

Обозначив коэффициенты при ДѲі и ДѲ2 через С{ и С2, под­ ставим (2.98) в (2.84) и получим новый вид уравнений поправок

102

с дополнительными неизвестными — систематическими

ошибками

аі и 02 регистрации моментов наблюдений на станциях

 

sin АрѴу — cos Ap.viV—c^-Oj + c2ja.2+ l\V — ej

(2.99)

с весом

 

Кроме рассмотренных нами способов уравнивания наблюдений, Е. Г. Бойко в [7] было предложено уравнивать измерения по ме­ тоду условий, где в качестве условия выдвигается требование компланарности всех векторов Wj в данной серии наблюдений. В этом методе для каждой тройки векторов составляются услов­ ные уравнения

(IP'i 1Ѵ2 IKg) =

О

( 2. 100)

(W2W3W4) -

О

 

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

§ П. ПРОЕКТИРОВАНИЕ НАБЛЮДЕНИЙ

Задача проектирования оптимальных наблюдений ПВЦ для определения направления хорды связана с необходимостью соблю­ дения ряда своеобразных, а подчас и противоречивых, требований. Но если при создании векторной сети не ставится каких-либо особых задач, то главными требованиями при организации наблю­ дений следует считать равноточность определения искомых направ­ ляющих углов ф и Л хорды при заданной точности т0 окончатель­ ного результата.

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

Ранее теоретически было доказано, что выгодно наблюдать ПВЦ на больших зенитных расстояниях, однако опыт показал, что вследствие искажающего влияния рефракции точность изме­

ренных направлений на ПВЦ падает

с увеличением зенитного

расстояния, и при z > 70° наблюдать

такие цели нежелательно.

 

ЮЗ

Следовательно, при проектировании наблюдений нужно искать компромиссное решение.

Как следует из соотношений (2.32), для каждой хорды суще­ ствуют две особые позиции синхронных плоскостей (0 = Ф и ij = = 90°), при которых направляющие углы хорды ф и Л определя­ ются с наивысшей точностью. Но осуществить наблюдения ПВЦ в этих позициях не всегда возможно. Известно также, что равноточность определения направляющих углов может быть достигну­ та только тогда, когда плоскости синхронизации пересекаются под утлом у^90°. А чтобы обеспечить гарантированную возмож­

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

 

шим (110—120°).

высоту

 

ИСЗ

 

Оптимальную

 

 

над поверхностью Земли некото­

 

рые

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

 

определяют

 

под условием минимума диспер­

 

сии

(2.45)

направления

хорды.

 

Но

такое

решение

не

всегда

 

обеспечивает

получение

равно­

Рмс. 40. Геометрия синхронных

точных значений

углов

т|)

и Л.

наблюдений

По-впднмому,

чтобы

достигнуть

 

выполнения

этого

главного

тре­

бования, оптимальную высоту ИСЗ нужно рассчитать из условия построения веера плоскостей с размахом у = 2е= от ПО до 120°, при наблюдении его на границе зоны видимости на предельно допустимых зенитных расстояниях.

На рис. 40 показана принципиальная схема синхронных наблю­ дений спутника С с концевых пунктов хорды 1 и 2 шарообразной Земли. На схеме приняты следующие обозначения: 102 — плос­ кость геоцентрического сечения хорды, 12С— плоскость синхрони­ зации, R — радиус Земли, Я — высота спутника над поверхностью Земли, zt — зенитное расстояние спутника в точке 1, s — горизон­ тальная угловая дальность до спутника, а — дуга большого круга, стягивающего хорду L12, ßi,2, с — углы в треугольнике 12С, е —• угол между плоскостью синхронизации и плоскостью геоцентри­

ческого сечения хорды.

между

Элементы геоцентрического тетраэдра 12 СО связаны

собой такими соотношениями:

 

длина хорды

 

LV1=-■2R sin - j - ,

(2.101)

зенитное расстояние спутника в пункте наблюдения 1

 

cos zx = cos вcos — sin j?x — sin — cos ßlf

(2,102)

104

или при симметричной схеме наблюдения со станции

cos 2, =

cos е cos — cos

2

— sin — sin

2

(2.103)

1

2

2

4

'

горизонтальная дальность s до спутника при заданном значении z

 

 

sin (2 ■— s) —

R

sin 2,

 

(2.104)

 

 

R + H

 

 

 

 

 

 

 

 

зенитное расстояние спутника при заданном s

 

 

 

 

tgz =

(R +

Н) sin s

 

(2.105)

 

 

(R + Н) cos s R

 

 

 

 

 

 

дальность г1С до спутника

 

 

 

 

 

гіс

{ЯЛ-И) г 2

 

4 (R + H).R sin2 y

+ H2’

(2.106)

 

L\c- + Я2 =

высота спутника

над поверхностью Земли при ßi = ß2= ß

 

Н

= 1А 5

т~ sec2 ß 4- RL12sec ß cos 2 — R

(2.107)

 

Ч—

или

 

 

 

 

 

 

 

H ~ R ^ \ / r * + sin2 -^-sec2ß Ч- 2 s i n s e c ß cos 2

— l ^.

(2.108)

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

на предельных зенитных расстояниях г0 ^s=

е= 0°^

Я0 ------- ----------------R.

(2.109)

5іп(г° ~ т )

 

Оптимальная высота ИСЗ для заданных

значений г, е и L

(а) может быть найдена по формулам (2.102) и (2.108); причем для решения трансцедентного уравнения (2.102) относительно ß

целесообразно воспользоваться

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

sin ß

1 - tg 2

cos ß =

: +tg2--—

 

^ е 2

 

105

и привести его к виду

 

 

 

( cos г — sin-y-j X2 — 2 cose cos у

л'-f- f cosz + sin у Л =

0,

(2.110)

где

 

 

 

х--= tg

р

 

 

Па рис. 41 показан график изменения оптимальной высоты

ПВЦ (ИСЗ) в зависимости от длины хорды при z0 = 70°

и размахе

веера плоскостей у = 2е=120°.

 

 

 

В конкретных случаях практики для проектирования наблюде­

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

оценки

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

мамп

характеристик

z0,

т\ѵ

Ар. (ф)

подспутниковых точек, со-

 

 

 

 

 

ставленных для наблюдения того или

 

 

 

 

 

иного ИСЗ.

 

картограмм

(см.

 

 

 

 

 

 

Расчет

таких

 

 

 

 

 

рис. 42 и 43) должен быть выполнен

 

 

 

 

 

заранее для всей зоны наблюдений

 

 

 

 

 

спутника,

ограниченной

линиями s0

 

 

 

 

 

предельных

горизонтальных дально­

 

 

 

 

 

стей до подспутниковых точек, опре­

 

 

 

 

 

деляемых

по средней высоте Н

ИСЗ

 

 

 

 

 

и максимальному зенитному расстоя­

 

 

 

 

 

нию г0.

Для

построения

изолиний

Рис.

41.

Зависимость

опти­ т\,

и Ар.

следует

выбрать в преде­

мальной

высоты ИСЗ

от

лах

зоны

видимости 11СЗ некоторое

 

длины хорды

 

 

 

 

 

 

 

число

равномерно

расположенных

 

 

 

 

 

подспутниковых точек и

рассчитать

для каждой из них (по приближенным кординатам)

численные

значения этих характеристик.

 

 

 

 

 

 

 

При этом ошибки р измеренных направлений на ИСЗ с на­ блюдательных станций должны быть известны a priori, а в случае большого диапазона зенитных расстояний ИСЗ в пределах зоны наблюдений они должны быть подсчитаны по формуле (2.40). В дальнейшем расчете дальности до ИСЗ определяются по фор­

муле (2.106),

а углы ß

— графически.

Значения

углов Лр

отстояния плоскостей синхронизации от

часового круга могут

быть найдены двумя способами. В первом

из них для каждой расчетной подспутниковой точки составляются уравнения плоскостей синхронизации

( R - R i) ( R . - R i ) (R* — Ді) = Ах + By + Cz + D = 0

(2.111)

и вычисляются углы наклона этих плоскостей к экватору

 

С

(2.112)

COS l j =

(42+ БЧ -С 2)'7“ ’

106

затем по формуле /4„ .=

cosі .

 

= arcsin ------определяют-

cos

i|>„

 

ся искомые

значения

этих

углов.

 

 

Второй

способ преду­

сматривает

вычисление

по

прибл ижениым коордииа-

там расчетных точек векто­ ров Wj с последующим оп­

ределением

углов

A Vj

по

формулам

(2.82)

способа

Вяйсяля.

 

 

 

 

На рис. 42 и 43 пред­

ставлены

две схематичес­

кие картограммы

характе­

ристик

Sq,

гп\ѵ , Ар под­

спутниковых

точек

для

круговой

 

орбиты

ИСЗ

«Эксплорер»

(а= 12 230 км)

при наблюдениях его с кон­ цевых пунктов хорды Ба­

мако—Ашхабад

(фи = 22°,

Z m , i x

= 75°).

Причем

на

пер­

вой

картограмме

 

(см.

рис.

42)

изолинии т\ѵ под­

считаны

при

условии,

что

дисперсии

измеренных

на­

правлений

на

спутник

со

станций А и Б равны еди­ нице и не зависят от зенит­ ного расстояния. Изолинии

т\Ѵ. второй картограммы

(см. рис. 43) вычислены с учетом падения точности измеренных направлений (2.40) по мере роста зе­ нитного расстояния. Из со­ поставления рисунков сле­ дует, что здесь па величину ошибки т \ заметно вли­

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

Рис. 42. Схематическая картограмма ха­ рактеристик подспутниковых точек для круговой орбиты ИСЗ

Рис. 43. Схематическая картограмма ха­ рактеристик подспутниковых точек, постро­ енная с учетом падения точности измеряе­ мых направлений на ИСЗ с ростом зенит­

ного расстояния

107-

Приведем несколько примеров практического использования этих картограмм.

Из рассмотрения рис. 42 можно заключить, во-первых, что по всей области наблюдении спутника значение ошибки /?г21Г почти одинаково и составляет примерно 3,0 ц2. Исходя из этого, можно сделать вывод, что при лучшей геометрической схеме наблюдений, когда углы между плоскостями синхронизации равны 90°, мы по­

лучим среднее значение ошибки направления

хорды (при необ­

ходимом минимуме измерений)

не лучше, чем

±1,7 ц. II, во-вто­

рых, если бы орбита спутника

имела наклон

к экватору около

22°, то при наблюдениях НСЗ с концевых пунктов хорды более точно определился угол хр, нежели А; при этом районами наблю­ дений в целом могла оказаться только область расположения точек апекса п район, расположенный на южной оконечности зо­ ны видимости ГІСЗ.

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

1. Определим по картограмме значения ш21Г и А р для каждой подспутниковой точки и вычислим веса плоскостей синхронизации по формуле

 

 

 

 

Рі = - 4 - .

(2-113.)

 

 

 

 

mwj

 

2.

Вычислим вспомогательные величины

 

 

 

2Q = [р];

2q=[p]2A

[р cos 2А ;

 

 

[р sin 2/1

 

 

 

 

 

 

 

2q = )/ [р cos 2Ар]2+

[р sin 2А р]2 .

3.

Определим элементы

эллипса

ошибок

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

 

аг _

С

Ь2 =

С

 

[р sin 2Ар\

 

Q — q

Q + q

tg 2фо “

[р cos 2.4 J

4. Вычисляем средние квадратические ошибки направляющих углов хорды

/Пф = а2 cos2 ф0 + Ьгsin2 ф0

(2.114)

ml cos2 ф = a2 sin2 ф0 -f- b2 cos2 %

Данные конкретного примера сведены в табл. 5.

108

Т а б л и ц а 5

 

10

А*

ЗЛ°

cos 2А р

sin 2 .4 р

р cos 2.-1 р

р sіп 2 .■4

шМ\”

 

Р

п

 

 

 

 

1

3,5

40°

80°

+ 0 ,1 7 4

+ 0 ,9 8 5

4 0,61

4-3,45

2

3,7

60

120

— 0,500

+ 0 ,8 6 6

— 1,85

4-3,20

3

3,5

90

180

— 1,000

0

— 3,50

0

4

3,1

ПО

220

— 0,765

— 0,644

— 2,37

— 2,00

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

таты:

 

 

 

 

 

 

 

 

 

 

2Q =

[p] = 13,8;

Q = 6,9;

[pcos27l/;] =

— 7,11;

[p sin 2AP] ~

4,65;

2^ =

]/50,3+21,6 = 8,5;

g -■ 4,25;

lg2cp0 =

— 0,654;

2cp =

147 '■>

cp0 =

73,5°;

a2 = 3,78;

ö2 = 0,90;

m2

=

1,23;

rn\ cos2 ф = 3,55.

Вторая

задача — проектирование

необходимых

наблюдений

при

заданной

точности

окончательного

результата

т.0 = піф —

= т Л cos,!, — решается исходя из условий

 

 

 

 

 

 

 

 

= 2q = 0 и

о

2с

или

 

 

 

 

 

ms —

----

 

 

 

 

/ [/? cos 2ЛР]

\ = 0

 

 

[д]

 

 

 

 

 

 

[р] =

 

(2.115)

 

 

V[р sin 2Ар]

И

- .

 

 

 

]

 

 

т5

 

 

 

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

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

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

Ар= 104-30° и Ар=

= 1004-120° (практически размеры района

можно увеличить на

10—15° без заметного искажения равноточности углов ф и А). Если в конкретных условиях поставленной задачи веса наблю­

денных позиций ИСЗ подсчитать по формуле

р = ----, то легко

убедиться, что для соблюдения условий (2.115),

m-w

при пг0= I, нужно

будет выполнить по 9—10 наблюдений спутника в .северной и южной частях зоны.

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

109

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