WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


 

На правах рукописи

  Медведев Павел Александрович 

Теория и методология повышения эффективности

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

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

25.00.32 – «Геодезия»

Автореферат

диссертации на соискание ученой степени

доктора технических наук

Омск – 2011

Работа выполнена в Федеральном государственном образовательном учреждении высшего профессионального образования «Омский государственный аграрный университет» (ФГОУ ВПО ОмГАУ)

Официальные оппоненты:

доктор физико-математических наук, профессор

Бордовицына Татьяна Валентиновна

доктор технических наук, профессор Маркузе Юрий Исидорович

доктор технических наук Мазуров Борис Тимофеевич

Ведущая организация:

Центральный научно-исследовательский институт геодезии, аэросъемки и картографии (ЦНИИГАиК), г. Москва

Защита  диссертации состоится  22 сентября  2011 г.  в  14-00  часов  на  заседании диссертационного совета Д 212.251.02 при Сибирской государственной  геодезической  академии (СГГА) по  адресу:  630108,  г. Новосибирск,  ул.  Плахотного, 10,  ауд.  403.

С диссертацией можно ознакомиться в библиотеке СГГА.

Автореферат разослан  «  »  _________________20 г.

Ученый секретарь

диссертационного совета                                Середович В.А.

Общая характеристика работы

Актуальность темы. В настоящее время для решения основной задачи геодезии по изучению фигуры Земли и ее внешнего гравитационного поля применяются астрономо-геодезический, гравиметрический и космический методы, дополняющие и контролирующие друг друга.

Так как за математическую модель Земли принимается эллипсоид вращения, то его параметры и координатная поверхность в явном или неявном видах используется при математической обработке всех видов измерений.

         Классическая редукционная задача по приведению выполненных измерений к поверхности эллипсоида решается: при построении сетей с помощью хорд эллипсоида; при трансформировании базисных линий, полученных системой GPS, в референцную систему координат; при обработке спутниковых измерений с целью вычисления наклонных расстояний между пунктами базиса; при редуцировании наклонной дальности, измеренной от спутника до объекта, для решения комбинированных засечек.

Проблема определения взаимного положения двух точек на эллипсоиде вращения, на поверхности Земли и в околоземном пространстве, получившая название «решения главных геодезических задач» (ГГЗ), под влиянием научно-технического прогресса меняет только аспекты своего решения, оставаясь актуальной длительный период времени. Формулы решения прямой геодезической задачи (ПГЗ) и обратной геодезической задачи (ОГЗ) между точками на физической поверхности Земли и в околоземном пространстве применяются: при обработке пространственных геодезических сетей;  при решениях  разнообразных геодезических  засечек; для определения уклонения отвесной линии по GPS-измерениям.

В связи с освоением шельфа и богатств Мирового океана широкое распространение при геодезическом обеспечении работ получили наземные навигационные и радионавигационные системы (типа Лоран, Омега, Селедис). С их применением координаты объектов определяются засечками (азимутальной, линейной, гиперболической), использующими формулы решения ПГЗ и ОГЗ на поверхности эллипсоида. Точность этих результатов существенно повышается при совместной обработке наземных радиогеодезических и спутниковых измерений. Решения ГГЗ на поверхности эллипсоида выполняются: для уточнения фундаментальных параметров земного эллипсоида; при установлении единой координатной системы; при ориентировке референц-эллипсоида; при уравнивании геодезических сетей сгущения; при исследованиях горизонтальных движений земной коры; при определениях уклонения отвесной линии на морской поверхности по альтиметрическим измерениям; при подготовке высокоточных маршрутов движения морских судов и воздушных объектов, при запусках ракет и ИСЗ, в целях  слежения за управляемыми ракетами; для определения промежуточных точек геодезической линии.

Теория определения взаимного положения точек на поверхности эллипсоида разработана Эйлером Л. еще в 1753 г.  Выведенные им дифференциальные уравнения геодезической линии (ДУГЛ) составили теоретическую основу построения математических моделей решения ГГЗ. Так как интегралы этих уравнений не выражаются в конечном виде через элементарные функции, то для приближенных решений применяют различные методы аппроксимации первообразных.

В течение более двух столетий, начиная с Эйлера Л., Лежандра А., Ориани Б., как отечественными, так и зарубежными учеными разработано значительное количество способов решения ГГЗ. Методика построения и практическая реализация математических моделей решения ГГЗ всегда определялись уровнем развития вычислительной техники, а их точность обусловливалась дорогостоящими технологиями производства геодезических работ. Для исключения вычислительных погрешностей к выводимым формулам предъявляется требование, чтобы точность определяемой по ним величины была на один-два порядка выше точности, вызванной погрешностями измерений.

Общий вывод о том, что для каждого конкретного случая следует выбирать соответствующие формулы, был, безусловно, верным при ручных способах счета. При использовании современной вычислительной техники в арсенале практических приложений должно остаться наименьшее число наиболее универсальных методов. Было бы идеальным, при любых расстояниях и с любой точностью решать задачи одним способом. Но, ввиду большого спектра приложений решения ГГЗ и значительного отличия решений по эффективности, используются разные методы построения математических моделей при больших и малых расстояниях.

При малых расстояниях приближенные методы решения ДУГЛ представляются, в основном, в форме отрезка ряда Тейлора. Внедрение в геодезическое производство радиоэлектронных средств в комплексе с ЭВМ способствовало разработке теоретических вопросов решения ГГЗ с привлечением нового математического аппарата, исходя из проблемы больших расстояний. Для решения ДУГЛ, кроме разложения в ряд Фурье, стали применять численные методы, разные аппроксимации, построенные с помощью рациональных дробей, полиномов Чебышева, методом экономизации.

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

Кроме этого, погрешности методов решения задач устанавливаются при условии, что исходные данные являются точными. При вычислениях с приближенными числами образуется неустранимая погрешность, обусловленная неточностью исходных данных, которая существенно влияет на точность результата при большом объеме вычислительных операций. Поэтому решения, приводящие к выполнению большого числа действий, образуют новую, очень важную проблему обоснования достоверности полученного результата. Эти положения не учитывались: при решениях ДУГЛ методом Рунге-Кутта и его модификациями, которые получили широкое распространение в практике геодезических вычислений; в методе решения ОГЗ при предельных расстояниях способом определения промежуточной точки.

Развитие спутниковых технологий позволило реализовать миллиметровый уровень точности при измерениях до тысячи и более километров, а  Федеральным законом «О геодезии и картографии» планируется создание геодезической сети на качественно новом, более высоком уровне точности, обеспечивающем решение фундаментальных и перспективных задач в области геодезии, геофизики, геодинамики, космонавтики, экономики и обороноспособности РФ.

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

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

При разработке математических моделей решения ГГЗ на эллипсоиде используются элементы вспомогательного сферического треугольника. С целью обоснования общего подхода к решению задач Морозов В.П. вывел дифференциальные уравнения, связывающие элементы сферического и сфероидического треугольников. Но в процессе исследований им не получена общая математическая модель, объединяющая всевозможные решения этих уравнений. В связи с этим не выполнен сравнительный анализ и не выведены оптимальные математические модели решения ГГЗ.

Разработанные в последние десятилетия простые и компактные формулы решения ГГЗ или не соответствуют по точности миллиметровым погрешностям измерений, или рассчитаны на применение при вычислениях логарифмов, разных таблиц и других вспомогательных средств. Необходимость получения более простых моделей решения ГГЗ с высокой точностью обусловлена как востребованностью в производственных вычислительных процессах, так и их методическим значением в процессе обучения студентов и аспирантов по геодезическим специальностям.

Кроме этого, существующие обзоры по решению ГГЗ не носят объективного и целостного подхода в связи с исключением из анализа публикаций Лежандра и Баховена, предложивших способы решения задач, известные в литературе под именами Бесселя и Мак-Коу. Не рассматриваются исследования Ориани, составившего полную сфероидическую тригонометрию и разработавшего методы решения ГГЗ.

Объект исследования – теория и методы определения взаимного положения точек на поверхности эллипсоида, физической поверхности Земли и в околоземном пространстве.

Предмет исследования – теоретические обоснование и систематизация разных принципов и подходов к решению ГГЗ; разработка методики сравнительного анализа математических моделей решения ДУГЛ по кратчайшей линии; развитие теории и разработка методологии совершенствования известных, построения новых математических моделей эффективного решения ГГЗ  на поверхности эллипсоида и в пространстве с использованием современной вычислительной техники.

Цель исследования –  совершенствование и разработка оптимальных математических моделей решения ГГЗ на поверхности эллипсоида и в пространстве, обеспечивающих точность результатов на один-два порядка выше точности, обусловленной миллиметровыми погрешностями измерений.

