Адаптация параметров модели ионосферного радиоканала к текущим условиям без организации специализированных измерений тема автореферата и диссертации по физике, 01.04.03 ВАК РФ
Голыгин, Виктор Александрович
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Иркутск
МЕСТО ЗАЩИТЫ
|
||||
2012
ГОД ЗАЩИТЫ
|
|
01.04.03
КОД ВАК РФ
|
||
|
На правах рукописи
005017080
Голыгин Виктор Александрович
Адаптация параметров модели ионосферного радиоканала к текущим условиям без организации специализированных измерений
Специальность 01.04.03-радиофизика
" АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата физико-математических наук
1 о;.;:.. л ?С12
ИРКУТСК-2012
005017080
Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования "Иркутский государственный университет"
Научный руководитель: доктор физико-математических наук, доцент, Сажин Виктор Иванович
Официальные оппоненты:
Куркин Владимир Иванович, доктор физико-математических наук, зам. директора федерального государственного бюджетного учреждения науки "Институт солнечно-земной физики" Сибирского отделения РАН.
Кулижский Андрей Владимирович, кандидат физико-математических наук, старший научный сотрудник лаборатории распространения радиоволн НИИПФ ФГБОУ ВПО "Иркутский государственный университет".
Ведущая организация: Научно-исследовательский институт физики федерального государственного автономного образовательного учреждения высшего профессионального образования "Южный федеральный университет" (г. Ростов-на-Дону).
Защита состоится 28.05.2012 в 12:00 на заседании диссертационного совета Д.212.074.04 при ФГБОУ ВПО "Иркутский государственный университет", по адресу: 664003, г. Иркутск, бульвар Гагарина, 20.
С диссертацией можно ознакомиться в научной библиотеке Иркутского государственного Университета.
Автореферат разослан 20.04.2012.
Ученый секретарь диссертационного совета кандидат физико-математических наук, доцент
Мангазеев Б.В.
Актуальность темы
Ионосферное распространение декаметровых радиоволн по-прежнему привлекает внимание исследователей. Резервная связь, зондирование удалённых участков земной и морской поверхностей, местоопределение положения передатчиков, исследования ионосферы- эти основные задачи определяют сферу его использования. На такое распространение влияет множество факторов, значение которых в конкретной ситуации определить чрезвычайно трудно. Поэтому важной и актуальной задачей является адаптация, хотя бы частичная, параметров рассматриваемого радиоканала к текущей ситуации.
В настоящее время разработаны некоторые способы частичной адаптации параметров радиоканала за счёт применения определённых технических устройств. Так, для повышения уровня принимаемого сигнала используется ретрансляторы и приём в двух разнесенных по пространству точках [1]. Применяются также устройства, реализующие автоматический выбор частоты из заданного набора частот [2], системы с адаптацией антенных устройств к некоторым особенностям принимаемого сигнала [3]. Развивается и другая группа методов, использующая модель параметров радиоканала. В одном из методов этого направления осуществляется прямой перенос характеристик модели радиоканала, определяемых по измерениям на опорных радиолиниях, на близкорасположенную трассу. Данный метод обладает хорошей точностью в небольшой окрестности опорных линий, т.к. "переносятся" значения непосредственно характеристик сигнала, измеренные в данный момент времени. Однако метод требует развертывания достаточно широкой сети опорных линий с проведением на них сложных специализированных измерений. В России развернута такая сеть станций, оборудованных ЛЧМ зондами, и ведутся работы групп исследователей в Институте солнечной и земной физики СОР АН (ИСЗФ), Марийском государственном университете, Научно-исследовательском радиофизическом институте (НИРФИ), Ростовском государственном университете (РГУ) [4] -[6]. Естественно, развитие и расширение такой сети требует существенных затрат.
В другом, менее затратном методе этого направления адаптации, используется модель радиоканала, входной частью которой является среднемесячная модель ионосферы, а выходной- блок расчёта характеристик сигнала. Так как среднемесячные вариации основных параметров ионосферы описываются в настоящее время достаточно удовлетворительно, то задачей адаптации является, в данном случае, уточнение значений параметров модели на конкретную ситуацию, обеспечивающее повышение точности соответствия рассчитываемых характеристик сигнала измеряемым. Т.к. станции, непосредственно измеряющие параметры ионосферы, имеются лишь в небольшом числе пунктов земной поверхности, то разрабатываются также способы получения оперативной информации об ионосфере на основе при-
ближённых решений обратной задачи распространения радиоволн на наклонных трассах. Здесь можно указать работы, проводимые в Институте земного магнетизма, ионосферы и распространения радиоволн (ИЗМИРАН) [7], ИСЗФ [8], Иркутском государственном университете (ИГУ) [9], РГУ [10], [11], Санкт-Петербургском государственном университете (СПБГУ) [12] и работы, выполненные в некоторых зарубежных учреждениях [13], [14]. Точность такого подхода (назовём его модельным) уступает точности метода прямого моделирования, однако адаптация может быть выполнена для более протяженного пространственного региона, примыкающего к опорной трассе, за счёт возможностей уточнённой модели ионосферы.
В то же время трудность широкого применения этого метода связана с тем, что для получения корректирующей информации требуется организация на опорной трассе специализированных измерений некоторых характеристик сигнала при работе передатчика на ряде частот.
В настоящей работе предлагаются способы адаптации, которые не требуют для своего применения организации специализированных измерений, что сделает адаптацию более доступной широкому кругу пользователей. В первом способе используется наличие в декаметровом диапазоне большого числа произвольных радиостанций. Это и вещательные станции, и станции точного времени, и радиолюбительские маяки, и некоторые другие станции. Они работают на определённом наборе частот в режиме непрерывного сигнала и известно их местоположение.
Кроме того, сейчас на земной поверхности развернута большая сеть специализированных приёмников сигналов глобальной спутниковой радионавигационной системы GPS, для которых определяются и размещаются для свободного доступа в сети Internet данные по полному электронному содержанию ионосферы (ПЭС). Во втором способе адаптации рассматривается использование данных по ПЭС для уточнения входных параметров модели радиоканала. Далее, на основе предлагаемых способов показывается возможность адаптации параметров модели радиоканала к текущей ситуации в некоторых задачах связи и местоопределения.
В соответствии с вышеизложенным формулируются цель работы и задачи исследования.
Цель работы п задачи исследования
Целями настоящей работы являются исследование некоторых вопросов адаптации параметров модели ионосферного радиоканала к текущим условиям на основе коррекции входной её части- блока задания среды распространения сигнала и разработка способов её выполнения по оперативной информации, для получения которой не требуется организации специализированных измерений.
Для достижения данной цели необходимо было решить следующие задачи:
1. Разработать способ получения корректирующей информации, для получения которой не требовалось бы организации специализированных измерений
2. Разработать способ адаптации параметров модели канала к текущей ситуации с использованием этой корректирующей информации
3. Разработать способ уточнения входных параметров модели радиоканала в текущей ситуации по величине полного электронного содержания ионосферы
4. Выполнить оценку эффективности разработанных способов для различных трасс и условий
5. Предложить способы адаптации входных параметров модели ионосферного радиоканала к текущим условиям в практических задач связи и местоопределения
Научная новизна работы заключается в следующем:
1. Разработан способ адаптации модели ионосферного радиоканала к текущим условиям по значениям МПЧ, определяемым на радиолиниях с постоянной рабочей частотой и известным местоположением передатчика в моменты радиовосхода и радиозахода на них.
2. Предложен способ значительного уточнения модельных значений критической частоты ионосферы в произвольном протяженном пространственном регионе по величине полного электронного содержания ионосферы.
3. Предложены способы адаптации входных параметров модели ионосферного радиоканала к текущим условиям в задачах связи и местоопределения, обеспечивающие повышение энергетического потенциала радиолинии, точности расчёта углов места, модовой структуры и дальности распространения сигнала в канале.
На защиту выносятся следующие положения:
1. Разработанный способ уточнения входной части модели ионосферного радиоканала по наблюдениям за сигналами реперных радиостанций существенно снижает погрешность задания критической частоты ионосферы в текущих условиях.
2. Предложенный способ коррекции исходных значений критической частоты ионосферы в модели радиоканала по данным о полном электронном содержании ионосферы, позволяет значительно уточнить текущие значения критической частоты ионосферы в определённом пространственном регионе и вдоль протяженных радиолиний.
3. Предлагаемые способы адаптации входных параметров модели ионосферного радиоканала к текущим условиям в задачах связи и местоопределения обеспечивают существенное повышение энергетического потенциала радиолинии, точности расчёта углов места, модовой структуры и дальности распространения сигнала.
Достоверность результатов, представленных в диссертационной работе обусловлена использованием физических обоснованных методов, про-
веренных большим объёмом численного моделирования с использованием хорошо апробированной модели радиоканала. Основные результаты работы находятся в количественном и качественном согласии с результатами экспериментальных данных, полученных в ИСЗФ СО РАН и ИГУ, а также ряда исследований, опубликованных в периодической печати.
Практическая ценность работы заключается в том, что применение разработанных способов и полученных результатов может повысить эффективность решения некоторых задач связи и местоопределения. Кроме того, уточняемая модель радиоканала может быть использована также при радиозондировании ионосферы.
Внедрение результатов
Результаты, полученные в диссертации использовались и используются при выполнении следующих тем и программ:
1. Госбюджетная научно-исследовательская тема "Исследование эффектов различных неоднородностей при распространении радиосигнала через многомасштабную неоднородную среду", выполняемая ФГОУ ВПО "Иркутский государственный университет" в 2008-2011гг.
2. Проект 14.740.11.0078 федеральной целевой программы "Научные и научно-педагогические кадры инновационной России" на 2009-2013гг.
Личный вклад автора
Автором создан программный комплекс, разработана методика и проведён весь объём численного моделирования, представленный в работе, внесён весомый вклад в реализацию установки, организацию и проведение наблюдений на ней. Автор также проводил сравнение данных численного моделирования с экспериментальными данными, определяющим образом участвовал в формулировании результатов и выводов работы.
Апробация работы
Основные результаты и выводы, приведенные в диссертации, докладывались и обсуждались на конференциях молодых учёных в рамках Байкальской научной молодежной школы по фундаментальной физике БШФФ-2003, БШФФ-2004, БШФФ-2005, БШФФ-2006, БШФФ-2007, БШФФ-2009, БШФФ-2011, Иркутск; на конференции аспирантов ИГУ, Иркутск; на семинарах кафедры радиофизики ИГУ, Иркутск.
Публикации
По теме диссертации опубликовано 17 научных работ: 3 - в журнале "Геомагнетизм и аэрономия"(входит в перечень ВАК), 1 - в электронном журнале "Исследовано в России", 13 - в трудах и научных сборниках Байкальской молодежной научной школы по фундаментальной физике.
Структура и объём диссертации
Диссертация состоит из введения, четырех глав, заключения и списка использованной литературы. Работа изложена на 115 листах, включает 32 рисунка и 6 таблиц.
КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ
В главе 1 рассматривается модель ионосферного радиоканала, реализованного в работе. Рассмотрим её кратко: она содержит два основных блока - исходный вариант и блок коррекции. Входная часть исходного варианта модели (см. слайд-рисунок №1), представляет среднемесячную модель ионосферы, в которой выходная часть- блок расчёта траекторных характеристик сигнала. В качестве модели ионосферы выбрана полуэмпирическая модель ионосферы (ПЭМИ), разработанная в ИГУ под руководством В.М.Полякова [15].
ИСХОДНЫЙ ВАРИАНТ -
МОДИФИЦИРОВАННАЯ ПЭМП
I
РАСЧЁТ ТРАЕКТОРНЫХ ХАРАКТЕРИСТИК СИГНАЛА
СРАВНЕНИЕ С ЭКСПЕР. ДАННЫМИ
БЛОК КОРРЕКЦИИ
ОПРЕДЕЛЕНИЕ ГоГ2
УТОЧНЕНИЕ Н1П
^ | ^
ПОЛУЧЕНИЕ ПОЛУЧЕНИЕ Го12
ДАННЫХ ПО ДАННЫХ ПО 1, ! МПЧ ПЭС'
Рисунок 1. Блок-схема модели радиоканала
Поскольку основной целью работы является разработка методов адаптации параметров модели радиоканала к текущей ситуации, то эта цель может быть реализована при использовании достаточно простого исходного варианта выходной части модели радиоканала, описывающего лишь основные, определяющие вариации параметров сигнала. Поэтому в выходной части модели радиоканала для расчета параметров сигнала использован известный метод характеристик (дающий строгое численное решение в рамках исходного приближения геометрической оптики).
Рассматривается распространение радиоволн в двумерной неоднородной изотропной ионосфере без учёта наличия в фоновой регулярной ионосфере неоднородностей различного масштаба. Численная реализация метода выполнена в соответствии с алгоритмами, хорошо апробированными в ИГУ [16]. Лучевые уравнения и уравнение для группового пути с выбором в качестве независимой переменной центрального угла 9, преобразованные из сферической системы координат в систему "выпрямленных" декартовых координат имеют следующий вид (ниже приведены используемые в уравнениях обозначения с поясняющим форму траектории сигнала рисунком):
г]и/ 1 г!г.
д£__
Ло 2 £ дг '
Рисунок 2. К лучевым уравнениям Я, 0 - сферические координаты, Я0 — радиус Земли, г = Я - Яо, х = Яо * е,
Р — групповой путь сигнала, е - диэлектрическая проницаемость ионосферы.
Дополнениями к известным реализациям метода характеристик, выполненными в настоящей работе, являются алгоритмы более оперативного расчёта МПЧ при односкачковом распространении на основе использования сплайн-интерполяции и расчёт МПЧ двухскачкового распространения, с включением косвенного способа учёта влияния случайных неоднородно-стей и шероховатостей земной поверхности в области отражения сигнала. Сопряжение с моделью ионосферы выполнено на основе кубической сплайн-интерполяции. Программа расчёта реализована в системе Делфи.
Ранее была проведена широкая проверка точности рассматриваемой модели декаметрового радиоканала [9]. Показано , что для периодов слабых и умеренных продольных градиентов на трассах средней протяженности в спокойных геомагнитных условиях модель хорошо применима для описания усреднённых за длительный временной интервал значений характеристик сигнала. Следовательно, она может быть эффективна при коррекции на текущую ситуацию.
В главе 2 рассматривается адаптация, используя данные наблюдений за сигналами реперных радиостанций. Мы обратили внимание, что в КВ диапазоне работает достаточно много радиостанций о которых известно, как правило, только точное местоположение и постоянная рабочая частота.
Такие линии мы рассматриваем как реперные. По наблюдениям за сигналами таких радиолиний, можно фиксировать значения МПЧ в моменты радиовосхода и радиозахода. Для фиксации МПЧ нами была реализована установка, блок-схема которой показана на рис.3.
Рисунок 3. Блок-схема установки.
Установка включает антенну, имеющую изотропную диаграмму направленности в горизонтальной плоскости, практически неизменную с частотой, устройство сопряжения, перестраиваемый по частоте приёмник и сопрягаемый с ним ПК с программой обработки.
Эта программа - Веасоп8ее [17] представляет собой анализатор сигнала, подаваемого по соединительному тракту на вход звуковой карты ПК. Обработанные данные выводятся в окна на экране (в каждом окне своя станция наблюдения). На рисунке 4 показан вид экрана при наблюдениях с помощью данной программы. Наличие сигнала соответствует засветке экрана, а его интенсивность- яркости засветки.
-1-1-1-=-
14:00 15:30 ЬТ
Рисунок 4. Вид экрана программы ВеасопБее
Для определения значения МПЧ реперной радиолинии на основе анализа временного поведения амплитуды сигнала по теоретическим расчётам работы [18] и результатов измерений, приводимых в [19], разработана следующая методика. Момент времени, соответствующий условию равенства
9
рабочей частоты и МПЧ трассы, фиксируется для середины участка резкого спада амплитуды от максимальных значений, соответствующих рефракционному прохождению сигнала, до минимальных, отвечающих рассеянию на мелкомасштабных неоднородностях на частотах выше МПЧ.
При наблюдениях на нашей установке, фиксируемый момент перехода через МПЧ соответствует резкому изменению яркости и интенсивности полос спектра сигнала, вплоть до пропадания засветки экрана. Полученные нами оценки длительности интервала перехода в сопоставлении с данными измерений МПЧ на среднеширотных радиолиниях показывают, что точность определения МПЧ по нашей методике составляет от +-0,3 до +-0,5МГц. Эта величина сопоставима с точностью измерения МПЧ при НЗ [20].
Получаемые значения МПЧ реперных радиолиний используются для уточнения в модели значений критических частот ионосферы. Теперь предстояло разработать метод использования этой информации для коррекции входной части модели радиоканала- модели ионосферы. Рассмотрим вначале методику уточнения по данным, получаемым на односкачковых радиолиниях.
Ранее в работе [21] указывалось на преобладающее влияние вариаций й)Р2 в области отражения на значения МПЧ в ПЭМИ для трасс средней протяженности в отсутствие продольных градиентов. Мы провели численное моделирование в большом объёме для ряда трасс различной протяженности и наличия умеренных продольных градиентов для оценки этого влияния. Получено, что наиболее существенное влияние на МПЧ радиолинии оказывают изменение критических частот в области отражения сигнала. На рисунке 5 представлены результаты оценки вариаций МПЧ в зависимости от вариаций параметров ПЭМИ в точке отражения для односкач-ковой трассы длиной 2500км с небольшим градиентом й>Е2, составляющем приблизительно изменение в 0,1МГц на 100км.
Видно, что влияние вариаций 1ЪЕ2 в средней точке трассы более чем в два раза превышает влияние вариации Ът¥2 и эффективной полутолщины слоя ¥2. В ходе моделирования получено, что при наличии умеренных градиентов для трасс от 1200 до 2700км степень влияние вариаций АэР2 остаётся преобладающей. Таким образом, подбирая в средней точке трассы значение й>Е2, обеспечивающее совпадение МПЧр с МПЧи, можно уточнить модель на текущую ситуацию в средней точке трассы. После этого по определённой методике выполняется коррекция значений foF2 во всей области ионосферного участка трассы.
Для дополнительного обоснования методики нами проведено численное моделирование, показавшее, что изменение положения области вершины скачка в исходной и скорректированной модели, составляет величину порядка 15-20 км, что соответствует изменению {о¥2 на величину не более 0,1 Мгц.
Рисунок 5. Влияние вариаций ионосферных параметров на значения МПЧ для трассы длиной 2500км.
Перейдём теперь к описанию методики коррекции модели ионосферы по данным двухскачковых радиолиний.
При расчёте МПЧ двухскачкового распространения используется метод косвенного учёта влияния отражающих свойств земной поверхности при расчёте МПЧ многоскачкового распространения, который в соответствии с работами [22, 23] заключается во введении при траекторных расчётах принципа равенства между собой МПЧ отдельных скачков и в целом МПЧ трассы.
Определение областей отражения на каждом из скачков при двухскач-ковом распространении выполняется на основе решения прямой задачи по расчёту МПЧ двухскачкового распространения (МПЧ2) в исходной модели без коррекции.
Далее их положение фиксируется, т.е. дальности скачков закрепляются, и в этих областях отражения независимо варьируются значения критических частот fm и /02, даваемых ПЭМИ. Вариации f0 проводятся до тех пор, пока в пределах погрешности измерений расчетные значения МПЧр на каждом скачке не станут равными измеренному значению МПЧ2„ радиолинии.
При этом пренебрегаем возможными изменениями геометрии скачков, а, следовательно, и изменением положения областей отражения сигнала в ионосфере на каждом скачке для скорректированной модели по сравнению с исходной. Степень этих приближений проверялась численным моделированием и в итоге получено, что они являются достаточно обоснованными.
Оценка эффективности коррекции входной части модели по данным о МПЧ реперных радиолиний выполнена с использованием данных ВЗ, полученных для пункта г.Иркутск в ИСЗФ, для других пунктов- из базы данных, размещаемой в сети Интернет.
Для односкачковой трассы обработаны данные наблюдений более чем 50дней МПЧ на трассе Новосибирск-Иркутск в невозмущенных магнитных условиях различных сезонов с 2002 по 2004гг. Улучшение точности описания критических частот в текущей ситуации оценивалось на основе сравнения с данными ВЗ в одном или другом конечных пунктах трассы. Данные ВЗ переносились с помощью ПЭМИ в область отражения сигнала на трассе.
Получено, что в ряде случаев коррекция модели позволяет значительно уменьшить расхождения с данными ВЗ. Вместе с тем, в утренний и вечерний периоды, когда наблюдаются более сильные продольные градиенты на реперной трассе, снижение отличий не такое существенное. Усредненные оценки эффективности уточнения критических частот показали, что рассмотренная коррекция снижает различия модельных значений /о с наблюдаемыми в текущей ситуации примерно на 70%. Эти результаты свидетельствуют о достаточно высокой эффективности методики коррекции и возможности ее практического использования.
Такую коррекцию можно проводить в определённом пространственном регионе с учетом радиусов пространственной корреляции отклонений величины критических частот от модельных значений. Как показывают наши оценки, эффективность коррекции сохраняется, если точка отражения сигнала на реперной трассе отстоит от точки отражения на "обслуживаемой" трассе (той, для которой выполняется коррекция) на небольшие расстояния, приблизительно до 300-400км. Это соответствует отличиям для азимутальных углов трасс в несколько градусов (до 10 градусов).
При оценке эффективности с использованием данных, получаемых на двухскачковых реперных радиолиниях, была выбрана система из двух трасс с близкими значениями азимутов, на одной из которых выполняется односкачковое распространение (Новосибирск — Иркутск), а на другой — многоскачковое (Москва — Иркутск). На рисунке 6 показана примерная геометрия скачков. Так как азимуты трасс почти совпадают, то в некоторый утренний период примерная геометрия скачков на них такова, что область отражения сигнала на трассе Новосибирск - Иркутск близко расположена к области отражения сигнала на втором скачке трассы Москва - Иркутск, что позволило для оценки эффективности уточнения значений ГоР2 по МПЧ2 использовать данные ВЗ в пунктах Иркутск и Новосибирск. Кроме того, измерения на такой системе трасс позволяют сопоставить степень уточнения ¿ЬР2 по МПЧ односкачкового распространения — МПЧ1 и МПЧ2. Получено, что абсолютное уточнение й>Р2 хотя и достаточно большое, но несколько уступает уточнению получаемому при использовании данных МПЧ1, составляя примерно 80% от последнего. Однако, важно подчеркнуть, что с использованием МПЧ2 можно уточнить значения й>Р2 в двух точках трассы, корректируя таким образом ход продольных градиентов в текущей ситуации.
z.
Рисунок 6. Приближенная геометрия скачков на трассах Москва-Иркутск и Новосибирск-Иркутск в утренний период
При разработке второго способа адаптации входной части модели радиоканала (глава 3) было обращено внимание на доступность и большую зону охвата измерений ПЭС в системе GPS. Была поставлена задача использовать эту информацию для адаптации модели.
Для того чтобы ПЭМИ отражала вклад в ПЭС участка профиля N(h) соответствующего высотам значительно больше максимума, в ПЭМИ включено описание области максимума параболой, с четверти полутолщины которой за ним используется экспоненциальная зависимость, идущая до больших высот. Часть параметров зависимости определяется при сопряжении с параболой, значения других подобраны на основе сравнения расчётных значений ПЭС по модели с значениями, получаемыми при усреднении данных измерений в системе GPS, размещаемых в сети Интернет. В силу того, что невозможно подобрать одно значение свободного параметра экспоненциальной зависимости, обеспечивающее удовлетворительное совпадение среднемесячных значений ПЭС с наблюдаемыми для всех гелиофи-зических условий, было выполнено разделение общей ситуации по сезонам (зима, лето, равноденствие), с выделением для каждого сезона дневного, ночного и переходного периодов. Таким образом, было подобрано 24 значения свободного параметра экспоненциальной зависимости.
Выполненное моделирование вариаций в ПЭМИ критической частоты fo и высоты максимума ионосферы показало, что вариации fo определяющим образом влияют на изменение ПЭС, рассчитываемого по модели (см рис 7). Степень изменения ПЭС составляет при этом в различных геофизических условиях до 40%-50%.
Также для исследования связи вариаций значений ПЭС с изменениями критических частот были рассмотрены временные месячные для одного часа наблюдений изменения значений ПЭС и fo в пункте г. Иркутск. Проведенное сравнение зависимостей показало их высокую степень коррелированное™. В целом получено, что коэффициенты корреляции fo и ПЭС составляют 0,8- 0.9, лишь в единичных случаях уменьшаясь до 0,7. На такую же высокую корреляцию для суточных и широтных зависимостей ПЭС указывается в работах [24, 25] соответственно.
\
40,035,0 -30,0-
• 25'° : / - - — 1er?
^ 20,0
\ 1 ; / \ю,о
- -5,0 ,
^ \ /
--•.•"O^HS- -- -
А/
мененпя функций *Ь
25 35
Рисунок 7. Оценка степени вариаций ПЭС при изменениях параметров модели
Используя, опять-таки, приближенный метод решения обратной задачи, появляется возможность, вычисляя значение ПЭС по модели в определённом пункте, подобрать значение foF2 для неё, пока рассчитываемое значение ПЭС не совпадёт с измеряемым. Для оценки эффективности такого, предложенного нами, способа уточнения foF2, проведена его апробация для некоторых среднеширотных пунктов российского региона. В качестве критерия использовано сравнение уточняемых значений foF2 с данными ВЗ, имеющимися для этих пунктов.
В целом, после усреднения по всем массиву из нескольких десятков выполненных уточнений fo получено, что в среднем снижение отличий в значениях foF2 по сравнению с ВЗ составило около 2 раз. Этот результат указывает на перспективность практического использования предлагаемого способа уточнения значений fo. Такое снижение абсолютного значения различий несколько уступает эффективности уточнения foF2 по наблюдениям за сигналами РРС. Вместе с тем, достоинством способа является его простота и доступность получения исходных данных, поэтому он может рассматриваться как дополнительный к методу определения foF2 по сигналам РРС.
Эффективность способа может существенно повысить использование для уточнения foF2 данных по ПЭС, получаемых на станциях с двухчас-тотными приёмниками непосредственно в регионе, для которого будет использоваться уточненная модель. В случае же использования данных по ПЭС из сети Интернет, надо учитывать задержку во времени получения данных. Показано, что при небольших временных задержках (до 3 часов) можно использовать экстраполяцию вперед значений ПЭС за предыдущие периоды.
Для протяженных трасс, коррекция по ПЭС может быть выполнена в нескольких точках, так чтобы восстановить ход foF2 вдоль трассы. Это даёт возможность уточнять ход градиентов foF2 для дальних трасс, а также отдельных участков на сверхдальних радиолиниях, если при этом отсутствуют подходящие реперные радиостанции. На рисунке 9 для примера показано сравнение эффективности уточнения foF2 по двум способам на протяженной трассе Москва-Иркутск. Видно, что эффективность уточнения по ПЭС хотя и ниже, чем для уточнения по РРС, но тем не менее достаточно высока.
Г... МГц .........
6,5
G 5.5
4.5 Л 3,5 3 2,5
Рисунок 8. Сравнение эффективности двух способов уточнения АэК2 на протяженной трассе (Москва- Иркутск).
В главе 4 описывается, как на основе уточнения параметров входной части модели радиоканала, может быть выполнена адаптация выходных параметров модели, в задачах связи и местоопределения. Рассмотрим вначале адаптацию в задачах связи.
В настоящее время имеются специальные системы КВ-связи, обеспечиваемые и поддерживаемые специализированным оборудованием. Но есть и потребители не имеющие таких возможностей, особенно в регионах с неразвитой инфраструктурой- месторождения нефти и полезных ископаемых, турбазы, экспедиции. Спутниковая связь в таких случаях не всегда доступна и значительно более дорогостояща. Кроме того КВ-диапазон используется при загоризонтной радиолокации морских поверхностей. Решаются в этом диапазоне и специальные задачи по обнаружению положения удалённого передатчика, т.н. задачи местоопределения. Для повышения энергетического потенциала линий предлагается использование на приёме остронаправленных антенн, с более точным нацеливанием их главного лепестка по углу места приходящего сигнала. Для этого значения углов вычисляются по корректируемой на текущую ситуацию модели.
В проведённых нами для различных трасс расчётах получено, что в некоторых случаях значения углов прихода, найденные по скорректированной на текущую ситуацию модели, могут отличаться от даваемых исходной
15
моделью на 7-8 градусов. Такое уточнение угла места соответствует, как показывают проведенные нами оценки, повышению энергетического потенциала линии до Юдб.
Как показывает численное моделирование, значения углов прихода в вертикальной плоскости существенно зависят и от вариаций Ьт в области отражения. Поэтому для повышения эффективности коррекции модели важно также при её проведении уточнить и значения Ит в текущей ситуации. По данным ВЗ это можно сделать, только если оно проводится в данной точке трассы, "перенос" значений Ьт из другой точки неэффективен, т.к. радиусы пространственной корреляции отклонений ЬтР2 от среднемесячных значений невелики. В такой ситуации для уточнения Ьт можно использовать значение МПЧ, определяемое ещё на одной реперной радиолинии. Беря значение й>Р2 из данных по МПЧ первой радиолинии, вариацией значений Ьт на второй подбираем лучшее соответствие для неё МПЧр с МПЧи.
Результат уточнения углов места максимален при совпадении азимутов реперной и обслуживаемой линий, он будет снижаться при различиях в несколько градусов для трасс средней протяженности.
На линиях ионосферной связи с передачей больших объёмов данных, коррекция модели ионосферы может быть использована для определения временного периода суток с одномодовым распространением. На рис. 9 для примера приведены данные расчёта такого периода односкачкового распространения для трассы Магадан-Иркутск в отдельный день эксперимента, когда на ней проводилось наклонное зондирование на скользящей частоте [20]. Результаты измерений представлены на рисунке сплошной линией с прямоугольниками, результаты расчёта по исходной модели- пунктирной линией, по модели с коррекцией по {о¥2 из данных ВЗ в Иркутске- более тонкой сплошной линией. Видно, что коррекция существенно улучшает описание периода одномодового распространения. Проделанные расчёты и сравнения с измерениями для 10 дней показали, что ошибка в определении продолжительности данного периода для скорректированной модели не превосходит 30 минут, а ошибка определения частоты одномодового распространения для отдельного часа снижается почти вдвое. Таким образом, эти данные демонстрируют эффективность уточнения модовой структуры на трассе.
Рисунок 9. Период одномодового распространения в отдельный день на трассе Магадан-Иркутск
¡08.11.33)
МПЧ1-МПЧ2, МГц
Д Экспериментальный участок цциомодоттого распр я ■ ^ ■ • Расчетный участок едномодоеого распр я до коррекции модели —— Расчетный участок одномодового распря после коррекции недели
Рассмотрим теперь адаптацию параметров ионосферного радиоканала в задачах местоопределения.
В частности, при ионосферном декаметровом распространении радиоволн в специализированных задачах находится положение радиопередатчика по измерению угла прихода в вертикальной плоскости сигнала от него и дальнейшему расчёту траектории, с использованием модели ионосферы, для определения дальности распространения. С целью более точного местоопределения передатчика нами предлагается адаптация входной части модели радиоканала на текущую ситуацию. В качестве примера эффективности подобной адаптации в таблице 1 представлены, в сравнении с данными измерений, приведённых в [АЗ], результаты расчёта углов на трассе Новосибирск-Иркутск в отдельные моменты времени для исходной модели ионосферы и с коррекцией по нашей методике значений £эР2. Коррекция проводилась для тех моментов, когда расчётные значения углов заметно отличались от измеренных. Ошибки в определении дальности получены расчётом траектории распространения сигнала из Иркутска в Новосибирск под углом выхода, равным измеренному значению угла прихода сигнала при распространении из Новосибирска.
Таблица 1. Сравнение результатов расчётов углов прихода в вертикальной плоскости и данных измерений на трассе Новосибирск-Иркутск
Время, Данные измерений Исходная модель Скорректированнаяя модель
ит (Зи' Р' Д1, мс ДО, км Р' мс ДО, км
02-00 12,1 12,1 0 0
03-00 10,5 10,6 0,1 18
0Ф00 13,4 12,4 -0,3 -107 12,9 -0,2 -51
05-00 12,7 11,2 -0,5 -169 12 -0,3 -87
06-00 12,4 10,7 -0,6 -195 11 -0,3 -110
Видно, что адаптация входной части модели радиоканала к текущим условиям приводит к существенному уточнению выходного параметра-дальности положения передатчика — заключающемуся в снижении ошибки до двух раз.
В заключение сформулированы основные результаты, полученные в работе:
Выполнена реализация исходной модели радиоканала, описывающей средние вариации дальности и времени распространения, углов прихода сигнала, максимально-применимых частот ионосферных радиолиний для двумерно-неоднородной изотропной регулярной ионосферы. Модель имеет хорошую применимость для описания средних значений характеристик сигнала при ионосферном распространении и поэтому может быть эффективна при использовании в задаче адаптации её параметров к текущим условиям.
Предложен способ фиксации МПЧ реперной радиолинии в моменты радиовосхода и радиозахода на ней по наблюдению временного поведения амплитуды сигнала.
В рамках численного моделирования определена связь вариаций параметров ионосферы и некоторых выходных характеристик канала- максимально-применимых частот радиолинии, углов прихода сигнала в вертикальной плоскости, времени распространения, а также связь вариаций параметров максимума электронной концентрации с вариациями полного электронного содержания ионосферы.
Разработан способ адаптации параметров модели канала к текущей ситуации по наблюдениям за сигналами реперных радиостанций на одно-скачковых и двухскачковых радиолиниях, в котором для получения корректирующей информации не требуется организации и проведения специализированных измерений
Разработан способ уточнения текущих значений йэР2 в протяженном пространственном регионе по данным о полном электронном содержании ионосферы. Способ приближенно учитывает средние вариации формы профиля N(11) в области максимума и на участке больших высот выше максимума ионосферы.
Выполнена оценка эффективности разработанных способов, показавшая их хорошую применимость для спокойных ионосферных условий. Для метода, использующего реперные трассы, применимость показана при наличии на них слабых и умеренных продольных градиентов электронной концентрации
Предложены способы адаптации входных параметров модели ионосферного радиоканала к текущим условиям в задачах связи и местоопреде-ления, существенно повышающие энергетический потенциал линии и точность определения углов места, модовой структуры сигнала и дальности его распространения.
Цитируемая литература:
1. Головин О.В. Декаметровая радиосвязь. М. : Радио и связь, 1990. 240 с.
2. Кузьмин Б.И. Адаптивные и автоматизированные системы связи. М. : Знание, 1994. 64 с.
3. Robert A. Monzingo, Thomas W. Miller // Introduction to Adaptive Arrays. Sci-Tech Publishing Inc, 2004.
4. Метод оперативной диагностики KB радиоканала / В.И. Кур-кин, В.Е. Носов, С.Н. Пономарчук, С.С. Савков, J1.B. Чистякова // Исследования по геомагнетизму, аэрономии и физике солнца. М. : Наука, 1993. Вып. 100. С. 168-188.
5. Использование JI4M ионозоида в адаптивной системе KB радиосвязи / В.А. Иванов, Н.М. Богута, Ю.В. Нога, С.А. Терехов, В.П. Урядов, В.В. Шумаев // Радиотехника. 1993. № 4. С. 77-79.
6. Особенности распространения KB сигналов на среднеширотных трассах в условиях геомагнитных возмущений / В.П. Урядов, В.И. Куркин, Г.Г. Вертоградов, В.Г. Вертоградов, А.А. Понятов, С.Н. Пономарчук // Известия вузов. Радиофизика. 2004. T.XLVII, №12. С. 1041-1056.
7. Погрешности прогнозирования ионосферного прохождения радиоволн на основе глобальной ионосферной модели / И.В. Крашенинников, И.Б. Егоров, О.П. Коломийцев, Ю.Н. Черкашин // Геомагнетизм и аэрономия. 2004. Т. 44, №2. С. 221-226.
8. Котович Г.В., Михайлов С.Я. Адаптация возможности модели IRI в прогнозировании характеристик декаметровых радиотрасс // Геомагнетизм и аэрономия. 2003. Т.43, №1. С.88-92.
9. Сажин В.И. Гибридное моделирование распространения декаметровых радиоволн: автореф. дис. ... д-ра физ.-мат. наук. Иркутск: ИГУ, 2001.30 с.
10. Барабашов Б.Г., Вертоградов Г.Г. Динамическая адаптивная модель связного декаметрового канала // Радиотехника. 1995. №12. С. 29-32.
11. Мальцева О.А., Родионова В.Т., Шлюпкин А.С. Использование полного электронного содержания для текущей диагностики состояния ионосферного канала // Геомагнетизм и аэрономия. 2005. Т.45, №4. С. 480486.
12. Ортиков М.Ю., Троицкий Б.В. Применение ионосферного индекса солнечной активности для расчёта параметров ионосферы при обеспечении КВ-радиосвязи // Геомагнетизм и аэрономия. 2000. Т.40, №4. С. 135-136.
13. McNamara L.F. The use of total electron density measurements to validate empirical models of the ionosphere // Adv. Space Res. 1985. Vol. 5, N 7. P. 81-90.
14. Komjathy A., Langley R.B. Improvement of a global ionospheric model to provide ionospheric range error corrections for single - frequency GPS users
// Presented at the ION 52nd Annual Meeting. Cambridge, Ma, 19-21 June 1996, 10 p.
15. Полуэмпирическая модель ионосферы для широкого диапазона геофизических условий / В.М. Поляков, В.Е. Суходольская, М.К. Ивель-ская, Г.Е. Сутырина, Г.В. Дубовская, М.Ю. Бузунова. М.: МЦЦ-Б, 1986. 136 с.
16. Сажин В.И. Компьютерное моделирование распространения радиоволн в регулярной ионосфере. Иркутск : Изд-во Иркут. гос. ун-та, 2010.
17. http://sapp.telepac.pt/coaa/beaconsee.htm
18. Афанасьев Н.Т., Грозов В.П., Тинин М.В. Эффекты ионосферных неоднородностей при наклонном распространении декаметровых радиоволн вблизи и выше МПЧ // Исслед. по геомагнетизму, аэрономии и физике Солнца. М.: Наука, 1987. Вып.78. С. 190-197.
19. Паньков JI.B., Унучков В.Е. Определение максимально-применимых частот по измерениям статистических характеристик сигналов // Тезисы докладов XIV Всесоюзной конференции по распространению радиоволн. М.: Наука, 1984. С. 161-162
20. Автоматизированный комплекс средств прогнозирования условий распространения декаметровых радиоволн на базе ионозонда с ЛЧМ / В.И. Алтынцева, И.Г. Брынько, И.А. Галкин, В.П. Грозов, Н.И. Двинских, Г.В. Котович, В.И. Куркин, С.М. Матюшонок, В.Е. Носов, С.С. Савков, В.В. Хахинов // Исследования по геомагнетизму, аэрономии и физике Солнца. М.: Наука, 1990. Вып. 92. С. 141-152.
21. Определение критических частот ионосферы по данным наклонного зондирования / Т. П. Зимнюхова, М. К. Ивельская, В.И. Сажин, Д. Ю. Свиридов, В. Е. Суходольская // Тезисы докладов Международной конференции «Физика ионосферы и атмосферы Земли». Иркутск, ИСЗФ СО РАН,
1998. С. 69-70.
22. Агарышев А.И. Метод расчёта МПЧ многоскачковых трасс // Радиотехника. 1985. №4. С. 67-70.
23. On some problems of the theory of radio wave propagation in a randomly inhomogeneous ionosphere / M.V. Tinin, N.T. Afanasyev S.M., Mikheev, A.P. Pobedina, O.V. Fridman // Radio Sei. 1992. Vol. 27, N 2. P. 245-255.
24. Полякова A.C. Исследование суточного хода "Вертикального" ПЭС в спокойных геомагнитных условиях // Труды 11 конференции молодых ученых «Гелио- и геофизические исследования». Иркутск: ИСЗФ СО РАН,
1999. С. 174-177.
25. Ma G., Maruyama Т. Construction of an empirical model relating TEC to foF2 // IEEE Xplore Digital Library, 2004
Основные результаты диссертации опубликованы в работах:
Издания, рекомендованные ВАК РФ:
AI. Голыгнн В.А., Сажнн В.И., Унучкоп В.Е. Адаптация к текущим условиям параметров ионосферного радиоканала по наблюдениям за сигналами реперных радиостанций // Геомагнетизм и Аэрономия. 2007. №47. С. 71-75.
А2.. Адаптация среднемесячной модели ионосферы к текущим условиям по данным о максимально - применимых частотах двухскачковых реперных радиолиний / В.А. Голыгин, Е.М. Вдовин, В.И. Сажин, В.Е. Унучков // Геомагнетизм и Аэрономия. 2009. Т. 49, № 3. С. 387-392.
A3. Модель наземного ионосферного канала в задачах связи и местоопределения / В.А. Голыгин, М.К. Ивельская, В.И. Сажин, В.Е. Унучков, C.B. Унучков // Геомагнетизм и Аэрономия. №49. С. 1-3. 2009
Другие издания:
A4. Аппроксимация формы профиля электронной концентрации для высот выше максимума в полуэмпирической модели ионосферы / В.А. Голыгин, М.И. Ивельская, В.И. Сажин, В.Е. Суходольская, A.B. Тащилин // Труды 6 сессии молодых ученых "Волновые процессы в проблеме космической погоды". Иркутск: ИСЗФ СО РАН, 2003. С. 72-74.
А5. Адаптация модели ионосферы к текущим условиям по критическим частотам без привлечения специализированных измерений / В.А. Голыгин, В.П. Грозов, В.А. Сажин, В.Е. Унучков // Труды 6 сессии молодых ученых "Волновые процессы в проблеме космической погоды". Иркутск: ИСЗФ СО РАН, 2003. С. 69-71.
А6. Коррекция критических частот в модели ионосферы по данным о полном электронном содержании / В.А. Голыгин, В.П. Грозов, О.С. Лесю-та, В.И. Сажин // Труды 6 сессии молодых ученых "Волновые процессы в проблеме космической погоды", Иркутск, ИСЗФ СО РАН, 2003г, с.56-57.
А7. Региональный контроль ионосферной обстановки без организации специализированных измерений / В.А. Голыгин, В.П. Грозов, В.И. Сажин, В.Е. Унучков // Труды 7 конференции молодых ученых "Взаимодействие полей и излучения с веществом", Иркутск: ИСЗФ СО РАН, 2004. С. 70-71.
А8. Коррекция формы профиля электронной концентрации ионосферы в области выше максимума / Б.А. Балаганский, В.А. Голыгин, М.К. Ивельская, А.О. Рупасов, В.И. Сажин, В.Е. Суходольская, A.B. Тащилин // Труды 7 конференции молодых ученых "Взаимодействие полей и излучения с веществом". Иркутск: ИСЗФ СО РАН, 2004. С. 117-118.
А9. Прогнозирование ионосферных условий для радиотрасс в глобальной сети / В.А. Голыгин, М.К. Ивельская, В.И. Сажин, П.П. Шутов // Труды 7 конференции молодых ученых "Взаимодействие полей и излучения с веществом". Иркутск: ИСЗФ СО РАН, 2004. С. 140-142.
А10. Голыгин В.А., Кинсар С.А., Сажин В.И. Уточнение критических частот ионосферы по данным о максимально-наблюдаемых частотах двух-
скачкового распространения декаметровых радиоволн // Труды 8 конференции молодых ученых «Астрофизика и физика околоземного космического пространства». Иркутск: ИСЗФ СО РАН, 2005. С. 66-68.
All. Голыгин В.А., Сажин В.И., Унучков В.Е. Оценка точности определения максимально-применимых частот радиолиний по наблюдениям за сигналами реперных радиостанций // Труды 9 конференции молодых ученых «Физические процессы в космосе и околоземной среде». Иркутск: ИСЗФ СО РАН, 2006. С. 256-258.
А12. Голыгин В.А., Сажин В.И., Унучков В.Е. Коррекция модели ионосферы по данным о максимально - применимых частотах реперных радиолиний // Электронный журнал "Исследовано в России". 255. С. 2463-2474. 2006 http://zhurnal.ape.relarn.ru/articles/2006/255.pdf
А13. Оценка интегральных возможностей модели ионосферы / Е.М. Вдовин, В.А. Голыгин, К.С. Паламарчук, В.И. Сажин // Труды 10 конференции молодых ученых «Современные проблемы в астрофизике и физике космической плазмы». Иркутск: ИСЗФ СО РАН,2007. С. 99-101. 2007.
А14. Адаптация параметров модели наземного ионосферного радиоканала к текущим условиям / Е.М. Вдовин, В.А. Голыгин, В.И. Сажин, В.Е. Унучков, C.B. Унучков // Труды 11 конференции молодых ученых «Гелио- и геофизические исследования». Иркутск: ИСЗФ СО РАН, 2009. С. 63-64.
AI5. Коррекция параметров модели трансионосферного радиоканала / Е.М. Вдовин, В.А. Голыгин, М.К. Ивельская, В.И. Сажин., O.A. Ситникова // Труды 11 конференции молодых ученых «Гелио- и геофизические исследования». Иркутск: ИСЗФ СО РАН, 2009. С. 65-67.
А16. Модель фоновой ионосферы с возможностью коррекции на текущую ситуацию в нескольких высотных областях / Е.М. Вдовин, В.А. Голыгин, М.К. Ивельская, В.И. Сажин // Труды 12 конференции молодых ученых "Взаимодействие полей и излучения с веществом". Иркутск: ИСЗФ СО РАН, 2011. С. 37.
А17. Голыгин В.А., Вдовин Е.М., Сажин В.И., Уточнение модельных значений критических частот ионосферы в текущей ситуации по величине полного электронного содержания // Труды 12 конференции молодых ученых "Взаимодействие полей и излучения с веществом". Иркутск: ИСЗФ СО РАН. с. 40. 2011
Научное издание
Голыгин Виктор Александрович
АДАПТАЦИЯ ПАРАМЕТРОВ МОДЕЛИ ИОНОСФЕРНОГО РАДИОКАНАЛА К ТЕКУЩИМ УСЛОВИЯМ БЕЗ ОРГАНИЗАЦИИ СПЕЦИАЛИЗИРОВАННЫХ ИЗМЕРЕНИЙ
АВТОРЕФЕРАТ
Подписано в печать 10.04,2012. Формат 60x90 1/16 Усл. печ. л. 1,5. Тираж 100 экз. Заказ 26
Издательство ИГУ 664003, Иркутск, бульвар Гагарина, 36; тел. (3952) 24-14-36
Введение. Актуальность темы, постановка задачи, цель и задачи работы.
Глава 1. Исходная модель среднего состояния ионосферного декаметрового радиоканала.
1.1 Получение системы лучевых уравнений в рамках геометро-оптического приближения.
1,20перативный расчёт максимально-применимой частоты (МПЧ) односкачкового и двухскачкового распространения.
1.33адание среднемесячных ионосферных условий для радиоканала.
1,4Разработка и реализация компьютерной программы модели радиоканала.
Глава 2. Уточнение модельных значений критической частоты ионосферы в текущих условиях по наблюдениям за прохождением сигнала на реперных радиолиниях.
2.1 Фиксация значений МПЧ в моменты радиовосхода и радиозахода.
2.2 Установка для фиксации МПЧ и оценка точности определения её величины.
2.3 Разработка способа уточнения критической частоты по наблюдениям за сигналами на односкачковых радиолиниях.
2.4 Разработка способа уточнения критической частоты по наблюдениям за сигналами на двухскачковых радиолиниях.
Глава 3. Уточнение модельных значений критической частоты ионосферы в текущих условиях по значениям полного электронного содержания ионосферы (ПЭС).
3.1 Развитие полуэмпирической модели ионосферы ИГУ (ПЭМИ) для описания верхнего участка.
3.2 Разработка способа уточнения критической частоты по значениям
Глава 4. Адаптация параметров модели радиоканала к текущим условиям в задачах связи и местоопределения.
4.1 Адаптация параметров модели радиоканала к текущим условиям в задачах связи.
4.2 Адаптация параметров модели радиоканала к текущим условиям в задачах местоопределения.
Ионосферное распространение декаметровых радиоволн по-прежнему привлекает внимание исследователей. Резервная связь, зондирование удалённых участков земной и морской поверхностей, местоопределение положения передатчиков, исследования ионосферы- эти основные задача определяют сферу его использования. На такое распространение влияет множество факторов, значение которых в конкретной ситуации определить чрезвычайно трудно. Поэтому важной и актуальной задачей является адаптация, хотя бы частичная, параметров рассматриваемого радиоканала к текущей ситуации.
В настоящее время разработаны некоторые способы частичной адаптации параметров радиоканала за счёт применения определённых технических устройств. Так, для повышения уровня принимаемого сигнала используются ретрансляторы и приём в двух разнесенных по пространству точках [1]. Применяются также устройства, реализующие автоматический выбор частоты из заданного набора частот [2], системы с адаптацией антенных устройств к некоторым особенностям принимаемого сигнала [3]. Развивается и другая группа методов, использующая модель параметров радиоканала. В одном из методов этого направления осуществляется прямой перенос характеристик модели радиоканала, определяемых по измерениям на опорных радиолиниях, на близкорасположенную трассу. Данный метод обладает хорошей точностью в небольшой окрестности опорных линий, т.к. "переносятся" значения непосредственно характеристик сигнала, измеренные в данный момент времени. Однако метод требует развертывания достаточно широкой сети опорных линий с проведением на них сложных специализированных измерений. В России развернута такая сеть станций, оборудованных ЛЧМ зондами, и ведутся работы групп исследователей в Институте солнечной и земной физики СОР АН (ИСЗФ), Марийском государственном университете, Научно-исследовательском радиофизическом институте (НИРФИ), Ростовском государственном университете (РГУ) [4]
6]. Естественно, развитие и расширение такой сети требует существенных затрат.
В другом, менее затратном методе этого направления адаптации, используется модель радиоканала, входной частью которой является среднемесячная модель ионосферы, а выходной- блок расчёта характеристик сигнала. Т.к. среднемесячные вариации основных параметров ионосферы описываются в настоящее время достаточно удовлетворительно, то задачей адаптации является, в данном случае, уточнение значений параметров модели на конкретную ситуацию, обеспечивающее повышение точности соответствия рассчитываемых характеристик сигнала измеряемым. Т.к. станции, непосредственно измеряющие параметры ионосферы, имеются лишь в небольшом числе пунктов земной поверхности, то разрабатываются также способы получения оперативной информации об ионосфере на основе приближённых решений обратной задачи распространения радиоволн на наклонных трассах. Здесь можно указать работы, проводимые в Институте земного магнетизма, ионосферы и распространения радиоволн (ИЗМИРАН)
7], ИСЗФ [8], Иркутском государственном университете (ИГУ) [9], РГУ [10], [11], Санкт-Петербургском государственном университете (СПБГУ) [12] и работы, выполненные в некоторых зарубежных учреждениях [13], [14]. Точность такого подхода (назовём его модельным) уступает точности метода прямого моделирования, однако адаптация может быть выполнена для более протяженного пространственного региона, примыкающего к опорной трассе, за счёт возможностей уточнённой модели ионосферы.
В то же время трудность широкого применения этого метода связана с тем, что для получения корректирующей информации требуется организация на опорной трассе специализированных измерений некоторых характеристик сигнала при работе передатчика на ряде частот.
В настоящей работе предлагаются способы адаптации, которые не требуют для своего применения организации специализированных измерений, что сделает адаптацию более доступной широкому кругу пользователей. В первом способе используется наличие в декаметровом диапазоне большого числа произвольных радиостанций. Это и вещательные станции, и станции точного времени, и радиолюбительские маяки, и некоторые другие станции. Они работают на определённом наборе частот в режиме непрерывного сигнала и известно их местоположение. Можно ли, ограничиваясь этой информацией о станциях, получить сведения, позволяющие выполнить адаптацию параметров модели радиоканала? Если да, то какова будет эффективность этой адаптации? На эти основные вопросы и отвечает данная работа.
Кроме того, сейчас на земной поверхности развернута большая сеть специализированных приёмников сигналов глобальной спутниковой радионавигационной системы GPS, для которых определяются и размещаются для свободного доступа в сети Internet данные по полному электронному содержанию ионосферы (ПЭС)- количеству электронов в объёме вертикального параллелепипеда с площадью основания в 1 квадратный метр. Посредством интерполяции строятся карты для любого пункта земной поверхности. Данные появляются с небольшой задержкой по времени, на которую вполне можно выполнить экстраполяцию. Решая приближенно обратную задачу, опять-таки можно попытаться получить информацию для адаптации параметров модели радиоканала. Это возможное направление также рассматривается в настоящей работе. Наконец, на основе предлагаемых способов, показывается возможность адаптации параметров модели радиоканала к текущей ситуации в некоторых задачах связи и местоопределения. В соответствии с вышеприведённым и были определены цель работы и задачи исследования.
Целями настоящей работы являются исследование некоторых вопросов адаптации параметров модели ионосферного радиоканала к текущим условиям на основе коррекции входной её части- блока задания среды распространения сигнала и разработка способов её выполнения по оперативной информации, для получения которой не требуется организации специализированных измерений.
Для достижения данной цели необходимо было решить следующие задачи:
1. Разработать способ получения корректирующей информации, для получения которой не требовалось бы организации специализированных измерений
2. Разработать способ адаптации параметров модели канала к текущей ситуации с использованием этой корректирующей информации
3. Разработать способ уточнения входных параметров модели радиоканала в текущей ситуации по величине полного электронного содержания ионосферы
4. Выполнить оценку эффективности разработанных способов для различных трасс и условий
5. Предложить способы адаптации входных параметров модели ионосферного радиоканала к текущим условиям в практических задач связи и местоопределения
Диссертационная работа состоит из введения, четырех глав, заключения и списка цитируемой литературы.
Заключение.
Подытожим основные результаты, полученные в работе.
Выполнена реализация исходной модели радиоканала, описывающей средние вариации дальности и времени распространения, углов прихода сигнала, максимально-применимых частот ионосферных радиолиний для двумерно-неоднородной изотропной регулярной ионосферы. Модель имеет хорошую применимость для описания средних значений характеристик сигнала при ионосферном распространении и поэтому может быть эффективна при использовании в задаче адаптации её параметров к текущим условиям.
Предложен способ фиксации МПЧ реперной радиолинии в моменты радиовосхода и радиозахода на ней по наблюдению временного поведения амплитуды сигнала.
В рамках численного моделирования определена связь вариаций параметров ионосферы и некоторых выходных характеристик канала- максимально-применимых частот радиолинии, углов прихода сигнала в вертикальной плоскости, времени распространения, а также связь вариаций параметров максимума электронной концентрации с вариациями полного электронного содержания ионосферы.
Разработан способ адаптации параметров модели канала к текущей ситуации по наблюдениям за сигналами реперных радиостанций на односкачковых и двухскачковых радиолиниях, в котором для получения корректирующей информации не требуется организации и проведения специализированных измерений
Разработан способ уточнения текущих значений в протяженном пространственном регионе по данным о полном электронном содержании ионосферы. Способ приближенно учитывает средние вариации формы профиля N(h) в области максимума и на участке больших высот выше максимума ионосферы.
Выполнена оценка эффективности разработанных способов, показавшая их хорошую применимость для спокойных ионосферных условий. Для метода, использующего реперные трассы, применимость показана при наличии на них слабых и умеренных продольных градиентов электронной концентрации
Предложены способы адаптации входных параметров модели ионосферного радиоканала к текущим условиям в задачах связи и местоопределения, существенно повышающие энергетический потенциал линии и точность определения углов места, модовой структуры сигнала и дальности его распространения.
В заключение, выражаю искреннюю и глубокую признательность:
- научному руководителю диссертации Сажину В.И.,
- соавторам отдельных работ Унучкову В.Е., Ивельской М.К., Вдовину Е.М., Унучкову C.B.
- за полезное обсуждение и замечания по работе Тинину М.В., Иванову В.Б., Афанасьеву Н.Т., Суходольской В.Е.,
- руководству и коллективам физического факультета ИГУ, отдела физики околоземного космического пространства ИСЗФ СО РАН.
1. Головин О.В. Декаметровая радиосвязь. М. : Радио и связь, 1990. 240 с.
2. Кузьмин Б.И. Адаптивные и автоматизированные системы связи. М. : Знание, 1994. 64 с.
3. Robert A. Monzingo, Thomas W. Miller // Introduction to Adaptive Arrays. Sci-Tech Publishing Inc, 2004.
4. Метод оперативной диагностики KB радиоканала / В.И. Куркин, В.Е. Носов, С.Н. Пономарчук, С.С. Савков, J1.B. Чистякова//Исследования по геомагнетизму, аэрономии и физике солнца. М. : Наука, 1993. Вып. 100. С. 168- 188.
5. Использование ЛЧМ ионозонда в адаптивной системе KB радиосвязи / В.А. Иванов, Н.М. Богута, Ю.В. Нога, С.А. Терехов, В.П. Урядов, В.В. Шумаев // Радиотехника. 1993. №4. С. 77-79.
6. Погрешности прогнозирования ионосферного прохождения радиоволн на основе глобальной ионосферной модели / И.В. Крашенинников, И.Б. Егоров, О.П. Коломийцев, Ю.Н. Черкашин // Геомагнетизм и аэрономия. 2004. Т. 44, №2. С. 221-226.
7. Котович Г.В., Михайлов С.Я. Адаптация возможности модели IRI в прогнозировании характеристик декаметровых радиотрасс // Геомагнетизм и аэрономия. 2003. Т.43, №1. С. 88-92.
8. Сажин В.И. Гибридное моделирование распространения декаметровых радиоволн. Автореферат диссертации на соискание ученой степени доктора физико-математических наук. Иркутск: ИГУ, 2001. 30 с.
9. Барабашов Б.Г., Вертоградов Г.Г. Динамическая адаптивная модель связного декаметрового канала // Радиотехника. 1995. №12. С. 29-32.
10. Мальцева О.А., Родионова В.Т., Шлюпкин А.С. Использование полного электронного содержания для текущей диагностики состояния ионосферного канала // Геомагнетизм и аэрономия. 2005. Т.45, N4. С. 480-486.
11. Ортиков М.Ю., Троицкий Б.В. Применение ионосферного индекса солнечной активности для расчёта параметров ионосферы при обеспечении КВ-радиосвязи // Геомагнетизм и аэрономия. 2000. Т.40, №4. С. 135-136.
12. McNamara L.F. The use of total electron density measurements to validate empirical models of the ionosphere // Adv. Space Res. 1985. Vol. 5, N7. P. 81-90.
13. Komjathy A., Langley R.B. Improvement of a global ionospheric model to provide ionospheric range error corrections for single frequency GPS users. Presented at the ION 52nd Annual Meeting. Cambridge, Ma, 19-21 June 1996, 10 p.
14. Кравцов Ю.А., Орлов Ю.И. Геометрическая оптика неоднородных сред. М. : Наука, 1980. 304 с.
15. Керблай Т.С., Ковалевская Е.М. О траекториях коротких радиоволн в ионосфере. М. : Наука, 1974. 165 с.
16. Егоров И.Б. О влиянии неоднородности и анизотропии ионосферы на характеристики распространения декаметровых волн. Автореферат дисс. кандидата физ.-мат. наук . М. : ИЗМИР АН, 1978. 45 с.
17. Унучков В.Е. Исследование угловых характеристик на трассах значительной протяженности. Автореферат дисс. кандидата физ.-мат. Наук. Горький, Изд-во ГГУ, 1979. 16 с.
18. Агафонников Ю.М., Лянной Б.Е., Лобачевская B.C. Влияние ионосферы на флуктуации азимутальных углов прихода радиоволн // Исследования распространения коротких радиоволн. М. : Наука, 1973. С. 47-56.
19. Бочаров В.И., Кольцов В.В., Лобачевский Л.А. Исследование распространения коротких радиоволн с боковым отклонением от дуги большого круга // Геомагнетизм и аэрономия. 1977. Т. 17. С. 50-56.
20. Наклонное зондирование ионосферы на радиотрассах различной протяженности / В.Б. Смирнов, В.М. Дриацкий, Ч.Л. Ходжа-Ахмедов, И.Н. Костерин // Труды ААНИИ. 1975. Т.32. С. 62-75.
21. Лукашкин В.М., Смирнов В.Б., Широчков A.B. Использование наклонного зондирования для изучения условий распространения радиоволн и ионосферного прогнозирования (обзор) // Труды ААНИИ, Ленинград. 1983. №390. С. 6-22.
22. Корка Н., Muller H.G. MUF Calculations Including the Earth's Magnetic Field // Radio Sei. 3 (New Series). 1968. №1. P. 53-56c.
23. Курант Р. Уравнения с частными производными. М. : Наука, 1964. 298 с.
24. Haselgrove J. Ray Theory and a New Method Ray-Tracing // Proc. Camb. Conf. Phys. Ionosphere (The Physical Society, Lond.). 1955. P. 355-364.
25. Казанцев A.H., Лукин Д.С., Спиридонов Ю.Г. Метод исследования распространения радиоволн в неоднородной магнитоактивной ионосфере // Космические исследования. 1967. Т. 5, №4. С. 593-600.
26. Терещенко В.Д. К вопросу о влиянии горизонтальных градиентов электромагнитной концентрации на величину МПЧ и траекторию распространения волн в ионосфере // Морфология и физика полярной ионосферы. Ленинград, Наука. 1971. С. 228-235.
27. Сажин В.И., Тинин М.В. Расчёт на ЭВМ траекторий и расходимости лучей близких к критическому для двумерно-неоднородной ионосферы // Исследования по геомагнетизму, аэрономии и физике Солнца. М. : Наука, 1975. Вып. 35. С. 83-85.
28. Сажин В.И. Моделирование на ЭВМ распространения радиоволн в регулярной ионосфере. Учебное пособие. Иркутск, Издательство ИГУ. 1993. 38 с.
29. Агарышев А.И. Алгоритм расчёта МПЧ односкачковых трасс // Исследования по геомагнетизму, аэрономии и физике Солнца. М. : Наука, 1981. Вып. 55. С. 168-172.
30. Калиткин Н.И. Численные методы. М. : Наука, 1978. 432 с.
31. Агарышев А.И. Метод расчёта МПЧ многоскачковых трасс // Радиотехника. 1985. №4. С. 67-70.
32. Полуэмпирическая модель ионосферы для широкого диапазона геофизических условий / В.М. Поляков, В.Е. Суходольская, М.К. Ивельская, Т.Е. Сутырина, Г.В. Дубовская, М.Ю. Бузунова. М. : МЦД-Б. 1986. 136 с.37. http://clust 1.wdcb.ru/spidr/index.jsp
33. Оценка применимости некоторых моделей ионосферы для расчёта углов прихода декаметровых радиоволн / А.И. Агарышев, Г.В. Дубовская, М.К.
34. Ивельская, C.B. Лопаткин, В.И. Сажин // Техника средств связи, серия «Системы связи». М. : Изд-во «Экое». 1982. Вып. 4. С. 3-10.
35. Михалевич Ю.И., Омельченко O.K. Процедуры кусочно-полиноминальных интерполяций функций одной и двух переменных. // Новосибирск. 1970. 64 с.
36. Кияновский М.П., Сажин В.И. К аналитическому представлению ионосферных данных при расчетах декаметровых радиоволн // Исследования по геомагнетизму, аэрономии и физике Солнца. М. : Наука, 1980. Вып. 57. С. 41-49.
37. Croft T.A., Hoogasin H. Extract ray calculations in a quasiparabolic ionosphere with no magnetic field // Radio Science. 1968. Vol. 3, №1. P. 6974.
38. Паньков Л.В., Унучков B.E. Определение максимально-применимых частот по измерениям статистических характеристик сигналов // Тез. докл. XIV Всесоюзной конференции по распространению радиоволн. М. : Наука, 1984. С. 161-162.
39. Голыгин В.А., Сажин В.И., Унучков В.Е. Коррекция модели ионосферы по данным о максимально применимых частотах реперных радиолиний // Электронный журнал "Исследовано в России". 255. 2006. С. 2463-2474. http://zhurnal.ape.relarn.ru/articles/2006/255.pdf
40. Афанасьев Н.Т., Грозов В.П., Тинин М.В. Эффекты ионосферных неоднородностей при наклонном распространении декаметровых радиоволн вблизи и выше МПЧ // Исследования по геомагнетизму, аэрономии и физике Солнца. М. : Наука, 1987. Вып. 78. С. 190-197.
41. Кочержевский Г.И., Ерохин Т.А., Козырев Н.Д. Антенно-Фидерные устройства. М. : Радио и Связь. 1999. 352 с.46. http://sapp.telepac.pt/coaa/beaconsee.htm
42. Troster J.G., Fabry R.S. The NCDXF/IARU International beacon network // QST. 1994. Part 1,№Ю. P. 31-37; Part 2, №11. P. 49-51.
43. Карякин В.И., Литвинцев О.Г., Сажин В.И. База данных радиозондирования ионосферы: реализация и некоторые применения // Радиофизика и электроника проблемы науки и обучения. Иркутск: Изд-воИГУ. 1995. С. 178-182.
44. Агарышев А.И., Дашеев H.H. Методы решения обратных задач наклонного зондирования в условиях горизонтально-неоднородной ионосферы // Исследования по геомагнетизму, аэрономии и физике Солнца. М . : Наука, 1988. №85. С. 155-166.
45. Диогенова Т.В., Ивельская М.К., Суходольская В.Е. Исследования пространственной корреляции критической частоты слоя F2 // Исследования по геомагнетизму, аэрономии и физике Солнца. М . : 1985. Вып. 71. С. 74-78.
46. Адаптивная модель ионосферы / И.А. Кринберг, В.И. Выборов, В.В. Кошелев, В.В. Попов, Н.А. Сутырин. М . : Наука, 1986. 132 с.57. ftp://cddis.gsfc.nasa.gov/pub/gps/data/daily/
47. Афраймович Э.Л., Перевалова Н.П. GPS-мониторинг верхней атмосферы Земли. Иркутск: ГУ НЦ ВСНЦ СО РАМН. 2006. 480 с.
48. Кринберг И.А., Тащилин А.В. Ионосфера и плазмосфера // М . : Наука, 1984. 188 с.
49. Tashchilin A.V., Romanova Е.В. Numerical modeling the high-latitude ionosphere // Proc. of COSPAR Colloquia Series. 2002. Vol. 1. P. 315-325.62. http://iszf.irk.ru/proiects.php
50. Оценка интегральных возможностей модели ионосферы / Е.М. Вдовин, В.А. Голыгин, К.С. Паламарчук, В.И. Сажин // Труды 10 конференции молодых ученых «Современные проблемы в астрофизике и физике космической плазмы». Иркутск: ИСЗФ СО РАН. 2007. С. 99-101.
51. Полякова A.C. Исследование суточного хода "Вертикального" ПЭС в спокойных геомагнитных условиях // Труды 11 конференции молодых ученых «Гелио- и геофизические исследования». Иркутск: ИСЗФ СО РАН. 2009. С. 174-177.
52. G. Ma and T. Maruyama. Construction of an empirical model relating TEC to foF2 // IEEE Xplore Digital Library, 2004
53. Агарышев А.И., Унучков В.Е. Диапазоны изменений углов места для различных способов распространения декаметровых радиоволн // M . : Наука, 1987. Исследования по геомагнетизму, аэрономии и физике солнца. Вып. 77. С. 60-65.
54. Гончаренко И.А. Компьютерное моделирование антенн. Всё о программе "MMANA" // M . : ИП Радиософт, ж. "Радио". 2002. 80 с.
55. Модель наземного ионосферного канала в задачах связи и местоопределения / В.А. Голыгин, М.К. Ивельская, В.И. Сажин, В.Е. Унучков, C.B. Унучков // Геомагнетизм и Аэрономия. 2009. Т. 49, №1. С. 1-3.
56. Расчёт статистических характеристик траекторных параметров KB сигнала на основе развития ПЭМИ / М.Ю. Бузунова, М.К. Ивельская, C.B. Лопаткин, X. Мяликгулиев, М.Д. Пежемская, В.И. Сажин, В.Е.
57. Суходольская // Исследования по геомагнетизму, аэрономии и физике Солнца. М . : Наука, 1991. Вып. 92. С. 46-50.
58. Сажин В.И., Унучков В.Е. Оценка эффективности коррекции параметров ионосферной модели по характеристикам углового спектра декаметрового сигнала // Радиофизика и электроника проблемы науки и обучения. Иркутск: изд-во ИГУ. 1995. С. 172-178.
59. Агарышев А.И., Алиев П.М, Труднев К.И. Системы коротковолновой радиосвязи с подавлением многолучёвости сигнала. Иркутск: Изд-во ИРГТУ. 2009. 160 с.
60. Антошкин Б.Н., Солоканов Е.Ю., Унучков В.Е. Установка для исследований угловой и пространственной структуры поля коротких радиоволн // Исследования по геомагнетизму, аэрономии и физике Солнца. М . : 1987. Вып. 76. С. 77-89.
61. Девис К. Радиоволны в ионосфере. М . : 1973. 447 с.
62. Голыгин В.А., Сажин В.И., Унучков В.Е. Адаптация к текущим условиям параметров ионосферного радиоканала по наблюдениям за сигналами реперных радиостанций // Геомагнетизм и Аэрономия. 2007. Т. 47, №1. С. 71-75.
63. Адаптация среднемесячной модели ионосферы к текущим условиям по данным о максимально применимых частотах двухскачковых реперных радиолиний / В.А. Голыгин, Е.М. Вдовин, В.И. Сажин, В.Е. Унучков // Геомагнетизм и Аэрономия. 2009. Т. 49, № 3. С. 387-392.
64. Региональный контроль ионосферной обстановки без организации специализированных измерений / В.А. Голыгин, В.П. Грозов, В.И. Сажин, В.Е. Унучков // Труды 7 конференции молодых ученых
65. Взаимодействие полей и излучения с веществом". Иркутск: ИСЗФ СО РАН. 2004. С. 70-71.
66. Прогнозирование ионосферных условий для радиотрасс в глобальной сети / В.А. Голыгин, М.К. Ивельская, В.И. Сажин, П.П. Шутов // Труды 7 конференции молодых ученых "Взаимодействие полей и излучения с веществом". Иркутск: ИСЗФ СО РАН. 2004. С. 140-142.
67. McNamara L.F. Ionospheric modeling in support of single station location of long range transmitters.// JATP. 1988. Vol. 50, №9. P. 781-795.
68. Коррекция параметров модели трансионосферного радиоканала / Е.М. Вдовин, В.А. Голыгин, М.К. Ивельская, В.И. Сажин, O.A. Ситникова // Труды 11 конференции молодых ученых «Гелио- и геофизические исследования». Иркутск: ИСЗФ СО РАН. 2009. С. 65-67.
69. Модель фоновой ионосферы с возможностью коррекции на текущую ситуацию в нескольких высотных областях / Е.М. Вдовин, В.А. Голыгин, М.К. Ивельская, В.И. Сажин // Труды 12 конференции молодых ученых
70. Взаимодействие полей и излучения с веществом". Иркутск: ИСЗФ СО РАН. 2011. С. 37.
71. Куркин В.И., Орлов И.И., Попов В.Н. Метод нормальных волн в проблеме коротковолновой радиосвязи. М . : Наука, 1981. 122 с.
72. Сажин В.И. Компьютерное моделирование распространения радиоволн в регулярной ионосфере: учебное пособие. Иркутск: Изд-во ИГУ. 2010. 77 с.
73. Керблай Т.С., Ковалевская Е.М. О траекториях коротких радиоволн в ионосфере. М . : Наука, 1974. 165 с.
74. Колесник А.Г., Колесник С.А., Нагорский П.М. Влияние ионосферной плазмы на распределение уровня космического радиоизлучения КВ-диапазона // Радиофизика и электроника: проблемы науки и обучения. Иркутск. 1995. С. 3-7.
75. Высокопотенциальный УКВ радар ИСЗФ СО РАН / Г.А. Жеребцов, A.B. Заворин, В.Е. Носов, А.П. Потехин // Тезисы докладов Российской научно-технической конференции по дифракции и распространению волн. Улан-Удэ: БНЦ СО РАН. 1996. С. 183-186.
76. Денисенко П.Ф. Энергетические потери декаметровых радиоволн в среднеширотной ионосфере. Автореферат дисс. на соискание учёной степени д.ф.-м.н. М . : ИЗМИР АН, 1990. 20 с.
77. Отчёт о НИР Совершенствование методов прогнозирования распространения декаметровых радиоволн. Чита: Читинский политехнический институт. 1989. №01890029725.
78. Черенкова Е.Л., Чернышев О.В. Распространение радиоволн. М . : Радио и связь. 1984. 271 с.
79. Глобальная эмпирическая модель распределения электронной концентрации, температуры и эффективной частоты соударений электронов в ионосфере / Ю.К. Часовитин, A.B. Широчков, A.C. Беспрозванная, Т.Л. Гуляева, П.Ф. Денисенко, O.A. Армейская, С.Е.
80. Иванова, А.И. Каширин, Н.М. Клюев, Е.А. Корякина, JI.C. Миронова, Т.Н. Сыкилинда, В.Б. Шушкова, В.И. Водолазкин, В.В. Соцкий, Н.Е. Шедаков // Ионосферные исслед. М . :МГК. 1988. №44. С. 6-13.
81. Брюннели Б.Е., Намгаладзе А.А. Физика ионосферы, М . : Наука, 1988. 528 с.
82. Афанасьев Н.Т., Тинин М.В., Унучков В.Е. О пространственной структуре волнового поля вблизи границы зоны тени // Геомагнетизм и аэрономия. 1998. №6. С. 835-841.
83. Stig Syndergaard. A new algorithm for retrieving GPS radio occupation total electron content // Geophysical Research Letters. 2002. Vol. 29, NO. 16, 1808, 10.1029/2001GL014478.
84. Raymond O., Davies K., Musman S. Comparisons of ionospheric total electron contents made at Boulder, Colorado, using GPS // Radio Sci. 1997. Vol. 32, №3., P. 1081-1089.
85. Суроткин B.A., Клименко B.B., Кореньков Ю.Н. Модельные расчёты полного электронного содержания для системы спутников GPS // Геомагнетизм и аэрономия. 2011. Т. 51, №5. С. 695-701.