По формулам со средними аргументами
При решении обратной геодезической задачи для расстояний между точками Q1 и Q2 не более 400 км наиболее удобной оказывается формула со средними аргументами.
Точности определения длины линии, прямого и обратного азимутов характеризуются следующими предельными погрешностями:
| S, км | δS, м | δА” |
| 0,01 0,1 1,0 | 0,02 0,1 0,5 |
Технологическая цепочка решения задачи и рабочие формулы:
1) по координатам точек Q1 и Q2 вычисляются значения
a Vm =√ (1 + e′ 2 cos 2 Bm), Nm = c / Vm , Mm = c / Vm 3 , c = ————
Контроль: Nm / Mm = V 2 .
2) вычисляются значения Q , Р , ∆A и Аm
ℓ 2 (2 + sin 2 Bm)
3∆B 2 + 2 ℓ 2 — 2(ℓ sin Bm) 2
3) вычисляются искомые значения прямого A12 и обратного A21 азимутов и расстояния S:
 |
S = P / sin Am = Q / cos Am = √(Q 2 + P 2 ).
Контроль – сравнение с данными, полученными при решении
прямой геодезической задачи.
Исходные данные – геодезические координаты точек 1 и 2 принять из решения прямой задачи по своему варианту.
Пример решения обратной геодезической задачи на
Дано: значения координат точек Q1 и Q2
B1 = 50 0 07′ 40.97”, B2 = 52 0 39′ 03,91”, (е′) 2 = 0,006 738525,
L1 = 23 45 13,43 , L2 = 24 00 25,46, с = 6 399 698,9018 м
Определить: значения прямого А12 и обратного А21 азимутов и
расстояния S между точками Q1 и Q2.
Решение выполняется в форме таблицы
| Обозначение | Значение | Обозначение | Значение |
| B2 рад | 0,918 934 806 | Q, м | 280 706,7048 |
| B1 рад | 0,874 899 472 | P, м | 17 636,3685 |
| ∆B рад | 0,044 035 334 | tg Am | 0,062 828 461 |
| ℓ рад | 0,004 421 597 | Am | 3 0 35′ 42,29” |
| Вm рад | 0,896 917 139 | ℓ sin Bm | 0,003 455 104 |
| sin Вm | 0,781 406 853 | ∆A“ | 712,79 |
| cos Вm | 0,624 021 898 | ∆A 0 /2 | 0 0 05′ 56,42” |
| (е′) 2 сos 2 Вm | 0,002 624 004 | A 0 12 | 3 29 45,85 |
| Vm | 1.001 311 142 | A 0 21 | 183 41 38,71 |
| Nm , м | 6 391 318,9751 | S=√ Q 2 + P 2 | 281 260,19 м |
| Mm , м | 6 374 592,0259 | S=P/sinAm=Q/cosAm | 281 260,19 м |
Вывод:результаты решения обратной геодезической задачи соответствуют результатам решения прямой задачи.
Рекомендуемая литература:
1. Бойко Е.Г. Сфероидическая геодезия, 2003 г. (глава 3).
2. Морозов В.П. Курс сфероидической геодезии, 1979 г. (глава IY)
3. Хаимов З.С. Основы высшей геодезии, 1984 г. (глава 5).
4. Практикум по высшей геодезии, 1982 г. (глава 17).
Источник
В ряды со средними аргументами
Разложение разностей широт, долгот и азимутов
Как отмечено, формулы для решения главной геодезической задачи, основанные на рядах с начальными аргументами, имеют ограниченные возможности. Рассмотрим один из возможных путей их усовершенствования, основанный на применении рядов со средними аргументами. При этом мы также будем удерживать в формулах для вычислений малые величины третьего порядка.
Пусть мы имеем некоторую геодезическую линию s12 ( рис. 6. 2 ), в середине которой имеем точку 0 с координатами B0 и L0, геодезические расстояния от которой до точек 1 и 2 будут соответственно –s / 2 и s / 2. Действуя аналогично тому, как это делали ранее, запишем следующие разложения для широты:
( 6. 9 )
Вычитая из первого уравнения второе, получим
, ( 6. 10 )
которое выражает разность широт как и первое уравнение ( 6. 4 ).
Аналогично получаем для разностей долгот и азимутов
( 6. 11 )
Отличие выражений ( 6. 10 ) – ( 6. 11 ) от ( 6. 4 ) состоит в том, что они содержат меньшее число слагаемых и проще для вычислений, но значения производных здесь следует вычислять по широте и азимуту середины линии. Возникает вопрос, быть может эти значения равны средним значениям, получаемым по формулам
Проверим, так ли это, для чего возьмем полу суммы уравнений ( 6. 9 ) и аналогичных им уравнений для долгот и азимутов, в результате получим
( 6. 12 )
Отсюда видно, что средние координаты отличаются от координат середины геодезической линии на малые величины второго порядка. При этом на данной геодезической линии существуют четыре различные точки со средними координатами. Это обстоятельство следует учитывать при дальнейшем выводе рабочих формул.
В формулах ( 6. 10 ) – ( 6. 11 ) коэффициенты разложений представлены производными, являющимися функциями широты и азимута середины геодезической линии, значения которых неизвестны. Мы можем вычислить средние широту и азимут, если известны их значения в двух точках. Поэтому перейдем в коэффициентах указанных формул к средним широтам и азимутам ( средним аргументам ). При этом будем иметь в виду порядок малых величин ( 6. 12 ). Для широт имеем
( 6. 13 )
( 6. 14 )
Учитывая значения разностей из ( 6. 12 ), выражения для производных, полученные ранее ( 6. 5 ) – ( 6. 8 ), а также вычисляя частные производные
, ( 6. 15 )
получаем для ( 6. 10 ) выражение в виде
( 6. 16 )
Действуя аналогично, получаем для разностей долгот и азимутов
( 6. 17 )
Заметим, что в правых частях полученных уравнений ( 6. 16 ) – ( 6. 17 ) дробные выражения, стоящие в скобках, являются малыми величинами третьего порядка, их можно для удобства вычислений с принятой точностью принять равными главным членам разложений и тогда можем записать
( 6. 18 )
Здесь вычисления ведутся в радианной мере, а в последнем уравнении, в поправочном члене принято
6. 5. Порядок решения прямой геодезической задачи
по формулам со средними аргументами
При решении прямой геодезической задачи известны следующие величины: B1, L1, S12 , A12 , требуется найти: B2, L2, A21. При этом очевидны уравнения:
При вычислении обратного азимута знак плюс берется, когда прямой азимут меньше π и знак минус – когда он больше π.
Особенностями применения формул ( 6. 18 ) при вычислениях искомых разностей является то, что нам неизвестны средние значения широт и азимутов, поэтому задача решается методом последовательных приближений. На первом этапе принимают ( Am )(1) = A12 ; ( Bm )(1) = B1, а значения разностей, стоящих в скобках уравнений ( 6. 18 ), равными нулю. Получают первоначальные значения разностей из этих уравнений b (1 ) , l (1) , a (1), , с учетом которых получают значения средних широт и долгот по формулам:
Для расстояний до 30 км второе приближение дает искомые разности с достаточной точностью, для контроля выполняют третье приближение. Получив значения разностей, искомые величины находят по формулам ( 6. 20 )
6. 6. Порядок решения обратной геодезической задачи
Здесь известными величинами являются: B1, B2 , L1, L2 , требуется определить: s12 , A12 , A21. Среднюю широту определяем по формуле Bm = ( B1 + B2 ) / 2 и разности широт и долгот b = B2 – B1; l = L2 – L1.
Далее замечаем, чт в ( 6. 18 ) первые два уравнения можно записать с учетом ( 6. 19 ) в виде
. ( 6. 21 )
Разделив второе уравнение на первое и ограничиваясь принятой точностью, получим после несложных преобразований уравнение для вычисления среднего азимута
Из любого из уравнений ( 6. 21 ) выражаем расстояние s и вычисляем его с контролем по формулам
Для вычисления разности азимутов используем третье уравнение ( 6. 18 ), затем вычисляем искомые прямой и обратный азимуты по формулам
6. 7. Способ Бесселя для решения главной геодезической задачи
Рассмотренный способ решения главной геодезической задачи пригоден для малых расстояний, мы привели формулы, обеспечивающие необходимую точность при расстояниях до 30 км.
Немецкий геодезист и астроном Ф. Бессель обратил внимание на то, что теорема Клеро для геодезической линии земного эллипсоида является аналогом теоремы синусов сферической тригонометрии для полярного сферического треугольника. Если в системе дифференциальных уравнений для геодезической линии земного эллипсоида ( 4. 40 )
( 6. 22 )
положить эксцентриситет, равный нулю, получим на сфере единичного радиуса
. ( 6. 23 )
Здесь α, φ, λ, σ – сферические: азимут, широта, долгота и расстояние.
Ф. Бессель рассмотрел полярные треугольники на сфере и эллипсоиде при условии, чтобы аргументы теоремы Клеро (азимуты и широты ) были бы равны на сфере и эллипсоиде. Тогда, разделив соответствующие уравнения ( 6. 22 ) и ( 6. 23 ) получим, при условии
α = A; φ = u ( 6. 24 )
следующие дифференциальные уравнения
( 6. 25 )
Интегрируя эти уравнения, получим для расстояний и долгот на эллипсоиде интегральные выражения:
.( 6. 26 )
Полученные уравнения вместе с ( 6. 24 ) определяют зависимости всех элементов соответствующего сфероидического и сферического полярных треугольников по Бесселю. Это позволяет применить формулы сферической тригонометрии для установления связи между величинами, входящими под знаки интегралов ( 6. 26 ).
На рис. 6. 3 имеем геодезическую линию, проходящую через некоторые две точки 1 и 2, пересекающую экватор в точке 0. Геодезические азимуты соответственно равны в этих точках: А1, 2π – А2 и А0 и широты: u1, u2 , u0 = 0. По условию Бесселя эти сфероидические величины сохраняются на сфере единичного радиуса, где длины больших кругов выражены в радианной мере и их значения показаны на рисунке. Точка Р – полюс, а дуги РО, Р1, Р2 – меридианы соответствующих точек.
Для дальнейших выводов введем условные долготы λ0 и сферические расстояния σ0, отсчитанные от экваториальной точки геодезической линии до текущей. При вычислении определенных интегралов ( 6. 26 ) участвуют их разности, поэтому искомые величины – расстояние и разность долгот получаются как разности
 |