Для достижения поставленной цели необходимо решить следу-ющие задачи

        1. Теоретически обосновать и усовершенствовать методологию вывода формул решения геодезических задач на малые расстояния путем преобразований разложений Лежандра без использования поверхности вспомогательной сферы. Разработать теорию и метод построения математических моделей и вывода на их основе формул с улучшенной сходимостью для эффективного решения ГГЗ. Исследовать точность формул со средними аргументами, используя принцип остаточных членов.
        2. Развить теорию и разработать метод исследований ДУГЛ, полученных с помощью теоремы Клеро. Выполнить теоретические исследования, применяемых в сфероидической геодезии биномиальных функций, и построить математическую модель обобщенного эллиптического интеграла для вычисления длины геодезической линии. На основе разложения обобщенного эллиптического интеграла в ряд с общим членом в замкнутой форме, провести классификацию формул по точности, структуре и способам их вывода. Установить перспективные направления исследований, приводящие к повышению эффективности и точности решения ГГЗ. Создать обобщенную математическую модель для вычисления геодезической долготы и выполнить ее анализ.
        3. Разработать теоретические основы и методы совершенствования известных и вывода новых математических моделей решения ГГЗ с высокой точностью на поверхностях эллипсоида и шара на большие расстояния, включая предельные. Выведенные математические модели в виде системы формул, приводящие к эффективному решению задач, удовлетворяют следующим требованиям: имеют логически обоснованную структуру с простой оценкой их достоверности и ясной геометрической интерпретацией; содержат компактные выражения, приведенные к удобному для вычислений виду, с наименьшим числом операций и промежуточных результатов; по значениям тригонометрических функций определяют аргументы без дополнительных исследований по установлению их квадрантов; обеспечивают взаимно однозначное соответствие по точности при решениях прямой и обратной геодезических задач; являются удобными для проведения дальнейших теоретических исследований и аналитического решения смежных задач.
        4. Выполнить  анализ решений ГГЗ, полученных численными методами, с помощью дробно-рациональных приближений, применением функций Якоби. Теоретически исследовать процесс экономизации степенных рядов на основе тождества, выражающего косинусы кратных аргументов через степени косинуса аргумента. Развить теорию и разработать обобщенный метод экономизации, используемых в сфероидической геодезии разложений эллиптических радикалов и интегралов, без применения полиномов Чебышева. Преобразовать формулы решения ГГЗ.
        5. Усовершенствовать математические модели решения ГГЗ в пространстве. Для этого разработать теорию и метод отделения корня трансцендентного уравнения для вычисления широты при переходе от пространственных прямоугольных координат к геодезическим .Используя принцип отделения корня, вывести разными способами формулы вычисления широты и закономерности для оценки ее погрешностей по итерациям. Обосновать метод дальнейшего сужения промежутка изоляции корня и вывести неитеративную формулу для определения широты с высокой точностью.

Методы исследований. Теоретическую основу исследований составили работы Лежандра, Ориани, Баховена и Беспалова с использованием разработок по сфероидической геодезии других авторов. Решение поставленных задач базировалось на методах математического моделирования, анализа и синтеза, сравнения, аналогий, обобщений и оценок с привлечением теории рядов, численных методов, математического анализа, аппарата цепных дробей, функций Якоби, полиномов Чебышева, аналитической и дифференциальной геометрии, высшей алгебры, теории погрешностей, дифференциальных уравнений и комбинаторики.

На защиту выносятся

  1. Теория и методология вывода обобщенных разложений Лежандра без использования поверхности вспомогательной сферы.  Теория и метод построения оптимальных моделей и вывода на их основе формул эффективного решения ГГЗ на малые и средние расстояния путем выделения из обобщенных разложений Лежандра сферических величин по закономерностям сферической тригонометрии. Усовершенствованный метод исследования точности формул со средними аргументами.
  2. Теория и результаты исследований ДУГЛ, полученных с помощью теоремы Клеро, на основе анализа и геометрической интерпретации подстановок, преобразующих несобственные интегралы в определенные. Совершенствование методов вывода формул решения ГГЗ с применением введенного понятия криволинейного вектора. Концепция и методология построения математической модели для вычисления долготы. Математическая модель разложения обобщенного эллиптического интеграла, способствующая изучению свойств и структуры формул для вычисления длин геодезических линий, выведенных по способам Лежандра и Баховена.
  3. Теоретические исследования и выведенные математические модели решения ГГЗ с высокой точностью на поверхностях эллипсоида и шара на большие расстояния, включая предельные, удовлетворяющие выдвинутым в задаче 3 требованиям.
  4. Теория и методология анализа и вывода формул решения геодезических задач в дробно-рациональном виде способом геометрической прогрессии. Теоретическое обоснование и обобщенный метод экономизации разложений эллиптических радикалов и интегралов на основе тригонометрических тождеств без использования полиномов Чебышева, его приложения к решению ГГЗ.
  5. Усовершенствованные математические модели решения геодезических задач в пространстве. Теория и метод отделения корня трансцендентного уравнения для вычисления широты при переходе от пространственных прямоугольных координат к геодезическим . Концепция и методы построения рабочих формул для определения широты и методология их сравнительного анализа. Теоретическое обоснование метода сужения промежутка изоляции корня и вывода неитеративной формулы высокой точности для определения широты с высокой точностью.

Научная новизна исследований

  1. Теоретически обоснована и усовершенствована методология вывода формул решения ГГЗ на малые расстояния путем преобразований рядов Лежандра без использования поверхности вспомогательной сферы. Такой подход позволил получить обобщенные математические модели и  привести в систему решения ГГЗ, полученные разными способами, осуществить их анализ более простым путем, не требующим учета погрешностей, вызванных заменой одной поверхности другой; разработать теорию и метод вывода математических моделей для эффективного решения ГГЗ на короткие и средние расстояния с выделенными из разложений сферическими членами путем преобразований обобщенных разложений Лежандра по схемам моделей сферических закономерностей.  Принципиальное отличие предложенного метода состоит в том, что при выводе преобразуются не формулы сферической тригонометрии с помощью зависимостей изображения эллипсоида на шаре, а обобщенные разложения Лежандра по схемам моделей сферических закономерностей. По предложенной методике выведены: оптимальные формулы решения ПГЗ; формулы  решения ОГЗ, в которых все члены имеют унифицированную и удобную для вычислений структуру.

Универсальность предложенного метода характеризуется следующими признаками: его применение не требует никакой проекции эллипсоида на шар; вывод формул осуществляется по единому логически обоснованному принципу как для решения ПГЗ, так и для решения ОГЗ; пригодностью для преобразований не только обобщенных разложений Лежандра, но и закономерностей изображения эллипсоида на шаре; возможностями использования для упрощения зависимостей картографических проекций эллипсоида на плоскость с помощью замкнутых формул проекций эллипсоида на шар.

Усовершенствованный метод исследования точности формул со средними аргументами для решения ОГЗ, позволяющий находить наибольшие погрешности определяемых величин,  реализован на сравнительном анализе выведенных формул с классическими  при сохранении в разложениях членов разных порядков, доказывая преимущества первых.

2. Теоретический анализ и разработанный метод исследований ДУГЛ, представленных по переменным , , обосновал перспективность вывода формул интегрированием уравнений по сферической дуге , так как интегрирование по другим переменным не приводит к новым результатам, а ведет только к сложным методам вывода формул.

       Усовершенствованные методы исследования четырех видов ДУГЛ с переменной , с геометрической интерпретацией  подстановок Ориани и Баховена позволили: установить зависимости между элементами сферических треугольников, с помощью которых дифференциальные уравнения одного из видов не только приводятся к  другим,  но и по решениям одного сфероидического треугольника определяются решения остальных трех видов без вывода дифференциальных уравнений и интегрирования.

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

       Теоретическими разработками установлено, что все формулы, выведенные по способу Лежандра для вычисления долготы, являются следствиями усовершенствованной общей модели Ориани. Для определения долготы предложена модель с точными выражениями постоянных величин , упрощающая ее структуру и процесс вычислений.

  1. Теоретически обоснована и разработана методика вывода рабочих формул эффективного решения ГГЗ на поверхностях эллипсоида и шара на большие расстояния. На ее основе построены математические модели в виде упорядоченных систем формул для решения прямой и обратной геодезических задач в функциях тангенсов половинных аргументов, не требующих дополнительных исследований по установлению квадрантов определяемых величин и  обеспечивающих числовую стабильность вычислительных операций с наименьшим числом десятичных разрядов.

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

       Разработана математическая модель решения ОГЗ на эллипсоиде, при антиподном расположении точек на вспомогательной сфере, по которой при малом объеме вычислительных операций, результаты определяются неитеративным способом с высокой точностью. Математическая модель решения  ОГЗ при точках, близких к антиподным, по сравнению с предложенной Боурингом, имеет более простые зависимости и обеспечивает квадратичную сходимость итерационного процесса при определении азимута. Итеративные формулы решения ОГЗ для достижения высокой точности результатов требуют менее десяти итераций, затрачивая на решения одну-две секунды машинного времени.

  1. Доказана тождественность дробно-рациональных приближений бинома , выполненных с помощью подходящих дробей цепной дроби Лагранжа и предложенным методом геометрической прогрессии, который не только приводит к тем же результатам более простым способом, но и имеет более широкую область применения. Теоретически обоснована эффективность построения моделей в дробно-рациональном виде методом геометрической прогрессии путем преобразования разложений эллиптических радикалов и интегралов, приводящих к компактным формулам с меньшей погрешностью без увеличения объема вычислительных операций. Этим методом приведены к дробно-рациональному виду, полученные автором зависимости для вычисления и с параметром и  формулы Ганзена, для сравнения с выведенными Беспаловым Н.А.

       Теоретическое обоснование процесса экономизации степенных рядов с использованием полиномов Чебышева на основе тригонометрического подхода, способствовало разработке обобщенного метода экономизации, базирующегося на тригонометрических тождествах. Его достоинства: экономизация этим методом выполняется без применения полиномов Чебышева и позволяет убедительно и просто установить равносильные между собой отрезки рядов, выявить причины, приводящие к неравноточным зависимостям, указать пути преобразования их к тождественному виду; обобщенный метод экономизации разложений эллиптических радикалов позволяет привести к экономизированному виду не только ряды, представленные по степеням , но и по степеням и ; тригонометрический принцип экономизации способствовал обоснованию и применению его к формулам, полученным с помощью обобщенного эллиптического интеграла и построению оптимальных рабочих формул решения ГГЗ; пригодностью для рядов с двумя переменными.

  1. Разработана теория и метод отделения корня трансцендентного уравнения для определения широты при переходе от прямоугольных координат к геодезическим и по предложенной методике установлен промежуток , внутри которого он содержится для точек с высотами . На основе свойств малой длины отрезка корня обоснована методология и получены рекуррентные зависимости, позволившие провести сравнительный анализ итерационных процессов вычисления геодезической широты и высоты как по выведенным, так и по рекомендованным для практического использования формулам. Выполнена классификация итерационных процессов определения широты по видам и скорости их сходимости, на основе которой уменьшен интервал корня и выведена неитеративная формула для вычисления геодезической широты с высокой точностью при .

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

