автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.01, диссертация на тему:Методы оценивания сигналов навигационных систем на основе многоальтернативного и неполного стохастического описания

[и

I

Кошаев Дмитрий Анатольевич

МЕТОДЫ ОЦЕНИВАНИЯ СИГНАЛОВ НАВИГАЦИОННЫХ СИСТЕМ НА ОСНОВЕ МНОГОАЛЬТЕРНАТИВНОГО И НЕПОЛНОГО СТОХАСТИЧЕСКОГО ОПИСАНИЯ

Специальность 05.13.01 - «Системный анализ, управление и обработка информации»

Автореферат диссертации на соискание ученой степени доктора технических наук

Санкт-Петербург 2010

2 2 НЮЛ ?010

004607385

Работа выполнена в ГНЦ РФ ОАО «КОНЦЕРН «ЦНИИ «ЭЛЕКТРОПРИБОР»

Научный консультант

доктор технических наук, профессор С.П. Дмитриев.

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

доктор технических наук, профессор А.В. Небылов,

доктор физико-математических наук, профессор А.Е. Барабанов,

доктор технических наук, профессор Р.И. Ивановский.

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

Институт проблем управления им. В.А. Трапезниш РАН (Москва).

Защита состоится 14 сентября 2010г. в 1500 на заседании диссертационно совета ДС 411.007.01 при ОАО «КОНЦЕРН «ЦНИИ «ЭЛЕКТРОПРИБОР» г адресу: 197046, С.-Петербург, ул. Малая Посадская, 30.

С диссертацией можно ознакомиться в научно-технической библиотеке ОА «КОНЦЕРН «ЦНИИ «ЭЛЕКТРОПРИБОР».

Автореферат разослан 2.3 июня 2010г.

Ученый секретарь диссертационного совета,

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

Н.В. Колесов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы. Существует большое число навигационных приложений, в которых удобные для реализации калмановские алгоритмы неприменимы в силу неопределенности стохастических моделей. В работе рассматривается два вида неопределенности, один из которых связан с многоальтернативностью, а второй с неполнотой стохастического описания. В обоих случаях для вектора состояния и вектора измерений можно записать следующие уравнения в дискретном времени:

Ч = фкЧ-1 + + Щ, Ч = Нкхк + G( рьЭ) + vk. (1)

Здесь к - номер дискретного момента времени tk = к At; At - интервал дискретности; \>к,Щ - чисто-случайные последовательности; F(), G(-)

- известные матрицы и вектор-функции; S- вектор известных или неизвестных параметров либо случайных величин. Величина р может иметь различную природу, которая и определяет тип неопределенности. В задаче требуется оценить значения функций F,G и вектора х.

Величина р может представлять собой случайный вектор и иметь дискретное распределение вероятности, например, при наличии дискретных помех измерений, нарушений или аномальных ошибок. В этом случае речь идет о многоальтернативном стохастическом описании. Классические варианты решения задач с применением такого описания сводятся к построению банка фильтров Калмана (ФК), каждый из которых настроен на одну из возможных гипотез о значении р. Впервые такое решение в случае отсутствия переходов между гипотезами предложил D.T. Magill. Впоследствии оно было детально проработано Д.Г. Лайниотисом. Многоальтернативное решение, предусматривающее переходы от одной гипотезы к другой (Interacting Multiple-Model), впервые представил Н.А.Р. Blom.

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

Другой рассматриваемый в диссертации тип неопределенности не предусматривает наличие стохастической модели для величины р. При этом Pit = pfe), где р(<) - функция, о которой известно лишь, что она непрерывно дифференцируема и для нее справедливы ограничения |p(i)| < 1, |р(<)| < р.

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

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

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

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

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

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

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

Основные задачи

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

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

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

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

5. Разработка метода оценивания непрерывно дифференцируемого сигнала с учетом ограничений на величину сигнала и его производную. Получение адекватной характеристики точности.

6. Подтверждение эффективности предложенных методов по результатам моделирования и обработки реальных данных навигационных систем.

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

Научные положения, выносимые на защиту

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

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

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

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

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

5. Метод приближенного вычисления среднеквадратической ошибки оценивания сигнала с учетом непрерывной дифференцируемости и ограничений.

Научная новизна

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

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

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

Практическая значимость

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

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

• Спроектирован алгоритм обнаружения и оценивания нарушений в навигационном комплексе, включающем две и более ИНС. С точки зрения качества получаемого решения алгоритм, не уступает классическому варианту на основе банка ФК для отдельных нарушений, но при этом требует значительно меньше вычислительных ресурсов. Эффективность алгоритма подтверждена путем моделирования задачи диагностики для навигационного комплекса типа «Аппассионата ЭКМ.1» и обработки реальных данных от ИНС «Дельта».

• Для навигационного комплекса из трех ИНС разработан алгоритм обнаружения и компенсации аномальных ошибок в показаниях курса без использования детального описания их динамики. Результаты стендовых испытаний навигационного комплекса типа «Ладога-1143 0», проведенные с имитацией различных типов нарушений, показывают, что алгоритм позволяет без значительного запаздывания обнаруживать и компенсировать дополнительную ошибку, вызванную нарушением, с точностью, сопоставимой с уровнем номинальных ошибок.

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

Применение результатов. К настоящему времени алгоритм обнаружения и оценивания аномальной ошибки внедрен в программное обеспечение навигационного комплекса «Ладога 11430» (2007г). Алгоритм исключения неоднозначности фазовых измерений ГНСС применялся для компенсации вертикальных ускорений в показаниях аэрогравиметра «Чекан-А» (2002-2003г). Разработанные алгоритмы определения курса и относительной траектории по фазовым измерениям ГНСС были апробированы в ходе проведения испытаний

в городских условиях и на оз. Ладога (тема «БАСН-М-ЭП», 2008-2009г). Сейчас алгоритм определения курса находятся на стадии внедрения в бортовую аппаратуру ГНСС-интерферометра с применением микромеханических инерциальных датчиков.

Апробация работы. Результаты диссертации докладывались на III, VI, IX, X, XI, XIV и XVI Санкт-Петербургской международной конференции по интегрированным навигационным системам (С.-Петербург, 1996, 1999, 2002, 2003, 2004, 2007, 2009); XX, XXI, XXII, XXIII, и XXV научно-технической межотраслевой конференции памяти H.H. Острякова (С.-Петербург, 1996,1998; 2000, 2002, 2006); I и II научно-технической конференции молодых ученых "Навигация и управление движением" (С.-Петербург, 1999, 2000); научно-технической конференции «Навигация и океанография» (С.-Петербург, 2001, 2007); форуме китайско-российских молодых ученых (Пекин, Китай, 2003); I и III Всероссийской конференции «Фундаментальное и прикладное координатно-временное обеспечение» (С.-Петербург, 2005, 2009); 5th International Conference on Differential Satellite Navigation Systems (С.-Петербург, 1996); Institute of Navigation 55th Annual Meeting (Boston, USA, 1999); 6th IFAC Symposium on Fault Detection, Supervision and Safety of Technical Processes (Beijing, PR China, 2006); 13th IEEE Conference "Methods and Models in Automation and Robotics" (Szczecin, Poland, 2007); 16th IEEE "Mediterranean Conference on Control and Automation" (Ajaccio, Corsica, France, 2008), на Общероссийском семинаре «Современные методы навигации и управления движением» (Москва, Институт проблем управления им. В.А. Трапезникова РАН, 2010).

Публикации. По материалам диссертации опубликовано 36 статей и докладов, включая 15 статей в журналах, рекомендованных ВАК.

Структура и объем работы. Диссертация состоит из пяти глав, введения, заключения и списка используемой литературы из 156 наименований. Общий объем работы составляет 242 страницы, включая 45 рисунков.

СОДЕРЖАНИЕ РАБОТЫ Первая глава посвящена оцениванию сигналов при наличии дискретно распределенной помехи измерений.

В задаче требуется оценить вектор состояния

хк=фкЧ+Щ> (2)

по измерениям вида

zk=Hkxk+AN + vk, (3)

где Ф к,Нк,А - известные матрицы; wk, v^ - чисто-случайные последовательности; N - случайный /и-мерный вектор, компоненты которого

принимают целочисленные значения. Отметим, что в отношении вектора N, по существу, применяется вариант многоальтернативного стохастического описания. При этом (2), (3) связаны с уравнениями (1) из общей постановки задачи следующим образом: F = О, р = N, Э - вектор, состоящий из элементов матрицы Л, G = AN.

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

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

чисел Ът, для которого апостериорная вероятность максимальна. Этот этап -наиболее трудоемкая часть решения. Здесь особенно остро стоит проблема сокращения вычислительных затрат. Решение задачи в два этапа является распространенным приемом, основанным на интуитивных соображениях. Его теоретическое обоснование представлено в диссертации.

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

N* может быть получено в результате минимизации

N* = arg min (N-N)T -N) = aig min q(N), (4) Ne Ът Neil"

где N, Ptf - оценка и ковариационная матрица ее ошибок, полученные на первом этапе. При наличии дополнительного измерения вида

(5)

где е - центрированная гауссовская погрешность в знании длины вектора х с

2 * дисперсией ст|, обосновано следующее выражение для N

V V

N* =arg min [(N - N)1 P^1 (N-N) + (Q- \ x{N) |)2 /a^] = arg min q(N),(6) Ne Zm Ne Ът

V

где x(N) - оценка вектора x, полученная по измерениям (3) в предположении,

V- V _ V v 2 2

что N является истинным значением, x(N) x(N)/\x(N)\ +ае -

условная по отношению кг и N дисперсия Рх/ дг - ковариационная матрица

V

ошибки оценки х(Щ.

Для вычисления минимизируемой в (4) квадратичной формы «/(ТУ) используется рекуррентный способ

Яj=Яj-l+(Hj-Nj/j-\)2/c'2¡/j_l, % =0. Ч = Ят> О)

где Щ/о =N1, ауо = 0^; Nj/j-l, Оу/у_1, ] =2,т - условные по отношению — т

к вектору Nj_^=(N^,...Nj_l) математическое ожидание и дисперсия компоненты Лу, которые определяются исходя т N и Р^, Тогда как для получения функции (¡(И), минимизируемой в (6), после (7) выполняется еще одна итерация = qm + (С,- |х(ЛГ) |) .

Поиск значения (4) производиться в области д(М)<д, а значения (6) - в области ¿¡(Щ<Ц, где ц - заданный порог. Этот целочисленный поиск и является источником вычислительных проблем. Основу предлагаемых методов оценивания составляют оригинальные приемы, позволяющие более рационально по сравнению с существующими методами использовать вычислительные ресурсы при выполнении целочисленного поиска

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

(Nj/j-\), при ,гу- =1;

.('/-о

(8)

Ы; 1 + (-1) 1 -1), при sJ■> 1 и -Л у/уч) <0;

Яр _1) + (~1р -1), при *у > 1 и (Л^ 0,

где ¿у =1,2,... номер значения Nj, <•> -

»г. . \ операция округления до ближайшего

} 1 \. иъ целого числа. С помощью этой нумерации значения компонент N '-' автоматически сортируются по мере их удаления от центра интервала перебора в котором Я) йч (см. рис. 1). Поскольку qJ■

