Описание программы исследований

Литература

1. Белоглазов И. Н., Джанджгава Г. И., Чигин Г. П. Основы навига­ции по геофизическим полям. - М.: Наука, 1985.-328 с.

2.Белоглазов И. Н., Тарасенко В. П. Корреляционно - экстремальные системы. - м.: Сов. радио, 1974. - 392 с.

{ESC} - Выход из программы или статистических исследований.

Вопросы для экпериментальных исследований

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

2. Используя первое окно интерактивного интерфейса приближенно подобрать дисперсию помехи ТИ (при размере фрагмента ЭИ - ширины участка поля для поиска в 20 элементов) для получения приблизительно половины успешных исходов экспериментов с правильным обнаружением за­даного фрагмента ЭИ. Зафиксировать числовые данные первого окна инте­рактивного интерфейса.

3. Для зафиксированных числовых данных первого окна интерактивно­го интерфейса п.1 и п.2 провести статистические исследования с исполь­зованием гистограмм при задании возможно наибольшего количества Nmax статистических экспериментов.

4. Исследовать сходимость результатов измерений гистограмм при возрастании количества Ni статистических экспериментов. Сходимость оп­ределить как среднеквадратическую разностьб² значений "соседних" по Ni (то есть для Ni и Ni+i, 1=1...n, Ni<Ni+i) гистограмм.

Содержание и порядок оформления отчета

Отчет оформляется в соответствии с СТП ГУАП и должен содержать следующие материалы:

-Теоретические сведения и алгоритмы работы поисковых КЭСН.

-План экспериментальных исследований по каждому пункту проведе­ния лабораторной работы.

-14-

 

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

Цель работы: изучение принципов работы и экспериментальное иссле­дование корреляционно-экстремальных систем (КЭСН) навигации по картам геофизических полей.

1. Теоретические сведения

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

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

В последнее время КЭНС находят широкое признание и применение. Они имеют разнообразную структуру и могут иметь различное предназначе­ние. Одним из удобных способов классификации КЭНС является классифика­ция по объему или характеру полезной информации, снимаемой с навигаци­онного поля бортовым датчиком в каждый момент времени или за короткий цикл сканирования. В соответствии с этим КЭНС делятся на три основных класса (рис. 1):

КЭНС-1 - в системах этого класса рабочая информация в текущий мо­мент времени снимается в "точке", т.е. выходной сигнал датчика поля

-3-

представляет скалярную величину;КЭНС-II - для систем данного класса характерным является то, что бортовой датчик за короткий цикл сканирования мгновенно снимает информацию о поле вдоль произвольных, но заранее выбранных, линий. Та­ковыми, в частности, могут служить линии, перпендикулярные к продоль­ной (в случае поперечного сканирования антенной зондируемой поверхнос­ти) или поперечной (при продольном сканирмовании) оси летательного аппарата;КЭНС-III - к этому классу относятся системы, в которых рабочая информация снимается с участка площади земной поверхности ("кадра"),просматриваемого бортовой системой в каждый момент времени или за короткий цикл сканирования (КЭНС с сопоставлением изображений).

По методу определения отклонения от экстремума КЭНС подразделяются на эвристические, поисковые ,беспоисковые оптимальные (субоптимальные) и комбинированные. Такая классификация в некоторой мере соответс­твует алгоритмам КЭНС, которые могут быть реализованы в КЭНС-I, КЭНС-П и КЭНС-III на базе аналоговой, аналого-цифровой и цифровой тех­ники.

В теоретических исследованиях и математическом моделировании на первом этапе в развитии теории и разрабоке КЭНС рассматривалась в ос­новном четырехточечная дифференциальная схема формирования сигналов коррекции отклонения от экстремума КФ[1].

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

Второе направление в развитии теории КЭНС составили поисковые ал­горитмы.Теория таких систем строится на основе широко известной теории статистических решений. Особенностью поисковых алгоритмов

-4-

 

линии полёта ЛА:”Ширина основного поля_(1000..10000)". Здесь и

далее (1000..10000) - размер поля в дискретных элементах. ЭИ задается в виде фиксированной реализации случайного процесса, типичной для изотропно отражающих земных поверхностей.Дисперсия реализации ЭИ нормирована к 1.

-Размер фрагмента ЭИ”Ширина участка полядля поиска_(10..300)”и начальную координату его расположения на ЭИ:”Начальное значение(координаты) поля для поиска_(1..900)”. Этот фрагмент и подлежит обнаружению в процессе выполнения лабораторной работы.

-Текущее изображение'(ТИ) как сумму выделенного фрагмента ЭИ и помехи. Интенсивность помехи задается в полях интерактивного интерфей­са: "Введите зашумленность в процентах (О..100)" или "Введите диспер­сию шума (О..1). Тип функции корреляции (ФК) помехи задается клавишами F1 - F3.