Научная и практическая ценность работы

       Теоретические положения и разработанные методы эффективного ре­шения ГГЗ на малые и большие расстояния позволили решение этой про­блемы выполнить на высоком научном уровне и привести в целостную и стройную систему. Полученные высокоточные формулы, иллюстрируе­мые в Приложениях диссертации числовыми примерами, целесообразно использовать при решениях фундаментальных перспективных задач в области геодезии, геофи­зики, геодинамики и космонавтики, обороноспособности РФ, а также при геодезическом обеспечении съемок акваторий морей и океанов, для удовле­творения возрастающих потребностей ракетной техники, морской, воздуш­ной и космической навигации.

       Достоверность полученных в работе научных результатов, выводов и математических моделей решения ГГЗ основывается на  использовании при исследованиях фундаментальных положений сфероидической геодезии и математических наук; обеспечивается строгими математическими преобразованиями при решениях уравнений и выводе формул; теоретическими обоснованиями и доказательствами сходимости итерационных процессов определяемых величин; анализом математических моделей по точности с помощью выведенных оценок и путем сравнения результатов взаимно обратных решений ГГЗ.

Реализация результатов работы. Итоги научной работы методологического характера были внедрены в учебный процесс кафедры высшей геодезии, фотограмметрии и географических информационных систем ФГОУ ВПО «Омский государственный аграрный университет» в лекционных и практических курсах по дисциплинам: «Математическое моделирование на ЭВМ», «Высшая геодезия», «Основы космической геодезии» по специальности 300100 – «Прикладная геодезия».

Усовершенствованные математические модели решения ГГЗ в пространстве внедрены в производственную практику Западно-Сибирского филиала ФГУП «Госземкадастрсъемка» – ВИСХАГИ для разработки методики создания планового геодезического обоснования с применением спутниковой системы GPS при межевании земель и использованы при камеральной обработке в процессе создания спутниковых геодезических сетей на территориях 13 областей и республик РФ.

Полученные в диссертации результаты составили теоретическую основу научных исследований по теме  «Решение главных геодезических задач на поверхности сфероида и в пространстве», зарегистрированной Всероссийским научно-техническим центром: 12.03.2003 № 01.9.90 000 949.

Апробация работы. Основные положения диссертации докладывались на научно-технических конференциях ФГОУ ВПО «Омский государственный аграрный университет» в период с 1980 по 2009 гг. (г. Омск); научных конференциях профессорско-преподавательского состава Сибирской автомобильно-дорожной академии (СибАДИ, г. Омск); научно-технической кон-ференции преподавателей ГОУ ВПО «Сибирская государственная геодезическая академия» (1998 г.); на третьем Сибирском конгрессе по прикладной и индустриальной математике (1998 г.) г. Новосибирск; научных семинарах кафедры высшей геодезии (1995–2000 гг.) ОмГАУ; на расширенных заседаниях кафедры прикладной геодезии, фотограмметрии и географических информационных систем ФГОУ ВПО ОмГАУ  (2009, 2010 гг.); на расширенном заседании кафедры высшей геодезии ГОУ ВПО СГГА (2009 г.).

Публикации. Основные результаты диссертационных исследований опубликованы в 30 работах, из которых девять статей (одна в соавторстве) – в журнале «Геодезия и картография», пять – в журнале «Известия вузов. Геодезия и аэрофотосъемка», входящих в Перечень  изданий, определенных ВАК  Минобрнауки РФ, в одной монографии, двух статьях в республиканском межведомственном сборнике «Геодезия, картография и аэрофотосъемка» (г. Львов) и научных трудах других учебных заведений.

Структура и объем работы. Диссертация состоит из введения, пяти разделов, заключения, библиографического списка и приложений. Общий объем составляет 426 страниц машинописного текста, 12 рисунков, 16 таблиц, 6 приложений. Список литературы включает 327 наименований, в том числе 115 на иностранных языках.

       

Основное содержание работы

Во введении обоснована актуальность рассматриваемой темы, определены объект и предмет исследования, поставлены цель и задачи, сформулированы положения, выносимые на защиту, научная новизна и практическая значимость полученных результатов.

В первом разделе рассмотрена теория и  методология повышения эффективности решения ГГЗ на поверхности эллипсоида при малых расстояниях и выведены оптимальные формулы решения прямой и обратной геодезических задач.

На основе математической модели, представленной системой ДУГЛ

  , ,  (1)

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

В качестве исходной взята математическая модель решения уравнений (1), построенная в форме рядов Тейлора способом Эйлера

  ,  ,  ,  (2)

с начальным условием прохождения кривой через точку в направлении . Впервые формулы для решения ПГЗ этим методом вывел в 1806 г. Лежандр, сохранив в разложениях (2) производные по третий порядок включительно. При выводе многих формул используется поверхность вспомогательной сферы. В связи с этим, сравнительный анализ и дальнейшее совершенствование методов решения задач приводят к сложным преобразованиям по учету погрешностей, вносимых заменой одной поверхности другой. Поэтому в работе дается вывод формул, получивших наибольшее распространение, без использования поверхности сферы путем непосредственного преобразования рядов Лежандра (2). Эти разложения рассматриваются как исходное целостное построение, на структурных и функциональных свойствах которого образуются другие модели, способствующие выводу новых формул решения ГГЗ, более ценных в практическом отношении.

Формулы Гаусса с сохранением членов третьего порядка непосредственно выводятся из (2). Для этого содержащиеся в них функции широты и азимута по формуле Тейлора преобразованы к средним аргументам и . В результате получены неполные степенные ряды, содержащие только нечетные степени , с коэффициентами в раза меньшими, которые символически представляются в виде

, ,   (3)

Если ПГЗ по зависимостям (3) решают способом последовательных приближений, то формулы, полученные обращением (3),

  ,  ,  ,  (4)

широко используют для неитеративного решения ОГЗ.

Разложения (2) существенно упрощаются при азимутах и , что приводит к образованию поверхностных прямоугольных координат  . С помощью формулы Тейлора и преобразованиями  рядов (2) получены для вычисления разностей координат и азимутов две группы отрезков рядов в функциях вспомогательной точки с широтой  , представленных по степеням , :

, ,,, , (5)

и по , :

,  ,  ,  ,  .  (6)

Все группы формул (2, 5, 6) являются равноточными, погрешности которых оцениваются отброшенными величинами шестого порядка малости. В работе системы (2–6) выведены без использования поверхности вспомогательной сферы. Поскольку отрезки рядов (3–5) являются следствиями (2), то равенства (2–6) названы обобщенными разложениями Лежандра.

       По разложениям (5) и (6), сохраняя в них величины по четвертый порядок включительно, путем введения вспомогательных величин получены нелогарифмические формулы, выведенные Шрейбером О., Изотовым А.А., а также автором. Выполнен их сравнительный анализ с использованием оценок погрешностей. Доказано, что по скорости сходимости и объему вычислительных операций они не являются оптимальными. К более эффективным результатам приводит структурно-функциональный принцип решения сфероидических треугольников, предложенный Эйлером Л. Деление отрезков рядов на сферическую и сфероидическую составляющие позволило сферическую часть решения задачи выполнить по строгим формулам сферической тригонометрии. Оставшаяся сфероидическая часть имеет не только более компактный вид, но и оценивается отброшенными членами сфероидического вида, меньшими на один-два порядка соответствующих сферических. Такой подход приводит к более простым формулам  и повышает их точность.

Этот способ в 1919 г. применил Крюгер Л. к разложениям (2). Решения ПГЗ на шаре он выполнил по формулам четвертого способа Гаусса. Наиболее сложная формула для вычисления широты автором преобразована к более простому виду. При  решениях этим способом Шедльбауэр А. и Крак К. применили теоремы косинусов, синусов, котангенсов и другие формулы сферической тригонометрии. Другой путь вывода формул, требующих при решении ОГЗ применения вспомогательной сферы, предложили Иордан В., а затем и Морозов В.П. На основе геодезического изображения эллипсоида на шаре они получили зависимости со средней широтой для вычисления и в форме рядов с поправочными членами сфероидического вида.