увеличивается с ростом номера \yi~~Vj aJ/j-lлjf~QJ значения И], первым рассматривается Рис. 1. Нумерация значений Nj

^Р, а затем - и т.д. до тех пор, пока д/\ ¿д. Введенная

4 4 у * ^ —¿v ^

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

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

При отыскании в ходе поиска каждого нового значения Ы, для которого й ц, порог д устанавливается равным . Это уже известный прием, но он более органично встраивается в предлагаемую схему перебора значений N на основе нумерации по правилу (8).

Достоверность значения Л^ оценивается с помощью апостериорной * * *

вероятности g =Р(ЛГ = ^/ | £) либо отношения g ко второй по величине

вероятности g . Отношение g /g вычисляется проще. Оно представляет собой верхнюю границу 1д2 наибольшей апостериорной

вероятности. В случае, если g* или

g / достигает порога, близкого

к 1, значение И* принимается за истинное, и для этого предположения вычисляется

оценка вектора состояния х. В предлагаемом методе уже по промежуточным результатам выявляется ситуация, когда соответствующая характеристика достоверности не сможет превысить заданный порог. Это позволяет избежать бесполезных уточняющих расчетов, связанных с * * . ** вычислением g и g I g .

На рис. 2 результаты использования методики для

представлены практического предложенной оценивания

относительной траектории судна по

3.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 дх, км

Рис. 2. Относительная траектория судна. Точки - положение с дискретностью 2с.

фазовым одночастотным измерениям ГНСС с частотой поступления 10Гц. В этой задаче подлежащий оцениванию вектор состояния х включает ДХ, ДУ, Ы, -прямоугольные координаты антенны бортового приемника относительно антенны опорного приемника вдоль географических осей. Опорный приемник необходим для организации дифференциального режима обработки. Он располагался на расстоянии порядка 20км от судна. Эксперимент проводился на оз. Ладога. В роли дискретно распределенных помех ДО в данной задаче, как уже отмечалось, выступают неизвестные заранее значения целого числа периодов фазовых измерений. Один период фазы соответствует 19см. Шумы фазовых измерений имеют значительно меньший уровень. Наиболее вероятное

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

Рис. 3. Результаты определения ориентации циркулирующего судна.

Рис. 3 демонстрирует результаты применения предложенной методики для оценивания курса и качки судна по фазовым одночастотным измерениям ГНСС, поступающим с частотой 10Гц. Суть этой задачи состоит в том, что на борту судна устанавливается две или более антенн (в данном эксперименте использовались две антенны, установленные в поперечном направлении на расстоянии 3,4м). Измеряя фазы спутниковых сигналов и располагая информацией о расположении антенн в связанной с объектом системе координат можно определить углы ориентации объекта относительно

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

значение вектора целого числа периодов определялось в соответствии с (6). При этом в качестве дополнительного измерения (5) выступало приближенно известное расстояние между антеннами, т.е. £=|х|+е. На верхнем графике показана реализация курса Атнсс+ммг, полученного в результате совместной обработки фазовых измерений и данных ММГ. Судно совершало циркуляцию в различных направлениях по траектории представленной на рис. 2. На втором графике показана угловая скорость АГгнсс+ммг, вычисленная по приращениям комплексной оценки курса. Третий график демонстрирует бортовую качку бгнсс, рассчитанную по фазовым измерениям. Отсутствие скачков в К шсс+ммг и Эпнсс говорит о правильном определении вектора целого числа периодов фазовых измерений N с помощью разработанной методики. Перепады по угловой скорости и по бортовой качке связаны не с ошибками, а с изменением направления циркуляции судна.

Вторая глава посвящена многоальтернативному методу обнаружения и оценивания нарушений в информационно-измерительных системах (ИИС) с помощью расширенного фильтра Калмана (РФК).

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

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

Ч = фкЧ-\ + Щ + ВкО(у)ц, гк = Нкхк + Ькщ +Ук+ СкБ(у)г\. (9)

1 X/ т

Здесь -л = (г|1,...т] ) - вектор всех возможных нарушений; -

известные матрицы; £>(у) - диагональная матрица с диагональю, образованной

из компонент вектора у = (у1,принимающих значения 0 или 1 (у' =1

передает наличие нарушения г|' в ИИК). При этом лишь одна из у' равна 1, а

остальные - 0. Априорные вероятности нарушений р' = Р(у-' = 5у, ] = 1,7/), где

5у - символ Кронекера, известны. Векторы хо, , уд., г] случайны и взаимно

независимы друг от друга и от у. Они имеют гауссовское центрированное распределение с известными ковариационными матрицами. Причем у^ образуют чисто-случайные последовательности.

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

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

= = 7 = I ) > локализовать действующее нарушение, оценить его

величину т|', а также вектор состояния х^ с учетом влияния т|'. Оценки т)' и хд. позволяют компенсировать нарушение и устранить последствия его влияния до начала компенсации, т.е. восстановить ИИС с нарушением.

Классическая схема решения такой задачи требует реализации банка N ФК, каждый из которых с учетом (9) оценивает х^ совместно с одним из нарушений

т|' в предположении, что именно оно действует в ИИК, т.е. =6у,_/ = 1,Лг.

В противовес такому трудоемкому решению предлагается применить один РФК, оценивающий вектор состояния увеличенной размерности ц г|Т)Т в предположении, что в ИИК одновременно действуют все нарушения, т.е. У'/=1 для j = l,N и, как следствие, £)(у) = Е. С точки зрения наблюдаемого процесса (см. (9)) события у7 =5у,у = 1,ЛГ и =0, где г|' -

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

й к 1 П (4Г1 (10)

где с . - априорная дисперсия г)7; ть , -оценка г] и ковариационная

V Щ

матрица ее ошибки, полученные РФК. Обоснована возможность получения