Переход к статистическим исследованиям (второму окну интерфейса) - клавиша F9. Во втором окне интерфейса можно задать количество ста­тистических экспериментов для заданной реализации ЭИ.

В этом разделе исследуются гистограммы

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

- вероятности правильного обнаружения фрагмента ЭИ от размера

ши­рины участка поля для поиска при ранее заданной и фиксированной дис­персии помехи, зашумляющей ТИ.

В программе также используются клавиши:

{TAB}-Переключение между полями ввода.

{F1}-{F3} -Выбор функции корреляции помехи ЭИ.

(F9) - Переход в режим статистических исследованийй.

{ENTER}-Начало поиска фрагмента ЭИ или расчета гистограмм.

-13-

 

Видно, что все функционалы имеют абсолютный экстремум для прове­ряемой ошибки 1Δ1, совпадающей с действительной ошибкой ГНС ΔX=X*N-XN. В этом случае I2i=I3i=0 и эти функционалы достигают абсолютного мини­мума. Функционал 1ц достигает на истинной гипотезе абсолютного макси­мума. Обоснованием функционала I1i является известное свойство корре­ляционной функции Rff(Δ) стационарного процесса f(x) достигать абсо­лютного максимума при нулевом сдвиге Δ=0. Действительно, если поле стационарное и эргодическое, а длина реализации NL достаточна для про­явления эргодичности, то

и абсолютный максимум I1i достигается при iΔ1=X*N-XN=ΔX. В лабораторной работе при экспериментальных исследованиях исполь­зуется алгоритм I1i.

Описание программы исследований

Программа Finder - эмулятор поисковой КЭНС для одномерного случая.

Программа позволяет задавать:

-Размер основного поля (эталонного изображения-ЭИ) вдоль

-12-

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

Однако реализация данных алгоритмов, особенно в случае решения задачи многопараметрического оценивания,требует значительной производительности БЦВМ. В данной лабораторной работе исследуется КЭНС-IIБ (см. рис. 1), которая является аналогомсистемы навигации по рельефу местности типа ТЕRСОМ. Метод навигации и наведения посредством сопоставления профилей рельефа местности ТЕRСОМ (Теrrian Сontour Matching) предназначен для применения на борту различных летательных аппаратов: пилотируемых и беспилотных самолетов, крылатых ракет, головных частей баллистических ракет. В основу метода построения системы ТЕRСОМ положен следующий принцип: географическое положение любой точки земной поверхности (су­ши) описывается единственным образом с помощью вертикальных профилей или топографией окружающей местности аналогично тому, как каждый чело­век обладает своими неповторимыми отпечатками пальцев. Как и при ана­лизе отпечатков пальцев, такая система требует предварительного кар­тографирования или иного определения характеристик профилей поверхнос­ти того района, над которым система будет использоваться. Это может быть сделано, например, по стереоскопическим аэроснимкам местности с помощью известных в настоящее время способов. Предварительно опреде­ленные данные о рельефе района полетов запоминаются в цифровом виде в бортовом запоминающем устройстве.

Во время полета система ТЕRСОМ измеряет вертикальный профиль рельефа вдоль истинной траектории полета с помощью радиовысотомера (для измерения геометрической высоты полета) и баровысотомера (для по­лучения уровня отсчета профиля).Вычитая текущую высоту,измеренную радиовысотомером,из барометрической высоты ,система определяет профиль рельефа вдоль траектории полёта

-5-

 

Корреляционно-экстремальные навигационные

системы (КЭНС)

По виду рабочей информации

 

КЭНС-I КЭНС-II КЭНС-III

(информация снимается ( информация снимается информация

в «точке») вдоль линии ) снимается с

кадра

 

По объёму исходной (априорной) информацией

 

КЭНС-IА КЭНС-IБ КЭНС-IIА КЭНС-IIБ КЭНС-IIIА КЭНС-II (системы (системы (системы системы системы системы

без памяти) с памятью) без памяти) с памятью без с па-

памяти) мятью

 

По способу хранения и обработки информации

Аналоговые Аналого-цифровые Цифровые

(непрерывные) (комбинированные)

 

По методу определения отклонения от экстремума

 

Эвристические Поисковые Беспоисковые Комбинированные

 

Рис.1

 

 

-6-

 

 

Могут использоваться различные виды функционалов Ii, но они долж­ны обладать общим обязательным свойством: в идеальных условиях, при отсутствии ошибок измерения поля и картографирования и других возмуще­ний, функционал должен достигать абсолютного (глобального) экстремума на истинной гипотезе, т.е. при ΔX=X*N –Xn. Тогда и в условиях неизбеж­ных возмущений, рассчитав значения функционала для всех возможных зна­чений ошибки ГНС и выбрав то значение ошибки Δх, при котором функцио­нал достиг экстремума, мы получим оценку местоположения, близкую к ис­тинной. Точность этой оценки будет тем выше, чем меньше ошибки измере­ния геофизического поля и больше длина обрабатываемой реализации (чис­ло измерений поля N).