Для более эффективного использования формул были выполнены разработки, исключающие при решении двойной переход с эллипсоида на шар и обратно. Буткевич А.В. с этой целью преобразовал правые части аналогий Не-пера к сфероидическим величинам с помощью полученных Иорданом В. зависимостей геодезического изображения эллипсоида на шаре. При таком подходе в системе предложенных логарифмических формул оказалась неисключенной сферическая величина . Формулы для вычисления геодези-ческих координат по заданным прямоугольным поверхностным координатам Гельмерт Ф. вывел путем преобразования формул четвертого способа Гаусса решения сферического треугольника к сфероидическим аргументам с по-мощью установленных им зависимостей геодезического изображения полярного сфероидического треугольника на шаре. В результате громоздких и сложных преобразований он получил отрезки рядов, поправочные члены в которых содержат только сфероидические величины. Если в равенствах, полученных Гельмертом Ф., первый эксцентриситет выразить через второй , то приходим к формулам Крюгера, выведенные им преобразованием (6).

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

Для построения таких формул необходимо сумму сферических членов ряда представить замкнутым выражением. Однако такой подход в общей постановке не приводит к цели. Решение этой задачи упрощается, если исходный ряд заменить равносильным ему с известной суммой. А это достигается только с помощью формул сферической тригонометрии. Но, если в предыдущих случаях они предназначались для решения задач на шаре, или служили основой для приведения их выражений к сфероидическим величинам с учетом условий изображения эллипсоида на шаре, то теперь они используются для построения математических моделей, по схемам которых выполняются преобразования рядов с целью выделения из них сферических величин. Для этого в формулах решения геодезических задач на шаре производится замена сферических элементов полярного треугольника на соответствующие сфероидические. В полученной таким способом модели ее левая часть указывает, какой функцией следует представить исходный ряд для определения искомой величины, а правая – сумму сферических членов преобразованного ряда.

По предложенной методике выделены сферические члены из нелогарифмических формул Шрейбера по закономерностям четвертого способа  Гаусса. Доказано, что относительная простота формул достигается при сохранении в разложениях величин не более четвертого порядка. С такой погрешностью получена зависимость для вычисления сфероидического избытка, а, выведенная система для решения ПГЗ, по объему вычислительных работ при одинаковой точности результатов, является более экономичной по сравнению с предложенной Гельмертом–Крюгером.

Доказано, что выделение сферических членов из разложений Лежандра (2) по зависимостям теоремы синусов и формулам косвенного метода решения задач на шаре, приводят к громоздким и сложным выражениям для вычисления и . Из всех формул решения ПГЗ с погрешностями на величины пятого порядка наиболее оптимальными являются, выведенные нами из (5) по закономерностям четвертого способа Гаусса:

,

,

,

,

,  (7)

В приведенных выражениях погрешности при и .

Для сравнительного анализа в работе получен нелогарифмический  вариант формул, состоящий из пяти уравнений со средними аргументами, предложенный Буткевичем А.В. для решения ОГЗ. Преобразованием формулы Иордана для определения по предложенной методике получена зависимость, с помощью которой исключается сферическая величина и система из пяти уравнений сводится к трем формулам с выделенными сферическими величинами

,

,

.  (8)

Однако равенства (8) эффективнее получаются из разложений (4). Преобразованием их выражений с учетом членов пятого порядка по аналогиям Непера автором получены формулы для функций и . Но такие построения в частных случаях приводят к неопределенностям вида . Формулы, свободные от этих недостатков, выведены путем преобразований (4) по закономерностям Деламбра. Для упрощения выражений вместо разности широт введена величина . Получены три варианта разложений в функциях , , . Для практического использования рекомендуются формулы в форме рядов по степеням :

, , , ,

,

,

  ,

  ,  , ,

  ,  ,  ,

, .  (9)

Если , то и

.

Отличительные особенности приведенных формул (9) заключаются в следующем: они не имеют аналогов по своей структуре; определяемые величины выражаются только через исходные данные; отсутствуют неопределенности  при любом расположении геодезической линии; главные и поправочные члены разложений имеют унифицированную структуру, способствующую эффективному решению ОГЗ.

Для исследования точности формул применен принцип остаточных членов. Дано обоснование, что при остаток ряда можно оценить по наибольшим величинам одного и того же порядка малости первых отброшенных членов. В этом случае величины и эксцентриситеты считаются малыми первого порядка.

Доказано, что для формул решения ОГЗ, представленных в функциях , , , дифференциальные зависимости между погрешностями определяемых величин , , и остаточными членами , , выражаются одинаковыми зависимостями:

  (10)

  (11)

  (12)

Для упрощения выводов величины с помощью главных членов разложений (4) преобразованы к переменным , , .

Исследованием на экстремум при выбранной широте устанавливается , при которой принимает наибольшее значение. После этого для погрешности с помощью (10) определяется предельная длина геодезической линии, характеризующая область применения анализируемых формул.

Аналогично, но с учетом полученного , на той же широте исследованием функций, содержащихся в (11) и (12), определяются наибольшие значения и . Погрешность прямого и обратного азимутов находится по соотношению

.

Анализ формул проводился на широтах при вычисленных длинах геодезических линий с погрешностями ; 0,001; 0,01; 0,1; 1 (м).

По этой методике выполнен сравнительный анализ классических формул (4) с полученными при сохранении в разложениях членов соответственно первого, третьего и пятого порядков. И хотя с учетом членов третьего порядка эти алгоритмы являются тождественными, однако по (9) результаты вычисляются с более высокой точностью, что увеличивает пределы их действия при и с 125 км до 440 км. Это обусловлено тем, что погрешности формул (9) без членов пятого порядка определяются величинами седьмого порядка малости, а не пятого, как в отрезках рядов Иордана-Эггерта.

Полные формулы (9) на широтах имеют следующие характеристики:

если

если

Сохранение в них первых членов на широтах и обеспечивают вычисление результатов с погрешностями , .

       Итак, принцип пересмотра известных методов решения ГГЗ с позиций построения их математических моделей путем преобразования рядов Лежандра без использования поверхности вспомогательной сферы, позволил эти решения не только обобщить и привести в систему, но и теоретически обосновать и разработать новый метод эффективного решения ПГЗ и ОГЗ.

       Второй раздел посвящен развитию теории и разработке методов построения обобщенных математических моделей решения ДУГЛ интегрированием их аппроксимаций, представленных в форме рядов.

Другой подход к построению модели взаимного положения двух точек на поверхности эллипсоида основан на свойствах общего решения третьего уравнения системы (1), полученного в 1733 г. Клеро А. в виде

    (13)

Методы решения задач во многом обусловлены начальными условиями интегрирования и видом представления равенства (13), допускающего различные аналитические преобразования и геометрические интерпретации.

ДУГЛ в функциях с геодезической широтой В, при любом ее расположении на поверхности сфероида, были выведены Эйлером Л. в 1753 г. Неопределенную постоянную интегрирования С он определил из условия прохождения кратчайшей линии через точку , т.е. при и . Но полученные им уравнения имеют сложную структуру, а результаты интегрирования, приведенные к рабочему виду, образуют громоздкие выражения при низкой точности определяемых величин.

Если для определения С использовать координаты , точки вертекса или координаты , точки пересечения линии с экватором, то уравнения Эйлера представляются в виде (14–16), полученном Баховеном Е. в 1865 г., и (17–19) с приведенной широтой, опубликованные Ориани Б. в 1806 г.:

(14)

(17)

,

(15)

(18)

(16)

(19)

Эти уравнения Ориани Б. и Баховен Е. преобразовали к более простой форме.

Преобразования Баховена Е.:

,

(20)

(23)

, ,

(21)

(24)

(22)

(25)

Преобразования Ориани Б.:

,

(26)

,

(27)

(28)

(29)

, ,

(30)

(31)

В геодезической литературе разные варианты полученных решений связывают с переменными, по которым выполнялось интегрирование. В связи с этим ДУГЛ были представлены в функциях переменных и др. Так как результаты интегрирования с помощью соответствующих подстановок снова выражаются через исходные, то для интегрирования применяют те из них, первообразные с которыми находятся наиболее простыми методами. Показано, что использование других переменных, кроме , для вывода формул не является рациональным.

С целью построения качественных, более простых и логически обоснованных решений четырех видов ДУГЛ, полученных Баховеном Е. и Ориани Б., необходимо выполнить их анализ с геометрической интерпретацией полученных результатов. Для этой цели применяется не общая теорема Клеро, представленная в форме теоремы синусов на сфере, а используются подстановки, преобразующие несобственные интегралы к определенным и приводящие к этим дифференциальным уравнениям. Подстановка Баховена (20) позволяет по катетам и построить полярный сферический треугольник с гипотенузой . В этом случае  основное уравнение геодезической линии (14) принимает вид

,  (32)

а сфероидическому полярному (рисунок А.7) ставится в соответствие сферический полярный (рисунок А.8),  решения которого используются для вычисления элементов сфероидического путем интегрирования соответствующих ему дифференциальных уравнений (21, 22) при начальных условиях, обусловленных координатами точки , принятой за начало отсчета длин.

