Аппаратура потребителей спутниковых радионавигационных систем (дисциплина) — различия между версиями

Материал из SRNS
Перейти к: навигация, поиск
(Курсовая работа)
(Этап 1. Использование сторонних средств)
Строка 420: Строка 420:
  
  
=== Этап 1. Использование сторонних средств ===
+
=== Этап 1. Обработка логов навигационного приемника ===
  
'''Дедлайн''': 03.03.21 (устранение замечаний и итоговый мёрдж до 09.03.21)
+
'''Дедлайн''': 04.03.21 (устранение замечаний и итоговый мёрдж до 13.03.21)
  
Конечная цель всего курсового проекта - получить библиотеку функций на Си++, позволяющую рассчитывать положение спутника Beidou по его эфемеридам.  
+
Конечная цель всего курсового проекта - получить библиотеку функций на Си++, позволяющую рассчитывать положение спутника GPS по данным с демодулятора его сигнала L1 C/A.  
На первом этапе подготовим вспомогательные данные для разработки: эфемериды и оценки положения спутника от сторонних сервисов (чтобы было с чем сравниваться на след. этапах)
+
На первом этапе реализуем модуль разбора навигационного сообщения до структуры эфемерид, сравним результаты со сторонней программой.
  
На крыше корпуса Е МЭИ установлена трехдиапазонная антенна [https://en.harxon.com/u_file/product/18_08_08/Harxon%20HX-CSX601A%20Brochure.pdf Harxon HX-CSX601A].  
+
В неизвестной локации установлен навигационный приемник, принимающий сигналы GPS L1C/A и логирующий результаты этого приема в формате NVS BINR [https://navis.ru/downloads/SOFT/NV08C_BINR_Protocol_Specification_V13_ENG_10082012.pdf NVS BINR].  
Она через 50-метровый кабель, [[Blog:Mikhaylova/25.06.2018_Антенный_пост_на_8_приемников|сплиттер, bias-tee и усилитель]] подключена к трем навигационным приемникам:
+
Собранный на пятиминутном интервале файл приложен в архиве под именем BINR.bin, см. таблицу вариантов.  
* Javad Lexon LGDD,
+
Файл содержит наблюдения псевдодальностей и прочих радионавигационных параметров, демодулированные и разобранные данные навигационного сообщения.  
* SwiftNavigation Piksi Multi,
+
* Clonicus разработки ЛНС МЭИ.  
+
  
Эти приемники осуществляют первичную обработку сигналов Beidou B1I, выдавая по интерфейсам соответствующие потоки данных - наблюдения псевдодальностей и эфемериды спутников.
+
Для удобства студентов  данные демодулятора продублированы в текстовый файл in.txt. Каждая строка файла содержит данные одного сабфрейма одного навигационного сигнала в формате:
Данные от приемника Clonicus, записанные вечером 16 февраля 2021 года, доступны в рабочем репозитории (директория logs) в нескольких форматах.  
+
<source lang="bash">
 +
1 0 013 0R GpsL1CA  # 13 212130404  29 125  53 100010111010101010101010001010010111000011001010010111111000101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100
 +
</source>
 +
где ''# 13'' - номер спутника
  
Во-первых, это дамп бинарного потока данных от приемника в формате [https://navis.ru/downloads/SOFT/NV08C_BINR_Protocol_Specification_V13_ENG_10082012.pdf NVS BINR].
+
Требуется:
 +
# Обработать
  
Во-вторых, текстовый файл данных пакета 0xF7, полученный из данного дампа - для тех, кому лень разбираться с бинарным файлом и протоколом.
 
Он получен подобным printf'ом для каждого спутника с периодом передачи эфемерид:
 
<source lang=C>
 
fprintf(fid, "%2u %9.3f %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %22.17e %u %u %u %u %u %u\n",
 
        pack->SatNum,      // 1
 
        pack->bds[0].toe,  // 2
 
        pack->bds[0].Crs,  // 3
 
        pack->bds[0].Dn,    // 4
 
        pack->bds[0].M0,    // 5
 
        pack->bds[0].Cuc,  // 6
 
        pack->bds[0].e,    // 7
 
        pack->bds[0].Cus,  // 8
 
        pack->bds[0].sqrtA, // 9
 
        pack->bds[0].Cic,  // 10
 
        pack->bds[0].Omega0,// 11
 
        pack->bds[0].Cis,  // 12
 
        pack->bds[0].i0,    // 13
 
        pack->bds[0].Crc,  // 14
 
        pack->bds[0].omega, // 15
 
        pack->bds[0].OmegaDot, // 16
 
        pack->bds[0].iDot,  // 17
 
        pack->bds[0].Tgd,  // 18
 
        pack->bds[0].toc,  // 19
 
        pack->bds[0].af2,  // 20
 
        pack->bds[0].af1,  // 21
 
        pack->bds[0].af0,  // 22
 
        pack->bds[0].URA,  // 23
 
        pack->bds[0].IODE,  // 24
 
        pack->bds[0].IODC,  // 25
 
        pack->bds[0].codeL2,// 26
 
        pack->bds[0].L2P,  // 27
 
        pack->bds[0].WN    // 28
 
        );
 
</source>
 
  
 
Для приемки первого этапа оформляется отчет по [https://mpei.ru/Science/ScienceDocuments/Documents/Polojenie_2018/0_GOST_7.32-2017.pdf ГОСТ 7.32], содержащий результаты работ:
 
Для приемки первого этапа оформляется отчет по [https://mpei.ru/Science/ScienceDocuments/Documents/Polojenie_2018/0_GOST_7.32-2017.pdf ГОСТ 7.32], содержащий результаты работ:

Версия 23:03, 15 февраля 2022

Дисциплина "Аппаратура потребителей спутниковых радионавигационных систем" входит в вариативную часть общенаучного цикла занятий учебного плана подготовки специалистов по направлению 11.05.01 "Радиоэлектронные системы и комплексы".

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

По дисциплине подготовлен электронный учебно-методический комплекс.

Содержание

Календарный план

Неделя Даты Лекция Лабораторная Контрольное мероприятие
2 07.09 Основы СРНС. Сегменты СРНС. Методы позиционирования. Ошибки позиционирования в городских условиях
3 13.09.21-19.09.21 Навигационные сигналы.
4 20.09.21-26.09.21 Антенна, МШУ КМ-1: Модель сигнала ГНСС
5 27.09.21-03.10.21 Радиочастотный блок
(Аналогово-цифровое преобразование. Синтезатор частот. План частот АП СРНС. )
6 04.10.21-10.10.21 Математическая модель коррелятора
(Статистический эквивалент коррелятора)
КМ-2: Ошибки позиционирования в городских условиях
7 11.10.21-17.10.21 Аппаратная реализация коррелятора
(Принципы построения и структура коррелятора. Цифровой генератор гармонического сигнала и дальномерного кода.)
8 18.10.21-24.10.21 Следящие системы
(Система слежения за фазой сигнала. Структурная схема ССФ, дискриминаторы, следящий фильтр.
Система слежения за частотой. Структурная схема ССЧ, дискриминаторы, следящий фильтр.
Система слежения за задержкой. Структурная схема ССЗ, дискриминаторы, следящий фильтр.
Система слежения за задержкой с поддержкой от ФАП. Комплексный фильтр.
)
GNSS Planning КМ-3: Статистический эквивалент коррелятора
9 25.10.21-31.10.21 Оценка SNR, детекторы захвата, детекторы CC
10 01.11.21-07.11.21 Обнаружение сигнала.
(Обнаружение сигнала. Алгоритмы, статистические характеристики обнаружения.)
11 08.11.21-14.11.21 Поиск сигнала
(Поиск сигнала по частоте и задержке. Алгоритмы поиска, быстрый поиск сигнала.
Блок быстрого поиска)
12 15.11.21-21.11.21 Символьная и оверлейная синхронизация
Исследование коррелятора АП СРНС с помощью имитационной модели КМ-4: Исследование коррелятора АП СРНС ГЛОНАСС с помощью имитационной модели
13 22.11.21-28.11.21 Демодуляция навигационных данных
14 29.11.21-05.12.21 Кадровая синхронизация, декодирование, парсинг КМ-5: Прогнозирование спутникового созвездия
15 08.12.21-12.12.21 Формирование измерений, расчет положения спутника
16 13.12.21-19.12.21 Модель наблюдений псевдодальности.
Вторичная обработка измерений. Одношаговый алгоритм решения навигационной задачи.
Контроль целостности
Решение навигационной задачи по псевдодальномерным измерениям КМ-6: Решение навигационной задачи по псевдодальномерным измерениям
17 20.12.21-26.12.21 Помехоустойчивость АП СРНС
18 27.12.21-02.01.21 Зачет

Журнал успеваемости.

Правила аттестации

Обязательным условием для получения допуска к экзамену является выполнение всех контрольных мероприятий: сдача двух контрольных и защита 4 лабораторных работ.

При написании отчетов следует опираться на правила оформления ГОСТ 7.32 (в части форматирования, оформления рисунков, таблиц, списков и т.д.).

Приемка осуществляется через репозиторий на github на базе проекта gnssrecbasics2021.

Демо модели

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

Производится обработка сигнала Gps L1 C/A, сформированного лабораторным имитатором сигналов, синхронизированным с приемником. Сигнал перенесен с помощью РЧБ NT1065 на промежуточную частоту f_{if} = 8.54 МГц и преобразован многоразрядным АЦП с частотой сэмплирования F_d = 99.375 МГц. Установленная на имитаторе мощность обеспечивает уровень отношения сигнал/шум 27 + N дБГц, где N - номер вариант по журналу.

В приемнике производтся корреляционная обработка сигнала, интервал накопления составляет T = 1 мс.

Сигнал имеет сдвиг частоты несущей f_d = 100 N Гц, при этом опорный сигнал коррелятора имеет аналогичный сдвиг f_{d, rep} = f_d. Задержка огибающей известна точно и используется в опорном сигнале \tau = \tau_{rep}.

Задание: Построить на плоскости I,Q вид множества точек корреляционных сумм, полученных в ходе большого (сотни) числа экспериментов.


ЧАП и ФАП на стат эквивалентах


Вопросы к экзамену

Экзамен проводится в письменной форме: тест, задача и один экзаменационный вопрос (экзаменационные вопросы). При подготовке рекомендуется использовать книги по тематике, конспекты лекций и презентации.

Курсовая работа

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

Курсовая работа по АП СРНС выполняется во втором семестре. В рамках КР студенты знакомятся с рядом инструментов и техник, используемых при разработке навигационных приемников.

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

Требования к разрабатываемому программному модулю:

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

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

  • разработка модуля разбора символов навигационного сообщения
  • расчет положения КА в Matlab/Python и его проверка сторонними сервисами;
  • реализация модуля расчета положения КА на С/С++ и его тестирование.

На каждом из этапов действуют следующие правила:

  • Взаимодействие осуществляется через github (пул реквесты, комментарии)
  • Этап сдан тогда, когда принят пул реквест


Варианты заданий:
Студент Номер спутника (варианта) Логи Примечание
Антропов Евгений Александрович 1 tar.gz
Бахолдин Никита Владимирович 10 tar.gz
Берестнев Владимир Сергеевич 13 tar.gz
Дубинин Леонид Андреевич 15 tar.gz
Еловский Никита Станиславович 17 tar.gz
Иванцова Дарья Николаевна 19 tar.gz
Капитонов Александр Ильич 24 tar.gz
Коробков Андрей Юрьевич 28 tar.gz
Куксин Владислав Александрович 30 tar.gz
Кусакин Павел Александрович
Ливенцов Василий Александрович
Михайлов Иван Олегович
Муратов Николай Сергеевич
Попов Максим Геннадьевич
Солнцева Милана Константиновна
Танкина Анна Михайловна
Цымбал Григорий Романович
Чепелев Игорь Игоревич


Этап 1. Обработка логов навигационного приемника

Дедлайн: 04.03.21 (устранение замечаний и итоговый мёрдж до 13.03.21)

Конечная цель всего курсового проекта - получить библиотеку функций на Си++, позволяющую рассчитывать положение спутника GPS по данным с демодулятора его сигнала L1 C/A. На первом этапе реализуем модуль разбора навигационного сообщения до структуры эфемерид, сравним результаты со сторонней программой.

В неизвестной локации установлен навигационный приемник, принимающий сигналы GPS L1C/A и логирующий результаты этого приема в формате NVS BINR NVS BINR. Собранный на пятиминутном интервале файл приложен в архиве под именем BINR.bin, см. таблицу вариантов. Файл содержит наблюдения псевдодальностей и прочих радионавигационных параметров, демодулированные и разобранные данные навигационного сообщения.

Для удобства студентов данные демодулятора продублированы в текстовый файл in.txt. Каждая строка файла содержит данные одного сабфрейма одного навигационного сигнала в формате:

1 0 013 0R GpsL1CA  # 13         212130404  29 125  53 100010111010101010101010001010010111000011001010010111111000101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100101010101010101010101010111100

где # 13 - номер спутника

Требуется:

  1. Обработать


Для приемки первого этапа оформляется отчет по ГОСТ 7.32, содержащий результаты работ:

  1. Определение формы орбиты и положения спутника на ней на начало рассматриваемого интервала времени по данным сервиса CelesTrak: общий вид + положение спутника на 18:00 МСК 16 февраля 2021, так, чтобы было видно подспутниковую точку и время (примерно так)
  2. Расчет графика угла места собственного спутника от времени по данным Trimble GNSS Planning Online на интервал времени с 18:00 МСК 16 февраля до 06:00 МСК 17 февраля 2021 года.
  3. Расчет диаграммы угла места и азимута спутника (SkyView, он же SkyPlot) по данным Trimble GNSS Planning Online на на интервал времени с 18:00 МСК 16 февраля до 06:00 МСК 17 февраля 2021 года.
  4. Формирование списка и описание параметров, входящих в состав эфемерид в сигнале B1I Beidou
  5. Формирование таблицы эфемерид собственного спутника (номер спутника Beidou см. в журнале выше) с подписанными размерностями (как тут, но только для одного спутника)

Работа осуществляется в репозитории на github на базе bdsephexercise. Для этого:

  1. Завести, если до сих пор нет, аккаунт на github
  2. Сделать форк проекта bdsephexercise (маленькая кнопка сверху справа в вэб-интерфейсе)
  3. Установить на рабочем компьютере git (для выполнения последующих этапов удобнее будет сделать это в ОС GNU/Linux, вероятно, на виртуальной машине)
  4. Клонировать свой форк на рабочий компьютер (git clone https://github.com/MySuperAccount/bdsephexercise.git)
  5. Промежуточные и итоговые результаты по первому этапу хранить в каталоге processing (git pull; git add <files>; git commit -m "Сделал то-то то-то"; git push)
  6. По результатам первого этапа завести Pull Request в исходный проект в ветку, название которой совпадает с вашей фамилией (кнопка наверху в вэб-интерфейсе)
  7. Устранить все замечания, добиться принятия Pull Request'а.

Этап сдан тогда, когда ваш ПР принят.

Этап 2. Моделирование

Дедлайн: 02.04.21 (исправления замечаний до 09.04.21)

Эфемериды - параметры некоторой модели движения спутника. В разных ГНСС эти модели разные, а значит отличается и формат эфемерид, и алгоритмы расчета положения спутника.

Одна из самых простых и удобных моделей - в системе GPS. Beidou наследует данную модель.

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

Построить трехмерные графики множества положений спутника Beidou с системным номером, соответствующим номеру студента по списку. Графики в двух вариантах: в СК ECEF WGS84 и соответствующей ей инерциальной СК. Положения должны соответствовать временному интервалу с 18:00 МСК 16 февраля до 06:00 МСК 17 февраля 2021 года. Допускается использовать одни и те же эфемериды на весь рассматриваемый интервал.

Построить SkyView за указанный временной интервал (напоминаю, антенна на крыше корпуса Е) и сравнить результат с Trimble GNSS Planning Online, полученный на прошлом этапе.

Оформить отчет по результатам этапа:

  1. Реализация в Matlab или Python (описание модели и её листинг)
  2. Таблица использованных эфемерид
  3. Трехмерные графики положений спутника в ECEF и ECI (не забудьте подписать оси, изобразите соответствующую Земле сферу в начале СК)
  4. Расчётный и полученный в GNSS Planing Online SkyView
  5. Сравнение результатов моделирования с результатами первого этапа

Работы по данному этапу осуществляются в каталоге simulation. Правила приемки этапа те же, что и на первом этапе.

Этап 3. Реализация

Дедлайн: 11.05.21 (исправления до 18.05.21)

Требуется разработать на языке С/С++ функцию расчета положения спутника Beidou на заданное время по шкале UTC, минимизируя время её исполнения и количество затрачиваемой оперативной памяти. Вызов функции не должен приводить к выбросу исключений или утечкам памяти при любом наборе входных данных.

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

Программный модуль должен сопровождаться unit-тестами под check:

  • Тесты функции решения уравнения Кеплера
  • Тест расчетного положения спутника в сравнении с Matlab/Python с шагом 0.1 секунды.

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

Требуется провести проверку на утечки памяти с помощью утилиты valgrind.

Оформить отчет по результатам курсового проекта. В качестве первых двух глав использовать отчёты с предыдущих этапов, в третьей главе отразить результаты этого этапа:

  1. Код реализации
  2. Вывод тестов, включая анализ времени исполнения
  3. Вывод valgrind
  4. Вывод по этапу
  5. Заключение по проекту

В репозитории расположен каркас программного модуля libbdssvpos, представленного в виде подключаемой библиотеки на базе системы сборки CMake. Настроена компиляция статической/динамической библиотек компилятором/кросскомпилятором и выполнение юнит-тестов (проверено в GNU/Linux системах). Для примера в библиотеке реализованы две функции: умножения и сложения.

В функцию сложения add внесена ошибка, на которую указывает результат соответствующего юнит-теста:

korogodin@KorPC:~/bdsephexercise/libbdssvpos$ mkdir build
korogodin@KorPC:~/bdsephexercise/libbdssvpos$ cd build
korogodin@KorPC:~/bdsephexercise/libbdssvpos/build$ cmake ..
-- The C compiler identification is GNU 5.4.0
-- The CXX compiler identification is GNU 5.4.0
-- Check for working C compiler: /usr/bin/cc
-- Check for working C compiler: /usr/bin/cc -- works
-- Detecting C compiler ABI info
-- Detecting C compiler ABI info - done
-- Detecting C compile features
-- Detecting C compile features - done
-- Check for working CXX compiler: /usr/bin/c++
-- Check for working CXX compiler: /usr/bin/c++ -- works
-- Detecting CXX compiler ABI info
-- Detecting CXX compiler ABI info - done
-- Detecting CXX compile features
-- Detecting CXX compile features - done
-- libbdssvpos version: cd76a7e
-- Found PkgConfig: /usr/bin/pkg-config (found version "0.29.1")
-- Checking for one of the modules 'check'
-- Configuring done
-- Generating done
-- Build files have been written to: /home/korogodin/bdsephexercise/libbdssvpos/build
korogodin@KorPC:~/bdsephexercise/libbdssvpos/build$ make
[ 33%] Built target bdssvpos-static
[ 66%] Built target bdssvpos
Scanning dependencies of target test_bdssvpos
[ 77%] Building CXX object test/CMakeFiles/test_bdssvpos.dir/check_position.cpp.o
[ 88%] Linking CXX executable test_bdssvpos
Running unit tests
Running suite(s): My test suite
50%: Checks: 2, Failures: 1, Errors: 0
/home/korogodin/bdsephexercise/libbdssvpos/test/check_position.cpp:11:F:Main tests:test_add:0: 1+2 should be 3
test/CMakeFiles/test_bdssvpos.dir/build.make:121: ошибка выполнения рецепта для цели «test/test_bdssvpos»
make[2]: *** [test/test_bdssvpos] Ошибка 1
make[2]: *** Удаляется файл «test/test_bdssvpos»
CMakeFiles/Makefile2:179: ошибка выполнения рецепта для цели «test/CMakeFiles/test_bdssvpos.dir/all»
make[1]: *** [test/CMakeFiles/test_bdssvpos.dir/all] Ошибка 2
Makefile:127: ошибка выполнения рецепта для цели «all»
make: *** [all] Ошибка 2

Работы по данному этапу осуществляются в каталоге libbdssvpos. Правила приемки этапа те же.

Курсовая прошлых лет

Задания в 2020, в 2021.

Персональные инструменты
Пространства имён

Варианты
Действия
SRNS Wiki
Рабочие журналы
Приватный файлсервер
QNAP Сервер
Инструменты