Применяя аналогии Непера к прямоугольному сферическому треугольнику О110 , запишем
( 6. 27 )
Дифференцируя первое из уравнений ( 6. 27 ) по переменным величинам ( азимут А0 – величина постоянная для данной геодезической линии ), получаем
( 6. 28 )
Тогда, с учетом второго и третьего уравнений ( 6. 27 ) можем записать
( 6. 29 )
Разложим подынтегральные выражения ( 6. 26 ) в биномиальный ряд по степеням малых величин e 2 cos 2 u ( e /2 sin 2 u ) ≤ 7*10 -3
( 6. 30 )
При этом заметим порядок величин членов разложений :
≤3,5*10 -3 ;
; 
; 
Следовательно, для обеспечения необходимой точности вычислений достаточно удерживать три члена разложения. Можно заметить то, что здесь степени разложения характеризуют малые величины определенного порядка, как и в способе со средними аргументами, но преимущество данных формул в том, что величины членов разложения ( 6. 30 ) практически не зависят от расстояния, а определяются степенями малой величины эксцентриситета меридианного эллипса.
С учетом изложенного запишем интегральные выражения ( 6. 26 ) в виде
( 6. 31 )
( 6. 32 )
С учетом третьего уравнения ( 6. 27 ) можем записать для ( 6. 31 ), введя обозначение для постоянной величины k = e / cosA0
( 6. 33 )
Интеграл долготы ( 6. 32 ) с учетом уравнений ( 6. 29 ) запишем в виде
,( 6. 34 )
где принято ω = λ2 – λ1 = λ02 – λ01 – разность сферических долгот.
Производя почленное интегрирование ( 6. 33 ) и ( 6. 34 ), получаем после тождественных преобразований с принятой точностью:
, ( 6. 35 )
где приняты обозначения постоянных коэффициентов:
6. 36 )
и для разности геодезических долгот
, ( 6. 37 )
где приняты обозначения постоянных коэффициентов:
( 6. 38 )
Таким образом получены основные формулы для решения главной геодезической задачи на большие расстояния. При вычислениях на ЭВМ можно заметить, что их объем здесь несущественно отличается от вычислений по формулам со средними аргументами. Поэтому в современных условиях лучше применять способ Бесселя для решения задачи на любые расстояния.
Источник