Аналогично, по второй подстановке Баховена (13) строится прямосторонний сферический треугольник . В этом случае вычисление сфероидических элементов выполняется по дифференциальным уравнениям (24, 25) при условии начала отсчета длин линий от экватора. Такое изображение сфероидического треугольника на сфере связывают с именем Мак-Коу, который получил уравнения Баховена в 1932 г., т.е. через 67 лет.

По зависимостям  (17, 26, 29) строится изображение геодезической линии на поверхности сферы способом Лежандра, который на 36 лет раньше Бесселя разработал и применил его для решения ПГЗ при обработке рядов триангуляции, проложенных на побережьях Франции и Англии через пролив Ла-Манш. В результате такого подхода каждому из двух видов сфероидических треугольников и (с началом отсчета длин линий соответственно от точки экватора и точки вертекса ) оказались поставленными в соответствие по два сферических прямосторонних  и прямоугольных  треугольников, построенных по способам Лежандра и Баховена.

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

,  ,  (33)

связывающие элементы прямоугольного и прямостороннего треугольников Лежандра. Выражения (33) при выбранном положительном направлении преобразуют формулы сферической тригонометрии для решения прямоугольного треугольника в формулы для решения прямостороннего сферического треугольника и наоборот. Кроме того, подстановка равенство (29) переводит в (26), а дифференциальные формулы Ориани (30, 31) приводит к (27, 28).

Аналогично устанавливаются соотношения между элементами треугольников Баховена с теми же свойствами, что и в способе Лежандра:

, . (34)

Зависимости между элементами сферических треугольников Лежандра и Баховена применяли Мак-Коу, Ганьшин и другие для вычислений величин одного треугольника по известным данным другого. Но эти результаты можно использовать и для обоснования других выводов. В связи с тем, что углы при вершине полюса в треугольниках Лежандра и Баховена равны, т.е. , , это свойство позволяет построить их на одной шаровой поверхности. Поэтому понятия «сфера Лежандра» или «сфера Баховена» теряют смысл.

       Зависимости

    (35)

(36)

преобразуют уравнение (30) к виду (24), а (31) – к (25), которые в силу (34) приводятся к (21, 22). Символически эти преобразования можно представить равенствами

  ,  (37)

  . (38)

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

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

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

Прямая задача: по координатам точки , длине вектора и его направлению в точке определяются координаты и направления вектора в точке .

Обратная задача: по координатам точек и определяются длина вектора и его направления в начальной и конечной точках.

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

Ограниченное количество показателей степени , используемых при построении формул для решения геодезических задач, позволяет образовать обобщенный степенной ряд

  . (39)

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

  ,  (40)

то степень достоверности приближенного равенства при любых знаках его членов оценивается по абсолютной величине первого отбрасываемого члена ряда (39). Исследования, проводимые на основе (39), дают результаты сразу для трех случаев и позволяют определить скорость сходимости ряда (39) в зависимости от величины . По закономерности между общими членами (39) при и при отрицательных устанавливаем, что вывод формул на основе ряда (39) с показателем является более перспективным, чем с . Например, при формула для вычисления длины геодезической линии, полученная по способу Лежандра интегрированием равенства (30), будет на два порядка точнее формулы, полученной по способу Баховена на основе (21).

Переменная в основных сфероидических функциях эллиптического вида представляется в формах и , где . Если , то обобщенное разложение (39) переходит в функциональный ряд

  (41)

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

  (42)

где 

,

    .  (43)

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

  (44)

с прежними выражениями коэффициентов. Этот ряд можно представить в других видах с помощью тождеств

  (45)

(46)

  ,  (47)

где – целая часть числа .

       Разложение (44) с учетом (45) преобразуется к виду

  (48)

На основе зависимости (48) образуется рекуррентная формула

    (49)

Доказано, что в результате применения подстановки в исходном разложении происходит простая замена переменных по соотношениям:

,

,

,

,

,

  (50)

при тех же выражениях коэффициентов.

В результате, как следствие (48), без выводов получается разложение:

.  (51)

Точно также на основе (44), имеем

  . (52)

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

Так, интегрированием уравнения (30) по правилу (44) при , , , , получается ряд для вычисления длины геодезической линии от точки экватора до текущей точки в форме:

(53)

  =  (54)

Аналогично, из (48) при , ,  следует:

  ,  (55)

    , (56)

с рекуррентной формулой

  ,  (57)

Из (53) при образуется отрезок ряда, который с параметром в коэффициентах k поистине стал классическим и приводится практически во всех учебных пособиях по сфероидической геодезии. Но никто не указывает, что его автором является Ориани Б. В лучшем случае, его относят к способу Бесселя и оставляют безымянным, в худшем –  приписывают Бесселю Ф.

Бессель Ф. при выводе формулы для вычисления S  исходит из тех же лежандровских условий, что и Ориани Б., но применяет другой модуль . Поэтому формулы Бесселя получаются из (53) при по разности , используя коэффициенты с параметром .

Множество введенных элементов для характеристики формы эллипсоида, привело к разным видам построения коэффициентов разложений. Так с помощью подстановок , равенство (55) приводится к зависимости:

  (58)

в которой коэффициенты формируются двумя способами. При группировке членов по степеням , получаются соотношения:

  ,  ,  (59)

а группировкой по степеням :

    (60)

Как видим, выражения (59, 60) имеют более сложную и громоздкую форму по сравнению с (56) и могут использоваться только при табличных вычислениях.

Аналогично, ряд (53) преобразуется к виду

(61)

Из (61) вытекают формулы Андреева М., которые он вывел интегрированием по переменной А.

Результаты Романовского В.А. следуют из (58) с коэффициентами (60), для вывода которых он применил две подстановки. Губений К. для вычисления коэффициентов (59) составил таблицы по параметрам эллипсоидов Бесселя, Хейфорда и Красовского. Если в (58) по зависимостям сферического треугольника Лежандра преобразовать функции к приведенной широте, то получим формулы Бахвалова, а также и Сержанова без использования аппроксимации поверхности сфероида шаровыми поясами.

Подстановка , выполненная по закономерностям (50), ряд (53) приводит к виду

  (62)

использованному Гельмертом Ф. с параметром при выводе формул для решения ПГЗ. Этим же способом на основе (55) определяем

.

Лежандр А. при выводе использует не параметр k, а второе сжатие . Его формулы вытекают из (62) при . Ганзен П. для вычисления коэффициентов применил подстановку и получил разложения по степеням . Основной недостаток формул, полученных Хандриковым М.Ф. и Вировцем А.М., заключается в том, что они положительному направлению геодезической линии на эллипсоиде поставили в соответствие отрицательное направление дуги большого круга на сфере. Указываются недостатки алгоритма Питмана, использующего переменную  .

Аналогично строятся формулы для вычисления длины геодезической линии по способу Баховена. Показано, что для вычисления длины геодезической линии наиболее перспективным является вывод формул по способу Лежандра с параметром .

Интегрированием уравнения (31) получена модель для определения долготы от точки до текущей точки 

  , 

  ,    (63)

Вычислить можно и по зависимостям

  ,

По этим соотношениям строятся рекуррентные формулы

. (64)

Введением в (63)  величины получаем равенство

  ,  , (65)

из которого следуют выражения, предложенные Ориани Б. для вычислений

,

  , , . 

Для определения Ориани не выводит общих формул (64, 65), а приводит выражения для при , являющиеся следствием (64), и их разложения по степеням эксцентриситетов.

Разность долгот между точками и вычисляется по правилу

,  , . 

Губений К. для использования формулы (63) коэффициенты вычислил по параметрам эллипсоидов Бесселя, Хайфорда и Красовского. Если в (63) отбросить члены с , а сферическую величину преобразовать к приведенной широте, то воспроизведем результаты Бахвалова и Сержанова.

Решение уравнения (31), представленное по синусам кратных дуг, имеет вид

.  (66)

Для вычисления справедливы все выведенные выше формулы.

       В отличие от (63) и (66), Ганьшин В.Н. общую модель формирует по другим закономерностям с введением большего числа вспомогательных величин, используя обычные и двойные факториалы, что усложняет не только процесс вычислений, но и структуру модели. И в данном случае с помощью соотношений и закономерностей (33, 50) из (66) получается формула для вычисления долготы от вертекса до текущей точки с прежними выражениями и :

  .  (67)

На основе (67) строятся формулы, выведенные Ориани Б., Ганзеном П., Лежандром А., Андреевым М., Ганьшиным В.Н., Гельмертом Ф., Хандриковым М., Вировцем А.М. и др. Из проведенного анализа следует, что все рассмотренные формулы для вычисления долготы являются следствиями обобщенной модели Ориани, записанной в виде (63) или (66).

Выведена общая модель для вычисления долготы по способу Баховена, из которой как следствие получаются формулы Баховена и Мак-Коу. Однако простое сравнение общих членов, выведенных по способам Лежандра и Баховена, показывает, что коэффициенты в разложениях Баховена в раз больше соответствующих коэффициентов в разложениях Лежандра, что при составляет величину десятичного разряда.

В связи с этим для определения долготы предлагается использовать разложение (66) с точными выражениями коэффициентов через сжатие со знаменателем , полученными по рекуррентным формулам (64)

, , 

, .  (68)

Значение можно приближенно определить по зависимости .

