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

На правах рукописи

□03053035

ЖИДКОВА НАТАЛЬЯ В А Л Е Р Ь Е В ....

АЛГОРИТМЫ ИДЕНТИФИКАЦИИ И КОМПЕНСАЦИИ ИНФОРМАЦИОННЫХ НАРУШЕНИЙ ИЗМЕРИТЕЛЬНОЙ СИСТЕМЫ

05.13.01 - Системный анализ, управление и обработка информации (по техническим наукам)

АВТОРЕФЕРАТ

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

Н. Новгород - 2007 г.

003053035

Работа выполнена в Арзамасском Политехническом Институте (филиале) Нижегородского Государственного Технического Университета.

Научный руководитель:

Доктор физико-математических наук, профессор Пакшин П.В.

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

Доктор технических наук, профессор Федосенко Ю.С.

Кандидат технических наук, доцент Волков Н.В.

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

ОАО АНПП «ТЕМП-АВИА» г. Арзамас

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

Защита состоится «.¿б чМйЛтп^- 2007 г. в диссертационного совета Д 212.16o.05 Нижегородского Государственного Технического Университета по адресу: 603600, Нижний Новгород, ул. Минина, 24, НГТУ, корпус / , аудитория .

С диссертацией можно ознакомиться в научно-технической библиотеке Нижегородского Государственного Технического Университета.

Автореферат разослан « » фЛ^аи-,2007 г.

Ученый секретарь ЬгАонУ к.т.н, Иванов А.П.

диссертационного совета ^»ХАШсг

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

Для достижения поставленной цели в диссертационной работе решались задачи:

- исследование моделей характерных информационных нарушений ИИС;

- разработка эффективных алгоритмов идентификации и компенсации постепенных информационных нарушений ИИС;

- разработка принципов построения цикла идентификации и компенсации информационных нарушений ИИС с использованием оптимизации вероятности ошибочных решений при обнаружении информационных нарушений;

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

Научная новизна. Научная новизна диссертационной работы заключается в следующем:

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

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

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

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

Также результаты работы могут быть использованы для измерительных информационных систем, в которых актуальным является высокодостоверное обнаружение и гарантированная компенсация постепенных отказов измерительных преобразователей ИИС. Методика экспериментальных исследований эффективности алгоритмов ИИС на основе пакета Matlab Simulink полезна для инженерного применения при проектировании ИИС.

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

Внедрение работы. Основные результаты диссертационной работы внедрены:

- ОАО АНПП «ТЕМП - АВИА» в составе специального программного обеспечения вычислителя блока | определения ориентации комплекса наземной отработки интегрированной навигационной системы;

- предприятием ООО «ЭЛЬСТЕР Газэлектроника» в программное обеспечение ИИС контура управления технологическим процессом испытательной температурной установки;

- в учебный процесс в Арзамасском политехническом институте филиале НГТУ на кафедре "Прикладная математика" по дисциплине "Теория навигационных систем" специальности 010501.65 "Прикладная математика и информатика".

Апробация работы. Основные положения и результаты диссертации опубликованы в статье журнала "Радиотехника" №6 2004, межвузовских сборниках "Прогрессивные технологии в машино- и приборостроении" Н.Новгород в 19992004гг.; доложены автором и обсуждены на 12 и 13-й Международных конференциях "Saint Peterburg International conference on integrated navigation systems", Санкт-Петербург, ЦНИИ «Электроприбор», май 2005, 2006 гг., 7-й конференции «Навигация и управление движением», Санкт-Петербург, ЦНИИ «Электроприбор», март 2005г.; 7-й ВНТК "Информационные технологии в промышленности и учебном процессе - 2004», Москва-Арзамас, сентябрь 2004г.; 10-й нижегородской сессии молодых ученых «Технические науки», Дзержинск, март 2005г, а также в ряде Всероссийских НТК, проводимых в НГТУ.

Структура и объем работы. Диссертация состоит из введения, пяти глав, выводов, списка литературы и приложения. Работа изложена на 148 страницах машинописного текста, содержит 40 рисунков и 110 наименований литературных источников.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

зующегося управления, развитой Саридисом Дж. и в исследованиях русских ученых - Цыпкина Я.З., Фельдбаума A.A.

Развитие методов идентификации параметров стохастических систем в реальном времени проведено группой ученых военно-воздушной инженерной академии им. Н.Е.Жуковского, Раменского приборостроительного конструкторского бюро: Чернодаровым A.B., Платоновым А.Ю., Будкиным B.J1., Голиковым В.П. и др.