входящих в (10) (-Р*, и | с использованием экономичного

преобразования (Р^)-1 и Оценка составного вектора = (хк ,т]')т и ковариационная матрица ее ошибки определяются по результатам РФК, скорректированным для случая г)1 = 0.

Таким образом, на основе РФК без потери оптимальности вычисляются те же характеристики состояния ИИК, которые определяются с помощью банка ФК для отдельных нарушений. Преимущество предложенного метода достаточно очевидно: вместо N ФК, оценивающих (л+1)-мерный вектор, применяется один РФК, оценивающий вектор размерности л-НУ. Чем выше размерность вектора состояния п, тем более существенное сокращение вычислительных затрат дает применение РФК. При малых и более экономичным может оказаться применение не одного, а нескольких РФК, каждый из которых предполагает действие своей группы нарушений.

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

номинальном состояниии - гипотеза А0: у = 0, либо в аномальном состоянии с

одним из N нарушений - гипотезы А': yJ j = Кроме того, с

известными вероятностями допускаются переходы А0 А', г= 1,7/ и обратно.

Традиционное решение такой задачи сводится к построению банка из 37/+1 взаимодействующих друг с другом ФК для всех возможных переходов между гипотезами на двух соседних шагах решения (см. рис. 4, слева). При этом на каждом к- м шаге вычисляются апостериорные вероятности гипотез

= а также оценки векторов = (хк,т\' с соответствующими

ковариационными матрицами для i = l,N. Характеристики для полных цепочек переходов между состояниями здесь не определяются, так как их число с течением времени неограниченно увеличивается.

Следуя традиционной схеме взаимодействия между фильтрами,

Л ф

предлагается использовать РФК и ограничиться всего двумя гипотезами А , А . Гипотеза А* предполагает наличие любой комбинации нарушений. Состав

шагЫ шаг к

ФК (А°-»А°), ФК (А°->А'), г = 1>, ФК(А'-»А'), г = 1777, ФК(А'->А°),г = йу

ФК (А°-»А°), РФК (А°->А*), РФК (А*-»А*), РФК(А*->А°)

о „>

£*>£*>**» "л

Если зд и 1, то вычисляются

■У 'шах р

А* ' х*1"'

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

При с помощью (10)

рассчитываются условные по

отношению к А* вероятности отдельных нарушений

Для нарушения г = /тах с наибольшей вычисляется

безусловная апостериорная вероятность

4=Р(у/=5 ^ = =

а при необходимости и оценка с ковариационной матрицей

ее ошибки Р . Такая схему

.1* .2*

для нескольких составных гипотез А , А

Рис. 4. Классическое решение на основе банка частных ФК (слева) и предлагаемое решение на основе РФК (справа).

решения можно использовать и о нарушениях.

На рис. 5 сопоставлены результаты предлагаемого и классического методов на примере навигационного комплекса типа «Аппассионата ЭКМ.1», включающего две платформенные ИНС, демпфированные по данным

относительного лага. Помимо гипотезы А0 об отсутствии нарушений в расчет принимаются гипотезы о нарушениях в уходах гироскопов , со^, соХ2, <ау2 и

смещениях нуля акселерометров , а^, аХ2, ау2 с горизонтальными осями

чувствительности первой и второй ИНС, а также нарушение в уходе азимутального гироскопа сог одной из двух ИНС. В качестве измерений использованы производные разностей показаний горизонтальных составляющих скорости и курса от первой и второй ИНС. В этом примере число измерений т=3, число нарушений N=9, размерность вектора состояния

для гипотез А', / = 1,9 о нарушениях и=6. При использовании РФК наиболее эффективным оказывается следующее распределение нарушений по гипотезам:

.1* , 2* гипотеза я отвечает за нарушения , со^, соХ2, <йуг, п - за

ах1>ау1>ах2>ауг'>1?* -засо2.

В данном примере для второй ИНС промоделировано нарушение в уходе гироскопа с горизонтальной осью чувствительности х на 0,2град/ч при номинальном уровне ухода 0,01 град/ч. Судя по графику апостериорных вероятностей, метод на основе РФК позволяет выявить промоделированное нарушение со^ лишь с немного бблыпим запаздыванием, чем классический

метод. Необходимые для восстановления ИНС оценка нарушения аХ2 и оценки

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

состояния я=6 здесь не велика, объем вычислений для реализации метода на основе РФК в 4 раза меньше по сравнению с классическим методом.

Оценка дополнительной ошибки вертикали а№

Апостериорные вероятности гипотез

вызванной действием нарушения о^.угп.с

0,3 ОД 0,1 0

1 2 Время действия нарушения, мин

Оценка нарушения тх град/ч

^Велич дна нару пения |

АА-Ч/*

Время действия нарушения, мин

Оценка дополнительной ошибки скорости 8 УЕ, вызванной действием нарушения оХ1, м/с

12 3 4 Время действия нарушения, мин

2 3 4 5

Время действия нарушения, мин

Рис. 5. Результаты, полученные с применением РФК (черные линии) и результаты классического решения (серые линии).

Помимо моделирования проводилась также экспериментальная проверка предложенного экономичного метода с использованием записанных на стенде реальных данных ИНС «Дельта». Ее результаты подтвердили сопоставимость классического и предложенного методов как с точки зрения задержки в определении типа нарушения, так и по точности вырабатываемых оценок.

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

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

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

одной ИИС. Исходя из этого, в задаче рассматривается гипотеза й° о номинальном состоянии ИИК и N гипотез И',1 = 1,Ы о нарушениях в какой-либо ИИС.

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

имеются априорные вероятности состояний ИИК р' = Р(Л'),1 = 0,# и переходов между ними

9" = А^), / = 1,7/. Важно отметить, что в задаче

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

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

Для гипотезы А0 используются вектор состояния , описывающий номинальные ошибки ИИС, либо их разности, входящие в измерения Для

каждой из гипотез А',¡ = \,И вводится расширенный вектор состояния

1 т —/ т т ■■■/

хк =(ХЛ>Т4 ) > включающий помимо хк подвектор т^, который в

соответствии с полиномиальной моделью служит для описания аномальных

ошибок /'-й ИИС. Компонентами г\'к являются сами аномальные ошибки, а

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

Полагаясь на введенные модели номинальных и аномальных ошибок и

априорные вероятности состояний ИИК />',г = 0,# и переходов между ними

«Д = \,И, получаем типичную байесовскую задачу

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

g1k = Р(й' 1/' = и если ^ «1 для некоторого />0, то вырабатывается

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

Полиномиальная модель аномальных ошибок может быть адекватной лишь на некотором ограниченном интервале

времени Т. Для решения задачи на сколь угодно большом времени можно регулярно с периодом Т обновлять полиномиальную модель

аномальных ошибок путем перезапуска банка ФК. Но при

перезапуске апостериорные Рис. 6. Временная диаграмма работы банков ФК. вероятности и оценки не смогут

сразу достоверно отражать реальное состояние ИИС из-за переходных процессов.

Чтобы преодолеть этот недостаток предлагается использовать два параллельно действующих банка ФК, которые периодически перезапускаются с разницей во времени, равной Г/2 их непрерывного функционирования (см. рис. 6). По структуре оба банка ФК идентичны. Для принятия решения о наличии аномальной ошибки и выработки оценки, компенсирующей ее влияние, используется тот банк, который отработал более половины времени после перезапуска Это «рабочий» интервал действия банка. На начальном «подготовительном» интервале после перезапуска результаты банка не участвуют в формировании окончательного решения. Такая схема позволяет всегда иметь результаты работы одного из банков без переходных процессов, возникающих после перезапуска (за исключением начального этапа решения).

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

Банк 1

Перезапуск

Банк 2

Банк 1

Банк 2

«Рабочий» интервал «Подготовительный» интервал

.ошио*.

^ Банк 1 &

X

■Оь Банк 2

£

Банк 1

"А

% : ^

Банк 2

хс

Комбинированная реализация £ служит дня принятия решения о наличии 1 аномальной ошибки _

чччч

Рис. 7. Комбинирование результатов двух банков ФК.

Достоинства представленного метода решения

проиллюстрированы на

рис. 7. Здесь предполагается, что в некоторый момент времени появилась аномальная ошибка в /-й ИИС. В этом случае первый банк вырабатывает повышенную

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

вероятности g', которая не позволит уверенно судить о наличии аномальной ошибки в 1-й ИИС. Оценка аномальной ошибки при использовании одного банка ФК также не может быть удачно использована для компенсации из-за ее недостаточной точности

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

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

корректируемых ИНС представлялись в виде марковских процессов 1-го порядка. Для описания аномальных ошибок

применялась

полиномиальная модель 3-го порядка. Период Т перезапуска банков ФК составляет 5ч. Для эксперимента

использовались реальные данные ИНС, записанные на стенде с имитацией рысканья, качек и циркуляции.

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

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

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

Рис. 8. Результаты стендовых испытаний с имитацией нарушений:

а) линейного тренда ухода гироскопа с горизонтальной осью чувствительности 0,06град/ч2;

б) скачкообразного изменения ухода азимутального гироскопа на 0,3град/ч;

в) скачкообразных изменений ошибки курса.

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

Четвертая глава посвящена методу оценивания сигналов при неполном стохастическом описании. Это второй из рассматриваемых в диссертации типов неопределенности моделей. Предлагаемый здесь метод основан на кусочно-полиномиальном представлении нестохастического сигнала с учетом ограничений и непрерывной дифференцируемое™.

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

собой /-мерный вектор 0(О = (0(')Т,с(Ор(О)Т» где (/-1)-мерный вектор 6(0

описывается дифференциальным уравнением

ё(0=м(/)9(0+Д0р(0, е(о)=е0. (и)

Здесь функция с(г) и матрицы М(1), Ц?) известны, а начальное условие 0о -неизвестно. Функция р(() непрерывно дифференцируема Для нее и ее производной справедливы ограничения

|р(0|£1, |р(фр- (12) Допускается любая размерность 9, в том числе /=1, когда 0(<) = с(()р(().

В задаче требуется получить оценки значений в для моментов времени = к А/ (А/ - интервал дискретности) по набору /и-мерных измерений вида

Ч=СкЩ) + нкхк+Ч> (13)

где хд. - и-мерный вектор, описывающий коррелированные во времени помехи с помощью стохастической модели

С]с,Н^,Фк - известные матрицы; ^к^к " чисто-случайные последовательности; х0> ^ представляют собой взаимно независимые гауссовские векторы с нулевым средним и известными ковариационными матрицами. Отметим, что первое слагаемое в (13) представляет собой функцию б в уравнении (1) для вектора измерений, общем для всех рассматриваемых в диссертации задач. При этом в качестве аргумента & этой функции выступает неизвестное начальное условие 9д.

Для определенности рассматривается задача сглаживания 9(^) по реализации г измерений г,- как для г <к, так и для ¡>к.

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

Предполагается, что для интервалов времени j/5-^,*^], где №=sT, s=l,2,... - номер интервала, Г - его длительность, для функции p(t) справедливо полиномиальное представление

p(t) = b^ +b¡s)(t-¿s-V)/T + (b¡s) /т2 . (14)

При этом связаны с коэффициентами на предыдущем интервале

соотношениями

Ъ^ = b[s), = 2¿j(í) - 2bW - 6W. (15)

Соблюдение (15) для каждого s обеспечивает непрерывность функции p(f) и ее производной на всей временной оси, в том числе, на границах интервалов. Поскольку непрерывная дифференцируемость р(г) соблюдается независимо от

значения коэффициента bfs+l\ предполагается, что этот коэффициент не

связан с коэффициентами на предыдущих.интервалах.

Представление функции p(í) в виде (14) позволяет свести задачу

оценивания вектора 0(/), порождаемого р(t), к оценке полиномиальных коэффициентов и (в случае />1) начального условия 6q. Причем, в свете соотношений (15) на S интервалах оцениванию подлежат только 2+S независимых коэффициента 6®, б® и b^-1, s =1,5.

Важно отметить, что полиномы второго порядка являются минимально достаточными для соблюдения заданных в задаче свойств функции р(/). Так, в случае представления p(l) с помощью полиномов первого порядка (кусочно-линейной функции) р(/) имеет разрывы на границах интервалов. Полиномы же третьего и более высоких порядков избыточны, имея сразу несколько «свободных» коэффициентов, связь которых с коэффициентами на предыдущем интервале по условию задачи никак не регламентируется.

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

В качестве критерия в данной задаче используется функция правдоподобия

измерений z относительно вектора Ь^ и 8q (для

задачи с многомерным полезным сигналом Q(í)). В этом случае задача сводится к определению значений

(^>Х) = аг8 шах /(г|6®,е0). (16)

где В^ - область, включающая те и только те значения вектора при которых для р(*) и р(0 с учетом (14), (15) выполняются неравенства (12).

Показано, что значения (16) могут быть получены в два этапа. На первом этапе с помощью ФК по измерениям г с учетом уравнений (11), (13), но без

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

ЬЧ8) =аг8_{ (ь^В*0 = в0(Г®), (17)

где /^(5) - полученная на первом этапе ковариационная матрица ошибки

оценки Ь^ \ ЫЬ*{3)) - условная по отношению к Ь^ -Ь*^ оценка 0о>

рассчитанная исходя из результатов оценивания Ь ^, 90 на первом этапе.

Центральное место в решении поисковой задачи (17) занимает нетривиальная проблема определения области

в аналитической форме.

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

коэффициент. В работе получено следующее определение области В^ из

От О —/С"\

(Л'+2)-мерного множества Е действительных значений Ьу ', при которых р(0 удовлетворяет заданным условиям:

В(5) = бЕ5+2 г = + (18)

Здесь (•) - нелинейные функции вида

-8,

¿(1) °0

-1+ . 2 / 6 = рт,

2(5+КТ

2 2 Ь^ - вектор, состоящий из коэффициентов Ь^, Ь^, связанных с

Ь^Ь?^ рекуррентными соотношениями (15); Ртах(^02^>

соответственно верхняя и нижняя границы коэффициента б^,

которые определяются по правилам

Ртах (¿>02

Ртт (¿02 ) =

ЙихРйгХ при Ь2 < 2(1-Ь0)-~8;

Ртах (¿02 )> при 2(\-Ь0)-Ь<Ь2<2{\-Ь0),

Ртах (¿02)> при Ь2> 2(1-Ь0)',

Р^пФй). при Ь2>-2(\ + Ь0) + Ь-

Р^п(¿02), при - 2(1 + Ь0) <Ьг <-2(1 + Ь0) + |б;

(19)

Кйп(йй2). при Ъ2 <-2(1 + £0);

Ртах(602) = 1"^

2 - Р?шп(А02) = 60 + ' 2 .

3Г2(1 -Ь0)-Ъ2 +-7(20~ ¿о)~ ¿2 + 35)2 -852

Ртп(602) = £

Ртах (¿02) = 1"

зГ2(1 + 60) + 62 + 7(2(1 + 60) + + 35)2 - 8821 -1,

(20)

(21)

(22)

-л/Г1^

. Ртт (¿02) =

+ ЛД + 6о

-1.(23)

С помощью (18)-(23) ограничения на полиномиальные коэффициенты накладываются последовательно: сначала на потом на и т.д.

до . При этом подкоренные выражения в (22), (23) всегда положительны.

Отметим, что при Ь^еВ^ коэффициенты ¿о, Ъ2 для первого и всех последующих интервалов оказываются заключенными в бочкообразную область, изображенную на рис. 9. На этом же рисунке отмечены зоны

ответственности функций Ршах (¿02) > 7 = 1»П,Ш за верхнюю границу коэффициента Ь\. Распределение зон для нижних границ Р^щ^ог), ] =1,П,Ш коэффициента Ь\ имеет вид, являющийся отображением зон для верхних границ справа налево и сверху вниз.

Благодаря определению области В® с помощью аналитических выражении (18)-(23) задача получения наиболее правдоподобной оценки вектора Ь^

(см. (17)) приобретает конкретные очертания. Для ее решения могут быть использованы стандартные методы условной

оптимизации.

После вычисления вектора оценок независимых

коэффициентов помощью (15) определяются оценки недостающих

коэффициентов для

5 = 2,5, затем по формуле (12) рассчитываются соответствующие оценки функции р(0, а в случае /> 1

и оценки полезного сигнала 0(г), пользуясь 0д и уравнением (И). Поскольку функции р(/) и 0(0 полностью определяются коэффициентами, входящими в Ь^, и начальным условием ©о, то рассчитанные с использованием наиболее правдоподобных оценок 0д значения полезного сигнала

также соответствуют критерию максимального

правдоподобия.

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

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

В настоящее главе в обозначении вектора

коэффициентов

полиномиального представления функции р(0 для упрощения записи опускается индекс 5. Трактуя Ь как случайный вектор, распределенный

Рис. 9. Область значений Ьд. ¿2 и определение верхней границы для Ь1 для всех 5 = 1,5.

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

точности оценки полезного сигнала 8(í) = ^9(í)T, c(í)p(0)T используются

апостериорные (условные по отношению к реализации измерений z) значения среднеквадратических ошибок (СКО) оценивания. Значения СКО определяются без учета погрешностей кусочно-полиномиальной аппроксимации (14).

Напомним, что при оценивании полезного сигнала вектор Ъ коэффициентов кусочно-полиномиального представления p(í) и начальное условие 6д рассматривались как неизвестные параметры, и для них не вводилась плотность распределения вероятности. Используемая здесь стохастическая интерпретация i и 0о служит исключительно для того, чтобы охарактеризовать точность

полученных оценок с помощью СКО. При этом нельзя не отметить, что

$ —*

наиболее правдоподобные оценки 0О,Ь , вообще говоря, не должны обеспечивать минимум СКО. Тем не менее, отдавая себе в этом отчет, в качестве меры точности полученных оценок полезного сигнала 0(/) решено воспользоваться такой распространенной характеристикой как СКО.

С учетом выражений (11), (14) СКО оценивания полезного сигнала 0(/) на любой момент времени нетрудно получить в результате линейного преобразования ковариационной матрицы ошибок

Т

f(Q0\b,z)dQ0

( T¿ 1 *Ь \ о _ f 00 f rb

(^е0)Т i Р% V I - J J -00 ,eo -0o, ,0O

-Гл

fib \z)db . (24)

Входящая в (24) плотность /(©о1 Ь,г) является гауссовской. Ее параметры

л. *

определяются по оценкам Ь, 8о и соответствующей ковариационной матрице, полученным с помощью ФК (первый этап решения, предложенного в четвертой главе). Плотность /(Ь\1) для Ье Е^ с точностью до нормировки

определяется как

/@|*)осехр| Аф-ъ?-Ъ)

Гауссовость

/(0q | b ,z) обеспечивает аналитическое вычисление интеграла по 0д.

Для вычисления в (24) интеграла по b применяется метод сеток, который в данном случае базируется на представлении плотности f(b\z) в виде

Kb^Ylf^W-b^),

где Ь<1> - узлы сетки значений вектора Ъ; </> - номер узла; м'</> осехр^---Ь)^Р^\ь<'>- веса, пропорциональные

значениям /(Ъ |г) при Ь =Ь<1> \ 5(-) - дельта-функция.

Непосредственный перебор всех узлов оказывается чрезмерно

трудоемким. Поэтому для вычисления ковариационной матрицы (24) используются значения Ь , удовлетворяющие двойному требованию

¿<'>еВ(5) Г! -Ъ?Р=^(Ь<1> -Ъ)<>Ц, (25)

где с/ - заранее заданный порог. Согласно (25) из области

отбираются

только те узлы ЬК1>, веса которых >р<'> существенно отличны от 0.

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

отбора узлов Ь<,:> в области (25). В ее основу положена представленная в первой главе процедура отбора значений вектора N с целочисленными

^ -Г 1 л _

компонентами в области ^-Щ Рм (М-Щ<д. В данной задаче она модифицирована с учетом первой части требования (24).

При значительном числе интервалов 5 кусочно-полиномиального представления функции р(<) для вычисления ковариационной матрицы (24) применяется упрощенная скользящая схема. Ее смысл заключается в том, что ограничения (12) учитываются не на всем времени решения задачи [0; 5Г], а на

более коротком отрезке т = включающем интервал [г^-1);/^]

кусочно-полиномиального представления функции р(*), для которого вычисляется СКО. Отрезки т состоят из фиксированного числа «У > 1 элементарных интервалов кусочно-полиномиального

представления функции р(<). Следствием учета ограничений на скользящем отрезке т является негладкий характер зависимости расчетных значений СКО от времени. Однако ввиду вспомогательного характера этих значений указанный недостаток не является принципиальным.

Разработанный метод вычисления СКО может применяться не только для

—* * —

наиболее правдоподобных значений Ь и 6о, но и для других оценок Ь и 9д.

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

В качестве простого примера на рис. 10 приведены результаты моделирования для задачи оценивания скалярного гармонического сигнала 0(*) = р(О по скалярным же измерениям г, в которых он смешивается с

коррелированной во времени помехой и шумом (см. верхний график). На среднем графике представлены ошибки оценок p(f) с использованием предложенного кусочно-полиномиального представления сигнала с учетом его непрерывной дифференцируемое™ и

ожидаемых ограничений на p(f) и p(í). Интервал кусочно-полиномиального представления Г был принят равным 1с. Для сравнения здесь также представлены ошибки оценок р((),

полученные из вектора Ь, вырабатываемого ФК без учета ограничений на р(Г) и p(f). На нижнем графике показаны рассчитанные с помощью описанного выше метода, условные по отношению к реализации измерений z, значения СКО.

Из представленных результатов видно многократное преимущество варианта решения, в котором учтены ограничения. Примечательно, что эффект от учета ограничений, с одной стороны, проявляется при приближении сигнала р(t) к границам ±1, а с другой - при прохождении p(f) через 0, когда производная p(f) достигает граничных значений ±р. В промежуточных точках каждый из этих факторов ослабляется, но зато они действуют совместно.

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

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

Рис. 10. Результаты оценивания функции р(/) по ее измерениям на фоне коррелированной и шумовой помех.

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

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

Заключение

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

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

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

вызванных действием нарушения. Полученные оценки открывают возможность восстановления системы с нарушением.

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

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

6. Обоснована двухэтапная схема оценивания сигнала на основе кусочно-полиномиальной модели с учетом ограничений и непрерывной дифференцируемости. На первом этапе полиномиальные коэффициенты оцениваются с помощью ФК. На втором - оценки полиномиальных коэффициентов уточняются с учетом заданных ограничений путем условной оптимизации с квадратичной целевой функцией.

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

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

9. Эффективность разработанных методов подтверждена результатами моделирования и обработки натурных данных. Так, методы оценивания сигналов при наличии дискретно распределенных помех и при неполном стохастическом описании были апробированы с использованием показаний приемников ГНС С и данных ММГ на автомобиле и на судне. Метод

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

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

Публикации по теме диссертации в журналах, рекомендованных ВАК

1. Степанов O.A., Кошаев Д.А. Решение задачи коррекции показаний навигационной системы по данным о береговой черте точечных ориентирах в рамках теории нелинейной фильтрации // Гироскопия и навигация. 1994. №3. С. 43-55.

2. Дмитриев С.П., Степанов O.A., Кошаев Д.А. Колеватов А.П. Эффективность использования инерциальных систем для исключения неоднозначности фазовых измерений спутниковых навигационных средств // Гироскопия и навигация. 1995. №4. С. 49-59.

3. Дмитриев С.П., Кошаев Д.А., Степанов O.A. Многоканальная фильтрация и ее применение для исключения неоднозначности при позиционировании объектов с помощью GPS // Известия РАН. Теория и системы управления. 1997. № 1.С. 65-70.

4. Степанов O.A., Кошаев Д.А. Применение неравенства Рао-Крамера в задачах нелинейного оценивания // Известия РАН. Теория и системы управления. 1997. №2. С. 65-72.

5. Кошаев Д.А. Сравнение нижних границ точности в задачах нелинейного оценивания // Известия РАН. Теория и системы управления. 1998. №2. С. 62-65.

6. Дмитриев С.П., Степанов O.A., Кошаев Д.А. Применение марковской теории нелинейной фильтрации при оптимизации алгоритмов навигации с использованием картографической информации о возможных траекториях движения объекта//Радиотехника, 1999. №10.

7. Степанов O.A., Кошаев Д.А. Исследование методов решения задачи ориентации с использованием спутниковых систем // Гироскопия и навигация. 1999. №2. С. 30-55.

8. Дмитриев С.П., Степанов O.A., Кошаев Д.А. Исследование способов комплексирования данных при построении инерциально-спутниковых систем //Гироскопия и навигация. 1999. №3. С. 36-52.

9. Дмитриев С.П., Степанов O.A., Ривкин Б.С., Кошаев Д.А., Чанг Д. Оптимальное решение задачи автомобильной навигации с использованием карты дорог //Гироскопия и навигация. 2000. №2. С.57-70.

10. Степанов O.A., Блажнов Б.А, Кошаев ДА. Исследование эффективности использования спутниковых измерений при определении ускорения силы тяжести на летательном аппарате //Гироскопия и навигация. 2002. №3.

П.Степанов О.А, Кошаев Д.А. Универсальные MATLAB-программы для анализа потенциальной точности и исследования чувствительности фильтров в интегрированных навигационных системах // Гироскопия и навигация. 2004. №2. С. 81-92.

12. Дмитриев С.П., Кошаев ДА. Информационный контроль и диагностика дублированных инерциальных систем // Гироскопия и навигация. 2005. №2. С. 30-46.

13. Кошаев ДА. Алгоритмы обнаружения и компенсации информационных нарушений в навигационных комплексах при большом числе типов нарушений и неопределенности их моделей // Гироскопия и навигация. 2007. №4. С. 47-59.

14. Блажнов Б.А., Кошаев Д.А. Определение относительной траектории движения и углов ориентации по фазовым спутниковым измерениям и данным микромеханического гироскопа // Гироскопия и навигация. 2009. №4. С. 15-34.

15. Кошаев Д.А. Многоальтернативный метод обнаружения и оценки нарушений на основе расширенного фильтра Калмана // Автоматика и Телемеханика 2010. №5. С. 70-83.

Основные публикации по теме диссертации в других изданиях

1. Кошаев Д.А. Экономичный алгоритм исключения неоднозначности фазовых спутниковых измерений // Тезисы докладов XX межотраслевой научно-технической конференции памяти H.H. Острякова Гироскопия и навигация. 1996. №4. С. 63-64.

2. Dmytriev S.P, Stepanov O.A., Koshaev D.A. Optimal Ambiguity Resolution and Efficiency in Using INS After Lock Loss // Proc. of 5th International Conference on Differential Satellite Navigation Systems. May 20-24, 1996, St.-Petersburg, Russia V. П. Paper 51.

3. Дмитриев С.П., Степанов О.А, Кошаев Д.А. Оптимальное разрешение неоднозначности фазовых измерений GPS с использованием ИНС И Труды Ш Санкт-Петербургской международной конференции по интегрированным навигационным системам. 28-29 мая 1996г., С.-Петербург, Россия.

4. Кошаев ДА. Экспериментальное исследование эффективности алгоритма исключения неоднозначности фазовых измерений спутниковой системы // Сборник докладов I научно-технической конференции молодых ученых "Навигация и управление движением". 6-7 апреля 1999г., С.-Петербург, Россия. С. 39-48.

5. Dmitriev S.P., Stepanov O.A., Rivkin B.S., Koshaev D.A, Chung D. Optimal Map-Matching For Car Navigation Systems // Proc. of 6th International Conference on Integrated Navigation Systems. May 24-26, 1999, St. Petersburg, Russia Paper 22.

6. Dmitriev S.P., Stepanov O.A., Rivkin B.S., Koshaev D.A., Chung D. Synthesis and Analysis of Map-Matching Algorithm for Car Navigation Systems // Proc. of Institute of Navigation 55th Annual Meeting. June 28-30,1999, Cambridge, MA.

7. Кошаев ДА. Анализ свойств ошибок измерений GPS на летательном аппарате // Сборник докладов II научно-технической конференции молодых ученых "Навигация и управление движением". 28-30 марта 2000г., С.-Петербург, Россия. С. 276-286.

8. Несенюк Л.П., Блажнов Б.А., Степанов О.А., Кошаев Д.А. Использование спутниковых измерений при определении ускорения силы тяжести на подвижном объекте // Сборник трудов научно-технической конференции «Навигация и океанография» (НО-2001). 2001, С.-Петербург, Россия.

9. Дмитриев С.П., Осипов А.В., Кошаев Д.А Контроль целостности интегрированных инерциально-спутниковых систем методами многоальтернативной фильтрации // Тезисы докладов Всероссийской конференции «Фундаментальное и прикладное координатно-временное и навигационное обеспечение» (КВО-2005), 11-15 апреля 2005г., С.-Петербург, Россия.

10.Dmitriev S.P., Osipov A.V., Koshaev D.A. Detection of information failures in marine navigation systems and their reconditioning // Proc. of 6th IF AC Symposium on Fault Detection, Supervision and Safety of Technical Processes. Aug. 30 - Sept. 1,2006, Beijing, P. R. China P. 1561-1566.

П.Дмитриев С.П., Осипов A.B., Кошаев Д.А. Обнаружение и устранение влияния информационных нарушений в навигационных системах // Тезисы докладов XXV межотраслевой научно-технической конференции памяти Н.Н. Острякова. Гироскопия и навигация. 2006. №4. С 96.

12. Дмитриев С.П., Осипов А.В., Кошаев Д. А. Учет неопределенности линейных моделей в задачах обработки навигационной информации // Сборник трудов VI Российской научно-технической конференции «Современное состояние и проблемы навигации и океанографии» (НО-2007). 23-27мая 2007г., С.-Петербург, Россия.

13.Кошаев Д.А. Экономичные алгоритмы обнаружения и устранения влияния информационных нарушений в навигационных комплексах // Труды 14-й Санкт-Петербургской международной конференции по интегрированным навигационным системам. 28-30 мая 2007г., С.-Петербург, Россия. С. 231-239.

14.Koshaev D.A. Efficient Fault Detection, Estimation and Recovery Using an Augmented Kalman Filter // Proc. of 13th IEEE/IF AC Conference on Methods and Models in Automation and Robotics. Aug. 27-30, 2007, Szczecin, Poland. P. 813-818.

15.Dmitriev S.P., Koshaev D.A. Filtering Algorithm with a Piecewise-Polynomial Description of Signals Taking Account of their Continuity and Level Constraint // Proc. of 16th IEEE Mediterranean Conference on Control and Automation. June 25-27,2008, Ajaccio, Corsica, France. P. 1646-1651.

Введение.

Глава 1. оценивание сигналов при наличии дискретно распределенных помех измерений.

1.1. Обзор существующих методов оценивания сигналов по измерениям с дискретно распределенными помехами.

1.2. Постановка и принципы решения задачи.

1.3. Рациональная организация отбора возможных значений дискретных помех.

1.4. Определение относительных координат по фазовым измерениям спутниковых систем.>.

1.5. Определение углов ориентации по фазовым измерениям спутниковых систем и данным микромеханического гироскопа.

Основные результаты первой главы.

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

2.1. Обзор классических многоальтернативных методов и пути снижения их трудоемкости.

2.2. Локализация и оценивание постоянно действующих нарушений.

2.3. Экономичная процедура пересчета результатов расширенного фильтра Калмана для получения многоальтернативного решения.

2.4. Обнаружение, локализация и оценивание нарушений с учетом возможности их возникновения и самоустранения.

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

Основные результаты второй главы.

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

3.1. Обзор методов обнаружение нарушений без использования стохастических моделей.

3.2. Постановка задачи обнаружения и оценивания аномальных ошибок • • • •

3.3. Полиномиальная модель для описания аномальной ошибки.

3.4. Решение на основе параллельно работающих банков фильтров Калмана •

3.5. Исследование эффективности метода на примере навигационного комплекса из трех инерциальных систем.

Основные результаты третьей главы.

Глава 4. метод оценивания сигналов на основе кусочнополиномиального представления с учетом ограничений и непрерывной дифференцируемости.

4.1. Обзор существующих методов оценивания нестохастических сигналов • •

4.2. Постановка задачи оценивания сигнала в условиях неполного стохастического описания.

4.3. Двухэтапная схема оценивания.

4.4. Определение области значений коэффициентов кусочно-полиномиального представления сигнала для одного интервала.

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

Основные результаты четвертой главы.

Глава 5. анализ точности оценивания сигналов при использовании кусочно-полиномиальной модели с соблюдением ограничений и непрерывной дифференцируемости.

5.1. Анализ априорного распределения значений кусочно-полиномиальной функции с учетом ограничений и непрерывной дифференцируемости •••

5.2. Определение среднеквадратических ошибок оценивания кусочно-полиномиальной функции с применением метода сеток.

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

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

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

Основные результаты пятой главы.

Актуальность темы

При обработке навигационной информации, как в реальном времени, так и в режиме постобработки, широко используются стохастические модели полезного сигнала и помех измерений. Наибольшее распространение получили модели в виде линейных уравнений динамики вектора состояния и линейных измерений, предпологающие гауссовский характер распределения вектора состояния и шумов измерений. Такие модели приводят к легко реализуемому фильтру Калмана (ФК) [2, 3, 9, 12, 16, 45, 56, 57, 71, 72, 77, 100, 102] с удобным рекуррентным правилом формирования оценки и ковариационной матрицы ее ошибки. Кроме того, имеются и хорошо проработанные методики анализа чувствительности ФК к разладкам между расчетной и действительной моделями вектора состояния и измерений [12, 57, 71].

Существует, однако, большое число задач, связанных с обработкой навигационной информации, решение которых не обеспечивается одними только калмановскими алгоритмами. К ним относятся задачи с нелинейной зависимостью измерений от оцениваемых величин и/или с нелинейными уравнениями, описывающими динамику вектора состояния [12, 17, 19, 21, 64, 66, 89, 90, 91]. Для оптимального решения таких нелинейных задач не существует столь же простого в вычислительном отношении метода как ФК, вследствие чего приходится разрабатывать субоптимальные методы, учитывающие специфику конкретной задачи. Здесь также возникают вычислительные трудности с анализом точности полученного решения, которые вынуждают прибегать к приближенным способам определения среднеквадратических ошибок оценивания [31, 63, 65, 66].

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

Ч = фкхк-1 + F(Pk >&) + wk> Ч = нкхк + >+ Ч • (В)

Здесь к - номер дискретного момента времени ^ = к At; At - интервал дискретности; Vk»wk " чисто-случайные последовательности; Ф^, F(-), G(-) - известные матрицы и вектор-функции; О - вектор известных или неизвестных параметров либо случайных величин. Величина р может иметь различную природу, которая и определяет тип неопределенности. В задаче требуется оценить значения функций F, G и/или вектора х.

Величина р может представлять собой случайный вектор и иметь t дискретное распределение вероятности, например, при наличии дискретных I помех измерений, нарушений или аномальных ошибок. В этом случае речь идет f о многоальтернативном стохастическом описании. Классические варианты решения таких задач сводятся к построению банка ФК, каждый из которых настроен на одну из возможных гипотез о значении р. Впервые такое решение в случае отсутствия переходов между гипотезами предложил D.T. Magill [130]. Впоследствии оно было детально проработано Д.Г. Лайниотисом [43]. Многоальтернативное решение, предусматривающее переходы от одной модели к другой (Interacting Multiple-Model) впервые представил Н.А.Р. Blom [81]. Сотрудниками ОАО «Концерн «ЦНИИ «Электроприбор» С.П. Дмитриевым, О.А. Степановым, Н.В. Колесовым, А.В. Осиновым также предложен ряд научно обоснованных вариантов многоальтернативных решений [20,22, 23].

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

Другой рассматриваемый в диссертации тип неопределенности не предусматривает наличие стохастической модели для величины р, входящей в уравнения (В). При этом рд. =р(^), где р(7) - функция, о которой известно лишь, что она непрерывно дифференцируема и для нее справедливы ограничения |р(/)| < 1, |p(i)| ^ р •

Следует отметить, что это достаточно типичная в навигационных задачах ситуация, когда* для ошибок измерительных средств, например, инерциальных навигационных систем (ИНС), приемников глобальной навигационной спутниковой системы (ГНСС), имеются достоверные стохастические модели, а детальная стохастическая модель подлежащих определению навигационных и динамических параметров ' неизвестна из-за недостатка статистического I материала. В то же время, для подвижного объекта, как правило, можно указать ряд ограничений, например, на составляющие линейных и угловых скоростей и ускорений.

На сегодняшний день существует несколько методов, позволяющих учесть неполноту стохастического описания. В частном случае, находят применение фильтры калмановского типа, учитывающие наличие неизвестного входного воздействия в уравнениях для вектора состояния. В англоязычной литературе их называют Unknown Input Kaiman Observers [118, 120, 137]. Этот тип фильтров полезен при полном отсутствии информации об ожидаемых свойствах возмущения.

В работах A.A. Красовского [38, 39] и О.С. Салычева [141] предлагаются робастные методы также на основе ФК, настроенного на квазидетерминированное описание оцениваемых процессов на ограниченных интервалах времени. Однако с помощью такого описания невозможно учесть 7 взаимосвязь квазидетерминированных функций на соседних участках и ограничения на оцениваемый процесс. Ожидаемый уровень процесса в фильтрах определяется с помощью задания дисперсий для параметров квазидетерминированного представления.

Для учета неполноты стохастической информации в задачах оценивания широко применяются гарантирующие методы, которые исходят из предположения о наихудшем для точности поведении сигнала при соблюдении некоторых детерминированных условий. Так существуют стационарные фильтры, построенные на основе критерия Н&, ориентированного на минимизацию наибольшей из возможных ошибок оценивания при заданных ограничениях на уровень или на дисперсии возмущений в уравнениях, описывающих динамику полезного сигнала [42, 54, 139]. В методе, предложенном A.B. Небыловым [40, 49], в качестве априорной информации используются предельные значения дисперсий производных от стационарного полезного сигнала. К сожалению, допущение о стационарности существенно ограничивает применимость таких методов в навигационных задачах, где зачастую приходится сталкиваться со сложным нестационарным поведением оцениваемого сигнала. Кроме того, для значительного числа навигационных приложений гарантирующий подход является излишне «жестким», поскольку применяемая в нем пессимистическая настройка, может не соответствовать физическому смыслу задачи.

Более общими являются методы гарантирующей фильтрации, предложенные А.Б. Куржанским [41] и А.И. Матасовым [44, 132], допускающие нестационарное описание вектора состояния и измерений. Вместо статистических данных о начальном условии для вектора состояния, о возмущениях, входящих в правую часть уравнения для вектора состояния, и о помехах измерений здесь используются ограничения. Однако и тут априорные представления об изменчивости оцениваемого процесса либо игнорируются, либо передаются излишне скупо, например, с помощью ограничения на производную от порождающего сигнал возмущения.

Для оценивания нестохастических сигналов с учетом их непрерывности и гладкости могут также применяться методы сплайн аппроксимации [1, 29, 78, 117], но они позволяют учесть ограничения на аппроксимируемый сигнал лишь косвенным путем — за счет ввода в минимизируемый критерий квадратов от производных аппроксимации с эмпирически выбранными весами.

Чтобы преодолеть недостатки существующих методов для нестохастических сигналов представляется целесообразным использовать кусочно-полиномиальное описание, коэффициенты которого находятся в области, обеспечивающей заданные свойства сигналов. Формализация такой задачи и доведение ее до реализуемого решения создает основу для построения эффективных алгоритмов комплексирования навигационной информации стохастического и детерминированного характера. Эти алгоритмы могут, например, использоваться для уточнения навигационных параметров подвижного объекта на основе очень ограниченных сведений о его динамических свойствах, таких как выдерживаемые при штатной эксплуатации пределы изменения линейной и угловой скоростей и ускорений, а также плавный характер изменения этих параметров. По существу это неинвариантное решение [18, 21, 72], которое не требует детальной стохастической модели динамики объекта.

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

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

Основные задачи

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

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

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

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

5. Разработка метода оценивания непрерывно дифференцируемого сигнала с учетом ограничений на величину сигнала и его производную. Получение адекватной характеристики точности.

6. Подтверждение эффективности предложенных методов по результатам моделирования и обработки реальных данных навигационных систем.

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

Научные положения, выносимые на защиту

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

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

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

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

5. Метод приближенного вычисления среднеквадратической ошибки оценивания сигнала, с учетом непрерывной' дифференцируемости и ограничений.

Научная новизна

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

Практическая значимость

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

Спроектирован алгоритм обнаружения и оценивания нарушений в навигационном комплексе, включающем две и более ИНС. С точки зрения качества получаемого решения алгоритм, не уступает классическому варианту на основе банка ФК для отдельных нарушений, но при этом требует значительно меньше вычислительных ресурсов. Эффективность алгоритма подтверждена путем моделирования задачи диагностики для навигационного комплекса типа «Аппассионата ЭКМ.1» и обработки реальных данных от ИНС «Дельта». ® Для навигационного комплекса из трех ИНС разработан алгоритм обнаружения и компенсации аномальных ошибок в показаниях курса без использования детального описания их динамики. Результаты стендовых испытаний навигационного комплекса типа «Ладога-11430», проведенные с имитацией различных типов нарушений, показывают, что алгоритм позволяет без значительного запаздывания обнаруживать и компенсировать дополнительную ошибку, вызванную нарушением, с точностью, сопоставимой с уровнем номинальных ошибок. ® Разработан алгоритм оценивания сигнала, модель которого опирается лишь на информацию о его непрерывности и гладкости, а также пределы изменения самого сигнала и его производной. При этом в полной мере учитываются корреляционные свойства ошибок измерителей: Алгоритм может использоваться в камеральном режиме для выработки эталонных данных при оценке точности навигационных средств, а также для высокоточной привязки измерений при съемке пространственных полей. Вместе с тем, он может применяться и в режиме реального времени на скользящем интервале для уточнения навигационных и» динамических параметров подвижных объектов. Результаты, моделирования- и обработки данных натурных испытаний демонстрируют существенный точностной . выигрыш за счет учета ограничений на динамические параметры объекта.

Применение результатов. К настоящему времени алгоритм обнаружения и оценивания аномальной ошибки внедрен в программное обеспечение навигационного комплекса «Ладога 11430» (2007г). Алгоритм исключения* неоднозначности фазовых измерений ГНСС применялся для компенсации вертикальных ускорений в показаниях аэрогравиметра «Чекан-А» (2002-2003г). Разработанные алгоритмы определения курса и относительной траектории по фазовым измерениям ГНСС были апробированы в ходе проведения испытаний в городских условиях и на оз. Ладога (тема «БАСН-М-ЭП», 2008-2009г).

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

Апробация работы. Результаты диссертации докладывались на III, VI, IX, X, XI, XIV и XVI Санкт-Петербургской международной конференции по интегрированным навигационным системам (С.-Петербург, 1996, 1999, 2002, 2003, 2004, 2007, 2009); XX, XXI, XXII, XXIII, и XXV научно-технической межотраслевой конференции памяти H.H. Острякова (С.-Петербург, 1996, 1998; 2000, 2002, 2006); I и II научно-технической конференции молодых ученых "Навигация и управление движением" (С.-Петербург, 1999, 2000); научно-технической конференции «Навигация и океанография» (С.-Петербург, 2001, 2007); форуме китайско-российских молодых ученых (Пекин, Китай, 2003); I и III Всероссийской конференции «Фундаментальное и прикладное координатно-временное обеспечение» (С.-Петербург, 2005, 2009); 5th International Conference on Differential Satellite Navigation Systems (С.-Петербург, 1996); Institute of Navigation 55th Annual Meeting (Boston, USA, 1999); 6th IFAC Symposium on Fault Detection, Supervision and Safety of Technical Processes (Beijing, PR China, 2006); 13th IEEE Conference "Methods and Models in Automation and Robotics" (Szczecin, Poland, 2007); 16th IEEE "Mediterranean Conference on Control and Automation" (Ajaccio, Corsica, France, 2008), на Общероссийском семинаре «Современные методы навигации и управления движением» (Москва, Институт проблем управления им. В.А. Трапезникова РАН, 2010).

Публикации. По материалам диссертации опубликовано 36 статей и докладов, включая 15 статей в журналах, рекомендованных ВАК.

Структура и объем работы. Диссертация состоит из пяти глав, введения, заключения и списка используемой литературы из 156 наименований. Общий объем работы составляет 242 страницы, включая 45 рисунков.

Основные результаты пятой главы

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

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

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

Заключение

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

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

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

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

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

6. Обоснована двухэтапная схема оценивания сигнала на основе кусочно

I \ полиномиальной модели с учетом ограничений и непрерывной дифференцируемости. На первом этапе полиномиальные коэффициенты оцениваются с помощью ФК. На втором - оценки полиномиальных коэффициентов уточняются с учетом заданных ограничений путем условной оптимизации с квадратичной целевой функцией.

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

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

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

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

1. Альберг Дж., Нильсон Э., Уолш Дж. Теория сплайнов и ее приложения. М.: Мир, 1972.

2. Анучин О.Н., Емельянцев Г.И. Интегрированные системы ориентации и навигации для морских подвижных объектов. Под общей ред. ПешехоноваВ.Г. С.-Петербург, ГНЦРФ ЦНИИ «Электроприбор», 1999.

3. Берман З.М., Вайсгант И.Б., Канушин В.М., Короленко A.B., Тупысев В.А., Шарыгин Б.Л. Преимущества инерциальной навигационной системы с фильтром калмановского типа в замкнутой схеме коррекции // Гироскопия и навигация. 1999. № 1. С. 48-56.

4. Бернелли-Заццера Ф., Молина М., Ванотти М. Недорогое оборудование для определения ориентации университетского спутника PalaMede по сигналам GPS // Гирскопия и навигация. 2001. №'4. С. 73-82.

5. Блажнов Б.А., Кошаев Д.А. Определение относительной траектории движения и углов ориентации по фазовым спутниковым измерениям иданным микромеханического гироскопа // Гироскопия и навигация. 2009. №4. с. 15-34.

6. Болотин Ю.А., Голован A.A., Парусников H.A. Информационное обеспечение задачи авиационной гравиметрии // Труды конференции по теории колебаний и управлению. 29 ноября 2000г. Москва, Россия.

7. Брайсон А., Хо Ю-Ши. Прикладная теория оптимального управления: М:: Мир, 1972.

8. Гантмахер Ф.Р. Теория матриц. М.: Наука, 1966.

9. Глобальная спутниковая радионавигационная система ГЛОНАСС / Под ред. Харисова В.Н., Перова А.И., Болдина В.А. М.: ИПРЖР, 1999.

10. Дмитриев С. П. Высокоточная морская навигация. СПб;: Судостроение, 1991.)

11. Дмитриев; С.П., Степанов О:А., Кошаев Д.А. Колеватов А.П. Эффективность (использования инерциальных систем' для исключения неоднозначности фазовых измерений спутниковых навигационных средств-//Гироскопия и навигация. 1995. №4. С. 49-59;:

12. Дмитриев С.П., Кошаев Д.А., Степанов O.A. Многоканальная фильтрация и ее' применение для исключения неоднозначности- при позиционировании объектов- с помощью GPS // Известия-РАН. Теория и системы управления. 1997. №1. С. 65-70.

13. Дмитриев С.П., Степанов O.A., Кошаев Д.А. Исследование способов комплексирования данных при ■ построении инерциально-спутниковых систем // Гироскопия и навигация. 1999. №3. С. 36-52.

14. Дмитриев С.П., Степанов O.A., Кошаев Д.А. Применение марковской теории нелинейной фильтрации при оптимизации алгоритмов навигации с использованием картографической информации о возможных траекториях движения объекта//Радиотехника, 1999. №10.

15. Дмитриев С.П., Степанов O.A. Неинвариантные алгоритмы обработки информации инерциальных навигационных систем // Гироскопия и навигация. 2000. №1. С. 24-38.

16. Дмитриев С.П., Степанов O.A., Ривкин Б.С., Кошаев Д.А., Чанг Д. Оптимальное решение задачи автомобильной навигации с использованием карты дорог // Гироскопия и навигация. 2000. №2. С.57-70.

17. Дмитриев С.П., Осипов A.B. Фильтрационный подход к задаче контроля целостности спутниковой навигационной системы // Радиотехника. 2002. №1. С. 39-47.

18. Дмитриев С.П., Пелевин А.Е. Задачи навигации и управления при стабилизации судна на траектории. СПб.: ГНЦ РФ ЦНИИ «Электроприбор», 2002.

19. Дмитриев С.П., Колесов Н.В., Осипов A.B. Информационная надежность, контроль и диагностика информационных систем. С.-Петербург, ГНЦ РФ ЦНИИ «Электроприбор», 2003.

20. Дмитриев С.П., Степанов O.A. Многоальтернативная фильтрация-в задачах обработки навигационной информации // Радиотехника. 2004. №7. С. 11-17.

21. Дмитриев С .П., Кошаев Д.А. Информационный контроль и диагностика дублированных инерциальных систем // Гироскопия и навигация. 2005. №2. С. 30-46.

22. Иванов Ю.П. Метод адаптивной комплексной оптимальной инвариантной! фильтрации сигналов // Известия ВУЗов. Приборостроение. 2003. Т. 46. №3. С. 3-8.

23. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука, 1968.

24. Кошаев Д.А. Экономичный алгоритм исключения неоднозначности фазовых спутниковых измерений // Тезисы докладов XX межотраслевой научно-технической конференции памяти H.H. Острякова. Гироскопия и навигация. 1996. №4. С. 63-64.

25. Кошаев Д.А. Сравнение нижних границ точности в задачах нелинейного оценивания // Известия РАН. Теория и системы управления. 1998. №2. С. 62-65.

26. Кошаев Д.А. Анализ свойств ошибок измерений GPS на летательном аппарате // Сборник докладов П научно-технической конференции молодых ученых "Навигация и управление движением". 28-30 марта 2000г., С.-Петербург, Россия. С. 276-286.

27. Кошаев Д.А. Исключение неоднозначности фазовых спутниковых измерений с использованием данных от инерциальных систем. Диссертация на соискание ученой степени к.т.н. ГНЦ РФ ФГУП ЦНИИ «Электроприбор», С.-Петербург, 2001.

28. Кошаев Д.А. Алгоритмы обнаружения и компенсации информационных нарушений в навигационных комплексах при большом числе типов нарушений и неопределенности их моделей // Гироскопия и навигация.2007. №4. С. 47-59.

29. Кошаев Д.А. Многоальтернативный метод обнаружения и оценки нарушений на основе расширенного фильтра Калмана // Автоматика и Телемеханика. 2010. №5. С. 70-83.

30. Красовский A.A. Адаптивный оптимальный регулятор с переменным порядком наблюдателя и временем экстраполяции // Автоматика и Телемеханика. 1994. № 11. С. 97-112.

31. Красовский A.A. Теория самоорганизующегося оптимального регулятора биномиального типа в детерминированно-стохастичеком приближении // Автоматика и Телемеханика. 1999. № 5. С. 97-113.

32. Кулакова В.И., Небылов A.B. Гарантрующее оценивание в системах с ограниченными дисперсиями производных // Автоматика и Телемеханика.2008. Хо 1.С. 83-96.41'. Куржанский А.Б., Управление и наблюдение в условиях неопределенности. М.: Наука, 1977.

33. Куркин О.М., Коробочкин Ю.Б., Шаталов С.А. Минимаксная обработка информации. М.: Энергоатомиздат, 1990.

34. Лайниотис Д.Г. Разделение единый метод построения адаптивных систем. I Оценивание. II Управление // Труды института инженеров по электротехнике и радиоэлектронике, 1976. Т.64. №8. С. 8-27, 74-94.

35. Матасов А.И. Метод гарантирующего оценивания. М.: Изд-во МГУ, 2009.

36. Медич Дж. Статистически оптимальные линейные оценки и управление М.: Энергия, 1973.

37. Миллер Г.Б., Панков А.Р. Минимаксная фильтрация в линейных неопределенно-стохастических дискретно-непрерывных системах // Автоматика-и Телемеханика. 2006. №3. С. 77-93.

38. Миллер Г.Б., Панков А.Р. Минимаксное управление процессом в линейной неопределенно-стохастической системе с неполными данными // Автоматика и Телемеханика. 2007. №11. С. 164-177.

39. Панков А.Р., Миллер Г.Б. Минимаксная линейная рекуррентная фильтрация неопределенно-стохастических последовательностей по интегральному критерию //Информ. процессы. 2001. Т. 1. №2. С. 150-166.

40. Пешехонов В.Г. Инерциальная навигация // Тезисы докладов; 3-й? Всероссийской конференции «Фундаментальное 1 и прикладное координатно-временноё и навигационное обеспечение (КВЫО-2009). 6-9 апреля 2009г., С.-Петербург, Россия. С. 29-32.

41. Поляк Б.Т. Введение в оптимизацию. М.: Наука, 1983.

42. Поляк Б.Т., Щербаков П.С. Робастная устойчивость и управление. М;: Наука, 2002.

43. Репин В.Г., Тартаковский Г.П. Статистический синтез при априорной неопределенности и адаптация* информационных систем. М.: Сов; Радио, 1977.

44. Ривкин С.С. Метод оптимальной фильтрации Калмана и его применение в инерциальных навигационных системах. JL: Судостроение, ч. 1, 1973, ч. 2, 1974. , ч

45. Ривкин С.С., Ивановский Р.ЮКостров A.B. Статистическая оптимизация; навигационных систем. JI.: Судостроение, 1976.

46. Рокафеллар Р. Выпуклый анализ. М:. Мир, 1973.

47. Сетевые спутниковые радионавигационные системы / Под ред. Шебшаевича B.C. М.: Радио и связь, 1993.

48. Соболь И.Н. Численные методы Монте-Карло. М.: Наука, 1973.

49. Соловьев Ю.А. Системы спутниковой навигации. Москва, ЭКО-ТРЕНЗ, 2003.

50. Справочник по теории? автоматического.- управления / Под ред. A.A. Красовского. М.: Наука, 1987.

51. Степанов O.A. Приближенные методы анализа потенциальной точности; в нелинейных навигационных задачах. Л.: ЦШ-1И «Румб», 1986.

52. Степанов O.A., Кошаев Д.А. Решение задачи коррекции показаний навигационной системы по данным- о береговой черте точечных ориентирах в рамках теории нелинейной фильтрации // Гироскопия и навигация. 1994. №3. С. 43-55.

53. Степанов O.A., Кошаев Д.А. Применение неравенства Рао-Крамера, в задачах нелинейного оценивания // Известия РАН. Теория и системы управления. 1997. №2. С. 65-72.

54. Степанов О. А. Применение теории нелинейной фильтрации в задачах обработки навигационной информации. СПб: ГНЦ РФ ЦНИИ' «Электроприбор», 1998.

55. Степанов O.A., Кошаев Д.А. Исследование методов решения задача ориентации с использованием спутниковых систем // Гироскопия и навигация. 1999. №2. С. 30-55.

56. Степанов O.A., Блажнов Б.А., Кошаев Д.А. Исследование эффективности использования спутниковых измерений' при определений ускорения силы а тяжести на летательном аппарате // Гироскопия и навигация. 2002. №3.

57. Степанов O.A., Кошаев Д.А. Определение ускорения силы тяжести на летательном аппарате с применением измерений спутниковой системы // Сборник докладов на форуме китайско-российских молодых ученых. Пекин. 2003.

58. Степанов O.A., Кошаев Д.А. Универсальные MATLAB-программы для анализа потенциальной точности и исследования чувствительности фильтров в интегрированных навигационных системах // Гироскопия и навигация. 2004. №2. С. 81-92.

59. Степанов O.A. Основы теории оценивания с приложениями к задачам обработки навигационной информации. Часть 1. Введение в теорию оценивания. СПб.: ГНЦ РФ ЦНИИ «Электроприбор», 2009.

60. Химмельблау Д. Прикладное нелинейное программирование. М.: Мир, 1975.

61. Хлебников М.В. Робастная фильтрация при неслучайных возмущениях: метод инвариантных эллипсоидов // Автоматика и Телемеханика. 2009. № 1.С. 147-161.

62. Цыпкин Я.З. Информационная теория идентификации. М.: Наука, 1995.

63. Челпанов И.Б. Оптимальная обработка сигналов в навигационных системах. М.: Наука, 1967.

64. Ярлыков М.С., Миронов М.А. Марковская теория оценивания случайных процессов. М.: Радио и связь, 1993.

65. Baang D., Stoev J., Choi J.Y., Park J. Simplified adaptive nonlinear observer using B-spline based approximators // Proc. of 16th IF AC World Congress. July 4-9, 2005, Prague, Czech Republic.

66. Bar-Shalom Y., Li X.R. Estimation and Tracking: Principles, Techniques, and Software. Boston. MA. Artech House. 1993 (Reprinted by YBS Publishing. 1998).

67. Bassevile. M., Nikiforov I. Detection of Abrupt Changes: Theory and Application. Prentice Hall. 1993.

68. Blom H.A.P. An Efficient Filter for Abruptly Changing Systems // Proc. 23nd IEEE Conference on Decision and Control. 1984, Las Vegas. P. 656-658.

69. Blom H.A.P., Bar-Shalom Y. The interacting multiple model algorithm for systems with Markovian switching coefficients // IEEE Transactions on Automatic Control. Aug. 1988. V. 33. No. 8. P. 780-783.

70. Buist P. The Baseline Constrained LAMBDA Method for Single Epoch, Single Frequency Attitude Determination Applications // Proc. of ION GNSS 20th International Meeting of the Satellite Division, Sept. 25-28, 2007, Fort Worth, TX. P. 2962-2973.

71. Chen D., Lachapelle G. A Comparison of the FASF and Least-Squares Search Algorithms for Ambiguity Resolution On The Fly // Proc. of the International234

72. Symposium on Kinematic Systems in Geodesy, Geomatics and Navigation. Aug. 30 Sept. 2, 1994, Banff, Canada. P. 241-253.

73. Crassidis J.L., Lightsey E.G., Markley F.L. Efficient and Optimal Attitude Determination Using Recursive Global Positioning System Signal Operations // Journal of Guidance, Control, and Dynamics. 1999. V. 22. No. 2. P. 193-201.

74. De Jong P.J., Tiberius C.C.J.M., Teunissen P.J.G. Computational Aspects of the LAMBDA Method for GPS Ambiguity Resolution // Proc. of ION GPS-96. Sept. 17-20, 1996, Kansas City, Missouri.

75. De Weerd E., Van Kampen E., Chu Q.P., Mulder J.A. New Approach for Integer Ambiguity Resolution using Interval Analysis // NAVIGATION: Journal of The Institute of Navigation. 2008. V. 55. No. 4. P. 293-307.

76. Dmytriev S.P, Stepanov O.A., Koshaev D.A. Optimal Ambiguity Resolution andth

77. Efficiency in Using INS After Lock Loss // Proc. of 5 International Conference on Differential Satellite Navigation Systems. May 20-24, 1996, St.-Petersburg, Russia. V. II. Paper 51.

78. Dmitriev S.P., Stepanov O.A., Rivkin B.S., Koshaev D.A., Chung D. Optimal Map-Matching For Car Navigation Systems // Proc. of 6th International Conference on Integrated Navigation Systems. May 24-26, 1999, St. Petersburg, Russia. Paper 22.

79. Dmitriev S.P., Stepanov O.A., Rivkin B.S., Koshaev D.A., Chung D. Synthesis and Analysis of Map-Matching Algorithm for Car Navigation Systems // Proc. of Institute of Navigation 55th Annual Meeting. June 28-30, 1999, Cambridge, MA.

80. Dmitriev S.P., Stepanov O.A., Koshaev D.A. Map Matching for Automobile Navigation. // GIM International. The world Magazine for Geomatics. 2000. V. 14. No. 7.

81. Dmitriev S.P., Osipov A.V. Koshaev D.A. Detection of information failures in marine navigation systems and their reconditioning // Proc. of 6th IFAC

82. Symposium on Fault Detection, Supervision1 and Safety of Technical Processes. Aug. 30 Sept 1, 2006, Beijing, P. R. China. P. 1561-1566.

83. Fisher K.A., Maybeck P.S.' Multiple Model Adaptive Estimation with Filter$

84. Spawning // IEEE Transactions on Aerospace and- Electronic Systems. 2002.1. V. 38, No. 3. P. 755-768.i

85. Fletcher R. Practical methods for optimization. Volume 2: Constrained optimization. Wiley & Sons Inc. 1981.

86. Fouladirad M., Nikiforov I. On Line Change Detection with Nuisance , Parameters // Proc. of 6th IF AG Symposium on Fault Detection, Supervision and Safety of Technical Processes. Aug 30 Sept 1, 2006, Beijing, P. R. China. j P. 289-294. j

87. Frank P.M. Fault Diagnosis in Dynamic Systems Using Analytical and Knowledge Based Redundancy A Survey and Some New Results // ' j Automatica. 1990. V. 26. P: 459-474.

88. Frei E. and' Beutler G. Rapid static positioning based on the1 fast ambiguity approach FARA: theory and first results // Manuscripta Geodetica, 1991. P.325-356.

89. Gelb A. Applied Optimal Estimation. M.I.T. Thirteenth printing. Cambrige. MA. 1994.1 101. Global Positioning System: Theory and Applications / Edited by

90. B. W. Parkinson, J. J. Spliker Jr. V. I, II. 1996.

91. Grewal M.S., Andrews A.P. Kaiman Filtering. Theory and Practice. Prentice Hall, New Jersy. 1993.

92. Gustafsson F. Statistical Signal Processing Approaches to Fault Detection // Proc. of 6th IF AC Symposium on Fault Detection, Supervision and Safety of Technical Processes. Aug 30 Sept 1, 2006, Beijing, P. R. China. P. 37-48.

93. Han S., Wong K., Rizos C. Instantaneous Ambiguity Resolution for Real-Time GPS Attitude Determination // Proc. of the International Symposium on Kinematic Systems in Geodesy, Geomatics and Navigation. June 3-6, 1997, Banff, Canada. P. 409-416.

94. Hanlon P.D., Maybeck P.S. Characterization of Kaiman Filter Residuals in the Presence of Mismodeling // IEEE Trans. Aerospace Electronic Syst. 2000. V. 36. No. l.P. 114-131.

95. Hanlon P.D., Maybeck P.S. Interrelationship of Single-Filter and Multiple-Model Adaptive Algorithms // IEEE Trans. Aerospace Electronic Syst. 1998. Y. 34. No. 3. P. 934-946.

96. Hanlon P.D., Maybeck P.S. Multiple-Model Adaptive Estimation Using a Residual Correlation Kaiman Filter Bank // IEEE Trans. Aerospace Electronic Syst. 2000. V. 36. No. 1. P. 393-406.

97. Hein G.W., Werner W. Comparison of Different On-The Fly Ambiguity Resolution Techniques // Proc. of ION GPS-95. Sept. 12-15, 1995, Palm Springs, California. P. 1137-1144.

98. Hirokawa R., Ebinuma T. A Low-Cost Tightly Coupled GPS/INS for Small UAVs Augmented with Multiple GPS Antennas // NAVIGATION: Journal of The Institute of Navigation. 2009. V. 56. No. 1. P. 35-44.

99. Jamouli H, Sauter D.A New adaptive Kalman estimator integrated in a fault-tolerant control system // Proc. of 15th IEEE Mediterranean Conference on control and Automation, July 27-29, 2007, Athens, Greece.

100. Jordan S.K. Self-consistent statistical models for gravity anomaly and undulation of the geoid. J. Geophys. Res. 1972. V. 77. No. 20. P. 3660-3670.

101. Kano Hi,. Fujioka H., Egerstedt M., Martin C.F. Optimal smoothing spline curves and contour synthesis // Proc. of 16th IF AC World Congress. July 4-9, 2005, Prague, Czech Republic.

102. Kitanidis P.K. Unbiased; minimum-variance linear state estimation // Automatica. 1987. V. 23. P. 775-778.

103. Koshaev D.A. Efficient Fault Detection; Estimation and; Recovery Using an Augmented Kalman Filter // Proc. of 13 th IEEE/IF AC Conference on Methods and Models in Automation and'Robotics. Aug. 27-30, 2007, Szczecin, Poland. P. 813-818.

104. Li L., Zhou D.H., Liu K.D. Robust Fault Diagnosis of Nonlinear Systems Based on an Unknown Input Extended Kalman Observer // Proc. of 16th IF AC World Congress. July 4-9, 2005, Prague, Czech Republic.

105. Landau H., Euler H.J On-the-fly ambiguity resolution for precise differential positioning//Proc. of ION GPS-92, Alexandrian'

106. Lau L., Mok E. Improvement of GPS Relative Carrier Phase Positioning Accuracy by SNR Weighted Ambiguity Function Technique // Proc. of Institute of Navigation 55th Annual Meeting. 28-30 June,1999, Cambridge, MA. P. 323-332.

107. Lau L., Mok E. Precise GPS Processing Algorithm for Short Observation TimeiL

108. Span in Urban Area // Proc. of Institute of Navigation 55 Annual Meeting. June 28-30,1999, Cambridge, MA. P. 533-540.

109. Li X.R., Jilkov V.P. Survey of maneuvering target tracking; Part V: : Multiple-model methods // IEEE Transactions on Aerospace and! Electronic Systems. 2005. V. 41. No. 4. P. 1255-1321. '

110. Li Y., Zhang K., Grenfell R. Improoved Knight Method Based on Narrowed. Search Space for Instaneous GPS Attitude Determination // NAVIGATION: Journal of The Institute ofNavigation. 2005. V. 52. No. 2. P. 111-119.

111. Lightsey E.G., Crassidis J.L., Markley F.L. Fast Integer Ambiguity Resolution for GPS Attitude Detenninatuon // Proc. of the AIA A Guidance, Navigation and; Control Conference. 1999, Portland, OR. V. 1. P. 115-123.

112. Liu. Z.W., Reed M.D., Lapucha D;R. GPS Gyro Integration for Airborne Attitude Reference // Proc. of the International Symposium on Kinematic Systems in Geodesy, Geomatics and Navigation; June 3-6, 199.7, Banff, Canada;. P. 215-221. ". ' . .

113. Macabiau C. Comparsion of the LS AST and MAP AS methods for ambiguity resolution on-the-fly II Proc. of the 5 International Conference on. Differential Satellite Navigation Systems (DSNS-96), 20-24; May,: 1996, St.-Petersburg, Russia. Paper№50:

114. Madcr G.L. Kinematic GPS Phase Initialisation Using Ambiguity Function // Proc. of the 6th International Symposium on Satellite Positioning: March 17-20, 1992, Columbus, Ohio.

115. Magill D.T. Optimal-adaptive estimation of sampled'stochastic processes //IEEE Transactions on Automatic Control. 1965. V. 10. Issue 4. P. 434-439.

116. Markley F.L. Attitude Determination Using Vector Observations: A Fast Optimal Matrix Algorithm II The Journal of the Astronautical Sciences. 1993. V. 41.No. 2. P. 261-280.

117. Matasov A.I. Estimators for Uncertain Dynamic Systems. Kluwer Academic Publishers. Dordrecht. The Netherlands. 1998.

118. Maybeck P.S., Kozak M.C., Smith B.D. Mixed-Model Multiple-Hypothesis Tracking of Targets in Clutter // IEEE Trans. Aerospace Electronic Syst. 2008. < V. 44. No. 4. P. 1402-1415.

119. Nadler A., Bar-Itzhack Y., Weiss H. Iterative Algorithm for Attitude Estimation Using Global Positioning System Phase Measurements // Journal of Guidance,i

120. Control, and Dynamics. 2001. V. 24. No. 5. P. 983-990.

121. Nguang S.K., Shi P., Ding S.X. Fault detection filter for uncertain fuzzy systems: an LMI approach // Proc. of 16th IFAC World Congress. Prague, July 4-8, 2005, Czech Republic.

122. Park K. J., Crassidis J. L. Attitude Determination Methods Using Pseudolite Signal Phase Measurements // NAVIGATION: Journal of The Institute of Navigation. 2006. V. 53. No. 2. P. 121-133.

123. Remondi B.W. Pseudo-kinematic GPS results using the ambiguity function method. National information center, Rockville, Maryland, NOOA Technical memorandum NOS NGS-52, 1990.

124. Rios B.A., Rivas-Echeverria F., Garcia G. Robust FDI with mixed H^/H«, criteria for discrete-time linear systems // Proc. of 16th IFAC World Congress. July 4-8, 2005, Prague, Czech Republic.

125. Rizos C., Shaowei H. A new method for constructing multi-satellite ambiguity combinations for improved ambiguity resolution // Proc. of ION GPS-95. Sept. 12-15, 1995, Palm Springs, California. P." 1145-1153.

126. Salychev O.S. Inertial Systems in Navigation and Geophysics, Bauman Moscow State Technical University Press, Moscow, Russia, 1988.

127. Scherzinger B.M. Precise Robust Positioning with Inertially Aided RTK // NAVIGATION: Journal of The Institute of Navigation. 2006. V. 53, No. 2. P. 73-83.

128. Shumsky A. Parity relation method and its application to fault detection in nonlinear dynamic systems // Autom. Remote Contr. 1998. V. 59. No. 9: P. 1327-1335.

129. Shuster D. M. Maximum Likelihood Estimation of Spacecraft Attitude // The Journal of the Astronautical Sciences. 1989. V. 37. No. 1. P. 79-88.

130. Teunissen P.J.G. A New Method for Fast Carrier Phase Ambiguity Estimation // Proc. of Position Location and Navigation Symposium. Las Vegas, Nevada, April 11-15, 1994.

131. Teunissen P.J.G., De Jouge P.J., Tiberius, C.C.J.M. The Volume of the GPS Ambiguity Search Space and its Relevance for Integer Ambiguity Resolution // Proc. of ION GPS-96, Kansas City, Missoury, September 17-20, 1996.

132. Teunissen P.J.G. The LAMBDA Method for the GNSS Compass // Artificial Satellites. 2006. V. 41. No. 3. P. 89-103.

133. Van Kampen E., De Weerd E., Chu Q.P., Mulder J.A. Applied Interval Based Integer Ambiguity Resolution // NAVIGATION: Journal of The Institute of Navigation. 2009. V. 56. No. 3. P. 205-219.

134. Vasquez J.R., Maybeck P.S. Enhanced Motion and Sizing of Bank in Moving-Bank MMAE // IEEE Trans. Aerospace Electronic Syst. 2000. V. 36. No. 1. P. 770-779.

135. Wahba G. A Least-Squares Estimate of Spacecraft Attitude // SIAM Review. 1965. V. 7. No. 3. P. 409.

136. Ward L.M., Axelrad P. A Combined Filter for GPS-Based Attitude and Baseline Estimation//Proc. of ION GPS-96, Sept. 17-20, 1996, Kansas City, Missouri.

137. Willsky A.S., Jones H.L. A Generalized Likelihood Ratio Approach to Detection and Estimation of Jumps in Linear Systems // IEEE Trans. Autom. Control. 1976. V.21.P. 108-112.

138. Zhong M., Ye H., Ma C., Wang G. Robust FDI for uncertain systems of the politopic type // Proc. of 16th IF AC World Congress. Prague, Czech Republic, July 4-8, 2005.