В качестве функционалов в эвристических поисковых алгоритмах испольэуются:

среднее произведение

средняя абсолютная разность

средний квадрат разности

Если ошибки картографирования отсутствуют, т.е. в одних и тех же точках fn(x)=f(x), где f(x) - истинное значение геофизического поля, а датчик поля производит измерения абсолютно точно ZN-k=f(XN-kL) , то выражения для функционалов принимают вид

 

-11-

 

шаг дискретизации по коор­динате Δ1 и в доверительном интервале с

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

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

условие может и не выполняться. Однако если величина Δ1 невелика, введение дискретизации при проверке гипотез не приво­дит к существенным ошибкам. В рассматриваемом случае гипотезами Hi являются возможные дискретные значения iΔ1 продольной ошибки ГНС Δх (i= -n, ...,n; n=Δmax/Δ1). Поскольку мы предположили отсутствие ошибки измерения скорости и постоянство Δх в процессе проведения измерений геофизического поля, то с точностью до шумов датчика поля и ошибок картографирования измеренная реализация совпадает с некоторым участком реализации навигационного поля, хранящейся в блоке памяти.

В процессе перебора различных возможных значений ошибок ГНС (т.е. проверки различных возможных гипотез Hi, состоящих в предположении, что ошибка Δх=iΔ1) из блока памяти извлекается для каждой проверяемой гипотезы соответствующая ей (гипотетическая) реализация поля

{fn[x*N-iΔ1-(N-1)L],fn[x*N-iΔ1-(N-2)L],…,fn(x*N-iΔ1-L),fn(x*N-iΔ1)}. Сопоставляется гипотетическая реализация поля с измеренной z1=f[XN-(N-1)L], ..., ZN-1=f(XN-L), ZN=f(XN), где XN - действительное местоположение ЛА в момент проведения последнего замера, L - расстоя­ние, пролетаемое ЛА между последовательными измерениями. В результате такого сопоставления получается некоторое число Ii, которое характеризует степень совпадения измеренной и проверяемой гипотетической реали­зации. Величина Ii=I(Hi) зависит от измеренной и гипотетической реализаций навигационного поля. Кроме того, Ii является функцией проверяе­мого значения Δx=iΔl ошибки навигационной системы.

-10-

 

организует поиск в памяти вы­числителя наиболее "похожего", заранее запомненного профиля с извест­ными координатами. Поиск осуществляется последовательным сопоставлени­ем сигнала ТИ датчика поля с цифровым образом ЭИ и просмотром сначала назад и вперед предполагаемого местоположения, а затем поперек траек­тории в пределах матрицы 64Х64 ячейки (при размере ячейки ∆ от 30 Х 30м до 240м Х 240м).

Процедура сопоставления представляет собой минимизацию абсолютной разности, приближенно определяемой выражением

где Нист -превышение рельефа вдоль истинной траектории полёта,

Hкарт- превышение рельефа ,выбранное из бортовой карты,Λi,Φi-

координаты летательного аппарата.

Конкретному исследованию в лабораторной работе подвергаются

качественные показатели обнаружения заданного участка ЭИ по предъявляемому ТИ в поисковой КЭНС.

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

-7-

Может применяться "слепой поиск", когда вводится дискретизация координат и значения функционала рассчитываются абсолютно для всех возможных.дискретных значений ошибок ГНС (главной навигационной, нап­ример инерциальной) системы. Такой метод поиска, являющийся в опреде­ленных условиях единственно возможным, надежно решает проблему ликви­дации больших начальных отклонений в КЭНС, хотя в некоторых случаях требует чрезмерно высокой производительности БЦВМ.Другая группа методов поиска связана с вычислением градиентов (или их аналогов) функционала, определяющих направление движения к экстремуму. Здесь требования к быстродействию БЦВМ обычно оказываются ограниченными, так как не возникает необходимость рассчитывать значения функционалов для всевозможных ошибок ГНС, но возникают ограничения сходимости.

Рассмотрим эвристические поисковые алгоритмы одномерной КЭНС для оценки местоположения. Будем рассматривать равномерный горизонтальный полет. Пусть к текущему моменту времени t, когда принимается решение, проведено N измерений геофизического поля Z1=f*1, ...,ZN=f*N (рис. 2,а) и пусть в памяти БЦВМ записана карта поля fn(x) (рис. 2,6).

В поисковых КЭНС производится запись в оперативное запоминающее устройство измеренных значений поля Zi=f*iи выдаваемых ГНС координат местоположения x*i в момент проведения измерений. Пусть x*n - коорди­ната местоположения в момент последнего замера.

Эвристический алгоритм поисковой КЭНС сводится к следующему. Наз­начается доверительный интервал 2Δmax с центром в точке х*N. где Δmах-возможная максимальная ошибка ГНС. Задается и

-8-