Таким образом, построенные математические модели для вычисления долготы и длины геодезической линии, способствовали проведению исследований методов решения ГГЗ на более высоком научном уровне.

В третьем разделе рассмотрены теория и методы построения математических моделей и вывода оптимальных формул высокой точности решения ГГЗ на поверхностях эллипсоида и шара, включая предельные расстояния.

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

Так как функция при не определена, то для исключения этого случая Бэр Х. предложил использовать половинные и четвертичные аргументы по схеме: если , , то . Показано, что вместо схемы Бэра удобнее использовать тождество , на основе которого по формулам Деламбра получены выражения для вычислений , , Этим способом выведены формулы  для решения прямой и обратной геодезических задач  на шаре.

Для вычисления длины геодезической линии между точками и Ориани построил формулу на основе (53) по зависимости . Так как при решении ПГЗ величина является заданной, то для определения , используя теорему Лагранжа, он этот ряд обратил, сохранив в разложениях величины по порядок включительно. Существенным недостатком этой неитеративной формулы является громоздкость и сложность ее выражений.

Впервые алгоритм для прямого вычисления с началом отсчета от вертекса опубликовал в 1789 г. Лежандр А., используя величины и второе сжатие для образования коэффициентов. И только в 1877 г. Винтерберг Н. предложил для обращения ряда использовать не разложение , а ряды с началом отсчета от экватора или вертекса . В этом случае, если исходный ряд имеет вид

  (69)

то и после обращения он сохраняет прежнюю структуру

    (70)

       

Отрезок ряда (53) автором представлен в степенной форме по основанию , а простые и компактные выражения коэффициентов (54) с параметром преобразованы к виду

(71)

позволяющему для любого эллипсоида производить вычисления “с конца к началу” без фиксации промежуточных результатов с сохранением в ячейках памяти наименьшего числа данных , обеспечивая погрешность , где 

Обращением (71) получена формула

  (72)

где с погрешностью .

Для вычисления Гельмерт Ф. по разложению (70) образовал разность , а затем применил формулу разности синусов. Ганьшин В.Н. для вычисления использует график и специально составленные таблицы. Кроме рассмотренных способов для определения применяются методы итераций и дифференциальных поправок. Из (70) при следует: , тогда Морозов В.П. вводит вспомогательную величину , которая не упрощает конечную формулу.

Установлено, если в формуле (66) отбросить член , то долгота определится с погрешностью , что при вычислении длины геодезической линии вызовет ошибку мм. Пренебрегая этим членом и вынося общий множитель в за скобку, получим равенство

, (73)

с прежними зависимостями для , но с величинами . На основе (73) образуется формула для вычисления разности долгот между точками и .

При решениях геодезических задач на поверхности эллипсоида и в пространстве используют постоянные величины, которые приведены с числовыми значениями по параметрам эллипсоида Красовского:

,  ,  ,  , , ,

1)   2)   3)  ,

4) , 5),

6) , .  (74)

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

С использованием полученных зависимостей выведена математическая модель решения ПГЗ на большие расстояния в виде системы формул, реализуемых в следующей последовательности:

  1.  
  2.   если   то
  3. ,
  4.   5) ,
    1.   7)
  1. ,

9)


10)  11) 

  1. , 13)   14) ,

15) 

16)  17) 

18)  ,

19) .  Если в 2) , то

20)       (75)

Для решения ОГЗ разность по разложению (53) с применением тригонометрических тождеств представлена в форме

,

где ,

Проанализированы способы преобразования этой формулы, примененные Райнсфордом Х., Боурингом Б., Морозовым В.П., указаны их недостатки. Для упрощения выражений коэффициентов используется зависимость (43) при . На основе связей между коэффициентами ,, получено более компактное удобное для вычислений выражение, которое рекомен-дуется для решения ОГЗ с погрешностью мм. Таким же способом преобразовано и уравнение (73) для разности долгот.

Для решения ОГЗ на большие расстояния предложена система формул:

1)

2)

3)

4) , 

5)

6)  

7)  

8)

9)

10)

11)

12)

13)

14)

15)

16)

17)

18)

19)

20)

21)

22) 23) 

24) ,

25) ,

26)   27)

28) ,

29)

30)

31)

32)

33) .(76)

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

Кроме метода итераций, начиная с 1950 г., было предложено несколько вариантов формул для неитеративного решения ОГЗ, которые выводятся обращением ряда для разности долгот, выполненным разными способами. Сравнительный анализ таких алгоритмов выполнен  Даскаловой М.. Доказано, что различие между формулами заключается только в использовании разных аналитических выражений для одних и тех же величин. При сохранении в разложениях величин порядка такие формулы принимают громоздкую форму и как итеративные не решают задачу при предельных расстояниях и  требуют специального исследования. К тому же, формулы сферической тригонометрии, применяемые для решения, при теряют смысл. Такая область, согласно исследованиям Гельмерта Ф., определяется условиями

  (77)

Доказано, что принципиальное отличие основных способов решения ОГЗ (76) от решения ОГЗ при предельных расстояниях состоит только в применении разных формул для вычисления азимута .

Основная концепция построения таких решений базируется на выводе формул для приближенного вычисления азимутов и разработке способов, приводящих их результаты к дальнейшему уточнению. Рассмотрены два случая, когда точки и на сфере, соответствующие точкам и на эллипсоиде, являются: 1) антиподными; 2) близкими к антиподным.

При антиподном расположении точек и на вспомогательной сфере и для определения применяют уравнение

  (78)

Решение уравнения (78) Морозов В.П. выполняет методом простой итерации, Зорский З. – методом Ньютона, используя коэффициенты , выраженные через первый эксцентриситет. Боуринг Б. для выражения коэффициентов применил сжатие , а для исключения итераций отрезок ряда (78) обратил. Полученная им формула не содержит информации о ее точнос-ти и не является простой для вычислений.

В работе приведен вывод формулы с учетом в уравнении (78) величин порядка . Для выражения коэффициентов использовано третье сжатие . В результате получено неитеративное, более простое решение ОГЗ, которое применяется при , и  обеспечивает вычисление результатов с погрешностями , :

  1) 2) 

  3) , 4)

5) ,

6) , ,    (79)

При точках, близких к антиподным, Шмидт Х. рекомендует использовать формулу Гельмерта, а Боуринг Б. предлагает зависимость

  (80)

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

Принимая (80) за приближенное выражение , разработана математическая модель решения ОГЗ:

1) ,  2),  3)

Для начального приближения:

4) , 5)

6) , 7) .

Для дальнейших итераций:

8) 9) 10) 11)   12) 13)

14) , 15) ,

16) , 17) ,

18)         19)  ,

20)   21)   

22)  ,  23)  

24) , 25) ,  , 

26) ,  , 

27) . (81)

Приведенная модель решения ОГЗ отличается от предложенной Боурингом Б.: начальным приближением (4, 5), уменьшающим число итераций при определении ; более простыми зависимостями (6, 7) или (14, 15) для вычисления азимута пересечения геодезической линии с экватором; формулами (16, 17), приводящими к квадратичной сходимости итерационного процесса; формулами (18, 19) для определения прямого и обратного азимутов и (23–26) для вычисления длины геодезической линии.

Достоинства этого способа решения ОГЗ проиллюстрированы примерами, вычисленными по (81), в сравнении с формулами, предложенными Морозовым В.П., Саито Т., Шмидтом Х., Русиным М.И. и др. Об эффективности способов решения ОГЗ, предложенных Саито Т. и Шмидтом Х., судить трудно, так как все приведенные ими примеры решаются по основным формулам (76). Выполнены исследования критериев по применению основных формул (76) и формул решения ОГЗ (81) при предельных расстояниях.

Итак, по разработанной методике выведены формулы высокой точности для эффективного решения ГГЗ на поверхностях эллипсоида и шара на большие расстояния, включая предельные.

В четвертом разделе рассмотрена теория и методология решения ГГЗ на основе аппроксимаций эллиптических  радикалов и интегралов с помощью рациональных дробей, полиномов Чебышева, функций Якоби, методом экономизации и численными методами.

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

В работе доказана тождественность дробно-рациональных приближений бинома , построенных с помощью подходящих дробей цепной дроби Лагранжа и методом геометрической прогрессии. Действительно, если принять, что в общем разложении

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

  (82)

из которой при , , получаем выражение второй подходящей дроби цепной дроби Лагранжа, которое Беспалов Н.А., Хитров Б.Ф., Куштин И.Ф. использовали при выводе формул для решения геодезических задач.

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

Беспалов Н.А. с помощью формулы Лагранжа полученное выражение привел к обычному отрезку ряда

(83)

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

Автором коэффициенты формулы Ганзена, имеющей вид (83), преобразованы способом геометрической прогрессии к дробному виду

мм,

с которыми разложение (83) имеет не только более компактную форму, но и на два порядка точнее. При этом, для вывода не требуются ни цепные дроби, ни интегрирование приближенного равенства (82), ни использование ряда Лагранжа. Аналогично, в качестве альтернативной, с помощью разложения Ганзена строится формула для вычисления разности долгот.

       Этим способом к дробно-рациональному виду приводятся не только коэффициенты ряда, но и сам ряд. Так преобразованием (71) получено

.

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

