2014-02-09
4409
Меридиан земного эллипсоида представляет собой эллипс, радиус кривизны которого определяется величиной М, зависящей от широты. Длина дуги любой кривой переменного радиуса может быть вычислена по известной формуле дифференциальной геометрии, которая применительно к меридиану имеет выражение
(4. 25)
Здесь В1 и В2 широты, для которых определяется длина меридиана. Интеграл не берется в замкнутом виде в элементарных функциях. Для его вычисления возможны лишь приближенные методы интегрирования. При выборе метода приближенного интегрирования обратим внимание на то, что значение эксцентриситета меридианного эллипса величина малая, поэтому здесь возможно применить метод, основанный на разложении в ряд по степеням малой величины (e /2 cos2B < 7*10-3) биномиального выражения, стоящего под знаком интеграла. Число членов разложения будет зависеть от необходимой точности вычисления длины дуги меридиана, а также от разности широт ее конечных точек.
В геодезической практике могут возникать различные случаи, чаще приходится производить вычисления для малых длин (до 60 км), но для специальных целей может возникнуть потребность вычислений дуг меридианов большой длины: от экватора до текущей точки (до 10 000 км), между полюсами (до 20 000 км). Необходимая точность вычислений может достигать величины в 0. 001 м. Поэтому мы рассмотрим вначале общий случай, когда разность широт может достигать 1800, а длина дуги 20 000 км.
|
|
|
Для разложения в ряд биномиального выражения применяем известную из математики формулу.
(4. 26)
Погрешность вычисления с удержанием m членов разложения здесь достаточно определить с помощью остаточного члена в форме Лагранжа, который не меньше по абсолютной величине суммы всех отброшенных членов разложения и вычисляется по формуле
, (4. 27)
как первый из отброшенных членов разложения, вычисленный при максимально возможном значении величины x.
В нашем случае имеем
Подставляя полученное выражение в уравнение (4. 25), получим
, (4. 28)
которое допускает почленное интегрирование с удержанием необходимого числа членов разложений. Предположим, что длина дуги меридиана может достигать величины 10 000 км (от экватора до полюса), что соответствует разности широт DВ = p / 2, при этом требуется ее вычислить с точностью до 0. 001 м, что будет соответствовать относительной величине 10 –10. Значение cosB в любом случае не превзойдет единицы. Если при вычислениях будем удерживать третьи степени разложения, то остаточный член в форме Лагранжа имеет выражение
Как видим, для достижения необходимой точности такого числа членов разложения недостаточно, необходимо удерживать четыре члена разложения и остаточный член в форме Лагранжа будет иметь выражение
|
|
|
Следовательно, при интегрировании необходимо удерживать в данном случае четыре степени разложения.
Почленное интегрирование (4. 28) не вызывает труда, если преобразовать четные степени в кратные дуги (cos2nB в Cos(2nB)), используя известную формулу косинуса двойного аргумента
; cos2 B = (1 + cos2B)/2,
последовательно применяя которую, получаем
;
Действуя таким образом до cos8B, получим после несложных преобразований и интегрирования
(4. 29)
Здесь разность широт берется в радианной мере и приняты следующие обозначения коэффициентов, имеющих постоянные значения для эллипсоида с данными параметрами.
;
;
;
.
Полезно запомнить, что длина дуги меридиана с разностью широт в один градус примерно равна 111 км, в одну минуту – 1. 8 км, в одну секунду – 0. 031 км.
В геодезической практике очень часто возникает необходимость вычисления дуги меридиана малой длины (порядка длины стороны треугольника триангуляции), в условиях Беларуси это значение не превзойдет величины в 30 км. В этом случае нет необходимости применять громоздкую формулу (4. 29), а можно получить более простую, но обеспечивающую такую же точность вычислений (до 0. 001 м).
Пусть широты конечных точек на меридиане будут B1 и B2 соответственно. Для расстояний до 30 км это будет соответствовать разности широт в радианной мере, не более 0. 27. Вычисляя среднюю широту Bm дуги меридиана по формуле Bm = (B1 + B2) / 2, принимаем дугу меридиана за дугу окружности радиусом
(4. 30)
и ее длину вычисляем по формуле длины дуги окружности
, (4. 31)
где разность широт берется в радианной мере.
