CC BY
68
УДК 527:6287
МЕТОДЫ УТОЧНЕНИЯ ИНТЕГРИРОВАННОГО НАЗЕМНО-КОСМИЧЕСКОГО МОНИТОРИНГА В СИСТЕМЕ ГЛОНАСС
(Работа выполнена при поддержке РФФИ, грант № 11-07-00046-а)
И.П. Болодурина, д.т.н. (Оренбургский государственный университет, prmat@mail.osu.ru); В.Н. Решетников, д.ф.-м..н. (МАТИ- Российский государственный технологический университет им.. К.Э. Циолковского, г. Москва, kt-mati@mail.ru); М.Г. Таспаева (Оренбургский государственный университет, povtas@iribox.ru)
Рассматривается проблема повышения точности местоопределения наземного подвижного объекта в системе ГЛОНАСС. Предлагается метод фильтрации данных, основанный на применении многослойных персептронов, для прогнозирования навигационных параметров на текущем шаге по измерениям предыдущих шагов.
Ключевые слова: система ГЛОНАСС, точность местоопределения подвижного объекта, система мониторинга «АвтоГРАФ», нейросетевой подход, комплексная обработка информации.
Система ГЛОНАСС представляет второе поколение спутниковых навигационных систем, основным назначением которых является глобальная оперативная навигация приземных подвижных объектов: наземных (сухопутных, морских, воздушных) и низкоорбитальных космических. Термин «глобальная оперативная навигация» означает, что подвижный объект, оснащенный навигационной аппаратурой потребителей, может в любом месте приземного пространства в любой момент времени определить параметры своего движения -три координаты и три составляющие вектора скорости [1].
В последние годы спутниковые навигационные системы стали неотъемлемой частью транспортной инфраструктуры Российской Федерации. Наибольшее распространение они получили в секторе частного транспорта. Пользовательский сегмент включает в себя оборудование, позволяющее определять местоположение, скорость движения, расстояние до объекта следования и время, за которое можно достигнуть цели при определенной скорости. Наряду с ОТО-навигаторами все большее распространение получает навигационная аппаратура потребителей (НАП) ГЛОНАСС и ГЛОНАССОТО (мультисистемный приемник).
Основная задача НАП - прием информации со спутников, ее интерпретация и вывод на дисплей либо в канал связи в надлежащем виде. Так, для этих целей специалистами ООО «ТехноКом» (г. Челябинск) разработана система «АвтоГРАФ» спутникового мониторинга, которая позволяет проконтролировать в режиме реального времени местонахождение транспортного средства (ТС), направление его движения, пробег и объем затрачиваемого горючего. Для этого на борту каждого ТС достаточно установить прибор, принимающий сигналы навигационных спутников, а на диспетчерском пункте - соответствующую программу. Устройство, размещенное на борту ТС, накапливает информацию о своем местоположении с заданным периодом либо, анализируя характер
движения ТС, записывает маршрут движения, что значительно экономит объем передаваемых данных: время хранения маршрута в приборе может доходить до нескольких лет.
На рисунке 1 изображена функциональная схема контроллера системы «АвтоГРАФ».
Модуль ГЛОНАСС с помощью внешней активной антенны принимает кодовые сигналы со спутников системы ГЛОНАСС и, используя внутренний вычислитель, определяет географические координаты местоположения приемника, а также точное время, скорость и направление движения. Полученные данные по протоколу ЫЫЕА поступают с выхода модуля в блок центрального процессора. Центральный процессор - это ядро контроллера «АвтоГРАФ», связывающее все компоненты системы воедино и обеспечивающее их взаимодействие по заложенной в него программе. В качестве процессора выступает быстродействующая однокристалльная микроЭВМ, обеспечивающая скорость и точность вычислений, достаточную для решения различных навигационных и сервисных задач.
Антенна ГЛОНАСС V
Цифровой вход 1 /
Цифровой вход 2
Цифровой вход 3
Цифровой вход 4
Аналоговый вход 1
Аналоговый вход 2
* \
Модуль ГЛОНАСС
Центральный процессор
память
Драйвер питания
Выход 1
Выход 2
Рис. 1
Существующая НАП определяет координаты объекта с погрешностью 0,1-30 м и его скорость -0,02-0,2 м/с в зависимости от режима работы аппаратуры. Показания могут колебаться под влиянием внешних факторов:
- прохождение ионосферного и тропосферного слоев, влияющих на скорость прохождения сигнала;
- наличие электромагнитных помех;
- инфраструктура местности (отражение, затенение сигнала космических навигационных систем, туннели), создающая многолучевость принимаемых сигналов.
Задачей комплексной обработки навигационной информации является совместная обработка данных навигационного счисления для определения основных навигационных параметров движущегося объекта с максимально возможной точностью, которая зависит от качества навигационных измерителей (датчиков навигационной информации) и алгоритмов обработки навигационных сигналов.
Объединение (комплексирование) устройств и систем обработки информации в единый функционально, структурно и конструктивно взаимосвязанный навигационный комплекс позволяет полнее использовать имеющуюся на борту ТС избыточность информации, благодаря чему появляется возможность повышения точности, помехоустойчивости, непрерывности и надежности навигационных определений.
В настоящей статье предлагается новый подход при комплексировании на уровне вторичной обработки информации спутниковой радионавигационной системы ГЛОНАСС.
Широкое распространение получило разделение обработки аппаратных измерений в навигационных системах на первичную и вторичную. Под первичной обработкой информации подразумевается наблюдение, которое заключается в оценке вектора состояния регулируемой системы. Разного рода возмущения и помехи на входе системы в этой задаче не учитываются. Под вторичной - выполняемая в ЭВМ обработка выходных данных самих измерителей, результаты которой используются для определения и уточнения координат и скорости движения, углов ориентации подвижного объекта и источников погрешностей измерителей [2].
Сигналы, передаваемые каждым космическим аппаратом (КА) системы ГЛОНАСС в составе оперативной информации, описывают положение фазового центра передающей антенны данного КА в связанной с Землей геоцентрической системе координат ПЗ-90.02, определяемой следующим образом:
- начало координат расположено в центре масс Земли;
- ось Z направлена на условный полюс Земли, как определено в рекомендации Международной службы вращения Земли (ШКБ);
- ось X направлена в точку пересечения плоскости экватора и начального (нулевого) меридиана, установленного Международным бюро времени (В/Я);
- ось Y дополняет геоцентрическую прямоугольную систему координат до правой.
В этой системе координат положение точки в пространстве определяется значениями координат X, Y, Z. Система «АвтоГРАФ» выдает значения координат и проекций скорости ТС в дискретные моменты времени через заданный интервал, рассчитывая их по псевдодальностям до КА. Псевдодальности рассчитываются по временным задержкам Т сигнала по трассе ьй КА-потребитель и известной скорости распространения радиоволн [3].
На рисунке 2 изображен график движения ТС на стоянке. Кривая отображает зависимость от времени I (секунды) средней для координат x, у среднеквадратической ошибки определения координат S (в метрах). Здесь наблюдается пробег из-за дрейфа координат, возникающего при переотражении сигналов.
и с
Рис. 2
Запишем формальную дискретную модель движения:
x(k+1)=A(k)x(k)+V(k), ^Ц,, (1)
y(k)=C(k)x(k)+w(k). (2)
Недоступный непосредственному наблюдению фазовый вектор х(к) имеет размерность п, измеряемый вектор у(к) - размерность т. Параметры системы А(к) и С(к) представляют собой известные матрицы соответствующей размерности, а вектор возмущений V(k) размерности п и ошибка измерений w(k) размерности т - векторные случайные процессы с дискретным временем.
В РЛС автоматического сопровождения и обработки сигналов в инерциальных навигационных системах наиболее распространенным алгоритмом вторичной обработки является фильтр Калмана.
В алгоритме типа многомерной калмановской фильтрации начальными условиями на каждом новом цикле служат оценка состояния системы и величина, характеризующая ее погрешность. Данный алгоритм последовательно обрабатывает за-
ново поступающие векторы измерений, учитывая при этом значения, вычисленные на предшествующем цикле. Эта черта отличает способ фильтра Калмана от нерекуррентных алгоритмов, которым ради работы требуется оберегать целый массив обрабатываемых данных. На следующем шаге с помощью обрабатываемых на данном цикле измерений уточняются начальные условия. По мере последовательной обработки новых измерений происходит накопление фильтром полезной информации. Если элементы вектора состояния твердо выражаются сквозь измеренные величины, суммарная погрешность оценок, как правило, должна снижаться [4].
Реализация алгоритма обработки измерений может основываться на использовании в рамках информационно-управляющей системы, работающей в режиме реального времени. В силу того, что в данной системе требуется в реальном масштабе времени отслеживать значительное число быстроменяющихся параметров, вырабатывать эффективные и своевременные решения по управлению, необходим математический аппарат, позволяющий обрабатывать большой поток сложноструктурированной информации о текущем состоянии производственной системы. Базой такого математического аппарата являются алгоритмы глобального распараллеливания вычислительных процессов, в частности, модели нейронных сетей. В реальных системах часто невозможно определить значения некоторого набора параметров или определить их с высокой долей погрешности, возможен также выход из строя тех или иных датчиков. Кроме того, отсутствие точной математической модели является предпосылкой для применения в системе аппарата нечетких множеств.
Основным принципом нейросетевого моделирования является принцип черного ящика. В этом случае моделируется внешнее функционирование системы, а не ее внутренняя структура, как при аналитическом подходе. Функционирование системы в рамках нейросетевого моделирования описывается информационно на основе наблюдений или данных экспериментов над реальной системой. Интерпретация получаемых с помощью таких моделей результатов более сложна, чем при использовании формальных математических моделей, однако несомненное достоинство - отсутствие ограничений на сложность моделируемых систем - определяет их важную практическую значимость.
В данной работе предложен подход к определению местоположения подвижного объекта в системе ГЛОНАСС, использующий нейронные сети.
В процессе функционирования многослойная нейронная сеть формирует выходной сигнал Х1 в соответствии с входными сигналами Х0 по формуле
X, = ^ (• Х,_! ) =
= ^ (™,Т_1 • (^Т_2 •... • ^ ^0Т • Х„ ) ...)) ,
где 1 - число слоев в сети; XI - выход нейронов 1-го слоя; Х0 - вход сети; - матрица весовых коэффициентов между нейронами (1-1) и 1 слоев; ^ - активационная функция нейронов 1-го слоя,
1 = 01, X = Х(1), у = х.
Выходной сигнал - это результат обработки входного сигнала при его прохождении по связям сети. Следовательно, нейронная сеть реализует функциональное соответствие между входом и выходом и может служить информационной моделью Р* (1) системы Р (2).
Предложенный алгоритм реализуется следующим образом. На вход сети подаются измерения в навигационной системе. В выходном слое два нейрона: один отвечает за формирование оценки плотности распределения вероятности входного сигнала при условии перехода из состояния Sj в состояние Si, а второй оценивает вероятность противоположного события. Скрытый слой содержит радиальные базисные функции с центрами в точках из обучающей выборки. На рисунке 3 представлена архитектура предиктора.
ГгЛ уО-42 (1) чЛО""*^ «,13 (1)
\г
Гг \ п<и
А г7 У/О"23 (^ Л/О—► Аф12 (1) Аф1з (1) ►А<ф23 (1)
Входной слой Скрытый слой Выходной слой
Рис. 3
Измерения подаются на входной слой сети и хранятся в элементах кратковременной памяти 5Г1. Сети для разных состояний системы имеют одинаковую структуру.
На вход каждой сети подается вектор измерений предыдущего шага у(1-1)={Аа2,12(1-1), Да2дз(1-1), Аа2,2э(1-1), Афи(1-1), Аф1з(^1), Аф23^-1)}, на выходе получается оценка измерений текущего шага у(1) = {Аа112(1), Аа113(1),
Аа «23(1), Аф 12 (1), Аф 13 (1), Аф 23 (1)} для каждого состояния системы Si, 1=0, 1, 2, 3. При обучении использованы алгоритм обратного распространения ошибки и функция тренировки.
Для ограничения пространства поиска при обучении ставится задача минимизации целевой функции ошибки нейронной сети, которая находится по методу наименьших квадратов:
1 р
E(w) = - 2 (yj - dj) :
2 j=l
(4)
где yj - значение ]-го выхода нейросети; dj - целевое значение j-го выхода; р - число нейронов в выходном слое.
Обучение нейросети осуществляется методом градиентного спуска, то есть на каждой итерации изменение веса производится по формуле
Aw = —Л'
8E öwij
(5)
где ^ - параметр, определяющий скорость обучения.
Введем вспомогательную переменную:
5(„) = ' dyj.
(6)
Теперь можно определить рекурсивную формулу для определения 5(п) п-го слоя, если известно 5("+1) следующего п+1-го слоя:
5(n) =
w(n+-)
dS;
(7)
Нахождение 5(п) для последнего слоя нейронной сети не представляет трудности, так как известен целевой вектор, то есть вектор тех значений, которые должна выдавать нейронная сеть при данном наборе входных значений:
Запишем формулу (5) в раскрытом виде:
5(N) = ( y(N) — d )•dyL. j vJi ч HS
(8)
Awj =—Л'5('
(n) .
(9)
Рассмотрим полный алгоритм обучения нейронной сети.
1. Подать на вход нейросети один из требуемых образов и определить значения выходов нейронов нейронной сети.
2. Рассчитать 5(К) для выходного слоя нейронной сети по формуле (8), а изменения весов Aw[jN) выходного слоя N - по формуле (9).
3. Рассчитать по формулам (7) и (9) 5(11) и Awjjn) соответственно для остальных слоев нейронной сети, п=^1, ..., 1.
4. Скорректировать все веса нейронной сети: ,,(п) I*\ — „,(п) о- _ 1 л А„Дп)
w(n)(t) = w(n)(t — 1) + Aw(n)(t).
(10)
5. Если ошибка существенна, повторить алгоритм с шага 1.
На этапе 2 сети поочередно в случайном порядке предъявляются векторы из обучающей последовательности.
Для повышения точности местоопределения наземного подвижного объекта предложен метод контроля и диагностики информационных нарушений в навигационных системах, использующий нейросетевые технологии. Перед традиционными вычислительными системами нейронные сети имеют ряд преимуществ [5].
• Используя способность обучения на множестве примеров, нейронная сеть может решать задачи, в которых неизвестны закономерности развития ситуации и зависимости между входными и выходными данными.
• Имея возможность работы при наличии большого числа неинформативных, шумовых входных сигналов, нет необходимости делать их предварительный отсев - нейронная сеть сама определит малопригодность для решения задачи и отбросит их.
• Нейронные сети потенциально отказоустойчивы: то есть при неблагоприятных условиях их производительность падает незначительно. Например, если поврежден какой-то нейрон или его связи, извлечение находящейся в памяти информации затрудняется. Однако, принимая во внимание распределенный характер хранения информации в нейронной сети, можно утверждать, что только серьезные повреждения структуры нейронной сети существенно повлияют на ее работоспособность. Поэтому снижение качества работы нейронной сети происходит медленно.
• Нейронные сети обладают потенциальным сверхвысоким быстродействием за счет использования массового параллелизма обработки информации в системах, которые сейчас активно развиваются.
В заключение необходимо отметить, что одной из важнейших задач глобальной навигации в системе ГЛОНАСС является точное определение пространственных координат наземных подвижных объектов. Эта задача должна решаться комплексно. При первичной обработке оценивается вектор состояния наблюдаемой системы. Вторичная обработка информации предполагает уточнение координат и скорости движения объекта.
Следует ожидать, что предложенный метод фильтрации данных в навигационных системах, основанный на применении многослойных пер-септронов для прогнозирования навигационных параметров на текущем шаге: по измерениям на предыдущих шагах может повысить эффективность управления наземным подвижным объектом, использующим систему ГЛОНАСС.
Литература
1. Решетников В.Н. Космические телекоммуникации (Начала). Тверь: Изд-во ЗАО НИИ ЦПС, 2009. 128 с.
2. Сучилин В.И., Волобуев Г.Б. Оценка возможностей повышения точности местоопределения наземного подвижного объекта путем вторичной обработки показаний аппаратуры
пользователя систем GPS NAVSTAR и/или ГЛОНАСС // Кибернетика и высокие технологии XXI века: сб. докл. VIII Меж-дунар. науч.-технич. конф. Воронеж, 2007. Т. 2. С. 1066-1073.
3. Глобальная навигационная спутниковая система ГЛОНАСС: Навигационный радиосигнал в диапазонах L1, L2 с открытым доступом и частотным разделением. [Интерфейсный контрольный документ]. М., 2008.
4. Браммер К., Зиффлинг Г. Фильтр Калмана-Бьюси. Детерминированное наблюдение и стохастическая фильтрация; [пер. с нем.]. М.: Наука. Глав. ред. Физматлит, 1982. 200 с.
5. Деева А.С. Математическое и алгоритмическое обеспечение диагностики информационных нарушений инерциаль-ных навигационных систем на основе нейросетевого подхода: дис. ... канд. техн. наук. Челябинск: Изд-во ЮУрГУ, 2010.
УДК 519.8
МЕТОД ПЛАНИРОВАНИЯ ИННОВАЦИОННОЙ ДЕЯТЕЛЬНОСТИ
А.С. Гейда, к.т.н.; И.В. Лысенко, д.т.н..; А.И. Птушкин, к.т.н.; Е.В. Седлов
(Санкт-Петербургский институт информатики и автоматизации РАН, geida@iias.spb.su, ilys@iias.spb.su, russell9@yandex.ru)
Рассматривается метод решения задачи планирования инновационной деятельности. На основе особенностей решаемой задачи сделан вывод о том, что она должна решаться как оптимизационная. Для ее представления в виде математической оптимизационной задачи введено понятие «инновационный потенциал», предложены показатель его оценивания, вербальная и формальная постановки задачи оптимального планирования инновационной деятельности. Предложен способ решения задачи, основанный на модификации метода случайного поиска.
Ключевые слова: новация, инновация, модернизация, инновационный потенциал, эффективность функционирования, планирование, оптимизация.
Инновационный и модернизационный пути развития экономики делают актуальным решение ряда задач, связанных с инновационной деятельностью (ИД). Одной из них является задача планирования ИД, решение которой - план ИД - это список инновационных модернизационных мероприятий (ИММ), реализуемых тем или иным способом. Инновационной назовем деятельность по созданию и реализации новации. При этом под новацией будем понимать новое знание, которое может быть внедрено в практику и использовано, а под инновацией - новацию, внедренную в практику. В работе [1] описаны особенности планируемой инновационной деятельности и показано, что задачу следует ставить и решать как оптимизационную задачу распределения ресурсов на графовых моделях. Для ее решения были введены новые показатели оценивания качества ИД на основе инновационного потенциала [2] модернизируемой производственной системы (МПрС). А именно, под целевым результатом ИД в задаче понимается улучшение эффектов модернизируемой производственной деятельности (МПД). МПД протекает в виде функционирования МПрС для достижения ею возможных целей. Улучшение эффектов МПД оценивается по показателям, характеризующим спектр процессов, соответствующих разным возможным целям функционирования МПрС. Качество функционирования МПрС оценивают по эффективности ее функционирования, то есть по свойству целенаправленного процесса функционирования системы (МПрС), характеризующему приспособленность функционирования системы к достижению заданной цели [3].
Свойство системы (МПрС), позволяющее использовать ее для достижения любой из возможных целей, назовем потенциалом системы (МПрС). Инновационный потенциал системы - потенциал системы, получаемый в результате ИД. Задача планирования ИД состоит в следующем.
Пусть заданы: ИММ, которые могут быть включены в план ИД; возможные способы реализации ИММ и прогнозируемые результаты ИММ (возможные затраты денежных средств и времени на производственные мероприятия) для каждого способа реализации; графики МПД; ресурсы, которые могут быть израсходованы при инновационной и производственной деятельности.
Требуется найти оптимальный план ИД. При этом под оптимальным следует понимать такой план ИД, который позволит добиться лучшего значения показателя качества инновационной и производственной деятельности в процессе достижения возможных целей МПрС. Этот показатель - мера возможности того, что после завершения ИД будут получены требуемые целевые результаты функционирования МПрС и при этом не перерасходуются ресурсы, предоставленные для реализации ИД и МПД.
Введем следующие обозначения: М1 - 1-е ИММ; M -< ^ : 1 = М >; М^р - это М1 при р1-м
способе его реализации, pi е 0,Р;: р1=0, если М;>р не реализуется; Б - {р :1 = 1,1}; М1Р -- {М,,р : Р1 = Ор); © г {м,>ч : 1 = 1,1}; ^ -
вектор характеристик способа реализации М;>р ;