Кроме дробно-рациональных приближений эллиптических радикалов, для их вычислений применяют преобразование, получившее название экономизации. Морозов В.П. и Беспалов Н.А. использовали его для вычисления основных сфероидических функций, радиусов кривизны и упрощения выражений коэффициентов. Сущность процесса экономизации степенных рядов с помощью полиномов Чебышева заключается в следующем. Так как чебышевские многочлены подстановкой сводятся к виду

(84)

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

    (85)

где – целая часть числа ,  с учетом (84) представляет общую закономерность для построения полиномов Чебышева. В связи с этим появляется возможность анализа процесса экономизации с позиций тригонометрии. Так как , то исходный отрезок степенного ряда представляет разложение по степеням

    (86)

Соотношение устанавливает тождественность разложений по полиномам Чебышева и косинусам кратных аргументов

    (87)

Но зависимость (87) получается из (86) с помощью тождеств, выража-ющих степени через , поэтому результаты по тождественным между собой формулам будут определяться с одинаковой точностью. Их погрешность можно оценить по остатку исходного ряда (86): . Но отрезок ряда (86) подстановкой выражения , полученного из (85), приводится к виду

(88)

Формула (88) является не только тождественной (87), но и имеет равный с ней последний член. Поэтому Эта зависимость указывает, что ряд (87) по сравнению с (86) сходится быстрее, а последний член в равенстве (87), представленном по кратным аргументам , в раз меньше последнего члена (86), представленного по степеням .

При построении рабочих формул последние члены в отрезках рядов (86) и (88) отбрасывают и используют их вместе с остатком ряда для оценки достоверности оставшихся выражений. В результате получаются неравноточные между собой формулы. Для вычислений используется более точная (88) без последнего члена, которая и называется экономизированной.

       Этот способ экономизации основан на тригонометрических свойствах и не требует использования символики полиномов Чебышева. Анализ процесса экономизации на основе тригонометрических тождеств приводит к его обобщениям, позволяет убедительно и просто установить равносильные между собой алгоритмы; выявить причины, приводящие к неравноточным выражениям; указать способы преобразования их к тождественному виду. Такие исследования с применением полиномов Чебышева не только усложняют этот процесс, но и приводят иногда к неверным заключениям.

Обобщение процесса экономизации разложений эллиптических радикалов обусловлено тем, что тригонометрические тождества

 

        

  , 

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

При современной вычислительной технике действие извлечения квадратного корня не представляет никакой проблемы. Поэтому для вычисления значений основных сфероидических функций и радиусов кривизны не требуется никакая их аппроксимация: ни разложение в ряд, ни экономизация полученного отрезка ряда, ни представление в дробно-рациональном виде. Однако такие преобразования будут востребованными и иметь научную и практическую ценность при упрощениях результатов интегрирования эллиптических радикалов. Возможность процесса экономизации обусловлена теми же свойствами разложений по тригонометрическим переменным, а его реализация осуществляется с помощью соотношений (45–47).

Например, если за сумму ряда (48) принять частичную сумму первых его членов

  (89)

где то погрешность формулы (89), можно оценить по приближенной величине остатка ряда

Для проведения экономизации из тождества (45) получаем

. (90)

Подставив (90) в последнее слагаемое равенства (89) и пренебрегая членом с , приходим к экономизированному отрезку ряда

с погрешностью экономизации

  .

Аналогично проводится экономизация обобщенного эллиптического интеграла с помощью тождеств (46, 47). Рассматриваются способы экономизации, которые применили при преобразованиях формул Морозов В.П., Винкетти Т., Вольфрум О. Выполняется экономизация формул Вировца, повышающая их точность на порядок при прежнем объеме вычислительных операций. Этот метод использован автором при построении рабочих формул решения ГГЗ (75): 6, 8, (76): 24, 33, (81): 26.

Выполнен анализ вывода формул для решения ГГЗ с помощью функций Якоби. Полученные Якоби результаты не содержат в выражениях введенных им эллиптических функций, кроме одного параметра , определяющегося приближенно, независимо от них. Поэтому аппарат эллиптических и тэта-функций использовался как средство для проведения математических преобразований. В работе показано, что такие результаты несопоставимо эффективнее достигаются с помощью биномиального ряда. Например, если в формулах Лежандра второе сжатие выразить через первое, то после элементарных преобразований сразу получаются формулы Форсайта вместо сложнейших выводов с применением функций Якоби. Аналогично, выразив коэффициенты формул Лежандра через параметр , приходим к результатам, полученным Винтербергом с помощью эллиптических функций Якоби.

Из проведенного анализа сделан вывод о нецелесообразности применения функций Якоби для вывода формул в форме рядов, сложных по своей структуре и ведущих к громоздким преобразованиям.

Таким образом, поиски альтернативных решений путем представления эллиптических радикалов в дробно-рациональном виде, с использованием функций Якоби, картографических проекций и численных методов анализа, не привели к новым, более эффективным алгоритмам по сравнению с полученными в форме рядов.

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

При обработке пространственных геодезических сетей, решении ГГЗ между точками в пространстве приходится от геодезических координат точки переходить к ее прямоугольным геоцентрическим координатам и обратно. Проблеме преобразования этих систем координат, начиная с 1958 г. и до настоящего времени, в нашей стране и за рубежом посвящены сотни  публикаций.        Основные сложности здесь возникают при установлении связей между пространственными геодезическими координатами и прямоугольными координатами . Для определения долготы чаще всего приводится формула без описания способа вычисления . Один из возможных вариантов определения долготы по значению предложен автором:

1) 2)  ,

3) . (91)

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

1)   2)  

3)если ; ,если.  (92)

Для вычисления широты обычно используют формулу 

, (93)

которая в функциях с приведенной широтой принимает более простой вид

, ,  , . (94)

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

Используя аналитические и геометрические методы, установлены границы корня уравнения (93) при и уравнения (94):

  , ,  . 

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

Так как на любом референц-эллипсоиде е2/2<0,004 , то у чисел Т1 и Т2, представленных приближенно с помощью биноминального ряда , после их округлений окажутся одинаковыми не менее двух первых значащих цифр. Данное свойство имеет место и для всех точек, заключенных между ними. Следовательно, с такой точностью определяется корень уравнения и про-межуточная точка в остаточном члене, что дает возможность получить хоро-шие оценки для анализа итерационных процессов, принимая за промежуточную точку один из концов отрезка. Благодаря этому свойству можно сравнивать между собой погрешности как одного, так и разных алгоритмов относительно любых точек сегмента [T1; T2], чего нельзя было сделать ранее без решения данного вопроса. На основе этих положений автором и проведены исследования.

Решение уравнения (94) выполняется методом простой итерации по рекуррентной формуле

,  , ,  ,  .

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

По этой методике проведены исследования итерационных процессов определения широты по уравнениям, полученным из (94): с помощью тождества , приводящего к двусторонней линейной сходимости; подстановкой , приводящей к квадратичной сходимости; путем введения вспомогательного угла в функциях с приведенной широтой и геодезической широтой, рекомендованными Госстандартом России. Определение промежутка изоляции корня  позволило для вычисления широты применить метод хорд.

С помощью полученных результатов выполнен сравнительный анализ как разных способов решения уравнения (94), так и формул, которые предложили для вычисления широты Морозов В.П., Урмаев М.С., Боуринг Б.Р., Герасимов А.П., Машимов М.М., Буткевич А.В., Бопп Г. и Краусс Г., Винкетти Т., Баландин В.Н., Сорокин Н.А. и др. Согласно отрезку , проведенные исследования являются обоснованными только для точек с высотами . В работе уменьшен промежуток изоляции корня, что позволило вывести неитеративную формулу определения геодезической широты для точек с высотами и рекомендовать следующий вычислительный алгоритм:

1)   2)   3) 

4)  5)  (95)

При положительных высотах наибольшая погрешность широты на высоте . Если, то . При отрицательных высотах по мере приближения к центру эллипсоида погрешность возрастает: при ; при ; при .

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

В качестве переменной t чаще всего используется тангенс геодези-ческой или приведенной широт  (Ганзен Г., Сорокин Н.А., Фарао Ромао М., Пауль М., Пенев П., Фукусима Т. и др.). Основным недостатком таких решений является громоздкость их выражений, приводящих не только к большому объему вычислений, но и к действиям над числами порядка . Кроме этого, при решениях алгебраических уравнений четвертой степени по способу Феррари в его формулах в том или ином виде содержится операция вычитания двух близких по величине чисел, ведущая к потери верных значащих цифр. Следовательно, гораздо эффективнее при вычислениях применить систему формул (95). Полученные автором алгоритмы использованы в  следующих усовершенствованных математических моделях решения ГГЗ в пространстве.

Математическая модель решения ПГЗ:

1) ,

3) ,

2)

4)  

 

(96)

Математическая модель решения ОГЗ:

1)  ,

2)  ,

3) ,

 

4)  а)  б)

в) если ,  то

  если ,  то (97)

