книги из ГПНТБ / Разумов, О. С. Пространственная геодезическая векторная сеть
.pdfгеоцентрического сечения хорды, а их максимальный угол раство ра у>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 = s0—si :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- H» |
|
(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 ф , |
/ |
h» |
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