Исследования, проводимые в Санкт-Петербурге Дмитриевым С.П., Колесовым Н.В., Осиповым A.B., Степановым O.A., Кошаевым Д.А., Амосовым О.С. позволяют сформулировать локальные требования к точности оценок и достоверности контроля.

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

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

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

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

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

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

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

В основу положена матричная модель e(t) информационных нарушений, которая является аддитивной составляющей уравнения состояния ДС.

X{t) = FX(t) + е, (t) + Gw[t), i = (VV.

При этом расширенный вектор состояния, X - (Хт ,ет)т имеет уравнение динамики

X(t) = F'X(t) + G'w(t), i = ÖJ4-

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

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

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

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

В сложных измерительных системах, где вероятность возникновения информационных нарушений достаточно высока (как показывают статистические наблюдения на уровне 10"3'), ошибки контроля и диагностики требуется обеспечить на малом уровне практически равном апостериорной вероятности возникновения нарушений (= 10").

Апостериорные вероятности ошибок контроля 1-го и 2-го рода

ас = Р(В'/Н0)= , ' , Рс=Я(Я/Я0) =-П"о>В)РУ»-.(1)

Р(Н0/В')Р(.В,) + Р{Н0/В)Р(В) с Р(Н0/В)Р(В) + Р(Н0/В')Р(В')

где Н0, Н'о - гипотезы об отсутствии и наличии отказа; В, В' - события наличия и отсутствия отказа в системе; Р(Н0/В), Р(Н0/В% Р(Н'0/В), Р(Н'0/В'), - условные априорные вероятности; Р(В), Р(В') - безусловные априорные вероятности. Отказ измерительной системы наступает при возникновении хотя бы одного элементарного информационного нарушения. Полнота контроля, в этом случае, имеет вид

№ = (2)

где 1,..., т - множество контролируемых информационных нарушений,

{Л,},7 = т + 1,..., п - множество неконтролируемых нарушений. Самоорганизующаяся ИИС оптимизируется по гарантированному по вероятности показателю качества. Используется функционал вероятности

/>Вс(н) = Р{Ф(к,0))<Рс} (3)

Ррс(и) - вероятность события, состоящего в том, что случайная величина Ф(и,со) не превышает заданного значения (Зс ((Зс. - вероятность необнаруженного отказа). Эта вероятность функционально зависит от стратегии управления и процессом компенсации. На основе функционала вероятности получен функционал квантили

Ф6с(м) = ш1п{вс.Р,с(й)>бс}, (4)

где Ос - заданная вероятность (вероятность ошибки контроля первого рода).

Функционал квантили обозначает минимальный уровень рс не превышае-

мый случайным значением Ф(г/,со) с вероятностью Ррс(и) > ас, т.е. это гарантированное по вероятности значение функционала Ф(г/,со).

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

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

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

Рисунок 1. - Функционирование ПИС при идентификации и компенсации.

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

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

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

Х(к+1) = Р{Х(к), К(А)Л[/:,и|(т1),г)2("г2),...г)л-(тл.)]}, (5)

где вектор у зависит от N изменяющихся параметров, величина изменения г'-го параметра и, и время его появления т, (г = являются случайными величинами.

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

У{к) = 1г \Х(к), Щк), у [к, и,(/»,), ъ2(т2),..., (6)

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

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

пенсацию и оценить расширенный вектор состояния системы Х(к).

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

1 Общий алгоритм идентификации и компенсации информационных нарушений.

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

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

На рисунке 2 представлена структурная схема контроля, самоорганизующегося управления и оценивания, реализующая разработанный алгоритм. Приняты следующие обозначения: щк) - заданное детерминированное воздействие; у(к) - вектор измерений; Е\к+ик) - обновляющий процесс; Л^вд, - детерминированная и стохастическая составляющие вектора состояния объекта наблюдения; Хл(к} - вектор оценок; Г, й, К, С -матрицы коэффициентов; £>ь £>2 -дисперсии шумов измерения контролируемого канала в текущем и

предыдущем цикле контроля; Ь -масштабирующий коэффициент допуска контроля; <7 - вероятность отсутствия информационного нарушения.

Схема содержит два основных блока: 1) блок оценивания, 2) блок контроля.

В блоке идентификации нарушений используется аппроксимация зависимостей параметров алгоритма контроля от вероятностей ошибочных решений ас и (Зс. Эти зависимости получены на этапе экспериментальных исследований ИИС в условиях изменяющейся помеховой обстановки путем многофакторного эксперимента. Применение методики ¡многофакторного эксперимента позволило получить апроксиммирующие зависимости с заданной точностью. Выходные параметры блока идентификации нарушений: Ог /£), (дисперсии шума в текущем состоянии к дисперсии шума предыдущего состояния), Ь и д, определяются оптимальным образом - с гарантированным уровнем значимости.

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

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

Рисунок 2. - Структурная схема контроля, управления и оценивания.

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

По результатам идентификации настраиваются параметры: с = <з\/а] , Ь, д и коэффициенты матриц оценивающего устройства. Оперативный блок контроля, представленный в левой части схемы реализует сравнение обновляющего процесса Е = У - СХ, поступающего с блока оценивания с вычисленным пороговым значением

£„„„ = о?, . 2 In

■Kv (7)

„ (1-£)[1-<?(!/,)]_

где q^/ß - вероятность отсутствия информационного нарушения на j-м цикле контроля реального времени; е - доверительная вероятность при допусковом контроле процесса; Ъw - масштабирующий коэффициент допуска контроля; D4j) = а?0), D2U) =^2(7) " Дисперсии шумов измерения контролируемого канала в текущем и предыдущем цикле контроля.