Заключение

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

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

       При больших расстояниях искомые величины сфероидического треугольника определяются с помощью рядов Фурье, построенных по элементам вспомогательного сферического треугольника. Коэффициенты в этих разложениях представляются по степеням параметра сжатия, что увеличивает скорость сходимости рядов. Эффективность решения ГГЗ с высокой точностью в этом случае достигается разными способами.

       На основе разложения обобщенного эллиптического интеграла обоснована эффективность построения математических моделей решения ГГЗ по способу Лежандра и выведены формулы в степенной форме для вычисления длины геодезической линии S  и длины дуги на сфере , вычислительные операции по которым осуществляются без введения вспомогательных коэффициентов с параметром и без фиксации промежуточных результатов. Следует отметить, что модель разложения обобщенного эллиптического интеграла можно использовать для вычислений не только длин геодезических линий, но и кривых, образуемых локсодромой, нормальным или геоцентрическими сечениями.

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

       Погрешность отрезка ряда уменьшается путем преобразования его к дробно-рациональному виду предложенным методом геометрической прогрессии и реализованным на упрощениях математических моделей решения ГГЗ.

       Применение разработанного и обобщенного тригонометрического принципа экономизации разложений эллиптических радикалов и интегралов для упрощения выражений математических моделей решения ГГЗ повышает точность формул, не увеличивая объема вычислительных операций. Этот метод можно использовать и для экономизации рядов с двумя переменными, которые нашли широкое применение в геодезии.

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

       Математические модели решения ГГЗ в пространстве усовершенствованы с помощью выведенной неитеративной формулы для вычисления геодезической широты с высокой точностью и предложенными формулами преобразования геодезических и полярных топоцентрических координат.

       Исследования альтернативных решений путем представления эллиптических радикалов в дробно-рациональном виде, с использованием функций Якоби, картографических проекций и численных методов не привели к новым, более эффективным решениям по сравнению с полученным в форме рядов.

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

       Выведенные математические модели решении ГГЗ можно использовать при геодезическом обеспечении обороны страны, решении фундаментальных и перспективных задач геодезии, геофизики, геодинамики, космонавтики и решениях других специальных задач. Теоретические и практические аспекты результатов диссертационной работы могут найти свое применение в задачах вычислительной и прикладной математики.

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

  1. Медведев П.А. Решение обратной геодезической задачи на большие расстояния по формулам со средней широтой и средним азимутом / П.А. Медведев // Науч. тр. Ом. с.-х. ин-та. – Омск, 1971. – Т. 80. – С. 66 – 75.
  2. Медведев П.А. О точности решения обратной геодезической задачи по формулам со средней широтой и средним азимутом / П.А. Медведев  // Землеустройство, планировка сельских населенных пунктов и геодезия: сб. науч. тр. / Белорус. с.-х. акад. – Горки, 1972. – Т. 86. – С. 238 – 244.
  3. Медведев П.А.  Развитие формул В.П. Морозова со  средними  аргументами / П.А. Медведев // Науч. тр. Ом. с.-х. ин-та. – Омск, 1972. – Т. 80. – С. 186 – 189.
  4. Медведев П.А. Применение рядов со средней широтой и средним ази-мутом при выводе формул для решения обратной геодезической задачи  / П.А. Медведев // Науч. тр. Ом. с.-х. ин-та. – Омск, 1972. – Т. 90. – С. 8 – 9.
  5. Медведев П.А. Исследование точности формул со средними аргументами К. Губения / П.А. Медведев // Науч. тр. Ом. с.-х. ин-та. – Омск, 1972. – Т. 105. – С. 64 – 69.
  6. Медведев П.А. Преобразование формул для решения обратной геодезической задачи к сфероидическим аргументам / П.А. Медведев // Науч. тр. Ом. с.-х. ин-та. – Омск, 1974. – Т. 120. – С. 3 – 8.
  7. Медведев П.А. Решение главной геодезической задачи с применением метода проектирования, предложенного Н.С. Боголюбовой / П.А. Медведев, А.А. Мурашов // Проблемы строительства в нефтегазоносных районах Тюменской области: науч. тр. Тюм. индустр. ин-та. – Тюмень, 1974. – Вып. 29. – С. 169 – 273.
  8. Медведев П.А. Улучшение сходимости рядов со средними аргументами для решения прямой геодезической задачи / П.А. Медведев // Науч. тр. Ом. с. х. ин-та. – Омск, 1975. – Т. 132. – С. 51 – 55.
  9. Медведев П.А. Исследование точности формул со средними аргументами для решения прямой геодезической задачи / П.А. Медведев // Науч. тр. Ом. с.-х. ин-та. – Омск, 1976. – Т. 154. – С. 56 – 60.
  10. Медведев П.А. Об эквивалентности формул со средними аргументами для решения обратной геодезической задачи / П.А. Медведев // Науч. тр. Ом. с.-х. ин-та. – Омск, 1978. – Т. 176. – С. 69 – 76.
  11. Медведев П.А. Совершенствование формул со средними аргументами для решения главной геодезической задачи / П.А. Медведев // Прикладная геодезия в сельском хозяйстве: сб. тр. / Ом. с.-х. ин-т. – Омск, 1979. – С. 20 – 25.
  12. Медведев П.А. Формулы для вычисления синусов окончательных результатов решения прямой геодезической задачи / П.А. Медведев // Исследование точности и математическая обработка сетей сгущения: сб. тр. / Ом.  с.-х. ин-т. – Омск, 1981. – С. 38 – 41.
  13. Медведев П.А. Вывод формул со средними аргументами для решения обратной геодезической задачи / П.А. Медведев // Изв. вузов. Геодезия и  аэрофотосъемка. – 1984. – № 5. – С. 20 – 27.
  14. Медведев П.А. Замечания к статье о решении прямой геодезической задачи на любые расстояния / П.А. Медведев // Геодезия, картография и  аэрофотосъемка: респ. межвед. науч.-техн. сб. – Львов, 1984. – Вып. 39. – С. 49 – 51.
  15. Медведев П.А. Формулы с улучшенной сходимостью для решения главных геодезических задач / П.А. Медведев // Геодезия и картография. – 1985. – № 1. – С. 14 – 16.
  16. Медведев П.А. О некоторых способах решения прямой и обратной геодезической задачи в функциях начальных аргументов / П.А. Медведев  // Геодезия, картография и аэрофотосъемка: респ. межвед. науч.-техн. сб. – Львов, 1990. – Вып. 52. – С. 42 – 46.
  17. Медведев П.А. Способ Бесселя или Ориани? / П.А. Медведев // Изв. вузов. Геодезия и аэрофотосъемка. – 1993. – № 1, 2. – С. 11 – 24.
  18. Медведев П.А. Определение широты точки способом решения алгебраического уравнения / П.А. Медведев // Изв. вузов. Геодезия и аэрофотосъемка. – 1993. – № 4. – С. 66 – 72.
  19. Медведев П.А. Исследование рекуррентных формул для вычисления широты при переходе от пространственных прямоугольных координат к геоде-зическим / П.А. Медведев // Геодезия и картография. – 1994. – № 6. – С. 8 – 14.
  20. Медведев П.А. Вычисление широты методом хорд при преобразованиях пространственных прямоугольных координат к геодезическим / П.А. Медведев // Изв. вузов. Геодезия и аэрофотосъемка. – 1994. – № 5. – С. 25 – 34.
  21. Медведев П.А. Исследование рекуррентных формул для вычисления широты при переходе от пространственных прямоугольных координат к геодезическим / П.А. Медведев // Геодезия и картография. – 1994. – № 7. – С. 8 – 12.
  22. Медведев П.А. Исследование рекуррентных формул для вычисления широты при переходе от пространственных прямоугольных координат к геодезическим / П.А. Медведев // Геодезия и картография. – 1994. – № 8. – С. 4 – 9.
  23. Медведев П.А.  Вычисление  широты  методом  касательных  при  переходе от пространственных прямоугольных координат к геодезическим  / П.А. Медведев // Геодезия и картография. – 1995. – № 12. – С. 4 – 6.
  24. Медведев П.А. О соответствиях между погрешностями аргумента и тригонометрических функций / П.А. Медведев // Геодезия и картография. –1997. – № 12. – С. 7 – 14.
  25. Медведев П.А. Анализ преобразований пространственных прямоугольных координат в геодезические: монография / П.А. Медведев. – Омск: Изд-во ОмГАУ, 2000. – 104 с.
  26. Медведев П.А. Решение прямой геодезической задачи по способу Ле-жандра / П.А. Медведев // Земельные ресурсы Сибири: изучение, управление, реформирование: сб. тр. / Ом. гос. аграр. ун-т. – Омск, 2002. – С. 203 – 213.
  27. Медведев П.А. Решение обратной геодезической задачи на сфероиде  при антиподном расположении точек на вспомогательной сфере / П.А. Медведев // Геодезия и картография. – 2008. – № 10. – С. 24 – 26.
  28. Медведев П.А. Вычисление длины дуги меридиана методом проективной геометрии / П.А. Медведев, А.П. Макаров // Геодезия и картография. – 2008. – № 11. – С. 16 – 23.
  29. Медведев П.А. Вычисление эллиптических интегралов с помощью дробно-рациональных приближений / П.А. Медведев // Изв. вузов. Геодезия и аэрофотосъемка. – 2008. – № 2. – С. 66 – 81.
  30. Медведев П.А. Определение погрешностей геодезической высоты, широты и долготы аналитическими методами / П.А. Медведев // Геодезия и картография. – 2009. – № 1. – С. 25 – 27.

Приложение А

(обязательное)

Рисунки к тексту диссертации

Рисунок А.7 – Сфероидические треугольники

Рисунок А.8 – Сферические треугольники Лежандра и Баховена







© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.