Управление параметром Ъ (масштабирующим коэффициентом порога контроля) основано на использовании аппроксимирущих таблиц вероятностей ошибок контроля (ас + ßr = f{a\/<s], b, q)).

Экспериментальным путем определяются критические значения Сц, = ст;/ст,2 , ниже которых нецелесообразно выполнять управление параметрами (неразличимые шумы): с < скр - нет эффективности или отрицательная эффективность, с > скр - есть эффективность от управления параметрами.

2 Идентификация и компенсация информационных нарушений в реальном времени.

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

В качестве критерия оптимальности используется функция p(ß), которая связана с оценкой JE, вектора ошибок ИИС х, в г-й момент времени соотношением

*, = argmin ¿P(ß,) (8)

где p(ß) = -Inf (ß);ЛР)" функция плотности вероятности; ß, = v,/a, - нормированная невязка; v, - невязка между реальным у„ и прогнозируемым у, = С,х, значениями наблюдений; С, - вектор состава измерений; а, - нормирующий коэффициент.

При динамике движения x^y-Fx(k-\)-Gu^-\)=Q решением задачи (8) является алгоритм робастного оценивания, реализуемый на основе аппроксимационных зависимостей. Для идентификации информационных нарушений необходимо знать условную плотность вероятности flX^/У*). Эта плотность вероятности вычисляется по формуле Байеса для законов распределения

^у^ЩШШШ^ (9)

ДГда/Г*-1)

Плотность распределения оценок экстраполяцииЛХ^/У ) при отсутствии информационных нарушений аппроксимирована гауссовской плотностью. При нарушении гауссовости невязок появляется несоответствие априорно принятых параметров оценивающего фильтра реальному измерительному процессу.

Для определения отклонений параметров уу и РУ от их реальных значений используется адаптивный алгоритм распознавания гипотез на основе вычисления гауссовских функций плотности вероятности

Вероятность вычисляется по формуле Байеса:

«,/,)= = , (10) ,у =!)?(,/*-1)

где 9(1«), (1 - (](\!к)) - вероятности соответственно наличия и отсутствия информационного нарушения на к-1 шаге.

В этом случае полагаем, что у(Ч неизвестная, но постоянная на интервале наблюдения величина. Априорная информация о величине уда задана ее распределением. Оптимальная оценка вектора состояния вычисляется в соответствии с соотношениями фильтра Калмана, апостериорная вероятность отсутствия информационного нарушения системы определяется по рекуррентной формуле (10).

Начальное условие для параметра уда при полной априорной неопределенности соответствует распределению параметра q^^\k.\) = <7 с равномерной плотностью на интервале 0,1. Математические преобразования дают конечное рекуррентное выражение для <7(]/ц:

— г\

__9(1/,ы)./а)_

9(11к)~- „1 п - 1/■' • С)

В реальном алгоритме (11) непрерывная плотность реализуется рядом дискретных точек. Алгоритм обнаружения с использованием вычисленной вероятности <7(|/<) удобно представить в виде последовательно выполняемых этапов.

1. При к = 1 начальное значение 9(1/определяются исходным <7;

2. Вычисляется д(1/ к) по формуле (11);

3. Сравнивая с пороговым значение Епор{к), определяется, что

- при Е(кц_ 1) < £пор((.) сбоя или внезапного отказа нет,

- при Е(ык_ 1) > Епор(ц сбой или внезапный отказ присутствует;

4. Если к < N,10 ожидание следующего дискретного шага к = к +1 (прерывания от таймера) и переход к циклу 2;

если к>Ы,ю окончание цикла контроля реального времени и выдачад^пу

Полученное значение вероятности отсутствия информационного нарушения 9(1/ к) используется далее в алгоритме реального времени для формирования порога контроля компаратора. Компаратор вырабатывает логический сигнал "0" или "1" в зависимости от результата сравнения отклонения контролируемого параметра с порогом.

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

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

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

где Л, - корреляционная функция, определяемая по экспериментальным данным: ту =у'Д А/, =- . Экспериментальные данные оценок ошибок ДПИ х, для определения RJ могут быть получены путем совместной обработки информации ИИС с внешним измерителем, например, спутниковой навигационной системой.

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

3 Идентификация и компенсация информационных нарушений на основе постданных.

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

При постобработке данных решаются следующие задачи:

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

- уточнение расширенного вектора ошибок ИИС путем адаптивно-робастной обработки зарегистрированных оценок в «обратном» времени;

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

(15)

(16)

(17)

В этом случае функция робастного правдоподобия аналогична (8)

4 Экспериментальные исследования эффективности алгоритма идентификации и компенсации информационных нарушений

Разработанный алгоритм реализован в системе управления технологическим процессом испытательной температурной установки предприятия ООО "ЭЛЬСТЕР Газэлектроника". При этом:

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

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

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

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

Моделировались характерные информационные нарушения типа скачкообразного изменения интенсивности измерительных шумов ДПИ.

Оценивание производилось с помощью настраиваемого фильтра Калмана Результаты моделирования предлагаемого алгоритма отказоустойчивого оценивания и контроля в сравнении с типовыми для измерительной части системы управления технологическим процессом представлены на рисунках 3, 4. На рисунке 3 представлены графики оценок температуры при информационном нарушении типа изменения шумов измерения от интенсивности 0,2 до 100. На рисунке 4 графики обозначены для. 1 - при оценивании процесса предлагаемым алгоритмом; 2 -фильтром Казаринова; 3 - обобщенным фильтром Калмана. Экспериментальные данные подтверждают наличие областей параметров устройства контроля, дающих максимальную достоверность контроля. Установлена зависимость суммы вероятностей ошибок контроля от порога контроля при рассматриваемых информационных нарушениях. Получены оптимальные значения порога контроля Епор.

Произведен сравнительный анализ эффективности предложенного варианта ИИС с устройством идентификации и компенсации информационных нарушений и вариантов: 1) обобщенного фильтра Калмана, 2) оказоустойчивого фильтра Ка-

/»г

«.I - - - ¡- - - - ; - - - ; — ----i----J-- — ! — -I - - - -j

0 100 21*1 10IÎ 40(1 '.ll'i (.1*1 7(H) К(Ю 41.)

В)

Рисунок 3. - Оценки температуры фильтрами: а) - с помощью ОФК; б) - с помощью отказоустойчивого фильтра Казаринова; в) - на основе предлагаемого алгоритма.

».со'

15 ----1---" ^----^

О 10« 200 300 400 500 <Î00~ 700 коо <М 'Л

Рисунок 4. - Дисперсия ошибок оценивания.

заринова Исследование эффективности проведено для различных вариантов показателя с = ст;/а-^ , где й\ = <з], В2 = а \.

Результаты подтверждают улучшение эффективности фильтрации в среднем на 17% по сравнению с обобщенным фильтром Калмана и на 5% по сравнению с отказоустойчивым фильтром Казаринова.

При экспериментальном исследовании алгоритма идентификации и компенсации информационных нарушений в системе управления технологическим процессом испытательной температурной установки предприятия ООО «ЭЛЬСТЕР Газэлектроника»

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

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

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

5 Математическое обеспечение оптимизации параметров схемы контроля

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

а. + Рг =0-<7(|Д))

¡Ме)с1е + [/¡(е)«/е

"<7(1,А) Ш6)^, (19)

где <7(и), - вероятность отсутствия информационного нарушения на А:-шаге.

При неизвестных плотностях вероятности/¡(е) и /2(е) минимизация суммы (ас + рс) сводится к построению условных аппроксимирующих зависимостей от параметров влияющих на пороговое значение Епор.

Для расчета и построения аппроксимирующих таблиц для с = а,2/а(2 ,Ъиц, соответствующих минимуму суммы (ас + рс) разработан алгоритм и программа, предназначенные для широкого класса линейных объектов наблюдения (различной размерности матриц коэффициентов).

В четвертой главе разработан алгоритм контроля, диагностики и компенсации информационных нарушений, возникающих в интегрированной навигационной системе (ИНС) в условиях изменяющейся помеховой обстановки и алгоритм расчета и оптимизации параметров схемы контроля. Модели информационных нарушений соответствуют Марковским моделям, рассмотренным в главах 2 и 3, для навигационной системы.

Для решения задачи обнаружения и компенсации информационных нарушений необходимо:

- установить событие информационного нарушения в ИНС;

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

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

- осуществить компенсацию возникших информационных нарушений управлением параметрами ИНС.

В комплексе наземной отработки при исследовании ИНС в динамике - оборудование размещается на автомобиле, который перемещается по заданному маршруту. Измеренные данные параметров движения объекта, полученные И11С в результате работы комплекса наземной отработки, записаны в файлы. Количество данных соответствует периоду времени 0 - 1692 е., т.е. по 16920 значений каждого параметра.

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

Вспомогательные измерения формируются в виде (ср - ф,) = -Дфс + Дф, (А. -XJ = = -АХС + АХ. Вводя векторы д: - ((л-) )', w -- ((w')\(w")')r, так чю х' -(Дфс, ЛХС)Г, х" (6ф, ЬХ, б К, 5К)' сформулирована задачу фильтрации шее in мерного вектора х, описываемого с помощью уравнений х = F(x, vv) цо двухмерным измерениям у ~ Нх + v.

Искомые координаты определяются с помощью оценок Дф(..Дл[ по формулам: Ф = <р, - Дф,, X = X., - АХt.

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

Вычисление Дф,,, ДА., производится в соответствии с уравнениями

„ cosAf - Ksin К • . .

Дф, =-8К--б.К + ^ДДф-Лф^ -5ф) +

+ ku (АХ - AXr -+ (ДV-5V) + к14 (АК -5К)

(20)

- -УътК „ бшА: - КсоэЛГ - , . „

АХЛ =-Дф4+--8У +-5АГ + ^,(Дф-Дф1 -5ф) +

Я.созр /г,соэср Я,соэф " (21)

+ к22 (АХ-АХ1 - 8Х) + кг,( А V - 8 V) + /I,,(АК-8К) Детерминированную часть соотношений вычисляют функциональные блоки. Обратную связь в фильтре А:п(Дф-Дфс-5ф) + А-,,(ДЛ-Д^ ~5Х), ¿2|(Дф-Дф( --5ф) + £22(ДА.-АХС-8Х) вычисляют блоки обновляющих процессов Дф--Лфс--6ф. ДХ.-ДХС -5Х и функциональные блоки (подробная схема вычислений приведена в диссертации). Входными переменными алгоритма контроля являются измерения (для канала широты - ф,„ , для канала долготы - >.,„). Этот блок идентифицирует значения дисперсий шумов (отдельно для широты и долготы).

При использовании разработанного алгоритма на основании данных, предоставленных ОЛО АНПП «ТЕМП - АВИА», проведены экспериментальные исследования посредством компьютерного моделирования и дан анализ эффективности предложенных решений в сравнении с традиционно-применяемыми.

Экспериментальное исследование эффекта фильтрации проводилось на основе численных оценок дисперсии ошибок. Результаты дисперсий ошибок оценивания при заданных дисперсиях ошибок измерения Д,фюм = 2,99-10'8 рад2 и Даим = 2,МО"7 рад2 составляют:

- для предлагаемого устройства оценивания - Ц,ф = 2,67-Ю"12 рад2, Д.>.= 5,9 10'4

. рад2;

- для фильтра Казаринова -ЛСФ= 7,45-10'12 рад2, Ол = 7,5-10"' рад2;

- для обобщенного фильтра Калмана - Д,ф = 9,7 10'" рад2, £)л= 1,52-10 рад2.

Эффект применения предлагаемого алгоритма контроля и адаптивного оценивания состоит в улучшении информационной точности:

- для ф по сравнению с отказоустойчивым фильтром Казаринова в 1,67 раза и обобщенным фильтром Калмана в 6,027 раза;

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

|—■—-:-1-:_и:с ринова в 1,127 раза и обоб-

I) 2<Ю 4ГН1 61Ю N110 ИНН) , ,г

„ „ ,,/-«- щенным фильтром Калмана в

Рисунок 8. - Дисперсия оценивании ср, к: 1 - обобщен- , г

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

- 1,5 10~6 рад' Оптимальные значения порога

г - вероятности ошибок контроля,

контроля Епор представлены пунктирнои линиеи. 1 1

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1 Проведен анализ моделей характерных информационных нарушений ИИС.

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

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

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

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

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

о, • ю рад-

Полу ченные, при моде-

Щ-

ОСНОВНЫЕ РАБОТЫ ПО ТЕМЕ ДИССЕРТАЦИИ

[1] Волкова (Жидкова) Н.В. Исследование датчиков первичной информации методом статистических испытаний Монте-Карло. - Сборник материалов Всероссийской молодежной НТК «Приборостроение в аэрокосмической технике». -Арзамас, 1999. с.24 - 28.

[2] Жидкова Н.В. Обзор исследований, связанных с моделированием состояния и обеспечением безотказной работы динамических систем. - Сборник статей по материалам ВНТК. Технологии в машино- и приборостроении. - Н. Новгород, 2003. с. 465 -469.

[3] Жидкова Н.В. Синтез оптимальной динамической системы по неявному функционалу качества. - Материалы 2-й региональной молодежной НТК «Будущее технической науки Нижегородского региона». - НГТУ, Н. Новгород, 2003.

[4] Жидкова Н.В., Волков В.Л. Достоверность контроля состояний нелинейной адаптивной измерительной системы - Прогрессивные технологии в машино-и приборостроении «ПТ-2004». - Материалы ВНТК Н.Новгород - Арзамас, 2004.

[5] Жидкова HB., Волков В.Л. Нелинейная адаптивная измерительная система с заданной достоверностью контроля состояний. - Сборник материалов седьмой научно-технической конференции «Информационные технологии в промышленности и учебном процессе». - М.: ООО Риасофт ЛТД, 2004, с,23 - 24.

[6] Жидкова Н.В , Волков В Л. Алгоритм адаптивной фильтрации измерительной системы с марковской моделью состояний. - Тесисы докладов ВНТК «Информационные системы и технологии ИСТ-2004». - НГТУ, Н.Новгород, 2004.

[7] Жидкова Н.В. Взаимосвязь точности оценок и полноты контроля состояний марковской адаптивной измерительной системы. - Тесисы докладов III Молодежной НТК «Будущее технической науки». - НГТУ, Н.Новгород, 2004. с. 94 -95.

[8] Жидкова Н.В. Математическое обеспечение для проектирования и исследования динамических систем. - Радиотехника, №6, 2004.

[9] Жидкова Н.В. Моделирование нелинейной адаптивной измерительной системы в Matlab Simulmk. - Сборник материалов седьмой научно-технической конференции «Информационные технологии в промышленности и учебном процессе». - М.: ООО Риасофт ЛТД, 2004. с. 24-28

[10] Жидкова HB., Волков В.Л Идентификация и компенсация отказов измерительной системы. — Гироскопия и навигация, № 4, 2005.

[11] Zhidkova N.V., Volkov V.L. The algorithm of self-organizing process for identification and compensation of measuring system failures. - 12th Saint Petersburg international conference on integrated navigation systems. - SPB.: Elektropribor, 2005. p. 126 - 127.

[12] Zhidkova N.V., Volkov V.L. The robust algorithm of identifications and compensation of information refusals of a navigation system. - 13th Saint Petersburg international conference on integrated navigation systems. — SPB.: Elektropribor, 2006. p. 131 - 135

Подписано в печать 01 02 2007 г Формат 60x84/16 Уел печ листов 1 Бумага офсетная Печать офсетная Гарнитура Times New Roman Заказ №711 Тираж 100 экз Отпечатано 8 ОАО «Арзамасская типография» 607220 г Арзамас Нижегородской области, ул Пландина, 8